Microscopic and soliton-like solutions of the Boltzmann Enskog and generalized Enskog equations for elastic and inelastic hard spheres

Transcription

1 arxiv: v2 [math-ph] 4 Mar 204 Microscopic and soliton-like solutions of the Boltzmann Enskog and generalized Enskog equations for elastic and inelastic hard spheres A. S. Trushechkin Steklov Mathematical Institute of the Russian Academy of Sciences Gubkina St. 8, 999 Moscow, Russia National Research Nuclear University MEPhI Kashirskoe Highway 3, 5409 Moscow, Russia trushechkin@mi.ras.ru Abstract N. N. Bogolyubov discovered that the Boltzmann Enskog kinetic equation has microscopic solutions. They have the form of sums of delta-functions and correspond to trajectories of individual hard spheres. But the rigorous sense of the product of the delta-functions in the collision integral was not discussed. Here we give a rigorous sense to these solutions by introduction of a special regularization of the delta-functions. The crucial observation is that the collision integral of the Boltzmann Enskog equation coincides with that of the first equation of the BBGKY hierarchy for hard spheres if the special regularization to the delta-functions is applied. This allows to reduce the nonlinear Boltzmann Enskog equation to the BBGKY hierarchy of linear equations in this particular case. Also we show that similar functions are exact smooth solutions for the recently proposed generalized Enskog equation. They can be referred to as particle-like or soliton-like solutions and are analogues of multisoliton solutions of the Korteweg de Vries equation. Introduction Kinetics of hard sphere gases is a classical field of physical and mathematical research and it is still actively developed [ 4]. One of kinetic equations in the focus of attention is the

2 Boltzmann Enskog kinetic equation. It is a generalization of the Boltzmann equation for the case of moderately dense (rather than dilute) gases. As the Boltzmann equation, it describes the time-irreversible behaviour and entropy production [5, 6]. However, N. N. Bogolyubov [7] (see also [8]) discovered that this equation has timereversible microscopic solutions as well. They have form of sums of delta-functions and correspond to trajectories of individual hard spheres. Another interest feature of these solutions (along with their reversibility) is that the Boltzmann Enskog kinetic equation is believed to give only an approximate aggregate description of the gas dynamics in terms of a single-particle distribution function. The Bogolyubov s result shows that the Boltzmann Enskog equation can describe the gas dynamics in terms of trajectories of individual hard spheres as well. Also the Vlasov kinetic equation is known to have microscopic solutions. They were discovered by A. A. Vlasov himself [9]. However, the justification of these solutions for the Boltzmann Enskog equation is performed by Bogolyubov on the physical level of rigour. In particular, Bogolyubov did not discuss the products of generalized functions in the collision integral. A way to give a rigorous sense to the microscopic solutions is the use of a regularization for delta-functions. A rather general type of regularization and a direct proof of the existence of the limit for the collision integral is proposed in [0]. But the proof is rather involved. Here we present a special kind of regularization. The crucial observation is that, under this regularization, the collision integral of the Boltzmann Enskog equation coincides with that of the first equation of the BBGKY hierarchy for hard spheres. This allows to reduce the nonlinear Boltzmann Enskog equation to the BBGKY hierarchy of linear equations in this particular case. This gives a rigorous sense to the microscopic solutions in a simpler and, in some sense, more natural and elegant way. Also we show that similar construction leads to exact smooth solutions for the recently proposed generalized Enskog equation []. In general, this equation describes the dynamics of infinite number of hard spheres. But our solutions are expressed in terms of the evolution operator of a finite number of hard spheres. These solutions can be referred to as particle-like (regularizations of microscopic solutions, which correspond to trajectories of particles) or soliton-like (evolution of well-localized bumps ) solutions and are analogues of multisoliton solutions of the Korteweg de Vries equation. The rest of the paper is organized as follows. In Section 2, we remind the Boltzmann Enskog equation. The notion of microscopic solution is reminded in Section 3. In Section 4, we construct a regularization for microscopic solutions of the Boltzmann Enskog equation for elastic hard spheres and so, give a rigorous sense to them. In Section 5, we do the same for the case of inelastic hard spheres. Finally, in Section 6 we remind the generalized Enskog equation and, in Section 7, construct explicit soliton-like solutions for it. 2

3 2 Boltzmann Enskog equation for elastic hard spheres The Boltzmann Enskog kinetic equation describes the dynamics of a gas of hard spheres. For elastic hard spheres, it has the form f t = v f + Q(f, f), () r Q(f, f)(r, v, t) = a 2 Ω v (v 2, σ)[f(r, v, t)f(r +aσ, v 2, t) f(r, v, t)f(r aσ, v 2, t)]dσdv 2. Here f = f(r, v, t) 0 is the density function of hard spheres, where r is the position of the center of a sphere, v is its velocity. To avoid the boundary effects, let us consider a gas in a three-dimensional torus T 3 = R 3 /(αz βz γz), where α, β, γ > a. The dynamics of hard spheres in a rectangular box with elastic reflection from the walls can be reduced to the dynamics on a torus. Then r T 3, v R 3, t R, positive constant (the diameter of a sphere). Further, v = v + (v 2, σ)σ, v 2 = v 2 (v 2, σ)σ, r is the gradient in r, a > 0 is a v 2 = v 2 v, σ S 2 (S 2 is a unit sphere in R 3 ), (, ) is a scalar product. v and v 2 have the meaning of the velocities of two hard spheres just before their collision, v and v 2 are the corresponding velocities just after the collision. In such interpretation, σ is a unit vector directed from the center of the second sphere to the center of the first sphere. The expression Q(f, f) is called the collision integral. The integration is performed over the region Ω v = {(σ, v 2 ) : σ S 2, v 2 R 3, (v 2, σ) 0} S 2 R 3. This equation differs from the Boltzmann kinetic equation for hard spheres by the terms ±aσ in the arguments of f in the collision integral. Thus, the Boltzmann equation assumes the size of spheres to be negligibly small (in comparison to the scale of spatial variation of f), while the Boltzmann Enskog equation takes the size of spheres into account. Formally, the Boltzmann equation for hard spheres is obtained from the Boltzmann Enskog equation by the Boltzmann Grad limit: a 0, f is rescaled according to a 2 f. Namely, one can show that a function of the form a 2 f(r, v, t) where f is a solution of the Boltzmann Enskog equation tends to a solution of the Boltzmann equation in the limit a 0 [, 2]. For recent investigations of the Boltzmann Enskog equation see [3]. 3 Microscopic solutions of the Boltzmann-Enskog equation Let (q, 0 w, 0..., qn 0, w0 N ) R6N and qi 0 qj 0 > a for all i j. Also let a system of hard spheres with initial positions q, 0..., qn 0 and velocities w0,..., wn 0 suffer only double collisions, don t suffer so called grazing collisions (with (v 2, σ) = 0 in (2)), and don t suffer infinite number of collisions in a finite time. The measure of the excluded phase points is zero [4]. 3 (2)

