[aspect-devel] Internal heating in aspect (Ludovic Jeanniot)

[Scott King]

Yes.  I’m aware of the manual but manuals and codes don’t always agree, even when they intend to.   Just sayin.   If you read Rene’s comment it sounds like the agreement for the ALA cases is not as good now that the flux is calculated more accurately and he asked if I had any magic.  It is the only thing I could think of...


> On Sep 12, 2018, at 8:50 AM, Timo Heister <heister at clemson.edu> wrote:
> 
> Scott,
> 
>> For the ALA, density should be \rho_{ref} everywhere, except for the buoyancy term.   That’s part of the ALA because the ALA assumes you only keep terms to first order in \mu=\alpha\Delta T.   I know ASPECT had/has the full density in the energy equation, not \rho_{ref}.  I don’t know what you actually use for density in the ALA benchmark problems.
> 
> We have implemented various formulations in ASPECT including ALA to
> mimic the same setup that is used in benchmarks like this. This was
> required to reproduce the results:
> 
> See manual section 2.11 and 2.11.5
> 
> http://www.math.clemson.edu/~heister/aspect-pdf-manual/web/viewer.html#subsection.2.11
> 
> http://www.math.clemson.edu/~heister/aspect-pdf-manual/web/viewer.html#sec:combined_formulations
> 
> Best,
> Timo
> 
> -- 
> Timo Heister
> http://www.math.clemson.edu/~heister/
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel