Fingering Instability Research Articles

Direct numerical simulation of droplet impact onto dry stationary and moving walls at low to high Weber numbers

The accurate prediction of the drop-wall impact dynamics and regime are of interest in many engineering applications. Considerable progress has been made in different elements of Direct Numerical Simulation (DNS) approaches. However, identifying the optimal combination of improved submodels that perform well under a wide range of operating conditions, covering low to high Weber numbers with varying surface wettabilities and wall motions, still requires further investigation. Here, we conduct a comparative study on the performance of key DNS elements, including the Interface Tracking (IT) algorithm for the widely-used volume-of-fluid approach, Contact Angle (CA), Contact-Line Velocity (CLV), Instability Actuation (IA), and numerical discretization. The study spans a broad range of Weber numbers, encompassing various impact phenomena such as drop oscillation, partial rebound, fingering, and splash. It is concluded that the Kistler CA model provides the most accurate predictions among the models considered here. In terms of IT, the Multidimensional Universal Limiter with Explicit Solution (MULES) algorithm is the most efficient one, especially for moderate to high Weber numbers. For high Weber numbers, involving the fingering instability and splash events, while applying the IA mechanism slightly improves the results, using a high-quality fine-enough grid and appropriate numerical discretization scheme to control the dispersion and dissipation errors are more important ingredients. It is shown that the recommended model combination is also able to predict the important features of drop impact on moving walls with reasonable accuracy.

Supramolecular Assemblies in Active Motor-Filament Systems: Micelles, Bilayers, and Foams

Active matter systems evade the constraints of thermal equilibrium, leading to the emergence of intriguing collective behavior. A paradigmatic example is given by motor-filament mixtures, where the motion of motor proteins drives alignment and sliding interactions between filaments and their self-organization into macroscopic structures. After defining a microscopic model for these systems, we derive continuum equations, exhibiting the formation of active supramolecular assemblies such as micelles, bilayers, and foams. The transition between these structures is driven by a branching instability, which destabilizes the orientational order within the micelles, leading to the growth of bilayers at high microtubule densities. Additionally, we identify a fingering instability, modulating the shape of the micelle interface at high motor densities. We study the role of various mechanisms in these two instabilities, such as contractility, active splay, and anchoring, allowing for generalization beyond the system considered here. Published by the American Physical Society 2024

Nonlinear flow phenomenon of a power-law non-Newtonian fluid falling down a cylinder surface

In this paper, we present a comprehensive study of the fingering phenomenon of a power-law non-Newtonian fluid falling down a cylinder surface. A theoretical analysis is firstly carried out and the governing equation describing the film thickness is established for the non-Newtonian fluid denoted by a power-law index n. Using the lubrication theory with dimensionless variables, the partial differential equation for the film thickness is derived, and both two-dimensional flow and three-dimensional flow are investigated. The traveling wave characteristic of the two-dimensional flow is displayed by using the finite difference scheme associated with a Newton iteration technique. The effect of changing the radius of the cylinder, precursor-layer thickness and power-law index is then considered through examining the variation of the capillary waves. Furthermore in three-dimensional complex flow, the fingering patterns for different parameters are simulated. The linear stability analysis based on the traveling wave solutions is given to elucidate the influence of different physical parameters, explaining the physical mechanism of nonlinear flow behaviors in two dimension and three dimension. Results from linear stability analysis show that the power-law index reinforces the fingering instability whether the fluid is shear-thinning or shear-thickening. Moreover, the numerical results illustrate that increasing the radii of the cylinder expands the unstable region. The flow profiles for different power-law coefficients are consistent with the result of the linear stability analysis, proving the unstable role of the power-law coefficient.

Emergence of lobed wakes during the sedimentation of spheres in viscoelastic fluids

The motion of rigid particles in complex fluids is ubiquitous in natural and industrial processes. The most popular toy model for understanding the physics of such systems is the settling of a solid sphere in a viscoelastic fluid. There is general agreement that an elastic wake develops downstream of the sphere, causing the breakage of fore-and-aft symmetry, while the flow remains axisymmetric, independent of fluid viscoelasticity and flow conditions. Using a continuum mechanics model, we reveal that axisymmetry holds only for weak viscoelastic flows. Beyond a critical value of the settling velocity, steady, non-axisymmetric disturbances develop peripherally of the rear pole of the sphere, giving rise to a four-lobed fingering instability. The transition from axisymmetric to non-axisymmetric flow fields is characterized by a regular bifurcation and depends solely on the interplay between shear and extensional properties of the viscoelastic fluid under different flow regimes. At higher settling velocities, each lobe tip is split into two new lobes, resembling fractal fingering in interfacial flows. For the first time, we capture an elastic fingering instability under steady-state conditions, and provide the missing information for understanding and predicting such instabilities in the response of viscoelastic fluids and soft media.

