[cig-commits] r14533 - in long/3D/Gale/trunk: . documentation

[walter at geodynamics.org]

Author: walter
Date: 2009-03-29 17:06:10 -0700 (Sun, 29 Mar 2009)
New Revision: 14533

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/documentation/gale.lyx
Log:
 r2623 at dante:  boo | 2009-03-29 17:06:14 -0700
 More manual updates



Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
   - 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2620
   + 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2623

Modified: long/3D/Gale/trunk/documentation/gale.lyx
===================================================================
--- long/3D/Gale/trunk/documentation/gale.lyx	2009-03-29 23:49:27 UTC (rev 14532)
+++ long/3D/Gale/trunk/documentation/gale.lyx	2009-03-30 00:06:10 UTC (rev 14533)
@@ -234,8 +234,13 @@
  Computing (VPAC) and Louis Moresi's group at Monash University.
  Walter Landry of CIG and Luke Hodkinson of VPAC have enhanced the code
  to satisfy the requirements of the long-term tectonics community.
- Roger Buck, Gus Correa, and Robert Bialas from the Lamont-Doherty Earth
- Observatory, Columbia University, provided valuable user testing.
+ Roger Buck, Gus Correa, Robert Bialas, Guillaume Duclaux, John Sheehan,
+ Garrett Ito, 
+\color none
+Noah Fay, Neil de Laplante, Matthieu Quinquis, Patrice Rey, Lara O'Dwyer,
+ Louise Kellogg, and Laetitia Le Pourhiet
+\color inherit
+ provided valuable user testing.
  The Gale team requests that in your oral presentations and in your papers
  that you indicate your use of this code and acknowledge the authors of
  the code, 
@@ -1131,6 +1136,11 @@
 \end_layout
 
 \begin_layout Subsubsection
+\begin_inset LatexCommand label
+name "sub:Hydrostatic-Pressure"
+
+\end_inset
+
 Hydrostatic Pressure
 \end_layout
 
@@ -1257,8 +1267,8 @@
 \end_layout
 
 \begin_layout Standard
-So if you are simulating a simple isoviscous field, you do not even have
- to know what the hydrostatic pressure and forces are.
+Note that if you are simulating a simple viscous problem in a box, you never
+ even use the complete pressure.
  In the general case such as with a yielding rheology, the viscosity may
  depend on the overall pressure.
 \end_layout
@@ -7635,7 +7645,7 @@
  and making sure that the plugin is compiled into the Gale executable.
 \end_layout
 
-\begin_layout Subsection
+\begin_layout Subsubsection
 EulerDeform
 \begin_inset LatexCommand label
 name "sub:EulerDeform"
@@ -15348,13 +15358,347 @@
 \end_layout
 
 \begin_layout Section
+\begin_inset LatexCommand label
+name "sec:HydrostaticTerm"
+
+\end_inset
+
 HydrostaticTerm
 \end_layout
 
 \begin_layout Standard
+This term subtracts out the hydrostatic part of the pressure as described
+ in Section 
+\begin_inset LatexCommand ref
+reference "sub:Hydrostatic-Pressure"
 
+\end_inset
+
+.
+ It can subtract out the pressure from a two layer model with a temperature
+ profile described by 
+\family typewriter
+the TemperatureProfile
+\family default
+ standard condition function (see Section 
+\begin_inset LatexCommand ref
+reference "sec:Standard-Condition-Functions"
+
+\end_inset
+
+).
+ Specifically, it subtracts out a density given by
 \end_layout
 
