Implementation and analysis of numerical solution of the population balance equation in CFD packages

Simulation of polydisperse flows must include the effects of particle–particle interaction, as breakage and aggregation, coupling the population balance equation (PBE) with the multiphase modelling. In fact, the implementation of efficient and accurate new numerical techniques to solve the PBE is ne...

Full description

Saved in:
Bibliographic Details
Published in:Computers & chemical engineering 2008-12, Vol.32 (12), p.2933-2945
Main Authors: Silva, L.F.L.R., Damian, R.B., Lage, P.L.C.
Format: Article
Language:English
Subjects:
ANSYS CFX
Disperse flow
DQMOM
Multiphase flow
OpenFOAM
Population balance
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Simulation of polydisperse flows must include the effects of particle–particle interaction, as breakage and aggregation, coupling the population balance equation (PBE) with the multiphase modelling. In fact, the implementation of efficient and accurate new numerical techniques to solve the PBE is necessary. The direct quadrature method of moments, known as DQMOM, is a moment-based method that uses an optimal adaptive quadrature closure and came into view as a promising choice for this implementation. In the present work, DQMOM was implemented in two CFD packages: the commercial ANSYS CFX, through FORTRAN subroutines, and the open-source OpenFOAM, by directly coding the PBE solution. Transient zero-dimensional and steady one-dimensional simulations were performed in order to explore the PBE solution accuracy using several interpolation schemes. Simulation cases with dominant breakage, dominant aggregation and invariant solution (equivalent breakage and aggregation) were simulated and validated against an analytical solution. The solution of the population balance equation was then coupled to the two-fluid model, considering that all particles classes share the same velocity field. Momentum exchange terms were evaluated using the local instantaneous Sauter mean diameter of the size distribution function. The two-dimensional tests were performed in a backward facing step geometry where the vortex zones traps the particles and provides high rates of breakage and aggregation.
ISSN:0098-1354
1873-4375
DOI:10.1016/j.compchemeng.2008.03.007