[cig-commits] r19995 - short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes
willic3 at geodynamics.org
willic3 at geodynamics.org
Wed Apr 25 16:16:53 PDT 2012
Author: willic3
Date: 2012-04-25 16:16:53 -0700 (Wed, 25 Apr 2012)
New Revision: 19995
Added:
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_20km.exo.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_20km.jou
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_6.7km.exo.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_6.7km.jou
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_6.7km.jou
Removed:
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_refined_20km.exo.gz
short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_refined_20km.jou
Log:
Added files for coarser meshes (net yet complete).
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_20km.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_20km.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/x-gzip
Added: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_20km.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_20km.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_20km.jou 2012-04-25 23:16:53 UTC (rev 19995)
@@ -0,0 +1,27 @@
+# The volume is meshed at an initial resolution of 20 km.
+#
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size to 20 km.
+# ----------------------------------------------------------------------
+#{dx=20.0*km}
+volume all size {dx}
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "hex8_20km.exo" dimension 3 overwrite
+
+
+
+
Copied: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_6.7km.exo.gz (from rev 19990, short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_refined_20km.exo.gz)
===================================================================
(Binary files differ)
Copied: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_6.7km.jou (from rev 19990, short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_refined_20km.jou)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_6.7km.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_6.7km.jou 2012-04-25 23:16:53 UTC (rev 19995)
@@ -0,0 +1,40 @@
+# The volume is meshed at an initial resolution of 20 km.
+#
+# An inner volume with dimensions of 480 km (x), 240 km (y), and 100 km (z)
+# is then put through 1 level of refinement.
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size to 20 km.
+# ----------------------------------------------------------------------
+#{dx=20.0*km}
+volume all size {dx}
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Refine region inside inner box.
+# ----------------------------------------------------------------------
+group "inner_box" add hex with x_coord > {-240.0*km} and with x_coord < {240.0*km} and with y_coord > {-120.0*km} and with y_coord < {120.0*km} and with z_coord > {-100.0*km}
+refine hex in inner_box depth 0
+
+# ----------------------------------------------------------------------
+# Smooth mesh.
+# ----------------------------------------------------------------------
+volume all smooth scheme condition number beta 2.2
+smooth volume all
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "hex8_refined_20km.exo" dimension 3 overwrite
+
+
+
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_refined_20km.exo.gz
===================================================================
(Binary files differ)
Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_refined_20km.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_refined_20km.jou 2012-04-25 23:14:41 UTC (rev 19994)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/hex8_refined_20km.jou 2012-04-25 23:16:53 UTC (rev 19995)
@@ -1,40 +0,0 @@
-# The volume is meshed at an initial resolution of 20 km.
-#
-# An inner volume with dimensions of 480 km (x), 240 km (y), and 100 km (z)
-# is then put through 1 level of refinement.
-# ----------------------------------------------------------------------
-# Generate geometry
-# ----------------------------------------------------------------------
-playback 'geometry.jou'
-
-# ----------------------------------------------------------------------
-# Set discretization size to 20 km.
-# ----------------------------------------------------------------------
-#{dx=20.0*km}
-volume all size {dx}
-mesh volume all
-
-# ----------------------------------------------------------------------
-# Refine region inside inner box.
-# ----------------------------------------------------------------------
-group "inner_box" add hex with x_coord > {-240.0*km} and with x_coord < {240.0*km} and with y_coord > {-120.0*km} and with y_coord < {120.0*km} and with z_coord > {-100.0*km}
-refine hex in inner_box depth 0
-
-# ----------------------------------------------------------------------
-# Smooth mesh.
-# ----------------------------------------------------------------------
-volume all smooth scheme condition number beta 2.2
-smooth volume all
-
-# ----------------------------------------------------------------------
-# Generate geometry
-# ----------------------------------------------------------------------
-playback 'createbc.jou'
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "hex8_refined_20km.exo" dimension 3 overwrite
-
-
-
Copied: short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_6.7km.jou (from rev 19990, short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_graded_6.7km.jou)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_6.7km.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/sceccrustdeform/savageprescott/meshes/tet4_6.7km.jou 2012-04-25 23:16:53 UTC (rev 19995)
@@ -0,0 +1,188 @@
+# -*- Python -*- (syntax highlighting)
+# The volume is meshed at a variable resolution, with a maximum resolution
+# of 6.7 km at the upper center of the mesh, and a minimum resolution of
+# 20 km at the outer x-edges of the mesh.
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size and scheme
+# ----------------------------------------------------------------------
+#{dx=20.0*km/3.0}
+
+# ----------------------------------------------------------------------
+# Set vertex spacing with increasing spacing away from fault
+# ----------------------------------------------------------------------
+#{bias_factor=1.02}
+
+# Function for computing discretization size at end of curve with bias
+#{fun_dx='dxStart*bias_factor**ceil( ln(1-curveL/dxStart*(1-bias_factor))/ln(bias_factor))'}
+
+# ----------------------------------------------------------------------
+# Compute sizes at curve endpoints
+# ----------------------------------------------------------------------
+# dxA - size at vertex 9 28 22 12 29 21
+#{dxStart=dx}{curveL=Length(27)}{execute("dxA="//fun_dx)}
+
+# dxB - size at vertex 27
+#{dxStart=dx}{curveL=Length(41)}{execute("dxB="//fun_dx)}
+
+# dxC - size at vertex 10 11
+#{dxStart=dxB}{curveL=Length(47)}{execute("dxC="//fun_dx)}
+
+# dxD - size at vertex 28 19 23 26
+#{dxStart=dx}{curveL=Length(25)}{execute("dxD="//fun_dx)}
+
+# dxE - size at vertex 17 20
+#{dxStart=dxD}{curveL=Length(43)}{execute("dxE="//fun_dx)}
+
+# dxF - size at vertex 1 2 3 4 13 14 15 16
+#{dxStart=dxD}{curveL=Length(26)}{execute("dxF="//fun_dx)}
+
+# dxG - size at vertex 5 6 7 8
+#{dxStart=dxF}{curveL=Length(16)}{execute("dxG="//fun_dx)}
+
+# ----------------------------------------------------------------------
+# Reset sizes
+# ----------------------------------------------------------------------
+curve all scheme default
+surface all sizing function none
+surface all scheme trimesh
+
+# Fault
+curve 56 57 size {dx}
+
+# fault to A
+curve 27 scheme bias fine size {dx} coarse size {dxA} start vertex 25
+curve 31 scheme bias fine size {dx} coarse size {dxA} start vertex 25
+curve 54 scheme bias fine size {dx} coarse size {dxA} start vertex 30
+curve 58 scheme bias fine size {dx} coarse size {dxA} start vertex 30
+curve 29 scheme bias fine size {dx} coarse size {dxA} start vertex 24
+curve 30 scheme bias fine size {dx} coarse size {dxA} start vertex 24
+
+# fault to B
+curve 41 scheme bias fine size {dx} coarse size {dxB} start vertex 24
+
+# fault to D
+curve 25 scheme bias fine size {dx} coarse size {dxD} start vertex 25
+curve 35 scheme bias fine size {dx} coarse size {dxD} start vertex 25
+curve 23 scheme bias fine size {dx} coarse size {dxD} start vertex 24
+curve 34 scheme bias fine size {dx} coarse size {dxD} start vertex 24
+
+# ----------------------------------------------------------------------
+# bias starting at A
+
+# A to C
+curve 13 scheme bias fine size {dxA} coarse size {dxC} start vertex 22
+curve 15 scheme bias fine size {dxA} coarse size {dxC} start vertex 21
+
+# A to F
+curve 1 scheme bias fine size {dxA} coarse size {dxF} start vertex 9
+curve 8 scheme bias fine size {dxA} coarse size {dxF} start vertex 9
+curve 2 scheme bias fine size {dxA} coarse size {dxF} start vertex 12
+curve 6 scheme bias fine size {dxA} coarse size {dxF} start vertex 12
+
+curve 10 scheme bias fine size {dxA} coarse size {dxF} start vertex 22
+curve 17 scheme bias fine size {dxA} coarse size {dxF} start vertex 22
+curve 9 scheme bias fine size {dxA} coarse size {dxF} start vertex 21
+curve 18 scheme bias fine size {dxA} coarse size {dxF} start vertex 21
+
+# A to A
+curve 53 59 52 55 size {dxA}
+
+# ----------------------------------------------------------------------
+# bias starting at B
+
+# B to C
+curve 44 scheme bias fine size {dxB} coarse size {dxC} start vertex 27
+curve 47 scheme bias fine size {dxB} coarse size {dxC} start vertex 27
+
+# B to E
+curve 42 scheme bias fine size {dxB} coarse size {dxE} start vertex 27
+curve 49 scheme bias fine size {dxB} coarse size {dxE} start vertex 27
+
+# ----------------------------------------------------------------------
+# bias starting at C
+
+# C to G
+curve 3 scheme bias fine size {dxC} coarse size {dxG} start vertex 10
+curve 7 scheme bias fine size {dxC} coarse size {dxG} start vertex 10
+curve 4 scheme bias fine size {dxC} coarse size {dxG} start vertex 11
+curve 5 scheme bias fine size {dxC} coarse size {dxG} start vertex 11
+
+# ----------------------------------------------------------------------
+# bias starting at D
+
+curve 26 scheme bias fine size {dxD} coarse size {dxF} start vertex 19
+curve 33 scheme bias fine size {dxD} coarse size {dxF} start vertex 19
+curve 37 scheme bias fine size {dxD} coarse size {dxF} start vertex 18
+curve 40 scheme bias fine size {dxD} coarse size {dxF} start vertex 18
+
+curve 32 scheme bias fine size {dxD} coarse size {dxE} start vertex 23
+curve 28 scheme bias fine size {dxD} coarse size {dxF} start vertex 23
+curve 38 scheme bias fine size {dxD} coarse size {dxF} start vertex 26
+curve 39 scheme bias fine size {dxD} coarse size {dxF} start vertex 26
+
+curve 43 scheme bias fine size {dxD} coarse size {dxF} start vertex 23
+curve 48 scheme bias fine size {dxD} coarse size {dxE} start vertex 26
+
+# D to D
+curve 24 36 size {dxD}
+
+# ----------------------------------------------------------------------
+# bias starting at E
+
+# E to G
+curve 45 scheme bias fine size {dxE} coarse size {dxG} start vertex 20
+curve 46 scheme bias fine size {dxE} coarse size {dxG} start vertex 20
+curve 50 scheme bias fine size {dxE} coarse size {dxG} start vertex 17
+curve 51 scheme bias fine size {dxE} coarse size {dxG} start vertex 17
+
+# ----------------------------------------------------------------------
+# bias starting at F
+
+curve 14 scheme bias fine size {dxF} coarse size {dxG} start vertex 15
+curve 16 scheme bias fine size {dxF} coarse size {dxG} start vertex 16
+curve 22 scheme bias fine size {dxF} coarse size {dxG} start vertex 14
+curve 21 scheme bias fine size {dxF} coarse size {dxG} start vertex 13
+
+# F to F
+curve 12 26 19 20 size {dxF}
+
+# ----------------------------------------------------------------------
+# Generate the mesh.
+# ----------------------------------------------------------------------
+mesh surface all
+surface all smooth scheme condition number beta 1.3
+smooth surface all
+
+volume all scheme tetmesh
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Smooth mesh to improve quality
+# ----------------------------------------------------------------------
+#{condnum=2.0}
+#{loop(4)}
+cleanup volume all
+volume all smooth scheme condition number beta {condnum} cpu 10
+smooth volume all
+#{condnum=condnum-0.1}
+#{endloop}
+
+# ----------------------------------------------------------------------
+# Setup boundary conditions.
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "tet4_6.7km.exo" dimension 3 overwrite
+
+
+
+
+
More information about the CIG-COMMITS
mailing list