[aspect-devel] Progress in writing the mantle convection code Aspect

[Magali Billen]

We are modeling sinking of the oceanic lithosphere into the mantle (subduction). The lithosphere is not just a thermal
boundary layer, but also has compositional layering: the basaltic crust, a harzburgitic layer (mostly olivine) beneath that 
and then normal mantle. Each of these compositions undergoes phase changes at different depths depending on its composition.
Also each layer has different rheologic properties (harzburgite is dryer than mantle so stronger, basalt has a lower viscosity, etc...).
So, this is what we'd like to be able model in Aspect. We've been doing this in CitcomS using tracer particles for each field but
it requires a lot of particles and is computationally expensive. 

It would great to see how the stabilization parameter affected this.
-Magali

On Oct 16, 2013, at 4:42 PM, Timo Heister wrote:

>> We're quantifying them for an AGU poster so we'll have more information then.  If I recall correctly, they are on the order of 5-10% of the gradient (so for a gradient of 10^4, the overshoot would be on the order of 10^2 to 10^3), versus up to 20-30% for CitcomS.  Finer resolution will make them more localized but generally worse in magnitude.
> 
> If you can set up simple to run examples, I am interested in seeing
> them. You can likely remove overshoots by increasing the stabilization
> parameter beta. What is your compositional field supposed to
> represent?
> 
> -- 
> Timo Heister
> http://www.math.clemson.edu/~heister/