4 Let q (t),..., q N (t) and w (t),..., w N (t) be positions and velocities of N hard spheres at time t provided that their initial state is (q 0, w 0,..., q 0 N, w0 N ). Consider the generalized function f(r, v, t) = N δ(r q i (t))δ(v w i (t)). (3) i= N. N. Bogolyubov proved that it satisfies the Boltzmann Enskog equation [7, 8]. He referred to solutions of this form as microscopic solutions, because they correspond to trajectories of individual hard spheres. The Boltzmann Enskog kinetic equation is believed to give only an approximate aggregate description of the gas dynamics in terms of a single-particle distribution function. The Bogolyubov s result shows that the Boltzmann Enskog equation can describe the gas dynamics in terms of trajectories of individual hard spheres as well. Another surprising fact is that the microscopic solutions (for the case of elastic collisions) are time-reversible, i.e., if we replace t t and v v, then a microscopic solution transforms into another microscopic solution. This is a consequence of the reversibility of N hard spheres dynamics. So, the Boltzmann Enskog kinetic equation, which is known to describe an irreversible dynamics and entropy production, has time-reversible solutions as well. However, Bogolyubov performed the proof just on the physical level of rigour. Namely, the formal substitution of function (3) into the Boltzmann Enskog equation yields products of delta-functions in the collision integral, which are not well-defined. The time derivative of this function in the left-hand side of the equation also yields such products because of discontinuity of the velocities w i (t). A way to give a rigorous sense to the notion of microscopic solutions is to introduce a regularization, i.e., to consider a family of smooth delta-like functions of the form f ε (r, v, t) = N δ ε (r q i (t), v w i (t)), i= ε > 0, where δ ε (r, v) δ(r, v) as ε 0, and to prove that the difference between leftand right-hand sides of () vanishes in some space of generalized functions as ε 0 if we substitute f ε instead of f. This is done in [0] under the assumption that the function δ ε is compactly supported and bounded. But the proof is rather involved. Here we consider a special and, in some sense, more natural type of regularization, which leads to a simpler and more elegant proof. 4 Regularization of microscopic solutions Let δ i ε(r, v) : T 3 R 3 R, be a family of functions, ε > 0, i =, 2,..., N. Let the following properties be satisfied for all ε > 0 and i =, 2,..., N: (i) δ i ε(r, v) is continuously differentiable, 4

5 (ii) δε(r, i v)drdv =, T 3 R 3 (iii) δ i ε(r, v) is compactly-supported. Moreover, for every neighbourhood of the origin, there exists ε 0 > 0 such that the supports of δ i ε lies in this neighbourhood for all ε < ε 0 and i =, 2,..., N. It follows from these properties that δ i ε(r, v) δ(r)δ(v) as ε 0, where δ( ) is the (three-dimensional) delta function. Let, further, (q 0, w 0,..., q 0 N, w0 N ) (T3 R 3 ) N and q 0 i q 0 j > a for all i j. As before, let a system of hard spheres with initial positions q 0,..., q 0 N and velocities w0,..., w 0 N suffer only double collisions, don t suffer so called grazing collisions, and don t suffer infinite number of collisions in a finite time. Consider the function f ε (r, v, t), ε > 0, defined as f ε (r, v, 0) f 0 ε (r, v ) = N δε(r i qi 0, v wi 0 ), (4) i= f ε (r, v, t) = f ε (r v (t T k ), v, T k ) + f ε (r v t v N (t t T k )σ, v, T k ) B + r,v,t T k for t (T k, T k+ ] and k = 0,, 2,.... Here f ε (r 2 v 2 t v 2(t t T k ), v 2, T k ) dr 2 dv 2, (5) B + r,v,t = {(r 2, v 2 ) : r 2 r a, t (r, v, r 2, v 2 ) < t}, t = t (r, v, r 2, v 2 ) is the time to collision of two hard spheres in their backward evolution from the phase points (r, v ) and (r 2, v 2 ) (t = if there is no collision of two hard spheres in their backward evolution in the absence of other spheres), v and v 2 are defined by (2) with σ = σ(r, v, r 2, v 2 ) = r 2 r (v 2 v )t (r, v, r 2, v 2 ) S 2 (if t (r 2 r, v 2 v ) < ). (6) Further, 0 = T 0 < T < T 2 <..., T k, are such that every hard sphere suffers no more than one collision on every interval [T k, T k+ ] provided that the initial state of N hard spheres is (q, 0 w, 0..., qn 0, w0 N ). Moreover, let there be no collisions at the instances T, T 2,.... Let us consider the space T = C0(T 3 R 3 [0, )) of functions ϕ(r, v, t) that are continuous, continuously differentiable in r and t and satisfy the property ϕ(r, v, t) = 0 for all (r, v) T 3 R 3 if t > T (T may be different for different functions). This is a space of test functions. The space of generalized functions T is the space of linear continuous functionals over T. 5

6 Theorem. as ε 0 in the space T. f ε t + v f ε + Q(f ε, f ε ) 0 (7) r Proof. Consider the first interval [0, T ] and the function { D ε (r, v,..., r N, v N, t) = N } N! S(N) t fε 0 (r i, v i ). i= Here, S (N) t is the operator of the shift along the phase space trajectory of N hard spheres: S (N) t ϕ(r, v,..., r N, v N ) = ϕ(r (t), V (t),..., R N (t), V N (t)), where ϕ is an arbitrary function and R (t), V (t),..., R N (t), V N (t) are the positions and velocities of N hard spheres at time t provided that their initial positions and velocities are r, v,..., r N, v N. The operator S (N) t is defined almost everywhere in G N, where G k = {(r N k+, v N k+,..., r N, v N ) (T 3 R 3 ) k : r i r j a, i, j =,..., N, i j}. (8) D ε can be regarded as a probability density function for a system of N hard spheres. Consider the single-particle and two-particle distribution functions: F ε (r, v, t) = N D ε (r, v,..., r N, v N, t) dr j dv j, G N j=2 F2 ε (r, v, r 2, v 2, t) = N(N ) D ε (r, v,..., r N, v N, t) dr j dv j. G N 2 They satisfy the first equation of the BBGKY hierarchy for elastic hard spheres (initially derived in [5], rigorously justified in [6 8], a derivation for the most general case is given in [2]): F ε (r, v, t) t F ε (r, v, t) = v + a 2 r In Lemma, we prove that j=3 Ω v (v 2, σ)[f ε 2 (r, v, r + aσ, v 2, t) (9) F ε 2 (r, v, r aσ, v, t)]dσdv 2. (0) F ε 2 (r, v, r aσ, v 2, t) = F ε 2 (r v t, v, r aσ v 2 t, v 2, 0) = f 0 ε (r v t, v )f 0 ε (r aσ v 2 t, v 2 ) () for all r, v, v 2, and σ such that (v 2 v, σ) 0, for all t [0, T ], and for sufficiently small ε. As a corollary, since (v 2 v, σ) = (v 2 v, σ) 0, we have F ε 2 (r, v, r +aσ, v 2, t) = F ε 2 (r v t, v, r +aσ v 2t, v 2, 0) = f 0 ε (r v t, v )f 0 ε (r +aσ v 2t, v 2) 6

