Comment on "Selection of the Saffman-Taylor finger width in the absence of surface tension: an exact result"

A Comment on the Letter by Mark Mineev-Weinstein, Phys. Rev. Lett. 80, 2113 (1998). The authors of the Letter offer a Reply.

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.

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.

Measurement of a Structured Backflow in an Open Small Channel Induced by Surface-Tension Gradients

We present experiments in which the laterally confined flow of a surfactant film driven by controlled surface tension gradients causes the subtended liquid layer to self-organize into an inner upstream microduct surrounded by the downstream flow. The anomalous interfacial flow profiles and the concomitant backflow are a result of the feedback between two-dimensional and three-dimensional microfluidics realized during flow in open microchannels. Bulk and surface particle image velocimetry data combined with an interfacial hydrodynamics model explain the dependence of the observed phenomena on channel geometry.

Pressure-dependent scaling scenarios in experiments of spontaneous imbibition

The scaling properties of the rough liquid-air interface formed in the spontaneous imbibition of a viscous liquid by a model porous medium are found to be very sensitive to the magnitude of the pressure difference applied at the liquid inlet. Interface fluctuations change from obeying intrinsic anomalous scaling at large negative pressure differences, to being super-rough with the same dynamic exponent z¿3 at less negative pressure differences, to finally obeying ordinary Family-Vicsek scaling with z¿2 at large positive pressure differences. This rich scenario reflects the relative importance on different length scales of capillary and permeability disorder, and the role of surface tension and viscous pressure in damping interface fluctuations.

Cavitation in 3He-4He liquid mixtures at low temperatures

Bubble formation in solutions of 3He and 4He is studied within a density-functional approach. In particular, the temperature dependence of the cavitation pressure for different 3He concentrations is calculated at low temperatures and compared to that of pure 4He. The presence of Andreev states lowers the surface tension and, consequently, nucleation barriers are drastically reduced. This fact means that even at low 3He concentrations the cavitation process takes place at higher pressures than the spinodal pressure, which is not the case for pure 4He.

Interfacial growth in driven diffusive systems

We have studied the growth of interfaces in driven diffusive systems well below the critical temperature by means of Monte Carlo simulations. We consider the region beyond the linear regime and of large values of the external field which has not been explored before. The simulations support the existence of interfacial traveling waves when asymmetry is introduced in the model, a result previously predicted by a linear-stability analysis. Furthermore, the generalization of the Gibbs-Thomson relation is discussed. The results provide evidence that the external field is a stabilizing effect which can be considered as effectively increasing the surface tension.

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.

Dynamical scaling of imbibition in columnar geometries

compatible with the usual nonlocal model governed by surface tension that results from a macroscopic description. To explore this discrepancy, we exhaustively analyze numerical integrations of a phase-field model with dichotomic columnar disorder. We find that two distinct behaviors are possible depending on the capillary contrast between the two values of disorder. In a high-contrast case, where interface evolution is mainly dominated by the disorder, an inherent anomalous scaling is always observed. Moreover, in agreement with experimental work, the interface motion has to be described through a local model. On the other hand, in a lower-contrast case, the interface is dominated by interfacial tension and can be well modeled by a nonlocal model. We have studied both spontaneous and forced-flow imbibition situations, giving a complete set of scaling exponents in each case, as well as a comparison to the experimental results.

Dynamics of the two-dimensional gonihedric spin model

In this paper, we study dynamical aspects of the two-dimensional (2D) gonihedric spin model using both numerical and analytical methods. This spin model has vanishing microscopic surface tension and it actually describes an ensemble of loops living on a 2D surface. The self-avoidance of loops is parametrized by a parameter ¿. The ¿=0 model can be mapped to one of the six-vertex models discussed by Baxter, and it does not have critical behavior. We have found that allowing for ¿¿0 does not lead to critical behavior either. Finite-size effects are rather severe, and in order to understand these effects, a finite-volume calculation for non-self-avoiding loops is presented. This model, like his 3D counterpart, exhibits very slow dynamics, but a careful analysis of dynamical observables reveals nonglassy evolution (unlike its 3D counterpart). We find, also in this ¿=0 case, the law that governs the long-time, low-temperature evolution of the system, through a dual description in terms of defects. A power, rather than logarithmic, law for the approach to equilibrium has been found.

Viscous fingering in liquid crystals: Anisotropy and morphological transitions

We show that a minimal model for viscous fingering with a nematic liquid crystal in which anisotropy is considered to enter through two different viscosities in two perpendicular directions can be mapped to a twofold anisotropy in the surface tension. We numerically integrate the dynamics of the resulting problem with the phase-field approach to find and characterize a transition between tip splitting and side branching as a function of both anisotropy and dimensionless surface tension. This anisotropy dependence could explain the experimentally observed (reentrant) transition as temperature and applied pressure are varied. Our observations are also consistent with previous experimental evidence in viscous fingering within an etched cell and simulations of solidification.

Interface dynamics in Hele-Shaw flows with centrifugal forces. Preventing cusp singularities with rotation

A class of exact solutions of Hele-Shaw flows without surface tension in a rotating cell is reported. We show that the interplay between injection and rotation modifies the scenario of formation of finite-time cusp singularities. For a subclass of solutions, we show that, for any given initial condition, there exists a critical rotation rate above which cusp formation is suppressed. We also find an exact sufficient condition to avoid cusps simultaneously for all initial conditions within the above subclass.

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.

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

Grain boundaries in vortex matter

We explore the statistical properties of grain boundaries in the vortex polycrystalline phase of type-II superconductors. Treating grain boundaries as arrays of dislocations interacting through linear elasticity, we show that self-interaction of a deformed grain boundary is equivalent to a nonlocal long-range surface tension. This affects the pinning properties of grain boundaries, which are found to be less rough than isolated dislocations. The presence of grain boundaries has an important effect on the transport properties of type-II superconductors as we show by numerical simulations: our results indicate that the critical current is higher for a vortex polycrystal than for a regular vortex lattice. Finally, we discuss the possible role of grain boundaries in vortex lattice melting. Through a phenomenological theory we show that melting can be preceded by an intermediate polycrystalline phase.