Ricardo J Alonso

AND Bertrand Lods
email: bertrand.lods@unito.it

Isabelle Tristani
email: isabelle.tristani@ens.fr

FROM BOLTZMANN EQUATION FOR GRANULAR GASES TO A MODIFIED NAVIER-STOKES-FOURIER SYSTEM

Keywords: Inelastic Boltzmann equation, Granular flows, Nearly elastic regime, Long-time asymptotic, Incompressible Navier-Stokes hydrodynamical limit, Knudsen number

 which provides the first rigorous derivation of hydrodynamic equations from the Boltzmann equation for inelastic hard spheres in 3D. In particular, we obtain a new system of hydrodynamic equations describing granular flows and prove existence of classical solutions to the aforementioned system. One of the main issue is to identify the correct relation between the restitution coefficient (which quantifies the rate of energy loss at the microscopic level) and the Knudsen number which allows us to obtain non trivial hydrodynamic behavior. In such a regime, we construct strong solutions to the inelastic Boltzmann equation, near thermal equilibrium whose role is played by the so-called homogeneous cooling state. We prove then the uniform exponential stability with respect to the Knudsen number of such solutions, using a spectral analysis of the linearized problem combined with technical a priori nonlinear estimates. Finally, we prove that such solutions converge, in a specific weak sense, towards some hydrodynamic limit that depends on time and space variables only through macroscopic quantities that satisfy a suitable modification of the incompressible Navier-Stokes-Fourier system.

INTRODUCTION

In this paper, we report on some recent results obtained in [START_REF] Alonso | Fluid dynamic limit of Boltzmann equation for granular hard-spheres in a nearly elastic regime[END_REF] about the problem of deriving rigorously some hydrodynamic limit from the Boltzmann equation for inelastic hard spheres with small inelasticity. Our aim here is to give an account of the main aspects of our work [START_REF] Alonso | Fluid dynamic limit of Boltzmann equation for granular hard-spheres in a nearly elastic regime[END_REF] in a shorterreader-friendly -version that includes the main results as well as the main ideas and arguments. We shall only sketch the proofs of our results, referring the reader to [START_REF] Alonso | Fluid dynamic limit of Boltzmann equation for granular hard-spheres in a nearly elastic regime[END_REF] for complete versions and details.

The problem.

The kinetic model. We consider here the (freely cooling) Boltzmann equation which provides a statistical description of identical smooth hard spheres suffering binary and inelastic collisions:

∂ t F + v • ∇ x F = Q α (F, F ) (1.1)
supplemented with initial condition F (0, x, v) = F in (x, v), where F = F (t, x, v) is the density of granular gases having position x ∈ T d ℓ and velocity v ∈ R d at time t 0. We consider here for simplicity the case of flat torus

T d ℓ = R d /(2π ℓ Z) d (1.2)
Acknowledgements. RA gratefully acknowledges the support from O Conselho Nacional de Desenvolvimento Científico e Tecnológico, Bolsa de Produtividade em Pesquisa -CNPq (303325/2019-4). BL gratefully acknowledges the financial support from the Italian Ministry of Education, University and Research (MIUR), "Dipartimenti di Eccellenza" grant 2018-2022. IT thanks the ANR EFI: ANR-17-CE40-0030 and the ANR SALVE: ANR-19-CE40-0004 for their support. Data sharing not applicable to this article as no datasets were generated or analysed during the current study.

Mathematics Subject Classification (2010): 76P05 Rarefied gas flows, Boltzmann equation [See also 82B40, 82C40, 82D05]; 76T25 Granular flows [See also 74C99, 74E20]; 47H20 Semigroups of nonlinear operators [See also 37L05, 47J35, 54H15, 58D07], 35Q35 PDEs in connection with fluid mechanics; 35Q30 Navier-Stokes equations [See also 76D05, 76D07, 76N10].

for some typical length-scale ℓ > 0. The so-called restitution coefficient α belongs to (0, 1] and the collision operator Q α is defined in weak form as

ˆRd Q α (g, f )(v) ψ(v) dv = 1 2 ˆR2d f (v) g(v * ) |v -v * |A α [ψ](v, v * ) dv * dv, (1.3)
where

A α [ψ](v, v * ) := ˆSd-1 (ψ(v ′ ) + ψ(v ′ * ) -ψ(v) -ψ(v * )) b(σ • q) dσ, (1.4)
and the post-collisional velocities (v ′ , v ′ * ) are given by

v ′ = v + 1 + α 4 (|q|σ -q), v ′ * = v * - 1 + α 4 (|q|σ -q),
where q = v -v * , q = q/|q|.

(1.5)

Here, dσ denotes the Lebesgue measure on S d-1 and the angular part b = b(σ • q) of the collision kernel appearing in (1.4) is a non-measurable mapping integrable over S d-1 . There is no loss of generality assuming

ˆSd-1 b(σ • q) dσ = 1, ∀ q ∈ S d-1 .
Notice that one can also give a strong formulation of the collision operator Q α (see [START_REF] Alonso | Fluid dynamic limit of Boltzmann equation for granular hard-spheres in a nearly elastic regime[END_REF]Appendix A]). This strong formulation is simpler in the elastic case (α = 1), we here give it for later use:

Q 1 (g, f )(v) = ˆRd ×S d-1 (g(v ′ * )f (v ′ ) -g(v * )f (v)) |v -v * | b(σ • q) dσ dv * . (1.6)
The true definition actually involves pre-collisional velocities and not post-collisional velocities v ′ and v ′ * but they match in the elastic case, which explains the formula (1.6).

The fundamental distinction between the classical elastic Boltzmann equation and the associated to granular gases lies in the role of the parameter α ∈ (0, 1), the coefficient of restitution that we suppose constant. This coefficient is given by the ratio between the magnitude of the normal component (along the line of separation between the centers of the two spheres at contact) of the relative velocity after and before the collision. The case α = 1 corresponds to perfectly elastic collisions where kinetic energy is conserved. However, when α < 1, part of the kinetic energy of the relative motion is lost since

|v ′ | 2 + |v ′ * | 2 -|v| 2 -|v * | 2 = - 1 -α 2 4 |q| 2 (1 -σ • q) 0.
Notice that the microscopic description (1.5) preserves the momentum Consequently, there is no loss of generality in assuming that

v ′ + v ′ * = v + v *
R(t) = R(0) = 1, U (t) = U (0) = 0, ∀ t 0.
The main contrast between elastic and inelastic gases is that in the latter the granular temperature,

T (t) := 1 |T d ℓ | ˆRd ×T d ℓ |v| 2 F (t, x, v) dv dx is constantly decreasing d dt T (t) = -(1 -α 2 )D α (F (t), F (t)) 0,
where D α ( • , • ) denotes the normalised energy dissipation associated to Q α , see [START_REF] Mischler | Cooling process for inelastic Boltzmann equations for hard-spheres, Part II: Self-similar solution and tail behavior[END_REF], given by

D α (g, g) := γ b 4 ˆTd ℓ dx |T d ℓ | ˆRd ×R d g(x, v)g(x, v * )|v -v * | 3 dv dv * , (1.7)
where γ b is a positive constant depending only on the angular kernel b.

The problem of hydrodynamic limits. To capture some hydrodynamic behaviour of the gas, we need to write the above equation in nondimensional form introducing the dimensionless Knudsen number which is proportional to the mean free path between collisions. We then introduce the classical Navier-Stokes rescaling of time and space (see [START_REF] Bardos | Fluid dynamic limits of kinetic equations. I. Formal derivations[END_REF]) to capture the hydrodynamic limit and introduce the particle density

F ε (t, x, v) := F t ε 2 , x ε , v , t 0. (1.8)
In this case, we choose for simplicity ℓ = ε in (1.2) which ensures now that F ε is defined on

R + × T d × R d with T d := T d 1 .
Under such a scaling, F ε satisfies the rescaled Boltzmann equation

ε 2 ∂ t F ε + ε v • ∇ x F ε = Q α (F ε , F ε ) on T d × R d , (1.9a)
supplemented with the initial condition

F ε (0, x, v) = F in ε (x, v) := F in ( x ε , v). (1.9b)
Conservation of mass and density is preserved under this scaling, if

F ε solves (1.9a), then d dt R ε (t) = d dt U ε (t) = 0
where R ε (t) := ´Td ×R d F ε (t, x, v) dv dx and U ε (t) := ´Td ×R d F ε (t, x, v)v dv dx, whereas the cooling of the granular gas is given by the equation

d dt T ε (t) = - 1 -α 2 ε 2 D α (F ε (t), F ε (t)), (1.10)
where

T ε (t) := 1 |T d | ´Td ×R d |v| 2 F ε (t, x, v
) dv dx and we recall that D α is defined in (1.7). The conservation properties of the equation imply that there is no loss of generality assuming that

R ε (t) = 1, U ε (t) = 0, ∀ ε > 0, t 0.
In order to understand the free-cooling inelastic Boltzmann equation (1.9a)-(1.9b), we perform a selfsimilar change of variables, which allows us to introduce an intermediate asymptotic and ensures that our equation has a non trivial steady state (see [START_REF] Mischler | Cooling process for inelastic Boltzmann equations for hard-spheres, Part I: The Cauchy Theory[END_REF][START_REF] Mischler | Cooling process for inelastic Boltzmann equations for hard-spheres, Part II: Self-similar solution and tail behavior[END_REF][START_REF] Mouhot | Stability, convergence to self-similarity and elastic limit for the Boltzmann equation for inelastic hardspheres[END_REF] for more details). After this change of variables, we are led to study the equation

ε 2 ∂ t f ε + εv • ∇ x f ε + (1 -α) ∇ v • (vf ε ) = Q α (f ε , f ε ), (1.11)
with initial condition

f ε (0, x, v) = F in ε (x, v
). Note that the drift term acts as an energy supply which prevents the total cooling down of the gas. It has been shown that there exists a spatially homogeneous steady state G α to (1.11). More specifically, there exists α 0 ∈ (0, 1) (where α 0 is an explicit threshold value) such that for α ∈ (α 0 , 1), there exists a unique distribution

G α = G α (v) satisfying (1 -α)∇ v • (v G α ) = Q α (G α , G α ) with ˆRd G α (v) 1 v dv = 1 0 . (1.12)
Moreover, there exists some constant C > 0 independent of α such that

G α -M L 1 v ( v 2 ) C(1 -α) (1.13)
where M is the Maxwellian distribution

M(v) := (2πϑ 1 ) -d/2 exp - |v| 2 2ϑ 1 , v ∈ R d , (1.14)
for some explicit temperature ϑ 1 > 0. The Maxwellian distribution M is a steady solution for α = 1 and its prescribed temperature ϑ 1 (which ensures (1.13) to hold) will play a role in the rest of the analysis.

It is important to emphasize that, in all the sequel, all the threshold values on ε and the various constants involved are actually depending only on this initial choice.

In order to reach some incompressible Navier-Stokes type equation in the limit ε → 0, we introduce the following fluctuation h ε around the equilibrium G α :

f ε (t, x, v) = G α (v) + ε h ε (t, x, v).

Our problem boils down to look at the following equation on h

