[aspect-devel] Is the correct Boussinesq approximation now the default in ASPECT?

John Naliboff jbnaliboff at ucdavis.edu
Sat Feb 4 09:10:18 PST 2017

Hi all,

Gerry - Clarification question. Your group was was using a modified form of the equations (relating to density approximation) at some point, yes?  

Did the set of runs from the summer use that custom approximation or was it using ASPECT’s default density approximation? Is it possible the new runs are using ASPECT’s default, while the old ones use your modified equations?

If both sets of models simply used ASPECT’s default for the given material model, from Timo’s email below the approximations should be the same.


John Naliboff
Assistant Project Scientist, CIG
Earth & Planetary Sciences Dept., UC Davis

> On Feb 4, 2017, at 6:11 AM, Timo Heister <timo.heister at gmail.com> wrote:
>>> How does one make certain that we are approximating solutions of the correct incompressible Mantle Convection Equations when running the current ASPECT master branch.
> As Max said, you need to check that
>> subsection Formulation
>>  set Formulation          = boussinesq approximation # default: custom
> is set. Note that we did not change the default behavior in ASPECT.
> The current plan is to release ASPECT 1.5 soon, which allows BA, ALA,
> etc. but doesn't use it by default. We will then switch to working on
> ASPECT 2.0 where we can introduce breaking changes and decide about
> better defaults, not only with respect to formulations.
>>> Also John Naliboff suggested that I also ask if the defaults are still no internal heating, adiabatic heating, etc. etc.
> Defaults did not change, but if you select the BA formulation and have
> incorrect heating settings, we will abort the program. See
> Simulator<dim>::check_consistency_of_formulation() or
> https://github.com/geodynamics/aspect/blob/2aac1f2342460dbce0f0f94180e79c1774cdf299/source/simulator/helper_functions.cc#L1382
>>> The genesis of this particular post is that we are getting startling different results with the same parameter file and presumably, the same code for a problem we're currently working on.
> Interesting. That would be worth exploring further.
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20170204/1026a92c/attachment.html>

More information about the Aspect-devel mailing list