Synergies of storing hydrogen at the crest of CO2${\rm CO}_{2}$ or other gas storage

AbstractThere are mutual benefits in storing in sedimentary reservoirs jointly with another gas serving as a cushion gas, such as the of a carbon capture and storage (CCS) operation or the natural gas of seasonal storage or left in a depleted hydrocarbon reservoir. When occupies the crest of the reservoir, the presence of either gas is beneficial to the other. reinforces the sealing efficiency of the caprock due to its very favorable interfacial properties with respect to brine and rock‐forming minerals. storage safety and capacity are also increased with cushion gases such as , which alleviate the buoyancy pressure at the top of the gas column. The potential drawback of this storage scheme is gas/gas mixing, which can, however, be strongly reduced if, by an appropriate choice of well completion and placement, is positioned in the upper zones of the reservoir, and its injection rate is kept below a critical value corresponding to the incipient fingering instability of the gas mixing zone. This value, which depends on reservoir permeability, the dip angle of the mixing front, and how density varies with viscosity in the mixing front, turns out to be well above practical injection rates. Therefore, dispersive mixing is the only cause of front spreading, which is acceptable for not‐too‐heterogeneous reservoirs. The mutual benefits identified in this study are the strongest when the cushion is made up of dense , which suggests that the crest of offshore storage reservoirs is a good candidate for storage. © 2024 The Authors. Greenhouse Gases: Science and Technology published by Society of Chemical Industry and John Wiley & Sons Ltd.

Study on two-phase displacement law considering mass transfer and diffusion

The influence of channeling on viscous fingering instability of CO2 miscible displacement is studied in the paper. Due to the viscous fingering not easily observed in porous media, channeling is used to simplify the viscous fingering instability. We adopted nonlinear simulation to investigate the development of viscous fingering instability during the displacement of Newtonian fluids in a channel by miscible fluids, and the influence of different Pe and different viscosity ratio R was studied. Under homogeneous conditions, when R is the same, the larger the Pe is, the more obvious the convection in the process of CO2 displacement is, the earlier the viscous fingering occurs, and the shorter the time for the finger to break through to the right boundary is. When Pe is the same, the larger the R is, the more unstable the contact area of miscible displacement becomes. More finger structures appear, and more complex fingering phenomena occur. The larger the R is, the earlier the fingering phenomenon appears, and the earlier it breaks through to the right boundary. In addition, we studied the change in Relative Mixing Length (RML) during the diffusion process quantitatively. Finally, our investigation delved into the impact of heterogeneity on viscous fingering. We observed that under significant heterogeneity, viscous fingering tends to manifest preferentially in the direction of increasing permeability.

Measurement of pressure gradients near the interface in the viscous fingering instability

The viscous fingering instability, which forms when a less-viscous fluid invades a more-viscous one within a confined geometry, is an iconic system for studying pattern formation. For both miscible and immiscible fluid pairs the growth dynamics change after the initial instability onset and the global structures, typical of late-time growth, are governed by the viscosity ratio. Here we introduce an experimental technique to measure flow throughout the inner and outer fluids. This probes the existence of a new length scale associated with the local pressure gradients around the interface and allows us to compare our results to the predictions of a previously proposed model for late-time finger growth. Published by the American Physical Society 2024

Trade Off between Hydrodynamic and Thermodynamic Forces at the Liquid-Liquid Interface.

Viscous fingering (VF) instability has been investigated in the case of a partially miscible binary system by nonlinear numerical simulations. Partially miscible fluid systems offer the possibility of phase separation coupled with VF instability. The thermodynamics of such systems are governed by the Margules parameter (interaction parameter) as well as the fluid concentrations. Kinetics of the decomposition is also influenced by dynamical parameters such as the viscosity of the fluid, which incidentally also affects the hydrodynamic forces. Here, we explore the effects of concentration and Margules parameter in order to ascertain the trade-offs incurred between hydrodynamic and thermodynamic effects at the interface as well as the thermodynamics of the bulk. Based on the Gibb's free energy versus concentration curve, we select concentrations (i) outside spinodal and binodal regions, (ii) within binodal but outside the spinodal, and (iii) within the spinodal curve. We solve the modified Cahn-Hilliard-Hele-Shaw equation employing the COMSOL Multiphysics software. Applying high-resolution numerical simulations, we show a strong dependence of the thermodynamic forces on the concentration of the mixtures. Rapid phase separation and hence a faster rate of droplet formation have been found when the concentration lies inside the spinodal region. Further, we have investigated the correlation between the fractal dimension and dynamics of the system. The spatiotemporal studies presented in this work clearly illustrate the competition between hydrodynamic and thermodynamic forces and provide insights on the kinetics of decomposition and growth of interfacial instabilities.

