Update stress build-up benchmark

geodynamics · Aug 21, 2017 · c41bf2c · c41bf2c
1 parent 2643312
commit c41bf2c
Showing 1 changed file with 21 additions and 24 deletions.
diff --git a/benchmarks/stress_build-up_maxwell_2d/stress_build-up_maxwell_2d.prm b/benchmarks/stress_build-up_maxwell_2d/stress_build-up_maxwell_2d.prm
@@ -1,30 +1,28 @@
 #  Global parameters
 set Dimension                              = 2
 set Start time                             = 0
-set End time                               = 50e2
+set End time                               = 1e6
 set Use years in output instead of seconds = true
 set Use direct solver for Stokes system    = false
 set Linear solver tolerance                = 1e-7
 set Nonlinear solver scheme                = IMPES
 set Nonlinear solver tolerance             = 1e-4
 set Max nonlinear iterations               = 10
-set Number of cheap Stokes solver steps    = 200
-set CFL number                             = 0.1
-set Maximum time step                      = 100
+set Number of cheap Stokes solver steps    = 500
+set CFL number                             = 0.5
+set Maximum time step                      = 1000
 set Output directory                       = output
 set Timing output frequency                = 1
 set Pressure normalization                 = surface
 set Surface pressure                       = 0.
 
-#### Parameters describing the model
-
-# Model geometry (100x100 km, 5 km spacing)
+# Model geometry (100x100 km, 1 km spacing)
 subsection Geometry model
   set Model name = box
 
   subsection Box
-    set X repetitions = 20
-    set Y repetitions = 20
+    set X repetitions = 100
+    set Y repetitions = 100
     set X extent      = 100e3
     set Y extent      = 100e3
   end
@@ -39,12 +37,11 @@ end
 
 # Boundary classifications
 subsection Model settings
-  set Include adiabatic heating               = false
-  set Include shear heating                   = false
-  set Fixed temperature boundary indicators   = bottom, top
-  set Prescribed velocity boundary indicators = bottom y:function, top y:function, left x:function, right x:function 
+  set Fixed composition boundary indicators   = bottom, top, left, right
+  set Fixed temperature boundary indicators   = bottom, top, left, right
+  set Prescribed velocity boundary indicators = bottom y:function, right x:function 
   set Free surface boundary indicators        = 
-  set Tangential velocity boundary indicators = 
+  set Tangential velocity boundary indicators = top, left 
   set Zero velocity boundary indicators       = 
   set Enable elasticity                       = true
 end
@@ -53,7 +50,7 @@ subsection Boundary velocity model
   subsection Function
     set Variable names      = x,y
     set Function constants  = cm=0.01, year=1
-    set Function expression = if (x<50e3 , -0.5*cm/year, 0.5*cm/year) ; if (y<50e3 , 0.5*cm/year, -0.5*cm/year) ; 
+    set Function expression = if (x<50e3 , -0.3154*cm/year, 0.); if (y<50e3 , -3.154*cm/year, 0.); 
   end
 end
 
@@ -74,6 +71,7 @@ subsection Initial composition model
 end
 
 # Composition boundary conditions
+# (Note: should change boundary values to be defined using analytical solution)
 subsection Boundary composition model
   set Model name = initial composition
 end
@@ -102,9 +100,9 @@ subsection Material model
 
   subsection Viscoelastic
     set Densities            =  2800
-    set Viscosities          = 1.e21
+    set Viscosities          = 1.e22
     set Elastic shear moduli = 1.e10
-    set Elastic time step = 1e2
+    set Elastic time step = 1e3
     set Viscosity averaging scheme = harmonic
   end
 
@@ -121,31 +119,30 @@ end
 # Element types
 subsection Discretization
   set Composition polynomial degree     = 2
-  set Stokes velocity polynomial degree = 3
+  set Stokes velocity polynomial degree = 2
   set Temperature polynomial degree     = 1
 end
 
 # Post processing
 subsection Postprocess
-  set List of postprocessors = velocity statistics, basic statistics, temperature statistics, visualization
-  set Run postprocessors on nonlinear iterations = true
+  set List of postprocessors = basic statistics, composition statistics, temperature statistics, velocity statistics, visualization
   subsection Visualization
     set List of output variables = material properties, strain rate 
 
     subsection Material properties
       set List of material properties = density, viscosity
     end
 
-    set Time between graphical output = 0
-    set Interpolate output = false
+    set Time between graphical output = 1e3
+    set Interpolate output = true
   end
 
 end
 
 # Termination criteria
 subsection Termination criteria
- set End step = 50
- set Termination criteria = end step
+ set End step = 1000
+ set Termination criteria = end time
 end