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

[walter at geodynamics.org]

Author: walter
Date: 2007-11-09 16:21:34 -0800 (Fri, 09 Nov 2007)
New Revision: 8272

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/documentation/gale.lyx
Log:
 r1975 at earth:  boo | 2007-11-09 16:03:53 -0800
 Beginnings of documentation for the circular inclusion benchmark



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

Modified: long/3D/Gale/trunk/documentation/gale.lyx
===================================================================
--- long/3D/Gale/trunk/documentation/gale.lyx	2007-11-10 00:21:31 UTC (rev 8271)
+++ long/3D/Gale/trunk/documentation/gale.lyx	2007-11-10 00:21:34 UTC (rev 8272)
@@ -13558,6 +13558,91 @@
  and subduction models.
 \end_layout
 
+\begin_layout Section
+Circular Inclusion
+\end_layout
+
+\begin_layout Standard
+To check the results yoursel, run the accompanying input file with 3 different
+ resolutions
+\end_layout
+
+\begin_layout LyX-Code
+Gale spherical_inclusion.xml --outputPath=64 -pc_type lu -ksp_type preonly
+\end_layout
+
+\begin_layout LyX-Code
+Gale spherical_inclusion.xml --elementResI=128 --elementResJ=128 --outputPath=128
+ -pc_type lu -ksp_type preonly
+\end_layout
+
+\begin_layout LyX-Code
+Gale spherical_inclusion.xml --elementResI=256 --elementResJ=256 --outputPath=256
+ -pc_type lu -ksp_type preonly
+\end_layout
+
+\begin_layout Standard
+The last run will require about 800M of memory and take about 30 minutes
+ on a reasonable machine.
+ Then run paraview and load in the state file
+\end_layout
+
+\begin_layout LyX-Code
+spherical_inclusion.pvsm
+\end_layout
+
+\begin_layout Standard
+That will give you a picture like Figure 
+\begin_inset LatexCommand \ref{fig:Circular Inclusion error}
+
+\end_inset
+
+.
+ The high resolution error (blue) is less than the medium resolution error
+ (green), which is less than the low resolution error (red).
+\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/circular_inclusion_error.png
+	width 80col%
+	keepAspectRatio
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Caption
+\begin_inset LatexCommand \label{fig:Circular Inclusion error}
+
+\end_inset
+
+Numerical error for the circular inclusion benchmark around the interface
+ between the low and high viscosity regision.
+ Blue is high resolution (256x256), green is medium resolution (128x128),
+ and red is low resolution (64x64).
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+Far away from the hole, the error becomes dominated by factors besides finite
+ resolution (e.g.
+ finite distance to the boundary and non-zero tolerance when solving the
+ stokes flow).
+\end_layout
+
 \begin_layout Chapter
 License 
 \end_layout