On high energy tails in inelastic gases arxiv:cond-mat/0510108v1 [cond-mat.stat-mech] 5 Oct 2005 R. Lambiotte a,b, L. Brenig a J.M. Salazar c a Physique Statistique, Plasmas et Optique Non-linéaire, Université Libre de Bruxelles, Campus Plaine, Boulevard du Triomphe, Code Postal 231, 1050 Bruxelles, Belgium b SUPRATECS, Université de Liège, B5 Sart-Tilman, B-4000 Liège, Belgium c LRRS UMR-5613 CNRS, Faculté des Sciences Mirande, 9 Av. Alain Savary 21078 Dijon Cedex, France. Abstract We study the formation of high energy tails in a one dimensional kinetic model for granular gases, the so-called Inelastic Maxwell Model. We introduce a timediscretized version of the stochastic process, and show that the original continuous time approach makes the particles have different mean energies. This feature is responsible for the high energy tails in the model. We verify this explanation by computer simulations and analytically, by introducing a linear Lorentz model. Key words: Granular models of complex systems, Random walks and Lévy flights, Kinetic theory PACS: 45.70.Vn, 05.40.Fb, 45.20.Dd 1 Introduction Lorentz gas models have been used historically in order to clarify the features of more complex non-linear kinetic equations. For instance, P. and T. Ehrenfest (1) introduced the so-called wind tree model, in order to discuss the Stosszahlansatz of the Boltzmann equation. In this paper, we will use a similar approach in order to illustrate a generic feature of inelastic gases, namely their anomalous velocity distribution. Inelastic gases have been introduced during the last decade in order to study the low-density limit of granular fluids, i.e. Corresponding author. Email addresses: rlambiot@ulb.ac.be (R. Lambiotte), lbrenig@ulb.ac.be (L. Brenig), jmarcos@u-bourgogne.fr (J.M. Salazar). Preprint submitted to Physica A 20 March 2008
large assemblies of moving macroscopic particles performing inelastic interactions. Due to their inelasticity, inelastic gases have a natural tendency to dissipate their kinetic energy, and to reach total rest state. In other words, their granular temperature, defined kinetically by T <V 2 > where V v u(r, t) and u(r, t) is the local mean velocity, vanishes asymptotically if there are not supplied by an external energy source. The need for a reference state in order to describe these systems has led many authors to focus on their asymptotic velocity distributions. One method consists in looking for self similar solutions, i.e. form preserving solutions whose time dependence occurs through the granular temperature, f(v; t) = 1 f S ( v ). Such a procedure has T(t) T(t) been applied to a large variety of kinetic models for inelastic gases, whose scaling velocity distributions have been shown to generically highlight non- Maxwellian features and overpopulated high energy tail (2; 3). Let us note that the specific shape of the tail may exhibit very different behaviours, depending on the details of the model and on the ways to heat it (4). In this paper, we show on simplified models a mechanism leading to the formation of high energy tails, namely an averaging process between particles which have different mean energies. To do so, we perform Direct Simulation Monte-Carlo (DSMC) simulations of the so-called Inelastic Maxwell Models, and introduce a discrete time version of it. Moreover, in order to perform analytical calculations up to the end, we focus on a linear Lorentz model that confirms the results. 2 Inelastic Maxwell Model In the following, we study scaling solutions in the context of the one-dimensional Inelastic Maxwell Model (IMM). This model (2) derives from a mathematical simplification of the non-linear Boltzmann equation for inelastic hard rods, assuming that the collision rate may be replaced by an average rate, independent of the relative velocity of the colliding particles. The resulting kinetic equation writes: f(v 1 ; t) t + f(v 1 ; t) = dv 2 ( 1 α )f(v 1; t)f(v 2; t) (1) where the primed velocities are the pre-collisional velocities, defined by the collision rule: v 1 = v 1 (1 + α) 2α v (1 + α) 12 v 2 = v 2 + 2α v 12 (2) 2
α is the so-called inelasticity parameter and belongs to the interval ]0:1]. It measures the dissipation of energy at each collision. The elastic limit corresponds to the case α = 1. This stochastic model is the simplest formulation for a non-linear inelastic gas. It is also the only one for which the exact form of the scaling solution has been found (5): f(v; t) = 1 1 v (3) π T(t) [1 + ( T(t) )2 ] 2 together with an exponential law for the granular temperature: T(t) = e (1 α)2 t 2. Recently, we have introduced (6) a time discretized version of (1): f N+1 (v 1 ; t) = dv 2 ( 1 α )f N(v 1 ; t)f N(v 2 ; t) (4) Argument based on the central limit theorem for the summation of random variables show that any initial velocity distribution converges towards a gaussian by the discrete time process (4): f N (v; t) = 1 (2πT N ) 1 2 e v 2 2T N (5) where T N is the temperature after N collision steps. Solution (5) is a scaling solution of the discrete time dynamics. Let us insist on the fact that (5) represents the first occurrence of a Maxwell-Boltzmann distribution in the context of inelastic gases, except in trivial models (7; 8). This remarkable feature has the further advantage to isolate a mechanism responsible for the emergence of high energy tail in these systems. Indeed, the transition from a power law tail (3) towards a Maxwell-Boltzmann distribution (5) is entirely due to the passage from a continuous time towards a discrete time description. In (4), every particle has suffered exactly N collisions at time N. Consequently, every particle has dissipated the same energy on the average at that time. In contrast, in (1), the constituents of the inelastic gas have performed different number of collisions N at a given time t, distributed according to the Poisson law tn e t. N! This implies that, on the average, the different particles are characterized by different quantities of energy. This behaviour suggests that the high energy tail should be due to a mixing between particles with different mean energies. We have performed DSMC simulations (9) of the kinetic equations (1) and (4), which are based on their interpretation as a stochastic process. The case (1) is well-known and consists in picking a given number of random collision pairs at each step. The case (4) is based on the same algorithm, except that at 3
0.045 0.04 (a) 0.035 Velocity distribution 0.03 0.025 0.02 0.015 0.01 0.005 0 0 0.5 1 1.5 2 2.5 3 3.5 4 v 0.1 (b) Velocity distribution 0.01 0.001 0.0001 0 0.5 1 1.5 2 2.5 3 3.5 4 Fig. 1. The system is composed by 150000 particles with α=0.5. The dashed line are a guide for the eye in order to isolate the initial conditions, that are uniform or Maxwell-Boltzmann distribution. In (a), we plot the velocity distributions in arithmetic and in logarithmic scale (inset) for the discrete time model. The solid line corresponds to the theoretical Maxwell-Boltzmann distribution. In (b), we plot the results obtained from the continuous time process, which highlights a convergence towards the solid line, i.e. the theoretical solution (3). In both cases, we plot the asymptotic velocity distribution with circles ( ). v 4
1 m1 m2 m3 m4 Velocity moments 0 10 20 30 40 50 Time steps Fig. 2. We study a system composed by 1000000 particles, with α=0.5. We plot in logarithmic scale the time evolution of the velocity moments m i(t) m i (0), i = 1,2,3,4. The solid lines correspond to the theoretical predictions. The different relaxation rates imply multiscaling, as discussed in the main text. each step, the N particles are distributed randomly into N collision couples. 2 This method ensures that each particle collides one and only one time at each step. At the end of each step, the velocities are rescaled in order to keep the total energy constant (10). The simulation results are in perfect agreement with the above discussion (fig.1). 3 One dimensional Lorentz model Lorentz-like models have already been applied to the study of inelastic gases, in order to highlight the non-equipartition of energy for instance (8). These systems are usually composed by two kinds of particles, the scatterers and the light particles, which are considered in physical limits such as the dynamics of light particles is determined only by their interactions with the scatterers. In this paper, we will consider a model which is very similar to the Ehrenfest wind tree model (1). Its kinetic equation writes: f(v; t) t = v ( 1 α 2f(v ; t) f(v; t)) (6) 5
where the pre-collisional is defined by v = 1 v. This kinetic equation has α already been considered in the literature. It has been shown (12; 13) that it is an alternative formulation of an infinite hierarchy of moments equations obtained from the inelastic Liouville equation (14). This equation is also related to population dynamics (12). In the following, we use an IMM-like approximation in order to treat this kinetic equation, i.e. we assume that the collision rate v may be replaced by a mean value T(t) and we rescale the time scale. We call the resulting model the Inelastic Lorentz Walk (ILW): f(v; t) t = 1 α f(v φ(k; t) ; t) f(v; t) α t = φ(αk; t) φ(k; t) (7) where the right side equation is the kinetic equation for the characteristic function of the velocity distribution φ(k; t) = 1 2π dvf(v, t)e ikv. Since the use of characteristic functions is particularly suited for the study of the tail properties of velocity distributions, we focus on the expressions in Fourier space since now on. The discrete time version of the process reads: φ N+1 (k; t) = φ N (αk; t) φ N (k; t) = φ 0 (α N k; t) (8) where we have solved the process recurrently. One should note that this stochastic process leaves the solution scale invariant, i.e. (8) is a self-similar solution whose form is that of the initial velocity distribution and whose temperature vanishes asymptotically. Obviously, the discrete process does not lead to the formation of high energy tails. This feature ceases to be true if only a fraction p of the discs performs a collision at each step. In this case, the state of the system after one step writes: φ 1 (k; t) = (1 p) φ 0 (k; t) + p φ 0 (αk; t) (9) Simple analytical calculations show that the kurtosis κ <v4 > <v 2 > 2 of the velocity distribution increases through (9) whatever the initial velocity distribution and proportion p 1. This is due to the fact that the superposition of identical distributions with different temperatures T 1, T 2 has a higher energy tail than the same distribution with the mean temperature T (1-p) T 1 + p T 2. This property applied to the Maxwell-Boltzmann distribution is well known and has been used by Beck (16) in order to construct the so-called superstatistics. Contrary to the non-linear IMM, the link between the discrete and the continuous time dynamics is straightforward. Methods from random walks theory (15) lead to: φ(k; t) = N=0 t N e t φ 0 (α N k; t) (10) N! 6
The system at time t is thus composed by different classes of particles which are characterized by their mean energy α N, or equivalently by the number N of interactions they have performed. The proportion of particles in these classes is the probability for a particle of having performed N collisions before time t tn e t N!. On should note that (10) consists in the superposition of the initial distribution at different energy, as in (9) and leads to the formation of high energy tails by the same mechanism. Assuming that the small k development of the initial characteristic function is φ 0 (k; t) = i=0 c 0i k i i!, the formal solution (10) reads: c 0i k i φ(k; t) = i! i=0 e (1 αi )t (11) Detailed analsis of this solution shows that the inelastic Lorentz walk exhibits multiscaling properties, i.e. the dimensionalized velocity moments m i (< (v 2 ) i >) 1 2i decay with different cooling rates: m i e λit. The cooling rates write λ i = 1 αi and are decreasing function of i: λ i i < λ j, i > j, i.e. the system is characterized by an infinite number of independent cooling rates. Consequently, the velocity moments m i, i 2 grow towards infinity, as compared to m 1. This behaviour is verified by computer simulations of the linear stochastic process (fig.2). Finally, let us stress that our computer simulations show that the ILW does not converge toward a scaling distribution, at long time, contrary to non-linear kinetic models. 4 Conclusions In this paper, we consider the formation of high energy tails of the scaling velocity distributions observed in inelastic gases. To do so, we focus on the one-dimensional Inelastic Maxwell Model, which consists in a mean field approximation of the Boltzmann equation for inelastic hard rods. By comparing a discrete and a continuous time version of the process, we discuss a mechanism responsible for the emergence of high energy tails. In a first approach, we perform DSMC simulations of the gas, thereby confirming theoretical predictions. Then, we introduce a simplified linear model in order to avoid technical difficulties, i.e. the Inelastic Lorentz Walk. Analytical arguments show that this effect is responsible for the broadening of the initial velocity distribution and for the formation of high energy tails in the Lorentz model. Moreover, we show and verify numerically that the system exhibits multiscaling of the velocity moments. Let us stress that this mechanism applies to inelastic gases, where energy is dissipated at each collision, and has no counterpart in elastic gases. It should also lead to strong consequences in inhomogeneous systems, where the particles average energies depend on space. 7
Acknowledgments R.L. would like to thank warmly M. Mareschal and M. Ausloos. This work has been done thanks to the financial support of FRIA, and ARC. References [1] P. and T. Ehrenfest, The conceptual foundations of the statistical approach in mexhanics, Cornell Univ. Press (1959) [2] E. Ben-Naim and P.L. Krapivsky Phys. Rev. E, 61(2000) R5 [3] M.H. Ernst and R.Brito, Europhys. Lett., 58 (2002) 182 [4] T. Biben, P.A. Martin, and J. Piasecki, Physica A, 310 (2002) 308 [5] A. Baldassarri, U.M.B. Marconi and A. Puglisi, Europhys. Lett., 58 (2002) 14 [6] R. Lambiotte and L. Brenig, Phys. Lett. A, 345 (2005) 309 [7] J. Brey et al., J. Stat. Phys., 97 (1999) 281 [8] P.A. Martin and J. Piasecki, Europhys. Lett., 46 (1999) 613 [9] G.A. Bird, Molecular Gas Dynamics, Clarendon Press (1976) [10] R. Soto and M. Mareschal, Phys. Rev. E, 63 (2001) 041303 [11] E. Wild Proc. Cambridge Philos. Soc., 47 (1951), 602 [12] Y. Elskens, J. Stat. Phys., 101 (2000) 397 [13] R. Lambiotte, phd thesis [14] J.M. Salazar and L. Brenig, Phys. Rev. E, 46 (1999) 2093 [15] R. Balescu, Statistical dynamics. Matter out of equilibrium Imperial college press (1997) [16] C. Beck, Phys. Rev. Letters, 87 (2001) 180601 8