[cig-commits] r18991 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205
brad at geodynamics.org
brad at geodynamics.org
Thu Sep 29 09:21:52 PDT 2011
Author: brad
Date: 2011-09-29 09:21:52 -0700 (Thu, 29 Sep 2011)
New Revision: 18991
Removed:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/geometry_refine.jou
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/createbc.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/geometry.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m.jou
Log:
Updated tet4 meshing procedure.
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/createbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/createbc.jou 2011-09-29 10:21:47 UTC (rev 18990)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/createbc.jou 2011-09-29 16:21:52 UTC (rev 18991)
@@ -8,41 +8,38 @@
# Create nodeset for fault
# ----------------------------------------------------------------------
group "fault" add node in fault_surface
-group "fault" remove node with x_coord > +15.0001e+3
-group "fault" remove node with x_coord < -15.0001e+3
-group "fault" remove node with z_coord < -15.0001e+3
-nodeset 10 group fault
-nodeset 10 name "fault"
+nodeset 20 group fault
+nodeset 20 name "fault"
# ----------------------------------------------------------------------
# Create nodeset for +x face
# ----------------------------------------------------------------------
-group "face_xpos" add node in surface 10
-group "face_xpos" add node in surface 17
-nodeset 11 group face_xpos
-nodeset 11 name "face_xpos"
+group "face_xpos" add node in surface 6
+nodeset 10 group face_xpos
+nodeset 10 name "face_xpos"
# ----------------------------------------------------------------------
# Create nodeset for -x face
# ----------------------------------------------------------------------
-group "face_xneg" add node in surface 12
-group "face_xneg" add node in surface 14
-nodeset 12 group face_xneg
-nodeset 12 name "face_xneg"
+group "face_xneg" add node in surface 4
+nodeset 11 group face_xneg
+nodeset 11 name "face_xneg"
# ----------------------------------------------------------------------
# Create nodeset for +y face
# ----------------------------------------------------------------------
-group "face_ypos" add node in surface 5
-nodeset 13 group face_ypos
-nodeset 13 name "face_ypos"
+group "face_ypos" add node in surface 12
+group "face_ypos" add node in surface 14
+nodeset 12 group face_ypos
+nodeset 12 name "face_ypos"
# ----------------------------------------------------------------------
# Create nodeset for -y face
# ----------------------------------------------------------------------
-group "face_yneg" add node in surface 3
-nodeset 14 group face_yneg
-nodeset 14 name "face_yneg"
+group "face_yneg" add node in surface 10
+group "face_yneg" add node in surface 17
+nodeset 13 group face_yneg
+nodeset 13 name "face_yneg"
# ----------------------------------------------------------------------
# Create nodeset for -z face
@@ -68,13 +65,5 @@
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 < -5.001e+3
-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 >= -12.001e+3 and x_max <= 12.001e+3 and z_min >= -7.501e+3 and y_coord = 3.000e+3
-nodeset 25 group off_fault
-nodeset 25 name "off_fault"
-
+nodeset 17 group face_zpos_subset
+nodeset 17 name "face zpos_subset"
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/geometry.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/geometry.jou 2011-09-29 10:21:47 UTC (rev 18990)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/geometry.jou 2011-09-29 16:21:52 UTC (rev 18991)
@@ -1,29 +1,58 @@
# ----------------------------------------------------------------------
# Create block
# ----------------------------------------------------------------------
-# Block is 40.8 km x 20.4 km x 20.4 km
-# -20.4 km <= x <= 20.4 km
-# -10.2 km <= y <= 10.2 km
-# -20.4 km <= z <= 0.0 km
+# Block is 48.0 km x 64.0 km x 32.0 km
+# -24.0 km <= x <= 24.0 km
+# -32.0 km <= y <= 32.0 km
+# -32.0 km <= z <= 0.0 km
#
+#{Units('si')}
+#
+#{blockWidth=48.0*km}
+#{blockLength=64.0*km}
+#{blockHeight=32.0*km}
+#
+#{faultLength=30.0*km}
+#{faultWidth=15.0*km}
+
reset
-brick x 40800 y 20400 z 20400
+brick x {blockWidth} y {blockLength} z {blockHeight}
# Translate block so the top is at z=0
-volume 1 move x 0 y 0 z -10200
+volume 1 move x 0 y 0 z {-0.5*blockHeight}
# ----------------------------------------------------------------------
# Create interface surfaces
# ----------------------------------------------------------------------
-create planar surface with plane yplane offset 0
-surface 7 name "fault_surface"
+create planar surface with plane xplane offset 0
+#{sI=Id("surface")}
# ----------------------------------------------------------------------
# Divide volumes using interface surfaces
# ----------------------------------------------------------------------
-webcut volume 1 with plane surface fault_surface
+webcut volume 1 with plane surface {sI}
# ----------------------------------------------------------------------
+# Split fault surface
+# ----------------------------------------------------------------------
+create vertex 0 {-0.5*faultLength} 0
+#{vAn=Id("vertex")}
+create vertex 0 {+0.5*faultLength} 0
+#{vAp=Id("vertex")}
+create vertex 0 {-0.5*faultLength} {-blockHeight}
+#{vBn=Id("vertex")}
+create vertex 0 {+0.5*faultLength} {-blockHeight}
+#{vBp=Id("vertex")}
+create vertex 0 {-0.5*faultLength} {-faultWidth}
+#{vCn=Id("vertex")}
+create vertex 0 {+0.5*faultLength} {-faultWidth}
+#{vCp=Id("vertex")}
+split surface 7 across location vertex {vAn} location vertex {vBn}
+split surface 18 across location vertex {vAp} location vertex {vBp}
+split surface 21 across location vertex {vCn} location vertex {vCp}
+surface 22 name "fault_surface"
+
+# ----------------------------------------------------------------------
# Imprint all volumes, merging surfaces
# ----------------------------------------------------------------------
imprint all
@@ -31,3 +60,5 @@
delete body 2
+
+
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/geometry_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/geometry_refine.jou 2011-09-29 10:21:47 UTC (rev 18990)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/geometry_refine.jou 2011-09-29 16:21:52 UTC (rev 18991)
@@ -1,34 +0,0 @@
-# ----------------------------------------------------------------------
-# Create block
-# ----------------------------------------------------------------------
-# Block is 128.0 km x 128.0 km
-# -64.0 km <= x <= 64.0 km
-# -64.0 km <= y <= 64.0 km
-# -64.0 km <= z <= 0.0 km
-#
-# We create the 2-D mesh by extracting the mid-surface from a 3-D mesh.
-reset
-brick x 128000 y 128000 z 64000
-
-# Translate block so the top is at z=0
-volume 1 move x 0 y 0 z -32000
-
-# ----------------------------------------------------------------------
-# Create interface surfaces
-# ----------------------------------------------------------------------
-create planar surface with plane yplane offset 0
-surface 7 name "fault_surface"
-
-# ----------------------------------------------------------------------
-# Divide volumes using interface surfaces
-# ----------------------------------------------------------------------
-webcut volume 1 with plane surface fault_surface
-
-# ----------------------------------------------------------------------
-# Imprint all volumes, merging surfaces
-# ----------------------------------------------------------------------
-imprint all
-merge all
-delete body 2
-
-
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m.jou 2011-09-29 10:21:47 UTC (rev 18990)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205/tet4_200m.jou 2011-09-29 16:21:52 UTC (rev 18991)
@@ -1,5 +1,5 @@
# ----------------------------------------------------------------------
-# Create tet4 mesh at 200m resolution.
+# Create tet4 mesh at 100m resolution.
# ----------------------------------------------------------------------
# ----------------------------------------------------------------------
@@ -10,18 +10,32 @@
# ----------------------------------------------------------------------
# Set discretization size and scheme
# ----------------------------------------------------------------------
-volume all size 200
+#{dx=200*m}
+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
# ----------------------------------------------------------------------
-# Improve element quality.
+# Smooth mesh to improve quality
# ----------------------------------------------------------------------
-playback 'tet4_smooth.jou'
+#{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.
@@ -33,7 +47,3 @@
# ----------------------------------------------------------------------
export mesh "tet4_200m.exo" dimension 3 overwrite
-
-
-
-
More information about the CIG-COMMITS
mailing list