[aspect-devel] convection in a (virtual) box

Juliane Dannberg judannberg at gmail.com
Thu Jan 4 13:25:15 PST 2018

Hi Cedric,

There was indeed a bug in this specific input file because we had set 
the thermal expansivity to 1, but did not use the Boussinesq 
approximation. This was only fixed in late September 
(https://github.com/geodynamics/aspect/pull/1930). So if the virtual box 
version is older, it just does not include this bugfix.

So I guess that just means we should update the virtual box soon.
For the meantime, you can simply update the code in your virtual box to 
the current developer version.


On 01/04/2018 10:02 PM, cedric thieulot wrote:
> hi all,
> for reasons too long to explain here I found myself this afternoon 
> quickly installing and running the VirtualBox software on my Mac.
> I quickly wanted to obtain the vrms and Nu values of the Blankenbach 
> benchmark (‘convection-box.prm’) for various values of Ra.
> Without even editing the prm file, I ran the model and plotted the 
> results, only to find out they are off with regards to the manual (and 
> the actual benchmark values from the paper.)
> And not half a percent off, a good 10-20% off :)
> I believe there probably is a simple explanation: since the thermal 
> expansion coefficient is set to 1 in the prm file the true Boussinesq 
> approximation should be used, yet it is absent from the prm file.
> I therefore suspect that the virtual box was updated in the Spring but 
> not correctly ? (the .prm file on github does contain a Formulation).
> Best,
> Cedric.
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel

Juliane Dannberg
Postdoctoral Fellow, Colorado State University

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20180104/47268b11/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: two.png
Type: image/png
Size: 908902 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20180104/47268b11/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: one.png
Type: image/png
Size: 679817 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20180104/47268b11/attachment-0003.png>

More information about the Aspect-devel mailing list