[cig-commits] r20693 - short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip
brad at geodynamics.org
brad at geodynamics.org
Thu Sep 6 08:09:36 PDT 2012
Author: brad
Date: 2012-09-06 08:09:35 -0700 (Thu, 06 Sep 2012)
New Revision: 20693
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/createbc.jou
short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/geometry.jou
short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/hex8_cubit.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/mat_elastic.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/mesh_hex8.jou
short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/risesliptime.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/slip.spatialdb
Removed:
short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/meshes/
short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/parameters/
short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/results/
Log:
Reorganized directory.
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/createbc.jou (from rev 17155, short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/meshes/createbc.jou)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/createbc.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/createbc.jou 2012-09-06 15:09:35 UTC (rev 20693)
@@ -0,0 +1,51 @@
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 volume elastic_xpos elastic_xneg
+block 1 name "elastic"
+
+# ----------------------------------------------------------------------
+# Create nodeset for fault
+# ----------------------------------------------------------------------
+group "fault" add node in fault_surface
+group "fault" add node in fault_surface at A
+group "fault" remove node with z_coord < -15.001e+3
+group "fault" remove node with y_coord < -15.001e+3
+group "fault" remove node with y_coord > +15.001e+3
+nodeset 10 group fault
+nodeset 10 name "fault"
+
+# ----------------------------------------------------------------------
+# Create nodesets for absorbing BC
+# ----------------------------------------------------------------------
+group "surf_xpos" add node in surface 6
+nodeset 11 group surf_xpos
+nodeset 11 name "surf_xpos"
+
+group "surf_xneg" add node in surface 4
+nodeset 12 group surf_xneg
+nodeset 12 name "surf_xneg"
+
+group "surf_ypos" add node in surface 12
+group "surf_ypos" add node in surface 14
+nodeset 13 group surf_ypos
+nodeset 13 name "surf_ypos"
+
+group "surf_yneg" add node in surface 10
+group "surf_yneg" add node in surface 17
+nodeset 14 group surf_yneg
+nodeset 14 name "surf_yneg"
+
+group "surf_zneg" add node in surface 11
+group "surf_zneg" add node in surface 15
+nodeset 15 group surf_zneg
+nodeset 15 name "surf_zneg"
+
+
+# ----------------------------------------------------------------------
+# Create nodeset for ground surface
+# ----------------------------------------------------------------------
+group "ground_surf" add node in surface 9
+group "ground_surf" add node in surface 16
+nodeset 16 group ground_surf
+nodeset 16 name "ground_surf"
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/geometry.jou (from rev 17155, short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/meshes/geometry.jou)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/geometry.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/geometry.jou 2012-09-06 15:09:35 UTC (rev 20693)
@@ -0,0 +1,41 @@
+# Geometry for simple dynamic benchmark with a kinematic fault rupture.
+#
+# Fault is 30km long and 15km wide
+# Fault runs north-south (y direction)
+#
+# Domain is 72.0 km by 36.0 km by 36.0 km
+# x: -18e+3 to +18e+3
+# y: -36e+3 to +36e+3
+# z: -36e+3 to 0
+
+# ----------------------------------------------------------------------
+# Create domain
+# ----------------------------------------------------------------------
+
+reset
+brick x 36.0e+3 y 72.0e+3 z 36.0e+3
+
+# Translate block so the top is at z=0
+volume 1 move x 0 y 0 z -18.0e+3
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+
+# Create fault (yz plane) at x = 0.0
+create planar surface with plane xplane offset 0
+surface 7 name "fault_surface"
+
+# ----------------------------------------------------------------------
+# Divide volumes using interface surfaces
+# ----------------------------------------------------------------------
+webcut volume 1 with plane surface fault_surface
+volume 1 name "elastic_xpos"
+volume 3 name "elastic_xneg"
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces
+# ----------------------------------------------------------------------
+imprint all with volume all
+merge all
+delete body 2
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/hex8_cubit.cfg (from rev 20123, short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/parameters/hex8_cubit.cfg)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/hex8_cubit.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/hex8_cubit.cfg 2012-09-06 15:09:35 UTC (rev 20693)
@@ -0,0 +1,83 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+meshiocubit = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1
+reader = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.reader]
+filename = ../meshes/halfspace_hex8_cubit.exo
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials.elastic]
+
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+label = surf_xpos
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.x_neg]
+label = surf_xneg
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_pos]
+label = surf_ypos
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_neg]
+label = surf_yneg
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.z_neg]
+label = surf_zneg
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+label = fault
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = ../results/pylith_hex8_cubit-volume.vtk
+
+[pylithapp.problem.formulation.output.groundsurface]
+writer.filename = ../results/pylith_hex8_cubit-ground.vtk
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = ../results/pylith_hex8_cubit-fault.vtk
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = ../results/pylith_hex8_cubit-statevars.vtk
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/mat_elastic.spatialdb (from rev 17906, short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/parameters/mat_elastic.spatialdb)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/mat_elastic.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/mat_elastic.spatialdb 2012-09-06 15:09:35 UTC (rev 20693)
@@ -0,0 +1,12 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = Density Vp Vs
+ value-units = kg/m**3 m/s m/s
+ num-locs = 1
+ data-dim = 0
+ cs-data = cartesian {
+ to-meters = 1.0
+ }
+}
+0.0 0.0 0.0 3050.0 7800.0 4500.0
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/mesh_hex8.jou (from rev 17155, short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/meshes/mesh_hex8.jou)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/mesh_hex8.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/mesh_hex8.jou 2012-09-06 15:09:35 UTC (rev 20693)
@@ -0,0 +1,29 @@
+# Mesh for simple dynamic benchmark with a kinematic fault rupture.
+#
+# Cell type: hex8
+# Cell size: 750m.
+
+# ----------------------------------------------------------------------
+# Generage geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+volume all size 750.0
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Generate blocks and nodesets
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "halfspace_hex8_cubit.exo" dimension 3 overwrite
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/pylithapp.cfg (from rev 20123, short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/parameters/pylithapp.cfg)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/pylithapp.cfg 2012-09-06 15:09:35 UTC (rev 20693)
@@ -0,0 +1,120 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+explicit = 1
+petsc = 1
+solverlinear = 1
+meshiolagrit = 1
+explicitelasticity = 1
+faultcohesivekin = 1
+fiatlagrange = 1
+meshimporter = 1
+
+[pylithapp.journal.debug]
+pylithapp = 1
+problem = 1
+explicit = 1
+
+# mesh generator
+[pylithapp.mesh_generator]
+#reorder_mesh = True
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 3
+
+formulation = pylith.problems.ExplicitLumped
+
+bc = [x_pos,x_neg,y_pos,y_neg,z_neg]
+bc.x_pos = pylith.bc.AbsorbingDampers
+bc.x_neg = pylith.bc.AbsorbingDampers
+bc.y_pos = pylith.bc.AbsorbingDampers
+bc.y_neg = pylith.bc.AbsorbingDampers
+bc.z_neg = pylith.bc.AbsorbingDampers
+bc.z_neg.up_dir = [1,0,0]
+
+interfaces = [fault]
+materials = [elastic]
+formulation.output = [domain,groundsurface]
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.25*s
+dt = 0.05*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials.elastic]
+
+label = Elastic material
+id = 1
+db_properties.iohandler.filename = mat_elastic.spatialdb
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+db.label = Absorbing BC +x
+db.iohandler.filename = mat_elastic.spatialdb
+
+[pylithapp.timedependent.bc.x_neg]
+db.label = Absorbing BC -x
+db.iohandler.filename = mat_elastic.spatialdb
+
+[pylithapp.timedependent.bc.y_pos]
+db.label = Absorbing BC +y
+db.iohandler.filename = mat_elastic.spatialdb
+
+[pylithapp.timedependent.bc.y_neg]
+db.label = Absorbing BC -y
+db.iohandler.filename = mat_elastic.spatialdb
+
+[pylithapp.timedependent.bc.z_neg]
+db.label = Absorbing BC -z
+db.iohandler.filename = mat_elastic.spatialdb
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+id = 100
+
+eq_srcs.rupture.slip_function = pylith.faults.LiuCosSlipFn
+
+[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
+slip.iohandler.filename = slip.spatialdb
+rise_time.iohandler.filename = risesliptime.spatialdb
+slip_time.iohandler.filename = risesliptime.spatialdb
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+log_summary = true
+
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.time_format = %05.2f
+vertex_data_fields = []
+
+[pylithapp.problem.formulation.output.groundsurface]
+skip = 1
+writer.time_format = %05.2f
+
+[pylithapp.timedependent.interfaces.fault.output]
+skip = 1
+writer.time_format = %05.2f
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.time_format = %05.2f
+cell_filter = pylith.meshio.CellFilterAvgMesh
+cell_data_fields = []
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/risesliptime.spatialdb (from rev 17906, short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/parameters/risesliptime.spatialdb)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/risesliptime.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/risesliptime.spatialdb 2012-09-06 15:09:35 UTC (rev 20693)
@@ -0,0 +1,13 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 2
+ value-names = rise-time slip-time
+ value-units = s s
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1
+ }
+}
+0.0 0.0 0.0 2.0 0.0
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/slip.spatialdb (from rev 17906, short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/parameters/slip.spatialdb)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/slip.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/performance/halfspace_strikeslip/slip.spatialdb 2012-09-06 15:09:35 UTC (rev 20693)
@@ -0,0 +1,13 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = left-lateral-slip reverse-slip fault-opening
+ value-units = m m m
+ num-locs = 1
+ data-dim = 0
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1
+ }
+}
+0.0 0.0 0.0 -1.0 0.0 0.0
More information about the CIG-COMMITS
mailing list