Yielding to percolation: a universal scale

A theoretical and computational study analysing the initiation of yield-stress fluids percolation in porous media is presented. Yield-stress fluid flows through porous media are complicated due to the non-linear rheological behaviour of this type of fluids, rendering the conventional Darcy type approach invalid. A critical pressure gradient must be exceeded to commence the flow of a yield-stress fluid in a porous medium. As the first step in generalising the Darcy law for yield-stress fluids, a universal scale based on the variational formulation of the energy equation is derived for the critical pressure gradient which reduces to purely geometrical feature of the porous media. The presented scaling is then validated by both exhaustive numerical simulations (using an adaptive finite element approach based on the augmented Lagrangian method), and also the previously published data. The considered porous media is constructed by randomised obstacles with various topologies; namely, square, circular and alternatively polygonal obstacles which are mimicked based on Voronoi tessellation of circular cases. Moreover, computations for the bi-dispersed obstacle cases are performed which further demonstrate the validity of the proposed universal scaling

General hydrodynamic features of elastoviscoplastic fluid flows through randomised porous media

A numerical study of yield-stress fluids flowing in porous media is presented. The porous media is randomly constructed by non-overlapping mono-dispersed circular obstacles. Two class of rheological models are investigated: elastoviscoplastic fluids (i.e. Saramito model) and viscoplastic fluids (i.e. Bingham model). A wide range of practical Weissenberg and Bingham numbers is studied at three different levels of porosities of the media. The emphasis is on revealing some physical transport mechanisms of yield-stress fluids in porous media when the elastic behaviour of this kind of fluids is incorporated. Thus, computations of elastoviscoplastic fluids are performed and are compared with the viscoplastic fluid flow properties. At a constant Weissenberg number, the pressure drop increases both with the Bingham number and the solid volume fraction of obstacles. However, the effect of elasticity is less trivial. At low Bingham numbers, the pressure drop of an elastoviscoplastic fluid increases compared to a viscoplastic fluid, while at high Bingham numbers we observe drag reduction by elasticity. At the yield limit (i.e. infinitely large Bingham numbers), elasticity of the fluid systematically promotes yielding: elastic stresses help the fluid to overcome the yield stress resistance at smaller pressure gradients. We observe that elastic effects increase with both Weissenberg and Bingham numbers. In both cases, elastic effects finally make the elastoviscoplastic flow unsteady, which consequently can result in chaos and turbulence. Keywords: Yield-stress fluids; Viscoplastic fluids; Elastoviscoplastic fluids; Porous medi

Yielding to percolation : a universal scale

A theoretical and computational study analysing the initiation of yield-stress fluid percolation in porous media is presented. Yield-stress fluid flows through porous media are complicated due to the nonlinear rheological behaviour of this type of fluid, rendering the conventional Darcy type approach invalid. A critical pressure gradient must be exceeded to commence the flow of a yield-stress fluid in a porous medium. As the first step in generalising the Darcy law for yield-stress fluids, a universal scale based on the variational formulation of the energy equation is derived for the critical pressure gradient which reduces to the purely geometrical feature of the porous media. The presented scaling is then validated by both exhaustive numerical simulations (using an adaptive finite element approach based on the augmented Lagrangian method), and also the previously published data. The considered porous media are constructed by randomised obstacles with various topologies; namely square, circular and alternatively polygonal obstacles which are mimicked based on Voronoi tessellation of circular cases. Moreover, computations for the bidispersed obstacle cases are performed which further demonstrate the validity of the proposed universal scaling

Sliding flows of yield-stress fluids

A theoretical and numerical study of complex sliding flows of yield-stress fluids is presented. Yield-stress fluids are known to slide over solid surfaces if the tangential stress exceeds the sliding yield stress. The sliding may occur due to various microscopic phenomena such as the formation of an infinitesimal lubrication layer of the solvent and/or elastic deformation of the suspended soft particles in the vicinity of the solid surfaces. This leads to a 'stick-slip' law which complicates the modelling and analysis of the hydrodynamic characteristics of the yield-stress fluid flow. In the present study, we formulate the problem of sliding flow beyond one-dimensional rheometric flows. Then, a numerical scheme based on the augmented Lagrangian method is presented to attack these kind of problems. Theoretical tools are developed for analysing the flow/no-flow limit. The whole framework is benchmarked in planar Poiseuille flow and validated against analytical solutions. Then two more complex physical problems are investigated: slippery particle sedimentation and pressure-driven sliding flow in porous media. The yield limit is addressed in detail for both flow cases. In the particle sedimentation problem, method of characteristics - slipline method - in the presence of slip is revisited from the perfectly plastic mechanics and used as a helpful tool in addressing the yield limit. Finally, flows through model and randomized porous media are studied. The randomized configuration is chosen to capture more sophisticated aspects of the yield-stress fluid flows in porous media at the yield limit - channelization

Stability of particles inside yield-stress fluid Poiseuille flows