Air invasion into three-dimensional foam induces viscous fingering instabilities

We conducted an experimental investigation to examine the immiscible radial displacement flows of air invading three-dimensional foam in a Hele-Shaw cell. Our study successfully identified three distinct flow regimes. In the initial regime, characterized by relatively low fingertip velocities, the foam underwent a slow displacement through plug flow. During this process, the three-phase contact lines slipped at the cell walls. Notably, we discovered that the air injection pressure exhibited a proportional relationship with the power of the fingertip velocity. This relationship demonstrated excellent agreement with a power law, where the exponent was determined to be 2/3. Transitioning to the second regime, we observed relatively high velocities, resulting in the displacement of the foam as a plug within single layers of foam bubbles. The movement of these bubbles near the cell walls was notably slower. Similar to the first regime, the behavior in this regime also adhered to a power law. In the third regime, which manifested at higher air injection pressures, the development of air fingers occurred through narrow channels. These channels had the potential to isolate the air fingers as they underwent a process of "healing." Furthermore, our results unveiled a significant finding that the width of the air fingers exhibited a continuous scaling with the air injection pressure, irrespective of the flow regimes being observed.

Interfacial instabilities in confined displacements involving non-Newtonian fluids

The growth of interfacial instabilities during fluid displacements can be driven by gradients in pressure, viscosity and surface tension, and by applying external fields. Since displacements of non-Newtonian fluids such as polymer solutions, colloidal and granular slurries are ubiquitous in natural and industrial processes, understanding the growth mechanisms and fully developed morphologies of interfacial patterns involving non-Newtonian fluids is extremely important. In this perspective, we focus on displacement experiments, wherein competitions between capillary, viscous, elastic and frictional forces drive the onset and growth of primarily viscous fingering instabilities in confined geometries. We conclude by highlighting several exciting open problems in this research area.

Controlling viscous fingering instabilities of complex fluids

Despite their aesthetic elegance, wavy or fingering patterns emerge when a fluid of low viscosity pushes another immiscible fluid of high viscosity in a porous medium, producing an incomplete sweep and hampering several crucial technologies. Some examples include chromatography, printing, coating flows, oil-well cementing, as well as large-scale technologies of groundwater and enhanced oil recovery. Controlling such fingering instabilities is notoriously challenging and unresolved for complex fluids of varying viscosity because the fluids’ mobility contrast is often predetermined and yet the predominant drive in determining a stable, flat or unstable, wavy interface. Here we show, experimentally and theoretically, how to suppress or control the primary viscous fingering patterns of a common type of complex fluids (of shear-thinning with a low yield stress) using a radially tapered cell of linearly varying gap thickness, h(r). Experimentally, we displace a complex viscous (PAA) solution with gas under a constant flow rate (Q), varied between 0.02 and 2 slpm (standard liter per minute), in a radially converging cell with a constant gap-thickness gradient, α=dh/dr<0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha = dh/dr < 0$$\\end{document}. A stable, uniform interface emerges at low Q and in a steeper cell (i.e., greater |α|\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$|\\alpha |$$\\end{document}) for the complex fluids, whereas unstable fingering pattern at high Q and smaller |α|\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$|\\alpha |$$\\end{document}. Our theoretical predictions with a simplified linear stability analysis show an agreeable stability criterion with experimental data, quantitatively offering strategies to control complex fluid-fluid patterns and displacements in microfluidics and porous media.

Spreading of volatile droplets in a humidity-controlled environment.

When a pure ethanol droplet is deposited on a dry, wettable and conductive substrate, it is expected to spread into a thin, uniform film. Here, we demonstrate that this uniform spreading behaviour can be altered significantly by controlling the ambient relative humidity. We show that higher relative humidity not only promotes faster spreading of the droplet, it also destabilizes the moving contact line, resulting in a fingering instability. We observe that these effects primarily emerge due to the hygroscopic nature of the pure droplet, which eventually leads to solutal-Marangoni effects. Additionally, heat transfer between the evaporating droplet and the underlying substrate also plays a crucial role in the overall dynamics. Thus, the overall spreading of a pure hygroscopic droplet is determined by a delicate interplay between solutal and thermal Marangoni effects.

