Viscous fingering instability in one-end-lifted Hele-Shaw cells for producing three-dimensional hierarchical structures.

Impact of an intermediate layer on immiscible viscous fingering instability in radial Hele-Shaw cell

Viscous fingering instabilities, common in confined environments such as porous media or Hele-Shaw cells, surprisingly also occur in unconfined, non-porous settings as revealed by recent experiments. These novel instabilities involve free-surface flows of dissimilar viscosity. We demonstrate that such a free-surface flow, involving a thin film of viscous fluid spreading over a substrate that is prewetted with a fluid of higher viscosity, is susceptible to a similar type of novel viscous fingering instability. Such flows are relevant to a range of geophysical, industrial and physiological applications from the small scales of thin-film coating applications and nasal drug delivery to the large scales of lava flows. In developing a theoretical framework, we assume that the intruding layer and the liquid film over which it flows are both long and thin, the effects of inertia and surface tension are negligible, and both layers are driven by gravity and resisted by viscous shear stress so that the principles of lubrication theory hold. We investigate the stability of axisymmetric similarity solutions, describing the base flow, by examining the growth of small-amplitude non-axisymmetric perturbations. We characterise regions of instability across parameter space and find that these instabilities emerge above a critical viscosity ratio. That is, a fluid of low viscosity intruding into another fluid of sufficiently high viscosity is susceptible to instability, akin to traditional viscous fingering in a porous medium. We identify the mechanism of instability, compare with other frontal instabilities and demonstrate that high enough density differences suppress the instability completely.

Viscous fingering and interface splitting instabilities in air-water-oil systems.

Viscous fingering instability has been analyzed through empirical studies using miscible flow displacement in fractured porous media. While significant research has been conducted on viscous fingering, limited information is available regarding its behavior in fractured porous structures. The experiments were conducted in rectangular porous models with fractures oriented at 0∘, 45∘, and 90∘, to investigate how fracture orientation influences fluid displacement, where both channeling and fingering mechanisms play significant roles. This paper, which is the second part of a previous study, places particular emphasis on the impact of the viscosity ratio, a crucial parameter in determining the complexity of the fingering patterns. Quantitative parameters such as sweep efficiency, tip location, and breakthrough time were evaluated and analyzed using image processing techniques. The results indicate that increasing the viscosity ratio leads to more complex finger formations. Additionally, as the injection rate increases, the size of the finger patterns slightly increases, while the channeling effect becomes less pronounced. Notably, fractures aligned at 0∘ had the most significant impact on the rate of sweep efficiency and tip location, increasing the tip velocity of the fingers by up to 90%.

A theoretical analysis of viscous fingering instability for a reactive system A + B → C with an infinitely fast reaction in a porous medium for a rectilinear flow is presented. By contrast to the traditional quasi-steady-state analysis (QSSA), a non-modal analysis (NMA) based on the fundamental matrix formulation is employed to study the reactive displacement, considering reactants and products with mismatched viscosities. This study investigates the transient growth of perturbations by analysing the singular values and singular vectors to address the optimal energy amplification. We illustrate that an increase in the viscosity contrast, | R b − R c | , resulting from a chemical reaction for a given endpoint viscosity contrast R b , leads to a more unstable system. However, there exist some reactions when R c > R b , the onset delays than the equivalent non-reactive case, R c = R b . It suggests that the stability of the flow is primarily influenced when instability develops downstream within the flow. Furthermore, R b is found to significantly affect the spatio-temporal evolution of perturbations and the underlying physical mechanism. It is demonstrated that the QSSA is inadequate to address the transient growth, and NMA is the most suitable approach to studying the underlying physical mechanism of instability. Furthermore, NMA results align more consistently with non-linear simulations compared with QSSA.

We experimentally demonstrate that the coacervation of a biopolymer can trigger a hydrodynamic instability when a coacervate is formed upon injection of a xanthan gum dispersion into a cationic surfactant (C14TAB) solution. The local increase of the viscosity due to the coacervate formation induces a viscous fingering instability. Three characteristic displacement regimes were observed: a viscous fingering dominated regime, a buoyancy-controlled "volcano" regime and a "fan"-like regime determined by the coacervate membrane dynamics. The dependence of the spatial properties of the viscous fingering pattern on the Péclet and Rayleigh numbers is investigated.

