[aspect-devel] Input parameters for the compressible mantle convection?
Juliane Dannberg
dannberg at gfz-potsdam.de
Wed Aug 31 16:18:59 PDT 2016
Hi Nan,
if you go to the benchmarks folder, there's the directory "tangurnis"
that contains input files for the different approximations tested in the
benchmark of Tan, Gurnis '07 (which is compressible).
For a time-dependent problem, we're still working on the King '10
benchmark cases. Timo has some preliminary results on one of his
branches (https://github.com/tjhei/aspect/tree/density_approximations);
but note that this branch also contains changes to the code itself, so
you will only be able to compare results if you checkout, compile and
run the models with his branch.
Best,
Juliane
On 08/31/2016 12:35 PM, Nan Zhang wrote:
> Hi all,
>
> I tried to do a calculation of a compressible mantle flow, and
> followed the user manual to get an input file like:
> ************************************************
> ...
> subsection Material model
> set Model name = simple compressible # simple
>
> subsection Simple compressible model
> set Reference density = 3300
> set Reference specific heat = 1250
> set Reference compressibility = 4e-12
> set Thermal expansion coefficient = 4e-5
> set Viscosity = 1e22
> end
> end
>
> subsection Heating model
> List of model names = adiabatic heating, constant heating, shear
> heating
>
> subsection Adiabatic heating
> set Use simplified adiabatic heating = true
> end
>
> subsection Constant heating
> set Radiogenic heating rate = 1e-12
> end
>
> subsection Shear heating
> end
> end
>
> subsection Model settings
> set Zero velocity boundary indicators = inner
> set Tangential velocity boundary indicators = outer
> set Prescribed velocity boundary indicators =
>
> set Fixed temperature boundary indicators = inner, outer
>
> set Include adiabatic heating = true
> set Include shear heating = true
> set Include latent heat = false
> end
> ...
> ************************************************
>
> This input can be run. My question is if any compressible benchmark
> case with ASPECT can be compared for testing my input parameter?
>
> Thanks!
> Nan
>
>
> _______________________________________________
> 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/20160831/83864d6d/attachment.html>
More information about the Aspect-devel
mailing list