The subdivision-based IGA-EIEQ numerical scheme for the Cahn–Hilliard–Darcy system of two-phase Hele-Shaw flow on complex curved surfaces

This paper investigates the numerical approximation of a two-phase Hele-Shaw fluid flow system, which consists of a phase-field Cahn–Hilliard model coupled with the Darcy flow, confined on a complex curved surface. To avoid geometric errors caused by the surface approximation and additional approximation errors introduced by numerical methods, for spatial discretization, we consider the recently developed strategy of subdivision-based isogeometric analysis (IGA), where we can accurately represent complex surfaces of arbitrary topology due to its high smoothness and hierarchical refinement properties. By incorporating multiple methods of temporal discretizations, including the Explicit-Invariant-Energy-Quadratization (EIEQ) approach for linearizing nonlinear potentials, the Zero-Energy-Contribution (ZEC) method for decoupling of nonlinear couplings, and the projection method for decoupling of linear couplings between the fluid velocity and pressure, and the finite-element method based on IGA for spatial discretizations, we arrive at a fully discrete scheme that possesses desirable properties such as decoupling, linearity, second-order accuracy in time, and unconditional energy stability. At each time step, only several elliptic equations with constant coefficients are needed to be solved. We also present the rigorous proof of unconditional energy stability, accompanied by numerous numerical examples to verify the stability and accuracy of the numerical scheme, such as the Saffman–Taylor fingering instability occurring on various complex curved surfaces induced by rotational forces.

Effects of viscosity ratio and surface wettability on viscous fingering instability in rectangular channel

This paper explores the impact of viscosity ratio and surface wettability on immiscible viscous fingering instability within a rectangular channel. Numerical investigations are conducted across a range of viscosity ratios (VR) from 0.0009 to 0.5 and wall wettability ([Formula: see text]) from 15° to 150°. The volume of fluid (VOF) model is employed to track the development of finger-shaped instability at the fluid interface. Our results indicate that higher viscosity ratios lead to increased displacement efficiency. Additionally, we find the formation of necking at low VR, which diminishes at higher VR values. The finger-shaped pattern splits into two parts at a wettability of 15°; beyond this threshold, no such splitting occurs. Furthermore, a transition from hydrophilic to superhydrophobic wettability abolishes necking, resulting in enhanced displacement efficiency. Notably, as wettability shifts from hydrophilic to super hydrophobic, instability shifts toward the left side. These findings hold relevance for applications in drug delivery, clinical processes and oil recovery.

Electroconvective viscous fingering in a single polyelectrolyte fluid on a charge selective surface

When a low-viscosity fluid displaces into a higher-viscosity fluid, the liquid-liquid interface becomes unstable causing finger-like patterns. This viscous fingering instability has been widely observed in nature and engineering systems with two adjoined fluids. Here, we demonstrate a hitherto-unrealizable viscous fingering in a single fluid-solid interface. In a single polyelectrolyte fluid on a charge selective surface, selective ion rejection through the surface initiates i) stepwise ion concentration and viscosity gradient boundaries in the fluid and ii) electroconvective vortices on the surface. As the vortices grow, the viscosity gradient boundary pushes away from the surface, resulting viscous fingering. Comparable to conventional one with two fluids, i) a viscosity ratio (M) governs the onset of this electroconvective viscous fingering, and ii) the boundary properties (finger velocity and rheological effects) - represented by M, electric Rayleigh ({{Ra}}_{E}), Schmidt ({Sc}), and Deborah ({De}) numbers - determine finger shapes (straight v.s. ramified, the onset length of fingering, and relative finger width). With controllable onset and shape, the mechanism of electroconvective viscous fingering offers new possibilities for manipulating ion transport and dendritic instability in electrochemical systems.

The thermo-viscous fingering instability of a cooling spreading liquid dome

