From: Wolfgang Bangerth <bangerth@colostate.edu>
Date: Mon, 19 Aug 2019 20:04:33 +0000 (-0600)
Subject: Further edits to the introduction of step-21.
X-Git-Tag: v9.2.0-rc1~1214^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=08ee1da5126a9cc629e8599cc5bca22ee674f21f;p=dealii.git

Further edits to the introduction of step-21.
---

diff --git a/examples/step-21/doc/intro.dox b/examples/step-21/doc/intro.dox
index 4fa98411d9..aeda981e08 100644
--- a/examples/step-21/doc/intro.dox
+++ b/examples/step-21/doc/intro.dox
@@ -77,21 +77,33 @@ similar). However, we have not said anything yet about the saturation, which
 of course is going to change as the fluids move around.
 
 The second part of the equations is the description of the
-dynamics of the saturation. The saturation equation for the displacing fluid (water) is:
+dynamics of the saturation, i.e., how the relative concentration of the
+two fluids changes with time. The saturation equation for the displacing 
+fluid (water) is given by the following conservation law:
 @f{eqnarray*}
   S_{t} + \nabla \cdot (F(S) \mathbf{u}) = q_{w},
-  \\
-  S_{t} + F(S) \nabla \mathbf{u} + \mathbf{u} \cdot \nabla F(S) = S_{t} + F(S) * q + \mathbf{u} \cdot \nabla F(S) = q_{w}.
 @f}
-where $q_{w}$ is the flow rate of the displacing fluid (water) and is related to the fractional flow F(S) through:
+which can be rewritten by using the product rule of the divergence operator
+in the previous equation:
+@f{eqnarray*}
+  S_{t} + F(S) \left[\nabla \cdot \mathbf{u}\right] 
+        + \mathbf{u} \cdot \left[ \nabla F(S)\right]
+  = S_{t} + F(S) q + \mathbf{u} \cdot \nabla F(S) = q_{w}.
+@f}
+Here, $q=\nabla\cdot \mathbf{u}$ is the total influx introduced
+above, and $q_{w}$ is the flow rate of the displacing fluid (water).
+These two are related to the fractional flow $F(S)$ in the following way:
+@f[
+  q_{w} = F(S) q,
+@f]
+where the fractional flow is often parameterized via the (heuristic) expression
 @f[
-  q_{w} = F(S) * q,
-  \\
   F(S)
   =
   \frac{k_{rw}(S)/\mu_{w}}{k_{rw}(S)/\mu_{w} + k_{ro}(S)/\mu_{o}}.
 @f]
-Thus, we obtain the saturation equation in the following advected form:
+Putting it all together yields the saturation equation in the following,
+advected form:
 @f{eqnarray*}
   S_{t} + \mathbf{u} \cdot \nabla F(S) = 0,
 @f}