[cig-commits] [commit] master: Add a test and a future cookbook for radiogenic heating. (9cdad95)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Fri May 16 18:39:34 PDT 2014
Repository : https://github.com/geodynamics/aspect
On branch : master
Link : https://github.com/geodynamics/aspect/compare/7b658caa05d489ae3f8d61b8049b015e7be94f77...cc43d3ac6f571573118eb54c2103bf92dfe8355f
>---------------------------------------------------------------
commit 9cdad95396a2d412fcc41a742285faa564002472
Author: Rene Gassmoeller <R.Gassmoeller at mailbox.org>
Date: Fri May 16 09:19:00 2014 -0500
Add a test and a future cookbook for radiogenic heating.
>---------------------------------------------------------------
9cdad95396a2d412fcc41a742285faa564002472
.../radiogenic_heating.prm} | 76 ++++----------
tests/radiogenic_heating.prm | 111 +++++++++++++++++++++
.../screen_output} | 28 +++---
tests/radiogenic_heating/statistics | 13 +++
4 files changed, 157 insertions(+), 71 deletions(-)
diff --git a/cookbooks/latent-heat.prm b/cookbooks/future/radiogenic_heating.prm
similarity index 50%
copy from cookbooks/latent-heat.prm
copy to cookbooks/future/radiogenic_heating.prm
index 58aa87e..9a48afc 100644
--- a/cookbooks/latent-heat.prm
+++ b/cookbooks/future/radiogenic_heating.prm
@@ -1,16 +1,23 @@
############### 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.
+# end time to approximately 2 million years.
set Dimension = 2
set Start time = 0
-set End time = 5e17
+set End time = 2e6
set Use years in output instead of seconds = false
set Output directory = output
+subsection Heating model
+ set Model name = constant heating
+
+ subsection Constant heating
+ set Radiogenic heating rate = 1
+ end
+end
subsection Geometry model
set Model name = box
@@ -32,12 +39,12 @@ end
subsection Model settings
- # As we only want to look at the effects of latent heating, we disable all
+ # As we only want to look at the effects of radiogenic, we disable all
# the other heating terms.
set Include adiabatic heating = false
- set Include latent heat = true
+ set Include latent heat = false
set Include shear heating = false
- set Radiogenic heating rate = 0
+ set Radiogenic heating rate = 1e5
# We only fix the temperature at the upper boundary, the other boundaries
# are isolating. To guarantuee a steady downward flow, we fix the velocity
@@ -60,7 +67,7 @@ end
# We prescribe a constant downward flow.
subsection Boundary velocity model
subsection Function
- set Function expression = 0;-2.1422e-11
+ set Function expression = 0;-1
set Variable names = x,y
end
end
@@ -75,74 +82,31 @@ 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 Model name = simple
+
+ subsection Simple model
set Reference density = 3400
set Reference specific heat = 1000
set Reference temperature = 1000
- set Thermal conductivity = 2.38
+ set Thermal conductivity = 0.0
- # 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
+ set Viscosity = 1e21
end
end
subsection Mesh refinement
set Initial adaptive refinement = 0
- set Initial global refinement = 7
+ set Initial global refinement = 6
set Time steps between mesh refinement = 0
end
-subsection Discretization
- subsection Stabilization parameters
- # The exponent $\alpha$ in the entropy viscosity stabilization. Units:
- # None.
- set alpha = 2
-
- # The $\beta$ 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
@@ -156,7 +120,7 @@ subsection Postprocess
# 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 Time between graphical output = 1e4
set List of output variables = density
end
end
diff --git a/tests/radiogenic_heating.prm b/tests/radiogenic_heating.prm
new file mode 100644
index 0000000..bca7cac
--- /dev/null
+++ b/tests/radiogenic_heating.prm
@@ -0,0 +1,111 @@
+############### Global parameters
+# This is a test for the radiogenic heating.
+
+set Dimension = 2
+
+set Start time = 0
+set End time = 1000
+set Use years in output instead of seconds = false
+
+set Output directory = radiogenic_heating
+
+subsection Heating model
+ set Model name = constant heating
+
+ subsection Constant heating
+ set Radiogenic heating rate = 1
+ end
+end
+
+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 radiogenic, we disable all
+ # the other heating terms.
+ set Include adiabatic heating = false
+ set Include latent heat = false
+ set Include shear heating = false
+
+ # 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;-1
+ 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 = simple
+
+ subsection Simple model
+
+ set Reference density = 3400
+ set Reference specific heat = 1000
+ set Reference temperature = 1000
+ set Thermal conductivity = 0.0
+
+ set Thermal expansion coefficient = 0.0
+
+ set Thermal viscosity exponent = 0.0
+ set Viscosity = 1e21
+ end
+end
+
+
+subsection Mesh refinement
+ set Initial adaptive refinement = 0
+ set Initial global refinement = 5
+ set Time steps between mesh refinement = 0
+
+end
+
+
+subsection Postprocess
+
+ set List of postprocessors = temperature statistics
+
+end
diff --git a/tests/compressibility_iterated_stokes/screen-output b/tests/radiogenic_heating/screen_output
similarity index 51%
copy from tests/compressibility_iterated_stokes/screen-output
copy to tests/radiogenic_heating/screen_output
index b1514b7..8196477 100644
--- a/tests/compressibility_iterated_stokes/screen-output
+++ b/tests/radiogenic_heating/screen_output
@@ -1,30 +1,28 @@
-----------------------------------------------------------------------------
--- This is ASPECT, the Advanced Simulator for Problems in Earth's ConvecTion.
--- . running in DEBUG mode
+-- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion.
+-- . running in OPTIMIZED mode
-- . running with 1 MPI process
-- . using Trilinos
-----------------------------------------------------------------------------
-Loading shared library <./libcompressibility_iterated_stokes.so>
-
Number of active cells: 1,024 (on 6 levels)
Number of degrees of freedom: 13,764 (8,450+1,089+4,225)
*** Timestep 0: t=0 seconds
+ Solving temperature system... 0 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 25 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 1000 K, 1000 K, 1000 K
+
+*** Timestep 1: t=1000 seconds
+ Solving temperature system... 9 iterations.
Rebuilding Stokes preconditioner...
- Solving Stokes system... 26 iterations.
- Nonlinear Stokes residual: 21.2064
- Solving Stokes system... 21 iterations.
- Nonlinear Stokes residual: 0.293291
- Solving Stokes system... 16 iterations.
- Nonlinear Stokes residual: 0.0159069
- Solving Stokes system... 10 iterations.
- Nonlinear Stokes residual: 0.000575274
- Solving Stokes system... 5 iterations.
- Nonlinear Stokes residual: 1.44305e-05
+ Solving Stokes system... 0 iterations.
Postprocessing:
- Top/bottom flux: 1/1
+ Temperature min/avg/max: 1000 K, 1001 K, 1001 K
Termination requested by criterion: end time
diff --git a/tests/radiogenic_heating/statistics b/tests/radiogenic_heating/statistics
new file mode 100644
index 0000000..04c5ccd
--- /dev/null
+++ b/tests/radiogenic_heating/statistics
@@ -0,0 +1,13 @@
+# 1: Time step number
+# 2: Time (seconds)
+# 3: Number of mesh cells
+# 4: Number of Stokes degrees of freedom
+# 5: Number of temperature degrees of freedom
+# 6: Iterations for temperature solver
+# 7: Iterations for Stokes solver
+# 8: Time step size (seconds)
+# 9: Minimal temperature (K)
+# 10: Average temperature (K)
+# 11: Maximal temperature (K)
+0 0.0000e+00 1024 9539 4225 0 25 1.0000e+03 1.00000000e+03 1.00000000e+03 1.00000000e+03
+1 1.0000e+03 1024 9539 4225 9 0 1.5625e+04 1.00000000e+03 1.00099578e+03 1.00110647e+03
More information about the CIG-COMMITS
mailing list