We investigate a theoretical model of a molten viscous planar liquid dome spreading under gravity over an inclined substrate. The liquid in the dome cools as it spreads, losing its heat to the surrounding colder air and substrate. Coupled nonlinear evolution equations for the dome's thickness and temperature describing the spreading flow are derived employing the lubrication approximation. The coupling between the flow and cooling is via a temperature-dependent viscosity. For intermediate Péclet numbers, a new one-dimensional free surface shape is identified. In this solution, the hotter and more mobile liquid piles up behind the dome's colder and less mobile leading edge, forming a distinct elevated ridge at the flow front. The ridge solution is mapped in parameter space. The transverse stability of the one-dimensional ridge solution is investigated using linear stability analysis and numerical simulations. The existence of a thermo-viscous fingering instability is revealed. For this instability to occur, the presence of the ridge is shown to be necessary. Two-dimensional simulations confirm the stability analysis elucidating the underlying thermo-viscous mechanism.

Research on the Interfacial Instability of Non-Newtonian Fluid Displacement Using Flow Geometry

The variation of the classical viscous fingering instability is studied numerically in this work. An investigation of the viscous fingering phenomenon of immiscible displacement in the Hele–Shaw cell (HSC), where the displaced fluid is a shear-thinning fluid, was carried out numerically using the volume of fluid (VOF) method by adding a minor depth gradient or altering the geometry of the top plate in the HSC. The findings demonstrate how the presence of depth gradients can change the stability of the interface and offer a chance to regulate and adapt the fingering instability in response to the viscous fingering properties of air driving non-Newtonian fluids under various depth gradients. The relative breadth will shrink under the influence of the depth gradient, and the negative consequences of the gradient will be increasingly noticeable. Specifically, under different power-law indices, we found that with the enhancement of shear-thinning characteristics (lower power-law exponent n) in both positive and negative depth gradients, the fingers that protrude from the viscous fingers become shorter and thicker, resulting in higher displacement efficiency. Additionally, several modifications were performed to the upper plate’s design, and the findings revealed that the shape had no effect on the viscous fingering and only had an impact on the longitudinal amplitude. Based on the aforementioned traits, we may alter the HSC’s form or depth gradient to provide high-quality and effective work.

Effects of dilute low molecular weight poly(ethylene oxide) solutions in immiscible radial viscous fingering instabilities

We experimentally study the effects of normal stress differences in the immiscible radial viscous fingering instability in a Hele–Shaw cell. Dilute low molecular weight poly(ethylene oxide) (PEO) solutions are used as the displaced fluid to focus on elastic effects without shear-thinning behavior. The molecular weight of PEO is varied to investigate the effects of normal stress differences. The experimental observations reveal that nonmonotonic and opposing effects are evident depending on the molecular weight of the PEO and the stage of the radial viscous fingering evolution. Decreases in the PEO molecular weight reduce the number of fingers and widen the finger width in the early stage. However, the increase in the PEO molecular weight promotes tip splitting and narrows finger width in the early stage but suppresses tip splitting in the intermediate stage. Weissenberg numbers are estimated at different stages of radial viscous fingering instabilities. Tip splitting occurs at the highest Weissenberg number covered in this study and suppression of tip splitting is observed at intermediate Weissenberg numbers. At low Weissenberg numbers, we observe an increased finger width and a reduced number of fingers.

Elastic fingering in a rotating Hele-Shaw cell.

We consider the steady-state fingering instability of an elastic membrane separating two fluids of different density under external pressure in a rotating Hele-Shaw cell. Both inextensible and highly extensible membranes are considered, and the role of membrane tension is detailed in each case. Both systems exhibit a centrifugally driven Rayleigh-Taylor-like instability when the density of the inner fluid exceeds that of the outer one, and this instability competes with the restoring forces arising from curvature and tension, thereby setting the finger scale. Numerical continuation is used to compute not only strongly nonlinear primary finger states up to the point of self-contact, but also secondary branches of mixed modes and circumferentially localized folds as a function of the rotation rate and the externally imposed pressure. Both reflection-symmetric and symmetry-broken chiral states are computed. The results are presented in the form of bifurcation diagrams. The ratio of system scale to the natural length scale is found to determine the ordering of the primary bifurcations from the unperturbed circle state as well as the solution profiles and onset of secondary bifurcations.

Viscosity ratio across interfaces controls the stability and self-assembly of microrollers

We investigate the individual and collective dynamics of torque-driven particles, called microrollers, near fluid-fluid interfaces. We find that the viscosity ratio across the interface controls the speed and direction of the particles, their relative motion, the growth of a fingering instability, and the self-assembled motile structures that emerge from it. By combining theory and large scale numerical simulations, we show how the viscosity ratio across the interface governs the long-range hydrodynamic interactions between particles and thus their collective behavior.