Mathematical modelling and numerical simulation of two-phase gas-liquid flows in stirred-tank reactors

Abstract This paper presents the mathematical modelling and numerical simulation of the turbulent, two-phase flow of liquid and gas in a gas-induced agitated stirred-tank reactor, using Computational Fluid Dynamics (CFD) techniques. The reactor used as an application demonstration of the developed model is the ozone-induced one, first designed and modeled by Yang et al. (1999). A three-dimensional (3D), transient, Euler-Euler two-phase flow model is developed and used to investigate the turbulent flow and mixing of liquid and bubbles in the stirred-tank reactor, applying the sliding mesh approach. Turbulence is simulated by means of several available models, the Renormalization Group (RNG) k-e model being the one finally recommended as the most appropriate of the ones studied, for the present application. Two-way coupling between the two phases is modeled by means of appropriate inter-phase interaction relations. The study focused on bubbles of one size group (mean aerodynamic diameter of 2.5E-03 m), but it is easily extended to any number of sizes. It is concluded that the predicted overall flow field pattern and the mixing of both phases around the two blades of the simulated baffled stirred vessel, as well as inside and outside of the main tube of the reactor, are physically plausible, appear reasonably accurate, and are, therefore, satisfying.

[1]  Nikos C. Markatos.  Computer simulation of turbulent fluid flow in chemical reactors , 1983 .

[2]  D. P. Karadimou,et al.  A novel flow-oriented discretization scheme for reducing false diffusion in three-dimensional (3D) flows: An application in the indoor environment , 2012 .

[3]  Olivier Petit,et al.  The ERCOFTAC centrifugal pump OpenFOAM case-study , 2009 .

[4]  G. Eigenberger,et al.  Simulation of Buoyancy Driven Bubbly Flow: Established Simplifications and Open Questions (Journal Review) , 2004 .

[5]  Mamoru Ishii,et al.  Two-fluid model and hydrodynamic constitutive relations , 1984 .

[6]  S. Orszag,et al.  Renormalization group analysis of turbulence. I. Basic theory , 1986 .

[7]  N. C. Markatos,et al.  Mathematical Modelling of Single- and Two-Phase Flow Problems in the Process Industries , 1993 .

[8]  G. Barakos,et al.  Sliding mesh algorithm for CFD analysis of helicopter rotor–fuselage aerodynamics , 2008 .

[9]  Vivek V. Ranade,et al.  Computational Flow Modeling for Chemical Reactor Engineering , 2001 .

[10]  Yuri Bazilevs,et al.  3D simulation of wind turbine rotors at full scale. Part I: Geometry modeling and aerodynamics , 2011 .

[11]  V. Roig,et al.  Eulerian–Eulerian two-fluid model for turbulent gas–liquid bubbly flows , 2003 .

[12]  Christian B Allen,et al.  Development and Validation of Sliding and Non-Matching Grid Technology for Control Surface Representation , 2006 .

[13]  A. D. Gosman,et al.  Full flow field computation of mixing in baffled stirred vessels , 1993 .

[14]  A. I. Lygeros,et al.  Simulation and design of fluid-catalytic cracking riser-type reactors , 1996 .

[15]  J. T. Kuo,et al.  Flow of bubbles through nozzles , 1988 .

[16]  William Yang,et al.  Two-Phase CFD Model of the Bubble-Driven Flow in the Molten Electrolyte Layer of a Hall–Héroult Aluminum Cell , 2015, Metallurgical and Materials Transactions B.

[17]  N. Markatos,et al.  Modelling of two-phase transient flow and combustion of granular propellants , 1986 .

[18]  B. Launder,et al.  The numerical computation of turbulent flows , 1990 .

[19]  David L. Marcum,et al.  A sliding interface method for unsteady unstructured flow simulations , 2007 .

[20]  D. B. Spalding,et al.  A general purpose computer program for multi-dimensional one- and two-phase flow , 1981 .