The stability of neutrally and non-neutrally buoyant particles immersed in a plane Poiseuille flow of a yield-stress fluid (Bingham fluid) is addressed numerically. Particles being carried by the yield-stress fluid can behave in different ways: they might (i) migrate inside the yielded regions or (ii) be transported without any relative motion inside the unyielded region if the yield stress is large enough compared to the buoyancy stress and the other stresses acting on the particles. Knowing the static stability of particles inside a bath of quiescent yield-stress fluid (Chaparian & Frigaard, J.A Fluid Mech., vol.A 819, 2017, pp.A 311-351), we analyse the latter behaviour when the yield-stress fluid Poiseuille flow is host to two-dimensional particles. Numerical experiments reveal that particles lose their stability (i.e. break the unyielded plug and sediment/migrate) with smaller buoyancy compared to the sedimentation inside a bath of quiescent yield-stress fluid, because of the inherent shear stress in the Poiseuille flow. The key parameter in interpreting the present results is the position of the particle relative to the position of the yield surface in the undisturbed flow (in the absence of any particle): the larger the portion of a particle located inside the undisturbed sheared regions, the more likely is the particle to be unstable. Yet, we find that the core unyielded plug can grow locally to some extent to contain the particles. This picture holds even for neutrally buoyant particles, although they are strictly stable when they are located wholly inside the undisturbed plug. We propose scalings for all cases

Atypical plug formation in internal elastoviscoplastic fluid flows over a non-smooth topology

An experimental and computational investigation of the internal flow of elastoviscoplastic fluids over non-smooth topologies is presented in two complimentary studies. In the first study, we visualize the creeping flow of a Carbopol gel over a cavity embedded in a thin slot using Optical Coherence Tomography (OCT) and confocal microscopy. We measure the size and shape of the plug as a function of Bingham and Weissenberg numbers. An asymmetry in the plug shape is observed which is also evident in our second study -- numerical simulations using adaptive finite element method based upon an augmented Lagrangian scheme. We quantify the asymmetry and present the results as a function of the product of the Weissenberg and Bingham numbers which collapse onto a single curve for each of these geometries. These findings underscore the theoretical underpinnings of the synergy between elasticity and plasticity of these complex fluids

Clouds of bubbles in a viscoplastic fluid

Viscoplastic fluids can hold bubbles/particles stationary by balancing the buoyancy stress with the yield stress - the key parameter here is the yield number, the ratio of the yield stress to the buoyancy stress. In the present study, we investigate a suspension of bubbles in a yield-stress fluid. More precisely, we compute how much is the gas fraction that could be held trapped in a yield-stress fluid without motion. Here the goal is to shed light on how the bubbles feel their neighbours through the stress field and to compute the critical yield number for a bubble cloud beyond which the flow is suppressed. We perform two-dimensional computations in a full periodic box with randomized positions of the monosized circular bubbles. A large number of configurations are investigated to obtain statistically converged results. We intuitively expect that for higher volume fractions, the critical yield number is larger. Not only here do we establish that this is the case, but also we show that short-range interactions of bubbles increase the critical yield number even more dramatically for bubble clouds. The results show that the critical yield number is a linear function of volume fraction in the dilute regime. An algebraic expression model is given to approximate the critical yield number (semi-empirically) based on the numerical experiment in the studied range of 0 ≤ φ ≤ 0.31,, together with lower and upper estimates

Flow onset for a single bubble in a yield-stress fluid

We use computational methods to determine the minimal yield stress required in order to hold static a buoyant bubble in a yield-stress liquid. The static limit is governed by the bubble shape, the dimensionless surface tension and the ratio of the yield stress to the buoyancy stress . For a given geometry, bubbles are static for Y_c]]>, which we determine for a range of shapes. Given that surface tension is negligible, long prolate bubbles require larger yield stress to hold static compared with oblate bubbles. Non-zero increases and for large the yield-capillary number determines the static boundary. In this limit, although bubble shape is important, bubble orientation is not. Two-dimensional planar and axisymmetric bubbles are studied

Global H-NS counter-silencing by LuxR activates quorum sensing gene expression

Bacteria coordinate cellular behaviors using a cell-cell communication system termed quorum sensing. In Vibrio harveyi, the master quorum sensing transcription factor LuxR directly regulates \u3e100 genes in response to changes in population density. Here, we show that LuxR derepresses quorum sensing loci by competing with H-NS, a global transcriptional repressor that oligomerizes on DNA to form filaments and bridges. We first identified H-NS as a repressor of bioluminescence gene expression, for which LuxR is a required activator. In an hns deletion strain, LuxR is no longer necessary for transcription activation of the bioluminescence genes, suggesting that the primary role of LuxR is to displace H-NS to derepress gene expression. Using RNA-seq and ChIP-seq, we determined that H-NS and LuxR co-regulate and co-occupy 28 promoters driving expression of 63 genes across the genome. ChIP-PCR assays show that as autoinducer concentration increases, LuxR protein accumulates at co-occupied promoters while H-NS protein disperses. LuxR is sufficient to evict H-NS from promoter DNA in vitro, which is dependent on LuxR DNA binding activity. From these findings, we propose a model in which LuxR serves as a counter-silencer at H-NS-repressed quorum sensing loci by disrupting H-NS nucleoprotein complexes that block transcription