Fluctuating hydrodynamics for dilute granular gases

Lattice Model to Derive the Fluctuating Hydrodynamics of Active Particles with Inertia

We derive the hydrodynamic equations with fluctuating currents for the density, momentum, and energy fields for an active system in the dilute limit. In our model, nonoverdamped self-propelled particles (such as grains or birds) move on a lattice, interacting by means of aligning dissipative forces and excluded volume repulsion. Our macroscopic equations, in a specific case, reproduce a transition line from a disordered phase to a swarming phase and a linear dispersion law accounting for underdamped wave propagation. Numerical simulations up to a packing fraction ∼10% are in fair agreement with the theory, including the macroscopic noise amplitudes. At a higher packing fraction, a dense-diluted coexistence emerges. We underline the analogies with the granular kinetic theories, elucidating the relation between the active swarming phase and granular shear instability.

Power-law decay of the velocity autocorrelation function of a granular fluid in the homogeneous cooling state

The hydrodynamic part of the velocity autocorrelation function of a granular fluid in the homogeneous cooling state has been calculated by using mode-coupling theory for a finite system with periodic boundary conditions. The existence of the shearing instability, leading to a divergent behavior of the velocity flow fluctuations, is taken into account. A time region in which the velocity autocorrelation function exhibits a power-law decay, when time is measured by the number of collisions per particle, has been been identified. Also the explicit form of the exponential asymptotic long time decay has been obtained. The theoretical prediction for the power-law decay is compared with molecular dynamics simulation results, and a good agreement is found, after taking into account finite size corrections. The effects of approaching the shearing instability are also explored.

Nonlinear driven diffusive systems with dissipation: fluctuating hydrodynamics

We consider a general class of nonlinear diffusive models with bulk dissipation and boundary driving and derive its hydrodynamic description in the large size limit. Both the average macroscopic behavior and the fluctuating properties of the hydrodynamic fields are obtained from the microscopic dynamics. This analysis yields a fluctuating balance equation for the local energy density at the mesoscopic level, characterized by two terms: (i) a diffusive term, with a current that fluctuates around its average behavior given by nonlinear Fourier's law, and (ii) a dissipation term which is a general function of the local energy density. The quasielasticity of microscopic dynamics, required in order to have a nontrivial competition between diffusion and dissipation in the macroscopic limit, implies a noiseless dissipation term in the balance equation, so dissipation fluctuations are enslaved to those of the density field. The microscopic complexity is thus condensed in just three transport coefficients-the diffusivity, the mobility, and a new dissipation coefficient-which are explicitly calculated within a local equilibrium approximation. Interestingly, the diffusivity and mobility coefficients obey an Einstein relation despite the fully nonequilibrium character of the problem. The general theory here presented is applied to a particular albeit broad family of systems, the simplest nonlinear dissipative variant of the so-called KMP model for heat transport. The theoretical predictions are compared to extensive numerical simulations, and an excellent agreement is found.

Internal energy fluctuations of a granular gas under steady uniform shear flow

The stochastic properties of the total internal energy of a dilute granular gas in the steady uniform shear flow state are investigated. A recent theory formulated for fluctuations about the homogeneous cooling state is extended by analogy with molecular systems. The theoretical predictions are compared with molecular dynamics simulation results. Good agreement is found in the limit of weak inelasticity, while systematic and relevant discrepancies are observed when the inelasticity increases. The origin of this behavior is discussed.

Structure factors in granular experiments with homogeneous fluidization

Velocity and density structure factors are measured over a hydrodynamic range of scales in a horizontal quasi-2D fluidized granular experiment, with packing fractions φ ∈ [10%, 40%]. The fluidization is realized by vertically vibrating a rough plate, on top of which particles perform a Brownian-like horizontal motion in addition to inelastic collisions. On one hand, the density structure factor is equal to that of elastic hard spheres, except in the limit of large length-scales, as it occurs in the presence of an effective interaction. On the other hand, the velocity field shows a more complex structure which is a genuine expression of a non-equilibrium steady state and which can be compared to a recent fluctuating hydrodynamic theory with non-equilibrium noise. The temporal decay of velocity modes autocorrelations is compatible with linear hydrodynamic equations with rates dictated by viscous momentum diffusion, corrected by a typical interaction time with the thermostat. Equal-time velocity structure factors display a peculiar shape with a plateau at large length-scales and another one at small scales, marking two different temperatures: the "bath" temperature T(b), depending on shaking parameters, and the "granular" temperature T(g) < T(b), which is affected by collisions. The two ranges of scales are separated by a correlation length which grows with φ, after proper rescaling with the mean free path.

