[aspect-devel] Compositional fields with different densities and viscosities

[Glerum, A.C. (Anne)]

Hi all,

On your advice, I've produced my own material models with multiple materials. This all works fine, thank you. I was just wondering if you will make it possible to output parameters in the statistics file per composition (other than the min/max/global composition values)? For example, I would be interested in the Vrms and dissipation of energy of only the downgoing slab (where the slab is represented by one compositional field). This would greatly help with the benchmark model I am working on.

Kind regards,

Anne


Anne Glerum, MSc | PhD candidate | Department of Earth Sciences | Utrecht University|  Budapestlaan 4, 3584 CD Utrecht | Room Z.204 | A.C.Glerum at uu.nl<mailto:D.J.J.vanHinsbergen at uu.nl> |



On Dec 10, 2012, at 5:37 PM, Juliane Dannberg wrote:

Hi Anne,

you are right, you can only apply one density difference in the Simple material model.
Until now, nobody needed as many as 8 different compositional fields, so until now aspect does not have many complex material models.
But I think the easiest and fastest way for you would be to write your own material model. Just copy the Simple model and then include all the input parameters you would like to have. In every material model there is the function density(), and that is the place where you want to change something. As the function gets the compositional fields as input parameter, you can easily apply different densities for different compositional fields there.


On 12/10/2012 05:18 PM, Glerum, A.C. (Anne) wrote:
Hi all,

For my PhD project I am working on modelling subduction with Aspect. I would like to start with a simple linear viscous subduction benchmark including eight materials. The materials all have different densities and different viscosities. As far as I understand, with material model "Simple model" it is possible to define eight compositional fields, but not to give them all different densities and/or viscosities. In the input file one can set "Density differential for compositional field 1" and the Composition viscosity prefactor, but then they apply to all fields; you cannot specify them for each field individually.

Can you tell me whether my understanding is correct, and if so, whether you intend to make it possible to define multiple compositional fields coupled to different densities/viscosities from the input file? I would guess the alternative is to hardwire the materials into the material model myself.

Thank you so much for your help,

Anne Glerum

PS. Cedric Thieulot and I have looked into the Van Keken et al. 1997 benchmark some more after the stabilisation parameters for advection were adjusted; I would like to show you some results later.



_______________________________________________
Aspect-devel mailing list
Aspect-devel at geodynamics.org<mailto:Aspect-devel at geodynamics.org>
http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel

_______________________________________________
Aspect-devel mailing list
Aspect-devel at geodynamics.org<mailto:Aspect-devel at geodynamics.org>
http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/aspect-devel/attachments/20130122/7fc34da5/attachment.htm