ε :      ∂ t h ε + 1 ε v • ∇ x h ε = 1 ε 2 L α h ε + 1 ε Q α (h ε , h ε ) h ε (t = 0) = h in ε := 1 ε (F in ε -G α ), (1.15)
where L α is the linearized collision operator (local in the x-variable) defined as

L α h := Q α (G α , h) + Q α (h, G α ) -(1 -α)∇ v • (vh). (1.16)
We also denote by L 1 the linearized operator around G 1 = M, that is,

L 1 h := Q 1 (M, h) + Q 1 (h, M). (1.17)
From now on, we will always assume that

ˆTd ×R d F in ε (x, v)   1 v |v| 2   dv dx =   1 0 E ε   with E ε > 0 and E ε -dϑ 1 ε ---→ ε→0 0. (1.18)
The choice of prescribing as initial energy some constant E ε > 0 satisfying ε -1 (E ε -dϑ 1 ) → 0 as ε → 0 for our problem is natural because dϑ 1 is the energy of the Maxwellian M introduced in (1.14) and as we shall see later on, the restitution coefficient α is intended to tend to 1 as ε goes to 0 in our analysis (see (1.22)).

It is also worth noticing that assumption (1.18) and (1.12) result in

ˆTd ×R d h in ε (x, v) 1 v dv dx = 0 0 .
Moreover, equation (1.15) preserves mass and vanishing momentum since, if h ε solves (1.15), then one formally has

d dt ˆTd ×R d h ε (t, x, v)v dv dx = ˆTd ×R d ∇ v • (vh ε (t, x, v))v dv dx = - ˆTd ×R d h ε (t, x, v)v dv dx. (1.19)
Consequently, there is no loss of generality assuming that

ˆTd ×R d h ε (t, x, v) 1 v dv dx = 0 0 , ∀ t 0. (1.20)

Relation between the restitution coefficient and the

(ε) = 1 -ε 2 (λ 0 + η(ε))
(1.21) with λ 0 0 and some function η( • ) that tends to 0 as ε goes to 0. If λ 0 = 0, we assume furthermore that there exists ε ⋆ > 0 such that η( • ) is positive on (0, ε ⋆ ).

Notice that under this assumption, the hypothesis made on the energy of the initial data in (1.18)

implies that ˆTd ×R d h in ε (x, v) |v| 2 dv dx ---→ ε→0 0. (1.22)
Indeed, using (1.18) and Assumption 1.1 combined with (1.13), we obtain

ˆTd ×R d h in ε (x, v)|v| 2 dv dx = 1 ε ˆTd ×R d F in ε (x, v) -G α(ε) (v) |v| 2 dv dx = E ε -dϑ 1 ε + 1 ε ˆTd ×R d M(v) -G α(ε) (v) |v| 2 dv dx ---→ ε→0 0.
Still under Assumption 1.1, we formally obtain that if ε → 0 in (1.15), then h ε → h with h ∈ Ker L 1 where L 1 is defined in (1.17). We recall that when seeing L 1 as an operator acting only on velocity on the space

L 2 v (M -1/2 ), then Ker L 1 = Span{M, v 1 M, • • • v d M, |v| 2 M}
and the projection π 0 onto Ker L 1 is given by

π 0 (g) := d+2 i=1 ˆRd g Ψ i dv Ψ i M, (1.23)
where

Ψ 1 (v) := 1, Ψ i (v) := 1 √ ϑ 1 v i-1 , i = 2, . . . , d + 1 and Ψ d+2 (v) := |v| 2 -dϑ 1 ϑ 1 √ 2d . (1.24)
We deduce formally that h takes the following form

h(t, x, v) = ̺(t, x) + u(t, x) • v + 1 2 θ(t, x)(|v| 2 -dϑ 1 ) M(v) with ̺(t, x) := ˆTd ×R d h(t, x, v) dv, u(t, x) := 1 ϑ 1 ˆTd ×R d h(t, x, v)v dv, θ(t, x) := ˆTd ×R d h(t, x, v) |v| 2 -dϑ 1 ϑ 2 1 d dv. (1.25)
It is worth mentioning that a careful spectral analysis of the linearized collision operator L α defined in (1.16) shows that unless one assumes 1 -α at least of order ε 2 , the eigenfunction associated to the energy dissipation would explode and prevent some exponential stability for (1.15) to hold true (see Theorem 2.1). Actually, in our study, we will require λ 0 to be relatively small with respect to the spectral gap associated to the elastic linearized operator to ensure stability in the inelastic case. If one assumes λ 0 = 0 (for example, one could assume 1 -α of order ε q with q > 2), the effect of the inelasticity is too weak in the hydrodynamic scale and the expected model is the classical Navier-Stokes-Fourier system. In short, we are left with two cases: Case 1: If λ 0 = 0, the expected model is the classical Navier-Stokes-Fourier system. Case 2: If 0 < λ 0 < ∞ is small enough (compared to some explicit quantities), the cumulative effect of inelasticity is visible in the hydrodynamic scale and we expect a different model to the Navier-Stokes-Fourier system accounting for that. In this nearly elastic regime, the energy dissipation rate in the system happens in a controlled fashion since the inelasticity parameter is compensated accordingly to the number of collisions per time unit. Other regimes can be considered depending on the rate at which kinetic energy is dissipated; for example, an interesting regime is the mono-kinetic one which considers the extreme case of infinite energy dissipation rate. In this way, the limit is formally described by enforcing a Dirac mass solution in the kinetic equation yielding the pressureless Euler system (corresponding to sticky particles). Such a regime has been rigorously addressed in the one-dimensional framework in the interesting contribution [START_REF] Jabin | Hydrodynamic limit of granular gases to pressureless Euler in dimension 1[END_REF]. It is an open question to extend such analysis to higher dimensions since the approach of [START_REF] Jabin | Hydrodynamic limit of granular gases to pressureless Euler in dimension 1[END_REF] uses the so-called Bony functional which is a tool specifically tailored for 1D kinetic equations.

1.2. Notations and definitions. Let us introduce some useful notations for functional spaces. For any nonnegative weight function m : R d → R + , we define, for all p > 1 the space L p (m) through the norm

f L p (m) := ˆRd |f (ξ)| p m(ξ) p dξ 1/p ,
We also define, for p 1

W k,p (m) = f ∈ L p (m) ; ∂ β ξ f ∈ L p (m) ∀ |β| k
with the usual norm, i.e., for k ∈ N:

f p W k,p (m) = |β| k ∂ β ξ f p L p (m) .
For m ≡ 1, we simply denote the associated spaces by L p and W k,p . Notice that all the weights we consider here will depend only on velocity, i.e. m = m(v). We will also use the notation ξ := 1 + |ξ| 2 for ξ ∈ R d . On the complex plane, for any a ∈ R, we set

C a := {z ∈ C ; Re z > -a}, C ⋆ a := C a \ {0} (1.26)
and, for any r > 0, we set

D(r) = {z ∈ C ; |z| r}.
We also introduce the following notion of hypo-dissipativity in a general Banach space (X, • ). A closed (unbounded) linear operator A : D(A) ⊂ X → X is said to be hypo-dissipative on X if there exists a norm, denoted by ||| • |||, equivalent to the • -norm such that A is dissipative on the space (X, ||| • |||), that is,

|||(λ -A)h||| λ |||h|||, ∀ λ > 0, h ∈ D(A).
Given two Banach spaces X and Y , we denote with • X→Y the operator norm on the space of B(X, Y ) linear and continuous operators from X to Y . Note also that in what follows, for two positive quantities A and B, we denote by A B if there exists a universal positive constant C (which is in particular independent of the parameters α and ε) such that A CB.

Main results.

The main results are both about the solutions to (1.15). The first one is the following Cauchy theorem regarding the existence and uniqueness of close-to-equilibrium solutions to (1.15). The functional spaces at stake are

L 1 v L 2 x -based Sobolev spaces E 1 ֒→ E defined through E := W k,1 v W m,2 x ( v q ), E 1 := W k,1 v W m,2 x ( v q+1 ) with m > d, m -1 k 0, q 3. ( 1
h ε ∈ C [0, ∞); E ∩ L 1 [0, ∞); E 1
satisfying for any r ∈ (0, 1),

h ε (t) E C η 0 exp (-(1 -r)λ ε t) , ∀ t > 0
for some positive constant C = C(r) > 0 independent of ε and where λ ε ∼ ε→0 λ 0 + η(ε) with λ 0 and η = η(ε)

that have been introduced in Assumption 1.1.

Remark 1.3. It is worth pointing out that the close-to-equilibrium solutions we construct are shown to decay with an exponential rate as close as we want to

λ ε ∼ 1-α(ε) ε 2
(which is the energy eigenvalue of the linearized operator, see Theorem 2.1 hereafter). The rate of convergence can thus be made uniform with respect to the Knudsen number ε (notice that if λ 0 = 0, we obtain a rate of decay as close as we want to η(ε), we thus obtain a uniform bound in time but not a uniform rate of decay).

The estimates on the solution h ε provided by Theorem 1.2 are enough to prove that the solution h ε (t) converges towards some hydrodynamic solution h which depends on (t, x) only through macroscopic quantities (̺(t, x), u(t, x), θ(t, x)) which are solutions to a suitable modification of the incompressible Navier-Stokes system. This is done under an additional assumption on the initial datum that is lightly restrictive. Before stating our main convergence result, we introduce the notation

W ℓ := W ℓ,2 x T d d+2 , ℓ ∈ N
and we furthermore assume that in the definition of the functional spaces (1.27), the following conditions are satisfied:

m > d, m -1 k 1, q 5.
Theorem 1.4. We suppose that the assumptions of Theorem 1.2 are satisfied. We assume furthermore that there exists

(̺ 0 , u 0 , θ 0 ) ∈ W m such that lim ε→0 π 0 h ε in -h 0 L 1 v W m,2
x ( v q ) = 0, where we recall that π 0 is the projection onto the kernel of L 1 defined in (1.23) and

h 0 (x, v) := ̺ 0 (x) + u 0 (x) • v + 1 2 θ 0 (x)(|v| 2 -dϑ 1 ) M(v).
(1.28)

Then, for any T > 0, the family of solutions {h ε } ε constructed in Theorem 1.2 converges in some weak sense to a limit h = h(t, x, v) which is such that

h(t, x, v) = ̺(t, x) + u(t, x) • v + 1 2 θ(t, x)(|v| 2 -dϑ 1 ) M(v), (1.29)
where

(̺, u, θ) ∈ C ([0, T ] ; W m-1 ) ∩ L 2 ((0, T ) ; W m ) is solution to the following incompressible Navier-Stokes-Fourier system with forcing            ∂ t u -ν ϑ1 ∆ x u + ϑ 1 u • ∇ x u + ∇ x p = λ 0 u, ∂ t θ -γ ϑ 2 1 ∆ x θ + ϑ 1 u • ∇ x θ = λ 0 c 2(d + 2) √ ϑ 1 θ, div x u = 0, ̺ + ϑ 1 θ = 0, (1.30)
subject to initial conditions (̺ in , u in , θ in ) defined by

u in := Pu 0 , θ in := d d + 2 θ 0 - 2 (d + 2)ϑ 1 ̺ 0 , ̺ in := -ϑ 1 θ in (1.31)
where P is the Leray projection on divergence-free vector fields and (̺ 0 , u 0 , θ 0 ) have been introduced in (1.28).

The viscosity ν > 0 and heat conductivity γ > 0 are explicit and λ 0 > 0 is the parameter appearing in Assumption 1.1. The parameter c > 0 is depending on the collision kernel b( • ).

Remark 1.5. The data that we consider here are actually quite general. Indeed, the assumption that we make only tells that the macroscopic projection of h in ε converges towards some macroscopic distribution and we do not make any assumption on the macroscopic quantities of this distribution. Namely, we do not suppose that the divergence free and the Boussinesq relations are satisfied by (̺ 0 , u 0 , θ 0 ), the initial layer that could be created by such a lack of assumption is actually absorbed in our notion of weak convergence, the precise notion of which being very peculiar and strongly related to the a priori estimates used for the proof of Theorem 1.2 (see Theorem 4.2 for more details on the type of convergence).

To prove Theorem 1.4, our approach is reminiscent of the program established in [START_REF] Bardos | Fluid dynamic limits of kinetic equations. I. Formal derivations[END_REF][START_REF] Bardos | Fluid dynamic limits of kinetic equations. II. Convergence proofs for the Boltzmann equation[END_REF][START_REF] Golse | The Navier-Stokes limit of the Boltzmann equation for bounded collision kernels[END_REF][START_REF] Saint-Raymond | Hydrodynamic limits of the Boltzmann equation[END_REF] but simpler because our solutions are stronger than the renormalized ones that are used in [START_REF] Golse | The Navier-Stokes limit of the Boltzmann equation for bounded collision kernels[END_REF]. It is based on computations and compactness arguments that were already used in the elastic case. Let us point out that in our case, additional terms appear due to the inelasticity and they can be handled in the framework of Assumption 1.1. In Section 4, we present the proof but only mention its main steps and arguments (details can be found in [START_REF] Alonso | Fluid dynamic limit of Boltzmann equation for granular hard-spheres in a nearly elastic regime[END_REF]Section 6]).

