A scaling theory of 3D spinodal turbulence

A new scaling theory for spinodal decomposition in the inertial hydrodynamic regime is presented. The scaling involves three relevant length scales, the domain size, the Taylor microscale and the Kolmogorov dissipation scale. This allows for the presence of an inertial "energy cascade", familiar from theories of turbulence, and improves on earlier scaling treatments based on a single length: these, it is shown, cannot be reconciled with energy conservation. The new theory reconciles the t^{2/3} scaling of the domain size, predicted by simple scaling, with the physical expectation of a saturating Reynolds number at late times.Comment: 5 pages, no figures, revised version submitted to Phys Rev E Rapp Comm. Minor changes and clarification

Waves and Instabilities in Accretion Disks: MHD Spectroscopic Analysis

A complete analytical and numerical treatment of all magnetohydrodynamic waves and instabilities for radially stratified, magnetized accretion disks is presented. The instabilities are a possible source of anomalous transport. While recovering results on known hydrodynamicand both weak- and strong-field magnetohydrodynamic perturbations, the full magnetohydrodynamic spectra for a realistic accretion disk model demonstrates a much richer variety of instabilities accessible to the plasma than previously realized. We show that both weakly and strongly magnetized accretion disks are prone to strong non-axisymmetric instabilities.The ability to characterize all waves arising in accretion disks holds great promise for magnetohydrodynamic spectroscopic analysis.Comment: FOM-Institute for plasma physics "Rijnhuizen", Nieuwegein, the Netherlands 12 pages, 3 figures, Accepted for publication in ApJ

Human behavior as origin of traffic phases

It is shown that the desire for smooth and comfortable driving is directly responsible for the occurrence of complex spatio-temporal structures (``synchronized traffic'') in highway traffic. This desire goes beyond the avoidance of accidents which so far has been the main focus of microscopic modeling and which is mainly responsible for the other two phases observed empirically, free flow and wide moving jams. These features have been incorporated into a microscopic model based on stochastic cellular automata and the results of computer simulations are compared with empirical data. The simple structure of the model allows for very fast implementations of realistic networks. The level of agreement with the empirical findings opens new perspectives for reliable traffic forecasts.Comment: 4 pages, 4 figures, colour figures with reduced resolutio

Diffusion in a multi-component Lattice Boltzmann Equation model

Diffusion phenomena in a multiple component lattice Boltzmann Equation (LBE) model are discussed in detail. The mass fluxes associated with different mechanical driving forces are obtained using a Chapman-Enskog analysis. This model is found to have correct diffusion behavior and the multiple diffusion coefficients are obtained analytically. The analytical results are further confirmed by numerical simulations in a few solvable limiting cases. The LBE model is established as a useful computational tool for the simulation of mass transfer in fluid systems with external forces.Comment: To appear in Aug 1 issue of PR

Multi-component lattice-Boltzmann model with interparticle interaction

A previously proposed [X. Shan and H. Chen, Phys. Rev. E {\bf 47}, 1815, (1993)] lattice Boltzmann model for simulating fluids with multiple components and interparticle forces is described in detail. Macroscopic equations governing the motion of each component are derived by using Chapman-Enskog method. The mutual diffusivity in a binary mixture is calculated analytically and confirmed by numerical simulation. The diffusivity is generally a function of the concentrations of the two components but independent of the fluid velocity so that the diffusion is Galilean invariant. The analytically calculated shear kinematic viscosity of this model is also confirmed numerically.Comment: 18 pages, compressed and uuencoded postscript fil

3D Spinodal Decomposition in the Inertial Regime

We simulate late-stage coarsening of a 3D symmetric binary fluid using a lattice Boltzmann method. With reduced lengths and times l and t respectively (scales set by viscosity, density and surface tension) our data sets cover 1 < l 100 we find clear evidence of Furukawa's inertial scaling (l ~ t^{2/3}), although the crossover from the viscous regime (l ~ t) is very broad. Though it cannot be ruled out, we find no indication that Re is self-limiting (l ~ t^{1/2}) as proposed by M. Grant and K. R. Elder [Phys. Rev. Lett. 82, 14 (1999)].Comment: 4 pages, 3 eps figures, RevTex, minor changes to bring in line with published version. Mobility values added to Table

