# A cookbook for using adiabatic conditions and a gravity # profile generated by a mineral physics toolkit and read in # through an ascii data file. In this example, material properties # are based on a Birch-Murnaghan equation of state. set Dimension = 2 set Use years in output instead of seconds = true set End time = 2.63e9 set Output directory = output-Mariana-07-12 subsection Material model set Model name = ascii reference profile subsection Ascii reference profile set Thermal viscosity exponent = 10.0 set Viscosity prefactors = 1.0, 0.1, 1.0, 10.0 subsection Ascii data model set Data directory = $ASPECT_SOURCE_DIR/data/adiabatic-conditions/ascii-data/ set Data file name = BurnMan_profiles.txt end end set Material averaging = harmonic average end # The geometry is a spherical shell with the inner and # outer radius of the Earth's mantle. subsection Geometry model set Model name = spherical shell subsection Spherical shell set Inner radius = 3481000 set Outer radius = 6336000 set Opening angle = 90 end end subsection Adiabatic conditions model set Model name = ascii data subsection Ascii data model set Data directory = $ASPECT_SOURCE_DIR/data/adiabatic-conditions/ascii-data/ set Data file name = BurnMan_profiles.txt end end # The gravity model reads an ascii data file that contains a gravity # profile consistent with the material properties computed using the # Birch-Murnaghan equation of state. It automatically uses the same file # specified in the adiabatic conditions model. subsection Gravity model set Model name = ascii data subsection Ascii data model set Data directory = $ASPECT_SOURCE_DIR/data/adiabatic-conditions/ascii-data/ set Data file name = BurnMan_profiles.txt end end # The model uses the anelastic liquid approximation. subsection Formulation set Formulation = anelastic liquid approximation end subsection Boundary velocity model set Zero velocity boundary indicators = inner set Tangential velocity boundary indicators = top, left, right end subsection Boundary temperature model set Fixed temperature boundary indicators = top, bottom set List of model names = spherical constant subsection Spherical constant set Inner temperature = 3500 set Outer temperature = 273 end end # The initial temperature is an adiabatic profile, which is # taken from the adabatic conditions model. In this case, this is # is the ascii data model that loads a file containing profiles # corresponding to the Birch-Murnaghan equation of state. # In addition to the adiabatic profile, there are thermal boundary # layers at the sirface and the core-mantle boundary. subsection Initial temperature model set Model name = adiabatic subsection Adiabatic set Age top boundary layer = 5e7 end end subsection Mesh refinement set Initial global refinement = 5 set Initial adaptive refinement = 0 set Time steps between mesh refinement = 0 #this means mesh will never change end subsection Heating model set List of model names = adiabatic heating, shear heating subsection Adiabatic heating set Use simplified adiabatic heating = true end end subsection Postprocess set List of postprocessors = visualization, velocity statistics, temperature statistics, heat flux statistics, depth average subsection Visualization set Output format = vtu set Time between graphical output = 0 set Number of grouped files = 0 set List of output variables = material properties, gravity, adiabat, nonadiabatic temperature, nonadiabatic pressure end subsection Depth average set Time between graphical output = 1e6 set Number of zones = 150 end end subsection Checkpointing set Time between checkpoint = 1800 end subsection Solver parameters subsection Stokes solver parameters set Linear solver tolerance = 1e-7 set Maximum number of expensive Stokes solver steps = 10000 set GMRES solver restart length=10000 end end