Coalescence of drops in stirred dispersion. A white noise model for coalescence

Coalescence between two droplets in a turbulent liquid-liquid dispersion is generally viewed as a consequence of forces exerted on the drop-pair squeezing out the intervening continuous phase to a critical thickness. A new synthesis is proposed herein which models the film drainage as a stochastic process driven by a suitably idealized random process for the fluctuating force. While the true test of the model lies in detailed parameter estimations with measurement of drop-size distributions in coalescing dispersions, experimental measurements on average coalescence frequencies lend preliminary support to the model.

Breakage of viscoelastic drops in turbulent stirred dispersions

The existing models of drop breakage in stirred turbulent dispersions are applicable only to purely viscous dispersed phases. In their present form, they are found to underpredict the diameters of the largest stable drops formed when a viscoelastic fluid is dispersed into a Newtonian liquid. In purely viscous fluids, the turbulent stresses are opposed both by the stresses due to interfacial tension and the viscous stresses generated as the drop deforms. In viscoelastic fluids, drop deformation produces additional retractive elastic stresses which also oppose turbulent stresses. As the deformation rates are large, the retractive stresses can be large in magnitude. Assuming that these additional stresses decay with time, a model of viscoelastic drop breakage in turbulent stirred dispersions has been developed. The new model quantitatively predicts the dmax of viscoelastic fluids. The model, however, does not predict the observation that when the time constant of the fluid becomes large (λ > 0.5 s), the fluid can not be dispersed into droplets up to agitator speeds of about 10 rps in our equipment.

Coalescence of rigid droplets in a stirred dispersion-II. Band-limited force fluctuations

The coalescence of nearly rigid liquid droplets in a turbulent flow field is viewed as the drainage of a thin film of liquid under the action of a stochastic force representing the effect of turbulence. The force squeezing the drop pair is modelled as a correlated random function of time. The drops are assumed to coalesce once the film thickness becomes smaller than a critical thickness while they are regarded as separated if their distance of separation is larger than a prescribed distance. A semi-analytical solution is derived to determine the coalescence efficiency. The veracity of the solution procedure is established via a Monte-Carlo solution scheme. The model predicts a reversing trend of the dependence of the coalescence efficiency on the drop radii, the film liquid viscosity and the turbulence energy dissipation per unit mass, as the relative fluctuation increases. However, the dependence on physical parameters is weak (especially at high relative fluctuation) so that for the smallest droplets (which are nearly rigid) the coalescence efficiency may be treated as an empirical constant. The predictions of this model are compared with those of a white-noise force model. The results of this paper and those in Muralidhar and Ramkrishna (1986, Ind. Engng Chem. Fundam. 25, 554-56) suggest that dynamic drop deformation is the key factor that influences the coalescence efficiency.

Gas-liquid interfacial area in stirred vessels: The effect of an immiscible liquid phase

The presence of an inert immiscible organic phase in gas�liquid dispersions in stirred vessels influences the interfacial area in a more complex fashion than hitherto reported. As the organic phase fraction is increased, the interfacial area expressed on the basis of a unit volume of dispersion or aqueous phase, first increases, passes through a maximum and then decreases. This trend is observed irrespective of whether the area is determined by chemical means or by physical method. It is found that for low values of inert phase fraction, the average bubble size decreases whereas the gas holdup increases, resulting in increased interfacial area. The lower average bubble size is found to be due to partial prevention of coalescence as the bubbles size generated in the impeller region actually increases with the organic phase fraction. The actual values of interfacial areas depend on the nature of the organic phase. It is also found that the organic phase provides a parallel path for mass transfer to occur, when the solubility of gas in it is high.

Breakage of viscous and non-Newtonian drops in stirred dispersions

