Darcy's law of yield stress fluids on a treelike network.

Understanding the flow of yield stress fluids in porous media is a major challenge. In particular, experiments and extensive numerical simulations report a nonlinear Darcy law as a function of the pressure gradient. In this letter we consider a treelike porous structure for which the problem of the flow can be resolved exactly due to a mapping with the directed polymer (DP) with disordered bond energies on the Cayley tree. Our results confirm the nonlinear behavior of the flow and expresses its full pressure dependence via the density of low-energy paths of DP restricted to vanishing overlap. These universal predictions are confirmed by extensive numerical simulations.

The fate of shear-oscillated amorphous solids.

The behavior of shear-oscillated amorphous materials is studied using a coarse-grained model. Samples are prepared at different degrees of annealing and then subjected to athermal and quasi-static oscillatory deformations at various fixed amplitudes. The steady-state reached after several oscillations is fully determined by the initial preparation and the oscillation amplitude, as seen from stroboscopic stress and energy measurements. Under small oscillations, poorly annealed materials display shear-annealing, while ultra-stabilized materials are insensitive to them. Yet, beyond a critical oscillation amplitude, both kinds of materials display a discontinuous transition to the same mixed state composed of a fluid shear-band embedded in a marginal solid. Quantitative relations between uniform shear and the steady-state reached with this protocol are established. The transient regime characterizing the growth and the motion of the shear band is also studied.

Darcy's Law for Yield Stress Fluids.

Predicting the flow of non-Newtonian fluids in a porous structure is still a challenging issue due to the interplay between the microscopic disorder and the nonlinear rheology. In this Letter, we study the case of a yield stress fluid in a two-dimensional structure. Thanks to an efficient optimization algorithm, we show that the system undergoes a continuous phase transition in the behavior of the flow, controlled by the applied pressure difference. In analogy with studies of plastic depinning of vortex lattices in high-T_{c} superconductors, we characterize the nonlinearity of the flow curve and relate it to the change in the geometry of the open channels. In particular, close to the transition, a universal scale-free distribution of the channel length is observed and explained theoretically via a mapping to the Kardar-Parisi-Zhang equation.

Experimental Evidence for Three Universality Classes for Reaction Fronts in Disordered Flows.

Self-sustained reaction fronts in a disordered medium subject to an external flow display self-affine roughening, pinning, and depinning transitions. We measure spatial and temporal fluctuations of the front in 1+1 dimensions, controlled by a single parameter, the mean flow velocity. Three distinct universality classes are observed, consistent with the Kardar-Parisi-Zhang (KPZ) class for fast advancing or receding fronts, the quenched KPZ class (positive-qKPZ) when the mean flow approximately cancels the reaction rate, and the negative-qKPZ class for slowly receding fronts. Both qKPZ classes exhibit distinct depinning transitions, in agreement with the theory.

Moving line model and avalanche statistics of Bingham fluid flow in porous media.

In this article, we propose a simple model to understand the critical behavior of path opening during flow of a yield stress fluid in porous media as numerically observed by Chevalier and Talon (2015). This model can be mapped to the problem of a contact line moving in an heterogeneous field. Close to the critical point, this line presents an avalanche dynamic where the front advances by a succession of waiting time and large burst events. These burst events are then related to the non-flowing (i.e. unyielded) areas. Remarkably, the statistics of these areas reproduce the same properties as in the direct numerical simulations. Furthermore, even if our exponents seem to be close to the mean field universal exponents, we report an unusual bump in the distribution which depends on the disorder. Finally, we identify a scaling invariance of the cluster spatial shape that is well fit, to first order, by a self-affine parabola.

History effects on nonwetting fluid residuals during desaturation flow through disordered porous media.

We investigate experimentally the sweeping of a nonwetting fluid by a wetting one in a quasi-two-dimensional porous medium consisting of random obstacles. We focus primarily on the resulting phase distributions and the residual nonwetting phase saturation as a function of the normalized wetting fluid flow rate-the capillary number Ca-at steady state. The wetting liquid is then flowing in the medium partially saturated by immobile nonwetting liquid blobs. The decrease of the nonwetting saturation is an irreversible process that depends strongly on flow history and more specifically on the highest value of Ca reached in the past. At lower Ca values, when capillary forces are dominant, the residual steady state saturation depends significantly on the initial phase configuration. However, at higher Ca, the saturation becomes independent of the history and thus follows a master curve that converges to an asymptotic residual value. Blob sizes range over four orders of magnitude in our experimental domain, following a probability distribution function P that scales with the blob size s as P(s)∝s(-2) for blob sizes larger than the typical pore size. It also exhibits a maximum size cutoff s(max), that decreases as s(max)∝Ca(-1). To determine the flow properties, we have measured the pressure drop (B) versus the flow rate (Ca). In the ranges of low and high Ca values, the relationship between Ca and B is found to be linear, following Darcy's law (B∝Ca). In the intermediate regime, the progressive mobilization of blobs leads to a nonlinear dependence B∝Ca(0.65), due to an increase of the available flow paths.

