[cig-commits] commit 2516 by heister to /var/svn/dealii/aspect
dealii.demon at gmail.com
dealii.demon at gmail.com
Mon Apr 14 07:55:04 PDT 2014
Revision 2516
merge from mainline
_U branches/freesurface/
U branches/freesurface/doc/manual/parameters.tex
U branches/freesurface/doc/modules/0.1-to-0.2.h
U branches/freesurface/doc/modules/0.2-to-0.3.h
U branches/freesurface/doc/modules/changes.h
U branches/freesurface/doc/modules/groups.h
U branches/freesurface/doc/modules/main.h
U branches/freesurface/doc/modules/todo.h
U branches/freesurface/doc/release-tasklist
U branches/freesurface/doc/update_parameters.sh
U branches/freesurface/include/aspect/adiabatic_conditions.h
U branches/freesurface/include/aspect/boundary_composition/box.h
U branches/freesurface/include/aspect/boundary_composition/initial_composition.h
U branches/freesurface/include/aspect/boundary_composition/interface.h
U branches/freesurface/include/aspect/boundary_composition/spherical_constant.h
U branches/freesurface/include/aspect/boundary_temperature/box.h
U branches/freesurface/include/aspect/boundary_temperature/constant.h
U branches/freesurface/include/aspect/boundary_temperature/initial_temperature.h
U branches/freesurface/include/aspect/boundary_temperature/interface.h
U branches/freesurface/include/aspect/boundary_temperature/spherical_constant.h
U branches/freesurface/include/aspect/boundary_temperature/tan_gurnis.h
U branches/freesurface/include/aspect/compositional_initial_conditions/function.h
U branches/freesurface/include/aspect/compositional_initial_conditions/interface.h
U branches/freesurface/include/aspect/geometry_model/box.h
U branches/freesurface/include/aspect/geometry_model/interface.h
U branches/freesurface/include/aspect/geometry_model/sphere.h
U branches/freesurface/include/aspect/geometry_model/spherical_shell.h
U branches/freesurface/include/aspect/global.h
U branches/freesurface/include/aspect/gravity_model/interface.h
U branches/freesurface/include/aspect/gravity_model/radial.h
U branches/freesurface/include/aspect/gravity_model/vertical.h
U branches/freesurface/include/aspect/initial_conditions/adiabatic.h
U branches/freesurface/include/aspect/initial_conditions/box.h
U branches/freesurface/include/aspect/initial_conditions/function.h
U branches/freesurface/include/aspect/initial_conditions/harmonic_perturbation.h
U branches/freesurface/include/aspect/initial_conditions/interface.h
U branches/freesurface/include/aspect/initial_conditions/spherical_shell.h
U branches/freesurface/include/aspect/introspection.h
U branches/freesurface/include/aspect/material_model/duretz_et_al.h
U branches/freesurface/include/aspect/material_model/interface.h
U branches/freesurface/include/aspect/material_model/latent_heat.h
U branches/freesurface/include/aspect/material_model/simple.h
U branches/freesurface/include/aspect/material_model/steinberger.h
U branches/freesurface/include/aspect/material_model/table.h
U branches/freesurface/include/aspect/material_model/tan_gurnis.h
U branches/freesurface/include/aspect/mesh_refinement/composition.h
U branches/freesurface/include/aspect/mesh_refinement/density.h
U branches/freesurface/include/aspect/mesh_refinement/interface.h
U branches/freesurface/include/aspect/mesh_refinement/nonadiabatic_temperature.h
U branches/freesurface/include/aspect/mesh_refinement/temperature.h
U branches/freesurface/include/aspect/mesh_refinement/thermal_energy_density.h
U branches/freesurface/include/aspect/mesh_refinement/topography.h
U branches/freesurface/include/aspect/mesh_refinement/velocity.h
U branches/freesurface/include/aspect/mesh_refinement/viscosity.h
U branches/freesurface/include/aspect/particle/generator.h
U branches/freesurface/include/aspect/particle/integrator.h
U branches/freesurface/include/aspect/particle/output.h
U branches/freesurface/include/aspect/particle/particle.h
U branches/freesurface/include/aspect/particle/world.h
U branches/freesurface/include/aspect/plugins.h
U branches/freesurface/include/aspect/postprocess/composition_statistics.h
U branches/freesurface/include/aspect/postprocess/depth_average.h
U branches/freesurface/include/aspect/postprocess/duretz_et_al.h
U branches/freesurface/include/aspect/postprocess/dynamic_topography.h
U branches/freesurface/include/aspect/postprocess/heat_flux_statistics.h
U branches/freesurface/include/aspect/postprocess/interface.h
U branches/freesurface/include/aspect/postprocess/table_heat_flux_statistics.h
U branches/freesurface/include/aspect/postprocess/table_velocity_statistics.h
U branches/freesurface/include/aspect/postprocess/tan_gurnis.h
U branches/freesurface/include/aspect/postprocess/temperature_statistics.h
U branches/freesurface/include/aspect/postprocess/tracer.h
U branches/freesurface/include/aspect/postprocess/velocity_statistics.h
U branches/freesurface/include/aspect/postprocess/visualization/artificial_viscosity.h
U branches/freesurface/include/aspect/postprocess/visualization/density.h
U branches/freesurface/include/aspect/postprocess/visualization/dynamic_topography.h
U branches/freesurface/include/aspect/postprocess/visualization/error_indicator.h
U branches/freesurface/include/aspect/postprocess/visualization/friction_heating.h
U branches/freesurface/include/aspect/postprocess/visualization/melt_fraction.h
U branches/freesurface/include/aspect/postprocess/visualization/nonadiabatic_pressure.h
U branches/freesurface/include/aspect/postprocess/visualization/nonadiabatic_temperature.h
U branches/freesurface/include/aspect/postprocess/visualization/partition.h
U branches/freesurface/include/aspect/postprocess/visualization/seismic_anomalies.h
U branches/freesurface/include/aspect/postprocess/visualization/seismic_vp.h
U branches/freesurface/include/aspect/postprocess/visualization/seismic_vs.h
U branches/freesurface/include/aspect/postprocess/visualization/specific_heat.h
U branches/freesurface/include/aspect/postprocess/visualization/strain_rate.h
U branches/freesurface/include/aspect/postprocess/visualization/thermal_expansivity.h
U branches/freesurface/include/aspect/postprocess/visualization/thermodynamic_phase.h
U branches/freesurface/include/aspect/postprocess/visualization/viscosity.h
U branches/freesurface/include/aspect/postprocess/visualization/viscosity_ratio.h
U branches/freesurface/include/aspect/postprocess/visualization.h
U branches/freesurface/include/aspect/simulator.h
U branches/freesurface/include/aspect/simulator_access.h
U branches/freesurface/include/aspect/termination_criteria/end_step.h
U branches/freesurface/include/aspect/termination_criteria/end_time.h
U branches/freesurface/include/aspect/termination_criteria/interface.h
U branches/freesurface/include/aspect/termination_criteria/steady_rms_velocity.h
U branches/freesurface/include/aspect/termination_criteria/user_request.h
U branches/freesurface/include/aspect/velocity_boundary_conditions/function.h
U branches/freesurface/include/aspect/velocity_boundary_conditions/gplates.h
U branches/freesurface/include/aspect/velocity_boundary_conditions/interface.h
U branches/freesurface/source/boundary_composition/initial_composition.cc
U branches/freesurface/source/boundary_composition/interface.cc
U branches/freesurface/source/boundary_composition/spherical_constant.cc
U branches/freesurface/source/boundary_temperature/constant.cc
U branches/freesurface/source/boundary_temperature/initial_temperature.cc
U branches/freesurface/source/boundary_temperature/interface.cc
U branches/freesurface/source/geometry_model/box.cc
U branches/freesurface/source/geometry_model/interface.cc
U branches/freesurface/source/geometry_model/sphere.cc
U branches/freesurface/source/gravity_model/interface.cc
U branches/freesurface/source/gravity_model/radial.cc
U branches/freesurface/source/initial_conditions/adiabatic.cc
U branches/freesurface/source/initial_conditions/harmonic_perturbation.cc
U branches/freesurface/source/initial_conditions/interface.cc
U branches/freesurface/source/main.cc
U branches/freesurface/source/material_model/interface.cc
U branches/freesurface/source/material_model/latent_heat.cc
U branches/freesurface/source/material_model/simple.cc
U branches/freesurface/source/material_model/steinberger.cc
U branches/freesurface/source/mesh_refinement/composition.cc
U branches/freesurface/source/mesh_refinement/density.cc
U branches/freesurface/source/mesh_refinement/nonadiabatic_temperature.cc
U branches/freesurface/source/mesh_refinement/topography.cc
U branches/freesurface/source/mesh_refinement/velocity.cc
U branches/freesurface/source/mesh_refinement/viscosity.cc
U branches/freesurface/source/postprocess/composition_statistics.cc
U branches/freesurface/source/postprocess/depth_average.cc
U branches/freesurface/source/postprocess/dynamic_topography.cc
U branches/freesurface/source/postprocess/heat_flux_statistics.cc
U branches/freesurface/source/postprocess/table_heat_flux_statistics.cc
U branches/freesurface/source/postprocess/temperature_statistics.cc
U branches/freesurface/source/postprocess/visualization/density.cc
U branches/freesurface/source/postprocess/visualization/dynamic_topography.cc
U branches/freesurface/source/postprocess/visualization/friction_heating.cc
U branches/freesurface/source/postprocess/visualization/melt_fraction.cc
U branches/freesurface/source/postprocess/visualization/nonadiabatic_pressure.cc
U branches/freesurface/source/postprocess/visualization/nonadiabatic_temperature.cc
U branches/freesurface/source/postprocess/visualization/partition.cc
U branches/freesurface/source/postprocess/visualization/seismic_vp.cc
U branches/freesurface/source/postprocess/visualization/seismic_vs.cc
U branches/freesurface/source/postprocess/visualization/specific_heat.cc
U branches/freesurface/source/postprocess/visualization/strain_rate.cc
U branches/freesurface/source/postprocess/visualization/thermal_expansivity.cc
U branches/freesurface/source/postprocess/visualization/thermodynamic_phase.cc
U branches/freesurface/source/postprocess/visualization/viscosity.cc
U branches/freesurface/source/postprocess/visualization/viscosity_ratio.cc
U branches/freesurface/source/postprocess/visualization.cc
U branches/freesurface/source/simulator/assembly.cc
U branches/freesurface/source/simulator/checkpoint_restart.cc
U branches/freesurface/source/simulator/core.cc
U branches/freesurface/source/simulator/helper_functions.cc
U branches/freesurface/source/simulator/initial_conditions.cc
U branches/freesurface/source/simulator/nullspace.cc
U branches/freesurface/source/simulator/parameters.cc
U branches/freesurface/source/simulator/solver.cc
U branches/freesurface/source/termination_criteria/end_time.cc
U branches/freesurface/source/termination_criteria/interface.cc
U branches/freesurface/source/velocity_boundary_conditions/gplates.cc
A branches/freesurface/tests/maxtimestep/
D branches/freesurface/tests/maxtimestep/screen-output
A branches/freesurface/tests/maxtimestep/screen-output
D branches/freesurface/tests/maxtimestep/statistics
A branches/freesurface/tests/maxtimestep/statistics
A branches/freesurface/tests/maxtimestep.prm
A branches/freesurface/tests/refine_vel/
D branches/freesurface/tests/refine_vel/screen-output
A branches/freesurface/tests/refine_vel/screen-output
D branches/freesurface/tests/refine_vel/statistics
A branches/freesurface/tests/refine_vel/statistics
A branches/freesurface/tests/refine_vel.prm
http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=2516&peg=2516
Diff:
Index: branches/freesurface
===================================================================
--- branches/freesurface 2014-04-14 14:54:56 UTC (rev 2515)
+++ branches/freesurface 2014-04-14 14:55:02 UTC (rev 2516)
Property changes on: branches/freesurface
___________________________________________________________________
Modified: svn:mergeinfo
## -2,4 +2,4 ##
/branches/compositional:1141-1251
/branches/fully-nonlinear:542-728
/branches/j-dannberg:1679,1875-1882
-/trunk/aspect:2386-2478
+/trunk/aspect:2386-2513
\ No newline at end of property
Modified: branches/freesurface/doc/manual/parameters.tex
===================================================================
--- branches/freesurface/doc/manual/parameters.tex 2014-04-14 14:54:56 UTC (rev 2515)
+++ branches/freesurface/doc/manual/parameters.tex 2014-04-14 14:55:02 UTC (rev 2516)
@@ -16,7 +16,7 @@
{\it Description:} A list of names of additional shared libraries that should be loaded upon starting up the program. The names of these files can contain absolute or relative paths (relative to the directory in which you call ASPECT). In fact, file names that are do not contain any directory information (i.e., only the name of a file such as <myplugin.so> will not be found if they are not located in one of the directories listed in the LD\_LIBRARY\_PATH environment variable. In order to load a library in the current directory, use <./myplugin.so> instead.
-The typical use of this parameter is to so that you can implement additional plugins in your own directories, rather than in the ASPECT source directories. You can then simply compile these plugins into a shared library without having to re-compile all of ASPECT. See the section of the manual discussing writing extensions for more information on how to compile additional files into a shared library.
+The typical use of this parameter is so that you can implement additional plugins in your own directories, rather than in the ASPECT source directories. You can then simply compile these plugins into a shared library without having to re-compile all of ASPECT. See the section of the manual discussing writing extensions for more information on how to compile additional files into a shared library.
{\it Possible values:} [List list of <[FileName (Type: input)]> of length 0...4294967295 (inclusive)]
@@ -32,6 +32,7 @@
{\it Description:} In order to make the problem in the first time step easier to solve, we need a reasonable guess for the temperature and pressure. To obtain it, we use an adiabatic pressure and temperature field. This parameter describes what the `adiabatic' temperature would be at the surface of the domain (i.e. at depth zero). Note that this value need not coincide with the boundary condition posed at this point. Rather, the boundary condition may differ significantly from the adiabatic value, and then typically induce a thermal boundary layer.
+
For more information, see the section in the manual that discusses the general mathematical model.
@@ -86,13 +87,13 @@
\index[prmindex]{End time}
\index[prmindexfull]{End time}
-{\it Value:} 1e300
+{\it Value:} 5.6966627666142017e+300
-{\it Default:} 1e300
+{\it Default:} 5.6966627666142017e+300
-{\it Description:} The end time of the simulation. Units: years if the 'Use years in output instead of seconds' parameter is set; seconds otherwise.
+{\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.
{\it Possible values:} [Double -1.79769e+308...1.79769e+308 (inclusive)]
@@ -128,6 +129,21 @@
{\it Possible values:} [Integer range 0...2147483647 (inclusive)]
+\item {\it Parameter name:} { t Maximum time step}
+
+
+\index[prmindex]{Maximum time step}
+\index[prmindexfull]{Maximum time step}
+{\it Value:} 5.6966627666142017e+300
+
+
+{\it Default:} 5.6966627666142017e+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.
+
+
+{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
\item {\it Parameter name:} { t Nonlinear solver scheme}
@@ -143,6 +159,21 @@
{\it Possible values:} [Selection IMPES|iterated IMPES|iterated Stokes|Stokes only ]
+\item {\it Parameter name:} { t Nonlinear solver tolerance}
+
+
+\index[prmindex]{Nonlinear solver tolerance}
+\index[prmindexfull]{Nonlinear solver tolerance}
+{\it Value:} 1e-5
+
+
+{\it Default:} 1e-5
+
+
+{\it Description:} A relative tolerance up to which the nonlinear solver will iterate. This parameter is only relevant if Nonlinear solver scheme is set to 'iterated Stokes' or 'iterated IMPES'.
+
+
+{\it Possible values:} [Double 0...1 (inclusive)]
\item {\it Parameter name:} { t Number of cheap Stokes solver steps}
@@ -214,7 +245,7 @@
{\it Default:} 0
-{\it Description:} The start time of the simulation. Units: years if the 'Use years in output instead of seconds' parameter is set; seconds otherwise.
+{\it Description:} The start time of the simulation. Units: Years if the 'Use years in output instead of seconds' parameter is set; seconds otherwise.
{\it Possible values:} [Double -1.79769e+308...1.79769e+308 (inclusive)]
@@ -230,6 +261,7 @@
{\it Description:} The mathematical equations that describe thermal convection only determine the pressure up to an arbitrary constant. On the other hand, for comparison and for looking up material parameters it is important that the pressure be normalized somehow. We do this by enforcing a particular average pressure value at the surface of the domain, where the geometry model determines where the surface is. This parameter describes what this average surface pressure value is supposed to be. By default, it is set to zero, but one may want to choose a different value for example for simulating only the volume of the mantle below the lithosphere, in which case the surface pressure should be the lithostatic pressure at the bottom of the lithosphere.
+
For more information, see the section in the manual that discusses the general mathematical model.
@@ -298,6 +330,175 @@
+\subsection{Parameters in section t Boundary composition model}
+\label{parameters:Boundary_20composition_20model}
+
+egin{itemize}
+\item {\it Parameter name:} { t Model name}
+
+
+\index[prmindex]{Model name}
+\index[prmindexfull]{Boundary composition model!Model name}
+{\it Value:}
+
+
+{\it Default:}
+
+
+{\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''.
+
+
+{\it Possible values:} [Selection box|spherical constant|initial composition ]
+\end{itemize}
+
+
+
+\subsection{Parameters in section t Boundary composition model/Box}
+\label{parameters:Boundary_20composition_20model/Box}
+
+egin{itemize}
+\item {\it Parameter name:} { t Bottom composition}
+
+
+\index[prmindex]{Bottom composition}
+\index[prmindexfull]{Boundary composition model!Box!Bottom composition}
+{\it Value:}
+
+
+{\it Default:}
+
+
+{\it Description:} A comma separated list of composition boundary values at the bottom boundary (at minimal y-value in 2d, or minimal z-value in 3d). This list must have as many entries as there are compositional fields. Units: none.
+
+
+{\it Possible values:} [List list of <[Double -1.79769e+308...1.79769e+308 (inclusive)]> of length 0...4294967295 (inclusive)]
+\item {\it Parameter name:} { t Left composition}
+
+
+\index[prmindex]{Left composition}
+\index[prmindexfull]{Boundary composition model!Box!Left composition}
+{\it Value:}
+
+
+{\it Default:}
+
+
+{\it Description:} A comma separated list of composition boundary values at the left boundary (at minimal x-value). This list must have as many entries as there are compositional fields. Units: none.
+
+
+{\it Possible values:} [List list of <[Double -1.79769e+308...1.79769e+308 (inclusive)]> of length 0...4294967295 (inclusive)]
+\item {\it Parameter name:} { t Right composition}
+
+
+\index[prmindex]{Right composition}
+\index[prmindexfull]{Boundary composition model!Box!Right composition}
+{\it Value:}
+
+
+{\it Default:}
+
+
+{\it Description:} A comma separated list of composition boundary values at the right boundary (at maximal x-value). This list must have as many entries as there are compositional fields. Units: none.
+
+
+{\it Possible values:} [List list of <[Double -1.79769e+308...1.79769e+308 (inclusive)]> of length 0...4294967295 (inclusive)]
+\item {\it Parameter name:} { t Top composition}
+
+
+\index[prmindex]{Top composition}
+\index[prmindexfull]{Boundary composition model!Box!Top composition}
+{\it Value:}
+
+
+{\it Default:}
+
+
+{\it Description:} A comma separated list of composition boundary values at the top boundary (at maximal y-value in 2d, or maximal z-value in 3d). This list must have as many entries as there are compositional fields. Units: none.
+
+
+{\it Possible values:} [List list of <[Double -1.79769e+308...1.79769e+308 (inclusive)]> of length 0...4294967295 (inclusive)]
+\end{itemize}
+
+\subsection{Parameters in section t Boundary composition model/Initial composition}
+\label{parameters:Boundary_20composition_20model/Initial_20composition}
+
+egin{itemize}
+\item {\it Parameter name:} { t Maximal composition}
+
+
+\index[prmindex]{Maximal composition}
+\index[prmindexfull]{Boundary composition model!Initial composition!Maximal composition}
+{\it Value:} 1
+
+
+{\it Default:} 1
+
+
+{\it Description:} Maximal composition. Units: none.
+
+
+{\it Possible values:} [Double -1.79769e+308...1.79769e+308 (inclusive)]
+\item {\it Parameter name:} { t Minimal composition}
+
+
+\index[prmindex]{Minimal composition}
+\index[prmindexfull]{Boundary composition model!Initial composition!Minimal composition}
+{\it Value:} 0
+
+
+{\it Default:} 0
+
+
+{\it Description:} Minimal composition. Units: none.
+
+
+{\it Possible values:} [Double -1.79769e+308...1.79769e+308 (inclusive)]
+\end{itemize}
+
+\subsection{Parameters in section t Boundary composition model/Spherical constant}
+\label{parameters:Boundary_20composition_20model/Spherical_20constant}
+
+egin{itemize}
+\item {\it Parameter name:} { t Inner composition}
+
+
+\index[prmindex]{Inner composition}
+\index[prmindexfull]{Boundary composition model!Spherical constant!Inner composition}
+{\it Value:} 1
+
+
+{\it Default:} 1
+
+
+{\it Description:} Composition at the inner boundary (core mantle boundary). Units: K.
+
+
+{\it Possible values:} [Double -1.79769e+308...1.79769e+308 (inclusive)]
+\item {\it Parameter name:} { t Outer composition}
+
+
+\index[prmindex]{Outer composition}
+\index[prmindexfull]{Boundary composition model!Spherical constant!Outer composition}
+{\it Value:} 0
+
+
+{\it Default:} 0
+
+
+{\it Description:} Composition at the outer boundary (lithosphere water/air). Units: K.
+
+
+{\it Possible values:} [Double -1.79769e+308...1.79769e+308 (inclusive)]
+\end{itemize}
+
\subsection{Parameters in section t Boundary temperature model}
\label{parameters:Boundary_20temperature_20model}
@@ -315,14 +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.
-`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 'Sherical 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'.
-{\it Possible values:} [Selection box|Tan Gurnis|spherical constant ]
+{\it Possible values:} [Selection constant|box|initial temperature|Tan Gurnis|spherical constant ]
\end{itemize}
@@ -393,6 +600,63 @@
{\it Possible values:} [Double -1.79769e+308...1.79769e+308 (inclusive)]
\end{itemize}
+\subsection{Parameters in section t Boundary temperature model/Constant}
+\label{parameters:Boundary_20temperature_20model/Constant}
+
+egin{itemize}
+\item {\it Parameter name:} { t Boundary indicator to temperature mappings}
+
+
+\index[prmindex]{Boundary indicator to temperature mappings}
+\index[prmindexfull]{Boundary temperature model!Constant!Boundary indicator to temperature mappings}
+{\it Value:}
+
+
+{\it Default:}
+
+
+{\it Description:} A comma separated list of mappings between boundary indicators and the temperature associated with the boundary indicators. The format for this list is ``indicator1 : value1, indicator2 : value2, ...'', where each indicator is a valid boundary indicator and each value is the temperature of that boundary.
+
+
+{\it Possible values:} [Map map of <[Integer range 0...255 (inclusive)]:[Double -1.79769e+308...1.79769e+308 (inclusive)]> of length 0...4294967295 (inclusive)]
+\end{itemize}
+
+\subsection{Parameters in section t Boundary temperature model/Initial temperature}
+\label{parameters:Boundary_20temperature_20model/Initial_20temperature}
+
+egin{itemize}
+\item {\it Parameter name:} { t Maximal temperature}
+
+
+\index[prmindex]{Maximal temperature}
+\index[prmindexfull]{Boundary temperature model!Initial temperature!Maximal temperature}
+{\it Value:} 3773
+
+
+{\it Default:} 3773
+
+
+{\it Description:} Maximal temperature. Units: K.
+
+
+{\it Possible values:} [Double -1.79769e+308...1.79769e+308 (inclusive)]
+\item {\it Parameter name:} { t Minimal temperature}
+
+
+\index[prmindex]{Minimal temperature}
+\index[prmindexfull]{Boundary temperature model!Initial temperature!Minimal temperature}
+{\it Value:} 0
+
+
+{\it Default:} 0
+
+
+{\it Description:} Minimal temperature. Units: K.
+
+
+{\it Possible values:} [Double -1.79769e+308...1.79769e+308 (inclusive)]
+\end{itemize}
+
\subsection{Parameters in section t Boundary temperature model/Spherical constant}
\label{parameters:Boundary_20temperature_20model/Spherical_20constant}
@@ -497,16 +761,31 @@
\index[prmindex]{Data directory}
\index[prmindexfull]{Boundary velocity model!GPlates model!Data directory}
-{\it Value:} data/velocity-boundary-conditions/gplates/
+{\it Value:} \$ASPECT\_SOURCE\_DIR/data/velocity-boundary-conditions/gplates/
-{\it Default:} data/velocity-boundary-conditions/gplates/
+{\it Default:} \$ASPECT\_SOURCE\_DIR/data/velocity-boundary-conditions/gplates/
-{\it Description:} The path to the model data.
+{\it Description:} The name of a directory that contains the model data. This path may either be absolute (if starting with a '/') or relative to the current directory. The path may also include the special text '\$ASPECT\_SOURCE\_DIR' which will be interpreted as the path in which the ASPECT source files were located when ASPECT was compiled. This interpretation allows, for example, to reference files located in the 'data/' subdirectory of ASPECT.
{\it Possible values:} [DirectoryName]
+\item {\it Parameter name:} { t Interpolation width}
+
+
+\index[prmindex]{Interpolation width}
+\index[prmindexfull]{Boundary velocity model!GPlates model!Interpolation width}
+{\it Value:} 0
+
+
+{\it Default:} 0
+
+
+{\it Description:} Determines the width of the velocity interpolation zone around the current point. Currently equals the arc distance between evaluation point and velocity data point that is still included in the interpolation. The weighting of the points currently only accounts for the surface area a single data point is covering ('moving window' interpolation without distance weighting).
+
+
+{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
\item {\it Parameter name:} { t Point one}
@@ -542,16 +821,16 @@
\index[prmindex]{Time step}
\index[prmindexfull]{Boundary velocity model!GPlates model!Time step}
-{\it Value:} 3.1558e13
+{\it Value:} 1e6
-{\it Default:} 3.1558e13
+{\it Default:} 1e6
-{\it Description:} Time step between following velocity files. Default is one million years expressed in SI units.
+{\it Description:} Time step between following velocity files. Depending on the setting of the global 'Use years in output instead of seconds' flag in the input file, this number is either interpreted as seconds or as years. The default is one million, i.e., either one million seconds or one million years.
-{\it Possible values:} [Anything]
+{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
\item {\it Parameter name:} { t Velocity file name}
@@ -578,10 +857,10 @@
{\it Default:} 0.0
-{\it Description:} Time at which the velocity file with number 0 shall be loaded.Previous to this time, a no-slip boundary condition is assumed.
+{\it Description:} Time at which the velocity file with number 0 shall be loaded. Previous to this time, a no-slip boundary condition is assumed. Depending on the setting of the global 'Use years in output instead of seconds' flag in the input file, this number is either interpreted as seconds or as years.
-{\it Possible values:} [Anything]
+{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
\end{itemize}
\subsection{Parameters in section t Checkpointing}
@@ -882,8 +1161,12 @@
`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$.
-{\it Possible values:} [Selection box|spherical shell ]
+`sphere': Geometry model for sphere with a user specified radius.
+
+
+{\it Possible values:} [Selection box|spherical shell|sphere ]
\end{itemize}
@@ -922,6 +1205,21 @@
{\it Possible values:} [Bool]
+\item {\it Parameter name:} { t X repetitions}
+
+
+\index[prmindex]{X repetitions}
+\index[prmindexfull]{Geometry model!Box!X repetitions}
+{\it Value:} 1
+
+
+{\it Default:} 1
+
+
+{\it Description:} Number of cells in X direction.
+
+
+{\it Possible values:} [Integer range 1...2147483647 (inclusive)]
\item {\it Parameter name:} { t Y extent}
@@ -952,6 +1250,21 @@
{\it Possible values:} [Bool]
+\item {\it Parameter name:} { t Y repetitions}
+
+
+\index[prmindex]{Y repetitions}
+\index[prmindexfull]{Geometry model!Box!Y repetitions}
+{\it Value:} 1
+
+
+{\it Default:} 1
+
+
+{\it Description:} Number of cells in Y direction.
+
+
+{\it Possible values:} [Integer range 1...2147483647 (inclusive)]
\item {\it Parameter name:} { t Z extent}
@@ -982,8 +1295,44 @@
{\it Possible values:} [Bool]
+\item {\it Parameter name:} { t Z repetitions}
+
+
+\index[prmindex]{Z repetitions}
+\index[prmindexfull]{Geometry model!Box!Z repetitions}
+{\it Value:} 1
+
+
+{\it Default:} 1
+
+
+{\it Description:} Number of cells in Z direction.
+
+
+{\it Possible values:} [Integer range 1...2147483647 (inclusive)]
\end{itemize}
+\subsection{Parameters in section t Geometry model/Sphere}
+\label{parameters:Geometry_20model/Sphere}
+
+egin{itemize}
+\item {\it Parameter name:} { t Radius}
+
+
+\index[prmindex]{Radius}
+\index[prmindexfull]{Geometry model!Sphere!Radius}
+{\it Value:} 6371000
+
+
+{\it Default:} 6371000
+
+
+{\it Description:} Radius of the sphere. Units: m.
+
+
+{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
+\end{itemize}
+
\subsection{Parameters in section t Geometry model/Spherical shell}
\label{parameters:Geometry_20model/Spherical_20shell}
@@ -1056,10 +1405,12 @@
`radial earth-like': A gravity model in which the gravity direction is radially inward and with a magnitude that matches that of the earth at the core-mantle boundary as well as at the surface and in between is physically correct under the assumption of a constant density.
+`radial linear': A gravity model which is radially inward, where the magnitudedecreases linearly with depth, as you would get with a constantdensity spherical domain.
+
`vertical': A gravity model in which the gravity direction is vertically downward and at a constant magnitude by default equal to one.
-{\it Possible values:} [Selection radial constant|radial earth-like|vertical ]
+{\it Possible values:} [Selection radial constant|radial earth-like|radial linear|vertical ]
\end{itemize}
@@ -1073,18 +1424,39 @@
\index[prmindex]{Magnitude}
\index[prmindexfull]{Gravity model!Radial constant!Magnitude}
-{\it Value:} 30
+{\it Value:} 9.81
-{\it Default:} 30
+{\it Default:} 9.81
-{\it Description:} Magnitude of the gravity vector in $m/s^2$. The direction is always radially outward from the center of the earth.
+{\it Description:} Magnitude of the gravity vector in $m/s^2$. The direction is always radially inward towards the center of the earth.
{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
\end{itemize}
+\subsection{Parameters in section t Gravity model/Radial linear}
+\label{parameters:Gravity_20model/Radial_20linear}
+
+egin{itemize}
+\item {\it Parameter name:} { t Magnitude at surface}
+
+
+\index[prmindex]{Magnitude at surface}
+\index[prmindexfull]{Gravity model!Radial linear!Magnitude at surface}
+{\it Value:} 9.8
+
+
+{\it Default:} 9.8
+
+
+{\it Description:} Magnitude of the radial gravity vector at the surface of the domain, $m/s^2$
+
+
+{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
+\end{itemize}
+
\subsection{Parameters in section t Gravity model/Vertical}
\label{parameters:Gravity_20model/Vertical}
@@ -1139,8 +1511,10 @@
`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 ]
+
+{\it Possible values:} [Selection perturbed box|polar box|inclusion shape perturbation|mandelbox|adiabatic|spherical hexagonal perturbation|spherical gaussian perturbation|function|harmonic perturbation ]
\end{itemize}
@@ -1205,10 +1579,10 @@
{\it Default:} center
-{\it Description:} Where the initial temperature perturbation should be placed (in the center or at the boundary of the model domain).
+{\it Description:} Where the initial temperature perturbation should be placed. If 'center' is given, then the perturbation will be centered along a 'midpoint' of some sort of the bottom boundary. For example, in the case of a box geometry, this is the center of the bottom face; in the case of a spherical shell geometry, it is along the inner surface halfway between the bounding radial lines.
-{\it Possible values:} [Selection center|boundary ]
+{\it Possible values:} [Selection center ]
\item {\it Parameter name:} { t Radius}
@@ -1355,6 +1729,87 @@
{\it Possible values:} [Anything]
\end{itemize}
+\subsection{Parameters in section t Initial conditions/Harmonic perturbation}
+\label{parameters:Initial_20conditions/Harmonic_20perturbation}
+
+egin{itemize}
+\item {\it Parameter name:} { t Lateral wave number one}
+
+
+\index[prmindex]{Lateral wave number one}
+\index[prmindexfull]{Initial conditions!Harmonic perturbation!Lateral wave number one}
+{\it Value:} 3
+
+
+{\it Default:} 3
+
+
+{\it Description:} Doubled first lateral wave number of the harmonic perturbation. Equals the spherical harmonic degree in 3D spherical shells. In all other cases one equals half of a sine period over the model domain. This allows for single up-/downswings. Negative numbers reverse the sign of the perturbation but are not allowed for the spherical harmonic case.
+
+
+{\it Possible values:} [Integer range -2147483648...2147483647 (inclusive)]
+\item {\it Parameter name:} { t Lateral wave number two}
+
+
+\index[prmindex]{Lateral wave number two}
+\index[prmindexfull]{Initial conditions!Harmonic perturbation!Lateral wave number two}
+{\it Value:} 2
+
+
+{\it Default:} 2
+
+
+{\it Description:} Doubled second lateral wave number of the harmonic perturbation. Equals the spherical harmonic order in 3D spherical shells. In all other cases one equals half of a sine period over the model domain. This allows for single up-/downswings. Negative numbers reverse the sign of the perturbation.
+
+
+{\it Possible values:} [Integer range -2147483648...2147483647 (inclusive)]
+\item {\it Parameter name:} { t Magnitude}
+
+
+\index[prmindex]{Magnitude}
+\index[prmindexfull]{Initial conditions!Harmonic perturbation!Magnitude}
+{\it Value:} 1.0
+
+
+{\it Default:} 1.0
+
+
+{\it Description:} The magnitude of the Harmonic perturbation.
+
+
+{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
+\item {\it Parameter name:} { t Reference temperature}
+
+
+\index[prmindex]{Reference temperature}
+\index[prmindexfull]{Initial conditions!Harmonic perturbation!Reference temperature}
+{\it Value:} 1600.0
+
+
+{\it Default:} 1600.0
+
+
+{\it Description:} The reference temperature that is perturbed by theharmonic function. Only used in incompressible models.
+
+
+{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
+\item {\it Parameter name:} { t Vertical wave number}
+
+
+\index[prmindex]{Vertical wave number}
+\index[prmindexfull]{Initial conditions!Harmonic perturbation!Vertical wave number}
+{\it Value:} 1
+
+
+{\it Default:} 1
+
+
+{\it Description:} Doubled radial wave number of the harmonic perturbation. One equals half of a sine period over the model domain. This allows for single up-/downswings. Negative numbers reverse the sign of the perturbation.
+
+
+{\it Possible values:} [Integer range -2147483648...2147483647 (inclusive)]
+\end{itemize}
+
\subsection{Parameters in section t Initial conditions/Inclusion shape perturbation}
\label{parameters:Initial_20conditions/Inclusion_20shape_20perturbation}
@@ -1596,6 +2051,8 @@
`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.
@@ -1609,7 +2066,7 @@
`Inclusion': A material model that corresponds to the 'Inclusion' benchmark defined in Duretz et al., G-Cubed, 2011.
-{\it Possible values:} [Selection Steinberger|simple|Tan Gurnis|table|SolCx|SolKz|Inclusion ]
+{\it Possible values:} [Selection Steinberger|latent heat|simple|Tan Gurnis|table|SolCx|SolKz|Inclusion ]
\end{itemize}
@@ -1635,6 +2092,282 @@
{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
\end{itemize}
+\subsection{Parameters in section t Material model/Latent heat}
+\label{parameters:Material_20model/Latent_20heat}
+
+egin{itemize}
+\item {\it Parameter name:} { t Activation enthalpies}
+
+
+\index[prmindex]{Activation enthalpies}
+\index[prmindexfull]{Material model!Latent heat!Activation enthalpies}
+{\it Value:}
+
+
+{\it Default:}
+
+
+{\it Description:} A list of activation enthalpies for the temperature dependence of the viscosity of each phase. List must have one more entry than Phase transition depths. Units: $1/K$.
+
+
+{\it Possible values:} [List list of <[Double 0...1.79769e+308 (inclusive)]> of length 0...4294967295 (inclusive)]
+\item {\it Parameter name:} { t Composition viscosity prefactor}
+
+
+\index[prmindex]{Composition viscosity prefactor}
+\index[prmindexfull]{Material model!Latent heat!Composition viscosity prefactor}
+{\it Value:} 1.0
+
+
+{\it Default:} 1.0
+
+
+{\it Description:} A linear dependency of viscosity on composition. Dimensionless prefactor.
+
+
+{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
+\item {\it Parameter name:} { t Compressibility}
+
+
+\index[prmindex]{Compressibility}
+\index[prmindexfull]{Material model!Latent heat!Compressibility}
+{\it Value:} 5.124e-12
+
+
+{\it Default:} 5.124e-12
+
+
+{\it Description:} The value of the compressibility $\kappa$. Units: $1/Pa$.
+
+
+{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
+\item {\it Parameter name:} { t Corresponding phase for density jump}
+
+
+\index[prmindex]{Corresponding phase for density jump}
+\index[prmindexfull]{Material model!Latent heat!Corresponding phase for density jump}
+{\it Value:}
+
+
+{\it Default:}
+
+
+{\it Description:} A list of phases, which correspond to the Phase transition density jumps. The density jumps occur only in the phase that is given by this phase value. 0 stands for the 1st compositional fields, 1 for the second compositional field and -1 for none of them. List must have the same number of entries as Phase transition depths. Units: $Pa/K$.
+
+
+{\it Possible values:} [List list of <[Integer range 0...2147483647 (inclusive)]> of length 0...4294967295 (inclusive)]
+\item {\it Parameter name:} { t Density differential for compositional field 1}
+
+
+\index[prmindex]{Density differential for compositional field 1}
+\index[prmindexfull]{Material model!Latent heat!Density differential for compositional field 1}
+{\it Value:} 0
+
+
+{\it Default:} 0
+
+
+{\it Description:} If compositional fields are used, then one would frequently want to make the density depend on these fields. In this simple material model, we make the following assumptions: if no compositional fields are used in the current simulation, then the density is simply the usual one with its linear dependence on the temperature. If there are compositional fields, then the density only depends on the first one in such a way that the density has an additional term of the kind $+\Delta
ho \; c_1(\mathbf x)$. This parameter describes the value of $\Delta
ho$. Units: $kg/m^3/ extrm{unit change in composition}$.
+
+
+{\it Possible values:} [Double -1.79769e+308...1.79769e+308 (inclusive)]
+\item {\it Parameter name:} { t Phase transition Clapeyron slopes}
+
+
+\index[prmindex]{Phase transition Clapeyron slopes}
+\index[prmindexfull]{Material model!Latent heat!Phase transition Clapeyron slopes}
+{\it Value:}
+
+
+{\it Default:}
+
+
+{\it Description:} A list of Clapeyron slopes for each phase transition. A positive Clapeyron slope indicates that the phase transition will occur in a greater depth, if the temperature is higher than the one given in Phase transition temperatures and in a smaller depth, if the temperature is smaller than the one given in Phase transition temperatures. For negative slopes the other way round. List must have the same number of entries as Phase transition depths. Units: $Pa/K$.
+
+
+{\it Possible values:} [List list of <[Double -1.79769e+308...1.79769e+308 (inclusive)]> of length 0...4294967295 (inclusive)]
+\item {\it Parameter name:} { t Phase transition density jumps}
+
+
+\index[prmindex]{Phase transition density jumps}
+\index[prmindexfull]{Material model!Latent heat!Phase transition density jumps}
+{\it Value:}
+
+
+{\it Default:}
+
+
+{\it Description:} A list of density jumps at each phase transition. A positive value means that the density increases with depth. The corresponding entry in Corresponding phase for density jump determines if the density jump occurs in peridotite, eclogite or none of them.List must have the same number of entries as Phase transition depths. Units: $kg/m^3$.
+
+
+{\it Possible values:} [List list of <[Double 0...1.79769e+308 (inclusive)]> of length 0...4294967295 (inclusive)]
+\item {\it Parameter name:} { t Phase transition depths}
+
+
+\index[prmindex]{Phase transition depths}
+\index[prmindexfull]{Material model!Latent heat!Phase transition depths}
+{\it Value:}
+
+
+{\it Default:}
+
+
+{\it Description:} A list of depths where phase transitions occur. Values must monotonically increase. Units: $m$.
+
+
+{\it Possible values:} [List list of <[Double 0...1.79769e+308 (inclusive)]> of length 0...4294967295 (inclusive)]
+\item {\it Parameter name:} { t Phase transition temperatures}
+
+
+\index[prmindex]{Phase transition temperatures}
+\index[prmindexfull]{Material model!Latent heat!Phase transition temperatures}
+{\it Value:}
+
+
+{\it Default:}
+
+
+{\it Description:} A list of temperatures where phase transitions occur. Higher or lower temperatures lead to phase transition ocurring in smaller or greater depths than given in Phase transition depths, depending on the Clapeyron slope given in Phase transition Clapeyron slopes. List must have the same number of entries as Phase transition depths. Units: $K$.
+
+
+{\it Possible values:} [List list of <[Double 0...1.79769e+308 (inclusive)]> of length 0...4294967295 (inclusive)]
+\item {\it Parameter name:} { t Phase transition widths}
+
+
+\index[prmindex]{Phase transition widths}
+\index[prmindexfull]{Material model!Latent heat!Phase transition widths}
+{\it Value:}
+
+
+{\it Default:}
+
+
+{\it Description:} A list of widths for each phase transition. The phase functions are scaled with these values, leading to a jump betwen phases for a value of zero and a gradual transition for larger values. List must have the same number of entries as Phase transition depths. Units: $m$.
+
+
+{\it Possible values:} [List list of <[Double 0...1.79769e+308 (inclusive)]> of length 0...4294967295 (inclusive)]
+\item {\it Parameter name:} { t Reference density}
+
+
+\index[prmindex]{Reference density}
+\index[prmindexfull]{Material model!Latent heat!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!Latent heat!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!Latent heat!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!Latent heat!Thermal conductivity}
+{\it Value:} 2.38
+
+
+{\it Default:} 2.38
+
+
+{\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!Latent heat!Thermal expansion coefficient}
+{\it Value:} 4e-5
+
+
+{\it Default:} 4e-5
+
+
+{\it Description:} The value of the thermal expansion coefficient $eta$. Units: $1/K$.
+
+
+{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
+\item {\it Parameter name:} { t Thermal viscosity exponent}
+
+
+\index[prmindex]{Thermal viscosity exponent}
+\index[prmindexfull]{Material model!Latent heat!Thermal viscosity exponent}
+{\it Value:} 0.0
+
+
+{\it Default:} 0.0
+
+
+{\it Description:} The temperature dependence of viscosity. Dimensionless exponent.
+
+
+{\it Possible values:} [Double 0...1.79769e+308 (inclusive)]
+\item {\it Parameter name:} { t Viscosity}
+
+
+\index[prmindex]{Viscosity}
+\index[prmindexfull]{Material model!Latent heat!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)]
More information about the CIG-COMMITS
mailing list