Defect dynamics in viscous fingering

We study the dynamics of Staffman-Taylor fingering in terms of topological defects of the flow field. The defects are created and/or annihilated at the interface. The route towards the single-finger steady state is characterized by a detailed mechanism for defect annihilation. For small viscosity contrast this mechanism is impeded, and creation of new defects leads the system away from a single-finger solution. Strong evidence for a drastic reduction of the basin of attraction of the Saffman-Taylor finger is presented.

Effects of small surface tension in Hele-Shaw multifinger dynamics: an analytical and numerical study

We study the singular effects of vanishingly small surface tension on the dynamics of finger competition in the Saffman-Taylor problem, using the asymptotic techniques described by Tanveer [Philos. Trans. R. Soc. London, Ser. A 343, 155 (1993)] and Siegel and Tanveer [Phys. Rev. Lett. 76, 419 (1996)], as well as direct numerical computation, following the numerical scheme of Hou, Lowengrub, and Shelley [J. Comput. Phys. 114, 312 (1994)]. We demonstrate the dramatic effects of small surface tension on the late time evolution of two-finger configurations with respect to exact (nonsingular) zero-surface-tension solutions. The effect is present even when the relevant zero-surface-tension solution has asymptotic behavior consistent with selection theory. Such singular effects, therefore, cannot be traced back to steady state selection theory, and imply a drastic global change in the structure of phase-space flow. They can be interpreted in the framework of a recently introduced dynamical solvability scenario according to which surface tension unfolds the structurally unstable flow, restoring the hyperbolicity of multifinger fixed points.

Systematic weakly nonlinear analysis of interfacial instabilities in Hele-Shaw flows

We develop a systematic method to derive all orders of mode couplings in a weakly nonlinear approach to the dynamics of the interface between two immiscible viscous fluids in a Hele-Shaw cell. The method is completely general: it applies to arbitrary geometry and driving. Here we apply it to the channel geometry driven by gravity and pressure. The finite radius of convergence of the mode-coupling expansion is found. Calculation up to third-order couplings is done, which is necessary to account for the time-dependent Saffman-Taylor finger solution and the case of zero viscosity contrast. The explicit results provide relevant analytical information about the role that the viscosity contrast and the surface tension play in the dynamics of the system. We finally check the quantitative validity of different orders of approximation and a resummation scheme against a physically relevant, exact time-dependent solution. The agreement between the low-order approximations and the exact solution is excellent within the radius of convergence, and is even reasonably good beyond this radius.

Lateral instability in normal viscous fingers

We study a low-amplitude, long-wavelength lateral instability of the Saffman-Taylor finger by means of a phase-field model. We observe such an instability in two situations in which small dynamic perturbations are overimposed to a constant pressure drop. We first study the case in which the perturbation consists of a single oscillatory mode and then a case in which the perturbation consists of temporal noise. In both cases the instability undergoes a process of selection.

Phase-field model of Hele-Shaw flows in the high-viscosity contrast regime

A one-sided phase-field model is proposed to study the dynamics of unstable interfaces of Hele-Shaw flows in the high viscosity contrast regime. The corresponding macroscopic equations are obtained by means of an asymptotic expansion from the phase-field model. Numerical integrations of the phase-field model in a rectangular Hele-Shaw cell reproduce finger competition with the final evolution to a steady-state finger.

Experiments of periodic forcing of Saffman-Taylor fingers

We report on an experimental study of long normal Saffman-Taylor fingers subject to periodic forcing. The sides of the finger develop a low amplitude, long wavelength instability. We discuss the finger response in stationary and nonstationary situations, as well as the dynamics towards the stationary states. The response frequency of the instability increases with forcing frequency at low forcing frequencies, while, remarkably, it becomes independent of forcing frequency at large forcing frequencies. This implies a process of wavelength selection. These observations are in good agreement with previous numerical results reported in [Ledesma-Aguilar et al., Phys. Rev. E 71, 016312 (2005)]. We also study the average value of the finger width, and its fluctuations, as a function of forcing frequency. The average finger width is always smaller than the width of the steady-state finger. Fluctuations have a nonmonotonic behavior with a maximum at a particular frequency.

