Simulation of Elongated Bubbles in a Channel Using the Two-Fluid Model

Abstract

This paper investigates the capability of the two-fluid model to predict the bubble drift velocity of elongated bubbles in channels. The two-fluid model is widely used in the oil and gas industry for dynamic multiphase pipeline simulations. The bubble drift velocity is an important quantity in predicting pipeline flushing and slug flow. In this paper, it is shown that the two-fluid model in its standard form predicts a bubble drift velocity of (gH)^(1/2) (similar to the shallow water equations), instead of the exact value of 1/2(gH)^(1/2) as derived by Benjamin[1]. Modifying the two-fluid model with the commonly employed momentum correction parameter leads to a steady solution (in a moving reference frame), but still predicts an erroneous bubble drift velocity. To get the correct bubble drift velocity, it is necessary to include the pressure variation along the channel height due to both the hydrostatic component and the vertical momentum flux