We have studied viscous fingering instability patterns that frequently develop during thin-film detachment. An experimental investigation was conducted employing a thin layer of a simple yield stress fluid between two rigid circular parallel plates of a commercial rheometer operating in a probe-tack mode. We present a method to delay this type of interface instability by allowing an extra amount, a small external cap, of the mediating working fluid bulging out the contour boundary of the plates. When carefully prepared, this extra cap moves the air-fluid interface away from the high velocity flow field generated near the confines of the plates at the beginning of the detachment. A lower velocity flow field causes a delay on the instability growth which, within a given set of parameters, reduces or completely removes the fingering pattern. The effect of fingering instability in the presence of extra amounts of extra sample were investigated by monitoring the traction normal-force versus gap separation during detachment. Studies varying both the initial gap, the gap separation velocity, and the amount of a carefully prepared externally positioned sample layer were systematically conducted in order to determine the conditions where the instability delay was sufficient to remove fingering patterns. The force-decaying regimes were identified as well as their dependencies on the experiment control parameters. A proper adjustment of the extra cap resulted in complete removal of the fingering pattern while the peak of the normal traction force was unaffected. The absence of fingering both during and after the detachment was confirmed by videoing the process from underneath with a transparent glass bottom plate. Nevertheless, when the traction force profile was compared to the standard results obtained with trimmed samples, an additional stage with a levelled-off force interval, akin to those caused by cavitation, took place. Finally, we show the method is robust and may contribute as an alternative approach for delaying, and sometimes completely removing, the effect of fluid fingering on fluid-surface contact applications where one desires a final uniform sample distribution.

When a less-viscous solution of a reactant $A$ displaces a more-viscous solution of another reactant $B$ , a fast bimolecular $A + B \rightarrow C$ reaction decreasing locally the viscosity can influence the viscous fingering (VF) instability taking place between the two miscible solutions. We show both experimentally and numerically that, for monotonic viscosity profiles, this decrease in viscosity has opposite effects on the fingering pattern depending on the injection flow rate. For high flow rates, the reaction enhances the shielding effect, creating VF patterns with a lower surface density, i.e. thinner fingers covering a smaller area. In contrast, for lower flow rates, the reaction stabilises the VF dynamics, i.e. delays the instability and gives a less-deformed displacement, reaching in some cases an almost-stable displacement. Nonlinear simulations of reactive VF show that these opposite effects at low or high flow rates can only be reproduced if the diffusivity of $A$ is larger than that of $B$ , which favours a larger production of the product $C$ and, hence, a larger viscosity decrease. The analysis of one-dimensional viscosity profiles reconstructed on the basis of a one-dimensional reaction–diffusion–advection model confirms that the VF stabilisation at low Péclet number and in the presence of differential diffusion of reactants originates from an optimum reaction-driven decrease in the gradient of the monotonic viscosity profile.

We investigate the onset and growth of viscous fingering (VF) of miscible annulus in a radial Hele-Shaw cell. Systematic numerical study on a finite annulus domain is performed by employing finite element method solver in COMSOL Multiphysics software. We justify that concentration field analysis is not a good choice for dynamic study in radial flows. Instead, velocity magnitude is a better tool to understand the dynamics. Therefore, we propose velocity field analysis to better differentiate the stable and unstable states and present a new stability criterion using the velocity field method. Most interestingly, using the velocity field analysis and the new stability criterion, we show a restabilization of the VF at a critical time when the system becomes diffusion dominant and able to provide both the onset time, τon (time at which instability develops), and the time at which the interface returns to the stable state, τd. Furthermore, the study successfully suggests the critical values for several dimensionless parameters, the Péclet number (Pe), log-viscosity ratio (R), and volumetric ratio (ra) and time (τ), to induce instability. When Pe is higher than 103, the evolution of VF instability is no longer enhanced by Pe, and Rc converges to a certain value. In particular, for the transiently unstable system of low Pe, the restabilization of VF instability is identified even though R is higher than Rc. The unstable system (τ>τon) returns to the stable state as injection time increases further. Moreover, we obtained a critical value of the volumetric ratio (rc,a).

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.

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

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.

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.

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.

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.

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.

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 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.

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.