A Comparative Study of Closure Relations for CFD Modelling of Bubbly Flow in a Vertical Pipe

Abstract

Commercial code CFX was used to examine the performance of a two-fluid model to predict the details of upward isothermal bubbly flow of air and water in a vertical pipe. The model equations are volume-averaged Navier-Stokes equations that require closure models for interfacial forces and bubble-induced turbulence effects. Two-equation SST and k-epsilon RANS turbulence models were also used. A parametric study of closure models included both standard options in CFX and previously published novel closure models that were implemented with user-defined functions. The CFD simulations were compared with two cases from the MTLoop experiments by Lucas et al. at the Helmholtz-Zentrum Dresden Rossendorf: one with wall-peak void fraction profile (MT039), and another with a core-peak void fraction profile (MT118). The effect of changing the drag force closures was not significant for the set examined. Poor predictions were found when the lift force and wall lubrication models were incompatible in magnitude. There was no significant effect of changing the liquid phase turbulence model. Changing the bubble-induced turbulence models, however, had a significant impact on the radial void fraction profile. The novel wall force from Lubchenko et al. at the Massachusetts Institute of Technology significantly improved the prediction of the near wall void fraction in the wall peak profile.