STUDY OF THE KINETIC LINEARIZED PROBLEM

2.1. Main result on the linearized operator. The first step in the proof of Theorem 1.2 is the spectral analysis of the linearized problem associated to (1.15). To that end, we introduce

G α,ε h := - 1 ε v • ∇ x h + 1 ε 2 L α h.
We are going to state our main result on G α,ε in the space E defined in (1.27) but our analysis actually allows to treat even larger spaces (namely, we can obtain the same result under the softer constraints m k 0 and q > 2) but we only state the linear result in this case because it is the only one that will be used in the rest of the paper. Let us also recall that, in any reasonable space (in particular in E and Y j for j = -1, 0, 1 defined in (2.6)-(2.8)), the elastic operator has a spectral gap: there exists µ ⋆ > 0 such that

S(G 1,ε ) ∩ C µ⋆ = {0} (2.1)
where 0 is an eigenvalue of algebraic multiplicity d + 2 of G 1,ε associated to the eigenfunctions

{M, v 1 M, . . . , v d M, |v| 2 M}
(recall that C µ⋆ is defined in (1.26)). This can be proven by an enlargement argument due to [START_REF] Gualdani | Factorization for non-symmetric operators and exponential H-theorem[END_REF] based on the fact that in the Hilbert space

H := W m,2 x,v (M -1/2 ), m > d (2.2)
a result of hypocoercivity has been proven in [START_REF] Briant | From the Boltzmann equation to the incompressible Navier-Stokes equations on the torus: a quantitative error estimate[END_REF] (the constraint m 1 would actually be enough but we will only make use of this result for m > d in the sequel). More precisely, introducing the other Hilbert space

H 1 := W m,2 x,v (M -1/2 v 1/2 ), (2.3)
there exists µ ⋆ > 0 and a norm equivalent to the usual one uniformly in ε (we still denote it by • H and • , • H its associated scalar product to lighten the notations) such that

G 1,ε h, h H - µ ⋆ ε 2 (Id -π 0 )h 2 H1 -µ ⋆ h 2 H1 . (2.4)
As we shall see in the following result, the scaling (1.21) in Assumption 1.1 is precisely the one which allows to preserve exactly d + 2 eigenvalues in the neighborhood of zero for G α,ε . Let us now state our main spectral result (see Figure 1 for an illustration where we have denoted λ ε := -λ d+2 (ε)): Theorem 2.1. Assume that Assumption 1.1 is met. For µ close enough to µ ⋆ defined in (2.1) (in an explicit way), there are some explicit ε > 0 and λ > 0 depending only on χ := µ ⋆ -µ such that, for all ε ∈ (0, ε) and λ 0 ∈ [0, λ), the linearized operator G α(ε),ε has the following spectral property in E:

S(G α(ε),ε ) ∩ C µ = {λ 1 (ε), . . . , λ d+2 (ε)}, (2.5
)

with λ 1 (ε) = 0, λ j (ε) = 1 -α(ε) ε 2 , j = 2, . . . , d + 1, and
λ d+2 (ε) = - 1 -α(ε) ε 2 + O(ε 2 ) as ε → 0. O D(µ⋆ -µ) χ χ -µ ⋆ -µ Reλ Imλ -λε Cµ \ D(µ⋆ -µ) FIGURE 1. The set C µ \ D(µ ⋆ -µ) and the eigenvalue -λ ε . Remark 2.2.
It is worth noticing that the eigenvalue λ 1 (ε) = 0 corresponds with the property of mass conservation of the operator G α(ε),ε . Concerning the intermediate eigenvalues λ j (ε) for j = 2, . . . , d+ 1, as one can see on their definition, they may be positive, this is due to the fact that the collision operator Q α preserves momentum while the drift operator ∇ v (v • ) does not. However, using (1.19), one can prove that the vanishing momentum is preserved by the whole operator G α(ε),ε , consequently, those eigenvalues won't affect the long-time analysis of our problem. Finally, the eigenvalue λ d+2 (ε) is directly linked to the non-preservation of energy property of G α(ε),ε .

We are going to prove Theorem 2.1 in two stages. First, we perform a perturbative argument (reminiscent of [START_REF]Boltzmann equation for granular media with thermal force in a weakly inhomogeneous setting[END_REF]) in a L 2 v,x -based Sobolev space, namely in

Y := W s,2 v W ℓ,2 x ( v r ), ℓ ∈ N, s ∈ N * , ℓ s + 1, r > r ⋆ + κ + 2 (2.6)
where

r ⋆ := 4 σ 1 σ 0 + 3 2
with σ 0 and σ 1 defined in (2.9) and κ > d/2. The key point of our approach is to see G α,ε as a perturbation of the elastic linearized operator G 1,ε . We then use an enlargement argument (from [START_REF] Gualdani | Factorization for non-symmetric operators and exponential H-theorem[END_REF]) to extend the result from Y to the space E defined in (1.27). Several remarks are in order:

(i) First, let us remind that the global equilibrium of our equation G α defined in (1.12) has some exponential fat tail and in particular, decays more slowly than a standard Maxwellian distribution (see [START_REF] Mouhot | Stability, convergence to self-similarity and elastic limit for the Boltzmann equation for inelastic hardspheres[END_REF]). As a consequence, we can not rely on classical works on the elastic linearized operator which are developed in spaces of type L 2 v,x (M -1/2 ) with M defined in (1.14). To overcome this difficulty, we exploit results coming from [START_REF] Gualdani | Factorization for non-symmetric operators and exponential H-theorem[END_REF] in which an enlargement theory has been carried out. The results proven in [START_REF] Gualdani | Factorization for non-symmetric operators and exponential H-theorem[END_REF] include a spectral analysis of the elastic Boltzmann operator G 1,1 in larger spaces (in particular of type L 2 v,x ) with "soft weights" that can be polynomial or stretched exponential. In the same line of ideas, these results have been extended to the rescaled elastic operator G 1,ε in [START_REF] Briant | From Boltzmann to incompressible Navier-Stokes in Sobolev spaces with polynomial weight[END_REF]. (ii) Let us also point out that the perturbation at stake does not fall into the realm of the classical perturbation theory of unbounded operators as described in [START_REF] Kato | Perturbation theory for linear operators[END_REF] because the perturbation is not relatively bounded. Indeed, the domain of

G 1,ε in Y is given by W s+1,2 v W ℓ+1,2
x ( v r+1 ) while if one wants to be sharp in terms of rate, the best estimate in terms of functional spaces that we are able to get is

G α,ε -G 1,ε Yj →Yj-1 = 1 ε 2 L α -L 1 Yj →Yj-1 1 -α ε 2 , j = 0, 1, (2.7)
where the spaces Y j are defined through

Y -1 := W s-1,2 v W ℓ,2 x ( v r-κ-2 ), Y 0 := Y, Y 1 := W s+1,2 v W ℓ,2 x ( v r+κ+2 ) (2.8)
with κ > d/2. These estimates (whose proofs can be found in [START_REF] Alonso | Fluid dynamic limit of Boltzmann equation for granular hard-spheres in a nearly elastic regime[END_REF]Lemma 3.3]) are a generalization and optimisation of estimates obtained in [START_REF] Mouhot | Stability, convergence to self-similarity and elastic limit for the Boltzmann equation for inelastic hardspheres[END_REF] and are sharp in terms of rate, this sharpness being needed in our analysis since it allows us to deal with the case λ 0 > 0 in Assumption 1.1. (iii) As a consequence, we have to use refined perturbation arguments whose key insights come from [START_REF]Boltzmann equation for granular media with thermal force in a weakly inhomogeneous setting[END_REF].

Note however that we drastically simplify the analysis performed in [START_REF]Boltzmann equation for granular media with thermal force in a weakly inhomogeneous setting[END_REF] by remarking that the difference operator G α,ε -G 1,ε does not involve any spatial derivative and that we "only" need to develop a spectral analysis of G α,ε without being able to obtain decay properties on the associated semigroup. As a consequence, we don't need to use a spectral mapping theorem, nor do we need to use an iterated version of Duhamel formula and this is crucial in order to reach the optimal scaling (1.21) for our restitution coefficient. (iv) Let us finally mention that we perform our perturbative argument in Y which is a

L 2 v,x -based Sobolev space instead of performing it in E (which is L 1 v L 2 x -based) directly. This intermediate step seems necessary because even if L α -L 1 satisfies nice estimates in L 1
v , the use of Fubini theorem is actually crucial to get the rate (1 -α)/ε 2 in estimates of type (2.7).

Elements of proof of Theorem 2.1.

As mentioned above, the basis of the proof of this theorem is to see L α as a perturbation of L 1 .

We start by giving a splitting of it into two parts: one which has some good regularizing properties (in the velocity variable) and another one which is dissipative. For any δ > 0, one can write

L 1 = A (δ) + B (δ) 1
with A (δ) and B (δ) 1 defined through an appropriate mollification-truncation process (see [START_REF] Gualdani | Factorization for non-symmetric operators and exponential H-theorem[END_REF]Section 4.3.3] and [3, Section 2.2] for the details). The elastic collision operator L 1 writes (see the strong formulation of Q 1 in (1.6)):

L 1 g = ˆRd ×S d-1 b(σ • q)|v -v * |(M ′ * g ′ + M ′ g ′ * -Mg * ) dσ dv * - ˆRd M * |v -v * | dv * g
where we have used the shorthand notations g