7 as well. If we substitute function (5) into the Boltzmann Enskog equation, we will obtain (remind that we consider the time interval [0, T ]) f ε (r, v, t) t f ε (r, v, t) = v + a 2 r Ω v (v 2, σ)[f 0 ε (r v t, v )f 0 ε (r + aσ v 2t, v 2) f 0 ε (r v t, v )f 0 ε (r aσ v 2 t, v 2 )]dσdv 2. By view of (), the collision integral coincides with that of the first equation of the BBGKY hierarchy. This is the crucial observation. We can rewrite the last equation as f ε (r, v, t) t f ε (r, v, t) = v + a 2 r Ω v (v 2, σ)[f ε 2 (r, v, r + aσ, v 2, t) F ε 2 (r, v, r aσ, v 2, t)]dσdv 2. (2) Let us consider the subspace T T T of test functions ϕ(r, v, t) T with additional condition ϕ(r, v, t) = 0 for all r and v if t > T and the corresponding space of generalized functions T T. Let us show that (7) is satisfied in T T. Comparing (0) and (2), we see that, for this aim, it suffices to show that or f ε (r, v, t) t T 3 R 3 [0, ) f ε (r, v, t) + v F ε (r, v, t) F ε (r, v, t) v 0 in T, r t r [ ] [f ε (r, v, t) F ε ϕ(r, v, t) ϕ(r, v, t) (r, v, t)] + v dr dv dt 0 (3) t r as ε 0, for an arbitrary test function ϕ T T. According to Lemma 2, the function F ε (r, v, t) has the form F ε (r, v, t) = N f ε 0 (r v t, v ) fε 0 (r 2 v 2 t, v 2 ) dr 2 dv 2 Br,v,t + fε 0 (r v t v N (t t )σ, v ) fε 0 (r 2 v 2 t v 2(t t ), v 2) dr 2 dv 2, B + r,v,t where Br,v,t = {(r 2, v 2 ) : r 2 r a, t (r, v, r 2, v 2 ) > t}. So, it differs from f ε (r, v, t) only by the factor f N Br ε 0 (r 2 v 2 t, v 2 ) dr 2 dv 2 in the first,v,t term. The region of integration Br,v,t is the set of states (r 2, v 2 ) of the second sphere that do not lead to a collision with the first sphere in time t of the backward evolution given the (4) 7

8 state of the first sphere is (r, v ). Let us express the first term of the right-hand side of (4) as a sum: N f ε 0 (r v t, v ) fε 0 (r 2 v 2 t, v 2 ) dr 2 dv 2 = N B r,v,t N i= δ ε (i) (r v t qi 0, v wi 0 ) and consider the first term in this sum N δ(i) ε (r v t qi 0, v wi 0 ) B r,v,t B r,v,t B r,v,t f 0 ε (r 2 v 2 t, v 2 ) dr 2 dv 2 (5) f 0 ε (r 2 v 2 t, v 2 ) dr 2 dv 2. If sphere suffers no collisions on the interval [0, T ], then fε 0 (r 2 v 2 t, v 2 ) dr 2 dv 2 = fε 0 (r 2 v 2 t, v 2 ) dr 2 dv 2 =, (6) N N (T 3 \B a(r )) R 3 where B a (r ) R 3 is the ball with the center r and the radius a. Now let us assume that sphere collides with another sphere (without loss of generality let it have the number 2) and the moment of the collision is t (0, T ). Then for every t > 0, there exist ε > 0 such that all moments of collisions of spheres with initial phase points (r i, v i ) supp δ ε (i) (r i qi 0, v i wi 0 ), i =, 2, lies in the interval (t t, t + t). This means that N B r,v,t f 0 ε (r 2 v 2 t, v 2 ) dr 2 dv 2 = if t < t t (all phase points (r 2, v 2 ) do not lead to a collision in the backward evolution on time t) and δ ε (i) (r v t qi 0, v wi 0 ) = 0 if t > t + t (all phase points (r 2, v 2 ) lead to a collision, hence, there is no way for the first sphere to be in the phase point (r, v ) at the moment t by free propagation from (r v t, v ) ). So, N δ() ε (r v t q, 0 v w) 0 fε 0 (r 2 v 2 t, v 2 ) dr 2 dv 2 = δ ε () (r v t q, 0 v w) 0 B r,v,t Note that we use the absence of grazing collisions here. Otherwise, there are phase points (r, v ) such that the first sphere can get to them by two ways: by both free propagation and a collision with the second sphere (from different initial states of the two spheres within supp δ ε (i) (r i qi 0, v i wi 0 ), i =, 2). Then δ ε (i) (r v t qi 0, v wi 0) can be non-zero even if all phase points (r 2, v 2 ) lead to a collision. But if the collision of the hard spheres with the initial states (qi 0, w0 i ), i =, 2, is non-grazing, then the two sets {v r : (r, v ) supp δ () (r q 0, v w 0 )} and {v r, r 2, v 2 : (r, v ) supp δ (i) (r q 0 i, v w 0 i ), i =, 2} (v is the velocity of the first sphere after the collision of the spheres with the initial states (r i, v i ), i =, 2) do not intersect for sufficiently small ε, hence, there is at most one way to get to every phase point (r, v ): either by free propagation or by a collision (or neither). 8

9 if t < t t or t > t + t. Other terms in the sum of the right-hand side of (5) are analysed in the same way. Thus, the functions f ε and F ε coincide on [0, T ] apart from infinitesimal neighbourhoods of the moments of collisions. This proves limits (3) and (7) for T T. As ε 0, both f ε and F ε tend to the generalized function δ(r q(t), v w(t)) defined as (δ(r q(t), v w(t)), ϕ(r, v, t)) = 0 ϕ(q(t), w(t), t) dt. Since, by construction, there are no collisions at the moment T, the function f ε (r, v, T ) has the same form as fε 0 (r, v ) (according to Lemma 3) and we can analyse the time interval [T, T 2 ] in the same manner (regard T as the initial time instant). That is, for the interval [T, T 2 ], we set the function D ε to be { N } D ε (r, v,..., r N, v N, t) = N! S(N) (t T ) formulae (9), (0), and (2) go trough, formula () changes to i= f ε (r i, v i, T ) F ε 2 (r, v, r aσ, v 2, t) = F ε 2 (r v (t T ), v, r aσ v 2 (t T ), v 2, T ) = f ε (r v (t T ), v, T )f 0 ε (r aσ v 2 (t T ), v 2, T ). Again, F ε (r, v, t) and fε 0 (r, v, t) coincide on [T, T 2 ] apart from an infinitesimal time interval, which does not matter due to the integration with a test function. Hence, (7) is valid on T T 2. We can continue this procedure up to an arbitrarily large time. For an arbitrary interval [T k, T k+ ] we have { D ε (r, v,..., r N, v N, t) = N } N! S(N) (t T k ) f ε (r i, v i, T k ). (7) The theorem is proved. Lemma. Equality () is satisfied for all r, v, v 2, and σ such that (v 2 v, σ) 0, for all t [0, T ] and for sufficiently small ε. Proof. By construction, each hard sphere suffers no more than one collision on the interval [0, T ] provided that the initial state of N hard spheres is (q, 0 w, 0..., qn 0, w0 N ). Let ε be small such that the same is true for almost every initial state (r, v,..., r N, v N ) supp D ε (, 0) ( almost due to the excluded zero measure initial states, see above). Such ε exists since grazing collisions are prohibited, hence, the required set of initial states is open. The condition (v 2 v, σ) 0 means that two hard spheres are just before their collision. Since every hard sphere suffers no more than one collision, they moved free before the time i=, 9

