Collective dynamics of flowing colloids during pore clogging

Separation and Purification Technology, 2012

Prediction of pore fouling by microparticles is still challenging and remains a difficult step to optimize membrane and filtration processes. The scientific issue consists in determining the relevant operation parameters controlling the capture of particles and the clogging of the filter. In this study, we have developed for a dead-end and cross-flow filtration a poly-dimethylsiloxane (PDMS) device which allows direct observation of the clogging dynamics of microchannels (20 lm wide) by micrometric particles (5 lm diameter). The experiments highlight the formation of different 3D clogging patterns according to the filtration conditions (particle concentration, flowrate, particle flux density and physical-chemical conditions of suspension). Besides, we have determined under which specific conditions of filtration, the latex microparticles are captured and form arches, clusters or dendrites. For each type of structure, the temporal dynamics of the particle deposition are analyzed by means of the average thickness of deposit. The critical conditions for the formation of arches leading to deposit formation have been identified in term of a combination of operating conditions: the particle concentration and the particle velocity. A critical particle flux density yielding pore clogging is then observed and characterized. Studying these experimental results helps to identify pore clogging mechanisms: deposition, interception and bridging.

Journal of Colloid and Interface Science, 2000

The interaction between stable colloidal particles arriving at a pore entrance was studied using a numerical method for the case where the particle size is smaller than but of the same order as the pore size. The numerical method was adapted from a front-tracking technique developed for studying incompressible, multifluid flow by S. O. Unverdi and G. Tryggvason (J. Comp. Phys. 100, 25, 1992). The method is based on the finite difference solution of Navier-Stokes equation on a stationary, structured, Cartesian grid and the explicit representation of the particle-liquid interface using an unstructured grid that moves through the stationary grid. The simulations are in two dimensions, considering both deformable and nondeformable particles, and include interparticle colloidal interactions. The interparticle and particle-pore hydrodynamic interactions, which are very difficult to determine using existing analytical and semi-numerical, semi-analytical techniques in microhydrodynamics, are naturally accounted for in our numerical method and need not be explicity determined. Two-and three-particle motion toward a pore has been considered in our simulations. The simulations demonstrate how the competition between hydrodynamic forces and colloidal forces acting on particles dictate their flow behavior near the pore entrance. The predicted dependence of the particle flow behavior on the flow velocity and the ratio of pore size to particle size are qualitatively consistent with the experimental observations of V. Ramachandran and H. S. Fogler (J.

Physical Review Fluids, 2018

Physical Review E, 2014

We present numerical results of the effect that the driving force has on the clogging probability of inert particles passing through a bottleneck. When the driving force is increased by four orders of magnitude, the mean avalanche size remains almost unaltered (increases 1.6 times) while the flow rate and the avalanche duration display strong dependence on this magnitude. This indicates that in order to characterize the ability of a system to clog, the right variable to consider is the number of particles that pass through the outlet. The weak dependence of this magnitude on the driving force is explained in terms of the average kinetic energy of the flowing grains that has to be dissipated in order to get an arch stabilized.

Langmuir, 2007

Transport of colloidal particles in porous media is governed by the rate at which the colloids strike and stick to collector surfaces. Classic filtration theory has considered the influence of system hydrodynamics on determining the rate at which colloids strike collector surfaces, but has neglected the influence of hydrodynamic forces in the calculation of the collision efficiency. Computational simulations based on the sphere-in-cell model were conducted that considered the influence of hydrodynamic and Derjaguin-Landau-Verwey-Overbeek (DLVO) forces on colloid attachment to collectors of various shape and size. Our analysis indicated that hydrodynamic and DLVO forces and collector shape and size significantly influenced the colloid collision efficiency. Colloid attachment was only possible on regions of the collector where the torque from hydrodynamic shear acting on colloids adjacent to collector surfaces was less than the adhesive (DLVO) torque that resists detachment. The fraction of the collector surface area on which attachment was possible increased with solution ionic strength, collector size, and decreasing flow velocity. Simulations demonstrated that quantitative evaluation of colloid transport through porous media will require nontraditional approaches that account for hydrodynamic and DLVO forces as well as collector shape and size.

Transport in porous media, 2024

Colloidal transport and clogging in porous media is a phenomenon of critical importance in many branches of applied sciences and engineering. It involves multiple types of interactions that span from the sub-colloid scale (electrochemical interactions) up to the porescale (bridging), thus challenging the development of representative modelling. So far published simulation results of colloidal or particulate transport are based on either reduced set of forces or spatial dimensions. Here we present an approach enabling to overcome both computational and physical limitations posed by a problem of 3D colloidal transport in porous media. An adaptive octree mesh is introduced to a coupled CFD and DEM method while enabling tracking of individual colloids. Flow fields are calculated at a coarser scale throughout the domain, and at fine-scale around colloids. The approach accounts for all major interactions in such a system: elastic, electrostatic, and hydrodynamic forces acting between colloids, as well as colloids and the collector surface. The method is demonstrated for a single throat model made of four spherical segments, and the impact of clogging is reported in terms of the evolution of the critical path diameter for percolation and permeability. We identified four stages of clogging development depending on position and time of individual colloid entrapment, which in turn correlates to a cluster evolution and local transport.

Physical Review Letters, 2018

From observations of the progressive deposition of noncolloidal particles by geometrical exclusion effects inside a 3D model porous medium, we get a complete dynamic view of particle deposits over a full range of regimes from transport over a long distance to clogging and caking. We show that clogging essentially occurs in the form of an accumulation of elements in pore size clusters, which ultimately

Scientific reports, 2018

Blockage of pores by particles is found in many processes, including filtration and oil extraction. We present filtration experiments through a linear array of ten channels with one dimension which is sub-micron, through which a dilute dispersion of Brownian polystyrene spheres flows under the action of a fixed pressure drop. The growth rate of a clog formed by particles at a pore entrance systematically increases with the number of already saturated (entirely clogged) pores, indicating that there is an interaction or "cross-talk" between the pores. This observation is interpreted based on a phenomenological model, stating that a diffusive redistribution of particles occurs along the membrane, from clogged to free pores. This one-dimensional model could be extended to two-dimensional membranes.

Vadose Zone Journal, 2011

Values of the applied hydrodynamic torque (T applied) and the resisting adhesive torque (T adhesion) will determine whether a colloid will be immobilized (T applied £ T adhesion) or roll (T applied > T adhesion) on a solid water interface. Previous literature has demonstrated in 1-2 collector (grain) systems that the influence of T applied on colloid retention can be significant under unfavorable attachment conditions and that only a fraction of the solid surface may contribute to retention. However, many questions remain on how to obtain, analyze, and upscale information on the forces and torques that act on colloids near solid surfaces in porous media. To address some of these gaps in knowledge, high resolution pore-scale water flow simulations were conducted for sphere packs (25 spheres) over a range of Darcy velocities, grain sizes and distributions, and porosities. The spatial variability of T applied was calculated from this information, and successfully described using a lognormal cumulative density function (CDF). Linear interpolation and scaling techniques were subsequently used to predict the lognormal CDF of T applied for various colloid sizes, grain sizes and distributions, and water velocities. The lognormal CDF of T applied was then evaluated at select values of T adhesion (i.e, interaction energy) to quantify the fraction and locations on the solid surface that contributes to colloid retention (S f), and the theoretical maximum solid phase concentration of retained colloids (S max).

Soft Matter, 2013

In this paper, we used numerical simulations to investigate the flow properties of soft glassy materials. These systems display a mixed fluid-solid behavior whose theoretical description remains a challenging task. The molecular dynamic simulations exhibit non-local rheological behavior, in direct line with previous experimental results. The inverse viscosity of the material at a given point, denoted as fluidity, is not a local function of the local stress, but also depends on the state of the system in the neighborhood, with a spatial correlation length typically equal to a few particles. The fluidity is furthermore related directly to the velocity fluctuations and rate of plastic events in the form of a scaling function. Correlations are the signature of a cooperative process at the origin of the flow and of the non-local effects. We compare the obtained results with a scalar fluidity model and emphasize the similarities between the two approaches.