= g(v), g * = g(v * ), g ′ = g(v ′ ), g ′ * = g(v ′ * ). We define A (δ) g := ˆRd ×S d-1 Θ δ b(σ • q)|v -v * |(M ′ * g ′ + M ′ g ′ * -Mg * ) dσ dv * B (δ) 1 g := ˆRd ×S d-1 (1 -Θ δ ) b(σ • q)|v -v * |(M ′ * g ′ + M ′ g ′ * -Mg * ) dσ dv * -g ˆRd M * |v -v * | dv * where Θ δ = Θ δ (v, v * , σ) is an appropriate truncature function. The dissipativity property of B (δ)
1 comes from the fact that the truncature function Θ δ is defined so that the first term is small as δ goes to 0 and the fact that there exist σ 0 > 0 and σ 1 > 0 such that

σ 0 v ˆRd M * |v -v * | dv * σ 1 v , v ∈ R d . (2.9)
As a consequence, B

1 is going to be dissipative for δ small enough. This leads to the following decomposition of L α :

L α = B (δ) α + A (δ) , where B (δ) α := B (δ) 1 dissipative + [L α -L 1 ]
small as α → 1

(2.10) and then the following decomposition of G α,ε :

G α,ε = A (δ) ε + B (δ) α,ε , where
A (δ) ε := 1 ε 2 A (δ) , B (δ) α,ε := 1 ε 2 B (δ) α - 1 ε v • ∇ x . (2.11)
Our analysis of this splitting and then of the spectrum of G α,ε relies on several elements: the nice properties of the above-mentioned splitting of L 1 = A (δ) + B (δ) 1 coming from [START_REF] Gualdani | Factorization for non-symmetric operators and exponential H-theorem[END_REF], some refined bilinear estimates on the collision operator coming from [START_REF] Alonso | Convolution inequalities for the Boltzmann collision operator[END_REF], new estimates on G α -M that are reminiscent of estimates proven in [START_REF] Alonso | Uniqueness and regularity of steady states of the Boltzmann equation for viscoelastic hard-spheres driven by a thermal bath[END_REF] (see [START_REF] Alonso | Fluid dynamic limit of Boltzmann equation for granular hard-spheres in a nearly elastic regime[END_REF]Lemma 2.3]) and also new estimates on Q α -Q 1 (see [3, Lemmas 2.1 and 2.2]). Concerning the latter point, we exploit ideas developed in [START_REF] Mouhot | Stability, convergence to self-similarity and elastic limit for the Boltzmann equation for inelastic hardspheres[END_REF] but our situation is more involved because we work in polynomially weighted spaces whereas in [START_REF] Mouhot | Stability, convergence to self-similarity and elastic limit for the Boltzmann equation for inelastic hardspheres[END_REF], the authors were working with stretched exponential weights.

In the following lemma, we provide some regularization and hypodissipativity results on the splitting

G α,ε = A (δ) ε + B (δ)
α,ε (see [3, Lemma 2.7 and Proposition 2.9]): Lemma 2.3. There holds: [START_REF] Alonso | Convolution inequalities for the Boltzmann collision operator[END_REF] For any k ∈ N and δ > 0, there are two positive constants

C k,δ , R δ > 0 such that supp A (δ) g ⊂ B(0, R δ )
and

A (δ) g W k,2 v (R d ) C k,δ g L 1 v ( v ) , ∀ g ∈ L 1 v ( v ). ( 2

.12)

(2) There exist δ 0 , α 0 , ν 0 such that for all α ∈ (α 0 , 1) and δ ∈ (0, δ 0 ),

B (δ) α,ε + ε -2 ν 0 is hypo-dissipative in E and Y j , j = -1, 0, 1
, where we recall that the spaces E and Y j are respectively defined in (1.27) and (2.8).

In what follows, we suppose that Assumption 1.1 is satisfied. We introduce ε 0 which is such that α(ε 0 ) = α 0 (and thus α(ε) ∈ (α 0 , 1) for all ε ∈ (0, ε 0 )) and consider δ ∈ (0, δ 0 ), ε ∈ (0, ε 0 ). We will denote G ε := G α,ε as well as A ε := A (1) For any i, j ∈ {-1, 0, 1}, we have

A ε Yi→Yj 1 ε 2 .
(2) If ν > 0 is fixed, then for ε small enough (in terms of ν 0 and ν) and j = -1, 0, 1,

R(λ, B ε ) Yj→Yj 1 Re λ + ε -2 ν 0 ε 2 , ∀ Re λ > -ν .
The second keypoint to develop our perturbative argument is to have a good understanding of the spectrum of the operator G 1,ε . We here give some estimates on the associated resolvent that are a consequence of a result of decay of the associated semigroup (see [START_REF] Alonso | Fluid dynamic limit of Boltzmann equation for granular hard-spheres in a nearly elastic regime[END_REF]Theorem 2.12] which gives an improved version of [8, Theorem 2.1]): Lemma 2.5. There exists ε 1 ∈ (0, ε 0 ) such that for j = -1, 0, 1, for any λ ∈ C ⋆ µ⋆ and any ε ∈ (0, ε 1 ),

R(λ, G 1,ε ) Yj→Yj max 1 |λ| , 1 Reλ + µ ⋆
where µ ⋆ has been defined in (2.1).

Let us now explain how we develop our perturbative argument to prove Theorem 2.1. The following proposition (which is an adaptation of [START_REF]Boltzmann equation for granular media with thermal force in a weakly inhomogeneous setting[END_REF]Lemma 2.16]) is the first step in the development of the perturbative argument and its proof relies on Corollary 2.4 and Lemma 2.5.

Proposition 2.6. For all

λ ∈ C ⋆ µ⋆ , let J ε (λ) = (G ε -G 1,ε ) R(λ, G 1,ε )A ε R(λ, B ε ).
Then, for any µ ∈ (0, µ ⋆ ) and λ ∈ C µ \ D(µ ⋆ -µ), there exists ε 2 ∈ (0, ε 1 ) such that for any ε ∈ (0, ε 2 ),

J ε (λ) Y→Y 1 µ ⋆ -µ 1 -α(ε) ε 2 .
(2.13)

In addition, there exists ε 3 ∈ (0, ε 2 ) and λ 3 > 0 such that for ε ∈ (0, ε 3 ) and λ 0 ∈ [0, λ 3 ) (where λ 0 is defined in Assumption 1.1), Id -J ε (λ) and λ -G ε are invertible in Y with

R(λ, G ε ) = Γ ε (λ)(Id -J ε (λ)) -1 , λ ∈ C µ \ D(µ ⋆ -µ), (2.14)
where

Γ ε (λ) := R(λ, B ε ) + R(λ, G 1,ε )A ε R(λ, B ε ). Finally, we have for ε ∈ (0, ε 3 ), R(λ, G ε ) Y→Y 1 µ ⋆ -µ , λ ∈ C µ \ D(µ ⋆ -µ). (2.15)
Sketch of the proof. The estimate on J ε (λ) can be easily deduced from (2.7), Corollary 2.4 and Lemma 2.5. First, fix µ ∈ (0, µ ⋆ ) and notice that from Corollary 2.4, we clearly have that there exists ε 2 ∈ (0, ε 1 ) (which depends on ν 0 and µ) such that for any Re λ > -µ, we have:

A ε R(λ, B ε ) Y→Y1 1.
We can then deduce that for µ ∈ (0, µ ⋆ ), for any Re λ > -µ, |λ| µ ⋆ -µ,

J ε (λ) Y→Y G ε -G 1,ε Y1→Y R(λ, G 1,ε ) Y1→Y1 A ε R(λ, B ε ) Y→Y1 C 1 -α(ε) ε 2 max 1 |λ| , 1 Reλ + µ ⋆ C 1 -α(ε) ε 2 1 µ ⋆ -µ (2.16)
for some C > 0. Moreover, one can choose ε 3 ∈ (0, ε 2 ) and λ 3 > 0 depending on the difference χ = µ ⋆ -µ, so that if λ 0 ∈ [0, λ 3 ) (recall that λ 0 is defined in Assumption 1.1 and is such that

(1 -α(ε))ε -2 ∼ λ 0 + η(ε)) ρ(ε) := C µ ⋆ -µ 1 -α(ε) ε 2 < 1, ∀ ε ∈ (0, ε 3 ). (2.17)
Under such an assumption, one sees that, for all

λ ∈ C µ \ D(µ ⋆ -µ), Id -J ε (λ) is invertible in Y with (Id -J ε (λ)) -1 = ∞ p=0 [J ε (λ)] p , ∀ ε ∈ (0, ε 3 ). Let us fix then ε ∈ (0, ε 3 ) and λ ∈ C µ \ D(µ ⋆ -µ). The range of Γ ε (λ) is clearly included in D(B ε ) = D(G 1,ε ).
Then, writing G ε = A ε + B ε , we easily get that

(λ -G ε )Γ ε (λ) = Id -J ε (λ) i.e. Γ ε (λ)(Id -J ε (λ)) -1 is a right-inverse of (λ -G ε ).
To prove that λ -G ε is invertible, it is therefore enough to prove that it is one-to-one, which can be done up to reducing the value of ε 3 using (2.7), Lemmas 2.3 and 2.5. Thus, for ε ∈ (0, ε 3 ), C µ \ D(µ ⋆ -µ) is included into the resolvent set of G ε and this shows (2.14).

To estimate now R(λ, G ε ) B(Y) , one simply notices that

(Id -J ε (λ)) -1 Y→Y ∞ p=0 J ε (λ) p Y→Y 1 1 -ρ(ε) , ∀ λ ∈ C µ \ D(µ ⋆ -µ) (2.18)
from which, as soon as

λ ∈ C µ \ D(µ ⋆ -µ), R(λ, G ε ) Y→Y 1 1 -ρ(ε) Γ ε (λ) Y→Y .
One checks, using the previous computations, that for λ

∈ C µ \ D(µ ⋆ -µ), Γ ε (λ) Y→Y ε 2 + A Y→Y R(λ, G 1,ε ) Y→Y (2.19)
and deduces (2.15). This achieves the proof.

A first obvious consequence of Proposition 2.6 is that, for any µ ∈ (0, µ ⋆ ), there is

ε 3 > 0 depending only on χ = µ ⋆ -µ such that S(G ε ) ∩ C µ ⊂ D(µ ⋆ -µ), ∀ ε ∈ (0, ε 3 ).
We denote by P ε (resp. P 0 ) the spectral projection associated to the set

S(G ε ) ∩ C µ = S(G ε ) ∩ D(µ ⋆ -µ) (resp. S(G 1,ε ) ∩ C µ = {0}).
One can deduce then the following lemma whose proof is similar to [21, Lemma 2.17].

Lemma 2.7. For any µ ∈ (0, µ ⋆ ) such that C µ ⊂ D(µ ⋆ -µ), there exist ε 4 ∈ (0, ε 3 ) and λ 4 ∈ (0, λ 3 ) depending only on χ = µ ⋆ -µ such that if λ 0 ∈ [0, λ 4 ) (where λ 0 is defined in Assumption 1.1),

