Modelling of Bubbly Flow Using CFD-PBM Solver in OpenFOAM: Study of Local Population Balance Models and Extended Quadrature Method of Moments Applications

In order to optimize and design new bubbly flow reactors, it is necessary to predict the bubble behavior and properties with respect to the time and location. In gas-liquid flows, it is easily observed that the bubble sizes may vary widely. The bubble size distribution is relatively sharply defined, and bubble rises are uniform in homogeneous flow; however bubbles aggregate, and large bubbles are formed rapidly in heterogeneous flow. To assist in the analysis of these systems, the volume, size and other properties of dispersed bubbles can be described mathematically by distribution functions. Therefore, a mathematical modeling tool called the Population Balance Model (PBM) is required to predict the distribution functions of the bubble motion and the variation of their properties. In the present paper, two rectangular bubble columns and a water electrolysis reactor are modeled using the open-source Computational Fluid Dynamic (CFD) package OpenFOAM. Furthermore, the Method of Classes (CM) and Quadrature-based Moments Method (QBMM) are described, implemented and compared using the developed CFD-PBM solver. These PBM tools are applied in two bubbly flow cases: bubble columns (using a Eulerian-Eulerian two-phase approach to predict the flow) and a water electrolysis reactor (using a single-phase approach to predict the flow). The numerical results are compared with measured data available in the scientific literature. It is observed that the Extended Quadrature Method of Moments (EQMOM) leads to a slight improvement in the prediction of experimental measurements and provides a continuous reconstruction of the Number Density Function (NDF), which is helpful in the modeling of gas evolution electrodes in the water electrolysis reactor.

[1]  Pierre Proulx,et al.  Simulation of bubbly flows: Comparison between direct quadrature method of moments (DQMOM) and method of classes (CM) , 2010 .

[2]  P. Saffman,et al.  On the collision of drops in turbulent clouds , 1956, Journal of Fluid Mechanics.

[3]  H. Rusche.  Computational fluid dynamics of dispersed two-phase flows at high phase fractions , 2003 .

[4]  F. Laurent,et al.  An extended quadrature method of moments for population balance equations , 2003 .

[5]  W. Cochran.  The flow due to a rotating disc , 1934, Mathematical Proceedings of the Cambridge Philosophical Society.

[6]  Tron Solberg,et al.  Flow Generated by an Aerated Rushton Impeller: Two‐phase PIV Experiments and Numerical Simulations , 2008 .

[7]  Jos Derksen,et al.  Population Balance Modeling of Aerated Stirred Vessels Based on CFD , 2002 .

[8]  Pierre Proulx,et al.  Two-phase mass transfer coefficient prediction in stirred vessel with a CFD model , 2008, Comput. Chem. Eng..

[9]  D. Ramkrishna,et al.  ON THE SOLUTION OF POPULATION BALANCE EQUATIONS BY DISCRETIZATION--I . A FIXED PIVOT TECHNIQUE , 2003 .

[10]  Vivek V. Ranade,et al.  An efficient computational model for simulating flow in stirred vessels: a case of Rushton turbine , 1997 .

[11]  Zhengming Gao,et al.  Droplet breakage and coalescence in liquid–liquid dispersions: Comparison of different kernels with EQMOM and QMOM , 2017 .

[12]  Pierre Proulx,et al.  Three-dimensional mathematical modeling of dispersed two-phase flow using class method of population balance in bubble columns , 2008, Comput. Chem. Eng..

[13]  R. Fox,et al.  Application of the direct quadrature method of moments to polydisperse gas–solid fluidized beds , 2004 .

[14]  M. Leahy,et al.  Modeling Natural Convection in Copper Electrorefining: Describing Turbulence Behavior for Industrial-Sized Systems , 2011 .

[15]  Jay Sanyal,et al.  A generalized approach to model oxygen transfer in bioreactors using population balances and computational fluid dynamics , 2005 .

[16]  Daniele Marchisio,et al.  Implementation of the Quadrature Method of Moments in CFD codes for aggregation-breakage problems , 2003 .

[17]  R. I. Issa,et al.  Modelling of dispersed bubble and droplet flow at high phase fractions , 2004 .

[18]  Bjørn H. Hjertager,et al.  LDA measurements and CFD modelling of gas-liquid flow in a stirred vessel , 1996 .

[19]  M. Leahy,et al.  Experimental Validation of a Computational Fluid Dynamics Model of Copper Electrowinning , 2010 .

[20]  M. Leahy,et al.  Computational Fluid Dynamics Modelling of Natural Convection in Copper Electrorefining , 2007 .

[21]  R. Fox,et al.  Multiphase reacting flows: modelling and simulation , 2007 .

[22]  G. Litrico,et al.  Simulation of bubbly flow and mass transfer in a turbulent gas-liquid stirred tank with CFD-PBM solver in OpenFOAM: EQMOM application , 2018 .

[23]  H. Svendsen,et al.  Theoretical model for drop and bubble breakup in turbulent dispersions , 1996 .

