CFD modeling of the mixing of water in oil emulsions

Abstract A computational fluid dynamics (CFD) model was developed for the mixing of water in oil emulsion in a lab-scale mixing tank equipped with a Rushton turbine impeller. Multiple reference frames (MRF) technique, k – ɛ model, and Eulerian–Eulerian approach were employed to model the impeller rotation, turbulence, and multiphase flow, respectively. The droplet size distribution within the mixing tank was estimated by means of the population balance approach, which employs the discrete method to describe coalescence and breakage of water droplets. To validate the CFD model, the cumulative probability size distribution computed using the model was compared with the experimentally determined values reported in the literature. This validated CFD model was then utilized to explore the effects of the impeller speed, oil type, and volume fraction of water on the cumulative probability droplet size distribution, number density, and distribution of local volume fraction of the dispersed phase.

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

[2]  Alberto Brucato,et al.  Numerical prediction of flow fields in baffled stirred vessels: A comparison of alternative modelling approaches , 1998 .

[3]  H. Neumann,et al.  Petroleum emulsions, micro-emulsions, and micellar solutions , 1988 .

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

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

[6]  Jiri Blazek,et al.  Computational Fluid Dynamics: Principles and Applications , 2001 .

[7]  J. P. Gupta,et al.  A model for transitional breakage probability of droplets in agitated lean liquid-liquid dispersions , 1979 .

[8]  S. Elghobashi.  Particle-laden turbulent flows: direct simulation and closure models , 1991 .

[9]  Arsène Isambert,et al.  Use of the capillary number for the prediction of droplet size in membrane emulsification , 2008 .

[10]  B. Launder,et al.  The numerical computation of turbulent flows , 1990 .

[11]  Margaritis Kostoglou,et al.  On the steady-state size distribution of dispersions in breakage processes , 1997 .

[12]  Joseph A. Wolkan,et al.  Introduction to probability and statistics , 1994 .

[13]  Steffen Schütz,et al.  Modeling fluid behavior and droplet interactions during liquid-liquid separation in hydrocyclones , 2009 .

[14]  Stuart E. Rogers,et al.  Steady and unsteady computation of impeller‐stirred reactors , 1996 .

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

[16]  Philip K. Chan,et al.  Using electrical resistance tomography and computational fluid dynamics modeling to study the formation of cavern in the mixing of pseudoplastic fluids possessing yield stress , 2008 .

[17]  Juan C. Lasheras,et al.  On the breakup of an air bubble injected into a fully developed turbulent flow. Part 1. Breakup frequency , 1999, Journal of Fluid Mechanics.

[18]  A. D. Gosman,et al.  Full flow field computation of mixing in baffled stirred vessels , 1993 .

[19]  Ashwin W. Patwardhan,et al.  Relation between hydrodynamics and drop size distributions in pump–mix mixer , 2010 .

[20]  J. F. Mitre,et al.  Comparison of the accuracy and performance of quadrature-based methods for population balance problems with simultaneous breakage and aggregation , 2010, Comput. Chem. Eng..

[21]  D. Spalding,et al.  A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows , 1972 .

[22]  Lawrence L. Tavlarides,et al.  The effect of continuous-phase viscosity on the unsteady state behavior of liquid-liquid agitated dispersions , 1987 .

[23]  Kohei Ogawa,et al.  Drop Size Distribution in Liquid-Liquid Mixing , 2003 .

[24]  Carlos A. Dorao,et al.  Analysis of breakage kernels for population balance modelling , 2009 .

[25]  M. Hounslow,et al.  Adjustable discretized population balance for growth and aggregation , 1995 .

[26]  Doraiswami Ramkrishna,et al.  Droplet breakage in stirred dispersions. Breakage functions from experimental drop-size distributions , 1996 .

[27]  Thomas Danner,et al.  Emulsification in turbulent flow 1. Mean and maximum drop diameters in inertial and viscous regimes. , 2007, Journal of colloid and interface science.

[28]  A. Gosman,et al.  PREDICTION OF IMPELLER- INDUCED FLOW IN MIXING VESSELS USING MULTIPLE FRAMES OF REFERENCE , 1994 .

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

[30]  Eric Loth,et al.  Quasi-steady shape and drag of deformable bubbles and drops , 2008 .

[31]  Gerald Warnecke,et al.  Solving population balance equations for two-component aggregation by a finite volume scheme , 2007 .

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

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

