Onset of Patterns in an Ocillated Granular Layer: Continuum and Molecular Dynamics Simulations

We study the onset of patterns in vertically oscillated layers of frictionless dissipative particles. Using both numerical solutions of continuum equations to Navier-Stokes order and molecular dynamics (MD) simulations, we find that standing waves form stripe patterns above a critical acceleration of the cell. Changing the frequency of oscillation of the cell changes the wavelength of the resulting pattern; MD and continuum simulations both yield wavelengths in accord with previous experimental results. The value of the critical acceleration for ordered standing waves is approximately 10% higher in molecular dynamics simulations than in the continuum simulations, and the amplitude of the waves differs significantly between the models. The delay in the onset of order in molecular dynamics simulations and the amplitude of noise below this onset are consistent with the presence of fluctuations which are absent in the continuum theory. The strength of the noise obtained by fit to Swift-Hohenberg theory is orders of magnitude larger than the thermal noise in fluid convection experiments, and is comparable to the noise found in experiments with oscillated granular layers and in recent fluid experiments on fluids near the critical point. Good agreement is found between the mean field value of onset from the Swift-Hohenberg fit and the onset in continuum simulations. Patterns are compared in cells oscillated at two different frequencies in MD; the layer with larger wavelength patterns has less noise than the layer with smaller wavelength patterns.Comment: Published in Physical Review

Linear instability and statistical laws of physics

We show that a meaningful statistical description is possible in conservative and mixing systems with zero Lyapunov exponent in which the dynamical instability is only linear in time. More specifically, (i) the sensitivity to initial conditions is given by $\xi =[1+(1-q)\lambda_q t]^{1/(1-q)}$ with $q=0$; (ii) the statistical entropy $S_q=(1-\sum_i p_i^q)/(q-1) (S_1=-\sum_i p_i \ln p_i)$ in the infinitely fine graining limit (i.e., $W\equiv$ {\it number of cells into which the phase space has been partitioned} $\to\infty$), increases linearly with time only for $q=0$; (iii) a nontrivial, $q$-generalized, Pesin-like identity is satisfied, namely the $\lim_{t \to \infty} \lim_{W \to \infty} S_0(t)/t=\max\{\lambda_0\}$. These facts (which are in analogy to the usual behaviour of strongly chaotic systems with $q=1$), seem to open the door for a statistical description of conservative many-body nonlinear systems whose Lyapunov spectrum vanishes.Comment: 7 pages including 2 figures. The present version is accepted for publication in Europhysics Letter

Measuring non-extensitivity parameters in a turbulent Couette-Taylor flow

We investigate probability density functions of velocity differences at different distances r measured in a Couette-Taylor flow for a range of Reynolds numbers Re. There is good agreement with the predictions of a theoretical model based on non-extensive statistical mechanics (where the entropies are non-additive for independent subsystems). We extract the scale-dependent non-extensitivity parameter q(r, Re) from the laboratory data.Comment: 8 pages, 5 figure

Fluctuations in viscous fingering

Our experiments on viscous (Saffman-Taylor) fingering in Hele-Shaw channels reveal finger width fluctuations that were not observed in previous experiments, which had lower aspect ratios and higher capillary numbers Ca. These fluctuations intermittently narrow the finger from its expected width. The magnitude of these fluctuations is described by a power law, Ca^{-0.64}, which holds for all aspect ratios studied up to the onset of tip instabilities. Further, for large aspect ratios, the mean finger width exhibits a maximum as Ca is decreased instead of the predicted monotonic increase.Comment: Revised introduction, smoothed transitions in paper body, and added a few additional minor results. (Figures unchanged.) 4 pages, 3 figures. Submitted to PRE Rapi

Shocks in supersonic sand

We measure time-averaged velocity, density, and temperature fields for steady granular flow past a wedge and calculate a speed of granular pressure disturbances (sound speed) equal to 10% of the flow speed. The flow is supersonic, forming shocks nearly identical to those in a supersonic gas. Molecular dynamics simulations of Newton's laws and Monte Carlo simulations of the Boltzmann equation yield fields in quantitative agreement with experiment. A numerical solution of Navier-Stokes-like equations agrees with a molecular dynamics simulation for experimental conditions excluding wall friction.Comment: 4 pages, 5 figure

Probing structural relaxation in complex fluids by critical fluctuations

