Simulation of polydisperse multiphase systems using population balances and example application to bubbly flows

Abstract In this work the relationship between multiphase computational fluid dynamics models and population balance models is illustrated by deriving the main governing equations from the generalized population balance equation. The resulting set of equations, consisting of the well known two-fluid model coupled with a bivariate population balance model, is then implemented in the CFD code OpenFOAM. The implementation is used to simulate a particular multiphase problem: bubbly flow in a rectangular column. Results show that, although the different mesoscale models for drag force, coalescence, breakup and mass transfer, can be improved, the agreement with experiments is nevertheless good. Moreover, although the problem investigated is quite complex, as the evolution of bubbles is solved in real-space, time and phase-space (i.e. bubble size and composition) the resulting computational costs are reasonable. This is due to the fact that the bivariate population balance model is solved here with the so-called conditional quadrature method of moments, that very efficiently deals with these problems. The overall approach is demonstrated to be efficient and robust and is therefore suitable for the simulation of many polydisperse multiphase flows.

[1]  R. Clift,et al.  Bubbles, Drops, and Particles , 1978 .

[2]  G. Montante,et al.  Experimental Analysis and Computational Modelling of Gas–Liquid Stirred Vessels , 2007 .

[3]  C. Wen.  Mechanics of Fluidization , 1966 .

[4]  Daniele Marchisio,et al.  Computational Models for Polydisperse Particulate and Multiphase Systems , 2013 .

[5]  Tron Solberg,et al.  Large eddy simulation of the Gas-Liquid flow in a square cross-sectioned bubble column , 2001 .

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

[7]  A. Bakker,et al.  A COMPUTATIONAL MODEL FOR THE GAS-LIQUID FLOW IN STIRRED REACTORS , 1994 .

[8]  Alfredo Iranzo,et al.  Numerical simulation of the gas–liquid flow in a laboratory scale bubble column: Influence of bubble size distribution and non-drag forces , 2008 .

[9]  Robert F. Mudde,et al.  Two- and three-dimensional simulations of a bubble plume using a two-fluid model , 1999 .

[10]  D. Wright,et al.  Numerical advection of moments of the particle size distribution in Eulerian models , 2007 .

[11]  Robert McGraw,et al.  Chemically resolved aerosol dynamics for internal mixtures by the quadrature method of moments , 2003 .

[12]  Miguel A. Galán,et al.  Experimental study of the transition between unsteady flow regimes in a partially aerated two-dimensional bubble column , 2008 .

[13]  Donald S. Scott,et al.  An eddy cell model of mass transfer into the surface of a turbulent liquid , 1970 .

[14]  D. Marchisio,et al.  Bubble size distribution modeling in stirred gas–liquid reactors with QMOM augmented by a new correction algorithm , 2010 .

[15]  Massimo Morbidelli,et al.  Role of turbulent shear rate distribution in aggregation and breakage processes , 2006 .

[16]  Y. Liao,et al.  A literature review of theoretical models for drop and bubble breakup in turbulent dispersions , 2009 .

[17]  S. Becker,et al.  Modelling and simulation of the dynamic flow behaviour in a bubble column , 2001 .

[18]  Francis J. Doyle,et al.  Solution technique for a multi-dimensional population balance model describing granulation processes , 2005 .

[19]  Gerhart Eigenberger,et al.  Applicability of the standard k–ε turbulence model to the dynamic simulation of bubble columns: Part I. Detailed numerical simulations , 1999 .

[20]  Jitendra Kumar,et al.  The cell average technique for solving multi-dimensional aggregation population balance equations , 2008, Comput. Chem. Eng..

[21]  Daniele Marchisio,et al.  Simulation of coalescence, break-up and mass transfer in a gas-liquid stirred tank with CQMOM , 2013 .

[22]  Lawrence L. Tavlarides,et al.  Description of interaction processes in agitated liquid-liquid dispersions , 1977 .

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

[24]  Dongsheng Zhang,et al.  Eulerian modeling of reactive gas-liquid flow in a bubble column , 2007 .

[25]  Antonello Barresi,et al.  Implementation of the population balance equation in CFD codes for modelling soot formation in turbulent flames , 2006 .

[26]  Bengt Andersson,et al.  On the breakup of fluid particles in turbulent flows , 2006 .

[27]  Jayanta Chakraborty,et al.  A new framework for solution of multidimensional population balance equations , 2007 .

[28]  C. Yuan,et al.  Conditional quadrature method of moments for kinetic equations , 2011, J. Comput. Phys..

[29]  Rodney O. Fox,et al.  Optimal Moment Sets for Multivariate Direct Quadrature Method of Moments , 2009 .

[30]  S. Katz,et al.  Some problems in particle technology: A statistical mechanical formulation , 1964 .

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

[32]  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 .

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

[34]  Daniele Marchisio,et al.  Multivariate Quadrature-Based Moments Methods for turbulent polydisperse gas–liquid systems , 2013 .

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

[36]  M. Fialová,et al.  Duality of the gas-liquid flow regimes in bubble column reactors , 1997 .

[37]  Young Gul Hur,et al.  Origin of regime transition to turbulent flow in bubble column: Orifice- and column-induced transitions , 2013 .

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

[39]  Ville Alopaeus,et al.  Modelling local bubble size distributions in agitated vessels , 2007 .

[40]  Jam Hans Kuipers,et al.  Numerical simulation of the dynamic flow behavior in a bubble column : A study of closures for turbulence and interface forces , 2006 .

[41]  Ville Alopaeus,et al.  Validation of bubble breakage, coalescence and mass transfer models for gas–liquid dispersion in agitated vessel , 2006 .

[42]  R. G. Rice,et al.  Bubble size prediction for rigid and flexible spargers , 1991 .

[43]  Ville Alopaeus,et al.  Numerical solution of moment-transformed population balance equation with fixed quadrature points , 2006 .

[44]  D. E. Rosner,et al.  Multivariate Population Balances via Moment and Monte Carlo Simulation Methods: An Important Sol Reaction Engineering Bivariate Example and “Mixed” Moments for the Estimation of Deposition, Scavenging, and Optical Properties for Populations of Nonspherical Suspended Particles , 2003 .

[45]  A. Buffo,et al.  Multidimensional population balance model for the simulation of turbulent gas–liquid systems in stirred tank reactors , 2012 .

[46]  N. H. Thomas,et al.  Homogeneous–heterogeneous regime transition in bubble columns , 2001 .

[47]  Vanni,et al.  Approximate Population Balance Equations for Aggregation-Breakage Processes. , 2000, Journal of colloid and interface science.

[48]  Harvey D. Mendelson.  The prediction of bubble terminal velocities from wave theory , 1967 .

[49]  Gerhart Eigenberger,et al.  Applicability of the standard k–ε turbulence model to the dynamic simulation of bubble columns. Part II:: Comparison of detailed experiments and flow simulations , 1999 .

[50]  Vivek V. Ranade,et al.  Dynamics of gas-liquid flow in a rectangular bubble column: Experiments and single/multi-group CFD simulations , 2002 .

[51]  Elena María Cachaza Gianzo.  Hydrodynamics and mass transfer effects in bubble columns , 2011 .

[52]  M. Hounslow,et al.  A discretized population balance for nucleation, growth, and aggregation , 1988 .