[cig-commits] r13315 - in short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip: meshes parameters
brad at geodynamics.org
brad at geodynamics.org
Fri Nov 14 17:21:45 PST 2008
Author: brad
Date: 2008-11-14 17:21:45 -0800 (Fri, 14 Nov 2008)
New Revision: 13315
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/meshes/createbc.jou
short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/meshes/geometry.jou
short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/meshes/halfspace_hex8_cubit.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/meshes/mesh_hex8.jou
short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/hex8_cubit.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/tet4_lagrit.cfg
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/pylithapp.cfg
Log:
Added hex8 case.
Added: short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/meshes/createbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/meshes/createbc.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/meshes/createbc.jou 2008-11-15 01:21:45 UTC (rev 13315)
@@ -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"
Added: short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/meshes/geometry.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/meshes/geometry.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/meshes/geometry.jou 2008-11-15 01:21:45 UTC (rev 13315)
@@ -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
Added: short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/meshes/halfspace_hex8_cubit.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/meshes/halfspace_hex8_cubit.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/meshes/mesh_hex8.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/meshes/mesh_hex8.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/meshes/mesh_hex8.jou 2008-11-15 01:21:45 UTC (rev 13315)
@@ -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
Added: short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/hex8_cubit.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/hex8_cubit.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/hex8_cubit.cfg 2008-11-15 01:21:45 UTC (rev 13315)
@@ -0,0 +1,84 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+meshiocubit = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1
+importer = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.importer]
+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 = 11
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.x_neg]
+label = 12
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_pos]
+label = 13
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_neg]
+label = 14
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.z_neg]
+label = 15
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+label = 10
+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
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/pylithapp.cfg 2008-11-14 23:30:06 UTC (rev 13314)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/pylithapp.cfg 2008-11-15 01:21:45 UTC (rev 13315)
@@ -13,17 +13,11 @@
explicitelasticity = 1
faultcohesivekin = 1
-# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator]
-#debug = 1
-importer = pylith.meshio.MeshIOLagrit
+[pylithapp.journal.debug]
+pylithapp = 1
+problem = 1
+explicit = 1
-[pylithapp.mesh_generator.importer]
-filename_gmv = ../meshes/halfspace_tet4_lagrit.gmv
-filename_pset = ../meshes/halfspace_tet4_lagrit.pset
-
# ----------------------------------------------------------------------
# problem
# ----------------------------------------------------------------------
@@ -43,6 +37,7 @@
interfaces = [fault]
materials = [elastic]
+formulation.output = [domain,groundsurface]
[pylithapp.timedependent.formulation.time_step]
total_time = 15.0*s
@@ -57,68 +52,41 @@
id = 1
db.iohandler.filename = mat_elastic.spatialdb
quadrature = pylith.feassemble.quadrature.Quadrature3D
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = tetrahedron
-quadrature.cell.quad_order = 2
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
[pylithapp.timedependent.bc.x_pos]
-label = elemfaces_xmax
db.label = Absorbing BC +x
db.iohandler.filename = mat_elastic.spatialdb
quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 2
[pylithapp.timedependent.bc.x_neg]
-label = elemfaces_xmin
db.label = Absorbing BC -x
db.iohandler.filename = mat_elastic.spatialdb
quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 2
[pylithapp.timedependent.bc.y_pos]
-label = elemfaces_ymax
db.label = Absorbing BC +y
db.iohandler.filename = mat_elastic.spatialdb
quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 2
[pylithapp.timedependent.bc.y_neg]
-label = elemfaces_ymin
db.label = Absorbing BC -y
db.iohandler.filename = mat_elastic.spatialdb
quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 2
[pylithapp.timedependent.bc.z_neg]
-label = elemfaces_zmin
db.label = Absorbing BC -z
db.iohandler.filename = mat_elastic.spatialdb
quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 2
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
[pylithapp.timedependent.interfaces.fault]
id = 100
-label = elemfaces_fault
quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 2
mat_db.iohandler.filename = mat_elastic.spatialdb
eq_srcs.rupture.slip_function = pylith.faults.LiuCosSlipFn
@@ -147,18 +115,19 @@
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
-[pylithapp.problem.formulation.output.output]
+[pylithapp.problem.formulation.output.domain]
+skip = 99999
+writer.time_format = %05.2f
+
+[pylithapp.problem.formulation.output.groundsurface]
skip = 1
-writer.filename = ../results/pylith_tet4_lagrit.vtk
writer.time_format = %05.2f
[pylithapp.timedependent.interfaces.fault.output]
skip = 1
-writer.filename = ../results/pylith_tet4_lagrit-fault.vtk
writer.time_format = %05.2f
[pylithapp.timedependent.materials.elastic.output]
-skip = 19
-writer.filename = ../results/pylith_tet4_lagrit-statevars.vtk
+skip = 999999
writer.time_format = %05.2f
cell_filter = pylith.meshio.CellFilterAvg
Added: short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/tet4_lagrit.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/tet4_lagrit.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/halfspace_strikeslip/parameters/tet4_lagrit.cfg 2008-11-15 01:21:45 UTC (rev 13315)
@@ -0,0 +1,87 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+meshiolagrit = 1
+gmvfile = 1
+psetfile = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1
+importer = pylith.meshio.MeshIOLagrit
+
+[pylithapp.mesh_generator.importer]
+filename_gmv = ../meshes/halfspace_tet4_lagrit.gmv
+filename_pset = ../meshes/halfspace_tet4_lagrit.pset
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials.elastic]
+
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = tetrahedron
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+label = elemfaces_xmax
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.x_neg]
+label = elemfaces_xmin
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_pos]
+label = elemfaces_ymax
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_neg]
+label = elemfaces_ymin
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.z_neg]
+label = elemfaces_zmin
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.fault]
+label = elemfaces_fault
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = ../results/pylith_tet4_lagrit-volume.vtk
+
+[pylithapp.problem.formulation.output.groundsurface]
+writer.filename = ../results/pylith_tet4_lagrit-ground.vtk
+
+[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = ../results/pylith_tet4_lagrit-fault.vtk
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = ../results/pylith_tet4_lagrit-statevars.vtk
More information about the CIG-COMMITS
mailing list