<div dir="ltr">he Timo,<div><br></div><div>are you implementing a proper levelset (advection of a signed distance function) or some step function?</div><div>i would be very much in favor of implementing a proper levelset.</div>

<div>we have been doing some experiments with that in another code and it gives you a lot of extra functionality to know how far from and on what side of an interface you are. </div><div><br></div><div>cheers</div><div>Thomas</div>

<div> </div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Oct 17, 2013 at 8:33 PM, Timo Heister <span dir="ltr"><<a href="mailto:heister@clemson.edu" target="_blank">heister@clemson.edu</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Magali,<br>
<div class="im"><br>
> We are modeling sinking of the oceanic lithosphere into the mantle (subduction). The lithosphere is not just a thermal<br>
> boundary layer, but also has compositional layering: the basaltic crust, a harzburgitic layer (mostly olivine) beneath that<br>
> and then normal mantle. Each of these compositions undergoes phase changes at different depths depending on its composition.<br>
<br>
</div>So that means you have basically three different distinct states A,B,<br>
and C (basaltic, harzburgitic, normal) and not a contiguous range (for<br>
example an oxygen concentration)?<br>
If that is the case then we can do quite a bit better compared to what<br>
is implemented in Aspect right now. At the moment we only assume that<br>
the field is a continuous function that can assume any value. My<br>
student is working on taking advantage when the users tells us that a<br>
field is a levelset (encoding distinct states). When we get that<br>
working you can specify for each individual field, if it is a<br>
level-set (encoding yes/no) or a continuous function (like it is right<br>
now).<br>
<div class="im"><br>
> It would great to see how the stabilization parameter affected this.<br>
<br>
</div>We are also experimenting with one additional parameter right now that<br>
does increase the smoothing at very steep gradients. (still work in<br>
progress)<br>
<div class="HOEnZb"><div class="h5"><br>
--<br>
Timo Heister<br>
<a href="http://www.math.clemson.edu/~heister/" target="_blank">http://www.math.clemson.edu/~heister/</a><br>
_______________________________________________<br>
Aspect-devel mailing list<br>
<a href="mailto:Aspect-devel@geodynamics.org">Aspect-devel@geodynamics.org</a><br>
<a href="http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel</a><br>
</div></div></blockquote></div><br></div>