Yannick Hallez

Effects of channel geometry on buoyancy-driven mixing

The evolution of the concentration and flow fields resulting from the gravitational mixing of two interpenetrating miscible fluids placed in a tilted tube or channel is studied by using direct numerical simulation. Three-dimensional (3D) geometries, including a cylindrical tube and a square channel, are considered as well as a purely two-dimensional (2D) channel. Striking differences between the 2D and 3D geometries are observed during the long-time evolution of the flow. We show that these differences are due to those existing between the 2D and 3D dynamics of the vorticity field. More precisely, in two dimensions, the strong coherence and long persistence of vortices enable them to periodically cut the channels of pure fluid that feed the front. In contrast, in 3D geometries, the weaker coherence of the vortical motions makes the segregational effect due to the transverse component of buoyancy strong enough to preserve a fluid channel near the front of each current. This results in three different regim...

Interaction between two spherical bubbles rising in a viscous liquid

The three-dimensional flow around two spherical bubbles moving in a viscous fluid is nstudied numerically by solving the full Navier-Stokes equations. The study considers the ninteraction between two bubbles for moderate Reynolds numbers (50 ≤ Re ≤ 500, Re nbeing based on the bubble diameter) and for positions described by the separation S n(2.5 ≤ S ≤ 10, S being the distance between the bubble centres normalized by the nbubble radius) and the angle θ (0o ≤ θ ≤ 90o ) formed between the line of centre and the ndirection perpendicular to the direction of the motion. We provide a general description nof the interaction extending the results obtained for two bubbles moving side by side n(θ = 0o ) by Legendre, Magnaudet & Mougin 2003 (J. Fluid Mech., 497,133-166) and nfor two bubbles moving in line (θ = 90o ) by Yuan & Prosperetti 1994 (J. Fluid Mech., n278, 325-349). Simple models based on physical arguments are given for the drag and lift nforces experienced by each bubble. The interaction is the combination of three effects: na potential effect, a viscous correction (Moore correction) and a significant wake effect nobserved on both the drag and the transverse force of the second bubble when located nin the wake of the first one.

A numerical investigation of horizontal viscous gravity currents

We study numerically the viscous phase of horizontal gravity currents of immiscible fluids in the lock-exchange configuration. A numerical technique capable of dealing with stiff density gradients is used, allowing us to mimic high-Schmidt-number situations similar to those encountered in most laboratory experiments. Plane two-dimensional computations with no-slip boundary conditions are run so as to compare numerical predictions with the ‘short reservoir’ solution of Huppert (J. Fluid Mech., vol. 121, 1982, pp. 43–58), which predicts the front position lf to evolve as t1/5, and the ‘long reservoir’ solution of Gratton & Minotti (J. Fluid Mech., vol. 210, 1990, pp. 155–182) which predicts a diffusive evolution of the distance travelled by the front xf ~ t1/2. In line with dimensional arguments, computations indicate that the self-similar power law governing the front position is selected by the flow Reynolds number and the initial volume of the released heavy fluid. We derive and validate a criterion predicting which type of viscous regime immediately succeeds the slumping phase. The computations also reveal that, under certain conditions, two different viscous regimes may appear successively during the life of a given current. Effects of sidewalls are examined through three-dimensional computations and are found to affect the transition time between the slumping phase and the viscous regime. In the various situations we consider, we make use of a force balance to estimate the transition time at which the viscous regime sets in and show that the corresponding prediction compares well with the computational results.

Turbulence-induced secondary motion in a buoyancy-driven flow in a circular pipe

We analyze the results of a direct numerical simulation of the turbulent buoyancy-driven flow that sets in after two miscible fluids of slightly different densities have been initially superimposed in an unstable configuration in an inclined circular pipe closed at both ends. In the central region located midway between the end walls, where the flow is fully developed, the resulting mean flow is found to exhibit nonzero secondary velocity components in the tube cross section. We present a detailed analysis of the generation mechanism of this secondary flow which turns out to be due to the combined effect of the lateral wall and the shear-induced anisotropy between the transverse components of the turbulent velocity fluctuations.

Experimental and numerical investigations of flow structure and momentum transport in a turbulent buoyancy-driven flow inside a tilted tube

Buoyancy-driven turbulent mixing of fluids of slightly different densities [At=Δρ/(2⟨ρ⟩)=1.15×10−2] in a long circular tube tilted at an angle θ=15° from the vertical is studied at the local scale, both experimentally from particle image velocimetry and laser induced fluorescence measurements in the vertical diametrical plane and numerically throughout the tube using direct numerical simulation. In a given cross section of the tube, the axial mean velocity and the mean concentration both vary linearly with the crosswise distance z from the tube axis in the central 70% of the diameter. A small crosswise velocity component is detected in the measurement plane and is found to result from a four-cell mean secondary flow associated with a nonzero streamwise component of the vorticity. In the central region of the tube cross section, the intensities of the three turbulent velocity fluctuations are found to be strongly different, that of the streamwise fluctuation being more than twice larger than that of the sp...

Quantitative assessment of the accuracy of the Poisson-Boltzmann cell model for salty suspensions.

