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

[Eric Heien]

Hi Timo,

Sounds good.  Let me know if Ted or I can provide anything else.  I know Ted's been working on some possible fixes with an applied math professor here (Gerry Puckett) so we can discuss these sometime.

-Eric

On Oct 17, 2013, at 4:15 PM, Timo Heister wrote:

> Eric,
> 
> cool, thanks! I will look into this. I think we can fix the jump from
> the first timestep and that should improve the situation later on.
> 
> On Thu, Oct 17, 2013 at 6:46 PM, Eric Heien <emheien at ucdavis.edu> wrote:
>> 
>> 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.
>> 
>> 
>> _______________________________________________
>> Aspect-devel mailing list
>> Aspect-devel at geodynamics.org
>> http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
> 
> 
> 
> -- 
> Timo Heister
> http://www.math.clemson.edu/~heister/
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel