[cig-commits] r16748 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d
brad at geodynamics.org
brad at geodynamics.org
Wed May 19 09:53:17 PDT 2010
Author: brad
Date: 2010-05-19 09:53:17 -0700 (Wed, 19 May 2010)
New Revision: 16748
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m_refine.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m_refine.jou
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry_refine.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.jou
Log:
Improved refinement.
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry_refine.jou 2010-05-19 14:28:19 UTC (rev 16747)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry_refine.jou 2010-05-19 16:53:17 UTC (rev 16748)
@@ -1,13 +1,13 @@
# ----------------------------------------------------------------------
# Create block
# ----------------------------------------------------------------------
-# Block is 72.0 km x 144.0 km
-# -36.0 km <= x <= 36.0 km
-# -72.0 km <= y <= 72.0 km
+# Block is 128.0 km x 128.0 km
+# -64.0 km <= x <= 64.0 km
+# -64.0 km <= y <= 64.0 km
#
# We create the 2-D mesh by extracting the mid-surface from a 3-D mesh.
reset
-brick x 72000 y 144000 z 72000
+brick x 128000 y 128000 z 128000
create midsurface volume 1 surface 1 2
delete volume 1
@@ -16,7 +16,7 @@
# ----------------------------------------------------------------------
# Create fault (yz plane) at x = 0.0
-split surface 7 across location position 0 -72000 0 location position 0 72000 0
+split surface 7 across location position 0 -64000 0 location position 0 64000 0
curve 17 name "fault_surface"
# ----------------------------------------------------------------------
@@ -24,6 +24,3 @@
# ----------------------------------------------------------------------
imprint all
merge all
-
-
-
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.cfg 2010-05-19 16:53:17 UTC (rev 16748)
@@ -0,0 +1,39 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = quad4_200m_refine.exo
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+dt = 0.01666666666666666666666667*s
+
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.output]
+writer.filename = output/quad4_200m_refine.vtk
+
+# Give basename for VTK domain output of solution over Off-fault stations.
+#[pylithapp.problem.formulation.output.subdomain]
+#label = 15 ; nodeset for subdomain
+#vertex_data_fields=[displacement,velocity]
+#output_freq = time_step
+#time_step = 0.04999*s
+#writer.filename = output/quad4_200m_OffFault.vtk
+#writer.time_format = %05.2f
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/quad4_200m_refine-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/quad4_200m_refine-elastic.vtk
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/x-gzip
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.jou 2010-05-19 14:28:19 UTC (rev 16747)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m_refine.jou 2010-05-19 16:53:17 UTC (rev 16748)
@@ -20,9 +20,16 @@
# ----------------------------------------------------------------------
# Refine elements within 1/2 of rupture width
# ----------------------------------------------------------------------
-group "Refined" add node with x_min >= -15000 and x_max <= 15000 and y_min >= -30000 and y_max <= 30000
+group "refinecells1" add node with x_min >= -24.001e+3 and x_max <= 24.001e+3 and y_min >= -36.001e+3 and y_max <= 36.001e+3
+refine node in refinecells1 depth 1
+group "refinecells2" add node with x_min >= -12.001e+3 and x_max <= 12.001e+3 and y_min >= -24.001e+3 and y_max <= 24.001e+3
+refine node in refinecells2 depth 1
-refine node in Refined numsplit 2
+# ----------------------------------------------------------------------
+# Smooth mesh to improve quality
+# ----------------------------------------------------------------------
+surface all smooth scheme condition number beta 1.3 cpu 10
+smooth surface all
# ----------------------------------------------------------------------
# Setup boundary conditions.
@@ -33,6 +40,3 @@
# Export exodus file
# ----------------------------------------------------------------------
export mesh "quad4_200m_refine.exo" dimension 2 overwrite
-
-
-
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m_refine.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m_refine.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/x-gzip
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m_refine.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m_refine.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tri3_200m_refine.jou 2010-05-19 16:53:17 UTC (rev 16748)
@@ -0,0 +1,49 @@
+# ----------------------------------------------------------------------
+# Create tri3 mesh at 200m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry_refine.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size and scheme
+# ----------------------------------------------------------------------
+surface all size 800
+surface all scheme trimesh
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh surface all
+
+# ----------------------------------------------------------------------
+# Refine elements within 1/2 of rupture width
+# ----------------------------------------------------------------------
+group "refinecells1" add node with x_min >= -24.001e+3 and x_max <= 24.001e+3 and y_min >= -36.001e+3 and y_max <= 36.001e+3
+refine node in refinecells1 depth 1
+group "refinecells2" add node with x_min >= -12.001e+3 and x_max <= 12.001e+3 and y_min >= -24.001e+3 and y_max <= 24.001e+3
+refine node in refinecells2 depth 1
+
+# ----------------------------------------------------------------------
+# Smooth mesh to improve quality
+# ----------------------------------------------------------------------
+surface all smooth scheme laplacian
+smooth surface all
+surface all smooth scheme condition number beta 1.3 cpu 10
+smooth surface all
+
+# ----------------------------------------------------------------------
+# Setup boundary conditions.
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "tri3_200m_refine.exo" dimension 2 overwrite
+
+
+
+
More information about the CIG-COMMITS
mailing list