[cig-commits] r16828 - short/3D/PyLith/trunk/examples/3d/hex8

[brad at geodynamics.org]

Author: brad
Date: 2010-05-28 21:50:24 -0700 (Fri, 28 May 2010)
New Revision: 16828

Modified:
   short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg
   short/3D/PyLith/trunk/examples/3d/hex8/step12.cfg
Log:
Adjusted step12 (static friction w/vel BC).

Modified: short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg	2010-05-29 03:46:14 UTC (rev 16827)
+++ short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg	2010-05-29 04:50:24 UTC (rev 16828)
@@ -77,7 +77,8 @@
 
 # Convergence parameters.
 ksp_rtol = 1.0e-8
-ksp_max_it = 500
+ksp_atol = 1.0e-15
+ksp_max_it = 100
 ksp_gmres_restart = 50
 
 # Linear solver monitoring options.
@@ -86,6 +87,9 @@
 ksp_converged_reason = true
 
 # Nonlinear solver monitoring options.
+snes_rtol = 1.0e-8
+snes_atol = 1.0e-15
+snes_max_it = 100
 snes_monitor = true
 snes_view = true
 snes_converged_reason = true

Modified: short/3D/PyLith/trunk/examples/3d/hex8/step12.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/step12.cfg	2010-05-29 03:46:14 UTC (rev 16827)
+++ short/3D/PyLith/trunk/examples/3d/hex8/step12.cfg	2010-05-29 04:50:24 UTC (rev 16828)
@@ -34,11 +34,8 @@
 # problem
 # ----------------------------------------------------------------------
 [pylithapp.timedependent]
-# Set bc to an array of 4 boundary conditions: 'x_pos_x', 'x_pos_y',
-# 'x_neg', and 'z_neg'. Note that 'x_pos_x' and 'x_pos_y' refer to the same
-# nodeset, because we want to apply both Neumann (x) and Dirichlet (y)
-# boundary conditions to the same set of nodes.
-bc = [x_pos_x,x_pos_y,x_neg,z_neg]
+# Set bc to an array of 3 boundary conditions: 'x_pos', 'x_neg', and 'z_neg'.
+bc = [x_pos,x_neg,z_neg]
 
 # Set interfaces to an array of 1 fault: 'fault'.
 interfaces = [fault]
@@ -67,17 +64,13 @@
 # Set the parameters for boundary conditions applied on the
 # +x, -x, and -z faces of the box.
 #
-# We apply axial (x) tractions to the +x face, and fix the x degree of
-# freedom on the -x face. We apply velocities in the y-direction on the
-# +x and -x faces, and fix the z degree of freedom on the bottom (-z) face.
+# On the -x and +x faces, we fix the x degrees of freedom and apply
+# velocities in the y-direction. We fix the z degree of freedom on the
+# bottom (-z) face.
 #
-# We first need to set the BC type to Neumann for +x, which allows us to
-# apply normal tractions on the face. We then specify Dirichlet y-velocities
-# to the same set of nodes using a different name.
-
-# We use a UniformDB to apply the tractions and velocities, while retaining
-# the default ZeroDispDB for zero x-displacements on -x and for zero
-# z-displacements on -z.
+# We use a UniformDB to apply the displacements and velocities, while
+# retaining the default ZeroDispDB for zero x-displacements on the -x
+# face and and zero z-displacements on -z.
 #
 # Note that since the fault cuts through the base of the model (z_neg),
 # we can only constrain the portion of the bottom boundary that does not
@@ -87,32 +80,20 @@
 
 # The label corresponds to the name of the nodeset in CUBIT.
 
-# +x face -- Neumann
-[pylithapp.timedependent.bc]
-x_pos_x = pylith.bc.Neumann
-
-[pylithapp.timedependent.bc.x_pos_x]
+# +x face -- Dirichlet
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0,1]
 label = face_xpos
 
-# Give quadrature information for integration of tractions.
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-quadrature.cell.quad_order = 2
-
-# Use a UniformDB to specify tractions on the face.
 db_initial = spatialdata.spatialdb.UniformDB
-db_initial.label = Neumann BC on +x
-db_initial.values = [traction-shear-horiz,traction-shear-vert,traction-normal]
-db_initial.data = [0.0*MPa,0.0*MPa,-1.0*MPa]
+db_initial.label = Dirichlet BC on +x
+db_initial.values = [displacement-x,displacement-y]
+db_initial.data = [-1.0*m,0.0*m]
 
-# +x face -- Dirichlet
-[pylithapp.timedependent.bc.x_pos_y]
-bc_dof = [1]
-label = face_xpos
 db_rate = spatialdata.spatialdb.UniformDB
 db_rate.label = Dirichlet rate BC on +x
-db_rate.values = [displacement-rate-y,rate-start-time]
-db_rate.data = [1.0*cm/year,0.0*year]
+db_rate.values = [displacement-rate-x,displacement-rate-y,rate-start-time]
+db_rate.data = [0.0*cm/year,1.0*cm/year,0.0*year]
 
 # -x face
 [pylithapp.timedependent.bc.x_neg]
@@ -154,7 +135,7 @@
 friction.db_properties = spatialdata.spatialdb.UniformDB
 friction.db_properties.label = Static friction
 friction.db_properties.values = [friction-coefficient,cohesion]
-friction.db_properties.data = [0.6,0.0*Pa]
+friction.db_properties.data = [0.6,2.0*MPa]
 
 # ----------------------------------------------------------------------
 # PETSc settings