A study on the numerical stability of the two-fluid model near ill-posedness

The two-fluid model is widely used in studying gas–liquid flow inside pipelines because it can qualitatively predict the flow field at low computational cost. However, the two-fluid model becomes ill-posed when the slip velocity exceeds a critical value, and computations can be quite unstable before the flow reaches the ill-posed condition. In this work, computational stability of various convection schemes together with the Euler implicit method for the time derivatives in conjunction with the two-fluid model is analyzed. A pressure correction algorithm for the two-fluid model is carefully implemented to minimize its effect on numerical stability. von Neumann stability analysis shows that the central difference scheme is more accurate and more stable than the 1st-order upwind, 2nd-order upwind, and QUICK schemes. The 2nd-order upwind scheme is much more susceptible to instability than the 1st-order upwind scheme and is inaccurate for short waves. Excellent agreement is obtained between the predicted and computed growth rates of harmonic disturbances. The instability associated with the two-fluid model discretized system of equations is related to but quantitatively different from the instability associated with ill-posedness of the two-fluid model. When the computation becomes unstable due to the ill-posedness, the machine roundoff errors from a selected range of short wavelengths, which scale with the grid size, are amplified rapidly to render the computation of any targeted long wavelength variation useless. For the viscous two-fluid model with wall friction and interfacial drag, a small-amplitude long wavelength disturbance grows due to viscous Kelvin–Helmholtz instability without triggering the grid scale short waves when the system remains well posed. Under such a condition, central difference is found to be the most accurate discretization scheme among those investigated.