10 t, so, F2 ε (r, v, r aσ, v 2, t) { = fε 0 (r, v )fε 0 (r aσ, v 2 ) (N 2)! = (N 2)! S (N) t G N 2 } N fε 0 (r i, v i ) dr j dv j i=3 fε 0 (r v t, v )fε 0 (r aσ v 2 t, v 2 )S (N 2) t G N 2 = f 0 ε (r v t, v )f 0 ε (r aσ v 2 t, v 2 ), j=3 { N } N fε 0 (r i, v i ) dr j dv j q.e.d. Here, in the second equality, we have used the fact of free motion of the two hard spheres before time t for sufficiently small ε. i=3 j=3 Lemma 2. Formula (4) is satisfied for t [0, T ]. Proof. Since each hard sphere suffers no more than one collision for almost every initial state (r, v,..., r N, v N ) supp D ε (, 0) given ε is sufficiently small (see the beginning of the proof of the previous lemma), F ε (r, v, t) = N = N (T 3 \B a(r )) R 3 S (2) t F ε 2 (r, v, r 2, v 2, 0) dr 2 dv 2 (T 3 \B a(r )) R 3 S (2) t {f 0 ε (r, v )f 0 ε (r 2, v 2 )} dr 2 dv 2, where B a (r ) R 3 is a ball with center r and radius a. operator S (2) t gives exactly (4). (8) An explicit expression of the Lemma 3. Let t [0, T ] is such that q i (t) q j (t) > a for all i j. Then the function f ε (r, v, t) has the form f ε (r, v, t) = n N i= δ i,t ε (r q i (t), v w i (t)) where the functions δ i,t ε (r, v) obeys properties (i) (iii) of the functions δ i ε(r, v) stated at the beginning of this section. Proof. As we see in the main part of the proof of Theorem and in the proof of Lemma 2, f ε (r, v, t) coincides with the function F (r, v, t) defined by (9) for sufficiently small ε if 0

11 there are no collisions at the moment t. So, we have { N } N f ε (r, v, t) = S (N) t fε 0 (r i, v i ) dr j dv j (N )! G N i= j=2 = fε 0 (R i ( t), V i ( t)) dr j dv j (N )! Here = = = (N )! (N )! (N )! G N i= G N G N i= j=2 {fε 0 (p 0 i, u 0 i )δ(p 0 i R i ( t), u 0 i V i ( t)) dp 0 i du 0 i } dr j dv j J {fε 0 (p 0 i, u 0 i )δ(r i p i (t), v i u i (t)) dp 0 i du 0 i } dr j dv j G N G N i= j=2 {fε 0 (p 0 i, u 0 i )δ(r i p i (t), v i u i (t)) dp 0 i du 0 i } dr j dv j. G N G N i= R i ( t) = R i ( t, r, v,..., r N, v N ), V i ( t) = V i ( t, r, v,..., r N, v N ), i =, 2,..., N, are the positions and velocities of N hard spheres at time t provided that their initial positions and velocities are r, v,..., r N, v N ; p i (t) = p i (t, p 0, u 0,..., p 0 N, u 0 N) u i (t) = u i (t, p 0, u 0,..., p 0 N, u 0 N), i =,..., N, are positions and velocities of N hard spheres at time t provided that their initial positions and velocities are p 0, u 0,..., p 0 N, u0 N ; J = D(R ( t), V ( t),..., R N ( t), V N ( t)) D(r, v,..., r N, v N ) is not defined if r i r j < a for some i j, hence, the derivatives in J are not defined if r i r j = a for some i j. Functions (20) are continuously differentiable in a point (r, v,..., r N, v N ) whenever: is the Jacobian, which emerges due to the definition of a composition of a 6N-dimensional delta function with smooth functions (20) of arguments r,..., v N. According to the Liouville s theorem, J =, but we must prove the continuously differentiability of functions (20) with respect to initial positions and velocities r,..., v N. Their continuous differentiability is not obvious, since the velocities of the spheres are discontinuous at moments of collisions. Also, the operator S (N) t (i) r i r j > a for all i j; j=2 j=2 (9) (20)

12 (ii) there is no collision at the instant t, i.e., for all i j. The reasoning is as follows. R i ( t, r, v,..., r N, v N ) R j ( t, r, v,..., r N, v N ) > a (2) (a) The positions of hard spheres are continuous with respect to the initial state. Hence, for every δ t > 0, there exists δ > 0 such that r i q i (t) < δ t for all i whenever R i ( t) q 0 i < δ and V i ( t) w 0 i < δ for all i. (b) For every δ > 0, there exists ε δ > 0 such that the expression N i= f 0 ε (R i ( t), V i ( t)), ε < ε δ, is non-zero only if R i ( t) q 0 i < δ and V i ( t) w 0 i < δ for all i. (c) By the condition of the lemma, q i (t) q j (t) > a for all i j. (d) From (a), (b), and (c): for sufficiently small ε, we have r i r j > a for all i j, otherwise N i= f 0 ε (R i ( t), V i ( t)) = 0. We have obtained (i). (e) From (b) and from q 0 i q 0 j > a for all i j: R i ( t) R j ( t) > a for all i j for sufficiently small ε, otherwise N i= f 0 ε (R i ( t), V i ( t)) = 0. We have obtained (ii). Thus, for sufficiently small ε, we integrate over a region of initial states r, v,..., r N, v N in which functions (20) are continuously differentiable. Now, f ε (r, v, t) can be expressed as f ε (r, v, t) = = {fε 0 (p 0 i, u 0 i )δ(r i p i (t), v i u i (t)) dp 0 i du 0 i } (N )! G N G N i= δ(r p (t), v u (t)) fε 0 (p 0 i, u 0 i ) dp 0 i du 0 i (N )! G N = n (N )! G N δ(r p (t), v u (t)) i= (j,...,j N ) i= dr j dv j j=2 δ j i ε (p 0 i q 0 j i, u 0 i w 0 j i ) dp 0 i du 0 i, where δ(r, v) δ(r)δ(v), the last sum is performed over all permutations (j,..., j N ) of the indices,..., N. 2

13 For every permutation, let us renumber the pairs (p 0 i, u 0 i ) according to this permutation: f ε (r, v, t) = n δ(r p j (t), v u j (t)) (N )! G N (j,...,j N ) = n δ(r p j (t), v u j (t)) (N )! G N (j,...,j N ) N = n δ(r p i (t), v u i (t)) where G N = n N i= i= δ i,t ε (r q i (t), v w i (t)), j= δε i,t (r q i (t), v w i (t)) = δ(r p i (t), v u i (t)) G N i= δ j i ε (p 0 j i q 0 j i, u 0 j i w 0 j i ) dp 0 i du 0 i δε(p i 0 i qi 0, u 0 i wi 0 ) dp 0 i du 0 i i= δ j ε(p 0 j q 0 j, u 0 j w 0 j ) dp 0 jdu 0 j δε(p j 0 j qj 0, u 0 j wj 0 ) dp 0 jdu 0 j. In the third equality we used the fact the integral depends on j and does not depend on j 2,..., j N. The number of permutations with a fixed j is (N )!. Now the required properties of the functions δε i,t follow from the corresponding properties of the functions δε i and the last expression. Remark. If we define the space of test functions as C = C 0 (T 3 R 3 [0, )) of continuous (not necessarily differentiable in r and t) functions ϕ(r, v, t) such that ϕ(r, v, t) = 0 if t > T (again, T may be different for different functions), then the Boltzmann Enskog equation is asymptotically satisfied in the mild sense: lim ε 0 { f ε (r, v, t) f 0 ε (r v t, v ) t 0 j= } Q(f ε, f ε )(r v (t s), v, s) ds = 0 in C. Note that a global existence theorem for the mild solution of the Boltzmann Enskog equation is proved in [, 2]. 5 The case of inelastic hard spheres Analogous solutions can be constructed for the case of inelastic hard spheres as well. If a collision is inelastic, then formula (2) is modified as v,2 = v ± + µ(g) (v 2 v, σ)σ, 2 where µ(g) (0, ] is the coefficient of restitution (µ for elastic collisions), g = (v 2 v, σ). The Boltzmann Enskog equation for inelastic hard spheres has the form [3] 3

