[cig-commits] commit 2532 by heister to /var/svn/dealii/aspect
dealii.demon at gmail.com
dealii.demon at gmail.com
Mon Apr 14 15:06:02 PDT 2014
Revision 2532
update parameters
U branches/releases/aspect-1.0/doc/manual/parameters.tex
http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=2532&peg=2532
Diff:
Modified: branches/releases/aspect-1.0/doc/manual/parameters.tex
===================================================================
--- branches/releases/aspect-1.0/doc/manual/parameters.tex 2014-04-14 22:05:43 UTC (rev 2531)
+++ branches/releases/aspect-1.0/doc/manual/parameters.tex 2014-04-14 22:05:59 UTC (rev 2532)
@@ -87,10 +87,10 @@
\index[prmindex]{End time}
\index[prmindexfull]{End time}
-{\it Value:} 5.6966627666142017e+300
+{\it Value:} 5.69e+300
-{\it Default:} 5.6966627666142017e+300
+{\it Default:} 5.69e+300
{\it Description:} The end time of the simulation. The default value is a number so that when converted from years to seconds it is approximately equal to the largest number representable in floating point arithmetic. For all practical purposes, this equals infinity. Units: Years if the 'Use years in output instead of seconds' parameter is set; seconds otherwise.
@@ -134,10 +134,10 @@
\index[prmindex]{Maximum time step}
\index[prmindexfull]{Maximum time step}
-{\it Value:} 5.6966627666142017e+300
+{\it Value:} 5.69e+300
-{\it Default:} 5.6966627666142017e+300
+{\it Default:} 5.69e+300
{\it Description:} Set a maximum time step size for the solver to use. Generally the time step based on the CFL number should be sufficient, but for complicated models or benchmarking it may be useful to limit the time step to some value. The default value is a value so that when converted from years into seconds it equals the largest number representable by a floating point number, implying an unlimited time step.Units: Years or seconds, depending on the ``Use years in output instead of seconds'' parameter.
@@ -347,16 +347,16 @@
{\it Description:} Select one of the following models:
-`box': A model in which the composition is chosen constant on all the sides of a box.
-
`spherical constant': A model in which the composition is chosen constant on the inner and outer boundaries of a spherical shell. Parameters are read from subsection 'Sherical constant'.
`initial composition': A model in which the composition at the boundaryis chosen to be the same as given in the initialconditions.
Because this class simply takes what the initial composition had described, this class can not know certain pieces of information such as the minimal and maximal composition on the boundary. For operations that require this, for example in postprocessing, this boundary composition model must therefore be told what the minimal and maximal values on the boundary are. This is done using parameters set in section ``Boundary composition model/Initial composition''.
+`box': A model in which the composition is chosen constant on all the sides of a box.
-{\it Possible values:} [Selection box|spherical constant|initial composition ]
+
+{\it Possible values:} [Selection spherical constant|initial composition|box ]
\end{itemize}
@@ -516,20 +516,20 @@
{\it Description:} Select one of the following models:
-`constant': A model in which the temperature is chosen constant on a given boundary indicator. Parameters are read from the subsection 'Constant'.
-
-`box': A model in which the temperature is chosen constant on all the sides of a box.
-
`initial temperature': A model in which the temperature at the boundaryis chosen to be the same as given in the initialconditions.
Because this class simply takes what the initial temperature had described, this class can not know certain pieces of information such as the minimal and maximal temperature on the boundary. For operations that require this, for example in postprocessing, this boundary temperature model must therefore be told what the minimal and maximal values on the boundary are. This is done using parameters set in section ``Boundary temperature model/Initial temperature''.
-`Tan Gurnis': A model for the Tan/Gurnis benchmark.
+`constant': A model in which the temperature is chosen constant on a given boundary indicator. Parameters are read from the subsection 'Constant'.
`spherical constant': A model in which the temperature is chosen constant on the inner and outer boundaries of a spherical shell. Parameters are read from subsection 'Spherical constant'.
+`box': A model in which the temperature is chosen constant on all the sides of a box.
-{\it Possible values:} [Selection constant|box|initial temperature|Tan Gurnis|spherical constant ]
+`Tan Gurnis': A model for the Tan/Gurnis benchmark.
+
+
+{\it Possible values:} [Selection initial temperature|constant|spherical constant|box|Tan Gurnis ]
\end{itemize}
@@ -1157,16 +1157,16 @@
{\it Description:} Select one of the following models:
-`box': A box geometry parallel to the coordinate directions. The extent of the box in each coordinate direction is set in the parameter file. The box geometry labels its 2*dim sides as follows: in 2d, boundary indicators 0 through 3 denote the left, right, bottom and top boundaries; in 3d, boundary indicators 0 through 5 indicate left, right, front, back, bottom and top boundaries. See also the documentation of the deal.II class ``GeometryInfo''.
-
`spherical shell': A geometry representing a spherical shell or a pice of it. Inner and outer radii are read from the parameter file in subsection 'Spherical shell'.
The model assigns boundary indicators as follows: In 2d, inner and outer boundaries get boundary indicators zero and one, and if the opening angle set in the input file is less than 360, then left and right boundaries are assigned indicators two and three. In 3d, inner and outer indicators are treated as in 2d. If the opening angle is chosen as 90 degrees, i.e., the domain is the intersection of a spherical shell and the first octant, then indicator 2 is at the face $x=0$, 3 at $y=0$, and 4 at $z=0$.
+`box': A box geometry parallel to the coordinate directions. The extent of the box in each coordinate direction is set in the parameter file. The box geometry labels its 2*dim sides as follows: in 2d, boundary indicators 0 through 3 denote the left, right, bottom and top boundaries; in 3d, boundary indicators 0 through 5 indicate left, right, front, back, bottom and top boundaries. See also the documentation of the deal.II class ``GeometryInfo''.
+
`sphere': Geometry model for sphere with a user specified radius.
-{\it Possible values:} [Selection box|spherical shell|sphere ]
+{\it Possible values:} [Selection spherical shell|box|sphere ]
\end{itemize}
@@ -1495,6 +1495,12 @@
{\it Description:} Select one of the following models:
+`adiabatic': Temperature is prescribed as an adiabatic profile with upper and lower thermal boundary layers, whose ages are given as input parameters.
+
+`spherical hexagonal perturbation': An initial temperature field in which the temperature is perturbed following a six-fold pattern in angular direction from an otherwise spherically symmetric state.
+
+`spherical gaussian perturbation': An initial temperature field in which the temperature is perturbed by a single Gaussian added to an otherwise spherically symmetric state. Additional parameters are read from the parameter file in subsection 'Spherical gaussian perturbation'.
+
`perturbed box': An initial temperature field in which the temperature is perturbed slightly from an otherwise constant value equal to one. The perturbation is chosen in such a way that the initial temperature is constant to one along the entire boundary.
`polar box': An initial temperature field in which the temperature is perturbed slightly from an otherwise constant value equal to one. The perturbation is such that there are two poles on opposing corners of the box.
@@ -1503,18 +1509,12 @@
`mandelbox': Fractal-shaped temperature field.
-`adiabatic': Temperature is prescribed as an adiabatic profile with upper and lower thermal boundary layers, whose ages are given as input parameters.
+`harmonic perturbation': An initial temperature field in which the temperature is perturbed following a harmonic function (spherical harmonic or sine depending on geometry and dimension) in lateral and radial direction from an otherwise constant temperature (incompressible model) or adiabatic reference profile (compressible model).
-`spherical hexagonal perturbation': An initial temperature field in which the temperature is perturbed following a six-fold pattern in angular direction from an otherwise spherically symmetric state.
-
-`spherical gaussian perturbation': An initial temperature field in which the temperature is perturbed by a single Gaussian added to an otherwise spherically symmetric state. Additional parameters are read from the parameter file in subsection 'Spherical gaussian perturbation'.
-
`function': Temperature is given in terms of an explicit formula
-`harmonic perturbation': An initial temperature field in which the temperature is perturbed following a harmonic function (spherical harmonic or sine depending on geometry and dimension) in lateral and radial direction from an otherwise constant temperature (incompressible model) or adiabatic reference profile (compressible model).
-
-{\it Possible values:} [Selection perturbed box|polar box|inclusion shape perturbation|mandelbox|adiabatic|spherical hexagonal perturbation|spherical gaussian perturbation|function|harmonic perturbation ]
+{\it Possible values:} [Selection adiabatic|spherical hexagonal perturbation|spherical gaussian perturbation|perturbed box|polar box|inclusion shape perturbation|mandelbox|harmonic perturbation|function ]
\end{itemize}
@@ -2049,24 +2049,22 @@
{\it Description:} Select one of the following models:
-`Steinberger': lookup viscosity from the paper of Steinberger/Calderwood2006 and material data from a database generated by Perplex. The database builds upon the thermodynamic database by Stixrude 2011 and assumes a pyrolitic composition by Ringwood 1988.
-
`latent heat': A material model that includes phase transitions and the possibility that latent heat is released or absorbed when material crosses one of the phase transitions of up to two different materials (compositional fields). This model implements a standard approximation of the latent heat terms following Christensen \& Yuen, 1986. The change of entropy is calculated as $Delta S = \gamma rac{\Delta
ho}{
ho^2}$ with the Clapeyron slope $\gamma$ and the density change $\Delta
ho$ of the phase transition being input parameters. The model employs an analytic phase function in the form $X=0.5 \left( 1 + anh \left( rac{\Delta p}{\Delta p_0}
ight)
ight)$ with $\Delta p = p - p_transition - \gamma \left( T - T_transition
ight)$ and $\Delta p_0$ being the pressure difference over the width of the phase transition (specified as input parameter).
-`simple': A simple material model that has constant values for all coefficients but the density and viscosity. This model uses the formulation that assumes an incompressible medium despite the fact that the density follows the law $
ho(T)=
ho_0(1-eta(T-T_{ ext{ref}})$. The temperature dependency of viscosity is switched off by default and follows the formula$\eta(T)=\eta_0*e^{\eta_T*\Delta T / T_{ ext{ref}})}$.The value for the components of this formula and additional parameters are read from the parameter file in subsection 'Simple model'.
-
-`Tan Gurnis': A simple compressible material model based on a benchmark from the paper of Tan/Gurnis (2007). This does not use the temperature equation, but has a hardcoded temperature.
-
-`table': A material model that reads tables of pressure and temperature dependent material coefficients from files. The default values for this model's runtime parameters use a material description taken from the paper extit{Complex phase distribution and seismic velocity structure of the transition zone: Convection model predictions for a magnesium-endmember olivine-pyroxene mantle} by Michael H.G. Jacobs and Arie P. van den Berg, Physics of the Earth and Planetary Interiors, Volume 186, Issues 1-2, May 2011, Pages 36--48. See \url{http://www.sciencedirect.com/science/article/pii/S0031920111000422}.
-
`SolCx': A material model that corresponds to the 'SolCx' benchmark defined in Duretz et al., G-Cubed, 2011.
`SolKz': A material model that corresponds to the 'SolKz' benchmark defined in Duretz et al., G-Cubed, 2011.
`Inclusion': A material model that corresponds to the 'Inclusion' benchmark defined in Duretz et al., G-Cubed, 2011.
+`simple': A simple material model that has constant values for all coefficients but the density and viscosity. This model uses the formulation that assumes an incompressible medium despite the fact that the density follows the law $
ho(T)=
ho_0(1-eta(T-T_{ ext{ref}})$. The temperature dependency of viscosity is switched off by default and follows the formula$\eta(T)=\eta_0*e^{\eta_T*\Delta T / T_{ ext{ref}})}$.The value for the components of this formula and additional parameters are read from the parameter file in subsection 'Simple model'.
-{\it Possible values:} [Selection Steinberger|latent heat|simple|Tan Gurnis|table|SolCx|SolKz|Inclusion ]
+`Tan Gurnis': A simple compressible material model based on a benchmark from the paper of Tan/Gurnis (2007). This does not use the temperature equation, but has a hardcoded temperature.
+
+`Steinberger': lookup viscosity from the paper of Steinberger/Calderwood2006 and material data from a database generated by Perplex. The database builds upon the thermodynamic database by Stixrude 2011 and assumes a pyrolitic composition by Ringwood 1988.
+
+
+{\it Possible values:} [Selection latent heat|SolCx|SolKz|Inclusion|simple|Tan Gurnis|Steinberger ]
\end{itemize}
@@ -2641,481 +2639,6 @@
{\it Possible values:} [Anything]
\end{itemize}
-\subsection{Parameters in section t Material model/Table model}
-\label{parameters:Material_20model/Table_20model}
-
-egin{itemize}
-\item {\it Parameter name:} { t Composition}
-
-
-\index[prmindex]{Composition}
-\index[prmindexfull]{Material model!Table model!Composition}
-{\it Value:} standard
-
-
-{\it Default:} standard
-
-
-{\it Description:} The Composition of the model.
-
-
-{\it Possible values:} [Anything]
-\item {\it Parameter name:} { t Compressible}
-
-
-\index[prmindex]{Compressible}
-\index[prmindexfull]{Material model!Table model!Compressible}
-{\it Value:} true
-
-
-{\it Default:} true
-
-
-{\it Description:} whether the model is compressible.
-
-
-{\it Possible values:} [Bool]
-\item {\it Parameter name:} { t ComputePhases}
-
-
-\index[prmindex]{ComputePhases}
-\index[prmindexfull]{Material model!Table model!ComputePhases}
-{\it Value:} false
-
-
-{\it Default:} false
-
-
-{\it Description:} whether to compute phases.
-
-
-{\it Possible values:} [Bool]
-\item {\it Parameter name:} { t Gravity}
-
-
-\index[prmindex]{Gravity}
-\index[prmindexfull]{Material model!Table model!Gravity}
-{\it Value:} 30
-
-
-{\it Default:} 30
-
-
-{\it Description:} The value of the gravity constant.Units: $m/s^2$.
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Path to model data}
-
-
-\index[prmindex]{Path to model data}
-\index[prmindexfull]{Material model!Table model!Path to model data}
-{\it Value:} data/material-model/table/
-
-
-{\it Default:} data/material-model/table/
-
-
-{\it Description:} The path to the model data.
-
-
-{\it Possible values:} [DirectoryName]
-\item {\it Parameter name:} { t Reference density}
-
-
-\index[prmindex]{Reference density}
-\index[prmindexfull]{Material model!Table model!Reference density}
-{\it Value:} 3300
-
-
-{\it Default:} 3300
-
-
-{\it Description:} Reference density $
ho_0$. Units: $kg/m^3$.
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Reference specific heat}
-
-
-\index[prmindex]{Reference specific heat}
-\index[prmindexfull]{Material model!Table model!Reference specific heat}
-{\it Value:} 1250
-
-
-{\it Default:} 1250
-
-
-{\it Description:} The value of the specific heat $cp$. Units: $J/kg/K$.
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Reference temperature}
-
-
-\index[prmindex]{Reference temperature}
-\index[prmindexfull]{Material model!Table model!Reference temperature}
-{\it Value:} 293
-
-
-{\it Default:} 293
-
-
-{\it Description:} The reference temperature $T_0$. Units: $K$.
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Thermal conductivity}
-
-
-\index[prmindex]{Thermal conductivity}
-\index[prmindexfull]{Material model!Table model!Thermal conductivity}
-{\it Value:} 4.7
-
-
-{\it Default:} 4.7
-
-
-{\it Description:} The value of the thermal conductivity $k$. Units: $W/m/K$.
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Thermal expansion coefficient}
-
-
-\index[prmindex]{Thermal expansion coefficient}
-\index[prmindexfull]{Material model!Table model!Thermal expansion coefficient}
-{\it Value:} 2e-5
-
-
-{\it Default:} 2e-5
-
-
-{\it Description:} The value of the thermal expansion coefficient $eta$. Units: $1/K$.
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\end{itemize}
-
-
-
-\subsection{Parameters in section t Material model/Table model/Viscosity}
-\label{parameters:Material_20model/Table_20model/Viscosity}
-
-egin{itemize}
-\item {\it Parameter name:} { t Reference Viscosity}
-
-
-\index[prmindex]{Reference Viscosity}
-\index[prmindexfull]{Material model!Table model!Viscosity!Reference Viscosity}
-{\it Value:} 5e24
-
-
-{\it Default:} 5e24
-
-
-{\it Description:} The value of the constant viscosity. Units: $kg/m/s$.
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Viscosity Model}
-
-
-\index[prmindex]{Viscosity Model}
-\index[prmindexfull]{Material model!Table model!Viscosity!Viscosity Model}
-{\it Value:} Exponential
-
-
-{\it Default:} Exponential
-
-
-{\it Description:} Viscosity Model
-
-
-{\it Possible values:} [Anything]
-\item {\it Parameter name:} { t Viscosity increase lower mantle}
-
-
-\index[prmindex]{Viscosity increase lower mantle}
-\index[prmindexfull]{Material model!Table model!Viscosity!Viscosity increase lower mantle}
-{\it Value:} 1e0
-
-
-{\it Default:} 1e0
-
-
-{\it Description:} The Viscosity increase (jump) in the lower mantle.
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\end{itemize}
-
-
-
-\subsection{Parameters in section t Material model/Table model/Viscosity/Composite}
-\label{parameters:Material_20model/Table_20model/Viscosity/Composite}
-
-egin{itemize}
-\item {\it Parameter name:} { t Activation energy diffusion}
-
-
-\index[prmindex]{Activation energy diffusion}
-\index[prmindexfull]{Material model!Table model!Viscosity!Composite!Activation energy diffusion}
-{\it Value:} 335e3
-
-
-{\it Default:} 335e3
-
-
-{\it Description:} activation energy for diffusion creep
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Activation energy dislocation}
-
-
-\index[prmindex]{Activation energy dislocation}
-\index[prmindexfull]{Material model!Table model!Viscosity!Composite!Activation energy dislocation}
-{\it Value:} 540e3
-
-
-{\it Default:} 540e3
-
-
-{\it Description:} activation energy for dislocation creep
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Activation volume diffusion}
-
-
-\index[prmindex]{Activation volume diffusion}
-\index[prmindexfull]{Material model!Table model!Viscosity!Composite!Activation volume diffusion}
-{\it Value:} 4.0e-6
-
-
-{\it Default:} 4.0e-6
-
-
-{\it Description:} activation volume for diffusion creep
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Activation volume dislocation}
-
-
-\index[prmindex]{Activation volume dislocation}
-\index[prmindexfull]{Material model!Table model!Viscosity!Composite!Activation volume dislocation}
-{\it Value:} 14.0e-6
-
-
-{\it Default:} 14.0e-6
-
-
-{\it Description:} activation volume for dislocation creep
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Prefactor diffusion}
-
-
-\index[prmindex]{Prefactor diffusion}
-\index[prmindexfull]{Material model!Table model!Viscosity!Composite!Prefactor diffusion}
-{\it Value:} 1.92e-11
-
-
-{\it Default:} 1.92e-11
-
-
-{\it Description:} prefactor for diffusion creep (1e0/prefactor)*exp((activation\_energy+activation\_volume*pressure)/(R*temperature))
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Prefactor dislocation}
-
-
-\index[prmindex]{Prefactor dislocation}
-\index[prmindexfull]{Material model!Table model!Viscosity!Composite!Prefactor dislocation}
-{\it Value:} 2.42e-10
-
-
-{\it Default:} 2.42e-10
-
-
-{\it Description:} prefactor for dislocation creep (1e0/prefactor)*exp((activation\_energy+activation\_volume*pressure)/(R*temperature))
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Stress exponent}
-
-
-\index[prmindex]{Stress exponent}
-\index[prmindexfull]{Material model!Table model!Viscosity!Composite!Stress exponent}
-{\it Value:} 3.5
-
-
-{\it Default:} 3.5
-
-
-{\it Description:} stress exponent for dislocation creep
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\end{itemize}
-
-\subsection{Parameters in section t Material model/Table model/Viscosity/Diffusion}
-\label{parameters:Material_20model/Table_20model/Viscosity/Diffusion}
-
-egin{itemize}
-\item {\it Parameter name:} { t Activation energy diffusion}
-
-
-\index[prmindex]{Activation energy diffusion}
-\index[prmindexfull]{Material model!Table model!Viscosity!Diffusion!Activation energy diffusion}
-{\it Value:} 335e3
-
-
-{\it Default:} 335e3
-
-
-{\it Description:} activation energy for diffusion creep
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Activation volume diffusion}
-
-
-\index[prmindex]{Activation volume diffusion}
-\index[prmindexfull]{Material model!Table model!Viscosity!Diffusion!Activation volume diffusion}
-{\it Value:} 4.0e-6
-
-
-{\it Default:} 4.0e-6
-
-
-{\it Description:} activation volume for diffusion creep
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Prefactor diffusion}
-
-
-\index[prmindex]{Prefactor diffusion}
-\index[prmindexfull]{Material model!Table model!Viscosity!Diffusion!Prefactor diffusion}
-{\it Value:} 1.92e-11
-
-
-{\it Default:} 1.92e-11
-
-
-{\it Description:} prefactor for diffusion creep (1e0/prefactor)*exp((activation\_energy+activation\_volume*pressure)/(R*temperature))
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\end{itemize}
-
-\subsection{Parameters in section t Material model/Table model/Viscosity/Dislocation}
-\label{parameters:Material_20model/Table_20model/Viscosity/Dislocation}
-
-egin{itemize}
-\item {\it Parameter name:} { t Activation energy dislocation}
-
-
-\index[prmindex]{Activation energy dislocation}
-\index[prmindexfull]{Material model!Table model!Viscosity!Dislocation!Activation energy dislocation}
-{\it Value:} 335e3
-
-
-{\it Default:} 335e3
-
-
-{\it Description:} activation energy for dislocation creep
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Activation volume dislocation}
-
-
-\index[prmindex]{Activation volume dislocation}
-\index[prmindexfull]{Material model!Table model!Viscosity!Dislocation!Activation volume dislocation}
-{\it Value:} 4.0e-6
-
-
-{\it Default:} 4.0e-6
-
-
-{\it Description:} activation volume for dislocation creep
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Prefactor dislocation}
-
-
-\index[prmindex]{Prefactor dislocation}
-\index[prmindexfull]{Material model!Table model!Viscosity!Dislocation!Prefactor dislocation}
-{\it Value:} 1.92e-11
-
-
-{\it Default:} 1.92e-11
-
-
-{\it Description:} prefactor for dislocation creep (1e0/prefactor)*exp((activation\_energy+activation\_volume*pressure)/(R*temperature))
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Stress exponent}
-
-
-\index[prmindex]{Stress exponent}
-\index[prmindexfull]{Material model!Table model!Viscosity!Dislocation!Stress exponent}
-{\it Value:} 3.5
-
-
-{\it Default:} 3.5
-
-
-{\it Description:} stress exponent for dislocation creep
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\end{itemize}
-
-\subsection{Parameters in section t Material model/Table model/Viscosity/Exponential}
-\label{parameters:Material_20model/Table_20model/Viscosity/Exponential}
-
-egin{itemize}
-\item {\it Parameter name:} { t Exponential P}
-
-
-\index[prmindex]{Exponential P}
-\index[prmindexfull]{Material model!Table model!Viscosity!Exponential!Exponential P}
-{\it Value:} 1
-
-
-{\it Default:} 1
-
-
-{\it Description:} multiplication factor or Pressure exponent
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\item {\it Parameter name:} { t Exponential T}
-
-
-\index[prmindex]{Exponential T}
-\index[prmindexfull]{Material model!Table model!Viscosity!Exponential!Exponential T}
-{\it Value:} 1
-
-
-{\it Default:} 1
-
-
-{\it Description:} multiplication factor or Temperature exponent
-
-
-{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
-\end{itemize}
-
\subsection{Parameters in section t Material model/Tan Gurnis model}
\label{parameters:Material_20model/Tan_20Gurnis_20model}
@@ -3440,26 +2963,26 @@
The following criteria are available:
-`viscosity': A mesh refinement criterion that computes refinement indicators from a field that describes the spatial variability of the logarithm of the viscosity, $\log\eta$. (We choose the logarithm of the viscosity because it can vary by orders of magnitude.)Because this quantity may not be a continuous function ($\eta$ may be a discontinuous function along discontinuities in the medium, for example due to phase changes), we approximate the gradient of this quantity to refine the mesh. The error indicator defined here takes the magnitude of the approximate gradient and scales it by $h_K^{1+d/2}$ where $h_K$ is the diameter of each cell and $d$ is the dimension. This scaling ensures that the error indicators converge to zero as $h_K
ightarrow 0$ even if the energy density is discontinuous, since the gradient of a discontinuous function grows like $1/h_K$.
+`composition': A mesh refinement criterion that computes refinement indicators from the compositional fields. If there is more than one compositional field, then it simply takes the sum of the indicators computed from each of the compositional field.
-`velocity': A mesh refinement criterion that computes refinement indicators from the velocity field.
+`thermal energy density': A mesh refinement criterion that computes refinement indicators from a field that describes the spatial variability of the thermal energy density, $
ho C_p T$. Because this quantity may not be a continuous function ($
ho$ and $C_p$ may be discontinuous functions along discontinuities in the medium, for example due to phase changes), we approximate the gradient of this quantity to refine the mesh. The error indicator defined here takes the magnitude of the approximate gradient and scales it by $h_K^{1.5}$ where $h_K$ is the diameter of each cell. This scaling ensures that the error indicators converge to zero as $h_K
ightarrow 0$ even if the energy density is discontinuous, since the gradient of a discontinuous function grows like $1/h_K$.
-`temperature': A mesh refinement criterion that computes refinement indicators from the temperature field.
+`nonadiabatic temperature': A mesh refinement criterion that computes refinement indicators from the excess temperature(difference between temperature and adiabatic temperature.
-`density': A mesh refinement criterion that computes refinement indicators from a field that describes the spatial variability of the density, $
ho$. Because this quantity may not be a continuous function ($
ho$ and $C_p$ may be discontinuous functions along discontinuities in the medium, for example due to phase changes), we approximate the gradient of this quantity to refine the mesh. The error indicator defined here takes the magnitude of the approximate gradient and scales it by $h_K^{1+d/2}$ where $h_K$ is the diameter of each cell and $d$ is the dimension. This scaling ensures that the error indicators converge to zero as $h_K
ightarrow 0$ even if the energy density is discontinuous, since the gradient of a discontinuous function grows like $1/h_K$.
-
`topography': A class that implements a mesh refinement criterion, which always flags all cells in the uppermost layer for refinement. This is useful to provide high accuracy for processes at or close to the surface.
To use this refinement criterion, you may want to combine it with other refinement criteria, setting the 'Normalize individual refinement criteria' flag and using the 'max' setting for 'Refinement criteria merge operation'.
-`composition': A mesh refinement criterion that computes refinement indicators from the compositional fields. If there is more than one compositional field, then it simply takes the sum of the indicators computed from each of the compositional field.
+`velocity': A mesh refinement criterion that computes refinement indicators from the velocity field.
-`thermal energy density': A mesh refinement criterion that computes refinement indicators from a field that describes the spatial variability of the thermal energy density, $
ho C_p T$. Because this quantity may not be a continuous function ($
ho$ and $C_p$ may be discontinuous functions along discontinuities in the medium, for example due to phase changes), we approximate the gradient of this quantity to refine the mesh. The error indicator defined here takes the magnitude of the approximate gradient and scales it by $h_K^{1.5}$ where $h_K$ is the diameter of each cell. This scaling ensures that the error indicators converge to zero as $h_K
ightarrow 0$ even if the energy density is discontinuous, since the gradient of a discontinuous function grows like $1/h_K$.
+`density': A mesh refinement criterion that computes refinement indicators from a field that describes the spatial variability of the density, $
ho$. Because this quantity may not be a continuous function ($
ho$ and $C_p$ may be discontinuous functions along discontinuities in the medium, for example due to phase changes), we approximate the gradient of this quantity to refine the mesh. The error indicator defined here takes the magnitude of the approximate gradient and scales it by $h_K^{1+d/2}$ where $h_K$ is the diameter of each cell and $d$ is the dimension. This scaling ensures that the error indicators converge to zero as $h_K
ightarrow 0$ even if the energy density is discontinuous, since the gradient of a discontinuous function grows like $1/h_K$.
-`nonadiabatic temperature': A mesh refinement criterion that computes refinement indicators from the excess temperature(difference between temperature and adiabatic temperature.
+`viscosity': A mesh refinement criterion that computes refinement indicators from a field that describes the spatial variability of the logarithm of the viscosity, $\log\eta$. (We choose the logarithm of the viscosity because it can vary by orders of magnitude.)Because this quantity may not be a continuous function ($\eta$ may be a discontinuous function along discontinuities in the medium, for example due to phase changes), we approximate the gradient of this quantity to refine the mesh. The error indicator defined here takes the magnitude of the approximate gradient and scales it by $h_K^{1+d/2}$ where $h_K$ is the diameter of each cell and $d$ is the dimension. This scaling ensures that the error indicators converge to zero as $h_K
ightarrow 0$ even if the energy density is discontinuous, since the gradient of a discontinuous function grows like $1/h_K$.
+`temperature': A mesh refinement criterion that computes refinement indicators from the temperature field.
-{\it Possible values:} [MultipleSelection viscosity|velocity|temperature|density|topography|composition|thermal energy density|nonadiabatic temperature ]
+
+{\it Possible values:} [MultipleSelection composition|thermal energy density|nonadiabatic temperature|topography|velocity|density|viscosity|temperature ]
\item {\it Parameter name:} { t Time steps between mesh refinement}
@@ -3687,29 +3210,25 @@
The following postprocessors are available:
-`velocity statistics for the table model': A postprocessor that computes some statistics about the velocity field.
-
-`heat flux statistics for the table model': A postprocessor that computes some statistics about the heat flux across boundaries.
-
-`visualization': A postprocessor that takes the solution and writes it into files that can be read by a graphical visualization program. Additional run time parameters are read from the parameter subsection 'Visualization'.
-
`temperature statistics': A postprocessor that computes some statistics about the temperature field.
-`composition statistics': A postprocessor that computes some statistics about the compositional fields, if present in this simulation. In particular, it computes maximal and minimal values of each field, as well as the total mass contained in this field as defined by the integral $m_i(t) = \int_\Omega c_i(\mathbf x,t) \; dx$.
-
`heat flux statistics': A postprocessor that computes some statistics about the (conductive) heat flux across boundaries. For each boundary indicator (see your geometry description for which boundary indicators are used), the heat flux is computed in outward direction, i.e., from the domain to the outside, using the formula $\int_{\Gamma_i} k
abla T
More information about the CIG-COMMITS
mailing list