Phase-field model for Hele-Shaw flows with arabitrary contrast II: numerical approach

We present a phase-field model for the dynamics of the interface between two inmiscible fluids with arbitrary viscosity contrast in a rectangular Hele-Shaw cell. With asymptotic matching techniques we check the model to yield the right Hele-Shaw equations in the sharp-interface limit, and compute the corrections to these equations to first order in the interface thickness. We also compute the effect of such corrections on the linear dispersion relation of the planar interface. We discuss in detail the conditions on the interface thickness to control the accuracy and convergence of the phase-field model to the limiting Hele-Shaw dynamics. In particular, the convergence appears to be slower for high viscosity contrasts.

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.

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.

Ballistic and diffusive corrections to front propagation in the presence of multiplicative noise

We study the dynamics of reaction-diffusion fronts under the influence of multiplicative noise. An approximate theoretical scheme is introduced to compute the velocity of the front and its diffusive wandering due to the presence of noise. The theoretical approach is based on a multiple scale analysis rather than on a small noise expansion and is confirmed with numerical simulations for a wide range of the noise intensity. We report on the possibility of noise sustained solutions with a continuum of possible velocities, in situations where only a single velocity is allowed without noise.

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

Coarsening of solid-liquid mixtures in a random acceleration field

The effects of flow induced by a random acceleration field (g-jitter) are considered in two related situations that are of interest for microgravity fluid experiments: the random motion of isolated buoyant particles, and diffusion driven coarsening of a solid-liquid mixture. We start by analyzing in detail actual accelerometer data gathered during a recent microgravity mission, and obtain the values of the parameters defining a previously introduced stochastic model of this acceleration field. The diffusive motion of a single solid particle suspended in an incompressible fluid that is subjected to such random accelerations is considered, and mean squared velocities and effective diffusion coefficients are explicitly given. We next study the flow induced by an ensemble of such particles, and show the existence of a hydrodynamically induced attraction between pairs of particles at distances large compared with their radii, and repulsion at short distances. Finally, a mean field analysis is used to estimate the effect of g-jitter on diffusion controlled coarsening of a solid-liquid mixture. Corrections to classical coarsening rates due to the induced fluid motion are calculated, and estimates are given for coarsening of Sn-rich particles in a Sn-Pb eutectic fluid, an experiment to be conducted in microgravity in the near future.

Low viscosity contrast fingering in a rotating Hele-Shaw cell

We study the fingering instability of a circular interface between two immiscible liquids in a radial Hele-Shaw cell. The cell rotates around its vertical symmetry axis, and the instability is driven by the density difference between the two fluids. This kind of driving allows studying the interfacial dynamics in the particularly interesting case of an interface separating two liquids of comparable viscosity. An accurate experimental study of the number of fingers emerging from the instability reveals a slight but systematic dependence of the linear dispersion relation on the gap spacing. We show that this result is related to a modification of the interface boundary condition which incorporates stresses originated from normal velocity gradients. The early nonlinear regime shows nearly no competition between the outgrowing fingers, characteristic of low viscosity contrast flows. We perform experiments in a wide range of experimental parameters, under conditions of mass conservation (no injection), and characterize the resulting patterns by data collapses of two characteristic lengths: the radius of gyration of the pattern and the interface stretching. Deep in the nonlinear regime, the fingers which grow radially outwards stretch and become gradually thinner, to a point that the fingers pinch and emit drops. We show that the amount of liquid emitted in the first generation of drops is a constant independent of the experimental parameters. Further on there is a sharp reduction of the amount of liquid centrifugated, punctuated by periods of no observable centrifugation.

Bifurcations and chaos in single-roll natural convection with low Prandtl number

Convective flows of a small Prandtl number fluid contained in a two-dimensional cavity subject to a lateral thermal gradient are numerically studied by using different techniques. The aspect ratio (length to height) is kept at around 2. This value is found optimal to make the flow most unstable while keeping the basic single-roll structure. Two cases of thermal boundary conditions on the horizontal plates are considered: perfectly conducting and adiabatic. For increasing Rayleigh numbers we find a transition from steady flow to periodic oscillations through a supercritical Hopf bifurcation that maintains the centrosymmetry of the basic circulation. For a Rayleigh number of about ten times that of the Hopf bifurcation the system initiates a complex scenario of bifurcations. In the conductive case these include a quasiperiodic route to chaos. In the adiabatic one the dynamics is dominated by the interaction of two Neimark-Sacker bifurcations of the basic periodic solutions, leading to the stable coexistence of three incommensurate frequencies, and finally to chaos. In all cases, the complex time-dependent behavior does not break the basic, single-roll structure.

Theory for correlation functions of processes driven by external colored noise

We study steady-state correlation functions of nonlinear stochastic processes driven by external colored noise. We present a methodology that provides explicit expressions of correlation functions approximating simultaneously short- and long-time regimes. The non-Markov nature is reduced to an effective Markovian formulation, and the nonlinearities are treated systematically by means of double expansions in high and low frequencies. We also derive some exact expressions for the coefficients of these expansions for arbitrary noise by means of a generalization of projection-operator techniques.

Intensity correlation functions of dye lasers: Comparison of colored-gain-noise and colored-loss-noise models

The intensity correlation functions C(t) for the colored-gain-noise model of dye lasers are analyzed and compared with those for the loss-noise model. For correlation times ¿ larger than the deterministic relaxation time td, we show with the use of the adiabatic approximation that C(t) values coincide for both models. For small correlation times we use a method that provides explicit expressions of non-Markovian correlation functions, approximating simultaneously short- and long-time behaviors. Comparison with numerical simulations shows excellent results simultaneously for short- and long-time regimes. It is found that, when the correlation time of the noise increases, differences between the gain- and loss-noise models tend to disappear. The decay of C(t) for both models can be described by a time scale that approaches the deterministic relaxation time. However, in contrast with the loss-noise model, a secondary time scale remains for large times for the gain-noise model, which could allow one to distinguish between both models.

Perturbing Hele-Shaw flow with a small gap gradient

A controlled perturbation is introduced into the Saffman-Taylor flow problem by adding a gradient to the gap of a Hele-Shaw cell. The stability of the single-finger steady state was found to be strongly affected by such a perturbation. Compared with patterns in a standard Hele-Shaw cell, the single Saffman-Taylor finger was stabilized or destabilized according to the sign of the gap gradient. While a linear stability analysis shows that this perturbation should have a negligible effect on the early-stage pattern formation, the experimental data indicate that the characteristic length for the initial breakup of a flat interface has been changed by the perturbation.