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

[Eric Heien]

On Oct 17, 2013, at 11:40 AM, Timo Heister wrote:

> I think the defaults are pretty good as is and I am sure that the best
> parameters depend on the problem.
> It sounded like Eric made an "evil" example where a larger value is
> acceptable and that is why I mentioned the option to increase the
> parameters. I mean if the only measure is the amount of overshoot, we
> can just do an O(1) first order diffusion and smear everything out.
> :-)

I've attached a sample parameter file developed by Ted where the overshoot is visible from an inclusion.  You can see the results at:

http://heien.org/overshoot_inclusion.mp4

This includes the temperature profile through the vertical domain center over time.  I've also attached a graph of the max/min temperatures.  You can see the over/undershoot goes to around 8-10% in the first few steps and decreases soon after, but it never fully disappears.

-------------- next part --------------
A non-text attachment was scrubbed...
Name: inclusion.prm
Type: application/octet-stream
Size: 6457 bytes
Desc: not available
URL: <http://geodynamics.org/pipermail/aspect-devel/attachments/20131017/014a7154/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: PastedGraphic-1.pdf
Type: application/pdf
Size: 18879 bytes
Desc: not available
URL: <http://geodynamics.org/pipermail/aspect-devel/attachments/20131017/014a7154/attachment-0001.pdf>