Numerical solution of the bivariate population balance equation for the interacting hydrodynamics and mass transfer in liquid¿liquid extraction columns

A comprehensive model for predicting the interacting hydrodynamics and mass transfer is formulated on the basis of a spatially distributed population balance equation in terms of the bivariate number density function with respect to droplet diameter and solute concentration. The two macro- (droplet breakage and coalescence) and micro- (interphase mass transfer) droplet phenomena are allowed to interact through the dispersion interfacial tension. The resulting model equations are composed of a system of partial and algebraic equations that are dominated by convection, and hence it calls for a specialized discretization approach. The model equations are applied to a laboratory segment of an RDC column using an experimentally validated droplet transport and interaction functions. Aside from the model spatial discretization, two methods for the discretization of the droplet diameter are extended to include the droplet solute concentration. These methods are the generalized fixed-pivot technique (GFP) and the quadrature method of moments (QMOM). The numerical results obtained from the two extended methods are almost identical, and the CPU time of both methods is found acceptable so that the two methods are being extended to simulate a full-scale liquid–liquid extraction column. 2005 Elsevier Ltd. All rights reserved.

[1]  Roy G. Gordon,et al.  Error Bounds in Equilibrium Statistical Mechanics , 1968 .

[2]  M. Manninen,et al.  On the mixture model for multiphase flow , 1996 .

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

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

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

[6]  Paul I. Barton,et al.  A Differentiation Index for Partial Differential-Algebraic Equations , 1999, SIAM J. Sci. Comput..

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

[8]  D. E. Rosner,et al.  Bivariate Extension of the Quadrature Method of Moments for Modeling Simultaneous Coagulation and Sintering of Particle Populations. , 2001, Journal of colloid and interface science.

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

[10]  N. P. Wilburn Mathematical Determination of Concentration Profiles in Two-Phase Continuous Countercurrent Extractors , 1964 .

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

[12]  Michael J. Hounslow,et al.  Finite-element methods for steady-state population balance equations , 1998 .

[13]  Solution of the droplet breakage equation for interacting liquid–liquid dispersions: a conservative discretization approach , 2004 .

[14]  D. E. Rosner,et al.  Bivariate moment simulation of coagulating and sintering nanoparticles in flames , 2002 .

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

[16]  Swati Mohanty,et al.  MODELING OF LIQUID-LIQUID EXTRACTION COLUMN: A REVIEW , 2000 .

[17]  J. Degrève,et al.  The Interaction of solute transfer, contaminants and dro break-up in rotating disc contactors: Part II. The coupling of the mass transfer and breakage processes via interfacial tension , 1997 .

[18]  V. Prakash MODELING OF LIQUID - LIQUID EXTRACTION COLUMN , 2004 .

[19]  D. R. Lewin,et al.  Modeling, simulation and control of liquid-liquid extraction columns , 1998 .

[20]  A. Kumar,et al.  Correlations for Prediction of Mass Transfer Coefficients in Single Drop Systems and Liquid–Liquid Extraction Columns , 1999 .

[21]  Hans-Jörg Bart,et al.  Numerical solution of the spatially distributed population balance equation describing the hydrodynamics of interacting liquid–liquid dispersions , 2004 .

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

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

[24]  Margaritis Kostoglou,et al.  An assessment of low-order methods for solving the breakage equation , 2002 .

[25]  D. Ramkrishna,et al.  Beyond log‐normal distributions: Hermite spectra for solving population balances , 2003 .

[26]  Thomas Baron,et al.  Mass and heat transfer from drops in liquid‐liquid extraction , 1957 .

[27]  Hans-Jörg Bart,et al.  The Droplet Population Balance Model – Estimation of Breakage and Coalescence , 2003 .

[28]  R. Diemer,et al.  A moment methodology for coagulation and breakage problems: Part 3—generalized daughter distribution functions , 2002 .

[29]  The generalized approximation method for modeling coagulation kinetics—Part 1: justification and implementation of the method , 2002 .

[30]  M. J. Slater,et al.  Liquid-liquid extraction equipment , 1994 .

[31]  Margaritis Kostoglou,et al.  Evaluation of Zero Order Methods for Simulating Particle Coagulation , 1994 .

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

[33]  Y. C. Zhou,et al.  A model for liquid‐liquid extraction column performance — The influence of drop size distribution on extraction efficiency , 1985 .