twoWayGPBEFoam: An open-source Eulerian QBMM solver for monokinetic bubbly flows

Abstract twoWayGPBEFoam is an open-source mesoscopic Eulerian QBMM solver for monokinetic bubbly flows. The solver is implemented within the OpenFOAM software framework. Unlike the existing macroscopic two-fluid model (TFM) solver twoPhaseEulerFoam , it can predict the size segregation phenomenon and the size-conditional velocities of the disperse phase, although it will not be able to predict the particle trajectory crossing (PTC). On theoretical grounds, the evolution of the disperse phase in multiphase flows is dictated by the generalized population balance equation (GPBE), which can be transformed into moment transport equations and solved using the finite-volume method with higher-order realizable spatial-discretization schemes and time-integration schemes. In order to address the closure problem of the size-conditional spatial flux, the size-conditional velocities need to be modeled. In many previous works, these are assumed to be identical with the disperse phase velocity predicted by the TFM. In this work, the size-conditional velocities were modeled using the velocity polynomial approximation (VPA), for which the velocity polynomial coefficients (VPCs) can be obtained from the moments themselves. By carrying out several test cases with both one-way and two-way coupling, we show that the results predicted by our solver agree well with the analytical solutions and the existing experimental data. Program summary Program Title: twoWayGPBEFoam and oneWayGPBEFoam Program Files doi: http://dx.doi.org/10.17632/rzstnw9ytw.1 Licensing provisions: GNU General Public License 3 Programming language: C++ Nature of problem: twoWayGPBEFoam and oneWayGPBEFoam have been developed to help investigate multiphase flows using the Eulerian QBMM. It provides an easily extended, parallelized, Eulerian QBMM environment. Solution method: The continuous phase is solved by the Eulerian approach. The disperse phase is solved by the QBMM. These equations are one-way coupled in oneWayGPBEFoam and two-way coupled in twoWayGPBEFoam . Additional comments including restrictions and unusual features: All appropriate methodological references are contained in the section entitled References.

[1]  Bo Kong,et al.  A solution algorithm for fluid-particle flows across all flow regimes , 2017, J. Comput. Phys..

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

[3]  Magne Hillestad,et al.  A multifluid-PBE model for simulation of mass transfer limited processes operated in bubble columns , 2018, Comput. Chem. Eng..

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

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

[6]  D. Lucas,et al.  Large eddy simulations of the gas–liquid flow in a rectangular bubble column , 2016 .

[7]  Ng Niels Deen,et al.  Particle mixing rates using the two-fluid model , 2018 .

[8]  A. Buffo,et al.  Simulation of a reacting gas–liquid bubbly flow with CFD and PBM: Validation with experiments ☆ , 2017 .

[9]  Zhengming Gao,et al.  Simulation of droplet breakage in turbulent liquid–liquid dispersions with CFD-PBM: Comparison of breakage kernels , 2016 .

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

[11]  J. J. J. Chen,et al.  Effect of liquid depth on circulation in bubble columns: a visual study , 1989 .

[12]  M. Bertodano,et al.  Effect of Collision Force on Well-Posedness and Stability of the Two-Fluid Model for Vertical Bubbly Flows , 2016 .

[13]  R. I. Issa,et al.  Numerical aspects of an algorithm for the Eulerian simulation of two‐phase flows , 2003 .

[14]  Olivier Desjardins,et al.  Verification of Eulerian–Eulerian and Eulerian–Lagrangian simulations for turbulent fluid–particle flows , 2017 .

[15]  A. D. Gosman,et al.  The computation of compressible and incompressible recirculating flows by a non-iterative implicit scheme , 1986 .

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

[17]  Daniele Marchisio,et al.  Simulation of polydisperse multiphase systems using population balances and example application to bubbly flows , 2013 .

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

[19]  Zhi Jian Wang,et al.  Realizable high-order finite-volume schemes for quadrature-based moment methods , 2011, J. Comput. Phys..

[20]  Eckhard Krepper,et al.  A population balance approach considering heat and mass transfer—Experiments and CFD simulations , 2011 .

[21]  Gerhart Eigenberger,et al.  Dynamic Numerical Simulation of Gas-Liquid Two Phase Flows, Euler-Euler versus Euler-Lagrange , 1997 .

[22]  Y. Liao,et al.  Application of new closure models for bubble coalescence and breakup to steam–water vertical pipe flow , 2014 .

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

[24]  Akio Tomiyama,et al.  On the role of the lateral lift force in poly-dispersed bubbly flows , 2011 .

[25]  Daniele Marchisio,et al.  On the implementation of moment transport equations in OpenFOAM: Boundedness and realizability , 2016 .

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

[27]  Marc Massot,et al.  Size-velocity correlations in hybrid high order moment/multi-fluid methods for polydisperse evaporating sprays: Modeling and numerical issues , 2013, J. Comput. Phys..

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

[29]  M. Bertodano,et al.  Stability and convergence of computational eulerian two-fluid model for a bubble plume , 2017 .

[30]  J. Carneiro,et al.  Large Eddy Simulation of particle–laden swirling flow with a presumed function method of moments , 2012 .

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

[32]  R'bei Bel Fdhila.  Analyse expérimentale et modélisation d'un écoulement vertical à bulles dans un élargissement brusque , 1991 .

[33]  Doraiswami Ramkrishna,et al.  Population Balances: Theory and Applications to Particulate Systems in Engineering , 2000 .

[34]  D. Drew.  Mathematical Modeling of Two-Phase Flow , 1983 .

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

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

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

[38]  Matthias Ullrich,et al.  Second-moment closure modeling of turbulent bubbly flows within the two-fluid model framework , 2017 .

[39]  Dongyue Li,et al.  Simulation of bubbly flows with special numerical treatments of the semi-conservative and fully conservative two-fluid model , 2017 .

[40]  P. J. O'rourke,et al.  The multiphase particle-in-cell (MP-PIC) method for dense particulate flows , 1996 .

[41]  Jeffrey C. Heylmun,et al.  An open-source quadrature-based population balance solver for OpenFOAM , 2018 .

[42]  Eckhard Krepper,et al.  The inhomogeneous MUSIG model for the simulation of polydispersed flows , 2008 .

[43]  Songze Chen,et al.  dugksFoam: An open source OpenFOAM solver for the Boltzmann model equation , 2017, Comput. Phys. Commun..

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

[45]  Alejandro Clausse,et al.  Two-Fluid Model Stability, Simulation and Chaos , 2016 .

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

[47]  D. Lucas,et al.  Comparison of Eulerian QBMM and classical Eulerian–Eulerian method for the simulation of polydisperse bubbly flows , 2019, AIChE Journal.

[48]  Olivier Simonin,et al.  Comparison between Lagrangian and mesoscopic Eulerian modelling approaches for inertial particles suspended in decaying isotropic turbulence , 2008, J. Comput. Phys..

[49]  Paola Lettieri,et al.  New quadrature‐based moment method for the mixing of inert polydisperse fluidized powders in commercial CFD codes , 2012 .

[50]  Y. Tsuji,et al.  Discrete particle simulation of two-dimensional fluidized bed , 1993 .

[51]  R. Fox,et al.  A fully coupled quadrature-based moment method for dilute to moderately dilute fluid–particle flows , 2010 .

[52]  Th. Frank,et al.  Extension of the inhomogeneous MUSIG model for bubble condensation , 2011 .