[cig-commits] commit 2337 by dannberg to /var/svn/dealii/aspect
dealii.demon at gmail.com
dealii.demon at gmail.com
Wed Mar 19 07:30:45 PDT 2014
Revision 2337
add a material model that includes latent heat and a corresponding benchmark case in the manual
A trunk/aspect/cookbooks/latent-heat.prm
A trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/
A trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-results-1.png
A trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-results-2.png
A trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-setup.png
A trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-temperature.png
U trunk/aspect/doc/manual/manual.tex
U trunk/aspect/include/aspect/material_model/interface.h
A trunk/aspect/include/aspect/material_model/latent_heat.h
A trunk/aspect/source/material_model/latent_heat.cc
U trunk/aspect/source/simulator/assembly.cc
http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=2337&peg=2337
Diff:
Added: trunk/aspect/cookbooks/latent-heat.prm
===================================================================
--- trunk/aspect/cookbooks/latent-heat.prm (rev 0)
+++ trunk/aspect/cookbooks/latent-heat.prm 2014-03-19 14:30:43 UTC (rev 2337)
@@ -0,0 +1,162 @@
+############### Global parameters
+# We use a 2d setup. Since it takes some time for
+# the model to reach a steady state we set the
+# end time to approximately 15 billion years.
+
+set Dimension = 2
+
+set Start time = 0
+set End time = 5e17
+set Use years in output instead of seconds = false
+
+set Output directory = output
+
+
+subsection Geometry model
+ set Model name = box
+
+ subsection Box
+ set X extent = 1000000
+ set Y extent = 1000000
+ end
+end
+
+
+subsection Gravity model
+ set Model name = vertical
+ subsection Vertical
+ set Magnitude = 10.0
+ end
+end
+
+
+subsection Model settings
+
+ # As we only want to look at the effects of latent heating, we disable all
+ # the other heating terms.
+ set Include adiabatic heating = false
+ set Include latent heat = true
+ set Include shear heating = false
+ set Radiogenic heating rate = 0
+
+ # We only fix the temperature at the upper boundary, the other boundaries
+ # are isolating. To guarantuee a steady downward flow, we fix the velocity
+ # at the top and bottom, and set it to free slip on the sides.
+ set Fixed temperature boundary indicators = 3
+ set Prescribed velocity boundary indicators = 2:function, 3:function
+ set Tangential velocity boundary indicators = 0, 1
+end
+
+
+############### Boundary conditions
+# We set the top temperature to T1=1000K.
+subsection Boundary temperature model
+ set Model name = box
+ subsection Box
+ set Top temperature = 1000
+ end
+end
+
+# We prescribe a constant downward flow.
+subsection Boundary velocity model
+ subsection Function
+ set Function expression = 0;-2.1422e-11
+ set Variable names = x,y
+ end
+end
+
+subsection Initial conditions
+ set Model name = function
+ subsection Function
+ set Function expression = 1000.0
+ set Variable names = x,y
+ end
+end
+
+
+subsection Material model
+ set Model name = latent heat
+ subsection Latent heat
+
+ # The change of density across the phase transition. Together with the
+ # Clapeyron slope, this is what determines the entropy change.
+ set Phase transition density jumps = 115.6
+ set Corresponding phase for density jump = 0
+
+ # If the temperature is equal to the phase transition temperature, the
+ # phase transition will occur at the phase transition depth. However,
+ # if the temperature deviates from this value, the Clapeyron slope
+ # determines how much the pressure (and depth) of the phase boundary
+ # changes. Here, the phase transition will be in the middle of the box
+ # for T=T1.
+ set Phase transition depths = 500000
+ set Phase transition temperatures = 1000
+ set Phase transition Clapeyron slopes = 1e7
+
+ # We set the width of the phase transition to 5 km. You may want to
+ # change this parameter to see how latent heating depends on the width
+ # of the phase transition.
+ set Phase transition widths = 5000
+
+ set Reference density = 3400
+ set Reference specific heat = 1000
+ set Reference temperature = 1000
+ set Thermal conductivity = 2.38
+
+ # We set the thermal expansion amd the compressibility to zero, so that
+ # all temperature (and density) changes are caused by advection, diffusion
+ # and latent heating.
+ set Thermal expansion coefficient = 0.0
+ set Compressibility = 0.0
+
+ # Viscosity is constant.
+ set Thermal viscosity exponent = 0.0
+ set Viscosity = 8.44e21
+ set Viscosity prefactors = 1.0, 1.0
+ set Composition viscosity prefactor = 1.0
+ set Activation enthalpies = 3.9473e-3, 3.9473e-3
+ end
+end
+
+
+subsection Mesh refinement
+ set Initial adaptive refinement = 0
+ set Initial global refinement = 7
+ set Time steps between mesh refinement = 0
+
+end
+
+
+subsection Discretization
+ subsection Stabilization parameters
+ # The exponent $lpha$ in the entropy viscosity stabilization. Units:
+ # None.
+ set alpha = 2
+
+ # The $eta$ factor in the artificial viscosity stabilization. An
+ # appropriate value for 2d is 0.052 and 0.078 for 3d. Units: None.
+ set beta = 0.078
+
+ # The $c_R$ factor in the entropy viscosity stabilization. Units: None.
+ set cR = 0.5 # default: 0.11
+ end
+end
+
+
+subsection Postprocess
+
+ set List of postprocessors = visualization
+
+ subsection Visualization
+
+ set Number of grouped files = 0
+ set Output format = vtu
+
+ # We are only interested in the last timestep (when the system hast reached
+ # a steady state). For following the development of the system or checking
+ # if the solution already reached steady state, this parameter can be set
+ # to a smaller value.
+ set Time between graphical output = 5e17
+ set List of output variables = density
+ end
+end
Added: trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-results-1.png
===================================================================
(Binary files differ)
Index: trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-results-1.png
===================================================================
--- trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-results-1.png 2014-03-18 17:58:13 UTC (rev 2336)
+++ trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-results-1.png 2014-03-19 14:30:43 UTC (rev 2337)
Property changes on: trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-results-1.png
___________________________________________________________________
Added: svn:mime-type
## -0,0 +1 ##
+application/octet-stream
\ No newline at end of property
Added: trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-results-2.png
===================================================================
(Binary files differ)
Index: trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-results-2.png
===================================================================
--- trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-results-2.png 2014-03-18 17:58:13 UTC (rev 2336)
+++ trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-results-2.png 2014-03-19 14:30:43 UTC (rev 2337)
Property changes on: trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-results-2.png
___________________________________________________________________
Added: svn:mime-type
## -0,0 +1 ##
+application/octet-stream
\ No newline at end of property
Added: trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-setup.png
===================================================================
(Binary files differ)
Index: trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-setup.png
===================================================================
--- trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-setup.png 2014-03-18 17:58:13 UTC (rev 2336)
+++ trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-setup.png 2014-03-19 14:30:43 UTC (rev 2337)
Property changes on: trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-setup.png
___________________________________________________________________
Added: svn:mime-type
## -0,0 +1 ##
+application/octet-stream
\ No newline at end of property
Added: trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-temperature.png
===================================================================
(Binary files differ)
Index: trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-temperature.png
===================================================================
--- trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-temperature.png 2014-03-18 17:58:13 UTC (rev 2336)
+++ trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-temperature.png 2014-03-19 14:30:43 UTC (rev 2337)
Property changes on: trunk/aspect/doc/manual/cookbooks/benchmarks/latent-heat/latent-heat-temperature.png
___________________________________________________________________
Added: svn:mime-type
## -0,0 +1 ##
+application/octet-stream
\ No newline at end of property
Modified: trunk/aspect/doc/manual/manual.tex
===================================================================
--- trunk/aspect/doc/manual/manual.tex 2014-03-18 17:58:13 UTC (rev 2336)
+++ trunk/aspect/doc/manual/manual.tex 2014-03-19 14:30:43 UTC (rev 2337)
@@ -266,6 +266,7 @@
\
\label{eq:temperature}
ho C_p \left(rac{\partial T}{\partial t} + \mathbf u
More information about the CIG-COMMITS
mailing list