[cig-commits] [commit] master: add a test for latent heat (426fbfc)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Tue May 20 07:15:06 PDT 2014
Repository : https://github.com/geodynamics/aspect
On branch : master
Link : https://github.com/geodynamics/aspect/compare/5dd9f8e83db40d26a5659c449816fc72358f3746...426fbfcad3568313f61490f77127c5e1c97a005d
>---------------------------------------------------------------
commit 426fbfcad3568313f61490f77127c5e1c97a005d
Author: Juliane Dannberg <dannberg at gfz-potsdam.de>
Date: Mon May 19 22:43:20 2014 +0200
add a test for latent heat
>---------------------------------------------------------------
426fbfcad3568313f61490f77127c5e1c97a005d
cookbooks/latent-heat.prm => tests/latent_heat.prm | 32 +++++----------
tests/latent_heat/screen-output | 45 ++++++++++++++++++++++
.../{radiogenic_heating => latent_heat}/statistics | 4 +-
3 files changed, 56 insertions(+), 25 deletions(-)
diff --git a/cookbooks/latent-heat.prm b/tests/latent_heat.prm
similarity index 82%
copy from cookbooks/latent-heat.prm
copy to tests/latent_heat.prm
index 8288722..1ebd9fd 100644
--- a/cookbooks/latent-heat.prm
+++ b/tests/latent_heat.prm
@@ -1,16 +1,15 @@
-############### 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.
+#########################################################
+# This is a model setup to test the latent heat generation.
+# Material flows in from the top and crosses a phase
+# transition, where energy is released and the temperature
+# increases.
set Dimension = 2
set Start time = 0
-set End time = 5e17
+set End time = 1e15
set Use years in output instead of seconds = false
-set Output directory = output
-
subsection Geometry model
set Model name = box
@@ -95,7 +94,7 @@ subsection Material model
# 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 Phase transition widths = 20000
set Reference density = 3400
set Reference specific heat = 1000
@@ -120,7 +119,7 @@ end
subsection Mesh refinement
set Initial adaptive refinement = 0
- set Initial global refinement = 7
+ set Initial global refinement = 4
set Time steps between mesh refinement = 0
end
@@ -144,18 +143,5 @@ 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
+ set List of postprocessors = temperature statistics
end
diff --git a/tests/latent_heat/screen-output b/tests/latent_heat/screen-output
new file mode 100644
index 0000000..c8ccb7f
--- /dev/null
+++ b/tests/latent_heat/screen-output
@@ -0,0 +1,45 @@
+-----------------------------------------------------------------------------
+-- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion.
+-- . running in OPTIMIZED mode
+-- . running with 1 MPI process
+-- . using Trilinos
+-----------------------------------------------------------------------------
+
+Number of active cells: 256 (on 5 levels)
+Number of degrees of freedom: 3,556 (2,178+289+1,089)
+
+*** 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=1e+15 seconds
+ Solving temperature system... 85 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 20 iterations.
+
+ Postprocessing:
+ Temperature min/avg/max: 1000 K, 1001 K, 1003 K
+
+Termination requested by criterion: end time
+
+
++---------------------------------------------+------------+------------+
+| Total wallclock time elapsed since start | 0.563s | |
+| | | |
+| Section | no. calls | wall time | % of total |
++---------------------------------+-----------+------------+------------+
+| Assemble Stokes system | 2 | 0.107s | 19% |
+| Assemble temperature system | 2 | 0.084s | 15% |
+| Build Stokes preconditioner | 2 | 0.123s | 22% |
+| Build temperature preconditioner| 2 | 0.012s | 2.1% |
+| Solve Stokes system | 2 | 0.102s | 18% |
+| Solve temperature system | 2 | 0.0271s | 4.8% |
+| Initialization | 2 | 0.0411s | 7.3% |
+| Postprocessing | 2 | 0.00273s | 0.49% |
+| Setup dof systems | 1 | 0.0372s | 6.6% |
++---------------------------------+-----------+------------+------------+
+
diff --git a/tests/radiogenic_heating/statistics b/tests/latent_heat/statistics
similarity index 70%
copy from tests/radiogenic_heating/statistics
copy to tests/latent_heat/statistics
index 04c5ccd..ade5039 100644
--- a/tests/radiogenic_heating/statistics
+++ b/tests/latent_heat/statistics
@@ -9,5 +9,5 @@
# 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
+0 0.0000e+00 256 2467 1089 0 25 1.0000e+15 1.00000000e+03 1.00000000e+03 1.00000000e+03
+1 1.0000e+15 256 2467 1089 85 20 1.3258e+15 1.00000000e+03 1.00145182e+03 1.00259980e+03
More information about the CIG-COMMITS
mailing list