[cig-commits] [commit] master: Add a test that verifies that we can deal with a situation where there is really absolutely nothing to do since both temperature and velocity are and remain zero. (0a06f6d)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/aspect

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/a35e4f5e47a47903957f7c84c785299a49101c46...2df080f4905a6be43fd1729ea0d6a7a956cce070

>---------------------------------------------------------------

commit 0a06f6de7f8738b149dc301511c73ca5320e8751
Author: Wolfgang Bangerth <bangerth at math.tamu.edu>
Date:   Wed May 21 12:37:43 2014 -0500

    Add a test that verifies that we can deal with a situation where there is really absolutely nothing to do since both temperature and velocity are and remain zero.


>---------------------------------------------------------------

0a06f6de7f8738b149dc301511c73ca5320e8751
 ...ent-temperature-bc.prm => zero-temperature.prm} | 18 ++++----
 tests/zero-temperature/screen-output               | 51 ++++++++++++++++++++++
 tests/{latent_heat => zero-temperature}/statistics |  5 ++-
 3 files changed, 63 insertions(+), 11 deletions(-)

diff --git a/tests/time-dependent-temperature-bc.prm b/tests/zero-temperature.prm
similarity index 87%
copy from tests/time-dependent-temperature-bc.prm
copy to tests/zero-temperature.prm
index 5b07e30..e0ae480 100644
--- a/tests/time-dependent-temperature-bc.prm
+++ b/tests/zero-temperature.prm
@@ -1,14 +1,13 @@
-# originally taken from diffusion.prm, but intended to test that we can deal
-# with time dependent temperature boundary conditions
-
-
+# derived from diffusion.prm. we have no velocity here, and the temperature
+# boundary values are all zero as well, so there is nothing really to do.
+# check that this works.
 
 set Dimension = 2
 
 
-set CFL number                             = 0.01
+set CFL number                             = 1
 
-set End time                               = 4e7
+set End time                               = 1
 
 
 set Resume computation                     = false
@@ -25,11 +24,11 @@ set Nonlinear solver scheme                = IMPES
 
 
 subsection Boundary temperature model
-  set Model name = time_dep_box
+  set Model name = box
 
-  subsection Time_Dep_Box
+  subsection Box
     set Bottom temperature = 0
-    set Top temperature    = 1 # the model interprets this as T(t)=1*t
+    set Top temperature    = 0
   end
 end
 
@@ -141,3 +140,4 @@ end
 subsection Postprocess
   set List of postprocessors = temperature statistics
 end
+
diff --git a/tests/zero-temperature/screen-output b/tests/zero-temperature/screen-output
new file mode 100644
index 0000000..3c8947a
--- /dev/null
+++ b/tests/zero-temperature/screen-output
@@ -0,0 +1,51 @@
+-----------------------------------------------------------------------------
+-- 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: 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... 0 iterations.
+
+   Postprocessing:
+     Temperature min/avg/max: 0 K, 0 K, 0 K
+
+*** Timestep 1:  t=0.5 seconds
+   Solving temperature system... 0 iterations.
+   Solving Stokes system... 0 iterations.
+
+   Postprocessing:
+     Temperature min/avg/max: 0 K, 0 K, 0 K
+
+*** Timestep 2:  t=1 seconds
+   Solving temperature system... 0 iterations.
+   Solving Stokes system... 0 iterations.
+
+   Postprocessing:
+     Temperature min/avg/max: 0 K, 0 K, 0 K
+
+Termination requested by criterion: end time
+
+
++---------------------------------------------+------------+------------+
+| Total wallclock time elapsed since start    |     0.458s |            |
+|                                             |            |            |
+| Section                         | no. calls |  wall time | % of total |
++---------------------------------+-----------+------------+------------+
+| Assemble Stokes system          |         3 |    0.0874s |        19% |
+| Assemble temperature system     |         3 |     0.129s |        28% |
+| Build Stokes preconditioner     |         1 |    0.0779s |        17% |
+| Build temperature preconditioner|         3 |    0.0444s |       9.7% |
+| Solve Stokes system             |         3 |    0.0101s |       2.2% |
+| Solve temperature system        |         3 |   0.00318s |       0.7% |
+| Initialization                  |         2 |    0.0216s |       4.7% |
+| Postprocessing                  |         3 |   0.00506s |       1.1% |
+| Setup dof systems               |         1 |    0.0644s |        14% |
++---------------------------------+-----------+------------+------------+
+
diff --git a/tests/latent_heat/statistics b/tests/zero-temperature/statistics
similarity index 56%
copy from tests/latent_heat/statistics
copy to tests/zero-temperature/statistics
index ade5039..13a17af 100644
--- a/tests/latent_heat/statistics
+++ b/tests/zero-temperature/statistics
@@ -9,5 +9,6 @@
 # 9: Minimal temperature (K)
 # 10: Average temperature (K)
 # 11: Maximal temperature (K)
-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 
+0 0.0000e+00 1024 9539 4225 0 0 5.0000e-01 0.00000000e+00 0.00000000e+00 0.00000000e+00 
+1 5.0000e-01 1024 9539 4225 0 0 5.0000e-01 0.00000000e+00 0.00000000e+00 0.00000000e+00 
+2 1.0000e+00 1024 9539 4225 0 0 5.0000e-01 0.00000000e+00 0.00000000e+00 0.00000000e+00