1. Introduction
In environmental management, the removal of particulate solids from liquid process effluents is of great importance. However, when the sizes of solid particles diminish and reach micron and submicron range, the particles tend to remain in suspension and cannot be removed by gravity settling. In order to achieve an acceptable solid-liquid separation at a reasonable cost, the particles need to be agglomerated by flocculation followed by phase separation (sedimentation, filtration, centrifugation etc.) [1]
Flocculation plays a key role in several natural processes in the aquatic and marine environment such as sediment erosion, transport, aggregation and deposition and it is often employed in several industrial purification and separation processes [2–8]. Aggregation of cohesive suspended particles by size enlargement to form larger flocs is typically encountered in natural systems such as rivers, reservoirs, lakes and estuarine and in engineered systems such as bioreactors and agitation vessels [9–13]. It is widely employed in both upstream and downstream solid-liquid separation processes where the value of the individual phases are enhanced by destabilization and aggregation of the charged particles in suspension using high molecular weight synthetic polymers [14–18].
Typical mean shear rates in these diverse natural and engineered systems range from 2.5/s in natural systems to 5000/s especially in high shear bioreactors and mixing chambers [19, 20]. A number of interfacial forces and interactions that play key roles in this process are well-understood in the light of classical and extended DLVO theory [21–25], and provides a basis for any theoretical analysis of particle dispersion and colloidal stability. Therefore, the aim of this review paper is to provide an insight into the role of colloidal physicochemical process and hydrodynamics interactions in particle aggregation process—floc aggregation, growth and stability, and summarize the recent contributions in this field to the understanding of turbulence phenomenon and the direction of future research.
2. Structure Formation Processes in Dispersed Systems
2.1. Colloidal Stability and Interfacial Forces
Most colloidal particles in suspension (0.001 – 10 μm) tend to remain suspended and settle very slowly due to the presence of surface charge arising from isomorphic substitution, chemical reaction at the interface or preferential adsorption of ions on the particle surface from the surrounding medium [26, 27]. The charge on the particle surface is surrounded by excess of oppositely charged ions (counterions) in solution in order to maintain electrical neutrality. The combined interactions of this system of opposite charged cloud is known as the electric double layer [28]. The developed electrostatic repulsive forces prevent particle aggregation in colloidal system and hence contribute to their stability in the dispersed system [29]. The diffuse layer and the Stern layer of the charge cloud surrounding a particle together constitute the electric double layer [30].
A quantitative description of the colloidal stability and aggregation process was facilitated by the assumptions of the DLVO theory which considered interactions between particles as additive [28]. The DLVO theory was primarily based on two types of interactions between two compact spherical particles—electrostatic (including Born repulsion) [31] repulsive force due to the electrical double layer and Van der Waals attractive force of the dispersion.
The combination of these two interactions defines a total energy barrier VT or the so-called Gibbs interaction energy as shown in the energy profile in Fig. 1 [32, 28, 33–35, 25] which, if overcome by the kinetic energy of the moving particles, will result in aggregation. This is expressed mathematically in Eqs. (1–2), where VA, VR, ζ refers to the energy barrier due to Van der Waals attraction and electrostatic repulsion and zeta potential respectively. In the case of submicron particles, the kinetic energy is derived primarily from Brownian diffusion and less from hydrodynamic (fluid shear) or gravitational forces (differential settling). However, for larger or coarser particles, hydrodynamic interactions such as turbulence or fluid shear and gravitational forces imparts much of the kinetic energy needed to overcome the flocculation barrier [36, 33, 37]. A critical shear rate is required to overcome the energy barrier in order to allow flocculation and floc formation, and this is dependent on the surface charge and the size of particles. In general, the higher the charge on the particles, and the smaller the particle size, the higher the shear rates that are required to bring about flocculation [21, 35, 38].
In addition to the interactions described above, there exist additional interaction energies and structural forces which may arise for instance from the perturbation or re-arrangement of water molecules near the interface [4, 34, 39]. Non-DLVO interactions such as hydrophobic, hydration, oscillatory, structural, steric and electrosteric and hydrodynamic forces, which were previously not considered in colloidal interactions are now widely recognized to play important roles in colloidal aggregation [4, 21, 40–42]. These forces may be attractive, repulsive or oscillatory in nature and may be more pronounced or stronger than the DLVO forces [4]. In some cases where the energy barrier for aggregation exists at subnanometer distances, significant discrepancies between theoretical (DLVO-based) and experimental measurements of flocculation rates have been reported [13, 36]. For instance, in the case of colloidal dispersion containing clay particles, the theory provides a general conceptual model for the gross interactions of particles but fails to explain certain aspects of the interactions [27]. However, in most practical applications of flocculation, the dominant forces are electrostatic repulsion and London-Van der Waals attraction forces [43].
2.2. Kinetics of Fine Particle Aggregation
Flocculation kinetics deal with time-dependent changes in the dispersions or suspensions and provides information on the flocculation rate, dispersion stability and particle interactions which depends on the number and efficiency of the collisions [44]. The collisions between particles has been suggested to occur at a rate that depends on the transport mechanisms of Brownian diffusion, fluid shear and differential settling and it is assumed to be a second-order rate or reaction process [45]. The aggregation of dilute systems (TS ≤ 1%) with an initial concentration of primary particles n1 and fast aggregation constant kf can be represented by the expression in Eq. (3) [31, 45].
Furthermore, several studies have shown that aggregates formed from different transport mechanisms exhibit marked variations in their structural attributes. The aggregates formed by shear-induced collision are known to be stronger than those formed from other transport mechanisms, those from Brownian motion tend to be easily dispersed by shearing, while those from differential settling results in ragged, weak, and low density aggregates [46].
Smoluchowski using the “population balance” approach, expressed time evolution of the number density of discrete particles of size k as they aggregate with respect to time in terms of the collision frequency or rate function β originally for Brownian diffusion and laminar or uniform shear flow, and later for differential settling assuming binary collisions between the particles [40, 45, 47–49]. Although the Smoluchowski equation has been modified in line with recent findings to include several other parameters, the general form of the expression is presented in Eq. (4). The quantity βij is the collision frequency or rate function for collisions between ith and jth sized particles while αij is the dimensionless collision efficiency factor. The collision frequency depends on the transport mechanisms of Brownian motion, fluid shear and differential settling, whereas the collision efficiency is a function of the degree of particle destabilization and it gives the probability of collision leading to attachment with values ranging from 0 to 1 [47, 50–53].
The first term of the right hand side of Eq. 4 describes the increase in number of particles of size k by flocculation of two particles whose total volume is equal to the volume of a particle of size k, while the second term on the right hand side describes the loss of particles of size k due to aggregation with particles of other sizes. The general form of the equation expresses the rate of change in the number concentration of particles of size k. In arriving at Eq. (4), Smoluchowski made a number of key assumptions listed below:
2.3. Flocculation Transport Processes
The interparticle collisions are promoted by transport mechanisms of Brownian motion or perikinetic (for particles with diameters less than 1 μm), fluid shear or orthokinetic (for particles in the diameter range 1–40 μm) and differential settling (typically for particles with diameter larger than 40 μm). The collision frequency or rate function β for these three mechanisms can be expressed mathematically in Eqs. (5–7) where βBM, βSH, βDS refers to the collision frequency function for Brownian motion, fluid shear, and differential settling respectively. The total collision rate or frequency βT is the sum of contributions from each of the three transport mechanisms (Eq. (8)) [3, 5, 10, 46, 48, 51, 54–56].
2.3.1. Perikinetic particle aggregation
The random displacement of particles in Brownian motion as a result of the thermal energy of the system is termed perikinetic [21, 40, 42, 45]. In practical applications, perikinetic aggregation is only important for very small particles where they are continuously bombarded by the surrounding water molecules ( < 1 μm) [57]. Under this prevailing condition, the rate of floc growth cannot be sufficiently sustained for an effective phase separation. Hence, the importance of other transport mechanisms in promoting the floc growth kinetics. However, in the case of contact between a small particle and micro flocs, Brownian diffusion can still control the transport of small particles across the layer of fluid on the floc surface [58].
2.3.2. Classical orthokinetic aggregation
The agitation of suspended particles in liquid medium will lead to aggregation due to the inter particle collision induced by the fluid motion [37]. The velocity gradient or shear rate is dependent on the nature of the fluid flow. In uniform laminar shear flow, the velocity gradient remains constant in the entire flow field while in turbulent flow, there is a rapid fluctuation of the velocity gradient. Consequently, the velocity gradient is a function of both space and time (spatial and temporal) [59, 60]. Camp and Stein in their critique of the Smoluchowski approach to shear flocculation modelling [61], introduced the concept of root-mean-square (R.M.S.) or absolute velocity gradient Ḡ to account for the variations in the shear rate. [26, 46, 58, 62, 63]. Considering the angular distortion of an elemental volume of fluid arising from the application of tangential surface forces, Ḡ is defined according to Eqs. (9–10) as the R.M.S. velocity gradient in the mixing vessel, where u, v, w are the components of the fluctuating velocity while x, y, z refers to the 3-D Cartesian coordinate system, P is the power dissipated, V is the volume of the mixing vessel, ε is the energy dissipation rate per unit mass, μ s the dynamic viscosity and v is the kinematic viscosity [26, 64].
In practical applications, owing to the difficulties associated with the calculation of absolute velocity gradient Ḡ due to the fluctuations in the energy dissipation within the mixing vessel, an average velocity gradient Gave has often been used in place of the absolute value Ḡ [59, 60] (Eq. (11)). Korpijarvi et. al. [65], in their study of mixing and flocculation in a jar using computational fluid dynamics (CFD) reported a large variation in the local velocity gradient GL within the mixing vessel. In similar studies conducted by Kramer and Clark [26, 66] as well as Mühle [67], a lower estimate of the absolute velocity gradient was proposed (Eqs. (12–13)) to account for the variation of the velocity gradient throughout the fluid where quantities GL, Pave, ε, V, μ, v represent the local velocity gradient at different points within the mixing vessel, average power consumption, volume of the mixing vessel, kinetic energy dissipation rate, dynamic and kinematic viscosity respectively.
2.3.1.1. Particles in laminar shear
According to the rectilinear model of particle trajectory, two particles moving in a fluid on different streamlines will experience a velocity gradient which indicates the relative motion of particles and the possibility of collision and aggregation [45, 68]. In simple laminar or uniform shear with a well-defined flow field such as in pipe flow or in some Couette devices [19, 69, 70], the transport of particles by laminar shear due to the fluid motion can be characterized by a single value of the shear rate G [26]. Under such condition, the particles in laminar shear will exhibit rotational motion in the direction of travel with a constant angular velocity ω [71, 72].
2.3.1.2. Particles in turbulent shear
In most practical applications of flocculation, the aggregation processes typically occur under turbulent conditions [26, 68, 73–75]. The particles in suspension under the influence of turbulent flow will experience fluctuating motion of the fluid with the particles being transported by the fluid eddies and the flow vortex [26, 58]. Thus, small particles suspended in fluid exist in an environment of small energy-dissipating eddies in most typical sheared reactors [19, 76]. Under such condition, particle collisions are promoted by eddy size similar to that of the colliding particles [47, 77]. The mechanism of flocculation in turbulent shear has been suggested to be similar to that of laminar shear for particles smaller than the Kolmogorov microscale λ0 of turbulence, while for particle larger than this length scale, the flocculation mechanism is similar to that of Brownian diffusion [21, 26].
The formation of aggregates in any turbulence or shear-induced flocculation essentially consist of the following phases: destabilization; collision and adhesion; floc growth and deformation phases [78]. Ives [69] and Bergenstahl [79] expressed the kinetics of orthokinetic aggregation in terms of the change in the initial particle concentration, while Levich [62, 80] in an earlier study derived the collision rate equation based on the concept of locally isotropic turbulence assuming a viscous diffusive subrange (i.e. particle size smaller than Kolmogorov micro scale) for shear-induced turbulent flow where N, r, d represent the particle number concentration per unit volume, particle radius and diameter respectively (Eqs. (14–15)). Several other models of flocculation kinetics are available in literature [27, 81]. In most of these models, the shear rate G and the kinetic energy dissipation rate ɛ are shown as important drivers of the flocculation process. Typically, orthokinetic flocculation can be experimentally observed in a conventional stirred tank consisting of an axially mounted impeller and a circular or rectangular mixing vessel.
2.3.2. Extended orthokinetic aggregation
Pelleting flocculation is considered as an extension of the classical orthokinetic transport mechanism based on the “metastable state” concept [82, 83]. It has been shown by several studies that the efficiency of the floc formation process as well as the floc structural attributes (size, shape, porosity, density, etc.) can be significantly improved by the application of suitable mechanical energy [84–94]. This is accomplished by the application of a non-destructive uneven force system on the floc surface, by the action of the turbulent fluid motion or the so-called floccule mechanical synaeresis induced by rolling and collision mechanisms [95–99].
In the rolling mechanism, as the floc rolls along a plane surface, it is subjected to a pressure fluctuation which becomes greater at the front of the central point in the contact area and lesser at the rear. On the other hand, in the collision mechanisms, a floc experiences an impact due to the floc-floc or floc-surface collision. In order to withstand this impact, the floc must adhere strongly, in order words, the impact results in a more compact floc [83]. The resulting pressure induced in a curvilinear trajectory on the floc surface can be expressed as a function of the velocity gradient G where m, R, G and A represent the mass of the floc, radius of the flow trajectory, velocity gradient and the area of the floc projection onto a plane oriented normally to the vector of the resultant force respectively (Eq. (16)). In practice, pelleting flocculation is realized through the selection of an appropriate process engineering conditions and stirrer-vessel system (geometry and configuration) in order to obtain the necessary conditions suitable for rolling and collision-mediated floc pelletization. [83, 91].
The structure formation by pelleting flocculation yields “pellet flocs” with superior structural attributes when compared to the “random flocs” produced by the classical flocculation process which tend to be loose and bulky [84]. Previous studies have reported two distinct models of pelleting flocculation, namely parallel and series models resulting in onion and raspberry-like pellet floc structures as presented schematically in Fig. 2 [87, 90, 91, 100–103]. In the series model, orthokinetic flocculation occurs first where random, loose and bulky flocs are produced. These voluminous flocs shrink and densify by the application of suitable mechanical energy (mechanical synaeresis).
Conversely, in the parallel model, there is no time delay between the orthokinetic flocculation and mechanical synaeresis. The dispersed particles in suspension and micro flocs are transported by the fluid motion to the surface of the “mother seed” where they are attached by polymer bridging [95, 100, 102]. In terms of the flocculation efficiency and process relevant (phase separation) parameters, pellet flocs can be characterized by their size, shape, density, porosity, compressive strength and deformability.
2.3.3. Polymer-mediated interactions
The energy barrier to flocculation can be overcome in several other ways apart from the kinetic energy derived from the transport mechanisms. The electrostatic charge can be lowered with the aid of polymeric flocculants by double layer compression, charge neutralization, and polymer bridging between particles [32, 33, 104]. In practice, both of these approaches are employed either concurrently or sequentially. In addition to these interactions, several other mechanisms have also been attributed to the polymer flocculation of charged particles in suspension. Electrostatic charge-patch; depletion flocculation; network flocculation and polymer complex formation are few of the other relevant interactions [2, 4, 36, 38, 73, 105–108].
However, the dominant concepts in the case of poly-electrolyte-mediated destabilization are surface charge neutralization (ion exchange), charge patch formation and polymer bridging as illustrated in Fig. 3 [44, 109–115]. The important physicochemical properties that influence the flocculation process in the case of polymer-mediated interactions include, polymer molecular weight, charge density and concentration, zeta potential or surface charge of suspension, particle size and distribution, specific surface area, solution conductivity and pH [105].
In the case of pelleting flocculation, Higashitani and Kubota [98] in their study of the pelleting flocculation of polystyrene latex particles reported that the formation of pellet flocs depends on the suspension concentration, flocculant molecular weight and charge density as well as the intensity of mixing. They argued that pellet flocs are formed within a limited range of these parameters hence pelleting flocculation requires a high degree of control of different process conditions. On the basis of their findings, they expressed mathematically the volume fraction of particles Øpt at the critical particle concentration Nc above which pellet flocs are formed as a function of the initial particle concentration N0, polymer concentration Cp, charge density p, and the molecular weight M, where NA is the Avogadro constant (Eq. (17)).
3. Hydrodynamics of Fine Particle Aggregation
Hydrodynamically-induced turbulent shear is an important driver of the flocculation process especially in the case of orthokinetic aggregation of particles [33, 36, 116]. The floc growth and stability in any flocculated system has been suggested to be a function fluid-particle interactions and the intrinsic physicochemical properties of the floc [75, 77, 117]. The dynamics of these interactions affects all facets of the flocculation process and the degree of aggregation in any sheared system [118]. In the case of hydrodynamic interactions, induced velocity gradient promote the aggregation process but might also be responsible for floc breakage as a result of increased viscous shear stress [119]. Consequently, in the case of shear-induced collisions, the hydrodynamic effect can be very significant [57].
The chemical transport in natural and engineered systems is often controlled more by the fluid properties, in this case the mixing rate of the fluid, than by chemical properties such as the molecular diffusivity of the compound. Therefore, in the design of flocculation reactors, it is desirable to have uniform mixing intensities throughout the reactor [26]. Under the condition of shear-induced turbulence due to the use of flow inducers such as paddles, mixers, stirrers or bubbles (Fig. 4), the convective currents in the system are insignificant as compared to the advective currents generated in the system [19]. The numerical modelling of the phenomena of fluid-particle mixing, floc formation, growth and breakup is very complex, as such minimizing the energy dissipation rate (degree of turbulent mixing) in the reactor will lead to improvement in flocculation efficiency [26].
3.1. Fluid Mixing and Particle Dispersion
Turbulent mixing characterized by high Reynolds number (Re > 104) [120], can be viewed as a hierarchy of irregular, rotational and dissipative motion containing vorticity (curl or rotation of the velocity vector) on different scales or eddy sizes [25, 121]. In turbulent flows such as those encountered in pipe or tube flow, channel flow and stirred tanks or mixing chambers, energy transfer occurs on different eddy scales [64, 67]. Eddies are spatial recognizable flow patterns that exist in a turbulent flow for at least a short time in which there is a correlation between the velocities at two different points [26]. Fluid deformation causes vortices to stretch and vorticity and kinetic energy to be transported from larger to smaller eddies [25, 121]. The turbulent vortex in such fluid motion is generally propagated in the tangential and axial direction than in the radial direction with the vorticity increasing with decreasing eddy size [64, 121, 122].
In most practical applications there exist three scales of mixing namely: macro, meso and micro mixing. Macromixing refers to mixing that is driven by the largest scale of motion in the fluid Λ (integral length scale), meso mixing on the other hand involves mixing on a scale smaller than the bulk circulation but larger than the micro mixing, while micro mixing is the mixing on the smallest scale of fluid motion λ0 (Kolmogorov microscale) expressed in Eq. (18) and at the final scales of molecular diffusivity (Batchelor scale) [120, 123–126]. Micromixing describes the process homogenization of liquid balls with their surroundings on a molecular level. In typical mixing conditions, the dividing line between micro and macro scale is between 100 and 1000 μm respectively [127].
The largest eddies in turbulent dispersion that represent the macro scale of turbulence or the integral length scale (Λ) contains most of the energy, and are produced by the stirrer or agitator head. The size of this macroscale of turbulence is of the order of the diameter of the agitator [19, 53, 120]. The turbulent flow can be viewed as an eddy continuum, with their size ranging from the dimension of the turbulence generating device to the Kolmogorov length scale [53]. In between the energy-containing eddies of the integral length scale Λ at the upper end of the inertia subrange and the smallest eddies of Kolmogorov microscale λ0 at the beginning of viscous dissipation range [128], there exist many eddies of other scales smaller than the integral scale Λ that transfer kinetic energy continually through the other length scales. The Batchelor λB and Taylor λT scales expressed in Eqs. (19–20) are the examples of other important length scales where D is the diameter of the flocs while ε is the turbulent kinetic energy dissipation rate.
The Taylor scale λT is an intermediate length scale in the viscous subrange that is representative of the energy transfer from large to small scales, but it is not a dissipation scale and does not represent any distinct group of eddies [129]. Batchelor micro-scale on the other hand is a limiting length scale where the rate of molecular diffusion is equal to the rate of dissipation of turbulent kinetic energy and represent the size of the region within which a molecule moves due to diffusional forces [123, 125, 128, 130].
3.2. Energy Dissipation in Turbulent Flow
The power input into any agitated system will eventually dissipates as heat due to viscous forces [76]. The larger eddies of the order of the agitator diameter (Λ~dA) draw energy from the fluid motion, while the smaller eddies transfer that energy gradually and continually to the smallest eddy, where the energy is ultimately dissipated into heat by friction [27]. The turbulence parameters of interest in any agitated micro-system (with respect to the mixer capacity and performance) are the Kolmogorov micro scale λ0, turbulence kinetic energy dissipation rate ɛ, and agitator tip velocity Vtip[122, 124, 131].
The Kolmogorov length scale λ0 is an indicator of the rate of micro mixing and mode of agglomerate deformation due to turbulent shear while the energy dissipation rate ɛ is the rate of the dissipation of kinetic energy. The tip velocity Vtip is a measure of the tangential velocity imparted by the flow inducer, an indicator of the strength of the vortex generated by the flow. Peripheral velocity Vpkol on the other hand gives an indication of the velocity of the flow vortex at Kolmogorov micro-scale [76, 122, 124, 131]. In typical laboratory mixing experiments (D < 1m), the micro-scale of turbulence predominates. Under such conditions, particle collision is promoted by eddy size similar to those of the colliding particles [47].
4. Floc Stability in Sheared Systems
Floc stability under the influence of hydrodynamic force has been suggested to be a function of floc binding or cohesive force FB, and hydrodynamic breaking force FH [64, 81, 116, 132–135]. While the binding force is determined by the flocs’ structure and physicochemical attributes, flow turbulence is the principal factor in the case of hydrodynamic force [59, 75, 117]. Of these two governing factors of floc stability, turbulence is the least understood owing to its complex nature [45, 136]. Therefore, a detailed analysis of floc stability under turbulent conditions often encountered in natural and agitated systems is not only difficult but often time consuming [45].
Recent advances in computational fluid dynamics has vastly improved the understanding of turbulence phenomenon. However, numerical modelling of complex turbulent fluid-particle interactions remains a challenge in fluid mechanics and it is computationally intensive [122]. In addition, the assumption of the rectilinear motion of particles prior to collision in flocculation modelling has been supplanted by recent findings of curvi-Environmental Engineering Research 21(1) 1–14 linear particle trajectory prior to collisions [47, 63]. Particles have been shown to move around other particles in a curvilinear path thereby reducing the probability of collision and attachment [47, 57].
4.1. Aggregate Strength and Hydrodynamic Stress
Floc formation and growth process in turbulent flocculation comprises of the lag phase, swift growth phase, and steady state phase or breakup/restructuring phase [135, 137, 138]. The overall rate of floc growth at the initial formation phase up to the steady state phase has been suggested to be a balance between collision-induced particle aggregation and floc cohesive strength on one hand, and the rate of floc breakage due to hydrodynamic stress on the other hand [64, 81, 116, 132–134, 139, 140].
The aggregate cohesive strength τ is a function of the physicochemical conditions and floc properties while the turbulent hydrodynamic stress σ depends on the design of the aggregation unit (geometry) and the mixing intensity [27, 75, 81, 117, 132, 134, 141, 142]. A number of empirical models have been proposed for predicting the maximum hydrodynamic breaking force FH and the global hydrodynamic stress σ acting on a spherical aggregate in the inertia and viscous domain of turbulence (Table. 1).
He et. al. [135], in a study of floc strength under turbulent flow conditions derived a mathematical expression for an estimate of the floc binding or cohesive force FB using fractal dimension approach while Lu et al. [35] in a similar study presented a theoretical model for the aggregate binding force FB of spherical mono-disperse particles both of which are presented in Eq. (21–22). Yuan & Farnood [143], as well as Jarvis et. al. [144], in their review of aggregate strength and breakage gave an estimate of the floc strength τ and global hydrodynamic stress σ obtained from empirical studies of various types of flocs.
The global hydrodynamic stress σ due to the shearing action of the fluid motion on the floc as well as the overall mechanical strength of an aggregate τ assuming a uniform floc shape and constant porosity can be expressed mathematically in Eqs. (23–24) [77, 116]. In turbulent hydrodynamic flow, the fluctuating motions of the fluid are responsible for shear, tensile, and compressive stress (normal stresses) on the flocs, depending on the direction of the fluctuating velocities acting on the aggregates [67, 77, 147, 148]. The exact type and magnitude of the prevailing stress causing aggregate disruption depends mainly on the relation between the floc size and the eddy size or radius, rw [67, 149].
4.2. Fluid-Particle Interactions and Floc Stability
The size of aggregates varies from molecular dimensions to a range that is visible to the unaided eye, with the smaller sizes being associated with the primary particles of diameter dp, while the largest size dFmax is determined by the balance of floc growth and rupture within the fluid [64, 139, 150, 151]. The floc growth and breakage is known to occur simultaneously—growing flocs are subjected to breakage while fragments of broken flocs undergo re-agglomeration until a levelling off of the floc sizes at steady state when the maximum stable size dFmax is attained [54, 139, 152, 153].
There is an increase in floc size as long as the hydrodynamic force due to turbulent shear is less than the floc cohesive or binding force, and after an extended period of time, an equilibrium floc size distribution is reached. In this case either a continued particle or micro flocs attachment to the larger flocs is prevented, or floc breakup kinetics balances the turbulence-induced collision [27, 67]. Increasing the shear rate (agitation or mixing speed) or prolonged shearing beyond the steady state leads to a sharp reduction in the mean floc size [10, 27, 35, 139, 154, 155]. The steady-state phase is regarded as the balance between floc growth and breakage under a given shear condition. In the viscous subrange, Lu et. al. [35], expressed this phenomenon in terms of the kinetic equation of flocculation as presented in Eq. (25). The quantities NA, NB, σS, τS, represent change in the particle number concentration per unit volume for aggregation and breakage, the aggregate shear strength and the shearing stress respectively.
In the analysis of floc breakup under turbulent hydrodynamic conditions, three approaches are normally employed namely: the limiting strength approach, maximum strain rate approach and the maximum floc size approach [26]. The first approach to floc breakup analysis which is based on force balances [67], requires an accurate description of the flow field in order describe the stress acting on the model floc. In this approach for analyzing the breakage phenomenon, under a given agitation condition, the critical condition of floc breakage is attained when the hydrodynamic breaking force FH is greater than or equal to the floc binding or cohesive force FB (B ≤ 1).
Consequently, the conceptual form of the floc growth and breakage rate mechanism can be expressed mathematically in Eqs. (26–27) while Eqs. (28–29) expressed the critical condition of floc breakage for the inertia and viscous subrange respectively [59, 59, 60, 64, 81, 116, 132–135, 140, 147]. A number of empirical models have been developed for estimating the maximum aggregate size dFmax. Few of such expressions are presented for the inertia and viscous subrange of turbulence in Table 2 [143, 144].
The fluctuating instantaneous fluid velocities acting parallel to the surface of the floc will induce a local shear stress σs on the aggregate in the viscous subrange reaching a maximum value when the floc size is roughly equal to the eddy scale (Kolmogorov micro scale) [118]. Similarly, fluctuating fluid velocities normal to the floc surface or dynamic pressure fluctuations acting on opposite sides of an aggregate will results in normal or bulk pressure stress σt (tensile or compressive) [53]. In addition, turbulent drag forces FD acting on the surface of an aggregate which originates from the local motion of the fluid relative to the motion of the aggregates will results in instantaneous surface shear forces and shear stress respectively [77]. It has been shown that for flocs in the inertial subrange of turbulence (dF > λ), tensile stress will predominate causing wholesale fracture or fragmentation [35], while for those in the viscous subrange (dF < λ), shear stress will cause erosion of the particles, floc shell or floc surface as shown schematically in Fig. 5 [47, 53].
4.3. Mechanisms of Aggregate Disruption
Several mechanisms have been proposed to account for the disruption of aggregates in orthokinetic flocculation [26–28, 46, 47, 57, 67, 117, 119, 144, 146, 156–158]. In addition to floc splitting, fracture or bulk rupture and surface erosion mentioned earlier, several other mechanisms have been proposed to be responsible for floc breakage but with little experimental data or theoretical analysis. Such mechanisms include aggregate-aggregate collisions, aggregate-stirrer collisions, as well as collisions with baffles and the tank wall [35, 138].
A number of theoretical and empirical models have been developed to account for the hydrodynamic stress exerted on particle agglomerates and their cohesive strength. Considering a model floc undergoing rupture by fragmentation (σt ≥ τt) in the inertia domain of turbulence (dF > λ), subjected to the shearing action of fluid motion, the turbulent hydrodynamic tensile stress σt causing bulk rupture or fragmentation by floc splitting and the aggregate tensile strength τt resisting the fluctuating pressure on the floc may be expressed mathematically in Eqs. (30–31).
Similarly, for a floc rupture due to the erosion of primary particles, micro flocs, or floc shell (σs ≥ τs) in the viscous subrange (dF < λ); the turbulent hydrodynamic shear stress σs eroding the surface of the aggregates and the corresponding aggregate shear strength τs resisting the viscous shear forces can be expressed mathematically in Eqs. (32–33) [35].
In a related study, Mühle [67] presented floc breakup analysis on the basis of rheology. The strength of a model floc was described in terms of the surface shear yield strength τy resisting floc shell erosion due to pseudo-surface tension force γF eroding the particle chain of the outer floc surface in the viscous domain of turbulence (σs ≥ τy) (Eq. (34)). The yield stress approach was also presented by Liu et. al. [147] for calculating the maximum floc tensile yield stress σy at which breakage is likely to occur in the inertia subrange (σt ≥ σy) (Eq. (35)).
On a similar basis, Attia [118] presented models for predicting the critical fluid velocity Vc above which there will be floc breakage by estimating floc yield stress σy (tensile or compressive) resulting from dynamic pressure acting on the floc (Eqs. (36–37)). The kinetic equations for fragmentation and surface erosion mechanisms are presented in Eqs. 38–39 based on the empirical results from flocculation in stirred reactor. The quantities NA, NB, σs, τs, represent change in the particle number concentration per unit volume for aggregation and breakage, the aggregate shear strength and the shearing stress respectively [35, 49].
4. Practical Applications and Future Outlook
Over the past decades, significant progress has been made in understanding the complex turbulent interactions that influence and drive the flocculation process. Improvement in modelling and simulation tools has led to the use of computational techniques in turbulent flow investigations. Nowadays, computational fluid dynamics provide a comprehensive information for process engineers necessary to design, optimize or retrofit various unit processes involving fluid flow as well as improving our understanding of flows in natural systems. This is extremely useful especially in investigations involving turbulent flows in devices with complex geometries with several design parameters requiring optimization without any need to perform the physical experiments [64, 136].
However, the numerical study of turbulent flow and flocculation process is still being done independently and only a few studies [159–164] deal with the integration of these interdependent processes. This is partly driven by the evolution pattern of computational fluid modelling (CFD) which was specifically developed for other industries but is now widely adopted across many fields [64]. Therefore, future studies must strive to advance research in this direction by coupling the modelling of turbulent fluid flow and aggregation processes within the computational environment in order to build on the recent progress made in these fields [64]. This in turn will assist process engineers in developing new generation of phase separation reactors more efficiently and improve our understanding and ability to accurately predict flow conditions in natural systems that are very vital for many environmental processes.
5. Conclusions
Flocculation process has been shown to be an effective particle phase separation process in engineered unit processes and is an important driver of many environmental processes (sediment transport, contaminant migration etc.) in the natural systems. The importance of hydrodynamics and flow turbulence as the main drivers of this process has been highlighted. Although significant progress has been made in the understanding of physicochemical aspects of flocculation, however, turbulence as remains a poorly understood phenomenon despite its profound influence on many engineered and life processes. This study highlights the need for closer integration of modelling and simulation tools in order to improve our understanding and ability to influence and predict these processes much more accurately.