[cig-commits] r16579 - in short/3D/PyLith/trunk/playpen/friction: . SCEC_benchmarks
surendra at geodynamics.org
surendra at geodynamics.org
Sun Apr 25 20:41:45 PDT 2010
Author: surendra
Date: 2010-04-25 20:41:45 -0700 (Sun, 25 Apr 2010)
New Revision: 16579
Added:
short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/
short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004.cfg
short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004.exo
short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004.jou
short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004initTract.spatialdb
short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004matprops.spatialdb
Log:
Added a smaller version of TPV5 SCEC benchmark
Added: short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004.cfg (rev 0)
+++ short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004.cfg 2010-04-26 03:41:45 UTC (rev 16579)
@@ -0,0 +1,206 @@
+# -*- Python -*-
+
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+#absorbingdampers = 1
+#explicit = 1
+#petsc = 1
+#meshiocubit = 1
+#explicitelasticity = 1
+#quadrature2d = 1
+#fiatlagrange = 1
+#faultcohesivedyn = 1
+
+[pylithapp.journal.debug]
+#pylithapp = 1
+#problem = 1
+#explicit = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+reader = pylith.meshio.MeshIOCubit
+
+[pylithapp.mesh_generator.reader]
+filename = SCEC2004.exo
+coordsys.space_dim = 3
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+
+dimension = 3
+
+# Set interfaces to an array with 1 faults: 'fault'.
+interfaces = [fault]
+
+# Set materials to an array with 1 material 'elastic'.
+materials = [elastic]
+
+formulation = pylith.problems.Explicit
+formulation.solver = pylith.problems.SolverNonlinear
+#formulation = pylith.problems.ExplicitLumped
+
+# Nondimensionalize problem using wave propagation parameters.
+normalizer = spatialdata.units.NondimElasticDynamic
+#normalizer.mass_density = 1000.0*kg/m**3
+#normalizer.shear_wave_speed = 1.0*m/s
+#normalizer.wave_period = 1.0*s
+
+bc = [x_neg,x_pos,z_neg,y_neg,y_pos]
+
+# Change the DirichletPoints BC to the AbsorbingDampers BC
+bc.x_neg = pylith.bc.AbsorbingDampers
+bc.x_pos = pylith.bc.AbsorbingDampers
+bc.z_neg = pylith.bc.AbsorbingDampers
+bc.y_neg = pylith.bc.AbsorbingDampers
+bc.y_pos = pylith.bc.AbsorbingDampers
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.1*s
+dt = 0.01*s
+
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+
+materials = [elastic]
+materials.elastic = pylith.materials.ElasticIsotropic3D
+#materials.elastic = pylith.materials.Homogeneous
+
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
+id = 1
+db_properties.iohandler.filename = SCEC2004matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 3
+
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+label = 26
+db.label = Absorbing BC +x
+db.iohandler.filename = SCEC2004matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+
+[pylithapp.timedependent.bc.x_neg]
+label = 27
+db.label = Absorbing BC -x
+db.iohandler.filename = SCEC2004matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+
+[pylithapp.timedependent.bc.z_neg]
+label = 28
+up_dir = [0, 1, 0]
+db.label = Absorbing BC -z
+db.iohandler.filename = SCEC2004matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+
+[pylithapp.timedependent.bc.y_pos]
+label = 22
+db.label = Absorbing BC +y
+db.iohandler.filename = SCEC2004matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+
+[pylithapp.timedependent.bc.y_neg]
+label = 23
+db.label = Absorbing BC -y
+db.iohandler.filename = SCEC2004matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+
+# Provide information on the fault (interface).
+[pylithapp.timedependent.interfaces]
+
+fault = pylith.faults.FaultCohesiveDyn
+
+[pylithapp.timedependent.interfaces.fault]
+
+# Material id
+id = 100
+
+# Nodeset id from CUBIT
+label = 10
+
+# Specify the basis functions and quadrature:
+# 1-D Lagrange cell in 2-D space with 2nd order quadrature
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+quadrature.cell.quad_order = 2
+
+db_initial_tractions = spatialdata.spatialdb.SimpleDB
+db_initial_tractions.iohandler.filename = SCEC2004initTract.spatialdb
+db_initial_tractions.label = "Initial fault tractions"
+
+friction = pylith.friction.SlipWeakening
+
+friction.db_properties = spatialdata.spatialdb.UniformDB
+friction.db_properties.label = Slip weakening
+friction.db_properties.values = [static-coefficient,dynamic-coefficient,slip-weakening-parameter,cohesion]
+friction.db_properties.data = [0.677,0.525,0.4*m,0.0*Pa]
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_type = gmres
+pc_type = asm
+
+# Change the preconditioner settings.
+sub_pc_factor_shift_type = nonzero
+
+ksp_rtol = 1.0e-8
+ksp_atol = 1.0e-15
+ksp_max_it = 200
+ksp_gmres_restart = 250
+snes_max_it = 300
+
+#ksp_monitor = true
+#ksp_view = true
+ksp_converged_reason = true
+
+snes_rtol = 1.0e-8
+snes_atol = 1.0e-15
+snes_monitor = true
+#snes_view = true
+snes_converged_reason = true
+
+#log_summary = true
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.output]
+writer.filename = SCEC2004.vtk
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+vertex_info_fields = [strike_dir,normal_dir,initial_traction]
+writer.filename = SCEC2004-fault.vtk
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
+writer.filename = SCEC2004-elastic.vtk
Added: short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004.exo
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004.exo
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Added: short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004.jou
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004.jou (rev 0)
+++ short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004.jou 2010-04-26 03:41:45 UTC (rev 16579)
@@ -0,0 +1,134 @@
+# ----------------------------------------------------------------------
+# Create block
+# ----------------------------------------------------------------------
+
+# Block is 7 km x 2 km x 6 km
+# -3.5 km <= x <= 3.5 km
+# -1 km <= y <= 1 km
+# -6 km <= z <= 0 km
+reset
+brick x 7000 y 2000 z 6000
+
+# Translate block so the top is at z=0
+volume 1 move x 0 y 0 z -3000
+
+# ----------------------------------------------------------------------
+# 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 with volume all
+merge all
+
+
+
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+volume all size 100
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+volume all scheme submap
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 volume 1 3
+block 1 name "elastic"
+
+# ----------------------------------------------------------------------
+# Create nodesets for faults
+# ----------------------------------------------------------------------
+group "fault" add node in surface fault_surface
+nodeset 10 group fault
+nodeset 10 name "fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in surface 10
+group "face_xpos" add node in surface 17
+nodeset 20 group face_xpos
+nodeset 20 name "face xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in surface 14
+group "face_xneg" add node in surface 12
+nodeset 21 group face_xneg
+nodeset 21 name "face xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in surface 5
+nodeset 22 group face_ypos
+nodeset 22 name "face ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in surface 3
+nodeset 23 group face_yneg
+nodeset 23 name "face yneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +z face
+# ----------------------------------------------------------------------
+group "face_zpos" add node in surface 16
+group "face_zpos" add node in surface 9
+nodeset 24 group face_zpos
+nodeset 24 name "face zpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg" add node in surface 11
+group "face_zneg" add node in surface 15
+nodeset 25 group face_zneg
+nodeset 25 name "face zneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos_nofault" add node in face_xpos
+group "face_xpos_nofault" remove node in fault
+nodeset 26 group face_xpos_nofault
+nodeset 26 name "face xpos nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg_nofault" add node in face_xneg
+group "face_xneg_nofault" remove node in fault
+nodeset 27 group face_xneg_nofault
+nodeset 27 name "face xneg nofault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -z face
+# ----------------------------------------------------------------------
+group "face_zneg_nofault" add node in face_zneg
+group "face_zneg_nofault" remove node in fault
+nodeset 28 group face_zneg_nofault
+nodeset 28 name "face zneg nofault"
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "/Users/surendranadhsomala/src/cig/pylith/playpen/friction/SCEC_benchmarks/SCEC2004.exo" dimension 3 overwrite
+
+
Added: short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004initTract.spatialdb
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004initTract.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004initTract.spatialdb 2010-04-26 03:41:45 UTC (rev 16579)
@@ -0,0 +1,37 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 3
+ value-names = traction-shear-leftlateral traction-shear-updip traction-normal
+ value-units = MPa MPa MPa
+ num-locs = 16
+ data-dim = 2
+ space-dim = 3
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // x, y, z are in km
+ space-dim = 3
+ }
+}
+-1.6 0.0 -1.4 -70.0 0.0 -120.0
+-1.4 0.0 -1.4 -70.0 0.0 -120.0
++1.4 0.0 -1.4 -70.0 0.0 -120.0
++1.6 0.0 -1.4 -70.0 0.0 -120.0
+
+-1.6 0.0 -1.6 -70.0 0.0 -120.0
+-1.4 0.0 -1.6 -81.6 0.0 -120.0
++1.4 0.0 -1.6 -81.6 0.0 -120.0
++1.6 0.0 -1.6 -70.0 0.0 -120.0
+
+-1.6 0.0 -4.4 -70.0 0.0 -120.0
+-1.4 0.0 -4.4 -81.6 0.0 -120.0
++1.4 0.0 -4.4 -81.6 0.0 -120.0
++1.6 0.0 -4.4 -70.0 0.0 -120.0
+
+-1.6 0.0 -4.6 -70.0 0.0 -120.0
+-1.4 0.0 -4.6 -70.0 0.0 -120.0
++1.4 0.0 -4.6 -70.0 0.0 -120.0
++1.6 0.0 -4.6 -70.0 0.0 -120.0
+
+
+
+
+
Added: short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004matprops.spatialdb
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004matprops.spatialdb (rev 0)
+++ short/3D/PyLith/trunk/playpen/friction/SCEC_benchmarks/SCEC2004matprops.spatialdb 2010-04-26 03:41:45 UTC (rev 16579)
@@ -0,0 +1,25 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ // number of physical properties
+ num-values = 3
+
+ // Names and units of physical properties
+ value-names = density vs vp
+ value-units = kg/m**3 m/s m/s
+
+ // Number of locations where physical properties are given
+ // 1 location -> data is uniform -> data dimension is 0
+ num-locs = 1
+ data-dim = 0
+
+ // Problem is in 3-D
+ space-dim = 3
+
+ // Coordinate system (3-D with coordinates in meters)
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 3
+ }
+}
+// Columns are x, y, z, density, vs, vp
+0.0 0.0 0.0 2670.0 3464.0 6000.0
More information about the CIG-COMMITS
mailing list