14 Q(f, f)(r, v, t) = na 2 Here f t = v f + Q(f, f), r Ω v (v 2 v, σ)[χf(r, v, t)f(r + aσ, v 2, t) f(r, v, t)f(r aσ, v 2, t)]dσdv 2. v,2 = v ± + µ(g ) 2µ(g ) (v 2 v, σ)σ (22) are the velocities before the inverse collision (i.e., they transforms into v,2 after the collision), χ = µ(g ) D(v, v 2) D(v, v 2 ) (23) (χ for elastic collisions), g = (v 2 v, σ), D( ) denotes the Jacobian. D( ) The microscopic solutions are again defined by (3), but with the new definition of the time evolution q i (t), w i (t), i =,..., N. Let us give a rigorous sense to them by introduction of a similar regularization. Let δε(r, i v), i =,..., N, ε > 0, and fε 0 (r, v ) be the same functions as in the previous section. Let f ε (r, v, t) now be defined by formula f ε (r, v, t) = f ε (r v (t T k ), v, T k ) + χf ε (r v t v N (t t T k )σ, v, T k ) B + r,v,t T k f ε (r 2 v 2 t v 2(t t T k ), v 2, T k ) dr 2 dv 2, (24) for t (T k, T k+ ] and k = 0,, 2,.... Here B + r,v,t, t = t (r, v, r 2, v 2 ), and σ = σ(r, v, r 2, v 2 ) are defined as before. The spaces of test functions T and generalized functions T are also defined as before. The evolution operator S t is now defined as S (N) t ϕ(r, v,..., r N, v N ) = ϕ(r (t), V (t),..., R N (t), V N (t)) D(R (t), V (t),..., R N (t), V N (t)) D(r, v,..., r N, v N ) where, again, R (t), V (t),..., R N (t), V N (t) are the positions and velocities of N inelastic hard spheres at time t provided that their initial positions and velocities are r, v,..., r N, v N. Note that the Jacobian is not defined if r i r j = a for some i j, because R k (t) and V k (t), k =,..., N, are not defined for r i r j < a (i.e., in any neighbourhood of a point where r i r j = a). So, in contrast to the case of elastic spheres, the operator S (N) t is now defined only if r i r j > a for all i j (the inequality is now strict). 4,

15 However, the Jacobian is a constant function of time on time intervals between collisions and has jumps at times of collisions, i.e., if R i (t) R j (t) = a for some i j. So, it has a finite limit as (r,..., v N ) tends to the hypersurface G N = {(r, v,..., r N, v N ) : r i r j a for all i j, r k r l = a for some k l} from G N \ G N (remind that the definition of G N is given by (8)). By this reason, we can define the operator S (N) t on this hypersurface by continuity and consider it to be defined on all G N, as in the case of elastic collisions. Theorem 2. as ε 0 in the space T. f ε t + v f ε + Q(f ε, f ε ) 0 r The proof is completely analogous to that of Theorem (Remark is also valid here). Namely, the functions D ε, F ε, and F2 ε are constructed in the same manner (with the new definition of S (N) t ). Single- and two-particle distribution functions F ε and F2 ε satisfy the first equation of the BBGKY hierarchy for a system of inelastic hard spheres [9, 20]: F ε F ε = v + na 2 (v 2 v, σ)[χf2 ε (r, v t r, r + aσ, v 2, t) Ω v F ε 2 (r, v, r aσ, v 2, t)]dσdv 2. (25) Let us comment this. In [9,20], only a particular case µ = const is considered. However, the analysis can be easily generalized to µ as a function of g = (v 2 v, σ). The main property to be satisfied for the validity of (25) is the following boundary condition. Lemma 4. Let the function D(r, v,..., r N, v N, t) be defined as D(r, v,..., r N, v N, t) = S (N) t D 0 (r, v,..., r N, v N ), where D 0 is a continuous function equal to zero in some neighbourhood of forbidden configurations (i.e., there exists δ > 0 such that D 0 = 0 whenever r i r j < a + δ for some i j). Then the boundary condition D(r, v,..., r i, v i,..., r i aσ, v j..., r N, v N, t) = χd(r, v,..., r i, v i,..., r i aσ, v j..., r N, v N, t) (26) is satisfied for almost all r, v,..., r N, v N, and σ such that (v 2 v, σ) < 0, and for all t. Here v i, v j, and χ are defined by (22) and (23) correspondingly (with and 2 replaced by i and j). Proof. By definition of S (N) t, D(r, v,..., r N, v N, t) = D 0 (R ( t), V ( t),..., R N ( t), V N ( t)) D(R ( t), V ( t),..., R N ( t), V N ( t)) D(r, v,..., r N, v N ). 5

16 Phase trajectories are continuous with respect to time and initial phase point almost everywhere between the collisions. As we said above, the Jacobian is a constant function of time on time intervals between collisions and has jumps at times of collisions. A collision happens if R k ( t) R l ( t) = a for some k l, but the function D 0 is zero in the neighbourhood of these hypersurfaces. So, the function D is continuous almost everywhere in G N R. It is sufficient to consider the case r k r l > a for all pairs k l except the pair (i, j). If we prove the equality for this case, then, by continuity of D, the equality is valid for the general case as well. As we noted above, if r i r j = a, then the operator S (N) t is defined as a limit from G N \ G N. So, the left- and right-hand sides of equality (26) should be understood as limits and Now we use continuity in t: D(r, v,..., r N, v N, t) = lim D(r, v,..., r i, v i,..., r j, v j..., r N, v N, t), r j r i aσ, lim χd(r, v,..., r i, v r j r i aσ, i,..., r j, v j..., r N, v N, t). lim D(r, v,..., r N, v N, t) r j r i aσ, = lim lim D(r, v,..., r N, v N, t + t) t 0 r j r i aσ, = lim lim D(R ( t), V ( t),..., R N ( t), V N ( t), t) t 0 r j r i aσ, D(R ( t), V ( t),..., R N ( t), V N ( t)) D(r, v,..., r N, v N ) = D(r, v,..., r i, v i,..., r j σ, v j..., r N, v N, t) lim t 0 lim r j r i aσ, D(R ( t), V ( t),..., R N ( t), V N ( t)) D(r, v,..., r N, v N ) Obviously, for sufficiently small t and r j r i a, we have R k ( t) = r k v k t, V k ( t) = v k for all k i, j, lim V i( t) = v i, lim V j( t) = v j, r j r i aσ r j r i aσ R i ( t) = r i v i t v i ( t t ), R j ( t) = r j v j t v j ( t t ). Here t = t (r, r 2 ) (0, T ) is the moment of the collision of the spheres i and j in the backward propagation, i.e., R i ( t ) R j ( t ) = a.. 6