Generalization of Darcy's law for Bingham fluids in porous media: from flow-field statistics to the flow-rate regimes.

In this paper, we numerically investigate the statistical properties of the nonflowing areas of Bingham fluid in two-dimensional porous media. First, we demonstrate that the size probability distribution of the unyielded clusters follows a power-law decay with a large size cutoff. This cutoff is shown to diverge following a power law as the imposed pressure drop tends to a critical value. In addition, we observe that the exponents are almost identical for two different types of porous media. Finally, those scaling properties allow us to account for the quadratic relationship between the pressure gradient and velocity.

Analysis and improvement of Brinkman lattice Boltzmann schemes: bulk, boundary, interface. Similarity and distinctness with finite elements in heterogeneous porous media.

This work focuses on the numerical solution of the Stokes-Brinkman equation for a voxel-type porous-media grid, resolved by one to eight spacings per permeability contrast of 1 to 10 orders in magnitude. It is first analytically demonstrated that the lattice Boltzmann method (LBM) and the linear-finite-element method (FEM) both suffer from the viscosity correction induced by the linear variation of the resistance with the velocity. This numerical artefact may lead to an apparent negative viscosity in low-permeable blocks, inducing spurious velocity oscillations. The two-relaxation-times (TRT) LBM may control this effect thanks to free-tunable two-rates combination Λ. Moreover, the Brinkman-force-based BF-TRT schemes may maintain the nondimensional Darcy group and produce viscosity-independent permeability provided that the spatial distribution of Λ is fixed independently of the kinematic viscosity. Such a property is lost not only in the BF-BGK scheme but also by "partial bounce-back" TRT gray models, as shown in this work. Further, we propose a consistent and improved IBF-TRT model which vanishes viscosity correction via simple specific adjusting of the viscous-mode relaxation rate to local permeability value. This prevents the model from velocity fluctuations and, in parallel, improves for effective permeability measurements, from porous channel to multidimensions. The framework of our exact analysis employs a symbolic approach developed for both LBM and FEM in single and stratified, unconfined, and bounded channels. It shows that even with similar bulk discretization, BF, IBF, and FEM may manifest quite different velocity profiles on the coarse grids due to their intrinsic contrasts in the setting of interface continuity and no-slip conditions. While FEM enforces them on the grid vertexes, the LBM prescribes them implicitly. We derive effective LBM continuity conditions and show that the heterogeneous viscosity correction impacts them, a property also shared by FEM for shear stress. But, in contrast with FEM, effective velocity conditions in LBM give rise to slip velocity jumps which depend on (i) neighbor permeability values, (ii) resolution, and (iii) control parameter Λ, ranging its reliable values from Poiseuille bounce-back solution in open flow to zero in Darcy's limit. We suggest an "upscaling" algorithm for Λ, from multilayers to multidimensions in random extremely dispersive samples. Finally, on the positive side for LBM besides its overall versatility, the implicit boundary layers allow for smooth accommodation of the flat discontinuous Darcy profiles, quite deficient in FEM.

Strong pinning of propagation fronts in adverse flow.

Reaction fronts evolving in a porous medium exhibit a rich dynamical behavior. In the presence of an adverse flow, experiments show that the front slows down and eventually gets pinned, displaying a particular sawtooth shape. Extensive numerical simulations of the hydrodynamic equations confirm the experimental observations. Here we propose a stylized model, predicting two possible outcomes of the experiments for large adverse flow: either the front develops a sawtooth shape or it acquires a complicated structure with islands and overhangs. A simple criterion allows one to distinguish between the two scenarios and its validity is reproduced by direct hydrodynamical simulations. Our model gives a better understanding of the transition and is relevant in a variety of domains, when the pinning regime is strong and only relies on a small number of sites.

On the determination of a generalized Darcy equation for yield-stress fluid in porous media using a Lattice-Boltzmann TRT scheme.