Experiments in a rotating Hele-Shaw cell

We introduce a modification to Hele-Shaw flows consisting of a rotating cell. A viscous fluid (oil) is injected at the rotation axis of the cell, which is open to air. The morphological instability of the oil-air interface is thus driven by centrifugal force and is controlled by the density (not viscosity) difference. We derive the linear dispersion relation and verify the maximum growth rate selection of initial patterns within experimental uncertainty. The nonlinear growth regime is studied in the case of vanishing injection rate. Several characteristic lengths are studied to quantify the patterns obtained. Experimental data exhibit good collapse for two characteristic lengths, namely, the radius of gyration and the radial finger length, which in the nonlinear regime appear to grow linearly in time.

Phase-field approach to spatial perturbations in normal Saffman-Taylor fingers

We make a numerical study of the effect that spatial perturbations have in normal Saffman-Taylor fingers driven at constant pressure gradients. We use a phase field model that allows for spatial variations in the Hele-Shaw cell. We find that, regardless of the specific way in which spatial perturbations are introduced, a lateral instability develops on the sides of the propagating Saffman-Taylor finger. Moreover, the instability exists regardless of the intensity of spatial perturbations in the cell as long as the perturbations are felt by the finger tip. If, as the finger propagates, the spatial perturbations felt by the tip change, the instability is nonperiodic. If, as the finger propagates, the spatial perturbations felt by the tip are persistent, the instability developed is periodic. In the later case, the instability is symmetrical or asymmetrical depending on the intensity of the perturbation.

Dynamical Systems approach to Saffman-Taylor fingering. A Dynamical Solvability Scenario

A dynamical systems approach to competition of Saffman-Taylor fingers in a Hele-Shaw channel is developed. This is based on global analysis of the phase space flow of the low-dimensional ordinary-differential-equation sets associated with the classes of exact solutions of the problem without surface tension. Some simple examples are studied in detail. A general proof of the existence of finite-time singularities for broad classes of solutions is given. Solutions leading to finite-time interface pinchoff are also identified. The existence of a continuum of multifinger fixed points and its dynamical implications are discussed. We conclude that exact zero-surface tension solutions taken in a global sense as families of trajectories in phase space are unphysical because the multifinger fixed points are nonhyperbolic, and an unfolding does not exist within the same class of solutions. Hyperbolicity (saddle-point structure) of the multifinger fixed points is argued to be essential to the physically correct qualitative description of finger competition. The restoring of hyperbolicity by surface tension is proposed as the key point to formulate a generic dynamical solvability scenario for interfacial pattern selection.

Surface tension and dynamics of fingering patterns

We study the minimal class of exact solutions of the Saffman-Taylor problem with zero surface tension, which contains the physical fixed points of the regularized (nonzero surface tension) problem. New fixed points are found and the basin of attraction of the Saffman-Taylor finger is determined within that class. Specific features of the physics of finger competition are identified and quantitatively defined, which are absent in the zero surface tension case. This has dramatic consequences for the long-time asymptotics, revealing a fundamental role of surface tension in the dynamics of the problem. A multifinger extension of microscopic solvability theory is proposed to elucidate the interplay between finger widths, screening and surface tension.

Fluctuations in Saffman-Taylor fingers with quenched disorder

We make an experimental characterization of the effect that static disorder has on the shape of a normal Saffman-Taylor finger. We find that static noise induces a small amplitude and long wavelength instability on the sides of the finger. Fluctuations on the finger sides have a dominant wavelength, indicating that the system acts as a selective amplifier of static noise. The dominant wavelength does not seem to be very sensitive to the intensity of static noise present in the system. On the other hand, at a given flow rate, rms fluctuations of the finger width, decrease with decreasing intensity of static noise. This might explain why the sides of the fingers are flat for typical Saffman-Taylor experiments. Comparison with previous numerical studies of the effect that temporal noise has on the Saffman-Taylor finger, leads to conclude that the effect of temporal noise and static noise are similar. The behavior of fluctuations of the finger width found in our experiments, is qualitatively similar to one recently reported, in the sense that, the magnitude of the width fluctuations decays as a power law of the capillary number, at low flow rates, and increases with capillary number for larger flow rates.