A model of breakage of drops in a stirred vessel has been proposed to account for the effect of rheology of the dispersed phase. The deformation of the drop is represented by a Voigt element. A realistic description of the role of interfacial tension is incorporated by treating it as a restoring force which passes through a maximum as the drop deforms and eventually reaching a zero value at the break point. It is considered that the drop will break when the strain of the drop has reached a value equal to its diameter. An expression for maximum stable drop diameter, dmax, is derived from the model and found to be applicable over a wide range of variables, as well as to data already existing in literature. The model could be naturally extended to predict observed values of dmax when the dispersed phase is a power law fluid or a Bingham plastic.

A multi-stage model for drop breakage in stirred vessels

Existing models for dmax predict that, in the limit of μd → ∞, dmax increases with 3/4 power of μd. Further, at low values of interfacial tension, dmax becomes independent of σ even at moderate values of μd. However, experiments contradict both the predictions show that dmax dependence on μd is much weaker, and that, even at very low values of σ,dmax does not become independent of it. A model is proposed to explain these results. The model assumes that a drop circulates in a stirred vessel along with the bulk fluid and repeatedly passes through a deformation zone followed by a relaxation zone. In the deformation zone, the turbulent inertial stress tends to deform the drop, while the viscous stress generated in the drop and the interfacial stress resist deformation. The relaxation zone is characterized by absence of turbulent stress and hence the drop tends to relax back to undeformed state. It is shown that a circulating drop, starting with some initial deformation, either reaches a steady state or breaks in one or several cycles. dmax is defined as the maximum size of a drop which, starting with an undeformed initial state for the first cycle, passes through deformation zone infinite number of times without breaking. The model predictions reduce to that of Lagisetty. (1986) for moderate values of μd and σ. The model successfully predicts the reduced dependence of dmax on μd at high values of μd as well as the dependence of dmax on σ at low values of σ. The data available in literature on dmax could be predicted to a greater accuracy by the model in comparison with existing models and correlations.

A multi-stage model for drop breakage in stirred vessels

Existing models for dmax predict that, in the limit of μd → ∞, dmax increases with 3/4 power of μd. Further, at low values of interfacial tension, dmax becomes independent of σ even at moderate values of μd. However, experiments contradict both the predictions show that dmax dependence on μd is much weaker, and that, even at very low values of σ,dmax does not become independent of it. A model is proposed to explain these results. The model assumes that a drop circulates in a stirred vessel along with the bulk fluid and repeatedly passes through a deformation zone followed by a relaxation zone. In the deformation zone, the turbulent inertial stress tends to deform the drop, while the viscous stress generated in the drop and the interfacial stress resist deformation. The relaxation zone is characterized by absence of turbulent stress and hence the drop tends to relax back to undeformed state. It is shown that a circulating drop, starting with some initial deformation, either reaches a steady state or breaks in one or several cycles. dmax is defined as the maximum size of a drop which, starting with an undeformed initial state for the first cycle, passes through deformation zone infinite number of times without breaking. The model predictions reduce to that of Lagisetty. (1986) for moderate values of μd and σ. The model successfully predicts the reduced dependence of dmax on μd at high values of μd as well as the dependence of dmax on σ at low values of σ. The data available in literature on dmax could be predicted to a greater accuracy by the model in comparison with existing models and correlations.

Breakage and coalescence of drops in turbulent stirred dispersions

The various existing models for predicting the maximum stable drop diameterd max in turbulent stirred dispersions have been reviewed. Variations in the basic framework dictated by additional complexities such as the presence of drag reducing agents in the continuous phase, or viscoelasticity of the dispersed phase have been outlined. Drop breakage in the presence of surfactants in the continuous phase has also been analysed. Finally, the various approaches to obtaining expressions for the breakage and coalescence frequencies, needed to solve the population balance equation for the number density function of the dispersed phase droplets, have been discussed.

A new model for the breakage frequency of drops in turbulent stirred dispersions

