[cig-commits] r1357 - in trunk/aspect/tests: . no_flow
bangerth at dealii.org
bangerth at dealii.org
Fri Nov 9 23:13:26 PST 2012
Author: bangerth
Date: 2012-11-10 00:13:25 -0700 (Sat, 10 Nov 2012)
New Revision: 1357
Added:
trunk/aspect/tests/no_flow.prm
trunk/aspect/tests/no_flow/
trunk/aspect/tests/no_flow/screen-output
trunk/aspect/tests/no_flow/statistics
Log:
add
Added: trunk/aspect/tests/no_flow/screen-output
===================================================================
--- trunk/aspect/tests/no_flow/screen-output (rev 0)
+++ trunk/aspect/tests/no_flow/screen-output 2012-11-10 07:13:25 UTC (rev 1357)
@@ -0,0 +1,140 @@
+Number of active cells: 16 (on 3 levels)
+Number of degrees of freedom: 268 (162+25+81)
+
+*** Timestep 0: t=0 seconds
+ Solving temperature system... 1 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-no_flow/solution-00000
+
+*** Timestep 1: t=0.0625 seconds
+ Solving temperature system... 17 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-no_flow/solution-00001
+
+*** Timestep 2: t=0.125 seconds
+ Solving temperature system... 16 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-no_flow/solution-00002
+
+*** Timestep 3: t=0.1875 seconds
+ Solving temperature system... 15 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-no_flow/solution-00003
+
+*** Timestep 4: t=0.25 seconds
+ Solving temperature system... 14 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-no_flow/solution-00004
+
+*** Timestep 5: t=0.3125 seconds
+ Solving temperature system... 14 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-no_flow/solution-00005
+
+*** Timestep 6: t=0.375 seconds
+ Solving temperature system... 13 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-no_flow/solution-00006
+
+*** Timestep 7: t=0.4375 seconds
+ Solving temperature system... 13 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-no_flow/solution-00007
+
+*** Timestep 8: t=0.5 seconds
+ Solving temperature system... 12 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-no_flow/solution-00008
+
+*** Timestep 9: t=0.5625 seconds
+ Solving temperature system... 12 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-no_flow/solution-00009
+
+*** Timestep 10: t=0.625 seconds
+ Solving temperature system... 12 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-no_flow/solution-00010
+
+Number of active cells: 16 (on 3 levels)
+Number of degrees of freedom: 268 (162+25+81)
+
+*** Timestep 11: t=0.6875 seconds
+ Solving temperature system... 11 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-no_flow/solution-00011
+
+*** Timestep 12: t=0.75 seconds
+ Solving temperature system... 11 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-no_flow/solution-00012
+
+*** Timestep 13: t=0.8125 seconds
+ Solving temperature system... 10 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-no_flow/solution-00013
+
+*** Timestep 14: t=0.875 seconds
+ Solving temperature system... 10 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-no_flow/solution-00014
+
+*** Timestep 15: t=0.9375 seconds
+ Solving temperature system... 10 iterations.
+ Solving Stokes system... 0 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-no_flow/solution-00015
+
+
+
++---------------------------------------------+------------+------------+
+| Total wallclock time elapsed since start | 0.323s | |
+| | | |
+| Section | no. calls | wall time | % of total |
++---------------------------------+-----------+------------+------------+
+| Assemble Stokes system | 16 | 0.0178s | 5.5% |
+| Assemble temperature system | 16 | 0.116s | 36% |
+| Build Stokes preconditioner | 2 | 0.0311s | 9.6% |
+| Build temperature preconditioner| 16 | 0.0116s | 3.6% |
+| Solve Stokes system | 16 | 0.00568s | 1.8% |
+| Solve temperature system | 16 | 0.00687s | 2.1% |
+| Initialization | 2 | 0.00787s | 2.4% |
+| Postprocessing | 16 | 0.0798s | 25% |
+| Refine mesh structure, part 1 | 1 | 0.00386s | 1.2% |
+| Refine mesh structure, part 2 | 1 | 0.000103s | 0.032% |
+| Setup dof systems | 2 | 0.0208s | 6.4% |
++---------------------------------+-----------+------------+------------+
+
Added: trunk/aspect/tests/no_flow/statistics
===================================================================
--- trunk/aspect/tests/no_flow/statistics (rev 0)
+++ trunk/aspect/tests/no_flow/statistics 2012-11-10 07:13:25 UTC (rev 1357)
@@ -0,0 +1,25 @@
+# 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: Visualization file name
+ 0 0.0000e+00 16 187 81 1 0 6.2500e-02 output-no_flow/solution-00000
+ 1 6.2500e-02 16 187 81 17 0 6.2500e-02 output-no_flow/solution-00001
+ 2 1.2500e-01 16 187 81 16 0 6.2500e-02 output-no_flow/solution-00002
+ 3 1.8750e-01 16 187 81 15 0 6.2500e-02 output-no_flow/solution-00003
+ 4 2.5000e-01 16 187 81 14 0 6.2500e-02 output-no_flow/solution-00004
+ 5 3.1250e-01 16 187 81 14 0 6.2500e-02 output-no_flow/solution-00005
+ 6 3.7500e-01 16 187 81 13 0 6.2500e-02 output-no_flow/solution-00006
+ 7 4.3750e-01 16 187 81 13 0 6.2500e-02 output-no_flow/solution-00007
+ 8 5.0000e-01 16 187 81 12 0 6.2500e-02 output-no_flow/solution-00008
+ 9 5.6250e-01 16 187 81 12 0 6.2500e-02 output-no_flow/solution-00009
+10 6.2500e-01 16 187 81 12 0 6.2500e-02 output-no_flow/solution-00010
+11 6.8750e-01 16 187 81 11 0 6.2500e-02 output-no_flow/solution-00011
+12 7.5000e-01 16 187 81 11 0 6.2500e-02 output-no_flow/solution-00012
+13 8.1250e-01 16 187 81 10 0 6.2500e-02 output-no_flow/solution-00013
+14 8.7500e-01 16 187 81 10 0 6.2500e-02 output-no_flow/solution-00014
+15 9.3750e-01 16 187 81 10 0 6.2500e-02 output-no_flow/solution-00015
Added: trunk/aspect/tests/no_flow.prm
===================================================================
--- trunk/aspect/tests/no_flow.prm (rev 0)
+++ trunk/aspect/tests/no_flow.prm 2012-11-10 07:13:25 UTC (rev 1357)
@@ -0,0 +1,127 @@
+# A test that deals with the problem that if there is no flow then the
+# conduction time step is infinite. in that case, we need to limit the
+# time step size to something related to the conduction time step
+#
+# We make sure that there is no flow in this particular situation by using
+# zero boundary conditions and a zero thermal expansion coefficient so
+# that the right hand side of the velocity is constant.
+
+set Dimension = 2
+
+
+set CFL number = 1.0
+
+set Use conduction timestep = true
+
+set End time = 1
+
+set Output directory = output-no_flow
+
+set Resume computation = false
+
+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 Discretization
+ set Stokes velocity polynomial degree = 2
+
+ set Temperature polynomial degree = 2
+
+ set Use locally conservative discretization = false
+
+ subsection Stabilization parameters
+ set alpha = 2
+
+ set beta = 0.078
+
+ set cR = 0.5 # default: 0.11
+ end
+
+end
+
+
+subsection Geometry model
+ set Model name = box
+
+ subsection Box
+ set X extent = 1
+
+ set Y extent = 1
+
+ set Z extent = 1
+ end
+end
+
+
+subsection Gravity model
+ set Model name = vertical
+
+end
+
+
+subsection Initial conditions
+ set Model name = perturbed box
+
+end
+
+
+subsection Material model
+ set Model name = simple
+
+ subsection Simple model
+ set Thermal expansion coefficient = 0
+ set Reference specific heat = 1
+ set Thermal conductivity = 1
+ set Reference density = 1
+ end
+end
+
+
+subsection Mesh refinement
+ set Initial adaptive refinement = 0 # default: 2
+ set Initial global refinement = 2 # default: 2
+
+ set Strategy = Normalized density and temperature
+end
+
+
+subsection Model settings
+ set Include adiabatic heating = false
+ set Include shear heating = false # default: true
+
+ set Radiogenic heating rate = 0
+
+ set Fixed temperature boundary indicators = 0, 1
+
+ set Prescribed velocity boundary indicators =
+
+ set Tangential velocity boundary indicators = 0,1,2,3
+
+ set Zero velocity boundary indicators =
+end
+
+
+subsection Postprocess
+ set List of postprocessors = visualization
+
+ subsection Visualization
+ set Number of grouped files = 0
+
+ set Output format = gnuplot
+
+ set Time between graphical output = 0 # default: 1e8
+ end
+end
More information about the CIG-COMMITS
mailing list