[cig-commits] r16685 - short/3D/PyLith/trunk/examples/3d/hex8
brad at geodynamics.org
brad at geodynamics.org
Tue May 11 09:36:14 PDT 2010
Author: brad
Date: 2010-05-11 09:36:14 -0700 (Tue, 11 May 2010)
New Revision: 16685
Modified:
short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg
short/3D/PyLith/trunk/examples/3d/hex8/step01.cfg
Log:
Moved Maxwell spatial database info from step01.cfg to pylithapp.cfg. Cleanup of some annotation.
Modified: short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg 2010-05-11 00:00:46 UTC (rev 16684)
+++ short/3D/PyLith/trunk/examples/3d/hex8/pylithapp.cfg 2010-05-11 16:36:14 UTC (rev 16685)
@@ -16,8 +16,6 @@
fiatlagrange = 1
pylithapp = 1
materials = 1
-powerlaw3d = 1
-##implicit = 1
# ----------------------------------------------------------------------
# mesh_generator
@@ -59,11 +57,10 @@
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 3
-# We do not define the database info for the inelastic material, since
-# a different database type may be used for some examples.
[pylithapp.timedependent.materials.inelastic]
label = Inelastic material
id = 2
+db_properties.iohandler.filename = spatialdb/mat_maxwell.spatialdb
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 3
@@ -83,18 +80,18 @@
ksp_max_it = 500
ksp_gmres_restart = 50
-# Linear solution logging options.
+# Linear solver monitoring options.
ksp_monitor = true
ksp_view = true
ksp_converged_reason = true
-# Nonlinear solution logging options.
+# Nonlinear solver monitoring options.
snes_monitor = true
snes_view = true
snes_converged_reason = true
-# PETSc run summary -- useful for performance information.
+# PETSc summary -- useful for performance information.
log_summary = true
-# Hopefully not necessary.
+# Uncomment to launch gdb when starting PyLith.
# start_in_debugger = true
Modified: short/3D/PyLith/trunk/examples/3d/hex8/step01.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/step01.cfg 2010-05-11 00:00:46 UTC (rev 16684)
+++ short/3D/PyLith/trunk/examples/3d/hex8/step01.cfg 2010-05-11 16:36:14 UTC (rev 16685)
@@ -1,6 +1,17 @@
# -*- Python -*-
[pylithapp]
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+
+#
+# [ADD PHYSICAL DESCRIPTION OF THE PROBLEM BEING SOLVED]
+
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+
# This is not a self-contained simulation configuration file. This
# file only specifies parameters specific to tutorial step01.
# The general parameters are specificed in the pylithapp.cfg
@@ -9,7 +20,7 @@
# To run the simulation:
# pylith step01.cfg
#
-# Output will be directed to directory step01.
+# Output will be directed to the directory step01.
# ----------------------------------------------------------------------
# problem
@@ -23,36 +34,29 @@
# We will output the solution over the domain and the ground surface.
output = [domain,subdomain]
-# Set subdomain component to OutputSolnSubset (subset of domain).
+# Set subdomain component to OutputSolnSubset (boundary of the domain).
output.subdomain = pylith.meshio.OutputSolnSubset
# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-
-# We specifically define the database info for the inelastic material
-# for each model, since the type of spatialdb can be different.
-
-[pylithapp.timedependent.materials.inelastic]
-db_properties.iohandler.filename = spatialdb/mat_maxwell.spatialdb
-
-# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
+# Set the parameters for Dirichlet boundary conditions applied on the
+# +x, -x, and -z faces of the box.
+#
+# We fix the x and y degrees of freedom on the +x and -x faces, and
+# fix the z degree of freedom on the bottom (-z) face.
+#
# For x_pos and x_neg, we change the database from ZeroDispDB (which
-# specifies a zero value everywhere) to SimpleDB (which uses a data file
-# to specify a spatial variation).
+# specifies a zero value) to SimpleDB (which uses a data file to
+# specify a spatial variation).
+#
# For z_neg, we keep the default ZeroDispDB, since we just want zero
# z-displacements on that boundary.
+#
-# Set the parameters for Dirichlet boundary conditions applied on the
-# +x and -x faces of the box.
-
-# We fix the x and y degrees of freedom on the +x and -x faces, and
-# fix the z degree of freedom on the bottom face.
-
# The label corresponds to the nodeset ID in CUBIT.
+# +x face
[pylithapp.timedependent.bc.x_pos]
bc_dof = [0, 1]
label = 11
@@ -60,6 +64,7 @@
db_initial.label = Dirichlet BC on +x
db_initial.iohandler.filename = spatialdb/fixeddisp_axial_shear.spatialdb
+# -x face
[pylithapp.timedependent.bc.x_neg]
bc_dof = [0, 1]
label = 12
@@ -67,6 +72,7 @@
db_initial.label = Dirichlet BC on -x
db_initial.iohandler.filename = spatialdb/fixeddisp_axial_shear.spatialdb
+# -z face
[pylithapp.timedependent.bc.z_neg]
bc_dof = [2]
label = 15
@@ -81,11 +87,13 @@
# Give basename for VTK domain output of solution over ground surface.
[pylithapp.problem.formulation.output.subdomain]
-label = 17 ; nodeset for subdomain
+# nodeset for boundary of domain
+label = 17
writer.filename = step01/step01-groundsurf.vtk
# Give basename for VTK output of elastic state variables.
[pylithapp.timedependent.materials.elastic.output]
+# Average values over quadrature points.
cell_filter = pylith.meshio.CellFilterAvgMesh
writer.filename = step01/step01-elastic.vtk
@@ -95,5 +103,6 @@
# information, so we include it here.
cell_info_fields = [density,mu,lambda,maxwell_time]
cell_data_fields = [total_strain,stress,viscous_strain]
+# Average values over quadrature points.
cell_filter = pylith.meshio.CellFilterAvgMesh
writer.filename = step01/step01-inelastic.vtk
More information about the CIG-COMMITS
mailing list