3.1. Global balance in steady upward flow

Consider an upward two-phase flow in a vertical pipe (Fig. 3.1(a)). In a steady state the momentum balance is given by

(3.1)\[\begin{split}\begin{split} &( \rho_{G} v_{G}^{2} A_{G} + \rho_{L} v_{L}^{2} A_{L} )_{+} - ( \rho_{G} v_{G}^{2} A_{G} + \rho_{L} v_{L}^{2} A_{L} )_{-} \\ &= ( - p_{G} A_{G} -p_{L} A_{L} )_{+} - ( - p_{G} A_{G} -p_{L} A_{L} )_{-} - \tau_{G} S_{G} - \tau_{L} S_{L} - \rho_{G} g V_{G} - \rho_{L} g V_{L} \end{split}\end{split}\]

where \(A_{k}\) is the cross sectional area occupied by the phase \(k\), \(S_{k}\) is the area of the wall in contact with the phase \(k\), \(V_{k}\) is the volume of the phase \(k\), and \(g\) is the magnitude of the acceleration of gravity. With the volume fractions,

(3.2)\[\begin{split}\begin{split} &( \alpha_{G} \rho_{G} v_{G}^{2} + \alpha_{L} \rho_{L} v_{L}^{2} )_{+} - ( \alpha_{G} \rho_{G} v_{G}^{2} + \alpha_{L} \rho_{L} v_{L}^{2} )_{-} \\ &= ( - p_{G} \alpha_{G} -p_{L} \alpha_{L} )_{+} - ( - p_{G} \alpha_{G} -p_{L} \alpha_{L} )_{-} - \tau_{G} \frac{S_{G}}{A} - \tau_{L} \frac{S_{L}}{A} - \rho_{G} g \frac{V_{G}}{A} - \rho_{L} g \frac{V_{L}}{A} \end{split}\end{split}\]

By substituting \(A = \pi D^{2}/4\), \(S_{k} = Pe_{k} \Delta z\) and \(V_{k} = \alpha_{k} A \Delta z\), we obtain

(3.3)\[\begin{split}\begin{split} &( \alpha_{G} \rho_{G} v_{G}^{2} + \alpha_{L} \rho_{L} v_{L}^{2} )_{+} - ( \alpha_{G} \rho_{G} v_{G}^{2} + \alpha_{L} \rho_{L} v_{L}^{2} )_{-} \\ &= ( - p_{G} \alpha_{G} -p_{L} \alpha_{L} )_{+} - ( - p_{G} \alpha_{G} -p_{L} \alpha_{L} )_{-} - \frac{\tau_{G} Pe_{G} \Delta z}{A} - \frac{\tau_{L} Pe_{L} \Delta z}{A} - ( \alpha_{G} \rho_{G} + \alpha_{L} \rho_{L} ) g \Delta z \end{split}\end{split}\]

where \(Pe\) is the perimeter. Dividing the both sides by \(\Delta z\) and taking a limit of \(\Delta z \rightarrow 0\) yield

(3.4)\[\begin{split} \frac{\partial \alpha_{G} \rho_{G} v_{G}^{2}}{\partial z} + \frac{\partial \alpha_{L} \rho_{L} v_{L}^{2}}{\partial z} = - \frac{\partial \alpha_{G} p_{G}}{\partial z} - \frac{\partial \alpha_{L} p_{L}}{\partial z} - \frac{\tau_{G} Pe_{G}}{A} - \frac{\tau_{L} Pe_{L}}{A} - ( \alpha_{G} \rho_{G} + \alpha_{L} \rho_{L} ) g \end{split}\]

Employing a single-pressure assumption \(p_{G} = p_{L} = p\) gives

(3.5)\[\begin{split}\begin{split} - \frac{\partial p}{\partial z} &= \underbrace{\frac{\partial \alpha_{G} \rho_{G} v_{G}^{2}}{\partial z} + \frac{\partial \alpha_{L} \rho_{L} v_{L}^{2}}{\partial z}}_{\text{acceleration}} + \underbrace{\frac{\tau_{G} Pe_{G}}{A} + \frac{\tau_{L} Pe_{L}}{A}}_{\text{friction}} + \underbrace{( \alpha_{G} \rho_{G} + \alpha_{L} \rho_{L} ) g}_{\text{static}} \\ &= - \left. \frac{\partial p}{\partial z} \right|_{a} - \left. \frac{\partial p}{\partial z} \right|_{f} - \left. \frac{\partial p}{\partial z} \right|_{s} \end{split}\end{split}\]

where

(3.6)\[\begin{split}\begin{split} &- \left. \frac{\partial p}{\partial z} \right|_{a} = \frac{\partial \alpha_{G} \rho_{G} v_{G}^{2}}{\partial z} + \frac{\partial \alpha_{L} \rho_{L} v_{L}^{2}}{\partial z} \\ &- \left. \frac{\partial p}{\partial z} \right|_{f} = \frac{\tau_{G} Pe_{G}}{A} + \frac{\tau_{L} Pe_{L}}{A} \\ &- \left. \frac{\partial p}{\partial z} \right|_{s} = ( \alpha_{G} \rho_{G} + \alpha_{L} \rho_{L} ) g \end{split}\end{split}\]
../_images/two-phase-pressure-drop.png

Fig. 3.1 Vertical and horizontal two-phase pipe flows.