We consider a dilute gas of hard spheres in dimension d ≥ 2 that upon collision either annihilate with probability p or undergo an elastic scattering with probability 1 − p. For such a system neither mass, momentum, nor kinetic energy are conserved quantities. We establish the hydrodynamic equations from the Boltzmann equation description. Within the Chapman-Enskog scheme, we determine the transport coefficients up to Navier-Stokes order, and give the closed set of equations for the hydrodynamic fields chosen for the above coarse grained description (density, momentum and kinetic temperature). Linear stability analysis is performed, and the conditions of stability for the local fields are discussed.

Recently van der Meer et al. studied the breakdown of a granular cluster [Phys. Rev. Lett. 88, 174302 (2002)]. We reexamine this problem using an urn model, which takes into account fluctuations and finite-size effects. General arguments are given for the absence of a continuous transition when the number of urns (compartments) is greater than two. Monte Carlo simulations show that the lifetime of a cluster tau diverges at the limits of stability as tau approximately N(1/3), where N is the number of balls. After the breakdown, depending on the dynamical rules of our urn model, either normal or anomalous diffusion of the cluster takes place.

We investigate the problem of ballistically controlled reactions where particles either annihilate upon collision with probability p, or undergo an elastic shock with probability 1-p. Restricting to homogeneous systems, we provide in the scaling regime that emerges in the long time limit, analytical expressions for the exponents describing the time decay of the density and the root-mean-square velocity, as continuous functions of the probability p and of a parameter related to the dissipation of energy. We work at the level of molecular chaos (nonlinear Boltzmann equation), and using a systematic Sonine polynomials expansion of the velocity distribution, we obtain in arbitrary dimension the first non-Gaussian correction and the corresponding expressions for the decay exponents. We implement Monte Carlo simulations in two dimensions, which are in excellent agreement with our analytical predictions. For p<1, numerical simulations lead to the conjecture that unlike for pure annihilation (p=1), the velocity distribution becomes universal, i.e., does not depend on the initial conditions.

We consider the velocity distribution for a granular gas of inelastic hard spheres described by the Boltzmann equation. We investigate both the free of forcing case and a system heated by a stochastic force. We propose a new method to compute the first correction to Gaussian behavior in a Sonine polynomial expansion quantified by the fourth cumulant $a_2$. Our expressions are compared to previous results and to those obtained through the numerical solution of the Boltzmann equation. It is numerically shown that our method yields very accurate results for small velocities of the rescaled distribution. We finally discuss the ambiguities inherent to a linear approximation method in $a_2$.Comment: 9 pages, 8 eps figures include

We study the effects of an external electric field on both the motion of the reaction zone and the spatial distribution of the reaction product, C, in an irreversible A- + B+ -->C reaction-diffusion process. The electrolytes A identical with (A+,A-) and B identical with (B+,B-) are initially separated in space and the ion-dynamics is described by reaction-diffusion equations obeying local electroneutrality. Without an electric field, the reaction zone moves diffusively leaving behind a constant concentration of C's. In the presence of an electric field which drives the reagents towards the reaction zone, we find that the reaction zone still moves diffusively but with a diffusion coefficient which slightly decreases with increasing field. The important electric field effect is that the concentration of C's is no longer constant but increases linearly in the direction of the motion of the front. The case of an electric field of reversed polarity is also discussed and it is found that the motion of the front has a diffusive as well as a drift component. The concentration of C's decreases in the direction of the motion of the front, up to the complete extinction of the reaction. Possible application of the above results to the understanding of the formation of Liesegang patterns in an electric field is briefly outlined.