Viscous fingering as a paradigm of interfacial pattern formation: Recent results and new challenges

We review recent results on dynamical aspects of viscous fingering. The Saffman¿Taylor instability is studied beyond linear stability analysis by means of a weakly nonlinear analysis and the exact determination of the subcritical branch. A series of contributions pursuing the idea of a dynamical solvability scenario associated to surface tension in analogy with the traditional selection theory is put in perspective and discussed in the light of the asymptotic theory of Tanveer and co-workers. The inherently dynamical singular effects of surface tension are clarified. The dynamical role of viscosity contrast is explored numerically. We find that the basin of attraction of the Saffman¿Taylor finger depends on viscosity contrast, and that the sensitivity to this parameter is maximal in the usual limit of high viscosity contrast. The competing attractors are identified as closed bubble solutions. We briefly report on recent results and work in progress concerning rotating Hele-Shaw flows, topological singularities and wetting effects, and also discuss future directions in the context of viscous fingering

Three-dimensional aspects of fluid flows in channels. II. Effects of meniscus and thin film regimes on viscous fingers

We perform a three-dimensional study of steady state viscous fingers that develop in linear channels. By means of a three-dimensional lattice-Boltzmann scheme that mimics the full macroscopic equations of motion of the fluid momentum and order parameter, we study the effect of the thickness of the channel in two cases. First, for total displacement of the fluids in the channel thickness direction, we find that the steady state finger is effectively two-dimensional and that previous two-dimensional results can be recovered by taking into account the effect of a curved meniscus across the channel thickness as a contribution to surface stresses. Second, when a thin film develops in the channel thickness direction, the finger narrows with increasing channel aspect ratio in agreement with experimental results. The effect of the thin film renders the problem three-dimensional and results deviate from the two-dimensional prediction.

Three-dimensional aspects of fluid flows in channels. I. Meniscus and thin film regimes

We study the forced displacement of a fluid-fluid interface in a three-dimensional channel formed by two parallel solid plates. Using a lattice-Boltzmann method, we study situations in which a slip velocity arises from diffusion effects near the contact line. The difference between the slip and channel velocities determines whether the interface advances as a meniscus or a thin film of fluid is left adhered to the plates. We find that this effect is controlled by the capillary and Péclet numbers. We estimate the crossover from a meniscus to a thin film and find good agreement with numerical results. The penetration regime is examined in the steady state. We find that the occupation fraction of the advancing finger relative to the channel thickness is controlled by the capillary number and the viscosity contrast between the fluids. For high viscosity contrast, lattice-Boltzmann results agree with previous results. For zero viscosity contrast, we observe remarkably narrow fingers. The shape of the finger is found to be universal.

Dynamics of driven three-dimensional thin films: from hdrophilic to superhydrophobic substrates

We study the forced displacement of a thin film of fluid in contact with vertical and inclined substrates of different wetting properties, that range from hydrophilic to hydrophobic, using the lattice-Boltzmann method. We study the stability and pattern formation of the contact line in the hydrophilic and superhydrophobic regimes, which correspond to wedge-shaped and nose-shaped fronts, respectively. We find that contact lines are considerably more stable for hydrophilic substrates and small inclination angles. The qualitative behavior of the front in the linear regime remains independent of the wetting properties of the substrate as a single dispersion relation describes the stability of both wedges and noses. Nonlinear patterns show a clear dependence on wetting properties and substrate inclination angle. The effect is quantified in terms of the pattern growth rate, which vanishes for the sawtooth pattern and is finite for the finger pattern. Sawtooth shaped patterns are observed for hydrophilic substrates and low inclination angles, while finger-shaped patterns arise for hydrophobic substrates and large inclination angles. Finger dynamics show a transient in which neighboring fingers interact, followed by a steady state where each finger grows independently.