[cig-commits] r17860 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14
brad at geodynamics.org
brad at geodynamics.org
Mon Feb 14 11:48:56 PST 2011
Author: brad
Date: 2011-02-14 11:48:56 -0800 (Mon, 14 Feb 2011)
New Revision: 17860
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/gradient.jou
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/geometry.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/params.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tet4_200m.jou
Log:
A little more work on tpv14 tet4 mesh.
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/geometry.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/geometry.jou 2011-02-14 05:45:26 UTC (rev 17859)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/geometry.jou 2011-02-14 19:48:56 UTC (rev 17860)
@@ -13,21 +13,60 @@
# Create interface surfaces
# ----------------------------------------------------------------------
-# Create Main Fault
+# Create main fault
create planar surface with plane yplane offset 0
-surface 7 name "fault_surface_main"
-webcut volume 1 with plane surface fault_surface_main
+surface 7 name "temp_main"
+webcut volume 1 with sheet surface temp_main
+delete body 2
-# Create Branch Fault
+# Create branch fault
create planar surface with plane xplane offset {blockLength/2.0}
-Rotate body 4 about vertex 2 5 angle {branchAngle-90}
-move body 4 x {branchOffsetY-blockWidth/2.0} y {-blockWidth/2.0}
-webcut volume 3 with sheet body 4
-surface 19 name "fault_surface_branch"
-
+surface 18 name "temp_branch"
+rotate surface temp_branch about vertex 2 5 angle {branchAngle-90}
+move surface temp_branch x {branchOffsetY-blockLength/2.0} y {-blockWidth/2.0}
+webcut volume 3 with sheet surface temp_branch
+delete body 4
+
# ----------------------------------------------------------------------
+# Split surfaces involving faults into pieces, isolating rupture regions
+# ----------------------------------------------------------------------
+
+# Split main fault surface using volume 1
+create vertex {-mainFaultLength/2} 0 0
+create vertex {+mainFaultLength/2} 0 0
+create vertex {-mainFaultLength/2} 0 {-faultWidth}
+create vertex {+mainFaultLength/2} 0 {-faultWidth}
+create vertex {-mainFaultLength/2} 0 {-blockHeight}
+create vertex {+mainFaultLength/2} 0 {-blockHeight}
+create vertex {-blockLength/2} 0 {-faultWidth}
+create vertex {+blockLength/2} 0 {-faultWidth}
+
+split surface 8 across location vertex 39 location vertex 40
+split surface 30 across location vertex 33 location vertex 35
+split surface 29 across location vertex 35 location vertex 37
+split surface 32 across location vertex 34 location vertex 36
+split surface 34 across location vertex 36 location vertex 38
+delete vertex all
+
+# Split branch fault surface using volume 3
+vertex 27 copy move z {-faultWidth}
+vertex 25 copy move z {-faultWidth}
+create vertex on curve 48 distance {12*km} from end
+vertex 51 copy move z {-faultWidth}
+vertex 51 copy move z {-blockHeight}
+
+split surface 19 across location vertex 49 location vertex 50
+split surface 40 across location vertex 51 location vertex 52
+split surface 39 across location vertex 52 location vertex 53
+delete vertex all
+
+# ----------------------------------------------------------------------
# Imprint all volumes, merging surfaces
# ----------------------------------------------------------------------
imprint all
merge all
-delete body 2 4
+
+surface 45 name "fault_surface_mainB"
+surface 48 name "fault_surface_mainA"
+surface 41 name "fault_surface_branch"
+
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/gradient.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/gradient.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/gradient.jou 2011-02-14 19:48:56 UTC (rev 17860)
@@ -0,0 +1,29 @@
+#{bias_factor=1.2}
+
+volume all size {dx}
+#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
+curve 69 scheme bias fine size {dx} factor {bias_factor} start vertex 47
+curve 71 scheme bias fine size {dx} factor {bias_factor} start vertex 47
+
+curve 59 scheme bias fine size {dx} factor {bias_factor} start vertex 43
+curve 60 scheme bias fine size {dx} factor {bias_factor} start vertex 44
+curve 63 scheme bias fine size {dx} factor {bias_factor} start vertex 44
+
+curve 83 scheme bias fine size {dx} factor {bias_factor} start vertex 56
+curve 82 scheme bias fine size {dx} factor {bias_factor} start vertex 57
+curve 84 scheme bias fine size {dx} factor {bias_factor} start vertex 57
+
+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}
+
+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}
+
+
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/params.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/params.jou 2011-02-14 05:45:26 UTC (rev 17859)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/params.jou 2011-02-14 19:48:56 UTC (rev 17860)
@@ -6,68 +6,40 @@
# Parameters for geometry.
# ----------------------------------------------------------------------
#
-# Block is 64.0 km x 64.0 km x 32.0 km
+# Block is 64.0 km x 80.0 km x 40.0 km
# -32.0 km <= x <= 32.0 km
-# -32.0 km <= y <= 32.0 km
-# -32.0 km <= z <= 0 km
+# -40.0 km <= y <= 40.0 km
+# -40.0 km <= z <= 0 km
#
#{blockWidth=64*km}
-#{blockLength=64*km}
-#{blockHeight=32*km}
+#{blockLength=80*km}
+#{blockHeight=40*km}
#
-# Main fault is 28.0 km long
+# Faults have a common width of 15 km
+#{faultWidth=15.0*km}
#
-#{mainFaultLength=28.1*km}
-#{mainFaultDepth=15.1*km}
+# Main fault is 28.0 km long.
+#{mainFaultLength=28.0*km}
#
# Branch fault is 12.0 km long and intersects main fault at y=+2.0 km
# at angle of 30.0 degrees.
#
#{branchAngle=30.0}
-#{branchFaultLength=12.1*km}
-#{branchFaultDepth=15.1*km}
+#{branchFaultLength=12.0*km}
#{branchOffsetY=2*km}
#
-# Note: Fault lengths are slightly longer than necessary to insure that
-# we capture the nodes that span the rupture.
-#
# ----------------------------------------------------------------------
# Parameters for mesh.
# ----------------------------------------------------------------------
-#
#{dx=200*m}
#
# ----------------------------------------------------------------------
-# Parameters for refinement.
-# ----------------------------------------------------------------------
-#
-# Level 1
-#
-# Within 24 km of deeper than fault depth
-# Within 24 km of rupture perpendicular to strike
-# Within 24 km of rupture along strike
-#{refine1Depth=40.01*km}
-#{refine1Width=24.01*km}
-#{refine1Length=24.01*km}
-#
-# Level 2
-#
-# Within 12 km of deeper than fault depth
-# Within 12 km of rupture perpendicular to strike
-# Within 12 km of rupture along strike
-#{refine2Depth=28.81*km}
-#{refine2Width=12.81*km}
-#{refine2Length=12.81*km}
-#
-# ----------------------------------------------------------------------
# Parameters for off-fault stations.
# ----------------------------------------------------------------------
#
-# All off-fault stations are at earth's surface so Z-coord=0.0
+# All off-fault stations are at Earth's surface so Z-coord=0.0
#{offfaultMinX=-1.01*km}
#{offfaultMaxX=+15.01*km}
-#{offfaultMinY=-8.51*km}
-#{offfaultMaxY=+4.01*km}
-#{offfaultCoordZ=0.0*km}
-
-
+#{offfaultMinY=-10.01*km}
+#{offfaultMaxY=+10.01*km}
+#{offfaultMinZ=-10.0*m}
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-14 05:45:26 UTC (rev 17859)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv14/tet4_200m.jou 2011-02-14 19:48:56 UTC (rev 17860)
@@ -6,6 +6,7 @@
# Set parameters.
# ----------------------------------------------------------------------
playback 'params.jou'
+#{dx=800*m}
# ----------------------------------------------------------------------
# Generate geometry
@@ -15,8 +16,8 @@
# ----------------------------------------------------------------------
# Set discretization size and scheme
# ----------------------------------------------------------------------
-volume all size {16*dx}
-volume all scheme tetmesh
+volume all scheme tetmesh
+playback 'gradient.jou'
# ----------------------------------------------------------------------
# Generate the mesh
@@ -48,3 +49,4 @@
+
More information about the CIG-COMMITS
mailing list