[cig-commits] r13660 - doc/snac
echoi at geodynamics.org
echoi at geodynamics.org
Thu Dec 11 13:59:11 PST 2008
Author: echoi
Date: 2008-12-11 13:59:10 -0800 (Thu, 11 Dec 2008)
New Revision: 13660
Modified:
doc/snac/snac.lyx
Log:
1. Chapters on Benchmark problems have been combined.
2. An example input file added.
Modified: doc/snac/snac.lyx
===================================================================
--- doc/snac/snac.lyx 2008-12-11 21:18:17 UTC (rev 13659)
+++ doc/snac/snac.lyx 2008-12-11 21:59:10 UTC (rev 13660)
@@ -2889,7 +2889,7 @@
.
SNAC will not be confused even if both sets of parameters are present in
one input file.
- Finally,
+
\family typewriter
meshType
\family default
@@ -2907,7 +2907,44 @@
\family default
by default.
However, it is important to assign a correct value for remeshing.
-
+ Finally, two parameters for a parallel run: ``
+\family typewriter
+shadowDepth
+\family default
+
+\begin_inset Quotes srd
+\end_inset
+
+ and ``
+\family typewriter
+decompDims
+\family default
+
+\begin_inset Quotes srd
+\end_inset
+
+.
+ ``
+\family typewriter
+shadowDepth
+\family default
+
+\begin_inset Quotes srd
+\end_inset
+
+ determines the depth from the parallel boundary in terms of the number
+ of elements.
+ To run Snac in parallel, this parameter should be 1 or larger.
+ ``
+\family typewriter
+decompDims
+\family default
+
+\begin_inset Quotes srd
+\end_inset
+
+ determines whether parallel domain decomposition is performed in 1, 2 or
+ 3 dimensions.
\end_layout
\begin_layout Subsection
@@ -4040,9 +4077,947 @@
\end_layout
\begin_layout LyX-Code
+<?xml version="1.0"?>
+\end_layout
+\begin_layout LyX-Code
+<!DOCTYPE StGermainData SYSTEM "stgermain.dtd">
\end_layout
+\begin_layout LyX-Code
+<!-- StGermain-Snac input file -->
+\end_layout
+
+\begin_layout LyX-Code
+<StGermainData xmlns="http://www.vpac.org/StGermain/XML_IO_Handler/Jun2003">
+
+\end_layout
+
+\begin_layout LyX-Code
+<!-- StGermain simulation parameters -->
+\end_layout
+
+\begin_layout LyX-Code
+<param name="start"> 0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="outputPath">./test6-2</param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="dumpEvery"> 1000 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="maxTimeSteps"> 100001 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<!-- Snac variables -->
+\end_layout
+
+\begin_layout LyX-Code
+<param name="density"> 2700 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="gravity"> 9.8 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="demf"> 0.8 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="dtType"> constant </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="timeStep"> 6.3072e+07 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="forceCalcType"> complete </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="decomposedAxis"> 0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="storeForces"> no </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="forceCheckSum"> no </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="topo_kappa"> 0.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="alpha"> 0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<!-- Extension modules -->
+\end_layout
+
+\begin_layout LyX-Code
+<list name="extensions">
+\end_layout
+
+\begin_layout LyX-Code
+ <param> SnacTemperature </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param> SnacViscoPlastic </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param> SnacHydroStaticIC </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param> SnacVPSeeds </param>
+\end_layout
+
+\begin_layout LyX-Code
+</list>
+\end_layout
+
+\begin_layout LyX-Code
+<struct name="mesh">
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="shadowDepth"> 1 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="decompDims"> 2 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <!-- Mesh size -->
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="meshSizeI"> 51 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="meshSizeJ"> 11 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="meshSizeK"> 53 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <!-- Initial geometry -->
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="minX"> 0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="minY"> 0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="minZ"> 0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="maxX"> 100000 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="maxY"> 20000 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="maxZ"> 104000 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <!-- Remeshing -->
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="meshType"> cartesian </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="buildNodeNeighbourTbl"> True </param>
+\end_layout
+
+\begin_layout LyX-Code
+</struct>
+\end_layout
+
+\begin_layout LyX-Code
+<!-- Elastic material parameters -->
+\end_layout
+
+\begin_layout LyX-Code
+<param name="lambda"> 1.0e+10 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="mu"> 1.0e+10 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<!-- Viscous material parameters -->
+\end_layout
+
+\begin_layout LyX-Code
+<param name="refvisc"> 1.0e+20 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="reftemp"> 1400.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="activationE"> 45.0e+03 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="vis_min"> 1.0e+30 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="vis_max"> 1.0e+30 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="srexponent"> 1 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<!-- Plastic material parameters -->
+\end_layout
+
+\begin_layout LyX-Code
+<param name="yieldcriterion"> mohrcoulomb </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="nsegments"> 2 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="plstrain0"> 0.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="plstrain1"> 0.02 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="plstrain2"> 1000.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="frictionAngle0"> 30.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="frictionAngle1"> 30.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="frictionAngle2"> 30.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="dilationAngle0"> 5.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="dilationAngle1"> 5.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="dilationAngle2"> 5.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="cohesion0"> 4.0e+07 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="cohesion1"> 4.0e+05 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="cohesion2"> 0.0e+00 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="ten_off"> 1.0e+13 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<!-- Temperature variables -->
+\end_layout
+
+\begin_layout LyX-Code
+<param name="topTemp"> 0.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="bottomTemp"> 700.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="thermal_conduct"> 1.6 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="heatCapacity"> 1000.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<!-- Remesher info -->
+\end_layout
+
+\begin_layout LyX-Code
+<!-- Used only when the SnacRemesher plugin is loaded.
+ -->
+\end_layout
+
+\begin_layout LyX-Code
+<!-- <param name="remeshCondition"> onBothTimeStepLength </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="remeshCondition"> onTimeStep </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="remeshCondition"> onMinLengthScale </param> -->
+\end_layout
+
+\begin_layout LyX-Code
+<param name="remeshCondition"> onMinLengthScale </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="remeshTimeStepCriterion"> 15000 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="remeshLengthCriterion"> 0.7 </param>
+\end_layout
+
+\begin_layout LyX-Code
+<param name="bottomResotre"> on </param>
+\end_layout
+
+\begin_layout LyX-Code
+<!-- node ICs -->
+\end_layout
+
+\begin_layout LyX-Code
+<struct name="nodeICs">
+\end_layout
+
+\begin_layout LyX-Code
+ <list name="vcList">
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> AllNodesVC </param>
+
+\end_layout
+
+\begin_layout LyX-Code
+ <list name="variables">
+
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="name"> vx </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> double </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="value"> 0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="name"> vy </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> double </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="value"> 0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="name"> vz </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> double </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="value"> 0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ </list>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> AllNodesVC </param>
+
+\end_layout
+
+\begin_layout LyX-Code
+ <list name="variables">
+
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="name"> temperature </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> func </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="value"> SnacTemperature_Top2BottomSweep
+ </param> </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ </list>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ </list>
+\end_layout
+
+\begin_layout LyX-Code
+</struct>
+\end_layout
+
+\begin_layout LyX-Code
+<!-- element ICs -->
+\end_layout
+
+\begin_layout LyX-Code
+<struct name="elementICs">
+\end_layout
+
+\begin_layout LyX-Code
+ <list name="vcList">
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> AllElementsVC </param>
+
+\end_layout
+
+\begin_layout LyX-Code
+ <list name="variables">
+
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="name"> elementMaterial </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> int </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="value"> 0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ </list>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ </list>
+\end_layout
+
+\begin_layout LyX-Code
+</struct>
+\end_layout
+
+\begin_layout LyX-Code
+<!-- Velocity BCs -->
+\end_layout
+
+\begin_layout LyX-Code
+<struct name="velocityBCs">
+\end_layout
+
+\begin_layout LyX-Code
+ <list name="vcList">
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> WallVC </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="wall"> left </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <list name="variables">
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="name"> vx </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> double </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="value"> -3.17e-10 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ </list>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> WallVC </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="wall"> right </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <list name="variables">
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="name"> vx </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> double </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="value"> 3.17e-10 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ </list>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> WallVC </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="wall"> back </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <list name="variables">
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="name"> vz </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> double </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="value"> 0.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ </list>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> WallVC </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="wall"> front </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <list name="variables">
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="name"> vz </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> double </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="value"> 0.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ </list>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> WallVC </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="wall"> bottom </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <list name="variables">
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="name"> vx </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> func </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="value"> variableVelBC </param>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="name"> vy </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> double </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="value"> 0.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ </list>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ </list>
+\end_layout
+
+\begin_layout LyX-Code
+</struct>
+\end_layout
+
+\begin_layout LyX-Code
+<!-- Temperature BCs -->
+\end_layout
+
+\begin_layout LyX-Code
+<struct name="temperatureBCs">
+\end_layout
+
+\begin_layout LyX-Code
+ <list name="vcList">
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> WallVC </param>
+
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="wall"> top </param>
+
+\end_layout
+
+\begin_layout LyX-Code
+ <list name="variables">
+
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="name"> temperature </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> double </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="value"> 0.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ </list>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> WallVC </param>
+
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="wall"> bottom </param>
+
+\end_layout
+
+\begin_layout LyX-Code
+ <list name="variables">
+
+\end_layout
+
+\begin_layout LyX-Code
+ <struct>
+
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="name"> temperature </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="type"> double </param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="value"> 700.0 </param>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ </list>
+\end_layout
+
+\begin_layout LyX-Code
+ </struct>
+\end_layout
+
+\begin_layout LyX-Code
+ </list>
+\end_layout
+
+\begin_layout LyX-Code
+</struct>
+\end_layout
+
+\begin_layout LyX-Code
+</StGermainData>
+\end_layout
+
\begin_layout Chapter
Postprocessing and Graphics
\end_layout
@@ -4365,19 +5340,11 @@
\end_layout
\begin_layout Chapter
-SNAC's Plasticity Model (from TWiki)
+Benchmark problems
\end_layout
-\begin_layout Standard
-The yield criterion corresponds to a Mohr-Coulomb failure envelope with
- tension cut-off.
- A non-associated flow rule [is ? -- TODO] for shear failure, and an associated
- flow rule for tensional failure controls plastic flow.
-
-\end_layout
-
\begin_layout Section
-Benchmark problem I
+Oedometer Test
\end_layout
\begin_layout Standard
@@ -4557,7 +5524,7 @@
\end_layout
\begin_layout Section
-Benchmark Problem II
+Thick-cylinder with pressure on the inner wall I
\end_layout
\begin_layout Standard
@@ -4789,7 +5756,7 @@
\end_layout
\begin_layout Section
-Benchmark Problem III
+Thick-cylinder with pressure on the inner wall II
\end_layout
\begin_layout Standard
@@ -5224,12 +6191,12 @@
\end_layout
-\begin_layout Chapter
-ViscoElasticRheology
+\begin_layout Section
+Parallel-plate viscometer problem
\end_layout
\begin_layout Standard
-We solve a parallel-plate viscometer problem (
+A parallel-plate viscometer problem (
\begin_inset LatexCommand cite
key "Jaeger_1969"
@@ -5249,8 +6216,8 @@
New York: John Wiley & Sons, Inc.)
\end_layout
-\begin_layout Section
-Benchmark problem I
+\begin_layout Subsection
+Model Setup
\end_layout
\begin_layout Itemize
More information about the CIG-COMMITS
mailing list