Introduction

Understanding and predicting bubbly flows is of critical importance in a large number of industrial applications, including boiling heat transfer in power plants, various metallurgical processes, and in bubble columns in the chemical industry. In bubble columns, used for partial oxidation of ethylene to acetaldehyde, isobutene separation, wet oxidation of heavily polluted effluent, and the production of synthetic fuels, for example, gas is injected at the bottom and, as the bubbles rise, the gas diffuses into the liquid and reacts (Furusaki et al., 2001). Bubble columns ranging from tens to hundreds of cubic meters are common in the chemical industry and up to thousands of cubic meters in bioreactors, where longer process times are needed. The absence of any moving parts and their relatively simple construction makes bubble columns particularly attractive for large-scale operations (Deckwer, 1992). Their operation, however, is usually dependent on the size of the vessel and the difficulty of scaling up small pilot models makes numerical predictions important. Similar considerations apply to other bubble systems.

Computational modeling of industrial-size multiphase flow systems must by necessity rely on models of the average flow. Such models range from simple mixture models to more sophisticated two-fluid models, where separate equations are solved for the dispersed and the continuous phase. Since no attempt is made to resolve the unsteady motion of individual bubbles, closure relations are necessary for the unresolved motion and the forces between the bubbles and the liquid.