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

walter at geodynamics.org walter at geodynamics.org
Tue Jun 10 14:31:05 PDT 2008 

Author: walter
Date: 2008-06-10 14:31:04 -0700 (Tue, 10 Jun 2008)
New Revision: 12136

Added:
   long/3D/Gale/trunk/documentation/images/Extension_comparision.png
   long/3D/Gale/trunk/documentation/images/Extension_setup.png
   long/3D/Gale/trunk/documentation/images/Shortening_setup.png
   long/3D/Gale/trunk/documentation/images/extension256_5cm.png
   long/3D/Gale/trunk/documentation/images/extension512_2e3_5cm.png
Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/documentation/gale.lyx
Log:
 r2221 at earth:  boo | 2008-06-10 14:30:36 -0700
 Add complete extension benchmark and part of the shortening benchmark



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

Modified: long/3D/Gale/trunk/documentation/gale.lyx
===================================================================
--- long/3D/Gale/trunk/documentation/gale.lyx	2008-06-10 05:08:04 UTC (rev 12135)
+++ long/3D/Gale/trunk/documentation/gale.lyx	2008-06-10 21:31:04 UTC (rev 12136)
@@ -2299,6 +2299,11 @@
 \end_layout
 
 \begin_layout Section
+\begin_inset LatexCommand label
+name "sec:Gauging-Accuracy"
+
+\end_inset
+
 Gauging Accuracy
 \end_layout
 
@@ -16553,49 +16558,226 @@
 
 \end_inset
 
-: one was benchmark models extension, and the other models shortening.
- At this time, Gale has partially implemented the extension benchmark.
- The input file is in 
+: one benchmark models extension, and the other models shortening.
+ 
+\end_layout
+
+\begin_layout Subsection
+Extension
+\end_layout
+
+\begin_layout Standard
+This benchmark models a sandbox being extended as in Figure 
+\begin_inset LatexCommand ref
+reference "fig:Extension-model-setup"
+
+\end_inset
+
+.
+ The right side and half of the bottom are translated to the right.
+ This creates a velocity discontinuity at the center which is the initial
+ seed for localization.
+ Gale's implementation of this benchmark is in 
 \family typewriter
 input/benchmarks/extension.xml
 \family default
 .
- The result of running the benchmark, using a direct solver in parallel,
- is shown in Figure 
+ There are two simplifications made for the purpose of benchmarking:
+\end_layout
+
+\begin_layout Enumerate
+The boundaries are kinematic, rather than frictional.
+\end_layout
+
+\begin_layout Enumerate
+The velocity discontinuity is smoothed over 0.02 cm.
+\end_layout
+
+\begin_layout Standard
+These simplifications are shared by most of the codes participating in the
+ benchmark.
+ In addition, the exact background viscosity is not prescribed by the benchmark.
+ We have used 
+\begin_inset Formula $10^{12}Pa\cdot s$
+\end_inset
+
+, the same as I2ELVIS.
+ This value is somewhere in the middle of the range of what the various
+ codes used.
+\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/Extension_setup.png
+	lyxscale 50
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+\begin_inset Caption
+
+\begin_layout Standard
+\begin_inset LatexCommand label
+name "fig:Extension-model-setup"
+
+\end_inset
+
+Extension model setup
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+A comparison against the other codes is in Figure 
 \begin_inset LatexCommand ref
-reference "fig:Strain-Geomod"
+reference "fig:Comparison-extension"
 
 \end_inset
 
 .
+ While it is difficult to perform exact comparisons, Gale does produce a
+ similar pattern of faults.
+ It is interesting to note that this qualitative comparison holds true even
+ though the code is not convergent.
+ For example, Figure 
+\begin_inset LatexCommand ref
+reference "fig:Low-resolution-(256x32)"
+
+\end_inset
+
+ shows a low resolution run (256x32).
+ Moreover, Figure 
+\begin_inset LatexCommand ref
+reference "fig:High-accuracy-(nonlinear"
+
+\end_inset
+
+ shows the result of increasing the accuracy of the nonlinear solver (see
+ Section 
+\begin_inset LatexCommand ref
+reference "sec:Gauging-Accuracy"
+
+\end_inset
+
+) from the default of 
+\begin_inset Formula $10^{-2}$
+\end_inset
+
+ to 
+\begin_inset Formula $2\cdot10^{-3}$
+\end_inset
+
+.
+ In both cases, there is a difference in the details, but not in the character
+ of the solution.
+ However, changing the accuracy of the linear solver (see Section 
+\begin_inset LatexCommand ref
+reference "sec:Gauging-Accuracy"
+
+\end_inset
+
+) from the default of 
+\begin_inset Formula $10^{-5}$
+\end_inset
+
+ to 
+\begin_inset Formula $10^{-7}$
+\end_inset
+
+ produced only imperceptible changes in the solution.
 \end_layout
 
 \begin_layout Standard
-\noindent
+\begin_inset Float figure
+wide false
+sideways false
+status open
+
+\begin_layout Standard
 \align center
+\begin_inset Graphics
+	filename images/Extension_comparision.png
+	lyxscale 50
+	width 100col%
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+\begin_inset Caption
+
+\begin_layout Standard
+\begin_inset LatexCommand label
+name "fig:Comparison-extension"
+
+\end_inset
+
+Strain rate invariant for the numerical extension models after 5 cm of extension
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
 \begin_inset Float figure
