[cig-commits] r20236 - short/3D/PyLith/branches/v1.7-trunk/examples/3d/hex8

[brad at geodynamics.org]

Author: brad
Date: 2012-05-30 16:30:15 -0700 (Wed, 30 May 2012)
New Revision: 20236

Modified:
   short/3D/PyLith/branches/v1.7-trunk/examples/3d/hex8/step20.cfg
Log:
Updated settings for dike example.

Modified: short/3D/PyLith/branches/v1.7-trunk/examples/3d/hex8/step20.cfg
===================================================================
--- short/3D/PyLith/branches/v1.7-trunk/examples/3d/hex8/step20.cfg	2012-05-30 23:29:54 UTC (rev 20235)
+++ short/3D/PyLith/branches/v1.7-trunk/examples/3d/hex8/step20.cfg	2012-05-30 23:30:15 UTC (rev 20236)
@@ -6,7 +6,7 @@
 #
 # This is a pure elastic static problem that demonstrates how to use
 # initial tractions with the FaultCohesiveDyn object to create fault
-# opening consistent with a fluid intrusion, such as a dike.
+# opening consistent with fluid intrusion, such as a dike.
 #
 # ----------------------------------------------------------------------
 # RUNNING THE SIMULATON
@@ -26,9 +26,9 @@
 # problem
 # ----------------------------------------------------------------------
 [pylithapp.timedependent.formulation.time_step]
-# Define the total time for the simulation and the default time step size.
-total_time = 0.0*s ; total time of simulation
-# Define an appropriat time step for simulations. Important for
+# Define the total time for the simulation.
+total_time = 0.0*year
+# Define an appropriate time step for simulations. Important for
 # nondimensionalization of velocities and slip rates.
 dt = 1.0*year
 
@@ -95,6 +95,9 @@
 # The label corresponds to the name of the nodeset in CUBIT.
 label = fault
 
+# Apply tractions to fault surface even when it is open.
+open_free_surface = False
+
 # Use the static friction model.
 friction = pylith.friction.StaticFriction
 friction.label = Static friction
@@ -111,29 +114,15 @@
 friction.db_properties.values = [friction-coefficient,cohesion]
 friction.db_properties.data = [0.1,0.0*Pa]
 
-db_initial_tractions = spatialdata.spatialdb.SimpleDB
-db_initial_tractions.label = Initial fault tractions
-db_initial_tractions.iohandler.filename = spatialdb/tractions_opening.spatialdb
-db_initial_tractions.query_type = nearest
+# Prescribed tractions for fluid intrusion.
+traction_perturbation = pylith.faults.TractPerturbation
 
-# ----------------------------------------------------------------------
-# PETSc settings
-# ----------------------------------------------------------------------
-# NOTE: There are additional settings specific to fault friction.
-[pylithapp.petsc]
+[pylithapp.timedependent.interfaces.fault.traction_perturbation]
+db_initial = spatialdata.spatialdb.SimpleDB
+db_initial.label = Initial fault tractions
+db_initial.iohandler.filename = spatialdb/tractions_opening.spatialdb
+db_initial.query_type = nearest
 
-# Friction sensitivity solve used to compute the increment in slip
-# associated with changes in the Lagrange multiplier imposed by the
-# fault constitutive model.
-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
 # ----------------------------------------------------------------------
@@ -149,7 +138,7 @@
 
 # Give basename for VTK fault output.
 [pylithapp.problem.interfaces.fault.output]
-vertex_info_fields = [initial_traction]
+vertex_info_fields = [traction_initial_value]
 writer.filename = output/step20-fault.vtk
 
 # Give basename for VTK output of upper_crust state variables.
@@ -163,3 +152,23 @@
 # Average values over quadrature points.
 cell_filter = pylith.meshio.CellFilterAvgMesh
 writer.filename = output/step20-lower_crust.vtk
+
+
+# ----------------------------------------------------------------------
+# PETSc settings
+# ----------------------------------------------------------------------
+# NOTE: There are additional settings specific to fault friction.
+[pylithapp.petsc]
+
+# Friction sensitivity solve used to compute the increment in slip
+# associated with changes in the Lagrange multiplier imposed by the
+# fault constitutive model.
+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
+