[cig-commits] r16736 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d

[surendra at geodynamics.org]

Author: surendra
Date: 2010-05-17 17:58:40 -0700 (Mon, 17 May 2010)
New Revision: 16736

Added:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry_refine.jou
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.jou
Modified:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry.jou
Log:
Started working on refined mesh for tpv205-2d

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry.jou	2010-05-18 00:30:19 UTC (rev 16735)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry.jou	2010-05-18 00:58:40 UTC (rev 16736)
@@ -7,7 +7,7 @@
 #
 # We create the 2-D mesh by extracting the mid-surface from a 3-D mesh.
 reset
-brick x 74000 y 48000 z 24000
+brick x 24000 y 48000 z 24000
 create midsurface volume 1 surface 1 2
 delete volume 1
 
@@ -37,3 +37,4 @@
 merge all
 
 
+

Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry_refine.jou	                        (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry_refine.jou	2010-05-18 00:58:40 UTC (rev 16736)
@@ -0,0 +1,29 @@
+# ----------------------------------------------------------------------
+# Create block
+# ----------------------------------------------------------------------
+# Block is 72.0 km x 144.0 km
+# -36.0 km <= x <= 36.0 km
+# -72.0 km <= y <= 72.0 km
+#
+# We create the 2-D mesh by extracting the mid-surface from a 3-D mesh.
+reset
+brick x 72000 y 144000 z 72000
+create midsurface volume 1 surface 1 2
+delete volume 1
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+
+# Create fault (yz plane) at x = 0.0
+split surface 7 across location position 0 -15000 0 location position 0 15000 0
+curve 17 name "fault_surface"
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces
+# ----------------------------------------------------------------------
+imprint all
+merge all
+
+
+

Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.jou	                        (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.jou	2010-05-18 00:58:40 UTC (rev 16736)
@@ -0,0 +1,41 @@
+# ----------------------------------------------------------------------
+# Create quad4 mesh at 200m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry_refine.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+surface all size 600
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh surface all
+
+# ----------------------------------------------------------------------
+# Setup boundary conditions.
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Refine elements within 1/2 of rupture width
+# ----------------------------------------------------------------------
+group "Refined"" add node in surface 8
+group "Refined"" add node in surface 9
+group "Refined" copy scale 0.5
+#nodeset 21 group Refined
+#nodeset 21 name "Refined"
+
+refine node Refined numsplit=3
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "quad4_200m_refine.exo" dimension 2 overwrite
+
+