[cig-commits] commit 2490 by heister to /var/svn/dealii/aspect

dealii.demon at gmail.com dealii.demon at gmail.com
Sat Apr 12 05:53:18 PDT 2014


Revision 2490

add test for maximum timestep

A   trunk/aspect/tests/maxtimestep/
A   trunk/aspect/tests/maxtimestep/screen-output
A   trunk/aspect/tests/maxtimestep/statistics
A   trunk/aspect/tests/maxtimestep.prm


http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=2490&peg=2490

Diff:
Added: trunk/aspect/tests/maxtimestep/screen-output
===================================================================
--- trunk/aspect/tests/maxtimestep/screen-output	                        (rev 0)
+++ trunk/aspect/tests/maxtimestep/screen-output	2014-04-12 12:53:16 UTC (rev 2490)
@@ -0,0 +1,69 @@
+-----------------------------------------------------------------------------
+-- This is ASPECT, the Advanced Simulator for Problems in Earth's ConvecTion.
+--     . running in DEBUG 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... 20 iterations.
+
+   Postprocessing:
+
+*** Timestep 1:  t=42 seconds
+   Solving temperature system... 5 iterations.
+   Solving Stokes system... 0 iterations.
+
+   Postprocessing:
+
+*** Timestep 2:  t=84 seconds
+   Solving temperature system... 3 iterations.
+   Solving Stokes system... 0 iterations.
+
+   Postprocessing:
+
+*** Timestep 3:  t=126 seconds
+   Solving temperature system... 3 iterations.
+   Solving Stokes system... 0 iterations.
+
+   Postprocessing:
+
+*** Timestep 4:  t=168 seconds
+   Solving temperature system... 3 iterations.
+   Solving Stokes system... 0 iterations.
+
+   Postprocessing:
+
+*** Timestep 5:  t=210 seconds
+   Solving temperature system... 3 iterations.
+   Solving Stokes system... 0 iterations.
+
+   Postprocessing:
+
+Termination requested by criterion: end step
+*** Snapshot created!
+
+
+
++---------------------------------------------+------------+------------+
+| Total wallclock time elapsed since start    |      1.37s |            |
+|                                             |            |            |
+| Section                         | no. calls |  wall time | % of total |
++---------------------------------+-----------+------------+------------+
+| Assemble Stokes system          |         6 |     0.257s |        19% |
+| Assemble temperature system     |         6 |     0.595s |        43% |
+| Build Stokes preconditioner     |         1 |     0.173s |        13% |
+| Build temperature preconditioner|         6 |    0.0475s |       3.5% |
+| Solve Stokes system             |         6 |    0.0631s |       4.6% |
+| Solve temperature system        |         6 |    0.0117s |      0.85% |
+| Create snapshot                 |         1 |    0.0321s |       2.3% |
+| Initialization                  |         2 |    0.0243s |       1.8% |
+| Postprocessing                  |         6 |  0.000359s |     0.026% |
+| Setup dof systems               |         1 |     0.102s |       7.4% |
++---------------------------------+-----------+------------+------------+
+

Added: trunk/aspect/tests/maxtimestep/statistics
===================================================================
--- trunk/aspect/tests/maxtimestep/statistics	                        (rev 0)
+++ trunk/aspect/tests/maxtimestep/statistics	2014-04-12 12:53:16 UTC (rev 2490)
@@ -0,0 +1,14 @@
+# 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)
+0 0.0000e+00 1024 9539 4225 0 20 4.2000e+01 
+1 4.2000e+01 1024 9539 4225 5  0 4.2000e+01 
+2 8.4000e+01 1024 9539 4225 3  0 4.2000e+01 
+3 1.2600e+02 1024 9539 4225 3  0 4.2000e+01 
+4 1.6800e+02 1024 9539 4225 3  0 4.2000e+01 
+5 2.1000e+02 1024 9539 4225 3  0 4.2000e+01 

Copied: trunk/aspect/tests/maxtimestep.prm (from rev 2489, trunk/aspect/tests/box-end_time_1e7-terminate.prm)
===================================================================
--- trunk/aspect/tests/maxtimestep.prm	                        (rev 0)
+++ trunk/aspect/tests/maxtimestep.prm	2014-04-12 12:53:16 UTC (rev 2490)
@@ -0,0 +1,80 @@
+# a test that verifies that setting an end time earlier than when
+# we achieve steady state still terminates the program at the right
+# time
+
+set Dimension = 2
+set CFL number                             = 1.0
+set End time                               = 1e7
+set Start time                             = 0
+set Adiabatic surface temperature          = 0
+set Surface pressure                       = 0
+set Use years in output instead of seconds = false  # default: true
+set Nonlinear solver scheme                = IMPES
+
+
+subsection Boundary temperature model
+  set Model name = box
+end
+
+
+
+subsection Gravity model
+  set Model name = vertical
+end
+
+
+subsection Geometry model
+  set Model name = box
+
+  subsection Box
+    set X extent = 1.2 # default: 1
+    set Y extent = 1
+    set Z extent = 1
+  end
+end
+
+
+subsection Initial conditions
+  set Model name = perturbed box
+end
+
+
+subsection Material model
+  set Model name = simple
+
+  subsection Simple model
+    set Reference density             = 1    # default: 3300
+    set Reference specific heat       = 1250
+    set Reference temperature         = 1    # default: 293
+    set Thermal conductivity          = 1e-6 # default: 4.7
+    set Thermal expansion coefficient = 2e-5
+    set Viscosity                     = 1    # default: 5e24
+  end
+end
+
+
+subsection Mesh refinement
+  set Initial adaptive refinement        = 0
+  set Initial global refinement          = 5
+end
+
+
+subsection Model settings
+  set Fixed temperature boundary indicators   =
+  set Prescribed velocity boundary indicators =
+  set Tangential velocity boundary indicators = 1
+  set Zero velocity boundary indicators       = 0, 2, 3
+end
+
+subsection Postprocess
+  set List of postprocessors =
+end
+
+subsection Termination criteria
+  set Checkpoint on termination = true
+  set Termination criteria      = end step
+  set End step                  = 5
+
+end
+
+set Maximum time step                      = 42


More information about the CIG-COMMITS mailing list