P ε -P 0 Y→Y < 1, ∀ ε ∈ (0, ε 4 ).
In particular, dim Range(P ε ) = dim Range(P 0 ) = d + 2, ∀ ε ∈ (0, ε 4 ).

(2.20)

Sketch of the proof. Let µ ∈ (0, µ ⋆ ) be close enough to µ ⋆ so that D(µ ⋆ -µ) ⊂ C µ and 0 < r < χ = µ ⋆ -µ.

One has D(r) ⊂ C ⋆ µ . We set γ r := {z ∈ C ; |z| = r}. Recall that by definition

P ε := 1 2iπ ˛γr R(λ, G ε ) dλ, P 0 := 1 2iπ ˛γr R(λ, G 1,ε ) dλ. For λ ∈ γ r , set Z ε (λ) = R(λ, G 1,ε )A ε R(λ, B ε ) so that Γ ε (λ) = R(λ, B ε ) + Z ε (λ). Recall from (2.14) that, for λ ∈ γ r , R(λ, G ε ) = R(λ, B ε )(Id -J ε (λ)) -1 + Z ε (λ)(Id -J ε (λ)) -1 = R(λ, B ε ) + R(λ, B ε )J ε (λ)(Id -J ε (λ)) -1 + Z ε (λ)(Id -J ε (λ)) -1
where we wrote (Id -J ε (λ)) -1 = Id + J ε (λ)(Id -J ε (λ)) -1 to get the second equality. One also has

R(λ, G 1,ε ) = R(λ, B 1,ε ) + R(λ, G 1,ε )A ε [R(λ, B 1,ε ) -R(λ, B ε )] + Z ε (λ).
One can then obtain (see the proof of [3, Lemma 3.8] for the details)

P ε -P 0 = 1 2iπ ˛γr Γ ε (λ)J ε (λ)(Id -J ε (λ)) -1 dλ + 1 2iπ ˛γr R(λ, G 1,ε )A ε [R(λ, B ε ) -R(λ, B 1,ε )] dλ.
The first part is estimated thanks to (2.16), (2.18) and (2.19) combined with Lemma 2.5:

Γ ε (λ)J ε (λ)(Id -J ε (λ)) -1 Y→Y 1 r 2 1 1 -ρ(ε) 1 -α(ε) ε 2 .
For the second part, notice first that from Lemma 2.5,

R(λ, G 1,ε )A ε [R(λ, B ε ) -R(λ, B 1,ε )] Y→Y 1 r A ε R(λ, B ε ) -A ε R(λ, B 1,ε ) Y→Y .
Then, for λ ∈ γ r , we have

A ε R(λ, B ε ) -A ε R(λ, B 1,ε ) = A ε R(λ, B ε ) [B ε -B 1,ε ] R(λ, B 1,ε ) which implies that A ε R(λ, B ε ) -A ε R(λ, B 1,ε ) Y→Y A ε R(λ, B ε ) Y-1→Y B ε -B 1,ε Y→Y-1 R(λ, B 1,ε ) Y→Y 1 -α(ε) ε 2 .
Proceeding as in the proof of Lemma 2.6, one can conclude that for any 0 < r < χ = µ ⋆ -µ, With Lemma 2.7, we can now end the proof of Theorem 2.1. Sketch of the proof of Theorem 2.1. The structure of S(G ε ) ∩ C µ in the space Y comes directly from Lemma 2.7 together with Proposition 2.6. To describe more precisely the spectrum, one first remarks that

P ε -P 0 Y→Y C r 1 -α(ε) ε 2 1 r(1 -ρ(ε)) + 1 := ℓ(ε). ( 2
S(L α(ε) ) ∩ C µ ⊂ S(G ε ) ∩ C µ .
This comes from the fact that for each eigenvalue of L α(ε) , the eigenfunction depends only on v and thus remains an eigenfunction for the operator G ε . Since, for ε small enough, the same perturbative argument that we developed above implies that the spectral projection

Π L α(ε) associated to S(L α ) ∩ C µ satisfies dim(Range(Π L α(ε) )) = dim(Range(Π L1 )) = d + 2 = dim(Range(P ε )), we get that S(L α(ε) ) ∩ C µ = S(G ε ) ∩ C µ , (2.22)
that is, the eigenvalues λ j (ε) are actually eigenvalues of L α(ε) . The development of the energy eigenvalue λ d+2 (ε) comes from [START_REF] Mouhot | Stability, convergence to self-similarity and elastic limit for the Boltzmann equation for inelastic hardspheres[END_REF]. The conservation of mass gives us that 0 is an eigenvalue for our problem. The intermediate eigenvalues λ j (ε) for j = 2, . . . , d + 1 are obtained thanks to the fact that

ˆRd L α(ε) ϕ(v) v i dv = - 1 -α(ε) ε 2 ˆRd v i ∇ • (vϕ(v)) dv = 1 -α(ε) ε 2 ˆRd v i ϕ(v) dv.
Notice that all this allows us to find eigenfunctions (that depend only on v) in L 2 v,x ( v r ). Using once more the splitting L α = A (δ) + B δ α defined in (2.10) and the regularizing properties of A (δ) , one can actually prove that our eigenfunctions lie in Y, which yields the conclusion of Theorem 2.1 in the space Y. To extend the result to the space E, we use an enlargement argument coming from [START_REF] Gualdani | Factorization for non-symmetric operators and exponential H-theorem[END_REF], we omit the details here and just mention that this argument is based on the splitting G ε = A ε + B ε introduced in (2.11).

STUDY OF THE KINETIC NONLINEAR PROBLEM

Let us recall that the spaces E and E 1 are defined in (1.27). In this section, we assume that Assumption 1.1 is met and consider ε ∈ (0, ε), λ 0 ∈ 0, λ where ε and λ are defined in Theorem 2.1. As in Section 2, to lighten the notations, we write G ε = G α(ε),ε as well as B ε = B α(ε),ε .

3.1. Splitting of the nonlinear inelastic Boltzmann equation. Now that the spectral analysis of the linearized operator G ε in the space E has been performed, in order to prove Theorem 1.2, we are going to prove several a priori estimates for the solutions to (1.15). The crucial point in the analysis lies in the splitting of (1.15) into a system of two equations mimicking a spectral enlargement method from a PDE perspective (see [START_REF] Mouhot | Exponential stability of slowly decaying solutions to the kinetic-Fokker-Planck equation[END_REF]Section 2.3] and [START_REF] Briant | From Boltzmann to incompressible Navier-Stokes in Sobolev spaces with polynomial weight[END_REF] for pioneering ideas on such a method). More precisely, using (2.11), the splitting amounts to look for a solution of (1.15) of the form

h ε (t) = h 0 ε (t) + h 1 ε (t) with h 0 ε solution to            ∂ t h 0 ε = B ε h 0 ε + 1 ε Q α(ε) (h 0 ε , h 0 ε ) + 1 ε Q α(ε) (h 0 ε , h 1 ε ) + Q α(ε) (h 1 ε , h 0 ε ) + G ε h 1 ε -G 1,ε h 1 ε + 1 ε Q α(ε) (h 1 ε , h 1 ε ) -Q 1 (h 1 ε , h 1 ε ) , h 0 ε (0, x, v) = h in ε (x, v) ∈ E, (3.1)
and

h 1 ε solution to    ∂ t h 1 ε = G 1,ε h 1 ε + 1 ε Q 1 (h 1 ε , h 1 ε ) + A ε h 0 ε , h 1 ε (0, x, v) = 0. (3.2)
In order to lighten the notations, in this section, we will write h in , h, h 0 and h 1 instead of h in ε , h ε , h 0 ε and h 1 ε . The goal is to obtain nice nested a priori estimates on h 0 and h 1 . Notice first that our splitting is more complicated than the one of [START_REF] Briant | From Boltzmann to incompressible Navier-Stokes in Sobolev spaces with polynomial weight[END_REF] because it relies on perturbative considerations around the elastic case that come out in the equation satisfied by h 0 . As a consequence, our a priori estimates are more intricate and require the use of non standard Gronwall lemma. Notice also that since the initial datum of h 1 is vanishing, we can study the equation on h 1 in any functional space. In particular, we can study it in the Hilbert space H = W m,2

x,v M -1/2 in which we have a good understanding of the elastic linearized operator G 1,ε . Indeed, in this type of spaces, the symmetries of the collision operator Q 1 allow to get some nice hypocoercive estimates (see (2.4)).

Remark 3.1. In [START_REF] Gualdani | Factorization for non-symmetric operators and exponential H-theorem[END_REF], the authors treat the elastic case (α = 1) of the non-rescaled equation (ε = 1) and they do not resort to such a splitting method to study the nonlinear equation, their approach is based on the use of a norm which is equivalent to the usual one and is such that G 1,1 is dissipative in this norm in large spaces. Such an approach is no longer usable when one wants to deal with rescaled equations and obtain uniform in ε estimates. Indeed, the definition of the equivalent norm in [START_REF] Gualdani | Factorization for non-symmetric operators and exponential H-theorem[END_REF] does not take into account the different behaviors of microscopic and macroscopic parts of the solution with respect to ε: typically, the microscopic part of the solution vanishes as ε → 0 whereas the macroscopic one does not. Conversely, in the splitting method, the equation that defines h 1 is treated thanks to hypocoercivity tricks that allow to distinguish microscopic and macroscopic behaviors.

3.2. Estimating h 0 . Concerning h 0 , let us first mention that the dissipativity properties of B ε stated in Lemma 2.3 can actually be improved a bit. More precisely, one can show that there exist norms on the spaces E and E 1 that are equivalent to the standard ones (with multiplicative constants independent of ε) that we still denote • E and • E1 and that satisfy:

d dt S Bε (t)g E - ν 0 ε 2 S Bε (t)g E1 (3.3)
where we have denoted by (S Bε (t)) t 0 the semigroup generated by B ε and ν 0 is defined in Lemma 2.3. Let us also introduce the Banach space E 2

E 2 := W k+1,2 v W m,2 x (̟ q+2κ+2 ), κ > d 2 which satisfies the following continuous embeddings: H ֒→ E 2 ֒→ E 1 (recall that E 1 is defined in (1.27)).
Let us point out that the spaces E 1 and E 2 allow us to get the following estimates (see [START_REF] Alonso | Fluid dynamic limit of Boltzmann equation for granular hard-spheres in a nearly elastic regime[END_REF]Remark 3.5] and [START_REF] Alonso | Convolution inequalities for the Boltzmann collision operator[END_REF]2]):

(Q α -Q 1 )(g, f ) E (1 -α) g E2 f E2 and Q α (g, f ) E g E f E1 + g E1 f E (3.4)
where the multiplicative constants are uniform in α. One can then obtain the following proposition:

Proposition 3.2. Assume that h 0 ∈ E, h 1 ∈ H are such that sup t 0 h 0 (t) E + h 1 (t) H ∆ 0 < ∞.
For ν ∈ (0, ν 0 ) (where ν 0 is defined in Lemma 2.3), there exists an explicit ε 5 ∈ (0, ε) (where ε is defined in Theorem 2.1) such that:

h 0 (t) E h in E e -ν ε 2 t + λ ε ˆt 0 e -ν ε 2 (t-s) h 1 (s) H ds (3.5)
where we recall that

