[cig-commits] commit 2298 by bangerth to /var/svn/dealii/aspect

dealii.demon at gmail.com dealii.demon at gmail.com
Wed Feb 12 07:33:17 PST 2014


Revision 2298

Start adding something for the geophysical cookbooks section.

U   trunk/aspect/doc/manual/manual.tex


http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=2298&peg=2298

Diff:
Modified: trunk/aspect/doc/manual/manual.tex
===================================================================
--- trunk/aspect/doc/manual/manual.tex	2014-02-12 15:01:04 UTC (rev 2297)
+++ trunk/aspect/doc/manual/manual.tex	2014-02-12 15:33:15 UTC (rev 2298)
@@ -4110,12 +4110,62 @@
 
 \subsection{Geophysical setups}
 \label{sec:cookbooks-geophysical}
-\marginpar{To be written}
 
-Use cookbooks/periodic-box.prm
+Having gone through the ways in which one can set up problems in rectangular
+geometries, let us now move on to situations that are directed more towards the
+kinds of things we want to use spect{} for: the simulation of convection in
+the rocky mantles of planets or other celestial bodies.
 
-Include something that uses the GPlates interface
+To this end, we need to go through the list of issues that have to be described
+and that were outlined in Section~
ef{sec:cookbooks-overview}, and address them
+one by one:
+egin{itemize}
+  \item 	extit{What internal forces act on the medium (the equation)?}
+  	This may in fact be the most difficult to answer part of it all. The real
+  	material in Earth's mantle is certainly no Newtonian fluid where the stress
+  	is a linear function of the strain with a proportionality constant (the
+  	viscosity) $\eta$ that only depends on the temperature. Rather, the
+  	real viscosity almost surely also depends on the pressure and the strain
+  	rate. Because the issue is complicated and the exact material model not
+  	entirely clear, for the next few subsections we will therefore ignore the
+  	issue and start with just using the ``simple'' material model where the
+  	viscosity is constant and most other coefficients depend at most on the
+  	temperature.
+  	
+  \item 	extit{What external forces do we have (the right hand side)}
+  	There are of course other issues: for example, should the model include terms
+  	that describe shear heating? Should it be compressible? Adiabatic heating due
+  	to compression? Most of the terms that pertain to these questions appear on
+  	the right hand sides of the equations, though some (such as the
+  	compressibility) also affect the differential operators on the left. Either
+  	way, for the moment, let us just go with the simplest models and come back to
+  	the more advanced questions in later examples.
+  	
+  	One right hand side that will certainly be there is that due to gravitational
+  	acceleration. To first order, within the mantle gravity points radially
+  	inward and has a roughly constant magnitude. In reality, of course, the
+  	strength and direction of gravity depends on the distribution and density of
+  	materials in Earth -- and, consequently, on the solution of the model at
+  	every time step. We will discuss some of the associated issues in the
+  	examples below.
 
+  \item 	extit{What is the domain (geometry)?}
+  \item 	extit{What happens at the boundary for each variable involved
+  (boundary conditions)?}
+  \item 	extit{How did it look at the beginning (initial conditions)?}
+\end{itemize}
+
+units
+geometry
+
+\subsubsection{Simple convection in a quarter of an annulus}
+
+
+\subsubsection{Simple convection in a spherical shell}
+
+\marginpar{Use cookbooks/periodic-box.prm}
+\marginpar{Include something that uses the GPlates interface}
+
 \subsection{Benchmarks}
 \label{sec:cookbooks-benchmarks}
 


More information about the CIG-COMMITS mailing list