Scaling relations for interface motion through disordered media: Application to two-dimensional fluid invasion

Criticality in sheared, disordered solids. I. Rate effects in stress and diffusion

Rate effects in sheared disordered solids are studied using molecular dynamics simulations of binary Lennard-Jones glasses in two and three dimensions. In the quasistatic (QS) regime, systems exhibit critical behavior: the magnitudes of avalanches are power-law distributed with a maximum cutoff that diverges with increasing system size L. With increasing rate, systems move away from the critical yielding point and the average flow stress rises as a power of the strain rate with exponent 1/β, the Herschel-Bulkley exponent. Finite-size scaling collapses of the stress are used to measure β as well as the exponent ν which characterizes the divergence of the correlation length. The stress and kinetic energy per particle experience fluctuations with strain that scale as L^{-d/2}. As the largest avalanche in a system scales as L^{α}, this implies α<d/2. The diffusion rate of particles diverges as a power of decreasing rate before saturating in the QS regime. A scaling theory for the diffusion is derived using the QS avalanche rate distribution and generalized to the finite strain rate regime. This theory is used to collapse curves for different system sizes and confirm β/ν.

Criticality in sheared, disordered solids. II. Correlations in avalanche dynamics

Disordered solids respond to quasistatic shear with intermittent avalanches of plastic activity, an example of the crackling noise observed in many nonequilibrium critical systems. The temporal power spectrum of activity within disordered solids consists of three distinct domains: a novel power-law rise with frequency at low frequencies indicating anticorrelation, white-noise at intermediate frequencies, and a power-law decay at high frequencies. As the strain rate increases, the white-noise regime shrinks and ultimately disappears as the finite strain rate restricts the maximum size of an avalanche. A new strain-rate- and system-size-dependent theory is derived for power spectra in both the quasistatic and finite-strain-rate regimes. This theory is validated using data from overdamped two- and three-dimensional molecular dynamics simulations. We identify important exponents in the yielding transition including the dynamic exponent z which relates the size of an avalanche to its duration, the fractal dimension of avalanches, and the exponent characterizing the divergence in correlations with strain rate. Results are related to temporal correlations within a single avalanche and between multiple avalanches.

Anisotropic avalanches and critical depinning of three-dimensional magnetic domain walls

Simulations with more than 10^{12} spins are used to study the motion of a domain wall driven through a three-dimensional random-field Ising magnet (RFIM) by an external field H. The interface advances in a series of avalanches whose size diverges at a critical external field H_{c}. Finite-size scaling is applied to determine critical exponents and test scaling relations. Growth is intrinsically anisotropic with the height of an avalanche normal to the interface ℓ_{⊥} scaling as the width along the interface ℓ_{∥} to a power χ=0.85±0.01. The total interface roughness is consistent with self-affine scaling with a roughness exponent ζ≈χ that is much larger than values found previously for the RFIM and related models that explicitly break orientational symmetry by requiring the interface to be single-valued. Because the RFIM maintains orientational symmetry, the interface develops overhangs that may surround unfavorable regions to create uninvaded bubbles. Overhangs complicate measures of the roughness exponent but decrease in importance with increasing system size.

Universality of Critically Pinned Interfaces in Two-Dimensional Isotropic Random Media

Based on extensive simulations, we conjecture that critically pinned interfaces in two-dimensional isotropic random media with short-range correlations are always in the universality class of ordinary percolation. Thus, in contrast to interfaces in >2 dimensions, there is no distinction between fractal (i.e., percolative) and rough but nonfractal interfaces. Our claim includes interfaces in zero-temperature random field Ising models (both with and without spontaneous nucleation), in heterogeneous bootstrap percolation, and in susceptible-weakened-infected-removed epidemics. It does not include models with long-range correlations in the randomness and models where overhangs are explicitly forbidden (which would imply nonisotropy of the medium).

Verification of a Dynamic Scaling for the Pair Correlation Function during the Slow Drainage of a Porous Medium

