[cig-commits] r16598 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d
brad at geodynamics.org
brad at geodynamics.org
Wed Apr 28 15:39:01 PDT 2010
Author: brad
Date: 2010-04-28 15:39:00 -0700 (Wed, 28 Apr 2010)
New Revision: 16598
Added:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/createbc.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/friction.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/matprops.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/output/
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m.exo.gz
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tractions.spatialdb
Removed:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/TPV5_2D.cfg
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/TPV5_2D.exo
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/TPV5_2D.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/frictpar_TPV5_2D.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/initTract_TPV5_2D.spatialdb
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/matprops2D.spatialdb
Log:
Adjusted dimensions of domain. Limited fault size to region of interest. Adjusted parameters.
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/TPV5_2D.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/TPV5_2D.cfg 2010-04-28 22:09:59 UTC (rev 16597)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/TPV5_2D.cfg 2010-04-28 22:39:00 UTC (rev 16598)
@@ -1,207 +0,0 @@
-# -*- Python -*-
-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# journal
-# ----------------------------------------------------------------------
-[pylithapp.journal.info]
-timedependent = 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 = TPV5_2D.exo
-coordsys.space_dim = 2
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-
-dimension = 2
-
-# 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,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.y_neg = pylith.bc.AbsorbingDampers
-bc.y_pos = pylith.bc.AbsorbingDampers
-
-[pylithapp.timedependent.formulation.time_step]
-total_time = 1.0*s
-dt = 0.1*s
-
-
-# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-
-materials = [elastic]
-materials.elastic = pylith.materials.ElasticPlaneStrain
-
-[pylithapp.timedependent.materials.elastic]
-label = Elastic material
-id = 1
-db_properties.iohandler.filename = matprops2D.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_pos]
-label = 11
-db.label = Absorbing BC +x
-db.iohandler.filename = matprops2D.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
-
-[pylithapp.timedependent.bc.x_neg]
-label = 12
-db.label = Absorbing BC -x
-db.iohandler.filename = matprops2D.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
-
-[pylithapp.timedependent.bc.y_pos]
-label = 13
-db.label = Absorbing BC +y
-db.iohandler.filename = matprops2D.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
-
-[pylithapp.timedependent.bc.y_neg]
-label = 14
-db.label = Absorbing BC -y
-db.iohandler.filename = matprops2D.spatialdb
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 1
-quadrature.cell.quad_order = 2
-
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-
-# Provide information on the fault (interface).
-[pylithapp.timedependent.interfaces]
-
-fault = pylith.faults.FaultCohesiveDyn
-
-[pylithapp.timedependent.interfaces.fault]
-
-# Material id
-id = 10
-
-# 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 = 1
-quadrature.cell.quad_order = 2
-
-db_initial_tractions = spatialdata.spatialdb.SimpleDB
-db_initial_tractions.iohandler.filename = initTract_TPV5_2D.spatialdb
-db_initial_tractions.label = "Initial fault tractions"
-
-friction = pylith.friction.SlipWeakening
-
-friction.db_properties = spatialdata.spatialdb.SimpleDB
-friction.db_properties.iohandler.filename = frictpar_TPV5_2D.spatialdb
-friction.db_properties.label = Slip weakening
-
-# ----------------------------------------------------------------------
-# 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
-
-# Friction sensitivity solve
-friction_pc_type = asm
-friction_sub_pc_factor_shift_type = nonzero
-
-friction_ksp_max_it = 25
-friction_ksp_gmres_restart = 30
-
-#friction_ksp_monitor = true
-#friction_ksp_view = true
-friction_ksp_converged_reason = true
-
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
-[pylithapp.problem.formulation.output.output]
-writer.filename = output/TPV5_2D.vtk
-
-# Give basename for VTK fault output.
-[pylithapp.timedependent.interfaces.fault.output]
-vertex_info_fields = [strike_dir,normal_dir,initial_traction]
-writer.filename = output/TPV5_2D-fault.vtk
-
-# Give basename for VTK output of state variables.
-[pylithapp.timedependent.materials.elastic.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = output/TPV5_2D-elastic.vtk
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/TPV5_2D.exo
===================================================================
(Binary files differ)
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/TPV5_2D.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/TPV5_2D.jou 2010-04-28 22:09:59 UTC (rev 16597)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/TPV5_2D.jou 2010-04-28 22:39:00 UTC (rev 16598)
@@ -1,108 +0,0 @@
-# ----------------------------------------------------------------------
-# Create block
-# ----------------------------------------------------------------------
-# Bar is 20.4 km x 40.8 km
-# -10.2 km <= x <= 10.2 km
-# -20.4 km <= y <= 20.4 km
-#
-# We create the 2-D TPV5 by extracting the mid-surface from a 3-D TPV5.
-reset
-brick x 20400 y 40800 z 20400
-create midsurface volume 1 surface 1 2
-delete volume 1
-
-# ----------------------------------------------------------------------
-# Create interface surfaces
-# ----------------------------------------------------------------------
-
-# Create fault (yz plane) at x = 0.0
-split surface 7 across location position 0 -20400 0 location position 0 20400 0
-curve 17 name "fault_surface"
-
-# ----------------------------------------------------------------------
-# Imprint all volumes, merging surfaces
-# ----------------------------------------------------------------------
-imprint all
-merge all
-
-
-
-
-
-
-
-# ----------------------------------------------------------------------
-# Set discretization size
-# ----------------------------------------------------------------------
-surface all size 200
-
-# ----------------------------------------------------------------------
-# Generate the mesh
-# ----------------------------------------------------------------------
-mesh surface all
-
-# ----------------------------------------------------------------------
-# Create blocks for materials
-# ----------------------------------------------------------------------
-block 1 surface 8 9
-block 1 name "elastic"
-
-# ----------------------------------------------------------------------
-# Create nodeset for fault
-# ----------------------------------------------------------------------
-group "fault" add node in fault_surface
-nodeset 10 group fault
-nodeset 10 name "fault"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +x face
-# ----------------------------------------------------------------------
-group "face_xpos" add node in curve 14
-nodeset 11 group face_xpos
-nodeset 11 name "face xpos"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -x face
-# ----------------------------------------------------------------------
-group "face_xneg" add node in curve 16
-nodeset 12 group face_xneg
-nodeset 12 name "face xneg"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +y face
-# ----------------------------------------------------------------------
-group "face_ypos" add node in curve 18
-group "face_ypos" add node in curve 21
-nodeset 13 group face_ypos
-nodeset 13 name "face ypos"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -y face
-# ----------------------------------------------------------------------
-group "face_yneg" add node in curve 19
-group "face_yneg" add node in curve 20
-nodeset 14 group face_yneg
-nodeset 14 name "face yneg"
-
-# ----------------------------------------------------------------------
-# Create nodeset for +y face without fault
-# ----------------------------------------------------------------------
-group "face_ypos" add node in curve 18
-group "face_ypos" add node in curve 21
-group "face_ypos" remove node in fault
-nodeset 15 group face_ypos
-nodeset 15 name "face ypos"
-
-# ----------------------------------------------------------------------
-# Create nodeset for -y face without fault
-# ----------------------------------------------------------------------
-group "face_yneg" add node in curve 19
-group "face_yneg" add node in curve 20
-group "face_yneg" remove node in fault
-nodeset 16 group face_yneg
-nodeset 16 name "face yneg"
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "TPV5_2D.exo" dimension 2 overwrite
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/createbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/createbc.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/createbc.jou 2010-04-28 22:39:00 UTC (rev 16598)
@@ -0,0 +1,44 @@
+# ----------------------------------------------------------------------
+# Create blocks for materials
+# ----------------------------------------------------------------------
+block 1 surface 8 9
+block 1 name "elastic"
+
+# ----------------------------------------------------------------------
+# Create nodeset for fault
+# ----------------------------------------------------------------------
+group "fault" add node in fault_surface
+group "fault" remove node with y_coord > +15.0001e+3
+group "fault" remove node with y_coord < -15.0001e+3
+nodeset 10 group fault
+nodeset 10 name "fault"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +x face
+# ----------------------------------------------------------------------
+group "face_xpos" add node in curve 14
+nodeset 11 group face_xpos
+nodeset 11 name "face xpos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -x face
+# ----------------------------------------------------------------------
+group "face_xneg" add node in curve 16
+nodeset 12 group face_xneg
+nodeset 12 name "face xneg"
+
+# ----------------------------------------------------------------------
+# Create nodeset for +y face
+# ----------------------------------------------------------------------
+group "face_ypos" add node in curve 18
+group "face_ypos" add node in curve 21
+nodeset 13 group face_ypos
+nodeset 13 name "face ypos"
+
+# ----------------------------------------------------------------------
+# Create nodeset for -y face
+# ----------------------------------------------------------------------
+group "face_yneg" add node in curve 19
+group "face_yneg" add node in curve 20
+nodeset 14 group face_yneg
+nodeset 14 name "face yneg"
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/friction.spatialdb (from rev 16594, short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/frictpar_TPV5_2D.spatialdb)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/friction.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/friction.spatialdb 2010-04-28 22:39:00 UTC (rev 16598)
@@ -0,0 +1,17 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 4
+ value-names = static-coefficient dynamic-coefficient slip-weakening-parameter cohesion
+ value-units = none none m MPa
+ num-locs = 4
+ data-dim = 2
+ space-dim = 2
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // x, y are in km
+ space-dim = 2
+ }
+}
+0.0 -15.1 1000.0 0.525 0.4 0.0
+0.0 -14.9 0.677 0.525 0.4 0.0
+0.0 +14.9 0.677 0.525 0.4 0.0
+0.0 +15.1 1000.0 0.525 0.4 0.0
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/frictpar_TPV5_2D.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/frictpar_TPV5_2D.spatialdb 2010-04-28 22:09:59 UTC (rev 16597)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/frictpar_TPV5_2D.spatialdb 2010-04-28 22:39:00 UTC (rev 16598)
@@ -1,17 +0,0 @@
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 4
- value-names = static-coefficient dynamic-coefficient slip-weakening-parameter cohesion
- value-units = none none m MPa
- num-locs = 4
- data-dim = 2
- space-dim = 2
- cs-data = cartesian {
- to-meters = 1.0e+3 // x, y are in km
- space-dim = 2
- }
-}
-0.0 -15.1 1000.0 0.525 0.4 0.0
-0.0 -14.9 0.677 0.525 0.4 0.0
-0.0 +14.9 0.677 0.525 0.4 0.0
-0.0 +15.1 1000.0 0.525 0.4 0.0
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/geometry.jou 2010-04-28 22:39:00 UTC (rev 16598)
@@ -0,0 +1,26 @@
+# ----------------------------------------------------------------------
+# Create block
+# ----------------------------------------------------------------------
+# Block is 24.0 km x 48.0 km
+# -12.0 km <= x <= 12.0 km
+# -24.0 km <= y <= 24.0 km
+#
+# We create the 2-D mesh by extracting the mid-surface from a 3-D mesh.
+reset
+brick x 24000 y 48000 z 24000
+create midsurface volume 1 surface 1 2
+delete volume 1
+
+# ----------------------------------------------------------------------
+# Create interface surfaces
+# ----------------------------------------------------------------------
+
+# Create fault (yz plane) at x = 0.0
+split surface 7 across location position 0 -24000 0 location position 0 24000 0
+curve 17 name "fault_surface"
+
+# ----------------------------------------------------------------------
+# Imprint all volumes, merging surfaces
+# ----------------------------------------------------------------------
+imprint all
+merge all
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/initTract_TPV5_2D.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/initTract_TPV5_2D.spatialdb 2010-04-28 22:09:59 UTC (rev 16597)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/initTract_TPV5_2D.spatialdb 2010-04-28 22:39:00 UTC (rev 16598)
@@ -1,25 +0,0 @@
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 2
- value-names = traction-shear traction-normal
- value-units = MPa MPa
- num-locs = 12
- data-dim = 2
- space-dim = 2
- cs-data = cartesian {
- to-meters = 1.0e+3 // x, y, z are in km
- space-dim = 2
- }
-}
-0.0 -9.1 -70.0 -120.0
-0.0 -8.9 -78.0 -120.0
-0.0 -6.1 -78.0 -120.0
-0.0 -5.9 -70.0 -120.0
-0.0 -1.6 -70.0 -120.0
-0.0 -1.4 -81.6 -120.0
-0.0 +1.4 -81.6 -120.0
-0.0 +1.6 -70.0 -120.0
-0.0 +5.9 -70.0 -120.0
-0.0 +6.1 -62.0 -120.0
-0.0 +8.9 -62.0 -120.0
-0.0 +9.1 -70.0 -120.0
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/matprops.spatialdb (from rev 16594, short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/matprops2D.spatialdb)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/matprops.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/matprops.spatialdb 2010-04-28 22:39:00 UTC (rev 16598)
@@ -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 2-D
+ space-dim = 2
+
+ // Coordinate system (2-D with coordinates in meters)
+ cs-data = cartesian {
+ to-meters = 1.0
+ space-dim = 2
+ }
+}
+// Columns are x, y, density, vs, vp
+0.0 0.0 2670.0 3464.0 6000.0
Deleted: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/matprops2D.spatialdb
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/matprops2D.spatialdb 2010-04-28 22:09:59 UTC (rev 16597)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/matprops2D.spatialdb 2010-04-28 22:39:00 UTC (rev 16598)
@@ -1,25 +0,0 @@
-#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 2-D
- space-dim = 2
-
- // Coordinate system (2-D with coordinates in meters)
- cs-data = cartesian {
- to-meters = 1.0
- space-dim = 2
- }
-}
-// Columns are x, y, density, vs, vp
-0.0 0.0 2670.0 3464.0 6000.0
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/pylithapp.cfg (from rev 16594, short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/TPV5_2D.cfg)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/pylithapp.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/pylithapp.cfg 2010-04-28 22:39:00 UTC (rev 16598)
@@ -0,0 +1,188 @@
+# -*- Python -*-
+
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+explicit = 1
+petsc = 1
+meshiocubit = 1
+elasticityexplicit = 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]
+coordsys.space_dim = 2
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+
+dimension = 2
+
+# 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,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.y_neg = pylith.bc.AbsorbingDampers
+bc.y_pos = pylith.bc.AbsorbingDampers
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 1.00*s
+dt = 0.01*s
+
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+
+materials = [elastic]
+materials.elastic = pylith.materials.ElasticPlaneStrain
+
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
+id = 1
+db_properties.iohandler.filename = matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+label = 11
+db.label = Absorbing BC +x
+db.iohandler.filename = matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.x_neg]
+label = 12
+db.label = Absorbing BC -x
+db.iohandler.filename = matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_pos]
+label = 13
+db.label = Absorbing BC +y
+db.iohandler.filename = matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_neg]
+label = 14
+db.label = Absorbing BC -y
+db.iohandler.filename = matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+quadrature.cell.quad_order = 2
+
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+
+# Provide information on the fault (interface).
+[pylithapp.timedependent.interfaces]
+
+fault = pylith.faults.FaultCohesiveDyn
+
+[pylithapp.timedependent.interfaces.fault]
+
+# Material id
+id = 10
+
+# 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 = 1
+quadrature.cell.quad_order = 2
+
+db_initial_tractions = spatialdata.spatialdb.SimpleDB
+db_initial_tractions.iohandler.filename = tractions.spatialdb
+db_initial_tractions.label = "Initial fault tractions"
+
+friction = pylith.friction.SlipWeakening
+
+friction.db_properties = spatialdata.spatialdb.SimpleDB
+friction.db_properties.iohandler.filename = friction.spatialdb
+friction.db_properties.label = Slip weakening
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+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 = 100
+ksp_gmres_restart = 50
+snes_max_it = 200
+
+ksp_monitor = true
+ksp_view = true
+ksp_converged_reason = true
+
+snes_monitor = true
+snes_view = true
+snes_converged_reason = true
+
+log_summary = true
+
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+
+friction_ksp_monitor = true
+friction_ksp_view = true
+friction_ksp_converged_reason = true
+
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m.cfg 2010-04-28 22:39:00 UTC (rev 16598)
@@ -0,0 +1,35 @@
+# -*- Python -*-
+
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = quad4_200m.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.output]
+output_freq = time_step
+time_step = 0.05*s
+writer.filename = output/quad4_200m.vtk
+writer.time_format = %05.2f
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+vertex_info_fields = [strike_dir,normal_dir,initial_traction]
+output_freq = time_step
+time_step = 0.05*s
+writer.filename = output/quad4_200m-fault.vtk
+writer.time_format = %05.2f
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
+output_freq = time_step
+time_step = 1.0*s
+writer.filename = output/quad4_200m-elastic.vtk
+writer.time_format = %05.2f
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m.exo.gz
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m.exo.gz
___________________________________________________________________
Name: svn:mime-type
+ application/x-gzip
Added: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/quad4_200m.jou 2010-04-28 22:39:00 UTC (rev 16598)
@@ -0,0 +1,29 @@
+# ----------------------------------------------------------------------
+# Create quad4 mesh at 200m resolution.
+# ----------------------------------------------------------------------
+
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+surface all size 200
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh surface all
+
+# ----------------------------------------------------------------------
+# Setup boundary conditions.
+# ----------------------------------------------------------------------
+playback 'createbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "quad4_200m.exo" dimension 2 overwrite
+
Copied: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tractions.spatialdb (from rev 16594, short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/initTract_TPV5_2D.spatialdb)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tractions.spatialdb (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv205-2d/tractions.spatialdb 2010-04-28 22:39:00 UTC (rev 16598)
@@ -0,0 +1,25 @@
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 2
+ value-names = traction-shear traction-normal
+ value-units = MPa MPa
+ num-locs = 12
+ data-dim = 2
+ space-dim = 2
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // x, y, z are in km
+ space-dim = 2
+ }
+}
+0.0 -9.1 -70.0 -120.0
+0.0 -8.9 -78.0 -120.0
+0.0 -6.1 -78.0 -120.0
+0.0 -5.9 -70.0 -120.0
+0.0 -1.6 -70.0 -120.0
+0.0 -1.4 -81.6 -120.0
+0.0 +1.4 -81.6 -120.0
+0.0 +1.6 -70.0 -120.0
+0.0 +5.9 -70.0 -120.0
+0.0 +6.1 -62.0 -120.0
+0.0 +8.9 -62.0 -120.0
+0.0 +9.1 -70.0 -120.0
More information about the CIG-COMMITS
mailing list