The cell model is a ubiquitous, fast, and relatively easily implemented model used to estimate the osmotic pressure of a colloidal dispersion. It has been shown to yield accurate approximations of the pressure in dispersions with a low salt content. It is generally accepted that it performs well when long-ranged interactions are involved and the structure of the dispersion is solidlike. The aim of the present work is to determine quantitatively the error committed by assuming the pressure computed with the cell model is the real osmotic pressure of a dispersion. To this end, cell model pressures are compared to a correct estimation of the actual pressures obtained from Poisson-Boltzmann Brownian dynamics simulations including many-body electrostatics and the thermal motion of the colloids. The comparison is performed for various colloidal sizes and charges, salt contents, and volume fractions. It is demonstrated that the accuracy of the cell model predictions is a function of only the average intercolloid distance scaled by Debyes length κd̅ and the normalized colloidal charge. The cell model is accurate for κd̅ < 1 and not reliable for κd̅ > 5 independently of the colloidal charge. In the 1 < κd̅ < 5 range, covering a wide set of experimental conditions, the colloidal surface charge has a large influence on the error associated with the cell approximation. The results presented in this article should provide a useful reference to determine a priori if the cell model can be expected to predict accurately an equation of state for a given set of physicochemical parameters.

A Three-Step Scenario Involved in Particle Capture on a Pore Edge

A scenario is proposed to describe the capture of a spherical particle around a cylindrical pore. This geometry, ideal as far as the problem of particle capture on a filtration membrane is concerned, is clearly relevant in view of the pore-scale geometry of nucleopore or microsieve filtration membranes, and also of some microfluidic systems used to perform fluid-particle separation. The present scenario consists of three successive steps: particle deposition on the membrane away from the pore, subsequent reentrainment of some of the deposited particles by rolling on the membrane surface, and final arrest by a stabilizing van der Waals torque when the particle rolls over the pore edge. A modeling of these three steps requires the hydrodynamic and physicochemical particle-membrane interactions to be detailed close to the singular pore edge region and raises questions concerning the role of particle surface roughness. The relevance and robustness of such of a scenario for rough micrometer-sized latex particles is emphasized and comparisons are made with existing experimental data.

The continuous modeling of charge-stabilized colloidal suspensions in shear flows

Flows of concentrated colloidal suspensions may exhibit a rich set of behaviors due to both hydrodynamic and colloidal interactions between the particles. Colloidal flows are generally modeled with an effective Navier–Stokes equation and a mass balance for the solid phase involving a diffusion coefficient given by the generalized Stokes–Einstein relation. This picture corresponds to a near equilibrium regime in which entropic and colloidal forces dominate over hydrodynamic interactions, the latter being totally ignored. On the other hand, suspension flows far from equilibrium require the modeling of significant hydrodynamic stresses responsible in particular for shear-induced migration, a phenomenon known to occur in some industrial processes involving colloids, such as ultrafiltration. The choice of the proper model ingredients requires a knowledge of the domains in parameter space in which colloidal or hydrodynamic effects are dominant. In this article, such a phase diagram is established for a channel flow of charge-stabilized colloids with a version of the suspension balance model including both colloidal and hydrodynamic effects at the continuous level. It is shown that the classical Peclet number is not sufficient to characterize the flow regime. The phase boundary between the colloidal and hydrodynamic regimes exhibits an original shape explained by the dependence of electrostatic interactions with the colloidal surface charge, and in particular by the phenomenon of ionic condensation. We also show that the phase diagram can be predicted based on the knowledge of a rescaled Peclet number comparing the hydrodynamic stress scale to the bulk modulus of the suspension. The criterion determined here provides important guidelines for the efficient modeling of colloidal flows.

Fast, Robust Evaluation of the Equation of State of Suspensions of Charge-Stabilized Colloidal Spheres

Increasing demand is appearing for the fast, robust prediction of the equation of state of colloidal suspensions, notably with a view to using it as input data to calculate transport coefficients in complex flow solvers. This is also of interest in rheological studies, industrial screening tests of new formulations, and the real-time interpretation of osmotic compression experiments, for example. For charge-stabilized spherical particles, the osmotic pressure can be computed with standard liquid theories. However, this calculation can sometimes be lengthy and/or unstable under some physicochemical conditions, a drawback that precludes its use in multiscale flow simulators. As a simple, fast, and robust replacement, the literature reports estimations of the osmotic pressure that have been built by adding the Carnahan-Starling and the cell model pressures (CSCM model). The first contribution is intended to account for colloid-colloid contacts, and the second, for electrostatic effects. This approximation has not yet been thoroughly tested. In this work, the CSCM is evaluated by comparison with data from experiments on silica particles, Monte Carlo simulations, and solutions of the accurate Rogers-Young integral equation scheme with a hard-sphere Yukawa potential obtained from the extrapolated point-charge renormalization method for a wide range of volume fractions, surface charge densities, and interaction ranges. We find that the CSCM is indeed perfectly adequate in the electrostatically concentrated regime, where it can be used from vanishingly small to high surface charge because there is error cancellation between the Carnahan-Starling and cell model contributions at intermediate charge. The CSCM is thus a nice extension of the cell model to liquid-like dense suspensions, which should find application in the domains mentioned above. However, it fails for dilute suspensions with strong electrostatics. In this case, we show that, and explain why, perturbation methods and the rescaled mean spherical approximation are good alternatives in terms of precision, ease of implementation, computational cost, and robustness.

Drying colloidal systems: Laboratory models for a wide range of applications

Abstract.The drying of complex fluids provides a powerful insight into phenomena that take place on time and length scales not normally accessible. An important feature of complex fluids, colloidal dispersions and polymer solutions is their high sensitivity to weak external actions. Thus, the drying of complex fluids involves a large number of physical and chemical processes. The scope of this review is the capacity to tune such systems to reproduce and explore specific properties in a physics laboratory. A wide variety of systems are presented, ranging from functional coatings, food science, cosmetology, medical diagnostics and forensics to geophysics and art.Graphical abstract