[cig-commits] [commit] master: add a test for the composition reaction (76adf6a)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Fri May 16 18:31:08 PDT 2014
Repository : https://github.com/geodynamics/aspect
On branch : master
Link : https://github.com/geodynamics/aspect/compare/fec9fff95d1f00eabc29ab290d8294518b63b49a...7b658caa05d489ae3f8d61b8049b015e7be94f77
>---------------------------------------------------------------
commit 76adf6a11c50234352a7c7e4206c3168ddd4555b
Author: Juliane Dannberg <dannberg at gfz-potsdam.de>
Date: Sat May 17 03:27:44 2014 +0200
add a test for the composition reaction
>---------------------------------------------------------------
76adf6a11c50234352a7c7e4206c3168ddd4555b
.../composition_reaction.prm | 23 ++--
tests/composition_reaction/screen-output | 142 +++++++++++++++++++++
tests/composition_reaction/statistics | 28 ++++
3 files changed, 180 insertions(+), 13 deletions(-)
diff --git a/cookbooks/composition-reaction.prm b/tests/composition_reaction.prm
similarity index 81%
copy from cookbooks/composition-reaction.prm
copy to tests/composition_reaction.prm
index 6d5aabe..b012a8c 100644
--- a/cookbooks/composition-reaction.prm
+++ b/tests/composition_reaction.prm
@@ -1,15 +1,16 @@
-# A description of convection in a rectangular box where we prescribe
-# temporally variable boundary conditions for the velocity along the top
-# boundary. See the manual for more information.
+#########################################################
+# This is a variation of the composition-active.prm
+# parameter file, with the exception that we use the
+# material model composition reaction that accounts for
+# reactions between the two compositional fields.
############### Global parameters
set Dimension = 2
set Start time = 0
-set End time = 20
+set End time = 0.3
set Use years in output instead of seconds = false
-set Output directory = output
############### Parameters describing the model
@@ -88,7 +89,8 @@ subsection Compositional initial conditions
subsection Function
set Variable names = x,z
- set Function expression = if(z<0.5,1,0);0
+ set Function constants = pi=3.1415926
+ set Function expression = if(z<0.5+0.4*sin((x-0.5)*pi),1,0);0
end
end
@@ -111,18 +113,13 @@ end
# mesh is refined and what to do with the solution once computed
subsection Mesh refinement
set Initial adaptive refinement = 0
- set Initial global refinement = 5
+ set Initial global refinement = 4
set Time steps between mesh refinement = 0
end
subsection Postprocess
- set List of postprocessors = visualization, temperature statistics, heat flux statistics
-
- subsection Visualization
- set Time between graphical output = 0.1
- set List of output variables = density
- end
+ set List of postprocessors = composition statistics
end
diff --git a/tests/composition_reaction/screen-output b/tests/composition_reaction/screen-output
new file mode 100644
index 0000000..a54ad28
--- /dev/null
+++ b/tests/composition_reaction/screen-output
@@ -0,0 +1,142 @@
+-----------------------------------------------------------------------------
+-- 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: 5,734 (2,178+289+1,089+1,089+1,089)
+
+*** Timestep 0: t=0 seconds
+ Solving temperature system... 0 iterations.
+ Solving composition system 1... 11 iterations.
+ Solving composition system 2... 12 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 30+2 iterations.
+
+ Postprocessing:
+ Compositions min/max/mass: -0.2262/1.01/0.9695 // -0.03313/1.223/0.02968
+
+*** Timestep 1: t=0.03125 seconds
+ Solving temperature system... 11 iterations.
+ Solving composition system 1... 11 iterations.
+ Solving composition system 2... 12 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 18 iterations.
+
+ Postprocessing:
+ Compositions min/max/mass: -0.03882/1.044/0.969 // -0.02427/1.008/0.03015
+
+*** Timestep 2: t=0.0625 seconds
+ Solving temperature system... 12 iterations.
+ Solving composition system 1... 12 iterations.
+ Solving composition system 2... 12 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 29 iterations.
+
+ Postprocessing:
+ Compositions min/max/mass: -0.05886/1.072/0.9642 // -0.0257/0.9568/0.03493
+
+*** Timestep 3: t=0.09375 seconds
+ Solving temperature system... 12 iterations.
+ Solving composition system 1... 12 iterations.
+ Solving composition system 2... 13 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 28 iterations.
+
+ Postprocessing:
+ Compositions min/max/mass: -0.06562/1.075/0.9562 // -0.02472/0.9421/0.04302
+
+*** Timestep 4: t=0.125 seconds
+ Solving temperature system... 13 iterations.
+ Solving composition system 1... 12 iterations.
+ Solving composition system 2... 13 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 28 iterations.
+
+ Postprocessing:
+ Compositions min/max/mass: -0.06213/1.074/0.9474 // -0.03111/0.9764/0.05175
+
+*** Timestep 5: t=0.15625 seconds
+ Solving temperature system... 12 iterations.
+ Solving composition system 1... 12 iterations.
+ Solving composition system 2... 12 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 29 iterations.
+
+ Postprocessing:
+ Compositions min/max/mass: -0.05478/1.069/0.9394 // -0.01816/1/0.05976
+
+*** Timestep 6: t=0.1875 seconds
+ Solving temperature system... 12 iterations.
+ Solving composition system 1... 12 iterations.
+ Solving composition system 2... 12 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 27 iterations.
+
+ Postprocessing:
+ Compositions min/max/mass: -0.04742/1.058/0.9325 // -0.007577/0.9952/0.06664
+
+*** Timestep 7: t=0.21875 seconds
+ Solving temperature system... 12 iterations.
+ Solving composition system 1... 12 iterations.
+ Solving composition system 2... 11 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 27 iterations.
+
+ Postprocessing:
+ Compositions min/max/mass: -0.04069/1.041/0.9261 // -0.005582/0.9845/0.07301
+
+*** Timestep 8: t=0.25 seconds
+ Solving temperature system... 12 iterations.
+ Solving composition system 1... 12 iterations.
+ Solving composition system 2... 12 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 28 iterations.
+
+ Postprocessing:
+ Compositions min/max/mass: -0.03538/1.021/0.9199 // -0.003598/0.9756/0.07918
+
+*** Timestep 9: t=0.28125 seconds
+ Solving temperature system... 12 iterations.
+ Solving composition system 1... 11 iterations.
+ Solving composition system 2... 12 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 26 iterations.
+
+ Postprocessing:
+ Compositions min/max/mass: -0.03189/1.019/0.9144 // -0.002049/0.9688/0.08467
+
+*** Timestep 10: t=0.3 seconds
+ Solving temperature system... 10 iterations.
+ Solving composition system 1... 10 iterations.
+ Solving composition system 2... 12 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 28 iterations.
+
+ Postprocessing:
+ Compositions min/max/mass: -0.03066/1.019/0.9104 // -0.002088/0.9857/0.08875
+
+Termination requested by criterion: end time
+
+
++---------------------------------------------+------------+------------+
+| Total wallclock time elapsed since start | 4.77s | |
+| | | |
+| Section | no. calls | wall time | % of total |
++---------------------------------+-----------+------------+------------+
+| Assemble Stokes system | 11 | 1.13s | 24% |
+| Assemble composition system | 22 | 0.937s | 20% |
+| Assemble temperature system | 11 | 0.576s | 12% |
+| Build Stokes preconditioner | 11 | 0.874s | 18% |
+| Build composition preconditioner| 22 | 0.133s | 2.8% |
+| Build temperature preconditioner| 11 | 0.0629s | 1.3% |
+| Solve Stokes system | 11 | 0.719s | 15% |
+| Solve composition system | 22 | 0.101s | 2.1% |
+| Solve temperature system | 11 | 0.0467s | 0.98% |
+| Initialization | 2 | 0.053s | 1.1% |
+| Postprocessing | 11 | 0.0272s | 0.57% |
+| Setup dof systems | 1 | 0.0615s | 1.3% |
++---------------------------------+-----------+------------+------------+
+
diff --git a/tests/composition_reaction/statistics b/tests/composition_reaction/statistics
new file mode 100644
index 0000000..eecdcfe
--- /dev/null
+++ b/tests/composition_reaction/statistics
@@ -0,0 +1,28 @@
+# 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: Number of composition degrees of freedom
+# 7: Iterations for temperature solver
+# 8: Iterations for composition solver 1
+# 9: Iterations for composition solver 2
+# 10: Iterations for Stokes solver
+# 11: Time step size (seconds)
+# 12: Minimal value for composition 0
+# 13: Maximal value for composition 0
+# 14: Global mass for composition 0
+# 15: Minimal value for composition 1
+# 16: Maximal value for composition 1
+# 17: Global mass for composition 1
+ 0 0.0000e+00 256 2467 1089 1089 0 11 12 32 3.1250e-02 -2.26183349e-01 1.01030114e+00 9.69455003e-01 -3.31278643e-02 1.22330718e+00 2.96769411e-02
+ 1 3.1250e-02 256 2467 1089 1089 11 11 12 18 3.1250e-02 -3.88158899e-02 1.04404545e+00 9.68991014e-01 -2.42706699e-02 1.00778984e+00 3.01517042e-02
+ 2 6.2500e-02 256 2467 1089 1089 12 12 12 29 3.1250e-02 -5.88558597e-02 1.07195015e+00 9.64228221e-01 -2.56952565e-02 9.56836539e-01 3.49275891e-02
+ 3 9.3750e-02 256 2467 1089 1089 12 12 13 28 3.1250e-02 -6.56193323e-02 1.07483229e+00 9.56158422e-01 -2.47217533e-02 9.42149651e-01 4.30162373e-02
+ 4 1.2500e-01 256 2467 1089 1089 13 12 13 28 3.1250e-02 -6.21268685e-02 1.07382996e+00 9.47392405e-01 -3.11142057e-02 9.76351197e-01 5.17484148e-02
+ 5 1.5625e-01 256 2467 1089 1089 12 12 12 29 3.1250e-02 -5.47786350e-02 1.06889311e+00 9.39374743e-01 -1.81602956e-02 1.00043483e+00 5.97632415e-02
+ 6 1.8750e-01 256 2467 1089 1089 12 12 12 27 3.1250e-02 -4.74164607e-02 1.05753878e+00 9.32505293e-01 -7.57681282e-03 9.95185287e-01 6.66415171e-02
+ 7 2.1875e-01 256 2467 1089 1089 12 12 11 27 3.1250e-02 -4.06931970e-02 1.04120271e+00 9.26116868e-01 -5.58161493e-03 9.84505910e-01 7.30074549e-02
+ 8 2.5000e-01 256 2467 1089 1089 12 12 12 28 3.1250e-02 -3.53818226e-02 1.02113392e+00 9.19931737e-01 -3.59807045e-03 9.75580067e-01 7.91817573e-02
+ 9 2.8125e-01 256 2467 1089 1089 12 11 12 26 1.8750e-02 -3.18874841e-02 1.01880895e+00 9.14432953e-01 -2.04944019e-03 9.68804576e-01 8.46703424e-02
+10 3.0000e-01 256 2467 1089 1089 10 10 12 28 3.1250e-02 -3.06581313e-02 1.01947680e+00 9.10352847e-01 -2.08847675e-03 9.85720626e-01 8.87458550e-02
More information about the CIG-COMMITS
mailing list