[cig-commits] r17865 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14
brad at geodynamics.org
brad at geodynamics.org
Mon Feb 14 16:59:48 PST 2011
Author: brad
Date: 2011-02-14 16:59:48 -0800 (Mon, 14 Feb 2011)
New Revision: 17865
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/createbc.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/gradient.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/params.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tet4_200m.jou
Log:
Initial tet4 mesh.
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/createbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/createbc.jou 2011-02-15 00:41:26 UTC (rev 17864)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/createbc.jou 2011-02-15 00:59:48 UTC (rev 17865)
@@ -5,25 +5,20 @@
block 1 name "elastic"
# ----------------------------------------------------------------------
-# Create nodeset for fault
+# Create nodeset for main fault
# ----------------------------------------------------------------------
-group "main_fault" add node in fault_surface_main
-group "main_fault" add node in fault_surface_main at A
-group "main_fault" remove node with x_coord > {mainFaultLength/2.0}
-group "main_fault" remove node with x_coord < {-mainFaultLength/2.0}
-group "main_fault" remove node with z_coord < {-mainFaultDepth}
+group "main_fault" add node in fault_surface_mainA
+group "main_fault" add node in fault_surface_mainB
nodeset 100 group main_fault
nodeset 100 name "main_fault"
# ----------------------------------------------------------------------
-# Create nodeset for fault
+# Create nodeset for branch fault
# ----------------------------------------------------------------------
group "branch_fault" add node in fault_surface_branch
-group "branch_fault" remove node with x_coord > {branchFaultLength*sind(branchAngle)}
-group "branch_fault" remove node with z_coord < {-branchFaultDepth}
group "branch_fault" remove node in main_fault
-nodeset 200 group branch_fault
-nodeset 200 name "branch_fault"
+nodeset 101 group branch_fault
+nodeset 101 name "branch_fault"
# ----------------------------------------------------------------------
# Create nodeset for +x face
@@ -38,7 +33,7 @@
# Create nodeset for -x face
# ----------------------------------------------------------------------
group "face_xneg" add node in surface 12
-group "face_xneg" add node in surface 14
+group "face_xneg" add node in surface 15
nodeset 12 group face_xneg
nodeset 12 name "face_xneg"
@@ -77,21 +72,12 @@
# ----------------------------------------------------------------------
# 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 < -10.001e+3
-#nodeset 20 group face_zpos_subset
-#nodeset 20 name "face zpos_subset"
+group "face_zpos_subset" add node in group face_zpos
+group "face_zpos_subset" remove node with x_coord > {offFaultMaxX}
+group "face_zpos_subset" remove node with x_coord < {offFaultMinX}
+group "face_zpos_subset" remove node with y_coord > {offFaultMaxY}
+group "face_zpos_subset" remove node with y_coord < {offFaultMinY}
+nodeset 20 group face_zpos_subset
+nodeset 20 name "face zpos_subset"
-# ----------------------------------------------------------------------
-# Create nodeset for off-fault stations
-# ----------------------------------------------------------------------
-group "off_fault" add node with x_min >= {offfaultMinX} and x_max <= {offfaultMaxX} and y_min >= {offfaultMinY} and y_max <= {offfaultMaxY} and z_coord = {offfaultCoordZ}
-nodeset 25 group off_fault
-nodeset 25 name "off_fault"
-
-
-
-
+# End of file
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/gradient.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/gradient.jou 2011-02-15 00:41:26 UTC (rev 17864)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/gradient.jou 2011-02-15 00:59:48 UTC (rev 17865)
@@ -1,7 +1,11 @@
-#{bias_factor=1.2}
+# ----------------------------------------------------------------------
+# Set vertex spacing with increasing spacing away from fault
+# ----------------------------------------------------------------------
+#{bias_factor=1.01}
-volume all size {dx}
-#surface fault_surface_mainA fault_surface_mainB fault_surface_branch size {dx}
+curve all scheme default
+surface all sizing function none
+surface fault_surface_mainA fault_surface_mainB fault_surface_branch size {dx}
# Set bias on curves extending from faults
curve 70 scheme bias fine size {dx} factor {bias_factor} start vertex 46
@@ -18,12 +22,92 @@
curve 75 scheme bias fine size {dx} factor {bias_factor} start vertex 54
-# Extend to surfaces
-surface 31 sizing function type bias start curve 59 60 factor {bias_factor}
-surface 9 sizing function type bias start curve 59 70 factor {bias_factor}
+# ----------------------------------------------------------------------
+# Extend bias based on 1st set of bias values
-surface 10 sizing function type bias start curve 28 32 factor {bias_factor}
-surface 11 sizing function type bias start curve 29 32 factor {bias_factor}
+# Must be set manually.
+#{dxA=452.58} # end size of curve 59
+#{dxB=440.784} # end size of curve 63
+#{dxC=511.68} # end size of curve 83
+# A to A
+curve 56 57 size {dxA}
+curve 85 92 93 size {dxB}
+# C to C
+curve 77 size {dxC}
+# bias starting at A
+curve 21 scheme bias fine size {dxA} factor {bias_factor} start vertex 16
+curve 22 scheme bias fine size {dxA} factor {bias_factor} start vertex 14
+curve 29 scheme bias fine size {dxA} factor {bias_factor} start vertex 14
+
+# bias starting at B
+curve 64 scheme bias fine size {dxB} factor {bias_factor} start vertex 45
+curve 73 scheme bias fine size {dxB} factor {bias_factor} start vertex 48
+curve 86 scheme bias fine size {dxB} factor {bias_factor} start vertex 58
+
+# bias starting at C
+curve 49 scheme bias fine size {dxC} factor {bias_factor} start vertex 25
+
+# bias A to C
+curve 42 scheme bias fine size {dxA} coarse size {dxC} start vertex 16
+curve 42 scheme bias fine size {dxA} coarse size {dxC} start vertex 16
+
+
+# ----------------------------------------------------------------------
+# Extend bias based on 2nd set of bias values
+
+# Must be set manually.
+#{dxD=682.805} # end size of curve 64
+#{dxE=744.36} # end size of curve 86
+#{dxF=608.892} # end size of curve 49
+#{dxG=758.417} # end size of curve 22
+
+# bias starting at D
+curve 24 scheme bias fine size {dxD} factor {bias_factor} start vertex 13
+curve 31 scheme bias fine size {dxD} factor {bias_factor} start vertex 13
+curve 23 scheme bias fine size {dxD} factor {bias_factor} start vertex 15
+
+# bias starting at E
+curve 50 scheme bias fine size {dxE} factor {bias_factor} start vertex 26
+
+# bias A to D
+curve 54 scheme bias fine size {dxA} coarse size {dxD} start vertex 41
+curve 55 scheme bias fine size {dxA} coarse size {dxD} start vertex 42
+
+# bias D to E
+curve 41 scheme bias fine size {dxD} coarse size {dxE} start vertex 15
+
+# bias F to G
+curve 4 scheme bias fine size {dxF} coarse size {dxG} start vertex 1
+
+# bias C to E
+curve 76 scheme bias fine size {dxC} coarse size {dxE} start vertex 55
+
+# G to G
+curve 2 size {dxG}
+
+# ----------------------------------------------------------------------
+# Extend bias based on 3rd set of bias values
+
+# Must be set manually.
+#{dxH=804.251} # end size of curve 50
+#{dxI=975.139} # end size of curve 24
+
+# bias G to I
+curve 9 scheme bias fine size {dxG} coarse size {dxI} start vertex 4
+curve 11 scheme bias fine size {dxG} coarse size {dxI} start vertex 3
+curve 12 scheme bias fine size {dxG} coarse size {dxI} start vertex 2
+
+# bias F to H
+curve 10 scheme bias fine size {dxF} coarse size {dxH} start vertex 1
+
+# bias H to I
+curve 6 scheme bias fine size {dxH} coarse size {dxI} start vertex 6
+
+# I to I
+curve 8 size {dxI}
+
+
+# End of file
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/params.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/params.jou 2011-02-15 00:41:26 UTC (rev 17864)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/params.jou 2011-02-15 00:59:48 UTC (rev 17865)
@@ -34,12 +34,9 @@
#{dx=200*m}
#
# ----------------------------------------------------------------------
-# Parameters for off-fault stations.
+# Parameters for off-fault stations on ground surface
# ----------------------------------------------------------------------
-#
-# All off-fault stations are at Earth's surface so Z-coord=0.0
-#{offfaultMinX=-1.01*km}
-#{offfaultMaxX=+15.01*km}
-#{offfaultMinY=-10.01*km}
-#{offfaultMaxY=+10.01*km}
-#{offfaultMinZ=-10.0*m}
+#{offFaultMinX=-25.01*km}
+#{offFaultMaxX=+25.01*km}
+#{offFaultMinY=-20.01*km}
+#{offFaultMaxY=+20.01*km}
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tet4_200m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tet4_200m.jou 2011-02-15 00:41:26 UTC (rev 17864)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tet4_200m.jou 2011-02-15 00:59:48 UTC (rev 17865)
@@ -6,7 +6,7 @@
# Set parameters.
# ----------------------------------------------------------------------
playback 'params.jou'
-#{dx=800*m}
+#{dx=200*m}
# ----------------------------------------------------------------------
# Generate geometry
@@ -16,22 +16,30 @@
# ----------------------------------------------------------------------
# Set discretization size and scheme
# ----------------------------------------------------------------------
+surface all scheme trimesh
volume all scheme tetmesh
playback 'gradient.jou'
# ----------------------------------------------------------------------
# Generate the mesh
# ----------------------------------------------------------------------
+mesh surface all
+surface all smooth scheme laplacian
+smooth surface all
+surface all smooth scheme condition number beta 1.3 cpu 10
+smooth surface all
+
mesh volume all
# ----------------------------------------------------------------------
-# Refine cells near rupture
+# Smooth mesh to improve quality
# ----------------------------------------------------------------------
-playback 'refinecells.jou'
+cleanup volume all
+volume all smooth scheme condition number beta 2.0 cpu 10
+smooth volume all
+volume all smooth scheme condition number beta 1.9 cpu 10
+smooth volume all
-# Refine to get final resolution
-refine volume all numsplit 1
-
# ----------------------------------------------------------------------
# Setup boundary conditions.
# ----------------------------------------------------------------------
@@ -41,12 +49,3 @@
# Export exodus file
# ----------------------------------------------------------------------
export mesh "tet4_200m.exo" dimension 3 overwrite
-
-
-
-
-
-
-
-
-
More information about the CIG-COMMITS
mailing list