In this Letter we give experimental grounding for the remarkable observation made by Furuberg et al. [Phys. Rev. Lett. 61, 2117 (1988)PRLTAO0031-900710.1103/PhysRevLett.61.2117] of an unusual dynamic scaling for the pair correlation function N(r,t) during the slow drainage of a porous medium. Those authors use an invasion percolation algorithm to show numerically that the probability of invasion of a pore at a distance r away and after a time t from the invasion of another pore scales as N(r,t)∝r^{-1}f(r^{D}/t), where D is the fractal dimension of the invading cluster and the function f(u)∝u^{1.4}, for u≪1 and f(u)∝u^{-0.6}, for u≫1. Our experimental setup allows us to have full access to the spatiotemporal evolution of the invasion, which is used to directly verify this scaling. Additionally, we connect two important theoretical contributions from the literature to explain the functional dependency of N(r,t) and the scaling exponent for the short-time regime (t≪r^{D}). A new theoretical argument is developed to explain the long-time regime exponent (t≫r^{D}).

No-Enclave Percolation Corresponds to Holes in the Cluster Backbone

The no-enclave percolation (NEP) model introduced recently by Sheinman et al. can be mapped to a problem of holes within a standard percolation backbone, and numerical measurements of such holes give the same size-distribution exponent τ=1.82(1) as found for the NEP model. An argument is given that τ=1+d_{B}/2≈1.822 for backbone holes, where d_{B} is the backbone dimension. On the other hand, a model of simple holes within a percolation cluster yields τ=1+d_{f}/2=187/96≈1.948, where d_{f} is the fractal dimension of the cluster, and this value is consistent with the experimental results of gel collapse of Sheinman et al., which give τ=1.91(6). This suggests that the gel clusters are of the universality class of percolation cluster holes. Both models give a discontinuous maximum hole size at p_{c}, signifying explosive percolation behavior.

Effects of Pore-Scale Disorder on Fluid Displacement in Partially-Wettable Porous Media

We present a systematic, quantitative assessment of the impact of pore size disorder and its interplay with flow rates and wettability on immiscible displacement of a viscous fluid. Pore-scale simulations and micromodel experiments show that reducing disorder increases the displacement efficiency and compactness, minimizing the fluid-fluid interfacial area, through (i) trapping at low rates and (ii) viscous fingering at high rates. Increasing the wetting angle suppresses both trapping and fingering, hence reducing the sensitivity of the displacement to the underlying disorder. A modified capillary number Ca^* that includes the impact of disorder λ on viscous forces (through pore connectivity) is direct related to λ, in par with previous works. Our findings bear important consequences on sweep efficiency and fluid mixing and reactions, which are key in applications such as microfluidics to carbon geosequestration, energy recovery, and soil aeration and remediation.

Effect of inertia on sheared disordered solids: critical scaling of avalanches in two and three dimensions

Molecular dynamics simulations with varying damping are used to examine the effects of inertia and spatial dimension on sheared disordered solids in the athermal quasistatic limit. In all cases the distribution of avalanche sizes follows a power law over at least three orders of magnitude in dissipated energy or stress drop. Scaling exponents are determined using finite-size scaling for systems with 10(3)-10(6) particles. Three distinct universality classes are identified corresponding to overdamped and underdamped limits, as well as a crossover damping that separates the two regimes. For each universality class, the exponent describing the avalanche distributions is the same in two and three dimensions. The spatial extent of plastic deformation is proportional to the energy dissipated in an avalanche. Both rise much more rapidly with system size in the underdamped limit where inertia is important. Inertia also lowers the mean energy of configurations sampled by the system and leads to an excess of large events like that seen in earthquake distributions for individual faults. The distribution of stress values during shear narrows to zero with increasing system size and may provide useful information about the size of elemental events in experimental systems. For overdamped and crossover systems the stress variation scales inversely with the square root of the system size. For underdamped systems the variation is determined by the size of the largest events.

Avalanches in strained amorphous solids: does inertia destroy critical behavior?

Simulations are used to determine the effect of inertia on athermal shear of amorphous two-dimensional solids. In the quasistatic limit, shear occurs through a series of rapid avalanches. The distribution of avalanches is analyzed using finite-size scaling with thousands to millions of disks. Inertia takes the system to a new underdamped universality class rather than driving the system away from criticality as previously thought. Scaling exponents are determined for the underdamped and overdamped limits and a critical damping that separates the two regimes. Systems are in the overdamped universality class even when most vibrational modes are underdamped.

