[cig-commits] r14538 - in long/3D/Gale/trunk: . documentation documentation/images
walter at geodynamics.org
walter at geodynamics.org
Mon Mar 30 02:43:45 PDT 2009
Author: walter
Date: 2009-03-30 02:43:44 -0700 (Mon, 30 Mar 2009)
New Revision: 14538
Added:
long/3D/Gale/trunk/documentation/images/Paraview_hydro.png
Removed:
long/3D/Gale/trunk/documentation/images/Octave_fig3.png
Modified:
long/3D/Gale/trunk/
long/3D/Gale/trunk/documentation/gale.lyx
long/3D/Gale/trunk/documentation/images/Matlab_fig111.png
long/3D/Gale/trunk/documentation/images/Matlab_fig112.png
long/3D/Gale/trunk/documentation/images/Octave_fig1.png
long/3D/Gale/trunk/documentation/images/Octave_fig2.png
long/3D/Gale/trunk/documentation/images/Paraview_viscous_extension.png
long/3D/Gale/trunk/documentation/images/Sphere_Error.eps
long/3D/Gale/trunk/documentation/images/Sphere_Scaled_Error.eps
long/3D/Gale/trunk/documentation/images/tibet2D.png
long/3D/Gale/trunk/documentation/images/tibet3D.png
Log:
r2632 at dante: boo | 2009-03-30 02:37:30 -0700
Doc update
Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
- 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2629
+ 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2632
Modified: long/3D/Gale/trunk/documentation/gale.lyx
===================================================================
--- long/3D/Gale/trunk/documentation/gale.lyx 2009-03-30 08:13:11 UTC (rev 14537)
+++ long/3D/Gale/trunk/documentation/gale.lyx 2009-03-30 09:43:44 UTC (rev 14538)
@@ -1,4 +1,4 @@
-#LyX 1.5.1 created this file. For more info see http://www.lyx.org/
+#LyX 1.5.5 created this file. For more info see http://www.lyx.org/
\lyxformat 276
\begin_document
\begin_header
@@ -123,7 +123,7 @@
\end_layout
\begin_layout Standard
-A 3D simulation of a mid-ocean ridge.
+A 3D simulation of a mid-ocean ridge courtesy of Garrett Ito.
Shown is the strain rate, which highlights the region of active faulting
in the lower right.
\end_layout
@@ -1685,7 +1685,7 @@
\end_layout
\begin_layout Itemize
-Debian stable (x86 and AMD64), testing (x86), and unstable (x86)
+Debian stable (x86 and AMD64), testing (x86 and AMD64), and unstable (x86)
\end_layout
\begin_layout Standard
@@ -1809,7 +1809,6 @@
\begin_layout LyX-Code
\begin_inset LatexCommand htmlurl
-name "cig-long Archives"
target "geodynamics.org/pipermail/cig-long/"
\end_inset
@@ -1823,7 +1822,6 @@
\begin_layout LyX-Code
\begin_inset LatexCommand htmlurl
-name "Roundup"
target "geodynamics.org/roundup"
\end_inset
@@ -2854,23 +2852,6 @@
\end_inset
-\begin_inset VSpace defskip
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\noindent
-\align center
-\begin_inset Graphics
- filename images/Octave_fig3.png
- lyxscale 75
- scale 75
-
-\end_inset
-
-
\newpage
\end_layout
@@ -3527,7 +3508,7 @@
\begin_inset Graphics
filename images/Paraview_viscous_extension.png
lyxscale 75
- scale 70
+ width 100col%
\end_inset
@@ -3565,7 +3546,13 @@
\end_inset
).
- The magnitude of this structure is small, so we ignore it for now.
+ The magnitude of this structure is small, so you can ignore it for now.
+ Later, in Section , you will remove it by adding a
+\family typewriter
+HydrostaticTerm
+\family default
+.
+
\end_layout
\begin_layout Section
@@ -3807,7 +3794,7 @@
\begin_inset Graphics
filename images/Paraview_viscous_split.png
lyxscale 75
- scale 70
+ width 100col%
\end_inset
@@ -4220,7 +4207,7 @@
\begin_inset Graphics
filename images/Paraview_viscous_inflow.png
lyxscale 75
- scale 70
+ width 100col%
\end_inset
@@ -4480,7 +4467,7 @@
\align center
\begin_inset Graphics
filename images/Paraview_viscous_fixed.png
- scale 40
+ width 100col%
\end_inset
@@ -4882,7 +4869,7 @@
\begin_inset Graphics
filename images/Paraview_normal_stress.png
lyxscale 75
- scale 40
+ width 100col%
\end_inset
@@ -5081,7 +5068,7 @@
\begin_inset Graphics
filename images/Paraview_viscous_sinusoid.png
lyxscale 75
- scale 40
+ width 100col%
\end_inset
@@ -5111,6 +5098,222 @@
\end_layout
\begin_layout Section
+\begin_inset LatexCommand label
+name "sec:Hydrostatic-Term"
+
+\end_inset
+
+Hydrostatic Term
+\end_layout
+
+\begin_layout Standard
+This example will add a HydrostaticTerm to the simple extension model (Section
+
+\begin_inset LatexCommand ref
+reference "sec:Viscous-Material-in"
+
+\end_inset
+
+) to improve its accuracy.
+ So first copy
+\family typewriter
+myviscous_extension.xml
+\family default
+ to
+\family typewriter
+myhydrostatic.xml
+\family default
+.
+ Then add a
+\family typewriter
+HydrostaticTerm
+\family default
+ and
+\family typewriter
+StressBC
+\family default
+ just before the
+\family typewriter
+BuoyancyForceTerm
+\end_layout
+
+\begin_layout LyX-Code
+<struct name="hydrostaticTerm">
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="Type">HydrostaticTerm</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="upperDensity">1</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="height">maxY</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="gravity">gravity</param>
+\end_layout
+
+\begin_layout LyX-Code
+</struct>
+\end_layout
+
+\begin_layout LyX-Code
+<struct name="stressBC">
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="Type">StressBC</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="ForceVector">mom_force</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="Swarm">picIntegrationPoints</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="wall">top</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="y_type">HydrostaticTerm</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="y_value">hydrostaticTerm</param>
+\end_layout
+
+\begin_layout LyX-Code
+</struct>
+\end_layout
+
+\begin_layout Standard
+Then, in the
+\family typewriter
+BuoyancyForceTerm
+\family default
+, add an entry for
+\family typewriter
+HydrostaticTerm
+\family default
+ so that the
+\family typewriter
+BuoyancyForceTerm
+\family default
+ now looks like
+\end_layout
+
+\begin_layout LyX-Code
+<struct name="buoyancyForceTerm">
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="Type">BuoyancyForceTerm</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="ForceVector">mom_force</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="Swarm">picIntegrationPoints</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="gravity">gravity</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="HydrostaticTerm">hydrostaticTerm</param>
+\end_layout
+
+\begin_layout LyX-Code
+</struct>
+\end_layout
+
+\begin_layout Standard
+Figure
+\begin_inset LatexCommand ref
+reference "fig:Hydro_sri"
+
+\end_inset
+
+ shows the strain rate invariant with the same color scale as Figure
+\begin_inset LatexCommand ref
+reference "fig:Strain-viscous-extension"
+
+\end_inset
+
+.
+ With the
+\family typewriter
+HydrostaticTerm
+\family default
+, the artifacts seen at the top and bottom in Figure
+\begin_inset LatexCommand ref
+reference "fig:Strain-viscous-extension"
+
+\end_inset
+
+ disappear.
+\end_layout
+
+\begin_layout Standard
+\begin_inset Float figure
+wide false
+sideways false
+status open
+
+\begin_layout Standard
+\align center
+\begin_inset Graphics
+ filename images/Paraview_hydro.png
+ lyxscale 50
+ width 100col%
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+\begin_inset Caption
+
+\begin_layout Standard
+\begin_inset LatexCommand label
+name "fig:Hydro_sri"
+
+\end_inset
+
+Strain rate invariant for an extension model with the hydrostatic pressure
+ subtracted out.
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+A worked example is in
+\family typewriter
+input/cookbook/hydrostatic.xml
+\family default
+.
+\end_layout
+
+\begin_layout Section
Multiple Viscous Materials
\begin_inset LatexCommand label
name "sec:Multiple-Viscous-Materials"
@@ -5390,7 +5593,7 @@
\begin_inset Graphics
filename images/Paraview_multi_material.png
lyxscale 75
- scale 40
+ width 100col%
\end_inset
@@ -5440,7 +5643,7 @@
\begin_inset Graphics
filename images/Paraview_multi_viscosity.png
lyxscale 75
- scale 40
+ width 100col%
\end_inset
@@ -5720,7 +5923,7 @@
\begin_inset Graphics
filename images/Paraview_yielding.png
lyxscale 75
- scale 40
+ width 100col%
\end_inset
@@ -5775,7 +5978,7 @@
\begin_inset Graphics
filename images/Paraview_yielding_viscosity.png
lyxscale 75
- scale 40
+ width 100col%
\end_inset
@@ -5830,7 +6033,7 @@
\begin_inset Graphics
filename images/Paraview_yielding_yielding.png
lyxscale 75
- scale 40
+ width 100col%
\end_inset
@@ -6472,7 +6675,7 @@
\begin_inset Graphics
filename images/Paraview_thermal.png
lyxscale 75
- scale 40
+ width 100col%
\end_inset
@@ -6691,13 +6894,13 @@
\end_inset
- shows the strain rate invariant after the model has extended 30 km.
+ shows the strain rate invariant after the model has extended 5.4 km.
The resolution is 2048
\begin_inset Formula $\times$
\end_inset
512, and we used a direct solver (Mumps).
- The most prominent faults occur near where the crust thickens, although
+ The most prominent faults occur near the edge of the plateau, although
smaller faults occur throughout the crust.
The depth of the faults is limited by the relatively low viscosity deeper
in the crust.
@@ -6731,7 +6934,7 @@
\end_inset
-Strain rate invariant for the 2D Tibetan scale extension model after 30
+Strain rate invariant for the 2D Tibetan scale extension model after 5.4
km of extension.
\end_layout
@@ -6778,29 +6981,17 @@
.
It covers a region of the Tibetan plateau under extension.
+ The material properties are the same as in the 2D case.
\end_layout
\begin_layout Standard
-The material properties are the same as in the 2D case with one exception.
- The crust cohesion weakens to
-\begin_inset Formula $4\cdot10^{6}kg\, m^{-1}s^{-2}$
-\end_inset
-
- rather than
-\begin_inset Formula $4\cdot10^{5}kg\, m^{-1}s^{-2}$
-\end_inset
-
-.
-\end_layout
-
-\begin_layout Standard
Figure
\begin_inset LatexCommand ref
reference "fig:tibet3D"
\end_inset
- shows the strain rate invariant after the model has extended 24 km.
+ shows the strain rate invariant after the model has extended 16 km.
The resolution is 128
\begin_inset Formula $\times$
\end_inset
@@ -6841,7 +7032,7 @@
\end_inset
-Strain rate invariant for the 3D Tibetan scale extension model after 24
+Strain rate invariant for the 3D Tibetan scale extension model after 16
km of extension
\end_layout
@@ -7288,10 +7479,18 @@
\begin_layout LyX-Code
void StgFEM_StandardConditionFunctions_xSquared( Node_LocalIndex node_lI,
- Variable_Index var_I, void* _context, void* _result ) {
\end_layout
\begin_layout LyX-Code
+ Variable_Index var_I,
+\end_layout
+
+\begin_layout LyX-Code
+ void* _context, void* _result
+ ) {
+\end_layout
+
+\begin_layout LyX-Code
DiscretisationContext* context = (DiscretisationContext*)_co
ntext;
\end_layout
@@ -7325,11 +7524,15 @@
\end_layout
\begin_layout LyX-Code
- feVariable = (FeVariable*)FieldVariable_Register_GetByName( context->fie
-ldVariable_Register, "VelocityField" );
+ feVariable = (FeVariable*)FieldVariable_Register_GetByName(
\end_layout
\begin_layout LyX-Code
+ context->fieldVariable_Register, "VelocityFiel
+d" );
+\end_layout
+
+\begin_layout LyX-Code
mesh = feVariable->feMesh;
\end_layout
@@ -7365,10 +7568,13 @@
\begin_layout LyX-Code
condFunc = ConditionFunction_New( StgFEM_StandardConditionFunctions_Extension,
- "Velocity_Extension" );
\end_layout
\begin_layout LyX-Code
+ "Velocity_Extension" );
+\end_layout
+
+\begin_layout LyX-Code
ConditionFunction_Register_Add( context->condFunc_Register, condFunc );
\end_layout
@@ -7399,18 +7605,28 @@
\begin_layout LyX-Code
void StgFEM_StandardConditionFunctions_Extension( Node_LocalIndex node_lI,
- Variable_Index var_I, void* _context, void* _result ) ;
+ Variable_Index var_I,
\end_layout
+\begin_layout LyX-Code
+ void* _context, void*
+ _result ) ;
+\end_layout
+
\begin_layout Standard
add the line
\end_layout
\begin_layout LyX-Code
void StgFEM_StandardConditionFunctions_xSquared( Node_LocalIndex node_lI,
- Variable_Index var_I, void* _context, void* _result ) ;
+ Variable_Index var_I,
\end_layout
+\begin_layout LyX-Code
+ void* _context, void* _result
+ ) ;
+\end_layout
+
\begin_layout Standard
Now recompile the code, and you should be able to use your new function
from input files.
@@ -13668,8 +13884,17 @@
\begin_layout Standard
Finally, DruckerPrager requires a pressure.
- If you are using HydrostaticTerm (see Section ), you have to remember to
- give it that as well.
+ If you are using
+\family typewriter
+HydrostaticTerm
+\family default
+ (see Section
+\begin_inset LatexCommand ref
+reference "sec:HydrostaticTerm"
+
+\end_inset
+
+), you have to remember to give it that as well.
\end_layout
\begin_layout Standard
@@ -14279,7 +14504,7 @@
\end_layout
\begin_layout Standard
-Gale's computational domain is logically Eulerian.
+Gale's computational domain is logically Euclidean.
So in 2D, there are four boundaries:
\family typewriter
right
@@ -14720,7 +14945,7 @@
\end_layout
\begin_layout Standard
-If the nature of our problem is that stresses are specified on the boundary
+If the nature of your problem is that stresses are specified on the boundary
rather than velocities, you can specify those conditions using the
\family typewriter
StressBC
@@ -14904,6 +15129,63 @@
<param name="GaussianDim">0</param>
\end_layout
+\begin_layout Standard
+The type can also be
+\family typewriter
+HydrostaticTerm
+\family default
+.
+ If you are using a
+\family typewriter
+HydrostaticTerm
+\family default
+ component (see Section
+\begin_inset LatexCommand ref
+reference "sec:HydrostaticTerm"
+
+\end_inset
+
+), then you need a
+\family typewriter
+StressBC
+\family default
+ component on the top to act as a restoring force when the surface of the
+ material dips below equilibrium.
+ So it would be something like
+\end_layout
+
+\begin_layout LyX-Code
+<struct name="stressBCTop">
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="Type">StressBC</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="ForceVector">mom_force</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="Swarm">picIntegrationPoints</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="wall">top</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="y_type">HydrostaticTerm</param>
+\end_layout
+
+\begin_layout LyX-Code
+ <param name="y_value">hydrostaticTerm</param>
+\end_layout
+
+\begin_layout LyX-Code
+</struct>
+\end_layout
+
\begin_layout Subsection
Kinetic Friction Boundary Conditions
\end_layout
@@ -15424,6 +15706,16 @@
\end_inset
+) and a
+\family typewriter
+StressBC
+\family default
+ (Section
+\begin_inset LatexCommand ref
+reference "sub:Stress-Boundary-Conditions"
+
+\end_inset
+
) in order for it to take effect.
You also have to pass it to any
\family typewriter
@@ -17343,7 +17635,7 @@
\begin_layout Standard
\align center
\begin_inset Tabular
-<lyxtabular version="3" rows="14" columns="2">
+<lyxtabular version="3" rows="13" columns="2">
<features>
<column alignment="center" valignment="top" leftline="true" width="0">
<column alignment="center" valignment="top" leftline="true" rightline="true" width="0">
@@ -17392,14 +17684,8 @@
\end_inset
-,
+,
\begin_inset LatexCommand ref
-reference "sec:Shear"
-
-\end_inset
-
-,
-\begin_inset LatexCommand ref
reference "sec:Geomod-2004"
\end_inset
@@ -17420,7 +17706,7 @@
\end_inset
</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
\begin_inset Text
\begin_layout Standard
@@ -17460,14 +17746,8 @@
\end_inset
-,
+,
\begin_inset LatexCommand ref
-reference "sec:Shear"
-
-\end_inset
-
-,
-\begin_inset LatexCommand ref
reference "sec:Geomod-2004"
\end_inset
@@ -17488,7 +17768,7 @@
\end_inset
</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
\begin_inset Text
\begin_layout Standard
@@ -17534,14 +17814,8 @@
\end_inset
-,
+,
\begin_inset LatexCommand ref
-reference "sec:Shear"
-
-\end_inset
-
-,
-\begin_inset LatexCommand ref
reference "sec:Geomod-2004"
\end_inset
@@ -17599,7 +17873,7 @@
\end_inset
</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
+<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
\begin_inset Text
\begin_layout Standard
@@ -17666,12 +17940,6 @@
,
\begin_inset LatexCommand ref
-reference "sec:Shear"
-
-\end_inset
-
-,
-\begin_inset LatexCommand ref
reference "sec:Geomod-2004"
\end_inset
@@ -17707,13 +17975,7 @@
\end_inset
-,
-\begin_inset LatexCommand ref
-reference "sec:Shear"
-\end_inset
-
-
\end_layout
\end_inset
@@ -17750,7 +18012,7 @@
\end_inset
</cell>
</row>
-<row topline="true">
+<row topline="true" bottomline="true">
<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
\begin_inset Text
@@ -17775,37 +18037,6 @@
\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
-Friction boundary conditions
-\end_layout
-
-\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Standard
-\begin_inset LatexCommand ref
-reference "sec:Shear"
-
-\end_inset
-
-,
-\begin_inset LatexCommand ref
-reference "sec:Geomod-2004"
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-</cell>
-</row>
</lyxtabular>
\end_inset
@@ -17836,12 +18067,18 @@
\end_layout
\begin_layout Standard
-With the exception of the GeoMod 2004 benchmark (Section
+With the exception of the GeoMod benchmarks (Sections
\begin_inset LatexCommand ref
reference "sec:Geomod-2004"
\end_inset
+ and
+\begin_inset LatexCommand ref
+reference "sec:Geomod-2008"
+
+\end_inset
+
), the benchmarks can be carried out to high precision (~1%).
In particular, the error should follow the relation
\end_layout
@@ -18172,10 +18409,18 @@
\end_inset
.
- These were done with resolutions of 16
+ These were done with resolutions of 8
\begin_inset Formula $\times$
\end_inset
+16
+\begin_inset Formula $\times$
+\end_inset
+
+8, 16
+\begin_inset Formula $\times$
+\end_inset
+
32
\begin_inset Formula $\times$
\end_inset
@@ -18200,7 +18445,7 @@
\begin_inset Formula $h$
\end_inset
-) of 0.25, 0.125, and 0.0625.
+) of 0.5, 0.25, 0.125, and 0.0625.
Because of the symmetries of the problem we only have to simulate a quarter
of the domain.
As we increase the resolution (decrease
@@ -18307,8 +18552,7 @@
\end_inset
.
- So the medium-resolution error is multiplied by 2 and the high-resolution
- error is multiplied by 4.
+ So the higher resolution errors are multiplied by 2, 4, and 8.
\end_layout
\end_inset
@@ -18965,7 +19209,7 @@
\end_layout
\begin_layout Standard
-This benchmarks tests the implementation of the divergence term in equation
+This benchmark tests the implementation of the divergence term in equation
\begin_inset LatexCommand ref
reference "eq:divergence"
@@ -19042,6 +19286,7 @@
\align center
\begin_inset Graphics
filename images/divergence_v.png
+ width 100col%
\end_inset
@@ -19167,328 +19412,8 @@
\end_layout
-\begin_layout Section
-\begin_inset LatexCommand label
-name "sec:Shear"
-
-\end_inset
-
-Shear
-\end_layout
-
\begin_layout Standard
-This benchmark simulates a region filled with a viscous medium subjected
- to simple shear as in Figure
-\begin_inset LatexCommand ref
-reference "fig:shear_friction"
-\end_inset
-
-.
- Friction (static or kinetic) provides resistance on the bottom boundary.
- When the shear stress on the bottom boundary is less than the yield stress
- determined by static friction, then the material sticks to the bottom boundary.
- Otherwise, it starts to slide and kinetic friction determines the shear
- stress.
- Since this is a simple viscous material, the velocity profile is linear
- in both cases.
- This means that the strain rate invariant is constant everywhere.
-\end_layout
-
-\begin_layout Standard
-\begin_inset Float figure
-wide false
-sideways false
-status open
-
-\begin_layout Standard
-\align center
-\begin_inset Graphics
- filename images/shear.eps
- scale 50
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\begin_inset Caption
-
-\begin_layout Standard
-\begin_inset LatexCommand label
-name "fig:shear_friction"
-
-\end_inset
-
-Viscous medium under simple shear with frictional boundaries
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-The files
-\family typewriter
-input/benchmarks/shear_static_friction.xml
-\family default
- and
-\family typewriter
-input/benchmarks/shear_kinetic_
-\newline
-friction.xml
-\family default
- have the implementations of this benchmark.
- Figure
-\begin_inset LatexCommand ref
-reference "fig:Error-in-critical"
-
-\end_inset
-
- shows the error in the critical static friction coefficient vs.
- resolution.
- Figure
-\begin_inset LatexCommand ref
-reference "fig:Scaled-Error-in-critical"
-
-\end_inset
-
- shows the error scaled with resolution.
- The observed error decreases linearly with resolution, giving us confidence
- that the implementation of static friction is correct.
-\end_layout
-
-\begin_layout Standard
-\begin_inset Float figure
-wide false
-sideways false
-status open
-
-\begin_layout Standard
-\align center
-\begin_inset Graphics
- filename images/shear_static_error.eps
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\begin_inset Caption
-
-\begin_layout Standard
-\begin_inset LatexCommand label
-name "fig:Error-in-critical"
-
-\end_inset
-
-Error in critical friction coefficient
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\begin_inset Float figure
-wide false
-sideways false
-status open
-
-\begin_layout Standard
-\align center
-\begin_inset Graphics
- filename images/shear_static_scaled_error.eps
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\begin_inset Caption
-
-\begin_layout Standard
-\begin_inset LatexCommand label
-name "fig:Scaled-Error-in-critical"
-
-\end_inset
-
-Error in critical friction coefficient scaled with resolution
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-For pure kinetic friction, the analysis is more complicated.
- Figure
-\begin_inset LatexCommand ref
-reference "fig:Error-kinetic"
-
-\end_inset
-
- shows the error in the strain-rate invariant for three resolutions.
- As expected, the error decreases with resolution.
- Figure
-\begin_inset LatexCommand ref
-reference "fig:scaled-Error-kinetic"
-
-\end_inset
-
- shows the error scaled with resolution.
- We can see that the green medium-resolution error is always between the
- blue high-resolution error and the red low-resolution error.
- This means that the scaled high-resolution error is closer to the medium-resolu
-tion error than the low-resolution error.
- In other words, the error is following the expected relation
-\end_layout
-
-\begin_layout Standard
-\begin_inset Formula \[
-error\varpropto h+O(higher\, order\, terms\, in\, h).\]
-
-\end_inset
-
-This gives us confidence that the implementation of kinetic friction is
- correct.
-\end_layout
-
-\begin_layout Standard
-\begin_inset Float figure
-wide false
-sideways false
-status open
-
-\begin_layout Standard
-\align center
-\begin_inset Graphics
- filename images/shear_kinetic_error.png
- width 80col%
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\begin_inset Caption
-
-\begin_layout Standard
-\begin_inset LatexCommand label
-name "fig:Error-kinetic"
-
-\end_inset
-
-Error in the strain rate for the kinetic friction benchmark.
- Red is low resolution (16
-\begin_inset Formula $\times$
-\end_inset
-
-16), green is medium resolution (32
-\begin_inset Formula $\times$
-\end_inset
-
-32), and blue is high resolution (64
-\begin_inset Formula $\times$
-\end_inset
-
-64).
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\begin_inset Float figure
-wide false
-sideways false
-status open
-
-\begin_layout Standard
-\align center
-\begin_inset Graphics
- filename images/shear_kinetic_scaled_error.png
- width 80col%
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\begin_inset Caption
-
-\begin_layout Standard
-\begin_inset LatexCommand label
-name "fig:scaled-Error-kinetic"
-
-\end_inset
-
-Error in the strain rate scaled with resolution for the kinetic friction
- benchmark.
- Red is low resolution (16
-\begin_inset Formula $\times$
-\end_inset
-
-16), green is medium resolution (32
-\begin_inset Formula $\times$
-\end_inset
-
-32), and blue is high resolution (64
-\begin_inset Formula $\times$
-\end_inset
-
-64).
- The scaled high-resolution error is closer to the scaled medium-resolution
- error than the scaled low-resolution error, indicating proper linear convergenc
-e.
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-
\newpage
\end_layout
@@ -19687,7 +19612,7 @@
\begin_inset Graphics
filename images/Mohr_coulomb_resolutions.png
lyxscale 25
- scale 20
+ width 100col%
\end_inset
@@ -20263,13 +20188,18 @@
\end_layout
\begin_layout Section
+\begin_inset LatexCommand label
+name "sec:Geomod-2008"
+
+\end_inset
+
Geomod 2008
\end_layout
\begin_layout Standard
Using the lessons learned from the Geomod 2004 benchmarks, new benchmarks
were created that would make it easier to compare numerical experiments
- with each other and with analogue experiments
+ with each other and with analog experiments
\begin_inset LatexCommand cite
key "Geomod 2008"
@@ -20307,7 +20237,7 @@
has no cohesion.
So the force at the tip will go to zero as the thickness of the material
goes to zero.
- However, analogue experiments suggest a finite cohesion, so this benchmark
+ However, analog experiments suggest a finite cohesion, so this benchmark
specifies a boundary cohesion.
\end_layout
@@ -20355,9 +20285,9 @@
\end_inset
.
- If we do not weaken the cohesion, then when we try to model an unstable
- wedge by reducing the internal angle of friction, then the wedge will never
- collapse on itself.
+ If we do not weaken the cohesion, when we try to model an unstable wedge
+ by reducing the internal angle of friction, the wedge never collapses on
+ itself.
\end_layout
\begin_layout Standard
@@ -20938,8 +20868,9 @@
\end_layout
\begin_layout Standard
-This benchmark is very similar to unstable shortening, but the bottom right
- part translates along as well as shown in Figure
+This benchmark is very similar to unstable shortening.
+ The only difference is that part of the bottom is also moving along as
+ shown in Figure
\begin_inset LatexCommand ref
reference "fig:Brittle_setup"
@@ -22162,9 +22093,9 @@
Dohrmann, C., and
\color none
-Bochev, P.
- (2004), A stabilized finite element method for the Stokes problem based
- on polynomial pressure projections,
+P.
+ Bochev, (2004), A stabilized finite element method for the Stokes problem
+ based on polynomial pressure projections,
\shape italic
\color inherit
Int.
Modified: long/3D/Gale/trunk/documentation/images/Matlab_fig111.png
===================================================================
(Binary files differ)
Modified: long/3D/Gale/trunk/documentation/images/Matlab_fig112.png
===================================================================
(Binary files differ)
Modified: long/3D/Gale/trunk/documentation/images/Octave_fig1.png
===================================================================
(Binary files differ)
Modified: long/3D/Gale/trunk/documentation/images/Octave_fig2.png
===================================================================
(Binary files differ)
Deleted: long/3D/Gale/trunk/documentation/images/Octave_fig3.png
===================================================================
(Binary files differ)
Added: long/3D/Gale/trunk/documentation/images/Paraview_hydro.png
===================================================================
(Binary files differ)
Property changes on: long/3D/Gale/trunk/documentation/images/Paraview_hydro.png
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
Modified: long/3D/Gale/trunk/documentation/images/Paraview_viscous_extension.png
===================================================================
(Binary files differ)
Modified: long/3D/Gale/trunk/documentation/images/Sphere_Error.eps
===================================================================
--- long/3D/Gale/trunk/documentation/images/Sphere_Error.eps 2009-03-30 08:13:11 UTC (rev 14537)
+++ long/3D/Gale/trunk/documentation/images/Sphere_Error.eps 2009-03-30 09:43:44 UTC (rev 14538)
@@ -1,7 +1,7 @@
%!PS-Adobe-2.0 EPSF-2.0
%%Title: Sphere_Error.eps
%%Creator: gnuplot 4.2 patchlevel 2
-%%CreationDate: Thu Mar 20 13:17:23 2008
+%%CreationDate: Mon Mar 30 00:19:38 2009
%%DocumentFonts: (atend)
%%BoundingBox: 50 50 410 302
%%EndComments
@@ -48,7 +48,7 @@
/Author (Walter Landry,,,)
% /Producer (gnuplot)
% /Keywords ()
- /CreationDate (Thu Mar 20 13:17:23 2008)
+ /CreationDate (Mon Mar 30 00:19:38 2009)
/DOCINFO pdfmark
end
} ifelse
@@ -436,26 +436,18 @@
5849 0 R
-63 0 V
966 420 M
-( 4e-05) Rshow
+( 0) Rshow
1.000 UL
LTb
-1050 977 M
+1050 1056 M
63 0 V
5849 0 R
-63 0 V
-966 977 M
-( 6e-05) Rshow
-1.000 UL
-LTb
-1050 1533 M
-63 0 V
-5849 0 R
--63 0 V
-5933 0 R
-( 8e-05) Rshow
+( 5e-05) Rshow
1.000 UL
LTb
-1050 2090 M
+1050 1692 M
63 0 V
5849 0 R
-63 0 V
@@ -463,36 +455,36 @@
( 0.0001) Rshow
1.000 UL
LTb
-1050 2646 M
+1050 2328 M
63 0 V
5849 0 R
-63 0 V
-5933 0 R
-( 0.00012) Rshow
+( 0.00015) Rshow
1.000 UL
LTb
-1050 3203 M
+1050 2964 M
63 0 V
5849 0 R
-63 0 V
-5933 0 R
-( 0.00014) Rshow
+( 0.0002) Rshow
1.000 UL
LTb
-1050 3759 M
+1050 3600 M
63 0 V
5849 0 R
-63 0 V
-5933 0 R
-( 0.00016) Rshow
+( 0.00025) Rshow
1.000 UL
LTb
-1050 4316 M
+1050 4236 M
63 0 V
5849 0 R
-63 0 V
-5933 0 R
-( 0.00018) Rshow
+( 0.0003) Rshow
1.000 UL
LTb
1050 4872 M
@@ -500,7 +492,7 @@
5849 0 R
-63 0 V
-5933 0 R
-( 0.0002) Rshow
+( 0.00035) Rshow
1.000 UL
LTb
1050 420 M
@@ -508,10 +500,10 @@
0 4389 R
0 -63 V
0 -4529 R
-( 0.05) Cshow
+( 0) Cshow
1.000 UL
LTb
-2232 420 M
+2035 420 M
0 63 V
0 4389 R
0 -63 V
@@ -519,36 +511,44 @@
( 0.1) Cshow
1.000 UL
LTb
-3415 420 M
+3021 420 M
0 63 V
0 4389 R
0 -63 V
0 -4529 R
-( 0.15) Cshow
+( 0.2) Cshow
1.000 UL
LTb
-4597 420 M
+4006 420 M
0 63 V
0 4389 R
0 -63 V
0 -4529 R
-( 0.2) Cshow
+( 0.3) Cshow
1.000 UL
LTb
-5780 420 M
+4991 420 M
0 63 V
0 4389 R
0 -63 V
0 -4529 R
-( 0.25) Cshow
+( 0.4) Cshow
1.000 UL
LTb
+5977 420 M
+0 63 V
+0 4389 R
+0 -63 V
+0 -4529 R
+( 0.5) Cshow
+1.000 UL
+LTb
6962 420 M
0 63 V
0 4389 R
0 -63 V
0 -4529 R
-( 0.3) Cshow
+( 0.6) Cshow
1.000 UL
LTb
1.000 UL
@@ -575,27 +575,34 @@
1.000 UP
1.000 UL
LT0
-5780 4703 M
-0 -695 V
--31 695 R
+5977 3230 M
+0 1091 V
+5946 3230 M
62 0 V
--62 -695 R
+-62 1091 R
62 0 V
-2824 2449 M
-0 -556 V
--31 556 R
+3513 2181 M
+0 404 V
+-31 -404 R
62 0 V
--62 -556 R
+-62 404 R
62 0 V
-1346 1058 M
-0 -556 V
--31 556 R
+2282 1350 M
+0 476 V
+-31 -476 R
62 0 V
-1315 502 M
+-62 476 R
62 0 V
-5780 4355 Pls
-2824 2171 Pls
-1346 780 Pls
+1666 773 M
+0 461 V
+1635 773 M
+62 0 V
+-62 461 R
+62 0 V
+5977 3775 Pls
+3513 2383 Pls
+2282 1588 Pls
+1666 1004 Pls
1.000 UL
LTb
1050 4872 N
Modified: long/3D/Gale/trunk/documentation/images/Sphere_Scaled_Error.eps
===================================================================
--- long/3D/Gale/trunk/documentation/images/Sphere_Scaled_Error.eps 2009-03-30 08:13:11 UTC (rev 14537)
+++ long/3D/Gale/trunk/documentation/images/Sphere_Scaled_Error.eps 2009-03-30 09:43:44 UTC (rev 14538)
@@ -1,7 +1,7 @@
%!PS-Adobe-2.0 EPSF-2.0
%%Title: Sphere_Scaled_Error.eps
%%Creator: gnuplot 4.2 patchlevel 2
-%%CreationDate: Thu Mar 20 13:18:31 2008
+%%CreationDate: Mon Mar 30 00:19:38 2009
%%DocumentFonts: (atend)
%%BoundingBox: 50 50 410 302
%%EndComments
@@ -48,7 +48,7 @@
/Author (Walter Landry,,,)
% /Producer (gnuplot)
% /Keywords ()
- /CreationDate (Thu Mar 20 13:18:31 2008)
+ /CreationDate (Mon Mar 30 00:19:38 2009)
/DOCINFO pdfmark
end
} ifelse
@@ -436,34 +436,34 @@
5849 0 R
-63 0 V
966 420 M
-( 0.00016) Rshow
+( 0) Rshow
1.000 UL
LTb
-1050 865 M
+1050 1056 M
63 0 V
5849 0 R
-63 0 V
-966 865 M
-( 0.00017) Rshow
+-5933 0 R
+( 5e-05) Rshow
1.000 UL
LTb
-1050 1310 M
+1050 1692 M
63 0 V
5849 0 R
-63 0 V
-5933 0 R
-( 0.00018) Rshow
+( 0.0001) Rshow
1.000 UL
LTb
-1050 1756 M
+1050 2328 M
63 0 V
5849 0 R
-63 0 V
-5933 0 R
-( 0.00019) Rshow
+( 0.00015) Rshow
1.000 UL
LTb
-1050 2201 M
+1050 2964 M
63 0 V
5849 0 R
-63 0 V
@@ -471,52 +471,28 @@
( 0.0002) Rshow
1.000 UL
LTb
-1050 2646 M
+1050 3600 M
63 0 V
5849 0 R
-63 0 V
-5933 0 R
-( 0.00021) Rshow
+( 0.00025) Rshow
1.000 UL
LTb
-1050 3091 M
+1050 4236 M
63 0 V
5849 0 R
-63 0 V
-5933 0 R
-( 0.00022) Rshow
+( 0.0003) Rshow
1.000 UL
LTb
-1050 3536 M
-63 0 V
-5849 0 R
--63 0 V
--5933 0 R
-( 0.00023) Rshow
-1.000 UL
-LTb
-1050 3982 M
-63 0 V
-5849 0 R
--63 0 V
--5933 0 R
-( 0.00024) Rshow
-1.000 UL
-LTb
-1050 4427 M
-63 0 V
-5849 0 R
--63 0 V
--5933 0 R
-( 0.00025) Rshow
-1.000 UL
-LTb
1050 4872 M
63 0 V
5849 0 R
-63 0 V
-5933 0 R
-( 0.00026) Rshow
+( 0.00035) Rshow
1.000 UL
LTb
1050 420 M
@@ -524,10 +500,10 @@
0 4389 R
0 -63 V
0 -4529 R
-( 0.05) Cshow
+( 0) Cshow
1.000 UL
LTb
-2232 420 M
+2035 420 M
0 63 V
0 4389 R
0 -63 V
@@ -535,36 +511,44 @@
( 0.1) Cshow
1.000 UL
LTb
-3415 420 M
+3021 420 M
0 63 V
0 4389 R
0 -63 V
0 -4529 R
-( 0.15) Cshow
+( 0.2) Cshow
1.000 UL
LTb
-4597 420 M
+4006 420 M
0 63 V
0 4389 R
0 -63 V
0 -4529 R
-( 0.2) Cshow
+( 0.3) Cshow
1.000 UL
LTb
-5780 420 M
+4991 420 M
0 63 V
0 4389 R
0 -63 V
0 -4529 R
-( 0.25) Cshow
+( 0.4) Cshow
1.000 UL
LTb
+5977 420 M
+0 63 V
+0 4389 R
+0 -63 V
+0 -4529 R
+( 0.5) Cshow
+1.000 UL
+LTb
6962 420 M
0 63 V
0 4389 R
0 -63 V
0 -4529 R
-( 0.3) Cshow
+( 0.6) Cshow
1.000 UL
LTb
1.000 UL
@@ -591,27 +575,34 @@
1.000 UP
1.000 UL
LT0
-5780 1931 M
-0 -1113 V
--31 1113 R
+5977 1552 M
+0 1091 V
+5946 1552 M
62 0 V
-5749 818 M
+-62 1091 R
62 0 V
-2824 3353 M
-0 -1781 V
--31 1781 R
+3513 1979 M
+0 807 V
+-31 -807 R
62 0 V
-2793 1572 M
+-62 807 R
62 0 V
-1346 4504 M
-0 -3561 V
--31 3561 R
+2282 1803 M
+0 1905 V
+2251 1803 M
62 0 V
-1315 943 M
+-62 1905 R
62 0 V
-5780 1374 Pls
-2824 2462 Pls
-1346 2724 Pls
+1666 913 M
+0 3684 V
+1635 913 M
+62 0 V
+-62 3684 R
+62 0 V
+5977 2098 Pls
+3513 2383 Pls
+2282 2755 Pls
+1666 2755 Pls
1.000 UL
LTb
1050 4872 N
Modified: long/3D/Gale/trunk/documentation/images/tibet2D.png
===================================================================
(Binary files differ)
Modified: long/3D/Gale/trunk/documentation/images/tibet3D.png
===================================================================
(Binary files differ)
More information about the CIG-COMMITS
mailing list