[aspect-devel] Input parameters for the compressible mantle convection?

[Nan Zhang]

Hi Juliane (or Timo?),

I had a look at these two .prm files for Tan & Gurnis' 07 and King' 10
cases.  I have several questions.

(1) When Timo tested Tan&Gurnis's compressibility with the dissipation
number, what is the difference of Tan&Gurnis from "Simple compressible"? Is
the "Simple compressible" material model good for the Earth's case?

(2) The King2dcompressible test has a subsection like:
********************************************
...
subsection Formulation
set Formulation = EBA # custom
set Mass density approximation = incompressible#implicit adiabatic
set Temperature density approximation = adiabatic
set Buoyancy density approximation = adiabatic pressure
end
...
********************************************
I could not find it in the current ASPECT version. I wonder if this change
is necessary for a compressible convection?

(3) Once I put "shear heating" into the Heating model, is there any other
parameter companion with "shear heating"?

(4) I think that the "subsection Adiabatic conditions model" is necessary
for a compressible flow. However, I think that it is going to conflict to
the "subsection Initial condition", isn't it?

Bests,
Nan

On Thu, Sep 1, 2016 at 7:18 AM, Juliane Dannberg <dannberg at gfz-potsdam.de>
wrote:

> 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 listAspect-devel at geodynamics.orghttp://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/20160901/c9511e47/attachment.html>