17 Now, Further, lim lim t 0 r j r i aσ, = lim t 0 = lim t 0 The parameter lim r j r i aσ, lim r j r i aσ, lim lim t 0 r j r i aσ, D(R ( t), V ( t),..., R N ( t), V N ( t)) D(r, v,..., r N, v N ) D(R i ( t), v i, R j ( t), v D(r i, v i, r j, v j ) D(R i ( t), v i, R j ( t), v j ) D(r i, v i, r j, v j ) j ) D(R i ( t), v i, R j ( t), v j ) D(R i ( t), v i, R j ( t), v j ). D(R i ( t), v i, R j ( t), v j ) D(R i ( t), v i, R j ( t), v j ) = D(v i, v j ) D(v i, v j ). (27) σ = R i( t ) R j ( t ) a depends on both R i,j ( t) and v i,j (because R i,j ( t ) = R i,j ( t)+v i,j( t t )). However, by properties of Jacobians, we can consider R i,j ( t) to be constant in the right-hand side of (27). Also, in the limit, t 0, R i,j ( t) R i,j ( t ), and σ does not depend on v. To finish the proof, we only need to prove that lim r j r i aσ, D(R i ( t), R j ( t)) D(r i, r j ) = µ, where we consider v i,j to be constant in the calculation of the left-hand side. Also remind that µ = µ( (v j v i, σ) ). We have D(R i ( t), R j ( t)) D(r i, r j ) = D(R i ( t) R j ( t), R i ( t) + R j ( t)) D(r i r j, r i + r j ) = DR Dr, where r = r i r j, R = R i( t) R j ( t). (28) a a Here we used that, by the law of momentum conservation, R i ( t) + R j ( t) = r i + r j (v i + v j ) t, D(R i ( t) + R j ( t)) D(r i + r j ) Variables (28) correspond to a reference frame attached to the jth hard sphere. Let us introduce the cylindrical coordinates r = (ρ, ϕ, z) and choose z-axis along v i v j v, i.e., v = (0, 0, v ). Then σ = (ρ, ϕ, ρ 2 ) (see Figure, a = in our rescaled variables), (v v 2, σ) = ρ 2, v i v j v = v + µ µ =. (v, σ)σ = ( + µ µ v ρ ρ 2, ϕ+π, v + µ µ v ( ρ2 )) (v ρ, ϕ+π, v z ), 7

18 R = σ v ( t t ) = (ρ + v ρ( t t ), ϕ, ρ 2 v z ( t t )) (R ρ, ϕ, R z ), µ = µ( ρ 2 ). We have lim lim t 0 r j r i aσ, lim lim t 0 r j r i aσ, t = z ρ 2, v R ρ ρ = v ρ t ρ, R z ρ = ρ + ρ 2 v z (the terms ( v ρ,z/ ρ)( t t ) disappear in the limit). Note also that Rρ (the ratio of ρ Jacobians of transformations from Cartesian to the cylindrical coordinates in R and r) tends to unity in this limit. R ρ and R z depend on z only through t. So, DR lim t 0 Dr = v ρ t v ρ ρ t z ρ ρ + v 2 z t z t = t v z z t ρ v ρ z = ρ 2 µ. lim r j r i aσ, The lemma is proved. v ρ z Again, f ε (r, v, t) and F ε (r, v, t) coincide apart from infinitesimal neighbourhoods of the moments of collisions. The counterparts of Lemmas, 2, and 3 are proved in completely the same way as for the case of elastic hard spheres. In particular, calculations in (9) are performed in such a way (the notations are the same): { N } N f ε (r, v, t) = S (N) t fε 0 (r i, v i ) dr j dv j (N )! G N i= j=2 = = = = (N )! (N )! (N )! (N )! J fε 0 (R i ( t), V i ( t)) dr j dv j G N i= j=2 J {fε 0 (p 0 i, u 0 i )δ(p 0 i R i ( t), u 0 i V i ( t)) dp 0 i du 0 i } dr j dv j G N G N i= j=2 J J {fε 0 (p 0 i, u 0 i )δ(r i p i (t), v i u i (t)) dp 0 i du 0 i } dr j dv j G N G N i= j=2 {fε 0 (p 0 i, u 0 i )δ(r i p i (t), v i u i (t)) dp 0 i du 0 i } dr j dv j. G N G N i= t ρ j=2 8

19 v Ball, t = 0 r z v Ball, t = t Ball 2 σ ρ 2 v ρ Ball, t = t R Figure : To the proof of Lemma 4. Cylindrical coordinates (ρ, ϕ, z), the plane ϕ = const. The frame of reference is attached to the second ball. The first ball is showed in three different time instants. 9

20 6 Generalized Enskog equation Thus, we see that smooth functions of form (5) (or (24) for the case of inelastic hard spheres) satisfy the Boltzmann Enskog equation asymptotically as ε 0. It is interesting that similar construction leads to exact smooth solutions for the recently proposed generalized Enskog equation [], which is more general than the Boltzmann Enskog one. Here we again consider the case of elastic hard spheres, the generalization to the inelastic case is straightforward. The generalized Enskog equation has the form F t = v F r + Q GE (F, F ), (29) Q GE (F, F )(r, v, t) = na 2 Ω v (v 2, σ)[f 2 (r, v, r + aσ, v 2 F (t)) F 2 (r, v, r aσ, v 2 F (t))]dσdv 2. Here F = F (r, v, t) (now we suppose that r R 3 ) has the meaning of the single-particle density function, F 2 = F 2 (r, v, r 2, v 2 F (t)) is the two-particle density function. It depends on time only through F : F 2 (r, v, r 2, v 2 F (t)) = n=0 n+2 n! V(+n) t ({, 2}, 3,..., n + 2) F (r i, v i, t), (30) i= V (+n) t ({Y }, X\Y ) = n! n ( ) k k=0 n m = n m... m k... (n m... m k )! m k = Â(+n m... m k ) t ({Y }, s +, s + 2,..., s + n m... m k ) k j= m j k j 2 =0... k j n m... m j +s k j n m... m j +s =0 s+n m... m j i j = (k j n m... m j +s+ i j k j n m... m j +s+2 i j )! Â(+kj n m... m j +s+ i j k j n m... m j +s+2 i j ) t (i j, s + n m... m j + + k j s+n m... m j +2 i j,..., s + n m... m j + k j s+n m... m j + i j ), where we mean k j = m j, k j n m... m j+s+ = 0. Also, Â (+n) t The operator A (+n) t A (+n) t ({Y }, X\Y ) = s+n ({Y }, X\Y ) = A (+n) t ({Y }, X\Y )I s+n i= is defined as follows: ( ) P ( P )! P : ({Y },X\Y )= i X i 20 X i P A () t (i). S ( θ(x i) ) t (θ(x i )).

