[cig-commits] r18994 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205

[brad at geodynamics.org]

Author: brad
Date: 2011-09-29 12:48:51 -0700 (Thu, 29 Sep 2011)
New Revision: 18994

Removed:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells_hex8.jou
Modified:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/createbc.jou
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m.jou
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m.jou
Log:
Updated mesh generation.

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/createbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/createbc.jou	2011-09-29 18:35:57 UTC (rev 18993)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/createbc.jou	2011-09-29 19:48:51 UTC (rev 18994)
@@ -61,9 +61,10 @@
 # Create nodeset for subset of +z face
 # ----------------------------------------------------------------------
 group "face_zpos_subset" add node in group face_zpos
-group "face_zpos_subset" remove node with x_coord > +15.001e+3
-group "face_zpos_subset" remove node with x_coord < -15.001e+3
-group "face_zpos_subset" remove node with y_coord > +5.001e+3
-group "face_zpos_subset" remove node with y_coord < -5.001e+3
+group "face_zpos_subset" remove node with y_coord > +15.001e+3
+group "face_zpos_subset" remove node with y_coord < -15.001e+3
+group "face_zpos_subset" remove node with x_coord > +10.001e+3
+group "face_zpos_subset" remove node with x_coord < -10.001e+3
 nodeset 17 group face_zpos_subset
 nodeset 17 name "face zpos_subset"
+

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m.jou	2011-09-29 18:35:57 UTC (rev 18993)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/hex8_200m.jou	2011-09-29 19:48:51 UTC (rev 18994)
@@ -8,9 +8,11 @@
 playback 'geometry.jou'
 
 # ----------------------------------------------------------------------
-# Set discretization size
+# Set discretization size and scheme
 # ----------------------------------------------------------------------
-volume all size 200
+#{dx=200*m}
+volume all scheme map
+volume all size {dx*3}
 
 # ----------------------------------------------------------------------
 # Generate the mesh
@@ -18,6 +20,27 @@
 mesh volume all
 
 # ----------------------------------------------------------------------
+# Refine around fault
+# ----------------------------------------------------------------------
+#{refineWidth=18.001*km}
+#{refineLength=45.001*km}
+#{refineDepth=21.001*km}
+group "refinecells" add Hex with x_min >= {-0.5*refineWidth} and x_max <= {+0.5*refineWidth} and y_min >= {-0.5*refineLength} and y_max <= {+0.5*refineLength} and z_min >= {-refineDepth}
+refine hex in group refinecells depth 1
+delete group refinecells
+
+# ----------------------------------------------------------------------
+# Smooth to improve quality
+# ----------------------------------------------------------------------
+#{condnum=3.0}
+#{loop(5)}
+volume all smooth scheme condition number beta {condnum} cpu 10
+smooth volume all
+#{condnum=condnum-0.2}
+#{endloop}
+
+
+# ----------------------------------------------------------------------
 # Setup boundary conditions.
 # ----------------------------------------------------------------------
 playback 'createbc.jou'
@@ -28,3 +51,5 @@
 export mesh "hex8_200m.exo" dimension 3 overwrite
 
 
+
+

Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells_hex8.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells_hex8.jou	2011-09-29 18:35:57 UTC (rev 18993)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/refinecells_hex8.jou	2011-09-29 19:48:51 UTC (rev 18994)
@@ -1,8 +0,0 @@
-# ----------------------------------------------------------------------
-# Refine elements within 1/2 of rupture width
-# ----------------------------------------------------------------------
-group "refinecells1" add hex with x_min >= -36.001e+3 and x_max <= 36.001e+3 and y_min >= -24.001e+3 and y_max <= 24.001e+3 and z_min >= -36.001e+3
-refine hex in refinecells1 depth 1
-group "refinecells2" add hex with x_min >= -24.001e+3 and x_max <= 24.001e+3 and y_min >= -12.001e+3 and y_max <= 12.001e+3 and z_min >= -24.001e+3
-refine hex in refinecells2 depth 1
-

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m.jou	2011-09-29 18:35:57 UTC (rev 18993)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m.jou	2011-09-29 19:48:51 UTC (rev 18994)
@@ -47,3 +47,4 @@
 # ----------------------------------------------------------------------
 export mesh "tet4_200m.exo" dimension 3 overwrite
 
+