[cig-commits] [commit] master: Add another test that ensures the simple compressible material model gives the same result as the existing compressibility_iterated_stokes test. (1cdd9a1)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/aspect

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/4f4b96e0569b8be211d2d72b283391b5d08b4e1d...0586ab6f9cf6ec525c331a8c0113ca155a3f1f87

>---------------------------------------------------------------

commit 1cdd9a19532f515e68024ead8e822911d470275d
Author: Rene Gassmoeller <R.Gassmoeller at mailbox.org>
Date:   Tue May 20 15:01:18 2014 -0500

    Add another test that ensures the simple compressible material model gives the same result as the existing compressibility_iterated_stokes test.


>---------------------------------------------------------------

1cdd9a19532f515e68024ead8e822911d470275d
 tests/simple_compressibility_iterated_stokes.cc    |  2 ++
 ... => simple_compressibility_iterated_stokes.prm} | 33 +++++++++++-----------
 .../screen-output                                  |  2 +-
 3 files changed, 19 insertions(+), 18 deletions(-)

diff --git a/tests/simple_compressibility_iterated_stokes.cc b/tests/simple_compressibility_iterated_stokes.cc
new file mode 100644
index 0000000..82b6adc
--- /dev/null
+++ b/tests/simple_compressibility_iterated_stokes.cc
@@ -0,0 +1,2 @@
+#include "compressibility.cc"
+
diff --git a/tests/no_adiabatic_heating.prm b/tests/simple_compressibility_iterated_stokes.prm
similarity index 67%
copy from tests/no_adiabatic_heating.prm
copy to tests/simple_compressibility_iterated_stokes.prm
index 8a0d7a3..c45cb24 100644
--- a/tests/no_adiabatic_heating.prm
+++ b/tests/simple_compressibility_iterated_stokes.prm
@@ -1,15 +1,11 @@
-# This is a follow-on on the 'compressibility_iterated_stokes' testcase
-# where we choose
-#   u_z = 1+z/100
-# If we switch off adiabatic heating and friction heating, then there
-# should be no effect on the temperature field from the velocity, and
-# no effect from the temperature onto the velocity either. Verify that
-# this is so
+# A variation on the 'compressibility_iterated_stokes' testcase. 
+# The only variation is the used material model. This test adapts the
+# 'simple compressible' material model plugin to do the same job as the 
+# shared library material model from compressibility_iterated_stokes.cc
 
 set Dimension = 2
-
 set CFL number                             = 1.0
-set End time                               = 1
+set End time                               = 0
 set Start time                             = 0
 set Adiabatic surface temperature          = 0
 set Surface pressure                       = 0
@@ -61,14 +57,14 @@ end
 
 
 subsection Material model
-  set Model name = no adiabatic heating
+  set Model name = simple compressible
 
-  subsection Simple model
-    set Reference density             = 1    # default: 3300
+  subsection Simple compressible model
+    set Reference density             = 0.909090909090    # default: 3300
+    set Reference compressibility     = 0.01
     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
@@ -76,7 +72,7 @@ end
 
 subsection Mesh refinement
   set Initial adaptive refinement        = 0
-  set Initial global refinement          = 4
+  set Initial global refinement          = 5
 end
 
 
@@ -86,7 +82,6 @@ subsection Model settings
   set Zero velocity boundary indicators       =
   set Prescribed velocity boundary indicators = 2: function
   set Include shear heating = false
-  set Include adiabatic heating = false
 end
 
 subsection Boundary velocity model
@@ -97,9 +92,13 @@ subsection Boundary velocity model
 end
 
 subsection Postprocess
-  set List of postprocessors = temperature statistics
+  set List of postprocessors = compressibility
 
   subsection Visualization
-    set Time between graphical output = 0.1
+
+    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
diff --git a/tests/compressibility_iterated_stokes/screen-output b/tests/simple_compressibility_iterated_stokes/screen-output
similarity index 94%
copy from tests/compressibility_iterated_stokes/screen-output
copy to tests/simple_compressibility_iterated_stokes/screen-output
index b1514b7..94fcb2c 100644
--- a/tests/compressibility_iterated_stokes/screen-output
+++ b/tests/simple_compressibility_iterated_stokes/screen-output
@@ -5,7 +5,7 @@
 --     . using Trilinos
 -----------------------------------------------------------------------------
 
-Loading shared library <./libcompressibility_iterated_stokes.so>
+Loading shared library <./libsimple_compressibility_iterated_stokes.so>
 
 Number of active cells: 1,024 (on 6 levels)
 Number of degrees of freedom: 13,764 (8,450+1,089+4,225)