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

Max Rudolph maxwellr at gmail.com
Fri Feb 3 20:36:12 PST 2017

There is a pull request here containing a jupyter notebook that calculates
and compares critical rayleigh number using ASPECT and compares it with the
analytic solution from Turcotte and Schubert.


You can also verify which approximation is being used:

$ grep Formulation output/parameters.prm
subsection Formulation
  set Formulation          = boussinesq approximation # default: custom

Hope this helps.


On Fri, Feb 3, 2017 at 3:53 PM, Elbridge Gerry Puckett <
egpuckett at ucdavis.edu> wrote:

> Hi,
> Timo, Wolfgang Juliane, Rene, Max Rudolf and CIG researchers know what I
> am referring to.
> My understanding is that Timo and Juliane have fixed the problem with the
> original - incorrect - Boussinesq approximation.
> How does one make certain that we are approximating solutions of the
> correct incompressible Mantle Convection Equations when running the current
> ASPECT master branch.
> Also John Naliboff suggested that I also ask if the defaults are still no
> internal heating, adiabatic heating, etc. etc. In other words, the only
> source of heat is the boundary conditions on the top and bottom walls (in
> 2D) assuming no heat flux BCs on the side walls.
> 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.
> Thanks!
> - Gerry
> --
> _______________________________________________
> 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/20170203/330459ac/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: EGP_AVATAR.png
Type: image/png
Size: 4465 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20170203/330459ac/attachment.png>

More information about the Aspect-devel mailing list