[cig-commits] r14382 - short/2.5D/benchmarks/savageprescott/meshes

[willic3 at geodynamics.org]

Author: willic3
Date: 2009-03-18 13:16:21 -0700 (Wed, 18 Mar 2009)
New Revision: 14382

Added:
   short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded3_20km.exo.gz
   short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded5_20km.exo.gz
   short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded5_20km.jou
Modified:
   short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded3_20km.jou
Log:
Updated journal file and added mesh for graded3 mesh.
Added journal file and mesh for graded5 mesh, which has 1 level of
refinement all along the fault.



Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded3_20km.exo.gz
===================================================================
(Binary files differ)


Property changes on: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded3_20km.exo.gz
___________________________________________________________________
Name: svn:mime-type
   + application/octet-stream

Modified: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded3_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded3_20km.jou	2009-03-18 19:39:54 UTC (rev 14381)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded3_20km.jou	2009-03-18 20:16:21 UTC (rev 14382)
@@ -21,7 +21,7 @@
 # The volume is meshed at an initial resolution of 20 km.
 # An inner volume with dimensions of 500 km (x), 250 km (y), and 100 km (z)
 # is then put through 1 level of refinement.
-# The final mesh contains 247405 vertices and 235482 cells.
+# The final mesh contains 156633 vertices and 147110 cells.
 # ----------------------------------------------------------------------
 # Create outer block
 # ----------------------------------------------------------------------

Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded5_20km.exo.gz
===================================================================
(Binary files differ)


Property changes on: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded5_20km.exo.gz
___________________________________________________________________
Name: svn:mime-type
   + application/octet-stream

Added: short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded5_20km.jou
===================================================================
--- short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded5_20km.jou	                        (rev 0)
+++ short/2.5D/benchmarks/savageprescott/meshes/spbm_hex8_graded5_20km.jou	2009-03-18 20:16:21 UTC (rev 14382)
@@ -0,0 +1,190 @@
+## /Applications/Plotting/Cubit_10_2.app/Contents/MacOS/Cubit_10_2
+## Cubit Version 10.2
+## Cubit Build 8
+## Revised 10/17/2006 15:13:41 MST
+## Running 06/16/2008 02:47:47 PM
+## Command Options:
+## -warning = On
+## -information = On
+
+# Description:
+# This journal creates a mesh with dimensions of 2000 km (x), 500 km (y),
+# and 800 km (z).  The mesh is divided into 3 horizontal layers:
+# a surface layer from 0 to 20 km, marking the bottom of the locking depth;
+# a middle layer from 20 to 40 km, marking the creeping section;
+# a bottom layer from 40 to 800 km, which approximates the viscoelastic
+# half-space.
+# The mesh is cut vertically by the fault surface, which lies along the
+# x-axis, and by a symmetry plane, which lies along the y-axis.
+# In this version of the mesh, the fault is restricted to the upper two
+# (elastic) layers.
+# The volume is meshed at an initial resolution of 20 km.
+# An inner volume with dimensions of 500 km (x), 500 km (y), and 200 km (z)
+# is then put through 1 level of refinement.
+# Note that the entire region along the fault (y-axis) is refined.
+# The final mesh contains 307051 vertices and 294194 cells.
+# ----------------------------------------------------------------------
+# Create outer block
+# ----------------------------------------------------------------------
+
+# Block is 2000 km x 500 km x 800 km
+# -1000 km <= x <= 1000 km
+# -250 km <= y <= 250 km
+# -800 km <= z <= 0 km
+reset
+brick x 2000000 y 500000 z 800000
+
+# Translate block so the top is at z=0
+volume 1 move x 0 y 0 z -400000
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+create planar surface with plane xplane offset 0
+surface 7 name "fault_surface"
+create planar surface with plane zplane offset -40000
+surface 8 name "material_interface"
+create planar surface with plane yplane offset 0
+surface 9 name "symmetry_plane"
+create planar surface with plane zplane offset -20000
+surface 10 name "fault_bottom"
+
+# ----------------------------------------------------------------------
+# Cut volume with interfaces.
+# ----------------------------------------------------------------------
+compress all
+webcut volume all with plane fault_surface
+webcut volume all with plane material_interface
+webcut volume all with plane symmetry_plane
+webcut volume all with plane fault_bottom
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces.
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
+
+# ----------------------------------------------------------------------
+# Remove surface bodies.
+# ----------------------------------------------------------------------
+delete body 2 3 4 5 7 8 9 11 12 14 16 17 18 20 21 23 26 27 29 31
+imprint all with volume all
+merge all
+
+# ----------------------------------------------------------------------
+# Set interval size for outer surfaces to 20 km
+# ----------------------------------------------------------------------
+volume all size 20000.0
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Refine region inside inner box.
+# ----------------------------------------------------------------------
+group "inner_box" add hex with x_coord > -250000.0 and with x_coord < 250000.0 and with y_coord > -300000.0 and with y_coord < 300000.0 and with z_coord > -200000.0
+refine hex in inner_box depth 0
+
+# ----------------------------------------------------------------------
+# Smooth mesh.
+# ----------------------------------------------------------------------
+set debug 91 on
+volume all smooth scheme smart laplacian
+smooth volume all
+volume all smooth scheme condition number
+smooth volume all
+volume all smooth scheme mean ratio
+smooth volume all
+volume all smooth scheme condition number
+smooth volume all
+
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 volume 1 6 15 19 25 28 30 32
+block 1 name "elastic"
+block 2 volume 10 13 22 24
+block 2 name "viscoelastic"
+
+# ----------------------------------------------------------------------
+# Create fault nodeset
+# ----------------------------------------------------------------------
+group "fault" add node in fault_surface
+group "fault" add node in fault_surface at A
+group "fault" add node in fault_surface at D
+group "fault" add node in fault_surface at E
+nodeset 10 fault
+nodeset 10 name "fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 133 109 139 111 89 87
+nodeset 11 face_xpos
+nodeset 11 name "face_xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 121 147 127 149 97 104
+nodeset 12 face_xneg
+nodeset 12 name "face_xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 106 120 113 128 35 47
+nodeset 13 face_ypos
+nodeset 13 name "face_ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 144 132 151 140 49 33
+nodeset 14 face_yneg
+nodeset 14 name "face_yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 103 92 96 84
+nodeset 15 face_zneg
+nodeset 15 name "face_zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face without fault
+# ----------------------------------------------------------------------
+group "face_zneg_nofault" add node in face_zneg
+group "face_zneg_nofault" remove node in fault
+nodeset 16 face_zneg_nofault
+nodeset 16 name "face_zneg_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 78 62 70 54
+nodeset 17 face_zpos
+nodeset 17 name "face_zpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face without fault or x-boundaries
+# ----------------------------------------------------------------------
+group "face_yneg_nofault" add node in face_yneg
+group "face_yneg_nofault" remove node in fault
+group "face_yneg_nofault" remove node in face_xpos
+group "face_yneg_nofault" remove node in face_xneg
+nodeset 18 face_yneg_nofault
+nodeset 18 name "face_yneg_nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face without fault or x-boundaries
+# ----------------------------------------------------------------------
+group "face_ypos_nofault" add node in face_ypos
+group "face_ypos_nofault" remove node in fault
+group "face_ypos_nofault" remove node in face_xpos
+group "face_ypos_nofault" remove node in face_xneg
+nodeset 19 face_ypos_nofault
+nodeset 19 name "face_ypos_nofault"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "spbm_hex8_graded5_20km.exo" dimension 3 overwrite