Discontinuous percolation transitions in epidemic processes, surface depinning in random media, and Hamiltonian random graphs

Discontinuous percolation transitions and the associated tricritical points are manifest in a wide range of both equilibrium and nonequilibrium cooperative phenomena. To demonstrate this, we present and relate the continuous and first-order behaviors in two different classes of models: The first are generalized epidemic processes that describe in their spatially embedded version--either on or off a regular lattice--compact or fractal cluster growth in random media at zero temperature. A random graph version of these processes is mapped onto a model previously proposed for complex social contagion. We compute detailed phase diagrams and compare our numerical results at the tricritical point in d = 3 with field theory predictions of Janssen et al. [Phys. Rev. E 70, 026114 (2004)]. The second class consists of exponential ("Hamiltonian," i.e., formally equilibrium) random graph models and includes the Strauss and the two-star model, where "chemical potentials" control the densities of links, triangles, or two-stars. When the chemical potentials in either graph model are O(logN), the percolation transition can coincide with a first-order phase transition in the density of links, making the former also discontinuous. Hysteresis loops can then be of mixed order, with second-order behavior for decreasing link fugacity, and a jump (first order) when it increases.

Dynamics of drying in 3D porous media

The drying dynamics in three dimensional porous media are studied with confocal microscopy. We observe abrupt air invasions in size from single particle to hundreds of particles. We show that these result from the strong flow from menisci in large pores to menisci in small pores during drying. This flow causes air invasions to start in large menisci and subsequently spread throughout the entire system. We measure the size and structure of the air invasions and show that they are in accord with invasion percolation. By varying the particle size and contact angle we unambiguously demonstrate that capillary pressure dominates the drying process.

Using wavelets to characterize the wettability of porous materials

Visualization experiments of rate-controlled immiscible displacement of oil by water are performed on a model porous medium of controlled fractional wettability, fabricated by mixing water-wet glass microspheres with strongly oil-wet polytetrafluoroethylene microspheres and packing them between two transparent glass plates. The growth pattern is video recorded and the transient response of the pressure drop across the pore network is measured for various fractions of oil-wet particles. The space-averaged capillary pressure coincides to the pressure drop measured across the porous medium. The oscillating transient signal of the capillary pressure is analyzed with multilevel wavelets to produce the best level wavelet details or best level capillary pressure spectrum (BLCPS) by minimizing the "entropy" of wavelet approximation. Invasion of water in oil-wet areas is reflected in high-frequency and low-amplitude fluctuations of the BLCPS. Correspondingly, invasion of water in water-wet areas is associated with low-frequency and high-amplitude fluctuations of the BLCPS. The displacement growth pattern is reflected in the "energy" and "frequency" of the BLCPS which along with the time-averaged capillary pressure are correlated with two parameters quantifying the spatial variation of wettability over the porous medium: the frontal wettability of the active interface of the fluids before each invasion step and the regional wettability of the area invaded by the displacing fluid.

Self-affinity in the gradient percolation problem

We study the scaling properties of the solid-on-solid front of the infinite cluster in two-dimensional gradient percolation. We show that such an object is self-affine with a Hurst exponent equal to 23 up to a cutoff length approximately g{-4/7}, where g is the gradient. Beyond this length scale, the front position has the character of uncorrelated noise. Importantly, the self-affine behavior is robust even after removing local jumps of the front. The previously observed multiaffinity is due to the dominance of overhangs at small distances in the structure function. This is a crossover effect.

Cluster evolution in steady-state two-phase flow in porous media

We report numerical studies of the cluster development of two-phase flow in a steady-state environment of porous media. This is done by including biperiodic boundary conditions in a two-dimensional flow simulator. Initial transients of wetting and nonwetting phases that evolve before steady state has occurred, undergo a crossover where every initial pattern is broken up. For flow dominated by capillary effects with capillary numbers in order of 10(-5), we find that around a critical saturation of nonwetting fluid the nonwetting clusters of size have a power-law distribution ns similar to s(-tau)with the exponent tau=1.92+/-0.04 for large clusters. This is a lower value than the result for ordinary percolation. We also present scaling relation and time evolution of the structure and global pressure.