A model of drop breakage in turbulent stirred dispersions based on interaction of a drop with eddies of a length scale smaller than the drop diameter has been developed. It predicts that, unlike the equal breakage assumed by earlier models, a large drop reduces in size due to stripping of smaller segments off it through unequal breakage. It is only when the drop nears the value of the maximum stable drop diameter that it breaks into equal parts. This new model of drop breakage, coupled with the pattern of interaction of drops with eddies of different sizes existing in the vessel, has been used to evaluate not only the breakage frequency, but also the size distribution of the daughter droplets(which was hitherto assumed). The model has been incorporated in the population balance equation and the resulting cumulative size distributions compared with those availble in the literature.

A new model for coalescence efficiency of drops in stirred dispersions

A model for coalescence efficiency of two drops embedded in an eddy has been developed. Unlike the other models which consider only head-on collisions, the model considers the droplets to approach at an arbitrary angle. The drop pair is permitted to undergo rotation while they approach each other. For coalescence to occur, the drops are assumed to approach each other under a squeezing force acting over the life time of eddy but which can vary with time depending upon the angle of approach. The model accounts for the deformation of tip regions of the approaching drops and, describes the rupture of the intervening film, based on stability considerations while film drainage is continuing under the combined influence of the hydrodynamic and van der Waals forces. The coalescence efficiency is defined as the ratio of the range of angles resulting in coalescence to the total range of all possible approach angles. The model not only reconciles the contradictory predictions made by the earlier models based on similar framework but also brings out the important role of dispersed-phase viscosity. It further predicts that the dispersions involving pure phases can be stabilized at high rps values. Apart from explaining the hitherto unexplained experimental data of Konno et al. qualitatively, the model also offers an alternate explanation for the interesting observations of Shinnar.

Prediction of reaeration rates in square, stirred tanks

Experiments were conducted on the oxygen transfer coefficient, k(L)a(20), through surface aeration in geometrically similar square tanks, with a rotor of diameter D fitted with six flat blades. An optimal geometric similarity of various linear dimensions, which produced maximum k(L)a(20) for any rotational speed of rotor N by an earlier study, was maintained. A simulation equation uniquely correlating k = k(L)a(20)(nu/g(2))(1/3) (nu and g are kinematic viscosity of water and gravitational constant, respectively), and a parameter governing the theoretical power per unit volume, X = (ND2)-D-3/(g(4/3)nu(1/3)), is developed. Such a simulation equation can be used to predict maximum k for any N in any size of such geometrically similar square tanks. An example illustrating the application of results is presented. Also, it has been established that neither the Reynolds criterion nor the Froude criterion is singularly valid to simulate either k or K = k(L)a(20)/N, simultaneously in all the sizes of tanks, even through they are geometrically similar. Occurrence of "scale effects" due to the Reynolds and the Froude laws of similitude on both k and K are also evaluated.

Studies on transport phenomena during solidification of an aluminum alloy in the presence of linear electromagnetic stirring

Numerical and experimental studies on transport phenomena during solidification of an aluminum alloy in the presence of linear electromagnetic stirring are performed. The alloy is electromagnetically stirred to produce semisolid slurry in a cylindrical graphite mould placed in the annulus of a linear electromagnetic stirrer. The mould is cooled at the bottom, such that solidification progresses from the bottom to the top of the cylindrical mould. A numerical model is developed for simulating the transport phenomena associated with the solidification process using a set of single-phase governing equations of mass. momentum, energy. and species conservation. The viscosity variation of the slurry, used in the model, is determined experimentally using a rotary viscometer. The set of governing equations is solved using a pressure-based finite volume technique, along with an enthalpy based phase change algorithm. The numerical study involves prediction of temperature, velocity, species and solid fraction distribution in the mould. Corresponding solidification experiments are performed, with time-temperature history recorded at key locations. The microstructures at various temperature measurement locations in the solidified billet are analyzed. The numerical predictions of temperature variations are in good agreement with experiments, and the predicted flow field evolution correlates well with the microstructures observed at various locations.