λ ε ∼ ε→0 1-α(ε) ε 2 is defined in Theorem 2.1.
Sketch of the proof. Using (3.3) as well as (3.4) and recalling that h 0 solves (3.1), we can compute the evolution of h 0 (t) E and estimate it:

d dt h 0 (t) E - ν 0 ε 2 h 0 (t) E1 + C ε h 0 (t) E + h 1 (t) E1 h 0 (t) E1 + C 1 -α(ε) ε 2 h 1 (t) E2 + C 1 -α(ε) ε h 1 (t) 2 E2 . (3.6)
Using that the embedding E 2 ֒→ H is continuous, recalling that h 0 (0) = h in and choosing ε 5 small enough so that C ε 5 ∆ 0 ν 0 -ν, we obtain

h 0 (t) E h in E e -ν ε 2 t + λ ε ˆt 0 e -ν ε 2 (t-s) h 1 (s) H ds + ε λ ε ˆt 0 e -ν ε 2 (t-s) h 1 (s) 2 H ds.
We conclude to (3.5) by assuming furthermore that ε 5 ∆ 0 1.

3.3.

Estimating h 1 . We now comment and study the equation satisfied by h 1 . Let us point out that getting estimates on h 1 is trickier than in [START_REF] Briant | From Boltzmann to incompressible Navier-Stokes in Sobolev spaces with polynomial weight[END_REF], indeed, in the latter paper, the idea is to estimate separately P 0 h 1 and (Id -P 0 )h 1 where P 0 is the projector onto Ker(G 1,ε ) defined by

P 0 g := d+2 i=1 ˆTd ×R d g Ψ i dv dx Ψ i M (3.7)
where the functions Ψ i have been defined in (1.24), and thanks to the properties of preservation of mass, momentum and energy of the whole equation, one could write that P 0 h = 0 so that P 0 h 1 = -P 0 h 0 and directly get an estimate on P 0 h 1 from the one on h 0 . In our case, the energy is no longer preserved which induces additional difficulties. However, we keep the same strategy and start by estimating P 0 h 1 (see Remark 3.4 for a comment on this choice of strategy).

For the sequel, we also introduce

P 0 h = d+1 i=1 ˆTd ×R d h Ψ i dv dx Ψ i M , Π 0 h = ˆTd ×R d hΨ d+2 dv dx Ψ d+2 M. (3.8)
Notice that thanks to Cauchy-Schwarz inequality in velocity, one can easily prove that we have P 0 ∈ B(E, H).

One can then obtain the following proposition:

Proposition 3.3. Assume that h 0 ∈ E, h 1 ∈ H are such that sup t 0 h 0 (t) E + h 1 (t) H ∆ 0 < ∞. For ε ∈ (0, ε 5 ) (ε 5 is defined in Proposition 3.2), P 0 h 1 (t) E h 0 (t) E + h in E e -λεt + λ ε ˆt 0 e -λε(t-s) h 0 (s) E + (Id -P 0 )h 1 (s) H ds + ελ ε ˆt 0 e -λε(t-s) h 1 (s) H ds (3.9)
where

λ ε = λ ε + O(1 -α(ε)) with λ ε ∼ ε→0 1-α(ε) ε 2 defined in Theorem 2.1.
Sketch of the proof. Due to the properties of preservation of mass and vanishing momentum of our equation, we have P 0 h = 0 which implies that P 0 h 1 = -P 0 h 0 . Consequently, we easily get an estimate on P 0 h 1 using that P 0 ∈ B(E, H):

P 0 h 1 (t) H h 0 (t) E .
(3.10) It now remains to estimate Π 0 h 1 . To this end, we first notice that

Π 0 h 1 = P 0 h 1 -P 0 h 1 = P 0 h -P 0 h 0 -P 0 h 1 = Π 0 h -P 0 h 0 -P 0 h 1
where we used that P 0 h = Π 0 h due to the preservation of mass and vanishing momentum so, using (3.5) and (3.10), we only need to estimate Π 0 h to get an estimate on P 0 h 1 . To this end, we start by computing the evolution of Π 0 h:

∂ t (Π 0 h) = Π 0 (G ε h) + 1 ε Π 0 Q α(ε) (h, h).
By direct inspection, using the definition of Π 0 given in (3.8) and the dissipation of energy (1.7) (see [3, Lemmas 4.2 and 4.5]), we obtain: as ε → 0,

Π 0 (G ε h) = -λ ε Π 0 h + O 1 -α(ε) ε 2 (Id -P 0 )h E with λ ε = λ ε + O(1 -α(ε)) ∼ ε→0 λ ε .
Similarly, we have by direct computation that

|Π 0 Q α(ε) (h, h)| = (1 -α 2 )|D α(ε) (h, h)|Ψ d+2 M
so that, using Minkowski's inequality to estimate D α(ε) (h, h), we obtain

Π 0 Q α(ε) (h, h) E ε 2 h 2 E . (3.11)
Gathering previous estimates, we are able to deduce that

P 0 h 1 (t) E h 0 (t) E + h in E e -λεt + λ ε ˆt 0 e -λε(t-s) h 0 (s) E + (Id -P 0 )h 1 (s) H ds + ελ ε ˆt 0 e -λε(t-s) h 0 (s) 2 E + h 1 (s) 2 H ds.
With this, inequality (3.9) holds by using ε 5 ∆ 0 1 from the proof of Proposition 3.2.

Remark 3.4. A natural approach would have been to adapt the method of [START_REF] Briant | From Boltzmann to incompressible Navier-Stokes in Sobolev spaces with polynomial weight[END_REF] by applying P ε (the projector associated to the eigenvalues λ j (ε) for j = 1, . . . , d + 2 of G ε around 0 that have been exhibited in Theorem 2.1) to our equation instead of P 0 . It implies that one would have had to estimate Π ε h where Π ε is the projector associated to the energy eigenvalue -λ ε = λ d+2 (ε) defined in Theorem 2.1. On the one hand, it simplifies the approach because Π ε G ε h = -λ ε Π ε h by definition. On the other hand, this projector is not explicit contrary to Π 0 and when applying Π ε to the equation satisfied by h

∂ t h = G ε h + 1 ε Q α(ε) (h, h), nothing guarantees that Π ε ε -1 Q α(ε) (h, h
) remains of order 1 with respect to ε whereas we have seen in (3.11) that due to the dissipation of kinetic energy,

Π 0 ε -1 Q α(ε) (h, h
) is actually of order ε. This explains our choice of strategy.

Let us now focus on the estimate of (Id -P 0 )h 1 . We can proceed similarly as in [START_REF] Briant | From Boltzmann to incompressible Navier-Stokes in Sobolev spaces with polynomial weight[END_REF], using in particular that P 0 Q 1 = 0. Another crucial point is that the source term A ε h 0 can be bounded in H using the fact that A ε ∈ B(E, H) (see Lemma 2.3). Moreover, it is important to mention that the fact that the bound on A ε induces a rate of ε -2 will be counterbalanced by the fact that the semigroup associated with B ε has an exponential decay rate of type e -νt/ε 2 (see (3.3)). We recall that the Hilbert space H 1 is defined in (2.3) and is such that

Q 1 (g, g) H g H g H1 . (3.12)
Proposition 3.5. Assume that h 0 ∈ E, h 1 ∈ H are such that

sup t 0 h 0 (t) E + h 1 (t) H ∆ 0 < ∞.
For µ ∈ (0, µ ⋆ ) (where µ ⋆ is defined in (2.1)) and for ∆ 0 small enough, we have that:

(Id -P 0 )h 1 (t) 2 H ∆ 2 0 ˆt 0 e -µ(t-s) h 1 (s) 2 H ds + 1 ε 2 ˆt 0 e -µ(t-s) h 1 (s) H h 0 (s) E ds. (3.13)
Sketch of the proof. From (3.2), the fact that P 0 Q 1 (g, g) = 0 and the fact that P 0 commutes with G 1,ε , we can compute the evolution of Φ(t) := (Id -P 0 )h 1 :

∂ t Φ = G 1,ε Φ + 1 ε Q 1 (h 1 , h 1 ) + (Id -P 0 )A ε h 0 .
We now use the hypocoercive norm on H for G 1,ε introduced in (2.4) and also denote by Φ ⊥ the microscopic part of Φ, namely Φ ⊥ := (Idπ 0 )Φ where we recall that π 0 is the projection onto the kernel of L 1 that has been introduced in (1.23). We compute the evolution of

Φ(t) 2 H : 1 2 d dt Φ(t) 2 = G 1,ε Φ(t), Φ(t) H + 1 ε Q 1 (h 1 (t), h 1 (t)), Φ ⊥ (t) H + (Id -P 0 )A ε h 0 (t), Φ(t) H .
Notice that we have been able to replace Φ by Φ ⊥ in the second term due to the conservation laws satisfied by Q 1 and the fact that π 0 is orthogonal in H. Then, from the properties of the hypocoercive norm (see (2.4)), using (3.12) and the facts that P 0 ∈ B(H), A ∈ B(E, H) (from Lemma 2.3) as well as Cauchy-Schwarz inequality, we obtain that

1 2 d dt Φ(t) 2 H - µ ⋆ ε 2 Φ ⊥ (t) 2 H1 -µ ⋆ Φ(t) 2 H1 + C ε h 1 (t) H h 1 (t) H1 Φ ⊥ (t) H + C ε 2 h 0 (t) E Φ(t) H .
Making an appropriate use of Young inequality to treat the third term, we obtain that for µ ∈ (0, µ ⋆ ),

1 2 d dt Φ(t) 2 H - µ ε 2 Φ ⊥ (t) 2 H1 -µ ⋆ Φ(t) 2 H1 + C h 1 (t) 2 H h 1 (t) 2 H1 + C ε 2 h 0 (t) E Φ(t) H -µ ⋆ Φ(t) 2 H1 + C h 1 (t) 2 H h 1 (t) 2 H1 + C ε 2 h 0 (t) E Φ(t) H .
In the second term, we decompose h 1 into two parts: h 1 = P 0 h 1 + Φ and use that P 0 = P 2 0 together with the fact that P 0 ∈ B(E, H) to obtain

h 1 (t) 2 H h 1 (t) 2 H1 ∆ 2 0 h 1 (t) 2 H + Φ(t) 2 H1 .
We can thus conclude the proof by taking ∆ 0 small enough and integrating the above differential inequality. Notice that the inequality stated in the Proposition holds for the equivalent "hypocoercive norm" introduced above and thus also holds for the usual norm on H because of the equivalence (uniformly in ε) between those two norms. Combining estimates from Propositions 3.2 and 3.5, one can obtain that Corollary 3.6. Assume that h 0 ∈ E, h 1 ∈ H are such that

sup t 0 h 0 (t) E + h 1 (t) H ∆ 0 < ∞.
For µ ∈ (0, µ ⋆ ) (where µ ⋆ is defined in (2.4)), for ∆ 0 small enough and for any δ > 0, we have that:

(Id -P 0 )h 1 (t) 2 H 1 δ h in 2 E e -µt + (∆ 2 0 + δ + λ ε ) ˆt 0 e -µ(t-s) h 1 (s) 2 H ds. (3.14)
Remark 3.7. The fact that we are able to obtain a multiplicative constant that can be chosen small in front of the second term is very important to recover a decay for h 1 . Indeed, in Proposition 3.3, in the estimate of P 0 h 1 , the term λ ε ˆt 0 e -λε(t-s) (Id -P 0 )h 1 (s) H ds is problematic when applying Gronwall lemma if one hopes to recover some decay in time but the extra small constant that appears in the estimate of (Id -P 0 )h 1 in (3.14) allows us to circumvent this difficulty.

In the end, we are able to prove the following proposition:

Corollary 3.8. Let r ∈ (0, 1). Assume that h 0 ∈ E, h 1 ∈ H are such that sup t 0 h 0 (t) E + h 1 (t) H ∆ 0 < ∞
where ∆ 0 is small enough so that the conclusion of Corollary 3.6 holds. There exists ε 6 ∈ (0, ε 5 ) (where ε 5 is defined in Proposition 3.2) and λ 6 ∈ (0, λ 4 ) (where λ 4 is defined in Lemma 2.7) such that for any ε ∈ (0, ε 6 ) and any λ 0 ∈ [0, λ 6 ) (where λ 0 is defined in Assumption 1.1),

h 1 (t) H C h in H e -(1-r)λεt
where λ ε has been introduced in Proposition 3.3 and the constant C depends on r, ∆ 0 , µ ⋆ (defined in (2.4)) and ν 0 (defined in Lemma 2.3).

Estimates on the kinetic problem.

Combining the previous corollary with Proposition 3.2, we are able to get our final a priori estimates on h in the space E: Proposition 3.9. Let r ∈ (0, 1). Assume that h 0 ∈ E, h 1 ∈ H are such that

sup t 0 h 0 (t) E + h 1 (t) H ∆ 0 < ∞
where ∆ 0 is small enough so that the conclusion of Corollary 3.6 holds. There exists ε † ∈ (0, ε 6 ), λ † ∈ (0, λ 6 ) (where ε 6 and λ 6 are defined in Proposition 3.8) such that for any ε ∈ (0, ε † ) and any λ 0 ∈ [0, λ † ) (where λ 0 is defined in Assumption 1.1), Remark 3.10. Notice that for a fixed ε > 0, the second a priori estimate shows that h = h ε belongs to the space L 1 ([0, ∞), E 1 ). If one is interested in getting bounds on the family {h ε } ε , then we obtain that if λ 0 > 0 (in Assumption 1.1), then the family is bounded in L 1 ([0, ∞), E 1 ) and if λ 0 = 0, then for any T > 0, the family is bounded in L 1 ([0, T ), E 1 ).

h(t) E C h in E e -(
Thanks to the above a priori estimates, we can prove Theorem 1.2 by introducing a suitable iterative scheme that is stable and convergent. We refer to [START_REF] Alonso | Fluid dynamic limit of Boltzmann equation for granular hard-spheres in a nearly elastic regime[END_REF]Section 5] for the details of the proof. We can actually prove the following more precise estimates (which will be useful in what follows) on h 0 ε and h 1 ε that are respectively solutions to (3.1) and (3.2):

h 0 ε L ∞ ([0,∞) ; E) 1 and h 0 ε L 1 ([0,∞) ; E1) ε 2 (3.15)
as well as

h 1 ε L ∞ ([0,∞) ; H) 1 and h 1 ε L 2 ([0,∞) ; H1) 1 (3.16)
where we recall that the spaces H and H 1 are respectively defined in (2.2) and (2.3). Notice that in the previous inequalities, the multiplicative constants only involve quantities related to the initial data of the problem and are independent of ε.

DERIVATION OF THE FLUID LIMIT SYSTEM

The Cauchy theory developed in the previous results give all the a priori estimates that will allow to prove Theorem 1.4. To this end, we make additional assumptions in the definition of the spaces E and E 1 , namely, in this section, those spaces are defined through:

E := W k,1 v W m,2 x ( v q ), E 1 := W k,1 v W m,2 x ( v q+1 ) with m > d, m -1 k 1, q 5. (4.1)
We assume that Assumption 1.1 is met, consider ε, λ 0 and η 0 sufficiently small so that the conclusion of Theorem 1.2 holds in those spaces and consider {h ε } ε a family of solutions to (1.15) constructed in this theorem that splits as h ε = h 0 ε + h 1 ε with h 0 ε and h 1 ε defined in Section 3. We also fix T > 0 for the rest of the section.

4.1. Weak convergence. We start by the following lemma which in particular tells that the microscopic part of h ε vanishes in the limit ε → 0: Lemma 4.1. For any 0 t 1 t 2 T , there holds:

ˆt2 t1 (Id -π 0 )h ε (τ ) E dτ ε √ t 2 -t 1 , (4.2)
where we recall that π 0 is the projection onto the kernel of L 1 defined in (1.23).

Proof. We first remark that

ˆt2 t1 (Id -π 0 )h ε (τ ) E dτ ˆt2 t1 (Id -π 0 )h 0 ε (τ ) 2 E dτ 1/2 √ t 2 -t 1 + ˆt2 t1 (Id -π 0 )h 1 ε (τ ) 2 H1 dτ 1/2 √ t 2 -t 1 .
The first term is estimated thanks to (3.15), which gives:

ˆt2 t1 (Id -π 0 )h 0 ε (τ ) 2 E dτ (Id -π 0 )h 0 ε L ∞ ((0,T ) ; E) (Id -π 0 )h 0 ε L 1 ((0,T ) ; E1) ε 2 .
Concerning the second one, we perform similar computations as in the proof of Proposition 3.5. We recall that h 1 ε solves (3.2) and consider • H an hypocoercive norm on H (see (2.4)). We then have for µ ∈ (0, µ ⋆ ):

1 2 d dt h 1 ε (t) 2 H - µ ε 2 (Id -π 0 )h 1 ε (t) 2 H1 -µ ⋆ h 1 ε (t) 2 H1 + C h 1 ε (t) 2 H h 1 ε (t) 2 H1 + C ε 2 h 0 ε (t) E h 1 ε (t) H
from which we deduce that

1 ε 2 ˆt2 t1 (Id -π 0 )h 1 ε (τ ) 2 H1 dτ h 1 ε (t 1 ) 2 H + ˆt2 t1 h 1 ε (τ ) 2 H h 1 ε (τ ) 2 dτ + ε 2 ˆt2 t1 h 0 ε (τ ) E h 1 ε (τ ) H dτ 1
where we used (3.15) and (3.16) to get the last estimate. Therefore, as for h 0 ε one has

ˆt2 t1 (Id -π 0 )h 1 ε (τ ) 2 H1 dτ ε 2
and this allows to conclude to the wanted estimate.

Using estimates (3.15), (3.16) and (4.2), one can prove the following result of weak convergence (we refer to [START_REF] Alonso | Fluid dynamic limit of Boltzmann equation for granular hard-spheres in a nearly elastic regime[END_REF]Theorem 6.4] for more details on the proof):

Theorem 4.2. Up to extraction of a subsequence, one has

   h 0 ε ε converges to 0 strongly in L 1 ((0, T ) ; E 1 ), h 1 ε ε converges to h weakly in L 2 ((0, T ) ; H) , (4.3)
where h = π 0 (h). In particular, there exist

̺ ∈ L 2 (0, T ); W m,2 x (T d ) , u ∈ L 2 (0, T ); W m,2 x (T d ) d , θ ∈ L 2 (0, T ); W m,2 x (T d ) , such that h(t, x, v) = ̺(t, x) + u(t, x) • v + 1 2 θ(t, x)(|v| 2 -dϑ 1 ) M(v) (4.4)
where M is the Maxwellian distribution introduced in (1.14).

Remark 4.3. Recall that (̺, u, θ) can be expressed in terms of h through the following equalities:

̺(t, x) = ˆRd h(t, x, v) dv, u(t, x) = 1 ϑ 1 ˆRd h(t, x, v)v dv, θ(t, x) = ˆRd h(t, x, v) |v| 2 -dϑ 1 ϑ 2 1 d dv. (4.5)
4.2. Limit system. As mentioned in the introduction, the path that we use to derive the limit system follows the same lines as in the elastic case. The main idea is to write equations satisfied by averages in velocity of h ε and to study the convergence of each term. It is worth mentioning that with the notion of weak convergence at hand presented above, we can adopt an approach which is reminiscent of the program established in [START_REF] Bardos | Fluid dynamic limits of kinetic equations. I. Formal derivations[END_REF][START_REF] Bardos | Fluid dynamic limits of kinetic equations. II. Convergence proofs for the Boltzmann equation[END_REF] but simpler. In particular, we can adapt some of the main ideas of [START_REF] Golse | The Navier-Stokes limit of the Boltzmann equation for bounded collision kernels[END_REF] regarding the delicate convergence of nonlinear terms. The detailed computations and arguments are included in [3, Section 6], we only mention the main steps and keypoints of the proof hereafter. In what follows, we will use the following notation: for g = g(x, v),

g := ˆRd g( • , v) dv.
Local conservation laws. We introduce

A(v) := v ⊗ v - 1 d |v| 2 Id and p ε := 1 d |v| 2 h ε (4.6) so that v ⊗ v h ε = A h ε + p ε Id. We integrate in velocity equation (1.15) multiplied by 1, v i , 1 2 |v| 2 , to obtain ∂ t h ε + 1 ε div x v h ε = 0, (4.7a) ∂ t v h ε + 1 ε Div x A h ε + 1 ε ∇ x p ε = 1 -α(ε) ε 2 v h ε , (4.7b)
∂ t 1 2 |v| 2 h ε + 1 ε div x 1 2 |v| 2 v h ε = 1 ε 3 J α(ε) (f ε , f ε ) + 2(1 -α(ε)) ε 2 1 2 |v| 2 h ε , (4.7c)
where we recall that f ε = G α(ε) + εh ε and where we have introduced

J α (f, f ) := ˆRd [Q α (f, f ) -Q α (G α , G α )] |v| 2 dv.
The goal is to study the convergence of each term in (4.7a)-(4.7b)-(4.7c). A first important remark to address this point is that thanks to the estimates recalled in (3.15)-(3.16), one can prove that for any function ψ = ψ(v) satisfying the bound |ψ(v)| v q , we have the following convergence in the distributional sense:

ψ h ε ---→ ε→0 ψ h in D ′ t,x (4.8)
where h is defined in (4.4) (see [START_REF] Alonso | Fluid dynamic limit of Boltzmann equation for granular hard-spheres in a nearly elastic regime[END_REF]Lemma 6.6]). Roughly speaking, the convergence of the terms in the LHS of (4.7a)-(4.7b)-(4.7c) is treated as in the elastic case. The RHS is going to be handled as a source term which takes into account the drift term and the dissipation of kinetic energy at the microscopic level. In this regard, using (4.8), we first remark that under Assumption 1.1,

1 -α(ε) ε 2 v h ε ---→ ε→0 ϑ 1 λ 0 u in D ′ t,x , (4.9)
since λ 0 = lim ε→0 + ε -2 (1 -α(ε)) and from the definition of u in (4.5). We then present a result of convergence for ε -3 J α(ε) (f ε , f ε ) in the following lemma in the proof of which there are not major difficulties. The said proof is thus omitted, we just mention that it is based on Assumption 1.1, on the estimates on h ε coming from (3.15)-(3.16) and involves the dissipation of energy (1.7), we refer to [3, Lemma 6.9] for more details.