-placement H
 wide false
 sideways false
 status open
 
 \begin_layout Standard
+\align center
 \begin_inset Graphics
-	filename images/strainrate_1cm.png
-	scale 50
+	filename images/extension256_5cm.png
+	lyxscale 50
+	width 90col%
 
 \end_inset
 
 
+\end_layout
+
+\begin_layout Standard
 \begin_inset Caption
 
 \begin_layout Standard
 \begin_inset LatexCommand label
-name "fig:Strain-Geomod"
+name "fig:Low-resolution-(256x32)"
 
 \end_inset
 
-Strain rate invariant after 1cm of extension
+Strain rate invariant of a low resolution (256x32) run of the extension
+ model after 5 cm of extension
 \end_layout
 
 \end_inset
@@ -16609,18 +16791,156 @@
 \end_layout
 
 \begin_layout Standard
-While it is difficult to perform exact comparisons, Gale does reproduce
- the multiple interlaced faults seen in both numerical and analogue experiments.
- However, the faults that Gale produces are not straight.
- This is probably because the viscosity we used for sand was on the low
- end (
-\begin_inset Formula $5\cdot10^{7}$
+\begin_inset Float figure
+wide false
+sideways false
+status open
+
+\begin_layout Standard
+\align center
+\begin_inset Graphics
+	filename images/extension512_2e3_5cm.png
+	lyxscale 50
+	width 90col%
+
 \end_inset
 
-).
- Using higher viscosities caused the solver to stop converging.
+
 \end_layout
 
+\begin_layout Standard
+\begin_inset Caption
+
+\begin_layout Standard
+\begin_inset LatexCommand label
+name "fig:High-accuracy-(nonlinear"
+
+\end_inset
+
+Strain rate invariant of a high accuracy (nonlinear tolerance=
+\begin_inset Formula $2\cdot10^{-3}$
+\end_inset
+
+) run of the extension model after 5 cm of extension
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Subsection
+Shortening
+\end_layout
+
+\begin_layout Standard
+This benchmark simulates a sandbox being shortened as in Figure 
+\begin_inset LatexCommand ref
+reference "fig:Shortening-model-setup"
+
+\end_inset
+
+.
+ The right side is moved to the left, creating a velocity discontinuity
+ at the bottom right corner.
+ Gale's implementation of this benchmark is in 
+\family typewriter
+input/benchmarks/shortening.xml
+\family default
+.
+ In this case, friction is critical to the behavior of the model.
+ Again, the exact background viscosity is not prescribed by the benchmark,
+ so we have used 
+\begin_inset Formula $10^{12}Pa\cdot s$
+\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/Shortening_setup.png
+	lyxscale 50
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+\begin_inset Caption
+
+\begin_layout Standard
+\begin_inset LatexCommand label
+name "fig:Shortening-model-setup"
+
+\end_inset
+
+Shortening model setup
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+A comparison against the other codes is in Figure .
+ Many of the same comments about the extension benchmark carry over.
+ Gale does produce a similar pattern of faults, although it is still not
+ convergent.
+ Figure compares low (256x32) and medium (512x64) resolution runs, and Figure
+ shows the result of increasing the accuracy of the nonlinear solver (see
+ Section 
+\begin_inset LatexCommand ref
+reference "sec:Gauging-Accuracy"
+
+\end_inset
+
+) from the default of 
+\begin_inset Formula $10^{-2}$
+\end_inset
+
+ to 
+\begin_inset Formula $2\cdot10^{-3}$
+\end_inset
+
+.
+ Also as before, changing the accuracy of the linear solver (see Section
+ 
+\begin_inset LatexCommand ref
+reference "sec:Gauging-Accuracy"
+
+\end_inset
+
+) from the default of 
+\begin_inset Formula $10^{-5}$
+\end_inset
+
+ to 
+\begin_inset Formula $10^{-7}$
+\end_inset
+
+ produced only imperceptible changes in the solution.
+\end_layout
+
 \begin_layout Chapter
 License 
 \end_layout

Added: long/3D/Gale/trunk/documentation/images/Extension_comparision.png
===================================================================
(Binary files differ)


Property changes on: long/3D/Gale/trunk/documentation/images/Extension_comparision.png
___________________________________________________________________
Name: svn:mime-type
   + application/octet-stream

Added: long/3D/Gale/trunk/documentation/images/Extension_setup.png
===================================================================
(Binary files differ)


Property changes on: long/3D/Gale/trunk/documentation/images/Extension_setup.png
___________________________________________________________________
Name: svn:mime-type
   + application/octet-stream

Added: long/3D/Gale/trunk/documentation/images/Shortening_setup.png
===================================================================
(Binary files differ)


Property changes on: long/3D/Gale/trunk/documentation/images/Shortening_setup.png
___________________________________________________________________
Name: svn:mime-type
   + application/octet-stream

Added: long/3D/Gale/trunk/documentation/images/extension256_5cm.png
===================================================================
(Binary files differ)


Property changes on: long/3D/Gale/trunk/documentation/images/extension256_5cm.png
___________________________________________________________________
Name: svn:mime-type
   + application/octet-stream

Added: long/3D/Gale/trunk/documentation/images/extension512_2e3_5cm.png
===================================================================
(Binary files differ)


Property changes on: long/3D/Gale/trunk/documentation/images/extension512_2e3_5cm.png
___________________________________________________________________
Name: svn:mime-type
   + application/octet-stream

More information about the cig-commits mailing list