Lattice-Gas Simulations of Minority-Phase Domain Growth in Binary Immiscible and Ternary Amphiphilic Fluid

We investigate the growth kinetics of binary immiscible fluids and emulsions in two dimensions using a hydrodynamic lattice-gas model. We perform off-critical quenches in the binary fluid case and find that the domain size within the minority phase grows algebraically with time in accordance with theoretical predictions. In the late time regime we find a growth exponent n = 0.45 over a wide range of concentrations, in good agreement with other simluations. In the early time regime we find no universal growth exponent but a strong dependence on the concentration of the minority phase. In the ternary amphiphilic fluid case the kinetics of self assembly of the droplet phase are studied for the first time. At low surfactant concentrations, we find that, after an early algebraic growth, a nucleation regime dominates the late-time kinetics, which is enhanced by an increasing concentration of surfactant. With a further increase in the concentration of surfactant, we see a crossover to logarithmically slow growth, and finally saturation of the oil droplets, which we fit phenomenologically to a stretched exponential function. Finally, the transition between the droplet and the sponge phase is studied.Comment: 22 pages, 13 figures, submitted to PR

Three-dimensional lattice-Boltzmann simulations of critical spinodal decomposition in binary immiscible fluids

We use a modified Shan-Chen, noiseless lattice-BGK model for binary immiscible, incompressible, athermal fluids in three dimensions to simulate the coarsening of domains following a deep quench below the spinodal point from a symmetric and homogeneous mixture into a two-phase configuration. We find the average domain size growing with time as $t^\gamma$, where $\gamma$ increases in the range $0.545 < \gamma < 0.717$, consistent with a crossover between diffusive $t^{1/3}$ and hydrodynamic viscous, $t^{1.0}$, behaviour. We find good collapse onto a single scaling function, yet the domain growth exponents differ from others' works' for similar values of the unique characteristic length and time that can be constructed out of the fluid's parameters. This rebuts claims of universality for the dynamical scaling hypothesis. At early times, we also find a crossover from $q^2$ to $q^4$ in the scaled structure function, which disappears when the dynamical scaling reasonably improves at later times. This excludes noise as the cause for a $q^2$ behaviour, as proposed by others. We also observe exponential temporal growth of the structure function during the initial stages of the dynamics and for wavenumbers less than a threshold value.Comment: 45 pages, 18 figures. Accepted for publication in Physical Review

Hydrodynamic Spinodal Decomposition: Growth Kinetics and Scaling Functions

We examine the effects of hydrodynamics on the late stage kinetics in spinodal decomposition. From computer simulations of a lattice Boltzmann scheme we observe, for critical quenches, that single phase domains grow asymptotically like $t^{\alpha}$, with $\alpha \approx .66$ in two dimensions and $\alpha \approx 1.0$ in three dimensions, both in excellent agreement with theoretical predictions.Comment: 12 pages, latex, Physical Review B Rapid Communication (in press

Probability distribution of the free energy of a directed polymer in a random medium

We calculate exactly the first cumulants of the free energy of a directed polymer in a random medium for the geometry of a cylinder. By using the fact that the n-th moment of the partition function is given by the ground state energy of a quantum problem of n interacting particles on a ring of length L, we write an integral equation allowing to expand these moments in powers of the strength of the disorder gamma or in powers of n. For n small and n of order (L gamma)^(-1/2), the moments take a scaling form which allows to describe all the fluctuations of order 1/L of the free energy per unit length of the directed polymer. The distribution of these fluctuations is the same as the one found recently in the asymmetric exclusion process, indicating that it is characteristic of all the systems described by the Kardar-Parisi-Zhang equation in 1+1 dimensions.Comment: 18 pages, no figure, tu appear in PRE 61 (2000