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

[dealii.demon at gmail.com]

Revision 2024

New test. Does not currently work.

A   trunk/aspect/tests/compressibility_iterated_stokes/
A   trunk/aspect/tests/compressibility_iterated_stokes.cc
A   trunk/aspect/tests/compressibility_iterated_stokes.prm


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

Diff:
Copied: trunk/aspect/tests/compressibility_iterated_stokes.cc (from rev 2022, trunk/aspect/tests/compressibility.cc)
===================================================================
--- trunk/aspect/tests/compressibility_iterated_stokes.cc	                        (rev 0)
+++ trunk/aspect/tests/compressibility_iterated_stokes.cc	2013-11-20 23:11:04 UTC (rev 2024)
@@ -0,0 +1,3 @@
+// use the same postprocessing facilities as for the 'compressibility'
+// testcase
+#include "compressibility.cc"

Added: trunk/aspect/tests/compressibility_iterated_stokes.prm
===================================================================
--- trunk/aspect/tests/compressibility_iterated_stokes.prm	                        (rev 0)
+++ trunk/aspect/tests/compressibility_iterated_stokes.prm	2013-11-20 23:11:04 UTC (rev 2024)
@@ -0,0 +1,108 @@
+# A variation on the 'compressibility' testcase. here, we compute not
+# only the first iteration of a nonlinear fixed point scheme, but in
+# fact the first two. The second iteration should have a better guess
+# u^* for the continuity equation,
+#    div u = -1/rho drho/dp u^* . g
+# and so should lead to an upward expansion of the flux.
+
+set Dimension = 2
+set CFL number                             = 1.0
+set End time                               = 0
+set Output directory                       = output-compressibility_iterated_stokes
+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                = Stokes only
+set Max nonlinear iterations               = 2
+
+set Additional shared libraries = ./libcompressibility_iterated_stokes.so
+
+
+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
+    set Y extent = 1
+    set Z extent = 1
+  end
+end
+
+
+# temperature field doesn't matter. set it to zero
+subsection Initial conditions
+  set Model name = function
+  subsection Function
+    set Function expression = 0
+  end
+end
+
+
+# no gravity. the pressure will equal just the dynamic component
+subsection Gravity model
+  set Model name = vertical
+  subsection Vertical
+    set Magnitude = 1.0
+  end
+end
+
+
+subsection Material model
+  set Model name = compressibility
+
+  subsection Simple model
+    set Reference density             = 1    # default: 3300
+    set Reference specific heat       = 1250
+    set Reference temperature         = 0    # default: 293
+    set Thermal conductivity          = 1e-6 # default: 4.7
+    set Thermal expansion coefficient = 0
+    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 Tangential velocity boundary indicators = 0, 1
+  set Zero velocity boundary indicators       =
+  set Prescribed velocity boundary indicators = 2: function
+  set Include shear heating = false
+end
+
+subsection Boundary velocity model
+  subsection Function
+    set Variable names = x,y
+    set Function expression = 0;1
+  end
+end
+
+subsection Postprocess
+  set List of postprocessors = compressibility, visualization
+
+  subsection Visualization
+
+    set List of output variables      = density
+    set Number of grouped files       = 0
+    set Output format                 = vtu
+    set Time between graphical output = 0                                                                                # default: 1e8
+  end
+end