Fluctuating hydrodynamics in a vertically vibrated granular fluid with gravity

We investigate hydrodynamic fluctuations in a two-dimensional granular fluid excited by a vibrating base and in the presence of gravity, focusing on the transverse velocity modes. Since the system is inhomogeneous, we measure fluctuations in horizontal layers whose width is smaller than the characteristic hydrodynamic lengths: they can be considered as almost-homogeneous subsystems. The large time decay of the autocorrelations of modes is exponential and compatible with vorticity diffusion due to shear viscosity, as in equilibrium fluids. The velocity structure factor, which strongly deviates from the equilibrium constant behavior, is well reproduced by an effective fluctuating hydrodynamics described by two noise terms: the first associated with vorticity diffusion and the second with the local energy exchange, which have internal and external character, respectively.

Fluctuating Navier-Stokes equations for inelastic hard spheres or disks

Starting from the fluctuating Boltzmann equation for smooth inelastic hard spheres or disks, closed equations for the fluctuating hydrodynamic fields to Navier-Stokes order are derived. This requires deriving constitutive relations for both the fluctuating fluxes and the correlations of the random forces. The former are identified as having the same form as the macroscopic average fluxes and involving the same transport coefficients. On the other hand, the random force terms exhibit two peculiarities as compared with their elastic limit for molecular systems. First, they are not white but have some finite relaxation time. Second, their amplitude is not determined by the macroscopic transport coefficients but involves new coefficients.

Continuum simulations of shocks and patterns in vertically oscillated granular layers

We study interactions between shocks and standing-wave patterns in vertically oscillated layers of granular media using three-dimensional, time-dependent numerical solutions of continuum equations to Navier-Stokes order. We simulate a layer of grains atop a plate that oscillates sinusoidally in the direction of gravity. Standing waves form stripe patterns when the accelerational amplitude of the plate's oscillation exceeds a critical value. Shocks also form with each collision between the layer and the plate; we show that pressure gradients formed by these shocks cause the flow to reverse direction within the layer. This reversal leads to an oscillatory state of the pattern that is subharmonic with respect to the plate's oscillation. Finally, we study the relationship between shocks and patterns in layers oscillated at various frequencies and show that the pattern wavelength increases monotonically as the shock strength increases.

Fluctuating hydrodynamics for dilute granular gases: a Monte Carlo study

We investigate hydrodynamic noise in a dilute granular gas during the homogeneous cooling state, by means of a proper application of the direct simulation Monte Carlo (DSMC) algorithm. The DSMC includes a source of randomization which is not present in molecular dynamics (MD) for inelastic hard disks. Notwithstanding this difference, a fair quantitative agreement is found, including a violation of the fluctuation-dissipation relation for the noise amplitude of the same order observed in MD. This study suggests that deterministic collision dynamics is not an essential ingredient to reproduce, up to a good degree of approximation, hydrodynamic fluctuations in dilute granular gases.

Effects of thermal noise on pattern onset in continuum simulations of shaken granular layers

The author investigates the onset of patterns in vertically oscillated layers of dissipative particles using numerical solutions of continuum equations to Navier-Stokes order. Above a critical accelerational amplitude of the cell, standing waves form stripe patterns which oscillate subharmonically with respect to the cell. Continuum simulations neglecting interparticle friction yield pattern wavelengths consistent with experiments using frictional particles. However, the critical acceleration for standing-wave formation is approximately 10% lower in continuum simulations without added noise than in molecular-dynamics simulations. This Brief Report incorporates fluctuating hydrodynamics theory into continuum simulations by adding noise terms with no fit parameters; this modification yields a critical acceleration in agreement with molecular-dynamics simulations.