Anomalous roughening of Hele-Shaw flows with quenched disorder

The kinetic roughening of a stable oil-air interface moving in a Hele-Shaw cell that contains a quenched columnar disorder (tracks) has been studied. A capillary effect is responsible for the dynamic evolution of the resulting rough interface, which exhibits anomalous scaling. The three independent exponents needed to characterize the anomalous scaling are determined experimentally. The anomalous scaling is explained in terms of the initial acceleration and subsequent deceleration of the interface tips in the tracks coupled by mass conservation. A phenomenological model that reproduces the measured global and local exponents is introduced.

K-Rb Fermi-Bose mixtures: Vortex states and sag

We study a confined mixture of bosons and fermions in the quantal degeneracy regime with attractive boson-fermion interaction. We discuss the effect that the presence of vortical states and the displacement of the trapping potentials may have on mixtures near collapse, and investigate the phase stability diagram of the K-Rb mixture in the mean-field approximation supposing in one case that the trapping potentials felt by bosons and fermions are shifted from each other, as it happens in the presence of a gravitational sag, and in another case, assuming that the Bose condensate sustains a vortex state. In both cases, we have obtained an analytical expression for the fermion effective potential when the Bose condensate is in the Thomas-Fermi regime, that can be used to determine the maxima of the Fermionic density. We have numerically checked that the values one obtains for the location of these maxima using the analytical formulas remain valid up to the critical boson and fermion numbers, above which the mixture collapses.

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.

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.

Particle dispersion in synthetic turbulent flows

We study particle dispersion advected by a synthetic turbulent flow from a Lagrangian perspective and focus on the two-particle and cluster dispersion by the flow. It has been recently reported that Richardson¿s law for the two-particle dispersion can stem from different dispersion mechanisms, and can be dominated by either diffusive or ballistic events. The nature of the Richardson dispersion depends on the parameters of our flow and is discussed in terms of the values of a persistence parameter expressing the relative importance of the two above-mentioned mechanisms. We support this analysis by studying the distribution of interparticle distances, the relative velocity correlation functions, as well as the relative trajectories.

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.

Critical frequency for vortex nucleation in Bose-Fermi mixtures in optical lattices

We investigate within mean-field theory the influence of a one-dimensional optical lattice and of trapped degenerate fermions on the critical rotational frequency for vortex line creation in a Bose-Einstein condensate. We consider laser intensities of the lattice such that quantum coherence across the condensate is ensured. We find a sizable decrease of the thermodynamic critical frequency for vortex nucleation with increasing applied laser strength and suggest suitable parameters for experimental observation. Since 87Rb-40K mixtures may undergo collapse, we analyze the related question of how the optical lattice affects the mechanical stability of the system.

Langevin approach to generate synthetic turbulent flows

We present an analytical scheme, easily implemented numerically, to generate synthetic Gaussian turbulent flows by using a linear Langevin equation, where the noise term acts as a stochastic stirring force. The characteristic parameters of the velocity field are well introduced, in particular the kinematic viscosity and the spectrum of energy. As an application, the diffusion of a passive scalar is studied for two different energy spectra. Numerical results are compared favorably with analytical calculations.

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.

Elasticity and melting of vortex crystals in anisotropic superconductors: Beyond the continuum regime.

The elastic moduli of vortex crystals in anisotropic superconductors are frequently involved in the investigation of their phase diagram and transport properties. We provide a detailed analysis of the harmonic eigenvalues (normal modes) of the vortex lattice for general values of the magnetic field strength, going beyond the elastic continuum regime. The detailed behavior of these wave-vector-dependent eigenvalues within the Brillouin zone (BZ), is compared with several frequently used approximations that we also recalculate. Throughout the BZ, transverse modes are less costly than their longitudinal counterparts, and there is an angular dependence which becomes more marked close to the zone boundary. Based on these results, we propose an analytic correction to the nonlocal continuum formulas which fits quite well the numerical behavior of the eigenvalues in the London regime. We use this approximate expression to calculate thermal fluctuations and the full melting line (according to Lindeman's criterion) for various values of the anisotropy parameter.

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.