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.

Statistical mechanics in the extended Gaussian ensemble

The extended Gaussian ensemble (EGE) is introduced as a generalization of the canonical ensemble. This ensemble is a further extension of the Gaussian ensemble introduced by Hetherington [J. Low Temp. Phys. 66, 145 (1987)]. The statistical mechanical formalism is derived both from the analysis of the system attached to a finite reservoir and from the maximum statistical entropy principle. The probability of each microstate depends on two parameters ß and ¿ which allow one to fix, independently, the mean energy of the system and the energy fluctuations, respectively. We establish the Legendre transform structure for the generalized thermodynamic potential and propose a stability criterion. We also compare the EGE probability distribution with the q-exponential distribution. As an example, an application to a system with few independent spins is presented.

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.

Ground-state properties of a dilute homogeneous Bose gas of hard disks in two dimensions

The energy and structure of a dilute hard-disks Bose gas are studied in the framework of a variational many-body approach based on a Jastrow correlated ground-state wave function. The asymptotic behaviors of the radial distribution function and the one-body density matrix are analyzed after solving the Euler equation obtained by a free minimization of the hypernetted chain energy functional. Our results show important deviations from those of the available low density expansions, already at gas parameter values x~0.001 . The condensate fraction in 2D is also computed and found generally lower than the 3D one at the same x.

Time-dependent couplings and crossover length scales in nonequilibrium surface roughening

We show that time-dependent couplings may lead to nontrivial scaling properties of the surface fluctuations of the asymptotic regime in nonequilibrium kinetic roughening models. Three typical situations are studied. In the case of a crossover between two different rough regimes, the time-dependent coupling may result in anomalous scaling for scales above the crossover length. In a different setting, for a crossover from a rough to either a flat or damping regime, the time-dependent crossover length may conspire to produce a rough surface, although the most relevant term tends to flatten the surface. In addition, our analysis sheds light into an existing debate in the problem of spontaneous imbibition, where time-dependent couplings naturally arise in theoretical models and experiments.

Influence of disorder strength on phase-field models of interfacial growth

We study the influence of disorder strength on the interface roughening process in a phase-field model with locally conserved dynamics. We consider two cases where the mobility coefficient multiplying the locally conserved current is either constant throughout the system (the two-sided model) or becomes zero in the phase into which the interface advances (one-sided model). In the limit of weak disorder, both models are completely equivalent and can reproduce the physical process of a fluid diffusively invading a porous media, where super-rough scaling of the interface fluctuations occurs. On the other hand, increasing disorder causes the scaling properties to change to intrinsic anomalous scaling. In the limit of strong disorder this behavior prevails for the one-sided model, whereas for the two-sided case, nucleation of domains in front of the invading front are observed.

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.

Brownian motion in short range random potentials

A numerical study of Brownian motion of noninteracting particles in random potentials is presented. The dynamics are modeled by Langevin equations in the high friction limit. The random potentials are Gaussian distributed and short ranged. The simulations are performed in one and two dimensions. Different dynamical regimes are found and explained. Effective subdiffusive exponents are obtained and commented on.

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.

Side-branch growth in two-dimensional dendrits. Part II: Phase-field model

The development of side-branching in solidifying dendrites in a regime of large values of the Peclet number is studied by means of a phase-field model. We have compared our numerical results with experiments of the preceding paper and we obtain good qualitative agreement. The growth rate of each side branch shows a power-law behavior from the early stages of its life. From their birth, branches which finally succeed in the competition process of side-branching development have a greater growth exponent than branches which are stopped. Coarsening of branches is entirely defined by their geometrical position relative to their dominant neighbors. The winner branches escape from the diffusive field of the main dendrite and become independent dendrites.

Brownian pump powered by a white-noise flasing ratchet

A Brownian pump of particles powered by a stochastic flashing ratchet mechanism is studied. The pumping device is embedded in a finite region and bounded by particle reservoirs. In the steady state, we exactly calculate the spatial density profile, the concentration ratio between both reservoirs and the particle flux. We propose a simulation framework for the consistent evaluation of such observable quantities.

Distorted-wave calculation of cross sections for inner-shell ionization by electron and positron impact

The relativistic distorted-wave Born approximation is used to calculate differential and total cross sections for inner shell ionization of neutral atoms by electron and positron impact. The target atom is described within the independent-electron approximation using the self-consistent Dirac-Fock-Slater potential. The distorting potential for the projectile is also set equal to the Dirac-Fock-Slater potential. For electrons, this guarantees orthogonality of all the orbitals involved and simplifies the calculation of exchange T-matrix elements. The interaction between the projectile and the target electrons is assumed to reduce to the instantaneous Coulomb interaction. The adopted numerical algorithm allows the calculation of differential and total cross sections for projectiles with kinetic energies ranging from the ionization threshold up to about ten times this value. Algorithm accuracy and stability are demonstrated by comparing differential cross sections calculated by our code with the distorting potential set to zero with equivalent results generated by a more robust code that uses the conventional plane-wave Born approximation. Sample calculation results are presented for ionization of K- and L-shells of various elements and compared with the available experimental data.