[cig-commits] r13686 - in doc/snac: . figures
echoi at geodynamics.org
echoi at geodynamics.org
Sun Dec 14 17:26:44 PST 2008
Author: echoi
Date: 2008-12-14 17:26:44 -0800 (Sun, 14 Dec 2008)
New Revision: 13686
Added:
doc/snac/figures/example_plots.png
Modified:
doc/snac/snac.lyx
Log:
The Cookbook chapter has been incorporated into the "Running SNAC" chapter.
Added: doc/snac/figures/example_plots.png
===================================================================
(Binary files differ)
Property changes on: doc/snac/figures/example_plots.png
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Modified: doc/snac/snac.lyx
===================================================================
--- doc/snac/snac.lyx 2008-12-14 20:25:33 UTC (rev 13685)
+++ doc/snac/snac.lyx 2008-12-15 01:26:44 UTC (rev 13686)
@@ -1,4 +1,4 @@
-#LyX 1.5.1 created this file. For more info see http://www.lyx.org/
+#LyX 1.5.3 created this file. For more info see http://www.lyx.org/
\lyxformat 276
\begin_document
\begin_header
@@ -329,7 +329,7 @@
\begin_layout Standard
SNAC (StGermaiN Analysis of Continua) is an updated Lagrangian explicit
finite difference code for modeling a finitely deforming elasto-visco-plastic
- solid in 3D, released under the GNU General Public License (see Appendix
+ solid in 3D, released under the GNU General Public License (see Ch.
\begin_inset LatexCommand vref
reference "cha:License"
@@ -350,7 +350,13 @@
the standard FEM and subdivided into two layouts of five tetrahedral elements
for symmetric response.
To avoid the over-stiff response of tetrahedrons in the incompressible
- limit, mixed discretization is applied [Marti and Cundall, 1982].
+ limit, mixed discretization is applied [
+\begin_inset LatexCommand cite
+key "MartCund1982"
+
+\end_inset
+
+].
The mixed discretization relieves over-stiffness by replacing the first
invariant of tetrahedral strain-rate tensor with the one averaged over
a zone.
@@ -1042,12 +1048,13 @@
\begin_layout Standard
SNAC uses a general elasto-visco-plastic rheological model to update the
- Cauchy stress tensor, e.g., Albe, et al.
+ Cauchy stress tensor (e.g.,
\begin_inset LatexCommand cite
key "Albe_etal2000"
\end_inset
+).
First, the initial guess of stress is acquired by the Maxwell viscoelastic
constitutive law
\begin_inset LatexCommand cite
@@ -1232,8 +1239,13 @@
\begin_inset Formula $\sigma^{n+1}$
\end_inset
- onto the yield surface using a return-mapping algorithm [Simo and Hughes,
- 1999].
+ onto the yield surface using a return-mapping algorithm [
+\begin_inset LatexCommand cite
+key "SimoHugh2004"
+
+\end_inset
+
+].
\end_layout
@@ -2104,8 +2116,8 @@
\end_layout
\begin_layout Standard
-It is handy to define the following environment variables for both building
- and running SNAC.
+It is handy, although not necessary, to define the following environment
+ variables for both building and running SNAC:
\end_layout
\begin_layout Standard
@@ -2747,6 +2759,26 @@
: If ``yes,'' the sum of forces is checked if it is consistent.
\end_layout
+\begin_layout Description
+
+\family typewriter
+topo_kappa
+\family default
+: Parameter for topography smoothing.
+\end_layout
+
+\begin_layout Description
+
+\family typewriter
+alpha
+\family default
+: Volumetric thermal expansion coefficient in the unit of
+\begin_inset Formula $K^{-1}$
+\end_inset
+
+.
+\end_layout
+
\begin_layout Subsection
Plug-ins List
\end_layout
@@ -2809,7 +2841,6 @@
\family default
: for elastic, elasto-plastic, Maxwell viscoelastic, and elasto-visco-plastic
rheology, respectively.
- Note that ViscoPlastic and Maxwell requires to load SNAC Temperature.
\end_layout
@@ -3052,7 +3083,11 @@
\begin_inset Formula $T_{0}$
\end_inset
- is ``reftemp'' (
+ is ``
+\family typewriter
+reftemp
+\family default
+'' (
\begin_inset Formula $^{\circ}C$
\end_inset
@@ -3287,7 +3322,7 @@
\begin_deeper
\begin_layout Standard
-Only two parameters need to be defined: thermal conductivity and heat capacity.
+Values of thermal conductivity and heat capacity can be defined.
The parameter names are
\family typewriter
thermal_conduct
@@ -3297,6 +3332,28 @@
heatCapacity
\family default
, and appropriate values in MKS unit should be assigned.
+ The default values are 2.0 W/(m
+\begin_inset Formula $\cdot$
+\end_inset
+
+K) and 1000.0 J/(kg
+\begin_inset Formula $\cdot$
+\end_inset
+
+K), respectively.
+ To be used for initial and boundary conditions, temperature values in
+\begin_inset Formula $^{\circ}C$
+\end_inset
+
+ can be assigned to
+\family typewriter
+topTemp
+\family default
+ and
+\family typewriter
+bottomTemp
+\family default
+.
\end_layout
\end_deeper
@@ -4084,939 +4141,1540 @@
Complete Example XML file
\end_layout
+\begin_layout Subsection
+Problem
+\end_layout
+
+\begin_layout Standard
+A 3-D Cartesian block of elaso-visco-plastic material is extended by the
+ velocities applied to two side walls.
+ The domain has a size of 40
+\begin_inset Formula $\times$
+\end_inset
+
+10
+\begin_inset Formula $\times$
+\end_inset
+
+80 km and is discretized into 20
+\begin_inset Formula $\times$
+\end_inset
+
+5
+\begin_inset Formula $\times$
+\end_inset
+
+40 elements (Fig.
+
+\begin_inset LatexCommand ref
+reference "fig:example-problem"
+
+\end_inset
+
+a).
+\end_layout
+
+\begin_layout Standard
+The two side walls perpendicular to x-axis are pulled at a constant velocity
+ of 1 cm/yr.
+ The other two side walls have free-slip conditions.
+ The bottom surface is supported by the Winkler resotring forces.
+\end_layout
+
+\begin_layout Standard
+An initial temperature field is generated by a condition function that makes
+ a horizontally uniform temperature field with a constant vertical gradient.
+ The boundary conditions fix the top and bottom temperature at their initial
+ values (0 and 700
+\begin_inset Formula $^{\circ}C$
+\end_inset
+
+, respectively) to keep the steady state.
+\end_layout
+
+\begin_layout Standard
+The Mohr-Coulomb model is used with a two-step strain weakening.
+ Cohesion decreases from 40 MPa to 0.4 Mpa as the second invariant of plastic
+ strain grows up to 2 % but remains constant at 0.4 MPa thereafter.
+ Note that extremely large values of vis_min/max are chosen such that viscosity
+ is uniformly high and thus no viscous flow is allowed.
+\end_layout
+
+\begin_layout Standard
+A finite value of plastic strain (2 % in this problem) is assigned to a
+ set of elements.
+ These element have a lower cohesion according to the given strain weakening
+ rule and thus strain localization is initiated from them.
+ The selection of elements is randomly made within the region at a certain
+ distance from the left and right boundaries (Fig.
+
+\begin_inset LatexCommand ref
+reference "fig:example-problem"
+
+\end_inset
+
+a).
+\end_layout
+
+\begin_layout Subsection
+Results
+\end_layout
+
+\begin_layout Standard
+Continued extension builds up strain and the seed elements start to yield
+ first.
+ Initiated at the seed elements, localization of plastic strain propagates
+ in the direction normal to the applied velocities.
+ Eventually, a long basin bounded by conjugate normal faults is created
+ (Fig.
+
+\begin_inset LatexCommand ref
+reference "fig:example-problem"
+
+\end_inset
+
+b).
+\end_layout
+
+\begin_layout Standard
+\begin_inset Float figure
+wide false
+sideways false
+status open
+
+\begin_layout Standard
+\begin_inset Graphics
+ filename figures/example_plots.png
+ scale 90
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+\begin_inset Caption
+
+\begin_layout Standard
+\begin_inset LatexCommand label
+name "fig:example-problem"
+
+\end_inset
+
+(a) The mesh for the example problem with the locations of the ``seed
+\begin_inset Quotes srd
+\end_inset
+
+ elements.
+ (b) Plastic strain mapped on the deformed mesh after 2 km extension (0.1
+ My elapsed).
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Subsection
+Complete listing of the input XML file
+\end_layout
+
+\begin_layout Standard
+The problem described above can be reproduced using the input file,
+\begin_inset Quotes sld
+\end_inset
+
+Snac/examples/Cookbook1/cookbook1.xml
+\begin_inset Quotes srd
+\end_inset
+
+.
+ The full contents of this input file is listed below for reference.
+\end_layout
+
\begin_layout LyX-Code
+
+\size small
<?xml version="1.0"?>
\end_layout
\begin_layout LyX-Code
+
+\size small
<!DOCTYPE StGermainData SYSTEM "stgermain.dtd">
\end_layout
\begin_layout LyX-Code
+
+\size small
<!-- StGermain-Snac input file -->
\end_layout
\begin_layout LyX-Code
+
+\size small
<StGermainData xmlns="http://www.vpac.org/StGermain/XML_IO_Handler/Jun2003">
\end_layout
\begin_layout LyX-Code
+
+\size small
<!-- StGermain simulation parameters -->
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="start"> 0 </param>
\end_layout
\begin_layout LyX-Code
-<param name="outputPath">./test6-2</param>
+
+\size small
+<param name="outputPath">./data</param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="dumpEvery"> 1000 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="maxTimeSteps"> 100001 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<!-- Snac variables -->
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="density"> 2700 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="gravity"> 9.8 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="demf"> 0.8 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="dtType"> constant </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="timeStep"> 6.3072e+07 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="forceCalcType"> complete </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="decomposedAxis"> 0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="storeForces"> no </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="forceCheckSum"> no </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="topo_kappa"> 0.0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="alpha"> 0 </param>
\end_layout
\begin_layout LyX-Code
+
+\end_layout
+
+\begin_layout LyX-Code
+
+\size small
<!-- Extension modules -->
\end_layout
\begin_layout LyX-Code
+
+\size small
<list name="extensions">
\end_layout
\begin_layout LyX-Code
- <param> SnacTemperature </param>
+
+\size small
+ <param> SnacTemperature </param>
\end_layout
\begin_layout LyX-Code
- <param> SnacViscoPlastic </param>
+
+\size small
+ <param> SnacViscoPlastic </param>
\end_layout
\begin_layout LyX-Code
- <param> SnacHydroStaticIC </param>
+
+\size small
+ <param> SnacHydroStaticIC </param>
\end_layout
\begin_layout LyX-Code
- <param> SnacVPSeeds </param>
+
+\size small
+ <param> SnacVPSeeds </param>
\end_layout
\begin_layout LyX-Code
-</list>
+
+\size small
+ <param> SnacWinklerForce </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
+</list>
+\end_layout
+
+\begin_layout LyX-Code
+
+\size small
<struct name="mesh">
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="shadowDepth"> 1 </param>
\end_layout
\begin_layout LyX-Code
- <param name="decompDims"> 2 </param>
+
+\size small
+ <param name="decompDims"> 1 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<!-- Mesh size -->
\end_layout
\begin_layout LyX-Code
- <param name="meshSizeI"> 51 </param>
+
+\size small
+ <param name="meshSizeI"> 21 </param>
\end_layout
\begin_layout LyX-Code
- <param name="meshSizeJ"> 11 </param>
+
+\size small
+ <param name="meshSizeJ"> 6 </param>
\end_layout
\begin_layout LyX-Code
- <param name="meshSizeK"> 53 </param>
+
+\size small
+ <param name="meshSizeK"> 41 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<!-- Initial geometry -->
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="minX"> 0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="minY"> 0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="minZ"> 0 </param>
\end_layout
\begin_layout LyX-Code
- <param name="maxX"> 100000 </param>
+
+\size small
+ <param name="maxX"> 40000 </param>
\end_layout
\begin_layout LyX-Code
- <param name="maxY"> 20000 </param>
+
+\size small
+ <param name="maxY"> 10000 </param>
\end_layout
\begin_layout LyX-Code
- <param name="maxZ"> 104000 </param>
+
+\size small
+ <param name="maxZ"> 80000 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<!-- Remeshing -->
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="meshType"> cartesian </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="buildNodeNeighbourTbl"> True </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<!-- Elastic material parameters -->
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="lambda"> 1.0e+10 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="mu"> 1.0e+10 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<!-- Viscous material parameters -->
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="refvisc"> 1.0e+20 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="reftemp"> 1400.0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="activationE"> 45.0e+03 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="vis_min"> 1.0e+30 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="vis_max"> 1.0e+30 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="srexponent"> 1 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<!-- Plastic material parameters -->
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="yieldcriterion"> mohrcoulomb </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="nsegments"> 2 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="plstrain0"> 0.0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="plstrain1"> 0.02 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="plstrain2"> 1000.0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="frictionAngle0"> 30.0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="frictionAngle1"> 30.0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="frictionAngle2"> 30.0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="dilationAngle0"> 5.0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="dilationAngle1"> 5.0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="dilationAngle2"> 5.0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="cohesion0"> 4.0e+07 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="cohesion1"> 4.0e+05 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="cohesion2"> 0.0e+00 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="ten_off"> 1.0e+13 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<!-- Temperature variables -->
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="topTemp"> 0.0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="bottomTemp"> 700.0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="thermal_conduct"> 1.6 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="heatCapacity"> 1000.0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<!-- Remesher info -->
\end_layout
\begin_layout LyX-Code
+
+\size small
<!-- Used only when the SnacRemesher plugin is loaded.
-->
\end_layout
\begin_layout LyX-Code
+
+\size small
<!-- <param name="remeshCondition"> onBothTimeStepLength </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="remeshCondition"> onTimeStep </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="remeshCondition"> onMinLengthScale </param> -->
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="remeshCondition"> onMinLengthScale </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="remeshTimeStepCriterion"> 15000 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="remeshLengthCriterion"> 0.7 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="bottomResotre"> on </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<!-- node ICs -->
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct name="nodeICs">
\end_layout
\begin_layout LyX-Code
+
+\size small
<list name="vcList">
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> AllNodesVC </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<list name="variables">
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="name"> vx </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> double </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="value"> 0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="name"> vy </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> double </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="value"> 0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="name"> vz </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> double </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="value"> 0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
</list>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> AllNodesVC </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<list name="variables">
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="name"> temperature </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> func </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="value">SnacTemperature_Top2BottomSweep</param>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
</list>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
</list>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<!-- element ICs -->
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct name="elementICs">
\end_layout
\begin_layout LyX-Code
+
+\size small
<list name="vcList">
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> AllElementsVC </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<list name="variables">
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
- <paramname="name">elementMaterial</param>
+ <param name="name">elementMaterial</pa
+ram>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> int </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="value"> 0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
</list>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
</list>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<!-- Velocity BCs -->
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct name="velocityBCs">
\end_layout
\begin_layout LyX-Code
+
+\size small
<list name="vcList">
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> WallVC </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="wall"> left </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<list name="variables">
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="name"> vx </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> double </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="value"> -3.17e-10 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
</list>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> WallVC </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="wall"> right </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<list name="variables">
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="name"> vx </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> double </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="value"> 3.17e-10 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
</list>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> WallVC </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="wall"> back </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<list name="variables">
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="name"> vz </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> double </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="value"> 0.0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
</list>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> WallVC </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="wall"> front </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<list name="variables">
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="name"> vz </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> double </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="value"> 0.0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
</list>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> WallVC </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="wall"> bottom </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<list name="variables">
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
- <param name="name"> vx </param>
-\end_layout
-\begin_layout LyX-Code
- <param name="type"> func </param>
+\size small
+ <param name="name"> vy </param>
\end_layout
\begin_layout LyX-Code
- <param name="value"> variableVelBC </param>
-\end_layout
-\begin_layout LyX-Code
- </struct>
+\size small
+ <param name="type"> double </param>
\end_layout
\begin_layout LyX-Code
- <struct>
-\end_layout
-\begin_layout LyX-Code
- <param name="name"> vy </param>
+\size small
+ <param name="value"> 0.0 </param>
\end_layout
\begin_layout LyX-Code
- <param name="type"> double </param>
-\end_layout
-\begin_layout LyX-Code
- <param name="value"> 0.0 </param>
+\size small
+ </struct>
\end_layout
\begin_layout LyX-Code
- </struct>
-\end_layout
-\begin_layout LyX-Code
+\size small
</list>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
</list>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<!-- Temperature BCs -->
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct name="temperatureBCs">
\end_layout
\begin_layout LyX-Code
+
+\size small
<list name="vcList">
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> WallVC </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="wall"> top </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<list name="variables">
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="name"> temperature </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> double </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="value"> 0.0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
</list>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> WallVC </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="wall"> bottom </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<list name="variables">
\end_layout
\begin_layout LyX-Code
+
+\size small
<struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="name"> temperature </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="type"> double </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
<param name="value"> 700.0 </param>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
</list>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
</list>
\end_layout
\begin_layout LyX-Code
+
+\size small
</struct>
\end_layout
\begin_layout LyX-Code
+
+\size small
</StGermainData>
\end_layout
@@ -5380,41 +6038,6 @@
\end_layout
\begin_layout Chapter
-SNAC Cookbook
-\end_layout
-
-\begin_layout Section
-Introduction
-\end_layout
-
-\begin_layout Standard
-These cookbook examples are meant to serve as a guide to some of the different
- types of problems SNAC can solve.
- The Cookbook examples include both regional as well as problems in a full
- spherical shell.
- These examples address more or less traditional convection problems.
- For additional information and introduction to the code coupling features
- made possible by the Pyre framework, please refer to Chapter [TODO - which
- chapter to reference?].
-\end_layout
-
-\begin_layout Section
-Cookbook 1
-\end_layout
-
-\begin_layout Subsection
-Problem
-\end_layout
-
-\begin_layout Subsection
-Solution
-\end_layout
-
-\begin_layout Subsection
-Discussion
-\end_layout
-
-\begin_layout Chapter
Benchmark Problems
\end_layout
@@ -5436,11 +6059,11 @@
\end_layout
\begin_layout Itemize
-Bulk modulus = 2.0e8 Pa
+Bulk modulus = 200 MPa
\end_layout
\begin_layout Itemize
-Shear modulus = 2.0e8 Pa
+Shear modulus = 200 MPa
\end_layout
\begin_layout Itemize
@@ -5448,19 +6071,35 @@
\end_layout
\begin_layout Itemize
-friction angle = 10 degrees
+friction angle = 10
+\begin_inset Formula $^{\circ}$
+\end_inset
+
+
\end_layout
\begin_layout Itemize
-dilation angle = 10 degrees
+dilation angle = 10
+\begin_inset Formula $^{\circ}$
+\end_inset
+
+
\end_layout
\begin_layout Itemize
-tension cut-off = 5.67 Mpa
+tension cut-off = 5.67 MPa
\end_layout
\begin_layout Itemize
-Boundary Conditions: Vy = -1.0e-5 m/sec on the top surface.
+Boundary Conditions:
+\begin_inset Formula $V_{y}$
+\end_inset
+
+ =
+\begin_inset Formula $-10^{-5}$
+\end_inset
+
+ m/sec on the top surface.
All the other surfaces are confined to their normal directions.
\end_layout
@@ -5489,15 +6128,15 @@
\end_layout
\begin_layout Itemize
-Use this file to reproduce this result: input file
-\family typewriter
-plastic.xml
-\end_layout
+\begin_inset Formula $\sigma_{yy}$
+\end_inset
-\begin_layout Itemize
-S_yy component of stress and E_yy component of strain are monitored.
- Due to the transient waves, values were averaged over the top elements
- layer.
+ component of stress and
+\begin_inset Formula $\epsilon_{yy}$
+\end_inset
+
+ component of strain are monitored.
+ These values were averaged over the top elements layer.
\end_layout
@@ -5621,11 +6260,11 @@
\end_layout
\begin_layout Itemize
-Bulk modulus = 2.0e8 Pa
+Bulk modulus = 200 MPa
\end_layout
\begin_layout Itemize
-Shear modulus = 2.0e8 Pa
+Shear modulus = 200 MPa
\end_layout
\begin_layout Itemize
@@ -5633,15 +6272,23 @@
\end_layout
\begin_layout Itemize
-friction angle = 0 degrees
+friction angle = 0
+\begin_inset Formula $^{\circ}$
+\end_inset
+
+
\end_layout
\begin_layout Itemize
-dilation angle = 0 degrees
+dilation angle = 0
+\begin_inset Formula $^{\circ}$
+\end_inset
+
+
\end_layout
\begin_layout Itemize
-tension cut-off = 5.67e8 Pa
+tension cut-off = 567 MPa
\end_layout
\begin_layout Itemize
@@ -5650,8 +6297,12 @@
\end_layout
\begin_layout Itemize
-Boundary Conditions: P = 2.0e7 Pa at r = a, and P = 0.0 at r = b.
- Top and bttom: Vy=0.0, free-slip.
+Boundary Conditions: P = 20 MPa at r = a, and P = 0.0 at r = b.
+ Top and bttom:
+\begin_inset Formula $V_{y}$
+\end_inset
+
+=0.0, free-slip.
\end_layout
@@ -5815,11 +6466,11 @@
\end_layout
\begin_layout Itemize
-Bulk modulus = 2.0e8 Pa
+Bulk modulus = 200 MPa
\end_layout
\begin_layout Itemize
-Shear modulus = 2.0e8 Pa
+Shear modulus = 200 MPa
\end_layout
\begin_layout Itemize
@@ -5827,15 +6478,23 @@
\end_layout
\begin_layout Itemize
-friction angle = 10 degrees
+friction angle = 10
+\begin_inset Formula $^{\circ}$
+\end_inset
+
+
\end_layout
\begin_layout Itemize
-dilation angle = 10 degrees
+dilation angle = 10
+\begin_inset Formula $^{\circ}$
+\end_inset
+
+
\end_layout
\begin_layout Itemize
-tension cut-off = 5.67e8 Pa
+tension cut-off = 567 MPa
\end_layout
\begin_layout Itemize
@@ -5844,8 +6503,12 @@
\end_layout
\begin_layout Itemize
-Boundary Conditions: P = 2.0e6 Pa at r = a, and P = 0.0 at r = b.
- Top and bttom: Vy=0.0, free-slip.
+Boundary Conditions: P = 2 MPa at r = a, and P = 0.0 at r = b.
+ Top and bttom:
+\begin_inset Formula $V_{y}$
+\end_inset
+
+=0.0, free-slip.
\end_layout
@@ -6251,19 +6914,27 @@
\end_layout
\begin_layout Itemize
-Viscosity = 1.0e9 Pa.sec
+Viscosity =
+\begin_inset Formula $10^{9}$
+\end_inset
+
+ Pa.sec
\end_layout
\begin_layout Itemize
-Bulk modulus = 1.5e9 Pa
+Bulk modulus = 1.5 GPa
\end_layout
\begin_layout Itemize
-Shear modulus = 5.0e8 Pa
+Shear modulus = 500 MPa
\end_layout
\begin_layout Itemize
-V0 = 1.0e-4 m/sec
+V0 =
+\begin_inset Formula $10^{-4}$
+\end_inset
+
+ m/sec
\end_layout
\begin_layout Itemize
More information about the CIG-COMMITS
mailing list