21 The notation S (N) t (i, i 2,..., i k ) means that the phase shift operator acts on the variables with the numbers i,..., i k (a function may depend on larger number of pairs of variables (r i, v i )). Now we treat (S (k) t (i, i 2,..., i k )ϕ)(r, v,..., r N, v N ) = 0 whenever r i r j < a for some i j. If k = N, we will still write S (N) t ϕ, because, in this case i,..., i N are defined unambiguously: i =,..., i N = N. Further, {Y } is a set consisting of a single element Y = {, 2,..., s}, i.e., {Y } =, X = {, 2,..., s + n}. The sum is performed over all partitions P of the set {{Y }, X\Y } = {{Y }, s +,..., s + n} into P non-empty mutually disjoint subsets X i {{Y }, X\Y }. The mapping θ is the declasterization mapping defined by the formulae where θ({{,..., s}, i, i 2,..., i k }) = {,..., s, i, i 2,..., i k }, θ({i, i 2,..., i k }) = {i, i 2,..., i k }. The operator I s+n acts on functions f s+n (r, v,..., r s+n, v s+n ) as follows: X k = X k (r,..., r k ) = I s+n f s+n = X s+n f s+n, {, r i r j a for all i j, 0, otherwise. The operator A (n+) t is called (n + )-th order cumulant of operators S, the operator  (n+) t is called (n + )-th order scattering cumulant of operators S. The cumulants of the lowest orders have the form A () t ({, 2}) = S (2) t (, 2), A (2) t ({, 2}, 3) = S (3) t (, 2, 3) S (2) t (, 2)S () t (3). It can be shown [] that the operator V () t ({, 2}) is the identity operator, so, the term n = 0 in (30) corresponds to the Boltzmann Enskog equation. Thus, equation (29) is a generalization of the Boltzmann Enskog equation. The two-particle distribution function F 2 can be expressed in terms of the initial function F 0 (r, v) F (r, v, 0) as follows: { 2 } F 2 (r, v, r 2, v 2 F (t)) = S (2) t + n= i= F 0 (r i, v i ) n+2 A (+n) t ({, 2}, 3,..., n + 2)X n+2 F 0 (r i, v i ) dr 3 dv 3... dr n+2 dv n+2. (3) n! R 6n A solution of (29) is given by F (r, v, t) = F 0 (r v t, v ) n+ + A (+n) t (, 2,..., n + )X n+ F 0 (r i, v i ) dr 2 dv 2... dr n+ dv n+. (32) n! R 6n n= 2 i= i=

22 Note that equation (29) is derived in such a way that it satisfies the dynamics of infinite number of particles with the initial s-particle distribution functions F 0 s (r, v,..., r s, v s ) = s F 0 (r i, v i )X s. i= 7 Explicit soliton-like solutions of the generalized Enskog equation Let δε(r, i v) : R 6 R, i =,..., N, and fε 0 (r, v), ε > 0, be the same functions as in Section 4. Let, additionally, ε is such that δε(r, i v) = 0 if r > a for all i. Now consider the function f ε (r, v, t) defined as { N } N f ε (r, v, t) = S (N) t fε 0 (r i, v i ) dr j dv j (33) (N )! R 6(N ) i= j=2 Theorem 3. The function f ε (r, v, t) is a solution of the generalized Enskog equation (29) for t > 0. Proof. Under the stated conditions, the series in (32) has a maximal term n = N, because there are no N + points from the support of F 0 (r, v) f 0 ε (r, v) with the distance equal or greater than a from each other. So, (32) has the form F ε (r, v, t) = fε 0 (r v t, v ) N n+ + A (n+) t (, 2,..., n + )X n+ fε 0 (r i, v i ) dr 2 dv 2... dr n+ dv n+. n! R 6n n= But this coincides with (33). Indeed, if we substitute here the expression for A n+, then the only surviving term is S (N) t { N i= f ε 0 (r i, v i )}, all other terms cancel each other. Thus, we obtained a family of explicit smooth solutions for the generalized Enskog equation. In general, this equation describes the dynamics of an infinite number of particles. The crucial observation is that special single-particle distribution functions admits only a finite number of particles. In this case, the solution of this equation can be expressed in terms of the evolution operator of a finite number of particles. The constructed solutions can be compared to multisoliton (or particle-like) solutions of the Korteweg de Vries equation: it also describes the evolution of particle-like bumps. The difference is that our solutions maintain their particle-like (localized) form only for a finite time interval. Our solutions can be referred to as particle-like (regularizations of microscopic solutions, which correspond to trajectories of particles) or soliton-like solutions. Let us consider the simplest nontrivial case N = 2. Then equation (29) for f ε is reduced to i= 22

23 f ε t = v f ε + Q (2) GE r (f ε, f ε ), (34) Q (2) GE (f ε, f ε )(r, v, t) = na 2 Ω v (v 2, σ)f ε (r, v, t)f ε (r + aσ, v 2, t)ζ(r, v, t)ζ(r + aσ, v 2, t) where ζ(r, v, t) = f ε (r, v, t)f ε (r aσ, v 2, t)ζ(r, v, t)ζ(r aσ, v 2, t)]dσdv 2, ( B r,v,t f 0 ε (r 2 v 2 t, v 2 ) dr 2 dv 2 ) Indeed, since we deal with the unbounded configuration space, the two particles collide at most once. Hence, an explicit expression of (33) is f ε (r, v, t) = fε 0 (r v t, v ) + B + r,v,t B r,v,t f 0 ε (r 2 v 2 t, v 2 ) dr 2 dv 2 f 0 ε (r v t v (t t )σ, v ) f 0 ε (r 2 v 2 t v 2(t t ), v 2) dr 2 dv 2 (35) (cf. (4), the notations are the same). As wee just saw, in this particular case, formula (32) is reduced to the dynamics of two particles and, hence, (29) is reduced to the first equation of the BBGKY hierarchy. F ε (r, v, t) t F ε (r, v, t) = v + a 2 r Ω v (v 2, σ)[f ε 2 (r, v, r + aσ, v 2, t) F ε 2 (r, v, r aσ, v, t)]dσdv 2. (36) with F ε (r, v, t) = f ε (r, v, t) and F ε 2 (r, v, r 2, v 2, t) = S (2) t { 2 i= f 0 ε (r i, v i )}. Since the particles collide at most once, () is satisfied: F ε 2 (r, v, r aσ, v 2, t) = F ε 2 (r v t, v, r aσ v 2 t, v 2, 0) = f 0 ε (r v t, v )f 0 ε (r aσ v 2 t, v 2 ) for all r, v, v 2, and σ such that (v 2 v, σ) 0, for all t [0, T ]. Since we consider the moment before the collision, the second term in the right-hand side of (35) is zero, because its integral is taken over the states that lead to collision in the backward evolution (see more detailed explanations in the proof of Theorem ). So, F ε 2 (r, v, r aσ, v 2, t) = f 0 ε (r v t, v )f 0 ε (r aσ v 2 t, v 2 ) Substitution of (37) to (36) gives (34). = f ε (r v t, v, t)ζ(r, v, t)f ε (r aσ v 2 t, v 2, t)ζ(r aσ, v 2, t) (37) 23

