]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
State other possible avenues for future studies.
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 9 Jun 2008 22:12:59 +0000 (22:12 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 9 Jun 2008 22:12:59 +0000 (22:12 +0000)
git-svn-id: https://svn.dealii.org/trunk@16283 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-21/doc/intro.dox
deal.II/examples/step-21/doc/results.dox

index 07e2a4651cf019f96c9e2d27dfdf6a1ca3baa371..3d5d70a1811e61e538f50df9758dfb8adf9f9a45 100644 (file)
@@ -385,7 +385,9 @@ anything useful once we have negative saturations or ones larger than one.
 
 Note that we will have similar restrictions on the time step also in @ref
 step_23 "step-23" and @ref step_24 "step-24" where we solve the time dependent
-wave equation, another hyperbolic problem.
+wave equation, another hyperbolic problem. We will also come back to the issue
+of time step choice below in the section on <a href="#extensions">possible
+extensions to this program</a>.
 
 
 <h3>The test case</h3>
index 99f28c148af5237d5cceb26ddd0c2c0cdb480609..8e339da8b69d1a044d0493b3d63a87ecdf17e5c4 100644 (file)
@@ -156,6 +156,41 @@ Using such techniques, it can reasonably be expected that the solution process
 will be faster by at least an order of magnitude.
 
 
+<h4>Time stepping</h4>
+
+In the introduction we have identified the time step restriction 
+@f[
+  \triangle t_{n+1} \le \frac h{|\mathbf{u}^{n+1}(\mathbf{x})|}
+@f]
+that has to hold globally, i.e. for all $\mathbf x$. After discretization, we
+satisfy it by choosing
+@f[
+  \triangle t_{n+1} = \frac {\min_K h_K}{\max_{\mathbf{x}}|\mathbf{u}^{n+1}(\mathbf{x})|}.
+@f]
+
+This restriction on the time step is somewhat annoying: the finer we make the
+mesh the smaller the time step; in other words, we get punished twice: each
+time step is more expensive to solve and we have to do more time steps.
+
+This is particularly annoying since the majority of the additional work is
+spent solving the implicit part of the equations, i.e. the pressure-velocity
+system, whereas it is the hyperbolic transport equation for the saturation
+that imposes the time step restriction.
+
+To avoid this bottleneck, people have invented a number of approaches. For
+example, they may only re-compute the pressure-velocity field every few time
+steps (or, if you want, use different time step sizes for the
+pressure/velocity and saturation equations). This keeps the time step
+restriction on the cheap explicit part while it makes the solution of the
+implicit part less frequent. Experiments in this direction are
+certainly worthwhile; one starting point for such an approach is the paper by
+Zhangxin Chen, Guanren Huan and Baoyan Li: </i>An improved IMPES method for
+two-phase flow in porous media</i>, Transport in Porous Media, 54 (2004),
+pp. 361--376. There are certainly many other papers on this topic as well, but
+this one happened to land on our desk a while back.
+
+
+
 <h4>Adaptivity</h4>
 
 Adaptivity would also clearly help. Looking at the movies, one clearly sees

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.