[cig-commits] r16716 - short/3D/PyLith/trunk/examples/bar_shearwave/quad4
brad at geodynamics.org
brad at geodynamics.org
Sat May 15 11:25:37 PDT 2010
Author: brad
Date: 2010-05-15 11:25:37 -0700 (Sat, 15 May 2010)
New Revision: 16716
Modified:
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_ratestateageing.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_slipweakening.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_staticfriction.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_sliding.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_stick.cfg
Log:
Finished setting up .cfg files.
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault.cfg 2010-05-15 18:25:13 UTC (rev 16715)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/dynamicfault.cfg 2010-05-15 18:25:37 UTC (rev 16716)
@@ -12,6 +12,63 @@
faultcohesivedyn = 1
# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+
+# Switch to a lumped Jacobian matrix.
+formulation = pylith.problems.ExplicitLumped
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc]
+# Change the DirichletPoints BC to the AbsorbingDampers BC
+x_pos = pylith.bc.AbsorbingDampers
+x_neg = pylith.bc.AbsorbingDampers
+
+[pylithapp.timedependent.bc.x_pos]
+# Absorbing BC on +x face of bar
+
+# Nodeset id from CUBIT
+label = 11
+
+# Specify label and filename for spatial database with physical properties
+db.label = Absorbing BC +x
+db.iohandler.filename = matprops.spatialdb
+
+# 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
+
+[pylithapp.timedependent.bc.x_neg]
+# Absorbing BC on -x face of bar
+
+# Nodeset id from CUBIT
+label = 12
+
+# Specify label and filename for spatial database with physical properties
+db.label = Absorbing BC -x
+db.iohandler.filename = matprops.spatialdb
+
+# 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
+
+[pylithapp.timedependent.bc.all_nofault]
+# Dirichlet BC on all vertices except fault vertices
+
+# Pin x and z degrees of freedom
+bc_dof = [0]
+
+# Nodeset id from CUBIT
+label = 13
+
+# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
[pylithapp.timedependent.interfaces]
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_ratestateageing.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_ratestateageing.cfg 2010-05-15 18:25:13 UTC (rev 16715)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_ratestateageing.cfg 2010-05-15 18:25:37 UTC (rev 16716)
@@ -1,99 +1,25 @@
# -*- Python -*-
-
[pylithapp]
# ----------------------------------------------------------------------
-# problem
+# RUNNING THE SIMULATON
# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-formulation = pylith.problems.Explicit
-formulation.solver = pylith.problems.SolverNonlinear
-#formulation = pylith.problems.ExplicitLumped
+# This is not a self-contained simulation configuration file. This
+# file only specifies the specific parameters for an explicit
+# time-stepping problem using a dynamic fault interface condition with
+# static friction.
+#
+# To run the simulation:
+# pylith dynamicfault.cfg shearwave_staticfriction.cfg
+#
+# Output will be directed to the directory output.
-# 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_pos,x_neg,all_nofault]
-
-# Change the DirichletPoints BC to the AbsorbingDampers BC
-bc.x_pos = pylith.bc.AbsorbingDampers
-bc.x_neg = pylith.bc.AbsorbingDampers
-
-[pylithapp.timedependent.formulation.time_step]
-total_time = 12.0*s
-dt = 0.05*s
-
# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_pos]
-# Absorbing BC on +x face of bar
-
-# Nodeset id from CUBIT
-label = 11
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC +x
-db.iohandler.filename = matprops.spatialdb
-
-# 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
-
-[pylithapp.timedependent.bc.x_neg]
-# Absorbing BC on -x face of bar
-
-# Nodeset id from CUBIT
-label = 12
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC -x
-db.iohandler.filename = matprops.spatialdb
-
-# 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
-
-[pylithapp.timedependent.bc.all_nofault]
-# Dirichlet BC on all vertices except fault vertices
-
-# Pin x and z degrees of freedom
-bc_dof = [0]
-
-# Nodeset id from CUBIT
-label = 13
-
-# ----------------------------------------------------------------------
# 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 = 1
-quadrature.cell.quad_order = 2
-
db_initial_tractions = spatialdata.spatialdb.UniformDB
db_initial_tractions.label = Initial fault tractions
db_initial_tractions.values = [traction-shear,traction-normal]
@@ -114,14 +40,12 @@
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
+# Set filenames for output.
[pylithapp.problem.formulation.output.output]
writer.filename = output/shearwave-ratestateageing.vtk
-# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
writer.filename = output/shearwave-ratestateageing-fault.vtk
-# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.elastic.output]
writer.filename = output/shearwave-ratestateageing-statevars.vtk
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_slipweakening.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_slipweakening.cfg 2010-05-15 18:25:13 UTC (rev 16715)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_slipweakening.cfg 2010-05-15 18:25:37 UTC (rev 16716)
@@ -1,99 +1,25 @@
# -*- Python -*-
-
[pylithapp]
# ----------------------------------------------------------------------
-# problem
+# RUNNING THE SIMULATON
# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-formulation = pylith.problems.Explicit
-formulation.solver = pylith.problems.SolverNonlinear
-#formulation = pylith.problems.ExplicitLumped
+# This is not a self-contained simulation configuration file. This
+# file only specifies the specific parameters for an explicit
+# time-stepping problem using a dynamic fault interface condition with
+# linear slip-weakening friction.
+#
+# To run the simulation:
+# pylith dynamicfault.cfg shearwave_slipweakening.cfg
+#
+# Output will be directed to the directory output.
-# 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_pos,x_neg,all_nofault]
-
-# Change the DirichletPoints BC to the AbsorbingDampers BC
-bc.x_pos = pylith.bc.AbsorbingDampers
-bc.x_neg = pylith.bc.AbsorbingDampers
-
-[pylithapp.timedependent.formulation.time_step]
-total_time = 1.0*s
-dt = 0.05*s
-
# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_pos]
-# Absorbing BC on +x face of bar
-
-# Nodeset id from CUBIT
-label = 11
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC +x
-db.iohandler.filename = matprops.spatialdb
-
-# 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
-
-[pylithapp.timedependent.bc.x_neg]
-# Absorbing BC on -x face of bar
-
-# Nodeset id from CUBIT
-label = 12
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC -x
-db.iohandler.filename = matprops.spatialdb
-
-# 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
-
-[pylithapp.timedependent.bc.all_nofault]
-# Dirichlet BC on all vertices except fault vertices
-
-# Pin x and z degrees of freedom
-bc_dof = [0]
-
-# Nodeset id from CUBIT
-label = 13
-
-# ----------------------------------------------------------------------
# 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 = 1
-quadrature.cell.quad_order = 2
-
db_initial_tractions = spatialdata.spatialdb.UniformDB
db_initial_tractions.label = Initial fault tractions
db_initial_tractions.values = [traction-shear,traction-normal]
@@ -109,16 +35,20 @@
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
[pylithapp.problem.formulation.output.output]
writer.filename = output/shearwave-slipweakening.vtk
-# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
+writer.filename = output/shearwave-slipweakening-fault.vtk
+
+# Add friction model parameters to fault info file in addition to
+# default values.
vertex_info_fields = [strike_dir,normal_dir,initial_traction,static_coefficient,dynamic_coefficient,slip_weakening_parameter,cohesion]
+
+# Add output of state variables cumulative slip and previous slip to default
+# values.
vertex_data_fields = [slip,traction,cumulative_slip,previous_slip]
-writer.filename = output/shearwave-slipweakening-fault.vtk
-# Give basename for VTK output of state variables.
+
[pylithapp.timedependent.materials.elastic.output]
writer.filename = output/shearwave-slipweakening-statevars.vtk
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_staticfriction.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_staticfriction.cfg 2010-05-15 18:25:13 UTC (rev 16715)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/shearwave_staticfriction.cfg 2010-05-15 18:25:37 UTC (rev 16716)
@@ -16,63 +16,6 @@
# Output will be directed to the directory output.
# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-
-# Switch to a lumped Jacobian matrix.
-formulation = pylith.problems.ExplicitLumped
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc]
-# Change the DirichletPoints BC to the AbsorbingDampers BC
-x_pos = pylith.bc.AbsorbingDampers
-x_neg = pylith.bc.AbsorbingDampers
-
-[pylithapp.timedependent.bc.x_pos]
-# Absorbing BC on +x face of bar
-
-# Nodeset id from CUBIT
-label = 11
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC +x
-db.iohandler.filename = matprops.spatialdb
-
-# 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
-
-[pylithapp.timedependent.bc.x_neg]
-# Absorbing BC on -x face of bar
-
-# Nodeset id from CUBIT
-label = 12
-
-# Specify label and filename for spatial database with physical properties
-db.label = Absorbing BC -x
-db.iohandler.filename = matprops.spatialdb
-
-# 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
-
-[pylithapp.timedependent.bc.all_nofault]
-# Dirichlet BC on all vertices except fault vertices
-
-# Pin x and z degrees of freedom
-bc_dof = [0]
-
-# Nodeset id from CUBIT
-label = 13
-
-# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
[pylithapp.timedependent.interfaces.fault]
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_sliding.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_sliding.cfg 2010-05-15 18:25:13 UTC (rev 16715)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_sliding.cfg 2010-05-15 18:25:37 UTC (rev 16716)
@@ -10,11 +10,17 @@
# with fault slip controlled by a fault constitutive model.
#
# To run the simulation:
-# pylith dynamicfault.cfg static_sliding.cfg
+# pylith static_sliding.cfg
#
# Output will be directed to the directory output.
# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+faultcohesivedyn = 1
+
+# ----------------------------------------------------------------------
# problem
# ----------------------------------------------------------------------
[pylithapp.timedependent]
@@ -77,17 +83,32 @@
friction.db_properties.data = [0.6*m,0.0*Pa]
# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+# Uncomment to view details of friction sensitivity solve.
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
+# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
-
-# Give basename for VTK domain output of solution over domain.
+# Set filenames for output.
[pylithapp.problem.formulation.output.output]
writer.filename = output/static_sliding.vtk
-# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
writer.filename = output/static_sliding-fault.vtk
-# Give basename for VTK output of state variables.
+# Add slip rate to the default output fiels of slip and traction.
+vertex_data_fields=[slip,slip_rate,traction]
+
[pylithapp.timedependent.materials.elastic.output]
writer.filename = output/static_sliding-statevars.vtk
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_stick.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_stick.cfg 2010-05-15 18:25:13 UTC (rev 16715)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/static_stick.cfg 2010-05-15 18:25:37 UTC (rev 16716)
@@ -10,11 +10,17 @@
# with a locked fault controlled by a fault constitutive model.
#
# To run the simulation:
-# pylith dynamicfault.cfg static-stick.cfg
+# pylith static-stick.cfg
#
# Output will be directed to the directory output.
# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+faultcohesivedyn = 1
+
+# ----------------------------------------------------------------------
# problem
# ----------------------------------------------------------------------
[pylithapp.timedependent]
@@ -61,6 +67,10 @@
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces]
+# Change fault to dynamic fault interface.
+fault = pylith.faults.FaultCohesiveDyn
+
[pylithapp.timedependent.interfaces.fault]
db_initial_tractions = spatialdata.spatialdb.UniformDB
@@ -76,17 +86,32 @@
friction.db_properties.data = [0.6*m,0.0*Pa]
# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+# Uncomment to view details of friction sensitivity solve.
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
+# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
-
-# Give basename for VTK domain output of solution over domain.
+# Set filenames for output.
[pylithapp.problem.formulation.output.output]
writer.filename = output/static_stick.vtk
-# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
writer.filename = output/static_stick-fault.vtk
-# Give basename for VTK output of state variables.
+# Add slip rate to the default output fiels of slip and traction.
+vertex_data_fields=[slip,slip_rate,traction]
+
[pylithapp.timedependent.materials.elastic.output]
writer.filename = output/static_stick-statevars.vtk
More information about the CIG-COMMITS
mailing list