Complex fluids, such as polymer solutions and blends, colloids and gels, are of growing interest in fundamental and applied soft-condensed-matter science. A common feature of all such systems is the presence of a mesoscopic structural length scale intermediate between atomic and macroscopic scales. This mesoscopic structure of complex fluids is often fragile and sensitive to external perturbations. Complex fluids are frequently viscoelastic (showing a combination of viscous and elastic behaviour) with their dynamic response depending on the time and length scales. Recently, non-invasive methods to infer the rheological response of complex fluids have gained popularity through the technique of microrheology, where the diffusion of probe spheres in a viscoelastic fluid is monitored with the aid of light scattering or microscopy. Here we propose an alternative to traditional microrheology that does not require doping of probe particles in the fluid (which can sometimes drastically alter the molecular environment). Instead, our proposed method makes use of the phenomenon of "avoided crossing" between modes associated with the structural relaxation and critical fluctuations that are spontaneously generated in the system.Comment: 4 pages, 4 figure

Transport Coefficients for Granular Media from Molecular Dynamics Simulations

Under many conditions, macroscopic grains flow like a fluid; kinetic theory pred icts continuum equations of motion for this granular fluid. In order to test the theory, we perform event driven molecular simulations of a two-dimensional gas of inelastic hard disks, driven by contact with a heat bath. Even for strong dissipation, high densities, and small numbers of particles, we find that continuum theory describes the system well. With a bath that heats the gas homogeneously, strong velocity correlations produce a slightly smaller energy loss due to inelastic collisions than that predicted by kinetic theory. With an inhomogeneous heat bath, thermal or velocity gradients are induced. Determination of the resulting fluxes allows calculation of the thermal conductivity and shear viscosity, which are compared to the predictions of granular kinetic theory, and which can be used in continuum modeling of granular flows. The shear viscosity is close to the prediction of kinetic theory, while the thermal conductivity can be overestimated by a factor of 2; in each case, transport is lowered with increasing inelasticity.Comment: 14 pages, 17 figures, 39 references, submitted to PRE feb 199

Identification of electron and hole traps in lithium tetraborate (Li2B4O7) crystals: Oxygen vacancies and lithium vacancies

Electron paramagnetic resonance (EPR) and electron-nuclear double resonance (ENDOR) are used to identify and characterize electrons trapped by oxygen vacancies and holes trapped by lithium vacancies in lithium tetraborate (Li2B4O7) crystals. Our study includes a crystal with the natural abundances of 10B and 11B and a crystal highly enriched with 10B. The as-grown crystals contain isolated oxygen vacancies, lithium vacancies, and copper impurities, all in nonparamagnetic charge states. During an irradiation at 77 K with 60 kV x-rays, doubly ionized oxygen vacancies trap electrons while singly ionized lithium vacancies and monovalent copper impurities trap holes. The vacancies return to their preirradiation charge states when the temperature of the sample is increased to approximately 90 K. Hyperfine interactions with 10B and 11B nuclei, observed between 13 and 40 K in the radiation-induced EPR and ENDOR spectra, provide models for the two vacancy-related defects. The electron trapped by an oxygen vacancy is localized primarily on only one of the two neighboring boron ions while the hole stabilized by a lithium vacancy is localized on a neighboring oxygen ion with nearly equal interactions with the two boron ions adjacent to the oxygen ion

Continuum-type stability balloon in oscillated granular layers

The stability of convection rolls in a fluid heated from below is limited by secondary instabilities, including the skew-varicose and crossroll instabilities. We observe a stability boundary defined by the same instabilities in stripe patterns in a vertically oscillated granular layer. Molecular dynamics simulations show that the mechanism of the skew-varicose instability in granular patterns is similar to that in convection. These results suggest that pattern formation in granular media can be described by continuum models analogous to those used in fluid systems.Comment: 4 pages, 6 ps figs, submitted to PR

Four-phase patterns in forced oscillatory systems

We investigate pattern formation in self-oscillating systems forced by an external periodic perturbation. Experimental observations and numerical studies of reaction-diffusion systems and an analysis of an amplitude equation are presented. The oscillations in each of these systems entrain to rational multiples of the perturbation frequency for certain values of the forcing frequency and amplitude. We focus on the subharmonic resonant case where the system locks at one fourth the driving frequency, and four-phase rotating spiral patterns are observed at low forcing amplitudes. The spiral patterns are studied using an amplitude equation for periodically forced oscillating systems. The analysis predicts a bifurcation (with increasing forcing) from rotating four-phase spirals to standing two-phase patterns. This bifurcation is also found in periodically forced reaction-diffusion equations, the FitzHugh-Nagumo and Brusselator models, even far from the onset of oscillations where the amplitude equation analysis is not strictly valid. In a Belousov-Zhabotinsky chemical system periodically forced with light we also observe four-phase rotating spiral wave patterns. However, we have not observed the transition to standing two-phase patterns, possibly because with increasing light intensity the reaction kinetics become excitable rather than oscillatory.Comment: 11 page