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

[Juliane Dannberg]

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
> http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel

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