<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix">Hi Anne, <br>
<br>
you are right, you can only apply one density difference in the
Simple material model. <br>
Until now, nobody needed as many as 8 different compositional
fields, so until now aspect does not have many complex material
models. <br>
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. <br>
<br>
<br>
On 12/10/2012 05:18 PM, Glerum, A.C. (Anne) wrote:<br>
</div>
<blockquote
cite="mid:86D0572567671044A6E8A89314399D9B321647@ICTSC-W-S202.soliscom.uu.nl"
type="cite">
<meta http-equiv="Content-Type" content="text/html;
charset=ISO-8859-1">
Hi all,
<div><br>
</div>
<div>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
<i>all</i> different densities and/or viscosities. In the input
file one can set "<span class="Apple-style-span"
style="font-family: CMR8; ">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.</span></div>
<div><span class="Apple-style-span" style="font-family: CMR8; "><br>
</span></div>
<div><font class="Apple-style-span" face="CMR8">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.</font></div>
<div><font class="Apple-style-span" face="CMR8"><br>
</font></div>
<div><font class="Apple-style-span" face="CMR8">Thank you so much
for your help,</font></div>
<div><font class="Apple-style-span" face="CMR8"><br>
</font></div>
<div><font class="Apple-style-span" face="CMR8">Anne Glerum</font></div>
<span style="font-family: CMR8; "> <br>
</span>
<div><span style="font-family: CMR8; ">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.</span></div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
Aspect-devel mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Aspect-devel@geodynamics.org">Aspect-devel@geodynamics.org</a>
<a class="moz-txt-link-freetext" href="http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel">http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel</a></pre>
</blockquote>
<br>
</body>
</html>