24 We see, that the corrections to the Boltzmann Enskog equation can be expressed by additional factors ζ (depending on the initial distribution f 0 ). As we see in Section 4, if ε 0, then these factors matter only on infinitesimal time intervals. Thus, for this special form of initial conditions, the generalized Enskog equation and its solutions have especially simple forms (34) and (35) correspondingly. The function t (r, v, r 2, v 2 ) = t (r 2 r, v 2 v ) can be explicitly expressed: it is the minimal solution of the quadratic equation r 2 r (v 2 v )t = a (or if a real solution does not exist or if there are negative solutions), i.e., ( ) (v,r) t a v v 2 r 2 + (v,r)2, r a, a 2 r 2 + (v,r)2 0, (v, r) 0, v (r, v) = 2 v 2, otherwise, where r = r 2 r, v = v 2 v. Finally, let us note that the generalization of the results of this section to the inelastic hard spheres is straightforward: one need only to redefine the operator S (N) t again and to add a factor χ (see Section 5) to the gain term of the collision integrals. 8 Conclusions We have given a rigorous sense to the microscopic solutions discovered by N. N. Bogolyubov. They are time-reversible. So, the reversibility/irreversibility property of the Boltzmann Enskog equation depends on the considered class of solutions. If the considered solutions have the form of sums of delta-functions, then this equation is reversible. If the considered solutions belong to the class of continuously differentiable functions, then the Boltzmann- Enskog equation is irreversible. The qualitative differences in solutions of evolution equations depending on the considered functional space is analysed in [2 23], it is crucial in so called functional mechanics [24 3]. Also we have found exact solutions for recently proposed generalized Enskog equation. The construction of these solutions is similar to the construction of regularizations of microscopic solutions for the Boltzmann Enskog equation. These solutions are expressed in terms of the evolution operator for a finite number of hard spheres. In particular, for the case of two hard spheres, these solutions as well as the generalized Enskog equation itself have especially simple forms. The constructed solutions of the generalized Boltzmann Enskog equation can be referred to as soliton-like solutions and are analogues of multisoliton solutions of the Korteweg de Vries equation. A difference is that our solitons spread out over time. Also it is worthwhile to mention the particle-like (microscopic, singular) wave propagations in solutions of the Boltzmann equation for hard spheres, which dominate the short-time behaviour of solutions (see [32,33], a generalization to hard potential model is given in [34]). On the contrary, fluid-like waves reveal the long-time dissipative behaviour of the solutions. An investigation of relations between the microscopic solutions of the Boltzmann Enskog equation and these two types of waves in solutions of the Boltzmann equation for hard spheres can be of interest. 24

25 Acknowledgements The author is grateful for useful remarks and discussions to H. van Beijeren, T. V. Dudnikova, V. V. Kozlov, S. V. Kozyrev, A. I. Mikhailov, T. Monnai, A. N. Pechen, E. V. Radkevich, S. Sasa, O. G. Smolyanov, H. Spohn, H. Tasaki, D. V. Treshchev, V. V. Vedenyapin, I. V. Volovich, and E. I. Zelenov. This work was partially supported by the Russian Foundation for Basic Research (project mol-a). References [] Gerasimenko V. I., Gapyak I. V. Hard sphere dynamics and the Enskog equation // Kinet. Relat. Models V. 5, N. 3. P ; arxiv: [math-ph]. [2] S. Simonella, Evolution of correlation functions in the hard sphere dynamics, J. Stat. Phys. (204) (to be published); arxiv: [math-ph] [3] N. V. Brilliantov and T. Pöschel, Kinetic theory of granular gases (Oxford University Press, 2004). [4] M. S. Borovchenkova and V. I. Gerasimenko, On the non-markovian Enskog equation for granular gases, J. Phys. A: Math. Theor. 47(3) (204). [5] N. Bellomo and M. Lachowicz, On the asymptotic theory of the Boltzmann and Enskog equations: A rigorous H-theorem for the Enskog equation, Springer Lecture Notes in Mathematics: Mathematical Aspects of Fluid and Plasma Dynamics (99). [6] P. Resibois, H-theorem for the (modified) nonlinear Enskog equation J. Stat. Phys. 9(6) (978). [7] N. N. Bogolyubov, Microscopic solutions of the Boltzmann-Enskog equation in kinetic theory for elastic balls, Theor. Math. Phys. 24(2) (975). [8] N. N. Bogolubov and N. N. Bogolubov, Jnr., Introduction to Quantum Statistical Mechanics (Gordon and Breach, New York, 992). [9] A. A. Vlasov, Many-particle theory and its application to plasma (Gordon and Breach, New York, 96). [0] A. S. Trushechkin, Microscopic solutions of the Boltzmann Enskog equation and the irreversibility problem, Proc. Steklov Inst. Math. (204) (to be published). [] L. Arkeryd and C. Cercignani, On the convergence of solutions of the Enskog equation to solutions of the Boltzmann equation Comm. PDE 4(8 9) (989). [2] L. Arkeryd and C. Cercignani, Global existence in L for the Enskog equation and convergence of the solutions to solutions of the Boltzmann equation, J. Stat. Phys. 59(3/4) (990). 25

26 [3] S.-Y. Ha and S. E. Noh, Global weak solutions and uniform L p -stability of the Boltzmann Enskog equation, J. Differential Equations 25() 25 (20). [4] C. Cercignani, V. I. Gerasimenko, and D. Ya. Petrina, Many-Particle Dynamics and Kinetic Equations (Kluwer Academic Publishing, Dordrecht, 997). [5] C. Cercignani, Theory and Application of the Boltzmann Equation (Scottish Academic Press, Edinburgh, 975; Elsevier, New York, 975); C. Cercignani, On the Boltzmann equation for rigid spheres, Transport Theory and Statistical Physics 2(3) (972). [6] H. Spohn, On the integrated form of the BBGKY hierarchy for hard spheres (985); arxiv:math-ph/ [7] R. Illner and M. Pulvirenti, A derivation of the BBGKY hierarchy for hard sphere particle systems Trans. Theory Stat. Phys. 6(7) (987). [8] D. Ya. Petrina, V. I. Gerasimenko, Mathematical problems of statistical mechanics of a system of elastic balls, Russian Mathematical Surveys 45(3) 53 (990). [9] D. Ya. Petrina and G. L. Garaffini, Analog of the Liouville equation and BBGKY hierarchy for a system of hard spheres with inelastic collisions, Ukrainian Mathematical Journal 57(6) (2005). [20] D. Ya. Petrina and G. L. Garaffini, Solutions of the BBGKY hierarchy for a system of hard spheres with inelastic collisions, Ukrainian Mathematical Journal 58(3) (2006). [2] V. V. Kozlov, Thermal Equilibrium in the sense of Gibbs and Poincaré (Institute of Computer Research, Moscow Izhevsk, 2002) (in Russian); V. V. Kozlov, Thermal Equilibrium in the sense of Gibbs and Poincar e Dokl. Math. 65() (2002). [22] V. V. Kozlov and D. V. Treshchev, Weak convergence of solutions of the Liouville equation for nonlinear Hamiltonian systems, Theor. Math. Phys. 34:3 (2003) [23] V. V. Kozlov and D. V. Treshchev, Evolution of measures in the phase space of nonlinear Hamiltonian systems, Theor. Math. Phys. 36:3 (2003) [24] I. V. Volovich, The irreversibilty problem and functional formulation of classical mechanics, Vestnik Samara State Univ. 8 (2008), 35-55; arxiv: [cond-mat.statmech] [25] I. V. Volovich, Randomness in classical mechanics and quantum mechanics, Found. Phys. 4:3 (20), ; arxiv: [quant-ph] [26] I. V. Volovich, Bogolyubov equations and functional mechanics, Theoret. and Math. Phys. 64:3 (200),