+\begin_layout Standard
+\begin_inset Formula \[
+\begin{array}{cc}
+0 & x>height\\
+upperDensity\cdot\left(1-upperAlpha\cdot T\right) & materialBoundary>x>height\\
+lowerDensity\cdot\left(1-lowerAlpha\cdot T\right) & x>materialBoundary\end{array},\]
+
+\end_inset
+
+ where
+\end_layout
+
+\begin_layout Standard
+\align left
+\begin_inset Formula \[
+T=\left\{ \begin{array}{cc}
+T\_0 & x>height\\
+T\_0+linearCoefficient\left(height-x\right) & x<height\\
++exponentialCoefficient1\left(1-\exp\left(-exponentialCoefficient2\left(height-x\right)\right)\right)\end{array}\right..\]
+
+\end_inset
+
+ This component also computes a pressure by integrating this density profile
+ analytically and multiplying by 
+\begin_inset Formula $gravity.$
+\end_inset
+
+ Once you have created a 
+\family typewriter
+HydrostaticTerm
+\family default
+, you have to remember to pass it along to 
+\family typewriter
+BuoyancyForceTerm
+\family default
+ (Section 
+\begin_inset LatexCommand ref
+reference "sub:BouyancyForceTerm"
+
+\end_inset
+
+) in order for it to take effect.
+ You also have to pass it to any 
+\family typewriter
+DruckerPrager
+\family default
+ rheologies (Section 
+\begin_inset LatexCommand ref
+reference "sub:DruckerPrager"
+
+\end_inset
+
+) that you are using, so that the rheology will have the correct pressure.
+\end_layout
+
+\begin_layout Standard
+\align center
+\begin_inset Tabular
+<lyxtabular version="3" rows="12" columns="2">
+<features>
+<column alignment="center" valignment="top" leftline="true" width="0">
+<column alignment="center" valignment="top" leftline="true" rightline="true" width="0">
+<row topline="true" bottomline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+Defaults
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row topline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+upperDensity
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+0
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row topline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+upperAlpha
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+0
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row topline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+lowerDensity
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+0
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row topline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+lowerAlpha
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+0
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row topline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+height
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+0
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row topline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+materialBoundary
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+0
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row topline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+T_0
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+0
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row topline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+linearCoefficient
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+0
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row topline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+exponentialCoefficient1
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+0
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row topline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+exponentialCoefficient2
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+0
+\end_layout
+
+\end_inset
+</cell>
+</row>
+<row topline="true" bottomline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+gravity
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+0
+\end_layout
+
+\end_inset
+</cell>
+</row>
+</lyxtabular>
+
+\end_inset
+
+
+\end_layout
+
 \begin_layout Section
 Buoyancy Forces
 \begin_inset LatexCommand label
@@ -15449,6 +15793,18 @@
 \end_layout
 
 \begin_layout Standard
+If you have defined a 
+\family typewriter
+HydrostaticTerm
+\family default
+ (Section ), then you need to pass it in to the 
+\family typewriter
+BuoyancyForceTerm
+\family default
+ for it to take effect.
+\end_layout
+
+\begin_layout Standard
 \begin_inset VSpace defskip
 \end_inset
 
@@ -15459,7 +15815,7 @@
 \noindent
 \align center
 \begin_inset Tabular
-<lyxtabular version="3" rows="6" columns="2">
+<lyxtabular version="3" rows="7" columns="2">
 <features>
 <column alignment="center" valignment="top" leftline="true" width="0">
 <column alignment="center" valignment="top" leftline="true" rightline="true" width="0">
@@ -15583,6 +15939,26 @@
 \end_inset
 </cell>
 </row>
+<row bottomline="true">
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+HydrostaticTerm
+\end_layout
+
+\end_inset
+</cell>
+<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+\begin_inset Text
+
+\begin_layout Standard
+none
+\end_layout
+
+\end_inset
+</cell>
+</row>
 </lyxtabular>
 
 \end_inset
@@ -15658,6 +16034,18 @@
 \end_layout
 
 \begin_layout Standard
+Note that this component does not take a 
+\family typewriter
+HydrostaticTerm
+\family default
+, so you can not use this component if you are using 
+\family typewriter
+HydrostaticTerm
+\family default
+.
+\end_layout
+
+\begin_layout Standard
 \begin_inset VSpace defskip
 \end_inset