[34]  Brian Launder,et al.  SECOND-MOMENT CLOSURE AND ITS USE IN MODELLING TURBULENT INDUSTRIAL FLOWS , 1989 .

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

[36]  C. J. Hoogendoorn,et al.  Natural-convection flow in a square cavity calculated with low-Reynolds-number turbulence models , 1991 .

[37]  Haresh Patel,et al.  CFD analysis of mixing in thermal polymerization of styrene , 2010, Comput. Chem. Eng..

[38]  Hallvard F. Svendsen,et al.  A model for turbulent binary breakup of dispersed fluid particles , 2002 .

[39]  Kazumi Takagi,et al.  Analysis of water entrainment into dispersed W/O emulsion drops , 1988 .

[40]  Stefan Turek,et al.  Numerical aspects and implementation of population balance equations coupled with turbulent fluid dynamics , 2011, Comput. Chem. Eng..

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

[42]  Hiroshi Takeuchi,et al.  Interfacial Area of Liquid-Liquid Dispersion in a Mixer-Settler Extraction Columun , 1992 .

[43]  P. Becher,et al.  Encyclopedia of emulsion technology , 1983 .

[44]  Costas Kiparissides,et al.  Use of CFD in prediction of particle size distribution in suspension polymer reactors , 1998 .

[45]  Gad Hetsroni,et al.  Turbulence intensity in dilute two-phase flows—3 The particles-turbulence interaction in dilute two-phase flow , 1994 .

[46]  Hariprasad J. Subramani,et al.  Analysis of Formation of Water-in-Oil Emulsions , 2009 .

[47]  Jinfu Wang,et al.  Population Balance Model for Gas−Liquid Flows: Influence of Bubble Coalescence and Breakup Models , 2005 .

[48]  Michael F. Malone,et al.  Experimental studies and population balance equation models for breakage prediction of emulsion drop size distributions , 2009 .

[49]  Robert E. Treybal,et al.  Drop Formation in Two-Liquid-Phase Systems , 1950 .

[50]  Ka Ming Ng,et al.  New discretization procedure for the breakage equation , 1995 .

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

[52]  J. Alvarez,et al.  A population balance approach for the description of particle size distribution in suspension polymerization reactors , 1994 .

[53]  Jerzy Bałdyga,et al.  Drop break‐up in intermittent turbulence: Maximum stable and transient sizes of drops , 1998 .

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

[55]  Kazumi Takagi,et al.  Effects of viscosity on drop diameter of W/O emulsion dispersed in a stirred tank. , 1987 .

[56]  Costas Kiparissides,et al.  Prediction of particle size distribution in suspension polymerization reactors: effect of turbulence nonhomogeneity , 2000 .

[57]  W. Agterof,et al.  Prediction of emulsion particle sizes using a computational fluid dynamics approach , 2003 .

[58]  A. Randolph,et al.  Theory of Particulate Processes: Analysis and Techniques of Continuous Crystallization , 1971 .

[59]  Carolyn A. Koh,et al.  Measurement and Calibration of Droplet Size Distributions in Water-in-Oil Emulsions by Particle Video Microscope and a Focused Beam Reflectance Method , 2010 .

[60]  Luiz Fernando L. R. Silva,et al.  Implementation and analysis of numerical solution of the population balance equation in CFD packages , 2008, Comput. Chem. Eng..

[61]  M. C. Ruiz,et al.  Drop size distribution in a batch mixer under breakage conditions , 2002 .

[62]  Bakhtier Farouk,et al.  A numerical study of three-dimensional natural convection in a differentially heated cubical enclosure , 1991 .

[63]  James M. Church,et al.  Statistical Theories of Turbulence in Predicting Particle Size in Agitated Dispersions , 1960 .

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

[65]  Reuel Shinnar,et al.  On the behaviour of liquid dispersions in mixing vessels , 1961, Journal of Fluid Mechanics.

[66]  Mitsutoshi Nakajima,et al.  Effect of dispersed phase viscosity on maximum droplet generation frequency in microchannel emulsification using asymmetric straight-through channels , 2011 .

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

[68]  M. Jahoda,et al.  Comparison of Different Techniques for Modelling of Flow Field and Homogenization in Stirred Vessels , 2005 .

[69]  Haresh Patel,et al.  Computational Fluid Dynamics Study of a Styrene Polymerization Reactor , 2010 .

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

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