[cig-commits] r16945 - in short/3D/PyLith/trunk/doc/userguide: materials tutorials/3dhex8 tutorials/3dhex8/quasistatic tutorials/3dhex8/quasistatic/figs tutorials/3dhex8/static
willic3 at geodynamics.org
willic3 at geodynamics.org
Tue Jun 8 16:40:00 PDT 2010
Author: willic3
Date: 2010-06-08 16:40:00 -0700 (Tue, 08 Jun 2010)
New Revision: 16945
Added:
short/3D/PyLith/trunk/doc/userguide/tutorials/3dhex8/quasistatic/figs/step07-displ-vel-t300.jpg
Modified:
short/3D/PyLith/trunk/doc/userguide/materials/materials.lyx
short/3D/PyLith/trunk/doc/userguide/tutorials/3dhex8/3dhex8.lyx
short/3D/PyLith/trunk/doc/userguide/tutorials/3dhex8/quasistatic/quasistatic.lyx
short/3D/PyLith/trunk/doc/userguide/tutorials/3dhex8/static/static.lyx
Log:
More work on tutorials and made some changes related to power-law.
Modified: short/3D/PyLith/trunk/doc/userguide/materials/materials.lyx
===================================================================
--- short/3D/PyLith/trunk/doc/userguide/materials/materials.lyx 2010-06-08 23:34:40 UTC (rev 16944)
+++ short/3D/PyLith/trunk/doc/userguide/materials/materials.lyx 2010-06-08 23:40:00 UTC (rev 16945)
@@ -1,4 +1,4 @@
-#LyX 1.6.5 created this file. For more info see http://www.lyx.org/
+#LyX 1.6.6 created this file. For more info see http://www.lyx.org/
\lyxformat 345
\begin_document
\begin_header
@@ -2645,6 +2645,13 @@
\begin_layout Subsubsection
Power-Law Maxwell Viscoelastic Material
+\begin_inset CommandInset label
+LatexCommand label
+name "sub:Power-Law-Maxwell-Viscoelastic"
+
+\end_inset
+
+
\end_layout
\begin_layout Standard
@@ -3009,8 +3016,14 @@
.
The user then specifies either a reference stress or a reference strain
rate, and a database suitable for PyLith is generated.
- This utility is described more fully in the Appendix.
- ***Maybe describe in Tutorials instead?***
+ This utility is described more fully in
+\begin_inset CommandInset ref
+LatexCommand ref
+reference "sub:Tutorial-Step08-Power-law"
+
+\end_inset
+
+.
\end_layout
\begin_layout Standard
Modified: short/3D/PyLith/trunk/doc/userguide/tutorials/3dhex8/3dhex8.lyx
===================================================================
--- short/3D/PyLith/trunk/doc/userguide/tutorials/3dhex8/3dhex8.lyx 2010-06-08 23:34:40 UTC (rev 16944)
+++ short/3D/PyLith/trunk/doc/userguide/tutorials/3dhex8/3dhex8.lyx 2010-06-08 23:40:00 UTC (rev 16945)
@@ -1,4 +1,4 @@
-#LyX 1.6.4 created this file. For more info see http://www.lyx.org/
+#LyX 1.6.6 created this file. For more info see http://www.lyx.org/
\lyxformat 345
\begin_document
\begin_header
@@ -516,7 +516,7 @@
ElasticIsotropic3D
\family default
material type for both materials.
- This behavior will be overriden by example-specific
+ This behavior will be overriden by example-specific
\family typewriter
.cfg
\family default
Added: short/3D/PyLith/trunk/doc/userguide/tutorials/3dhex8/quasistatic/figs/step07-displ-vel-t300.jpg
===================================================================
(Binary files differ)
Property changes on: short/3D/PyLith/trunk/doc/userguide/tutorials/3dhex8/quasistatic/figs/step07-displ-vel-t300.jpg
___________________________________________________________________
Name: svn:mime-type
+ image/jpeg
Modified: short/3D/PyLith/trunk/doc/userguide/tutorials/3dhex8/quasistatic/quasistatic.lyx
===================================================================
--- short/3D/PyLith/trunk/doc/userguide/tutorials/3dhex8/quasistatic/quasistatic.lyx 2010-06-08 23:34:40 UTC (rev 16944)
+++ short/3D/PyLith/trunk/doc/userguide/tutorials/3dhex8/quasistatic/quasistatic.lyx 2010-06-08 23:40:00 UTC (rev 16945)
@@ -206,10 +206,13 @@
\end_layout
\begin_layout Standard
-The step04.cfg file describes a problem with x-displacements fixed at zero
- on the positive and negative x-faces while velocity boundary conditions
- are applied in the y-directions on the same faces, yielding a left-lateral
- sense of movement.
+The
+\family typewriter
+step04.cfg
+\family default
+ file describes a problem with x-displacements fixed at zero on the positive
+ and negative x-faces while velocity boundary conditions are applied in
+ the y-directions on the same faces, yielding a left-lateral sense of movement.
The bottom (minimum z) boundary is held fixed in the z-direction.
We also use a Maxwell viscoelastic material for the lower crust, and the
simulation is run for 200 years using a constant time step size of 20 years.
@@ -243,7 +246,7 @@
\end_layout
\begin_layout Standard
-We than change the material type of the lower crust, provide a spatial database
+We then change the material type of the lower crust, provide a spatial database
from which to obtain the material properties (using the default
\family typewriter
SimpleDB
@@ -300,10 +303,6 @@
output.cell_data_fields = [total_strain,stress,viscous_strain]
\end_layout
-\begin_layout LyX-Code
-
-\end_layout
-
\begin_layout Standard
Note that the default
\family typewriter
@@ -438,10 +437,6 @@
db_rate.data = [0.0*cm/year,-1.0*cm/year,0.0*year]
\end_layout
-\begin_layout LyX-Code
-
-\end_layout
-
\begin_layout Standard
Note that
\family typewriter
@@ -495,7 +490,7 @@
\end_layout
\begin_layout LyX-Code
-writer.filename = step04/step04.vtk
+writer.filename = output/step04.vtk
\end_layout
\begin_layout LyX-Code
@@ -539,7 +534,7 @@
\end_layout
\begin_layout LyX-Code
-writer.filename = step04/step04-groundsurf.vtk
+writer.filename = output/step04-groundsurf.vtk
\end_layout
\begin_layout LyX-Code
@@ -566,9 +561,13 @@
When we have run the simulation, the output VTK files will be contained
in
\family typewriter
-examples/3d/hex8/step04
+examples/3d/hex8/output
\family default
-, and the results may be visualized using a tool such as ParaView or mayavi2.
+ (all with a prefix of
+\family typewriter
+step04
+\family default
+), and the results may be visualized using a tool such as ParaView or mayavi2.
Results using ParaView are shown in
\begin_inset CommandInset ref
LatexCommand ref
@@ -782,7 +781,7 @@
\end_layout
\begin_layout LyX-Code
-db_rate.label = Dirichlet rate BC on +x
+db_rate.label = Dirichlet rate BC on -x
\end_layout
\begin_layout LyX-Code
@@ -799,9 +798,13 @@
When we have run the simulation, the output VTK files will be contained
in
\family typewriter
-examples/3d/hex8/step05
+examples/3d/hex8/output
\family default
-, and the results may be visualized using a tool such as ParaView or mayavi2.
+ (all with a prefix of
+\family typewriter
+step05
+\family default
+), and the results may be visualized using a tool such as ParaView or mayavi2.
Results using ParaView are shown in
\begin_inset CommandInset ref
LatexCommand ref
@@ -826,7 +829,7 @@
\begin_inset Caption
\begin_layout Plain Layout
-Displacement field for example step05 visualized using ParaView.
+Displacement field for example step05 at t = 40 years visualized using ParaView.
The mesh has been distorted by the computed displacements (magnified by
500), and the vectors show the computed displacements.
\begin_inset CommandInset label
@@ -882,6 +885,80 @@
\end_layout
\begin_layout Standard
+In this example and the remainder of the examples in this section, we change
+ the time stepping behavior from the default
+\family typewriter
+TimeStepUniform
+\family default
+ to
+\family typewriter
+TimeStepAdapt
+\family default
+.
+ For adaptive time stepping, we provide the maximum permissible time step
+ size, along with a stability factor.
+ The stability factor controls how the time step size is allowed to deviate
+ from the stable time step size provided by the different materials in the
+ model.
+ A
+\family typewriter
+stability_factor
+\family default
+ of 1.0 means we should use the stable time step size, while a
+\family typewriter
+stability_factor
+\family default
+ greater than 1.0 means we want to use a smaller time step size.
+ A
+\family typewriter
+stability_factor
+\family default
+ less than 1.0 allows time step sizes greater than the stable time step size,
+ which may provide inaccurate results.
+ The adaptive time stepping information is provided as:
+\end_layout
+
+\begin_layout LyX-Code
+# Change time stepping algorithm from uniform time step, to adaptive
+\end_layout
+
+\begin_layout LyX-Code
+# time stepping.
+\end_layout
+
+\begin_layout LyX-Code
+time_step = pylith.problems.TimeStepAdapt
+\end_layout
+
+\begin_layout LyX-Code
+
+\end_layout
+
+\begin_layout LyX-Code
+# Change the total simulation time to 700 years, and set the maximum time
+\end_layout
+
+\begin_layout LyX-Code
+# step size to 10 years.
+\end_layout
+
+\begin_layout LyX-Code
+[pylithapp.timedependent.implicit.time_step]
+\end_layout
+
+\begin_layout LyX-Code
+total_time = 700.0*year
+\end_layout
+
+\begin_layout LyX-Code
+max_dt = 10.0*year
+\end_layout
+
+\begin_layout LyX-Code
+stability_factor = 1.0 ; use time step equal to stable value from materials
+\end_layout
+
+\begin_layout Standard
Due to the simplicity of the boundary conditions, we are able to use the
default
\family typewriter
@@ -1039,7 +1116,23 @@
\end_layout
\begin_layout Standard
-Since there are multiple earthquake sources of different types, there are
+Note that for all earthquake sources we provide both an
+\family typewriter
+origin_time
+\family default
+ and a
+\family typewriter
+slip_function.slip_time
+\family default
+.
+ The first provides the starting time for the entire earthquake source,
+ while the second provides any spatial variation in the slip time with respect
+ to the
+\family typewriter
+origin_time
+\family default
+ (if any).
+ Since there are multiple earthquake sources of different types, there are
a number of additional fault information fields available for output.
We therefore request that these additional fields be output to the fault
information file:
@@ -1062,13 +1155,17 @@
\family default
.
It will contain final slip information for each earthquake source along
- with origin time information.
+ with slip time information.
When we have run the simulation, the output VTK files will be contained
in
\family typewriter
-examples/3d/hex8/step06
+examples/3d/hex8/output
\family default
-, and the results may be visualized using a tool such as ParaView or mayavi2.
+ (all with a prefix of
+\family typewriter
+step06
+\family default
+), and the results may be visualized using a tool such as ParaView or mayavi2.
Results using ParaView are shown in
\begin_inset CommandInset ref
LatexCommand ref
@@ -1093,8 +1190,8 @@
\begin_inset Caption
\begin_layout Plain Layout
-Displacement field for example step06 visualized using ParaView at t = 300
- years.
+Displacement field for example step06 at t = 300 years visualized using
+ ParaView.
The mesh has been distorted by the computed displacements (magnified by
500), and the vectors show the computed displacements.
\begin_inset CommandInset label
@@ -1116,5 +1213,639 @@
\end_layout
+\begin_layout Subsubsection
+Step07 - Dirichlet Velocity Boundary Conditions with Time-Dependent Kinematic
+ Fault Slip
+\end_layout
+
+\begin_layout Standard
+The
+\family typewriter
+step07.cfg
+\family default
+ file describes a problem nearly identical to example step06.
+ The only difference is that velocity boundary conditions in the positive
+ and negative y-directions on the positive and negative x-faces, so that
+ the external boundaries keep pace with the average fault slip.
+ This problem is nearly identical to the strike-slip fault model of
+\begin_inset CommandInset citation
+LatexCommand cite
+key "Savage:Prescott:1978"
+
+\end_inset
+
+, except that the fault creep extends through the viscoelastic portion of
+ the domain.
+\end_layout
+
+\begin_layout Standard
+We use the default
+\family typewriter
+ZeroDispBC
+\family default
+ for the initial displacements on the positive and negative x-faces, as
+ well as the negative z-face.
+ For the velocities on the positive and negative x-faces, we use a
+\family typewriter
+UniformDB
+\family default
+:
+\end_layout
+
+\begin_layout LyX-Code
+# +x face
+\end_layout
+
+\begin_layout LyX-Code
+[pylithapp.timedependent.bc.x_pos]
+\end_layout
+
+\begin_layout LyX-Code
+bc_dof = [0, 1]
+\end_layout
+
+\begin_layout LyX-Code
+label = face_xpos
+\end_layout
+
+\begin_layout LyX-Code
+db_initial.label = Dirichlet BC on +x
+\end_layout
+
+\begin_layout LyX-Code
+db_rate = spatialdata.spatialdb.UniformDB
+\end_layout
+
+\begin_layout LyX-Code
+db_rate.label = Dirichlet rate BC on +x
+\end_layout
+
+\begin_layout LyX-Code
+db_rate.values = [displacement-rate-x,displacement-rate-y,rate-start-time]
+\end_layout
+
+\begin_layout LyX-Code
+db_rate.data = [0.0*cm/year,1.0*cm/year,0.0*year]
+\end_layout
+
+\begin_layout LyX-Code
+
+\end_layout
+
+\begin_layout LyX-Code
+# -x face
+\end_layout
+
+\begin_layout LyX-Code
+[pylithapp.timedependent.bc.x_neg]
+\end_layout
+
+\begin_layout LyX-Code
+bc_dof = [0, 1]
+\end_layout
+
+\begin_layout LyX-Code
+label = face_xneg
+\end_layout
+
+\begin_layout LyX-Code
+db_initial.label = Dirichlet BC on -x
+\end_layout
+
+\begin_layout LyX-Code
+db_rate = spatialdata.spatialdb.UniformDB
+\end_layout
+
+\begin_layout LyX-Code
+db_rate.label = Dirichlet rate BC on +x
+\end_layout
+
+\begin_layout LyX-Code
+db_rate.values = [displacement-rate-x,displacement-rate-y,rate-start-time]
+\end_layout
+
+\begin_layout LyX-Code
+db_rate.data = [0.0*cm/year,-1.0*cm/year,0.0*year]
+\end_layout
+
+\begin_layout Standard
+The fault definition information is identical to example step06.
+ In previous examples, we have just used the default output for the domain
+ and subdomain (ground surface), which just provides the displacements.
+ In many cases, it is also useful to provide the velocities.
+ PyLith is able to provide this information, approximating the velocities
+ for the current time step as the difference between the current displacements
+ and the displacements from the previous time step, divided by the time
+ step size.
+ We can obtain this information by explicitly requesting it in
+\family typewriter
+vertex_data_fields
+\family default
+:
+\end_layout
+
+\begin_layout LyX-Code
+# Give basename for VTK domain output of solution over domain.
+\end_layout
+
+\begin_layout LyX-Code
+[pylithapp.problem.formulation.output.domain]
+\end_layout
+
+\begin_layout LyX-Code
+# We specify that output occurs in terms of a given time frequency, and
+\end_layout
+
+\begin_layout LyX-Code
+# ask for output every 50 years.
+ The time stamps of the output files are
+\end_layout
+
+\begin_layout LyX-Code
+# in years (rather than the default of seconds), and we give a format for
+\end_layout
+
+\begin_layout LyX-Code
+# the time stamp.
+\end_layout
+
+\begin_layout LyX-Code
+# We also request velocity output in addition to displacements.
+\end_layout
+
+\begin_layout LyX-Code
+vertex_data_fields = [displacement,velocity]
+\end_layout
+
+\begin_layout LyX-Code
+output_freq = time_step
+\end_layout
+
+\begin_layout LyX-Code
+time_step = 50.0*year
+\end_layout
+
+\begin_layout LyX-Code
+writer.filename = output/step07.vtk
+\end_layout
+
+\begin_layout LyX-Code
+writer.time_format = %04.0f
+\end_layout
+
+\begin_layout LyX-Code
+writer.time_constant = 1.0*year
+\end_layout
+
+\begin_layout LyX-Code
+
+\end_layout
+
+\begin_layout LyX-Code
+# Give basename for VTK domain output of solution over ground surface.
+\end_layout
+
+\begin_layout LyX-Code
+[pylithapp.problem.formulation.output.subdomain]
+\end_layout
+
+\begin_layout LyX-Code
+# Name of nodeset for ground surface.
+\end_layout
+
+\begin_layout LyX-Code
+label = face_zpos
+\end_layout
+
+\begin_layout LyX-Code
+# We also request velocity output in addition to displacements.
+\end_layout
+
+\begin_layout LyX-Code
+vertex_data_fields = [displacement,velocity]
+\end_layout
+
+\begin_layout LyX-Code
+# We keep the default output frequency behavior (skip every n steps), and
+\end_layout
+
+\begin_layout LyX-Code
+# ask to skip 0 steps between output, so that we get output every time step.
+\end_layout
+
+\begin_layout LyX-Code
+skip = 0
+\end_layout
+
+\begin_layout LyX-Code
+writer.filename = output/step07-groundsurf.vtk
+\end_layout
+
+\begin_layout LyX-Code
+writer.time_format = %04.0f
+\end_layout
+
+\begin_layout LyX-Code
+writer.time_constant = 1.0*year
+\end_layout
+
+\begin_layout Standard
+When we have run the simulation, the output VTK files will be contained
+ in
+\family typewriter
+examples/3d/hex8/output
+\family default
+ (all with a prefix of
+\family typewriter
+step07
+\family default
+), and the results may be visualized using a tool such as ParaView or mayavi2.
+ Results using ParaView are shown in
+\begin_inset CommandInset ref
+LatexCommand ref
+reference "fig:step07-displ-vel-t300"
+
+\end_inset
+
+.
+\begin_inset Float figure
+wide false
+sideways false
+status open
+
+\begin_layout Plain Layout
+\begin_inset Graphics
+ filename figs/step07-displ-vel-t300.jpg
+ width 10cm
+
+\end_inset
+
+
+\begin_inset Caption
+
+\begin_layout Plain Layout
+Displacement field (color contours) and velocity field (vectors) for example
+ step07 at t = 300 years visualized using ParaView.
+ The mesh has been distorted by the computed displacements (magnified by
+ 500), and the vectors show the computed velocities.
+\begin_inset CommandInset label
+LatexCommand label
+name "fig:step07-displ-vel-t300"
+
+\end_inset
+
+.
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Subsubsection
+Step08 - Dirichlet Velocity Boundary Conditions with Time-Dependent Kinematic
+ Fault Slip and Power-Law Rheology
+\begin_inset CommandInset label
+LatexCommand label
+name "sub:Tutorial-Step08-Power-law"
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+The
+\family typewriter
+step08.cfg
+\family default
+ file describes a problem that is identical to example step06, except the
+ the lower crust is composed of power-law viscoelastic material.
+ Since the material behavior is now nonlinear, we must use the nonlinear
+ solver:
+\end_layout
+
+\begin_layout LyX-Code
+[pylithapp.timedependent]
+\end_layout
+
+\begin_layout LyX-Code
+# For this problem we must switch to a nonlinear solver.
+\end_layout
+
+\begin_layout LyX-Code
+implicit.solver = pylith.problems.SolverNonlinear
+\end_layout
+
+\begin_layout Standard
+Although we have not discussed the PyLith PETSc settings previously, note
+ that the use of the nonlinear solver may require additional options if
+ we wish to override the defaults.
+ These settings are contained in
+\family typewriter
+pylithapp.cfg
+\family default
+:
+\end_layout
+
+\begin_layout LyX-Code
+[pylithapp.petsc]
+\end_layout
+
+\begin_layout LyX-Code
+# Nonlinear solver monitoring options.
+\end_layout
+
+\begin_layout LyX-Code
+snes_rtol = 1.0e-8
+\end_layout
+
+\begin_layout LyX-Code
+snes_atol = 1.0e-12
+\end_layout
+
+\begin_layout LyX-Code
+snes_max_it = 100
+\end_layout
+
+\begin_layout LyX-Code
+snes_monitor = true
+\end_layout
+
+\begin_layout LyX-Code
+snes_view = true
+\end_layout
+
+\begin_layout LyX-Code
+snes_converged_reason = true
+\end_layout
+
+\begin_layout Standard
+These settings are unused unless we are using the nonlinear solver.
+\end_layout
+
+\begin_layout Standard
+When using laboratory results for power-law viscoelastic materials, the
+ properties provided do not generally match those expected by PyLith (see
+
+\begin_inset CommandInset ref
+LatexCommand ref
+reference "sub:Power-Law-Maxwell-Viscoelastic"
+
+\end_inset
+
+).
+ A utility code,
+\family typewriter
+powerlaw_gendb.py
+\family default
+, has been provided to simplify the process of using laboratory results
+ with PyLith.
+ This utility code should be installed in the same location as PyLith.
+ An example of how to use it is in
+\family typewriter
+examples/3d/hex8/spatialdb/powerlaw
+\family default
+.
+ The user must provide a spatial database defining the spatial distribution
+ of laboratory-derived parameters (contained in
+\family typewriter
+powerlaw_params.spatialdb
+\family default
+), another spatial database defining the temperature field in degrees K
+ (contained in
+\family typewriter
+temperature.spatialdb
+\family default
+), and a set of points for which values are desired (
+\family typewriter
+powerlaw_points.txt
+\family default
+).
+ The parameters for the code are defined in
+\family typewriter
+powerlaw_gendb.cfg
+\family default
+.
+ The properties expected by PyLith are
+\family typewriter
+reference_strain_rate
+\family default
+,
+\family typewriter
+reference_stress
+\family default
+, and
+\family typewriter
+power_law_exponent
+\family default
+.
+ The user must specify either
+\family typewriter
+reference_strain_rate
+\family default
+ or
+\family typewriter
+reference_stress
+\family default
+ so that
+\family typewriter
+powerlaw_gendb.py
+\family default
+ can compute the other property.
+ Default values of 1.0e-6 1/s and 1 MPa are provided.
+ A spatial database must then be provided for each laboratory-derived parameter
+ and for the temperature.
+ In this example, the same database was used for all parameters, and a separate
+ database was used to define the temperature distribution.
+ In practice, the user can provide any desired thermal model to provide
+ the spatial database for the temperature.
+ In this example, a simple 1D (vertically-varying) distribution was used.
+\end_layout
+
+\begin_layout Standard
+The fault definition information is identical to example step06.
+ In previous examples, we have just used the default output for the domain
+ and subdomain (ground surface), which just provides the displacements.
+ In many cases, it is also useful to provide the velocities.
+ PyLith is able to provide this information, approximating the velocities
+ for the current time step as the difference between the current displacements
+ and the displacements from the previous time step, divided by the time
+ step size.
+ We can obtain this information by explicitly requesting it in
+\family typewriter
+vertex_data_fields
+\family default
+:
+\end_layout
+
+\begin_layout LyX-Code
+# Give basename for VTK domain output of solution over domain.
+\end_layout
+
+\begin_layout LyX-Code
+[pylithapp.problem.formulation.output.domain]
+\end_layout
+
+\begin_layout LyX-Code
+# We specify that output occurs in terms of a given time frequency, and
+\end_layout
+
+\begin_layout LyX-Code
+# ask for output every 50 years.
+ The time stamps of the output files are
+\end_layout
+
+\begin_layout LyX-Code
+# in years (rather than the default of seconds), and we give a format for
+\end_layout
+
+\begin_layout LyX-Code
+# the time stamp.
+\end_layout
+
+\begin_layout LyX-Code
+# We also request velocity output in addition to displacements.
+\end_layout
+
+\begin_layout LyX-Code
+vertex_data_fields = [displacement,velocity]
+\end_layout
+
+\begin_layout LyX-Code
+output_freq = time_step
+\end_layout
+
+\begin_layout LyX-Code
+time_step = 50.0*year
+\end_layout
+
+\begin_layout LyX-Code
+writer.filename = output/step07.vtk
+\end_layout
+
+\begin_layout LyX-Code
+writer.time_format = %04.0f
+\end_layout
+
+\begin_layout LyX-Code
+writer.time_constant = 1.0*year
+\end_layout
+
+\begin_layout LyX-Code
+
+\end_layout
+
+\begin_layout LyX-Code
+# Give basename for VTK domain output of solution over ground surface.
+\end_layout
+
+\begin_layout LyX-Code
+[pylithapp.problem.formulation.output.subdomain]
+\end_layout
+
+\begin_layout LyX-Code
+# Name of nodeset for ground surface.
+\end_layout
+
+\begin_layout LyX-Code
+label = face_zpos
+\end_layout
+
+\begin_layout LyX-Code
+# We also request velocity output in addition to displacements.
+\end_layout
+
+\begin_layout LyX-Code
+vertex_data_fields = [displacement,velocity]
+\end_layout
+
+\begin_layout LyX-Code
+# We keep the default output frequency behavior (skip every n steps), and
+\end_layout
+
+\begin_layout LyX-Code
+# ask to skip 0 steps between output, so that we get output every time step.
+\end_layout
+
+\begin_layout LyX-Code
+skip = 0
+\end_layout
+
+\begin_layout LyX-Code
+writer.filename = output/step07-groundsurf.vtk
+\end_layout
+
+\begin_layout LyX-Code
+writer.time_format = %04.0f
+\end_layout
+
+\begin_layout LyX-Code
+writer.time_constant = 1.0*year
+\end_layout
+
+\begin_layout Standard
+When we have run the simulation, the output VTK files will be contained
+ in
+\family typewriter
+examples/3d/hex8/output
+\family default
+ (all with a prefix of
+\family typewriter
+step07
+\family default
+), and the results may be visualized using a tool such as ParaView or mayavi2.
+ Results using ParaView are shown in
+\begin_inset CommandInset ref
+LatexCommand ref
+reference "fig:step07-displ-vel-t300-1"
+
+\end_inset
+
+.
+\begin_inset Float figure
+wide false
+sideways false
+status open
+
+\begin_layout Plain Layout
+\begin_inset Graphics
+ filename /Users/charlesw/geoframe/cig/short/3D/PyLith/trunk/doc/userguide/tutorials/3dhex8/quasistatic/figs/step07-displ-vel-t300.jpg
+ width 10cm
+
+\end_inset
+
+
+\begin_inset Caption
+
+\begin_layout Plain Layout
+Displacement field (color contours) and velocity field (vectors) for example
+ step07 at t = 300 years visualized using ParaView.
+ The mesh has been distorted by the computed displacements (magnified by
+ 500), and the vectors show the computed velocities.
+\begin_inset CommandInset label
+LatexCommand label
+name "fig:step07-displ-vel-t300-1"
+
+\end_inset
+
+.
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
\end_body
\end_document
Modified: short/3D/PyLith/trunk/doc/userguide/tutorials/3dhex8/static/static.lyx
===================================================================
--- short/3D/PyLith/trunk/doc/userguide/tutorials/3dhex8/static/static.lyx 2010-06-08 23:34:40 UTC (rev 16944)
+++ short/3D/PyLith/trunk/doc/userguide/tutorials/3dhex8/static/static.lyx 2010-06-08 23:40:00 UTC (rev 16945)
@@ -165,9 +165,12 @@
\end_layout
\begin_layout Standard
-The step01.cfg file describes a problem with pure Dirichlet (displacement)
- boundary conditions corresponding to compression in the x-direction and
- shear in the y-direction.
+The
+\family typewriter
+step01.cfg
+\family default
+ file describes a problem with pure Dirichlet (displacement) boundary conditions
+ corresponding to compression in the x-direction and shear in the y-direction.
The bottom (minimum z) boundary is held fixed in the z-direction.
On the positive and negative x-faces, compressional displacements of 1
m are applied in the x-direction and shear displacements yielding a left-latera
@@ -229,7 +232,7 @@
\end_layout
\begin_layout LyX-Code
-writer.filename = step01/step01-groundsurf.vtk
+writer.filename = output/step01-groundsurf.vtk
\end_layout
\begin_layout Standard
@@ -238,7 +241,7 @@
\family typewriter
bc_dof
\family default
-), we must provide a label that describes which Cubit nodeset to which to
+), we must provide a label that describes the Cubit nodeset to which to
apply the BC, and we must specify what type of spatial database is being
used to describe the boundary conditions.
For the x-faces, we use a
@@ -351,13 +354,17 @@
When we have run the simulation, the output VTK files will be contained
in
\family typewriter
-examples/3d/hex8/step01
+examples/3d/hex8/output
\family default
-, and the results may be visualized using a tool such as ParaView or mayavi2.
+ (all with a prefix of
+\family typewriter
+step01
+\family default
+), and the results may be visualized using a tool such as ParaView or mayavi2.
Results using ParaView are shown in
\begin_inset CommandInset ref
LatexCommand ref
-reference "fig:step04-displ-t200"
+reference "fig:step01-displ"
\end_inset
@@ -415,10 +422,14 @@
\end_layout
\begin_layout Standard
-The step02.cfg file describes a problem with Dirichlet (displacement) boundary
- conditions corresponding to zero x and y-displacements applied on the negative
- x-face and Neumann (traction) boundary conditions corresponding to normal
- compression and horizontal shear applied on the positive x-face.
+The
+\family typewriter
+step02.cfg
+\family default
+ file describes a problem with Dirichlet (displacement) boundary conditions
+ corresponding to zero x and y-displacements applied on the negative x-face
+ and Neumann (traction) boundary conditions corresponding to normal compression
+ and horizontal shear applied on the positive x-face.
The bottom (minimum z) boundary is held fixed in the z-direction.
The problem is similar to example step01, except that 1 MPa of normal compressi
on and 1 MPa of shear (in a left-lateral sense) are applied on the positive
@@ -542,13 +553,17 @@
When we have run the simulation, the output VTK files will be contained
in
\family typewriter
-examples/3d/hex8/step02
+examples/3d/hex8/output
\family default
-, and the results may be visualized using a tool such as ParaView or mayavi2.
+ (all with a prefix of
+\family typewriter
+step02
+\family default
+), and the results may be visualized using a tool such as ParaView or mayavi2.
Results using ParaView are shown in
\begin_inset CommandInset ref
LatexCommand ref
-reference "fig:step05-displ-t40"
+reference "fig:step02-displ"
\end_inset
@@ -600,9 +615,13 @@
\end_layout
\begin_layout Standard
-The step03.cfg file describes a problem with Dirichlet (displacement) boundary
- conditions corresponding to zero x and y-displacements applied on the negative
- and positive x-faces and a vertical fault with a combination of left-lateral
+The
+\family typewriter
+step03.cfg
+\family default
+ file describes a problem with Dirichlet (displacement) boundary conditions
+ corresponding to zero x and y-displacements applied on the negative and
+ positive x-faces and a vertical fault with a combination of left-lateral
and updip motion.
The left-lateral component of fault slip has a constant value of 2 m in
the upper crust, and then decreases linearly to zero at the base of the
@@ -710,8 +729,8 @@
\family typewriter
SimpleDB
\family default
- to define the time at which slip begins (the start time is 0.0 since this
- is a static problem):
+ to define the spatial variation in slip initiation times (the start time
+ is 0.0 everywhere since this is a static problem):
\end_layout
\begin_layout LyX-Code
@@ -741,7 +760,7 @@
\end_layout
\begin_layout LyX-Code
-writer.filename = step03/step03-fault.vtk
+writer.filename = output/step03-fault.vtk
\end_layout
\begin_layout Standard
@@ -761,13 +780,17 @@
When we have run the simulation, the output VTK files will be contained
in
\family typewriter
-examples/3d/hex8/step03
+examples/3d/hex8/output
\family default
-, and the results may be visualized using a tool such as ParaView or mayavi2.
+ (all with a prefix of
+\family typewriter
+step03
+\family default
+), and the results may be visualized using a tool such as ParaView or mayavi2.
Results using ParaView are shown in
\begin_inset CommandInset ref
LatexCommand ref
-reference "fig:step06-displ-t300"
+reference "fig:step03-displ"
\end_inset
More information about the CIG-COMMITS
mailing list