[cig-commits] commit 2051 by bangerth to /var/svn/dealii/aspect

[dealii.demon at gmail.com]

Revision 2051

Add a new test to verify the 'topography' refinement.

A   trunk/aspect/tests/refinement_topography/
U   trunk/aspect/tests/refinement_topography/screen-output
D   trunk/aspect/tests/refinement_topography/statistics
A   trunk/aspect/tests/refinement_topography.prm


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

Diff:
Modified: trunk/aspect/tests/refinement_topography/screen-output
===================================================================
--- trunk/aspect/tests/sol_cx_2/screen-output	2013-11-25 14:06:39 UTC (rev 2050)
+++ trunk/aspect/tests/refinement_topography/screen-output	2013-11-25 14:14:05 UTC (rev 2051)
@@ -7,8 +7,23 @@
    Rebuilding Stokes preconditioner...
    Solving Stokes system... 17 iterations.
 
+Number of active cells: 22 (on 4 levels)
+Number of degrees of freedom: 392 (238+35+119)
+
+*** Timestep 0:  t=0 seconds
+   Solving temperature system... 0 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 28 iterations.
+
+Number of active cells: 46 (on 5 levels)
+Number of degrees of freedom: 791 (482+68+241)
+
+*** Timestep 0:  t=0 seconds
+   Solving temperature system... 0 iterations.
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system... 29 iterations.
+
    Postprocessing:
-     Errors u_L1, p_L1, u_L2, p_L2: 7.964441e-05, 1.174943e-01, 1.069688e-04, 1.192654e-01
 
 Termination requested by criterion: end time
 

Copied: trunk/aspect/tests/refinement_topography.prm (from rev 2050, trunk/aspect/tests/sol_cx_2.prm)
===================================================================
--- trunk/aspect/tests/refinement_topography.prm	                        (rev 0)
+++ trunk/aspect/tests/refinement_topography.prm	2013-11-25 14:14:05 UTC (rev 2051)
@@ -0,0 +1,109 @@
+# Test that we indeed refine towards the surface when using the
+# 'topography' mesh refinement criterion
+#
+# As far as output is concerned, all we care about is the number
+# of cells, which is always printed by default.
+
+set Dimension = 2
+
+
+set CFL number                             = 1.0
+
+set End time                               = 0
+
+set Output directory                       = output-refinement_topography
+
+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 = SolCx
+end
+
+
+subsection Mesh refinement
+  set Initial adaptive refinement        = 2
+  set Initial global refinement          = 2
+
+  set Strategy                           = topography
+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 = 
+end
+