[cig-commits] r1303 - in trunk/aspect/tests: . stokes
bangerth at dealii.org
bangerth at dealii.org
Sun Oct 21 14:53:34 PDT 2012
Author: bangerth
Date: 2012-10-21 15:53:34 -0600 (Sun, 21 Oct 2012)
New Revision: 1303
Added:
trunk/aspect/tests/stokes.prm
trunk/aspect/tests/stokes/
trunk/aspect/tests/stokes/screen-output
trunk/aspect/tests/stokes/solution-00000.pvtu
trunk/aspect/tests/stokes/solution-00000.visit
trunk/aspect/tests/stokes/solution.pvd
trunk/aspect/tests/stokes/statistics
Log:
Add another test.
Added: trunk/aspect/tests/stokes/screen-output
===================================================================
--- trunk/aspect/tests/stokes/screen-output (rev 0)
+++ trunk/aspect/tests/stokes/screen-output 2012-10-21 21:53:34 UTC (rev 1303)
@@ -0,0 +1,67 @@
+Number of active cells: 64 (on 3 levels)
+Number of degrees of freedom: 3,770 (2,187+125+729+729)
+
+*** Timestep 0: t=0 seconds
+ Solving temperature system... 0 iterations.
+ Solving composition system 1... 0 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 30+9 iterations.
+
+Number of active cells: 120 (on 4 levels)
+Number of degrees of freedom: 7,248 (4,215+223+1,405+1,405)
+
+*** Timestep 0: t=0 seconds
+ Solving temperature system... 0 iterations.
+ Solving composition system 1... 0 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 30+10 iterations.
+
+Number of active cells: 120 (on 4 levels)
+Number of degrees of freedom: 7,248 (4,215+223+1,405+1,405)
+
+*** Timestep 0: t=0 seconds
+ Solving temperature system... 0 iterations.
+ Solving composition system 1... 0 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 30+10 iterations.
+
+ Postprocessing:
+
+ Reference density (kg/m^3): 3300
+ Reference gravity (m/s^2): 9.81
+ Reference thermal expansion (1/K): 2e-05
+ Temperature contrast accross model domain (K): 1
+ Model domain depth (m): 2.89e+06
+ Reference thermal diffusivity (m^2/s): 1.13939e-06
+ Reference viscosity (Pas): 1e+22
+ Ra number: 1371.62
+ k_value: 4.7
+ reference_cp: 1250
+ reference_thermal_diffusivity: 1.13939e-06
+
+ Writing graphical output: output-stokes/solution-00000
+ RMS, max velocity: 5.46e-11 m/s, 7.7e-10 m/s
+
+
+
++---------------------------------------------+------------+------------+
+| Total wallclock time elapsed since start | 6.03s | |
+| | | |
+| Section | no. calls | wall time | % of total |
++---------------------------------+-----------+------------+------------+
+| Assemble Stokes system | 3 | 0.649s | 11% |
+| Assemble composition system | 3 | 1.36s | 23% |
+| Assemble temperature system | 3 | 1.72s | 29% |
+| Build Stokes preconditioner | 3 | 0.379s | 6.3% |
+| Build composition preconditioner| 3 | 0.0403s | 0.67% |
+| Build temperature preconditioner| 3 | 0.0404s | 0.67% |
+| Solve Stokes system | 3 | 0.599s | 9.9% |
+| Solve composition system | 3 | 0.00294s | 0.049% |
+| Solve temperature system | 3 | 0.00287s | 0.048% |
+| Initialization | 4 | 0.0752s | 1.2% |
+| Postprocessing | 1 | 0.321s | 5.3% |
+| Refine mesh structure, part 1 | 2 | 0.311s | 5.2% |
+| Refine mesh structure, part 2 | 2 | 0.00497s | 0.082% |
+| Setup dof systems | 3 | 0.443s | 7.3% |
++---------------------------------+-----------+------------+------------+
+
Added: trunk/aspect/tests/stokes/solution-00000.pvtu
===================================================================
(Binary files differ)
Property changes on: trunk/aspect/tests/stokes/solution-00000.pvtu
___________________________________________________________________
Added: svn:mime-type
+ application/xml
Added: trunk/aspect/tests/stokes/solution-00000.visit
===================================================================
--- trunk/aspect/tests/stokes/solution-00000.visit (rev 0)
+++ trunk/aspect/tests/stokes/solution-00000.visit 2012-10-21 21:53:34 UTC (rev 1303)
@@ -0,0 +1,2 @@
+!NBLOCKS 1
+solution-00000.0000.vtu
Added: trunk/aspect/tests/stokes/solution.pvd
===================================================================
(Binary files differ)
Property changes on: trunk/aspect/tests/stokes/solution.pvd
___________________________________________________________________
Added: svn:mime-type
+ application/xml
Added: trunk/aspect/tests/stokes/statistics
===================================================================
--- trunk/aspect/tests/stokes/statistics (rev 0)
+++ trunk/aspect/tests/stokes/statistics 2012-10-21 21:53:34 UTC (rev 1303)
@@ -0,0 +1,16 @@
+# 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
+# 9: Iterations for Stokes solver
+# 10: Time step size (seconds)
+# 11: Visualization file name
+# 12: RMS velocity (m/s)
+# 13: Max. velocity (m/s)
+0 0.0000e+00 64 2312 729 729 0 0 39 4.0417e+14 "" 0.00000000e+00 0.00000000e+00
+0 0.0000e+00 120 4438 1405 1405 0 0 40 1.8789e+14 "" 0.00000000e+00 0.00000000e+00
+0 0.0000e+00 120 4438 1405 1405 0 0 40 1.8789e+14 output-stokes/solution-00000 5.45841472e-11 7.70027618e-10
Added: trunk/aspect/tests/stokes.prm
===================================================================
--- trunk/aspect/tests/stokes.prm (rev 0)
+++ trunk/aspect/tests/stokes.prm 2012-10-21 21:53:34 UTC (rev 1303)
@@ -0,0 +1,143 @@
+# this is a version of the stokes' law input file
+# but with less mesh refinement
+
+
+
+############### Global parameters
+# We use a 3d setup. Since we are only interested
+# in a steady state solution, we set the end time
+# equal to the start time to force a single time
+# step before the program terminates.
+
+set Dimension = 3
+
+set Start time = 0
+set End time = 0
+set Use years in output instead of seconds = false
+
+set Output directory = output-stokes
+
+
+############### Parameters describing the model
+# The setup is a 3d box with edge length 2890000 in which
+# all 6 sides have free slip boundary conditions. Because
+# the temperature plays no role in this model we need not
+# bother to describe temperature boundary conditions or
+# the material parameters that pertain to the temperature.
+
+
+subsection Geometry model
+ set Model name = box
+
+ subsection Box
+ set X extent = 2890000
+ set Y extent = 2890000
+ set Z extent = 2890000
+ end
+end
+
+
+subsection Model settings
+ set Tangential velocity boundary indicators = 0,1,2,3,4,5
+end
+
+
+subsection Material model
+ set Model name = simple
+
+ subsection Simple model
+ set Reference density = 3300
+ set Viscosity = 1e22
+ end
+end
+
+
+subsection Gravity model
+ set Model name = vertical
+
+ subsection Vertical
+ set Magnitude = 9.81
+ end
+end
+
+
+############### Parameters describing the temperature field
+# As above, there is no need to set anything for the
+# temperature boundary conditions.
+
+subsection Boundary temperature model
+ set Model name = box
+end
+
+subsection Initial conditions
+ set Model name = function
+
+ subsection Function
+ set Function expression = 0
+ end
+end
+
+############### Parameters describing the compositional field
+# This, however, is the more important part: We need to describe
+# the compositional field and its influence on the density
+# function. The following blocks say that we want to
+# advect a single compositional field and that we give it an
+# initial value that is zero outside a sphere of radius
+# r=200000m and centered at the point (p,p,p) with
+# p=1445000 (which is half the diameter of the box) and one inside.
+# The last block re-opens the material model and sets the
+# density differential per unit change in compositional field to
+# 100.
+
+subsection Compositional fields
+ set Number of fields = 1
+end
+
+subsection Compositional initial conditions
+ set Model name = function
+
+ subsection Function
+ set Variable names = x,y,z
+ set Function constants = r=200000,p=1445000
+ set Function expression = if(sqrt((x-p)*(x-p)+(y-p)*(y-p)+(z-p)*(z-p)) > r, 0, 1)
+ end
+end
+
+subsection Material model
+ subsection Simple model
+ set Density differential for compositional field 1 = 100
+ end
+end
+
+
+
+
+############### Parameters describing the discretization
+# The following parameters describe how often we want to refine
+# the mesh globally and adaptively, what fraction of cells should
+# be refined in each adaptive refinement step, and what refinement
+# indicator to use when refining the mesh adaptively.
+
+subsection Mesh refinement
+ set Initial adaptive refinement = 2
+ set Initial global refinement = 2
+ set Refinement fraction = 0.2
+ set Strategy = Velocity
+end
+
+
+############### Parameters describing the what to do with the solution
+# The final section allows us to choose which postprocessors to
+# run at the end of each time step. We select to generate graphical
+# output that will consist of the primary variables (velocity, pressure,
+# temperature and the compositional fields) as well as the density and
+# viscosity. We also select to compute some statistics about the
+# velocity field.
+
+subsection Postprocess
+ set List of postprocessors = visualization, velocity statistics
+
+ subsection Visualization
+ set List of output variables = density, viscosity
+ end
+end
More information about the CIG-COMMITS
mailing list