Simulating flow of a Bingham fluid in porous media still remains a challenging task as the yield stress may significantly alter the numerical stability and precision. We present a Lattice-Boltzmann TRT scheme that allows the resolution of this type of flow in stochastically reconstructed porous media. LB methods have an intrinsic error associated to the boundary conditions. Depending on the schemes this error might be directly linked to the effective viscosity. As for non-Newtonian fluids viscosity varies in space the error becomes inhomogeneous and very important. In contrast to that, the TRT scheme does not present this deficiency and is therefore adequate to be used for simulations of non-Newtonian fluid flow. We simulated Bingham fluid flow in porous media and determined a generalized Darcy equation depending on the yield stress, the effective viscosity, the pressure drop and a characteristic length of the porous medium. By evaluating the flow in the porous structure, we distinguished three different scaling regimes. Regime I corresponds to the situation where fluid is flowing in only one channel. Here, the relation between flow rate and pressure drop is given by the non-Newtonian Poiseuille law. During Regime II an increase in pressure triggers the opening of new paths and the relation between flow rate and the difference in pressure to the critical yield pressure becomes quadratic: [Formula: see text]. Finally, Regime III corresponds to the situation where all the fluid is flowing. In this case, [Formula: see text].

Low Reynolds number suspension gravity currents.

The extension of a gravity current in a lock-exchange problem, proceeds as square root of time in the viscous-buoyancy phase, where there is a balance between gravitational and viscous forces. In the presence of particles however, this scenario is drastically altered, because sedimentation reduces the motive gravitational force and introduces a finite distance and time at which the gravity current halts. We investigate the spreading of low Reynolds number suspension gravity currents using a novel approach based on the Lattice-Boltzmann (LB) method. The suspension is modeled as a continuous medium with a concentration-dependent viscosity. The settling of particles is simulated using a drift flux function approach that enables us to capture sudden discontinuities in particle concentration that travel as kinematic shock waves. Thereafter a numerical investigation of lock-exchange flows between pure fluids of unequal viscosity, reveals the existence of wall layers which reduce the spreading rate substantially compared to the lubrication theory prediction. In suspension gravity currents, we observe that the settling of particles leads to the formation of two additional fronts: a horizontal front near the top that descends vertically and a sediment layer at the bottom which aggrandises due to deposition of particles. Three phases are identified in the spreading process: the final corresponding to the mutual approach of the two horizontal fronts while the laterally advancing front halts indicating that the suspension current stops even before all the particles have settled. The first two regimes represent a constant and a decreasing spreading rate respectively. Finally we conduct experiments to substantiate the conclusions of our numerical and theoretical investigation.

Autocatalytic reaction fronts inside a porous medium of glass spheres.

We analyze experimentally chemical wave propagation in the disordered flow field of a porous medium. The reaction fronts travel at a constant velocity that drastically depends on the mean flow direction and rate. The fronts may propagate either downstream and upstream but, surprisingly, they remain static over a range of flow rate values. Resulting from the competition between the chemical reaction and the disordered flow field, these frozen fronts display a particular sawtooth shape. The frozen regime is likely to be associated with front pinning in low velocity zones, the number of which varies with the ratio of the mean flow and the chemical front velocities.

Permeability of self-affine aperture fields.

We introduce a model that allows for the prediction of the permeability of self-affine rough channels (one-dimensional fracture) and two-dimensional fractures over a wide range of apertures. In the lubrication approximation, the permeability shows three different scaling regimes. For fractures with a large mean aperture or an aperture small enough to the permeability being close to disappearing, the permeability scales as the cube of the aperture when the zero level of the aperture is set to coincide with the disappearance of the permeability. Between these two regimes, there is a third regime where the scaling is due to the self-affine roughness. For rough channels, the exponent is found to be 3-1/H, where H is the Hurst exponent. For two-dimensional fractures, it is necessary to introduce an equivalent aperture b(c) to make the scaling regime apparent. b(c) is defined as the hydraulic aperture of the most restrictive barrier crossing the fracture normal to the flow direction. This regime is characterized by an exponent higher than that for the one-dimensional case: it is 2.25 for H=0.8 and 2.16 for H=0.3.

Critical behavior of the banded-unbanded spherulite transition in a mixture of ethylene carbonate with polyacrylonitrile.

We investigate the transition from unbanded to banded spherulitic growth in mixtures of ethylene carbonate with polyacrylonitrile. By carefully considering systematic errors, we show that the band spacing diverges with a power-law form showing scaling over nearly two decades. We also observe that the bands disorder as the transition point is approached. The critical exponent is nonclassical. One possible explanation is that the nonequilibrium transition is actually weakly first order (subcritical).