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.

Slow Dynamics of Ising Models with Energy Barriers

Using Monte Carlo simulations we study the dynamics of three-dimensional Ising models with nearest-, next-nearest-, and four-spin (plaquette) interactions. During coarsening, such models develop growing energy barriers, which leads to very slow dynamics at low temperature. As already reported, the model with only the plaquette interaction exhibits some of the features characteristic of ordinary glasses: strong metastability of the supercooled liquid, a weak increase of the characteristic length under cooling, stretched-exponential relaxation, and aging. The addition of two-spin interactions, in general, destroys such behavior: the liquid phase loses metastability and the slow-dynamics regime terminates well below the melting transition, which is presumably related with a certain corner-rounding transition. However, for a particular choice of interaction constants, when the ground state is strongly degenerate, our simulations suggest that the slow-dynamics regime extends up to the melting transition. The analysis of these models leads us to the conjecture that in the four-spin Ising model domain walls lose their tension at the glassy transition and that they are basically tensionless in the glassy phase.

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.

Reaction-diffusion fronts under stochastic advection

We study front propagation in stirred media using a simplified modelization of the turbulent flow. Computer simulations reveal the existence of the two limiting propagation modes observed in recent experiments with liquid phase isothermal reactions. These two modes respectively correspond to a wrinkled although sharp propagating interface and to a broadened one. Specific laws relative to the enhancement of the front velocity in each regime are confirmed by our simulations.

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.

On the many-time formulation of classical particle dynamics

Starting from the standard one-time dynamics of n nonrelativistic particles, the n-time equations of motion are inferred, and a variational principle is formulated. A suitable generalization of the classical LieKnig theorem is demonstrated, which allows the determination of all the associated presymplectic structures. The conditions under which the action of an invariance group is canonical are studied, and a corresponding Noether theorem is deduced. A formulation of the theory in terms of n first-class constraints is recovered by means of coisotropic imbeddings. The proposed approach also provides for a better understanding of the relativistic particle dynamics, since it shows that the different roles of the physical positions and the canonical variables is not peculiar to special relativity, but rather to any n-time approach: indeed a nonrelativistic no-interaction theorem is deduced.

Higher order Lagrangian systems: Geometric structures, Dynamics, and Constraints

In order to study the connections between Lagrangian and Hamiltonian formalisms constructed from aperhaps singularhigher-order Lagrangian, some geometric structures are constructed. Intermediate spaces between those of Lagrangian and Hamiltonian formalisms, partial Ostrogradskiis transformations and unambiguous evolution operators connecting these spaces are intrinsically defined, and some of their properties studied. Equations of motion, constraints, and arbitrary functions of Lagrangian and Hamiltonian formalisms are thoroughly studied. In particular, all the Lagrangian constraints are obtained from the Hamiltonian ones. Once the gauge transformations are taken into account, the true number of degrees of freedom is obtained, both in the Lagrangian and Hamiltonian formalisms, and also in all the intermediate formalisms herein defined.

Front dynamics in turbulent media

A study of a stable front propagating in a turbulent medium is presented. The front is generated through a reaction-diffusion equation, and the turbulent medium is statistically modeled using a Langevin equation. Numerical simulations indicate the presence of two different dynamical regimes. These regimes appear when the turbulent flow either wrinkles a still rather sharp propagating interfase or broadens it. Specific dependences of the propagating velocities on stirring intensities appropriate to each case are found and fitted when possible according to theoretically predicted laws. Different turbulent spectra are considered.

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.

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.

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.

On the determination of the interface density of states in a-Si:H/a-SiC:H multilayers

This paper deals with the determination of the interface density of states in amorphous silicon-based multilayers. Photothermal deflection spectroscopy is used to characterize two series of aSi:H/aSi1-xCx:H multilayers, and a new approach in the treatment of experimental dada is used in order to obtain accurate results. From this approach, an upper limit of 10^10 cm-2 is determined for the interface density of states.

Statistical dynamics of dislocation systems: The influence of dislocation-dislocation correlations.

During plastic deformation of crystalline materials, the collective dynamics of interacting dislocations gives rise to various patterning phenomena. A crucial and still open question is whether the long range dislocation-dislocation interactions which do not have an intrinsic range can lead to spatial patterns which may exhibit well-defined characteristic scales. It is demonstrated for a general model of two-dimensional dislocation systems that spontaneously emerging dislocation pair correlations introduce a length scale which is proportional to the mean dislocation spacing. General properties of the pair correlation functions are derived, and explicit calculations are performed for a simple special case, viz pair correlations in single-glide dislocation dynamics. It is shown that in this case the dislocation system exhibits a patterning instability leading to the formation of walls normal to the glide plane. The results are discussed in terms of their general implications for dislocation patterning.