Roughness and growth in a continuous fluid invasion model

We have studied interface characteristics in a continuous fluid invasion model, first introduced by Phys. Rev. Lett. 60, 2042 (1988)]. In this model, the interface grows as a response to an applied quasistatic pressure, which induces various types of instabilities. We suggest a variant of the model, which differs from the original model by the order of instabilities treatment. This order represents the relative importance of the physical mechanisms involved in the system. This variant predicts the existence of a third, intermediate regime, in the behavior of the roughness exponent as a function of the wetting properties of the system. The gradual increase of the roughness exponent in this third regime can explain the scattered experimental data for the roughness exponent in the literature. The growth exponent in this model was found to be around zero, due to the initial rough interface.

Experimental study of the immiscible displacement of shear-thinning fluids in pore networks

The pore scale mechanisms and network scale transient pattern of the immiscible displacement of a shear-thinning nonwetting oil phase (NWP) by a Newtonian wetting aqueous phase (WP) are investigated. Visualization imbibition experiments are performed on transparent glass-etched pore networks at a constant unfavorable viscosity ratio and varying values of the capillary number (Ca), and equilibrium contact angle (theta(e)). Dispersions of ozokerite in paraffin oil are used as the shear-thinning NWP, and aqueous solutions of PEG colored with methylene blue are used as the Newtonian WP. At high Ca values, the tip splitting and lateral spreading of WP viscous fingers are suppressed; at intermediate Ca values, the primary viscous fingers expand laterally with the growth of smaller capillary fingers; at low Ca values, network spanning clusters of capillary fingers separated by hydraulically conductive noninvaded zones of NWP arise. The spatial distribution of the mobility of shear-thinning NWP over the pore network is very broad. Pore network regions of low NWP mobility are invaded through a precursor advancement/swelling mechanism even at relatively high Ca and theta(e) values; this mechanism leads to irregular interfacial configurations and retention of a substantial amount of NWP along pore walls; it becomes the dominant mechanism in displacements performed at low Ca and theta(e) values. The residual NWP saturation increases and the end WP relative permeability decreases as Ca increases and both become more sensitive to this parameter as the shear-thinning behavior strengthens. The shear-thinning NWP is primarily entrapped in individual pores of the network rather than in clusters of pores bypassed by the WP. At relatively high flow rates, the amplitude of the variations of pressure drop, caused by fluid redistribution in the pore network, increase with shear-thinning strengthening, whereas at low flow rates, the motion of stable and unstable menisci in pores is reflected in strong pressure drop fluctuations.

Analytic results for scaling function and moments for a different type of avalanche in the bak-sneppen evolution model

Starting from the master equation for the hierarchical structure of avalanches of a different kind within the frame of the Bak-Sneppen evolution model, we derive the exact formula of the scaling function describing the probability distribution of avalanches. The scaling function displays features required by the scaling ansatz and verified by simulations. Using the scaling function we investigate the avalanche moment, denoted by <S(k)>(Delta&fmacr;). It is found that for any non-negative integer k, <S(k)>(Delta&fmacr;) diverges as Delta&fmacr;(-k), which gives an infinite group of exact critical exponents. Simulation outcomes of avalanche moments with k=1,2,3, are found to be consistent with the corresponding analytical results.

Depinning transition and thermal fluctuations in the random-field Ising model

We analyze the depinning transition of a driven interface in the three-dimensional (3D) random field Ising model (RFIM) with quenched disorder by means of Monte Carlo simulations. The interface initially built into the system is perpendicular to the [111] direction of a simple cubic lattice. We introduce an algorithm which is capable of simulating such an interface independent of the considered dimension and time scale. This algorithm is applied to the 3D RFIM to study both the depinning transition and the influence of thermal fluctuations on this transition. It turns out that in the RFIM characteristics of the depinning transition depend crucially on the existence of overhangs. Our analysis yields critical exponents of the interface velocity, the correlation length, and the thermal rounding of the transition. We find numerical evidence for a scaling relation for these exponents and the dimension d of the system.