[24]  Hallvard F. Svendsen,et al.  A coalescence and breakup module for implementation in CFD-codes , 2000 .

[25]  W. Schröder,et al.  Modeling and simulation of dispersed two-phase flow transport phenomena in electrochemical processes , 2009 .

[26]  Stefano Bove.  Computional Fluid Dynamics of Gas-Liquid Flows including Bubble Population Balances , 2005 .

[27]  H. Deconinck,et al.  Simulation of the Two-Phase Flow Hydrodynamics in an IRDE Reactor , 2009 .

[28]  Zoltan K. Nagy,et al.  Modelling of mass transfer in gas–liquid stirred tanks agitated by Rushton turbine and CD-6 impeller: A scale-up study , 2009 .

[29]  Marc Massot,et al.  Numerical simulation of spray coalescence in an Eulerian framework: Direct quadrature method of moments and multi-fluid method , 2007, J. Comput. Phys..

[30]  Boris V. Balakin,et al.  Coupling STAR-CD with a population-balance technique based on the classes method , 2014 .

[31]  Claude Dupuis,et al.  Mixing Characteristics of an Axial-Flow Rotor: Experimental and Numerical Study , 2005 .

[32]  W. J. Turner,et al.  Applicability of the standard k-ϵ turbulence model to gas-stirred baths , 1988 .

[33]  Ankur Gupta,et al.  Euler–Euler simulation of bubbly flow in a rectangular bubble column: Experimental validation with Radioactive Particle Tracking , 2013 .

[34]  L. Liu.  “Advances in Industrial Mixing: A Companion to the Handbook of Industrial Mixing” , 2017 .

[35]  V. Ranade,et al.  Eulerian-Lagrangian simulations of unsteady gas-liquid flows in bubble columns , 2006 .

[36]  Luiz Fernando L. R. Silva,et al.  Development and implementation of a polydispersed multiphase flow model in OpenFOAM , 2011, Comput. Chem. Eng..

[37]  H. Wagner,et al.  Hydrodynamic simulations of laboratory scale bubble columns fundamental studies of the Eulerian–Eulerian modelling approach , 1999 .

[38]  Gian Piero Celata,et al.  Terminal velocity of single bubbles in surface tension force dominant regime , 2002 .

[39]  A. Hubin,et al.  Modeling of mass and charge transfer in an inverted rotating disk electrode (IRDE) reactor , 2008 .

[40]  François Puel,et al.  Simulation and analysis of industrial crystallization processes through multidimensional population balance equations. Part 1: a resolution algorithm based on the method of classes , 2003 .

[41]  Vivek V. Ranade,et al.  Gas–liquid flow in stirred reactors: Trailing vortices and gas accumulation behind impeller blades , 1999 .

[42]  Physics of Electrolytic Gas Evolution , 2013, 1304.2701.

[43]  Investigation of Discrete Population Balance Models and Breakage Kernels for Dilute Emulsification Systems , 2011 .

[44]  N. Zuber,et al.  Drag coefficient and relative velocity in bubbly, droplet or particulate flows , 1979 .

[45]  Jesse T. Pikturna,et al.  Quadrature method of moments for population‐balance equations , 2003 .

[46]  Mohsen Karimi,et al.  PUFoam : A novel open-source CFD solver for the simulation of polyurethane foams , 2017, Comput. Phys. Commun..

[47]  Alberto Passalacqua,et al.  An extended quadrature‐based mass‐velocity moment model for polydisperse bubbly flows , 2014 .

[48]  D. Ramkrishna,et al.  On the solution of population balance equations by discretization—II. A moving pivot technique , 1996 .

[49]  M. P. Schwarz.  Simulation of gas injection into liquid melts , 1996 .

[50]  Ng Niels Deen.  An Experimental and Computational Study of Fluid Dynamics in Gas-Liquid Chemical Reactors , 2001 .

[51]  A. Bannari,et al.  CFD MODELING OF GAS DISPERSION AND BUBBLE SIZE IN A DOUBLE TURBINE STIRRED TANK , 2006 .

[52]  Daniele Marchisio,et al.  Solution of population balance equations using the direct quadrature method of moments , 2005 .

[53]  Mamoru Ishii,et al.  Interfacial area transport equation: model development and benchmark experiments , 2002 .

[54]  R. D. Vigil,et al.  Quadrature method of moments for aggregation-breakage processes. , 2003, Journal of colloid and interface science.

[55]  Marc Massot,et al.  Solution of population balance equations in applications with fine particles: Mathematical modeling and numerical schemes , 2016, J. Comput. Phys..

[56]  Robert McGraw,et al.  Description of Aerosol Dynamics by the Quadrature Method of Moments , 1997 .

[57]  Daniele Marchisio,et al.  On the Comparison between Population Balance Models for CFD Simulation of Bubble Columns , 2005 .

[58]  A. Passalacqua,et al.  An extended quadrature-based moment method with log-normal kernel density functions , 2015 .

[59]  A. C. Zettlemoyer,et al.  Homogeneous Nucleation Theory , 1974 .