diff --git a/docs/source/FluidEquations.rst b/docs/source/FluidEquations.rst
index c13ed9d525b470585768b84666876c30202e163c..92684f10b55826a1e659db644f16426714c6e512 100644
--- a/docs/source/FluidEquations.rst
+++ b/docs/source/FluidEquations.rst
@@ -1,6 +1,3 @@
-Here we describe the fluid variables, the governing equations, and the time discretization
-of the fluid evolution.
-
 Fluid Variables
 ===============
 
@@ -35,12 +32,3 @@ where :math:`\sum_p \beta_p (V_p - U_g)` is the drag term in which :math:`V_p` r
 Conservation of fluid volume:
 
 .. math:: \frac{\partial \varepsilon_g}{\partial t} + \nabla \cdot (\varepsilon_g  U_g)  = 0
-
-Time Discretization
-===============
-
-In the absence of reactions, we assume that the fluid density is unchanged.
-
-We compute the fluid volume fraction directly from the particle locations.
-
-Thus here we focus on the discretization of the momentum equation
diff --git a/docs/source/FluidTimeDiscretization.rst b/docs/source/FluidTimeDiscretization.rst
new file mode 100644
index 0000000000000000000000000000000000000000..fef179d0b87af682c9fdeffac99cb71a4d3a0ae8
--- /dev/null
+++ b/docs/source/FluidTimeDiscretization.rst
@@ -0,0 +1,41 @@
+
+Time Discretization
+===============
+
+In the absence of reactions, we assume that the fluid density is unchanged.
+
+We compute the fluid volume fraction directly from the particle locations.
+
+Thus here we focus on the discretization of the momentum equation
+
+In the predictor
+
+#. Define :math:`U^{MAC}`, the face-centered (staggered) MAC velocity which is used for advection.
+
+#. Define an approximation to the new-time state,:math:`(\varepsilon_g \rho_g U)^* = (\varepsilon_g \rho_g U)^n +  
+                                                          \Delta t ( -\nabla \cdot (\varepsilon_g \rho_g U^{MAC} U_g) 
+                                                                    + \varepsilon_g \nabla {p_g}^{n-1/2} + \nabla \cdot \tau^n
+                                                                    + \sum_{part} \beta_p (V_p - {U_g}^*) + \rho_g g )`
+
+#. Project :math:`U^*` by solving
+    :math:`\nabla \cdot \frac{\varepsilon_g}{\rho_g} \nabla \phi = \nabla \cdot (\varepsilon_g  U)^*`
+    then defining
+         :math:(\varepsilon_g  U)^{**} = (\varepsilon_g  U)^{*} - \frac{\varepsilon_g}{\rho_g} \nabla \phi
+    and 
+         :math:`{p_g}^{n+1/2,*} = {p_g}^{n-1/2} + \phi`  
+
+In the corrector
+
+#. Define an approximation to the new-time state,:math:`(\varepsilon_g \rho_g U)^{***} = (\varepsilon_g \rho_g U)^n +  
+                                                          \Delta t ( (-1/2) \nabla \cdot (\varepsilon_g \rho_g U^{MAC} U_g)^n 
+                                                                     -(1/2) \nabla \cdot (\varepsilon_g \rho_g U^{MAC} U_g)^{**} 
+                                                                    + \varepsilon_g \nabla {p_g}^{n+1/2,*} 
+                                                                    + (1/2) \nabla \cdot \tau^n + (1/2) \nabla \cdot \tau^{**}
+                                                                    + \sum_{part} \beta_p (V_p - {U_g}^{**}) + \rho_g g )`
+
+#. Project :math:`U^{***}` by solving
+    :math:`\nabla \cdot \frac{\varepsilon_g}{\rho_g} \nabla \phi = \nabla \cdot (\varepsilon_g  U)^{***}`
+    then defining
+         :math:(\varepsilon_g  U)^{n+1} = (\varepsilon_g  U)^{***} - \frac{\varepsilon_g}{\rho_g} \nabla \phi
+    and 
+         :math:`{p_g}^{n+1/2} = {p_g}^{n-1/2} + \phi`  
diff --git a/docs/source/Fluids.rst b/docs/source/Fluids.rst
index b8b2248a0a00b97c97c92d7db25d8d065f4e3d90..84294e7a72111ebe5c045eb2badb03dce0e16d4b 100644
--- a/docs/source/Fluids.rst
+++ b/docs/source/Fluids.rst
@@ -1,5 +1,8 @@
 .. _Chap:Fluids
 
+Here we describe the fluid variables, the governing equations, and the time discretization
+of the fluid evolution.
+
 Solving the Fluid Equations
 ===========================
 
@@ -7,4 +10,4 @@ Solving the Fluid Equations
    :maxdepth: 1
 
    FluidEquations
-
+   FluidTimeDiscretization