Lemma 4.4. It holds that

1 ε 3 J α(ε) (f ε , f ε ) ---→ ε→0 J 0 in D ′ t,x , where J 0 (t, x) := -λ 0 c ϑ 3 2 1 ̺(t, x) + 3 4 ϑ 1 θ(t, x)
for some positive constant c depending only on the angular kernel b( • ) and d and where λ 0 is defined in Assumption 1.1.

Incompressibility condition and Boussinesq relation. Using (4.8) in the equations (4.7a)-(4.7b), using also that the restitution coefficient satisfies Assumption 1.1, we can easily obtain the incompressibility condition as well as the Boussinesq relation:

div x u = 0 and ∇ x (̺ + ϑ 1 θ) = 0 (4.10)
where we recall that ̺, u and θ are defined in (4.5). Using furthermore that the global mass of h ε vanishes (see (1.20)), we have that Equations of motion and temperature. In order to identify the equations satisfied by u and θ, as in the elastic case, we start by studying the convergence of quantities that are related to

0 = ˆTd ×R d h ε (t, x, v) dv dx ---→ ε→0 ˆTd ̺(t,
̺ ε (t, x) := ˆRd h ε (t, x, v) dv, u ε (t, x) := 1 ϑ 1 ˆRd h ε (t, x, v)v dv, θ ε (t, x) := ˆRd h ε (t, x, v) |v| 2 -dϑ 1 ϑ 2 1 d
dv.

(4.12)

More precisely, we inverstigate the convergence of

u ε := exp -t 1 -α(ε) ε 2 Pu ε and θ ε := 1 2 (|v| 2 -(d + 2)ϑ 1 )h ε
where P is the Leray projection on divergence-free vector fields. Notice that if we compare our approach to the elastic case, we have added the exponential term in the definition of u ε in order to absorbe the term in the RHS in (4.7b). We compute the evolution of u ε and θ ε (by applying the Leray projector P to (4.7b) and by making an appropriate linear combination of (4.7a) and (4.7c)) and obtain:

∂ t u ε = -exp -t 1 -α(ε) ε 2 P ϑ -1 1 Div x 1 ε A h ε (4. 13
)
where A is defined in (4.6) and where we used the incompressibility condition and the strong Boussinesq relation given in (4.10)- (4.11). We refer to [3, Lemma 6.10] for more details.

∂ t θ ε + 1 ε div x b h ε = 1 ε 3 J α(ε) (f ε , f ε ) + 2(1 -α(ε)) ε 2
About initial data. Recall that, in Theorem 4.2, the convergence of h ε to h given by (4.4) is known to hold only for a subsequence and, in particular, at initial time, different subsequences could converge towards different initial datum and therefore (̺, u, θ) could be different solutions to the same system. In Theorem 1.4, the initial datum is prescribed by ensuring the convergence of π 0 h ε in towards a single possible limit where π 0 is defined in (1.23) (recall that the initial data for (̺, u, θ) is defined in (1.31)).

Using Lemma 4.1, one can apply Arzelà-Ascoli theorem to get that Pu ε and θ ε converge in C [0, T ] ; W m-1,2

x towards respectively u and θ 0 defined in (4.15) that also belong to C [0, T ] ; W m-1,2

x . We refer to [START_REF] Alonso | Fluid dynamic limit of Boltzmann equation for granular hard-spheres in a nearly elastic regime[END_REF]Proposition 6.19] for more details.

Limit equations. To get the limit equations, we need to study the convergence of the terms ε -1 PDiv x A h ε and ε -1 div x b h ε in (4.13) and (4.14). To this end, our approach relies on arguments coming from [START_REF] Golse | The Navier-Stokes limit of the Boltzmann equation for bounded collision kernels[END_REF] (in particular, the tricky convergence of the nonlinear terms is treated thanks to a compensated compactness argument coming from [START_REF] Lions | Une approche locale de la limite incompressible[END_REF]), the main difference being that we force the elastic collision operator to appear in our computations, we thus introduce terms that involve differences between the elastic and the inelastic collision operators. Those remainder terms vanish in the limit ε → 0 thanks to Assumption 1.1. We refer to [3, Lemmas 6.12-6.13-6.14] for more details. In the end, writing PDiv x (u ⊗ u) = Div x (u ⊗ u) + ϑ -1 1 ∇ x p (see [START_REF] Majda | Vorticity and incompressible flow[END_REF]Proposition 1.6]), we obtain the following result: Proposition 4.6. There are some constants ν > 0 and γ > 0 such that the limit velocity u = u(t, x) in (4.4) satisfies ∂ t u -ν ϑ 1 ∆ x u + ϑ 1 Div x (u ⊗ u) + ∇ x p = λ 0 u (4.16)

where λ 0 is defined in Assumption 1.1, while the limit temperature θ = θ(t, x) in (4.4) satisfies

∂ t θ - γ ϑ 2 1 ∆ x θ + ϑ 1 u • ∇ x θ = 2 (d + 2)ϑ 2 1 J 0 + 2dλ 0 d + 2 E + 2 d + 2 d dt E, (4.17)
where we recall that J 0 is defined in Lemma 4.4 and E is defined in (4.11).

Remark 4.7. The viscosity and heat conductivity coefficients ν and γ are explicit and fully determined by the elastic linearized collision operator L 1 (see [START_REF] Alonso | Fluid dynamic limit of Boltzmann equation for granular hard-spheres in a nearly elastic regime[END_REF]Lemma C.1]). Notice also that, due to (4.10), Div x (u ⊗ u) = (u • ∇ x ) u and (4.16) is nothing but a reinforced Navier-Stokes equation associated to a divergence-free source term given by λ 0 u which can be interpreted as an energy supply/self-consistent force acting on the hydrodynamical system because of the self-similar rescaling.

To end the identification of the limit equations, we go back to the strong Boussinesq equation (4.11) and prove the following result: where E = E(t) is defined in (4.11). Consequently, the limiting temperature θ(t, x) in (4.4) satisfies Gathering the results we obtained in Propositions 4.6 and 4.8, we are able to end the proof of Theorem 1.4.

∂ t θ - γ ϑ 2 1 ∆ x θ + ϑ 1 u • ∇ x θ = λ 0 c 2(d + 2) ϑ 1 θ. ( 4

  and, taking ψ = 1 and then ψ = v in (1.3) yields the following conservation of macroscopic density and bulk velocity defined asR(t) := ˆTd ℓ ×R d F (t, x, v) dv dx and U (t) := ˆTd ℓ ×R d vF (t, x, v) dv dx,for some solution F (t, x, v) to (1.1):

Assumption 1 . 1 .

 11 Knudsen number. The central underlying assumption in our study is the following relation between the restitution coefficient and the Knudsen number. The restitution coefficient α( • ) is a continuously decreasing function of the Knudsen number ε satisfying the scaling behaviour α

Corollary 2 . 4 .

 24 but do not change the notations G 1,ε and B 1,ε . The following corollary states immediate consequences of the previous lemma (we denote by R( • , B ε ) the resolvent of the operator B ε ): There holds:

. 21 )

 21 Thanks to Assumption 1.1, one can find ε 4 and λ 4 depending only on χ such that ℓ(ε) < 1 for any ε ∈ (0, ε 4 ) and λ 0 ∈ [0, λ 4 ). In particular, we deduce (2.20) from[START_REF] Kato | Perturbation theory for linear operators[END_REF] Paragraph I.4.6].

1 2 |v| 2 h ε with b(v) := 1 2 |v| 2 -

 12 (d + 2)ϑ 1 . (4.14)The study of the limit ε → 0 in those equations is more favorable because compared to (4.7a)-(4.7b)-(4.7c), the gradient term in (4.7b) has been eliminated thanks to the Leray projector and also because A and b belong to the range of Idπ 0 so that thanks to Lemma 4.1, we know that the quantities ε-1 Div x A h ε and ε -1 div x b h ε are bounded in W m-1,2 x . Then, applying a precised version of Aubin-Lions lemma [20, Corollary 4], we are able to prove that up to the extraction of a subsequence, {u ε } ε and {θ ε } ε converge strongly in L 1 (0, T ) ; W m-1,2 x respectively towards Pu = u and θ 0 := 1 2 (|v| 2 -(d + 2)ϑ 1 )h =

Proposition 4 . 8 .

 48 It holds thatE(t) = 0, t ∈ [0, T ],

. 18 ) 2 x. 1 2

 1821 where γ is defined in Proposition 4.6, λ 0 in Assumption 1.1 and c in Lemma 4.[START_REF] Alonso | Uniqueness and regularity of steady states of the Boltzmann equation for viscoelastic hard-spheres driven by a thermal bath[END_REF]. Moreover, the strong Boussinesq relation holds true:̺ + ϑ 1 θ = 0 on [0, T ] × T d . (4.19)Proof. Using Lemma 4.4 and averaging in position the equation (4.17), it is easy to prove that d dtE(t) = c0 E(t)for some some constant c0 ∈ R. Moreover, on the one hand, from (1.31), we haveE(0) = ˆTd θ(0, x) dx = -1 ϑ 1 ˆTd ̺(0, x) dx. (4.20)On the other hand, from the definition of θ 0 in (4.15), we also have that θ ε converges towards θ 0 in C [0, T ] ; W m-1,Consequently, we deduce thatˆTd θ 0 (0, x) dx = lim ˆTd |v| 2 -(d+2)ϑ1 2 h ε (0, x) dx = lim ε→0 ˆTd |v| 2 h ε (0, x) dxwhere we used (1.20) to get the last equality. From (1.22), we deduce that ˆTd θ 0 (0, x) dx = 0.Coming back to (4.20)-(4.21), we deduce thatE(0) = -1 ϑ 1 ˆTd ̺(0, x) dx = 2 ϑ 1 d ˆTd ̺(0, x) dxwhich implies that E(0) = 0. This concludes the proof.

  -α(ε))/ε 2 has been defined in Theorem 2.1 and the constant C depends on r, ∆ 0 , µ ⋆ (defined in (2.1)) and ν 0 (defined in Lemma 2.3).

	1-r)λεt	and	ˆt 0	h(s) E1 ds C h in	E min 1 + t, 1 +	1 λ ε
	where λ ε ∼

ε→0

(1

  x) dx in D ′ t and thus that ´Td ̺(t, x) dx = 0. It implies that we have the following strengthened Boussinesq relation: for almost every (t, x) ∈ (0, T ) × T d , Notice here that the derivation of the strong Boussinesq relation ̺ + ϑ 1 θ = 0 is not as straightforward as in the elastic case. In the elastic case, the classical Boussinesq relation ∇ x (̺+ϑ 1 θ) = 0 straightforwardly implies the strong form of Boussinesq because the two functions ̺ and θ have zero spatial averages. This cannot be deduced directly in the granular context due to the dissipation of energy and we will see later on how to obtain it (seeProposition 4.8).

	ˆTd
	θ(t, x) dx.	(4.11)
	Remark 4.5.

̺ + ϑ 1 (θ -E) = 0 with E = E(t) :=