Energy equalities for compressible Navier-Stokes equations

0 downloads 0 Views 302KB Size Report
Oct 24, 2018 - Key words and phrases. ... Roughly speaking, Onsager's conjecture (for Euler equations) addresses the question how much ...... Department of Mathematics, New York University Abu Dhabi, Abu Dhabi, United Arab Emirates.
ENERGY EQUALITIES FOR COMPRESSIBLE NAVIER-STOKES EQUATIONS

arXiv:1810.10264v1 [math.AP] 24 Oct 2018

QUOC-HUNG NGUYEN, PHUOC-TAI NGUYEN, AND BAO QUOC TANG Abstract. The energy equalities of compressible and inhomogeneous incompressible NavierStokes equations are studied. By using a unified approach, we give sufficient conditions on the regularity of weak solutions for these equalities to hold. The method of proof is suitable for the case of periodic as well as homogeneous Dirichlet boundary conditions. In particular, by a careful analysis using homogeneous Dirichlet boundary condition, no boundary layer assumptions are required when dealing with bounded domains with boundary.

Contents 1. Introduction and Main Results 2. Preliminaries 3. Proof of Theorem 1.2 3.1. Estimate of (A) 3.2. Estimate of (B) 3.3. Estimate of (C) 3.4. Estimate of (D). 3.5. Conclusion of the Proof of Theorem 1.2 4. Proof of Theorem 1.5 4.1. Estimate of (E) 4.2. Estimate of (F ) 4.3. Estimate of (G) 4.4. Estimate of (H) 4.5. Conclusion of the Proof of Theorem 1.5 5. Proof of Theorem 1.7 6. Proof of Theorem 1.8 7. Proof of Theorem 1.9 References

1 7 11 11 12 13 15 15 15 16 16 17 19 19 19 20 21 22

1. Introduction and Main Results Let d = 2, 3 and Ω be either the torus Td or a bounded domain in Rd with C 2 boundary ∂Ω. This paper studies the energy equalities for compressible Navier-Stokes equation with degenerate viscosity ( ∂t ̺ + ∇ · (̺u) = 0, in Ω × (0, T ), (cNSd) γ ∂t (̺u) + ∇ · (̺u ⊗ u) + ∇(̺ ) − 2ν∇ · (̺Du) = 0, in Ω × (0, T ), Date: October 25, 2018. 2010 Mathematics Subject Classification. 35Q30, 76B03, 76D05, 76N10. Key words and phrases. Compressible Navier-Stokes equations; Inhomogeneous incompressible Navier-Stokes equations; Energy equalities; Onsager’s conjecture. 1

2

Q-H. NGUYEN, P-T. NGUYEN, AND B. Q. TANG

and the equation ( ∂t ̺ + ∇ · (̺u) = 0, ∂t (̺u) + ∇ · (̺u ⊗ u) + ∇(̺γ ) − 2ν∆u − λ∇(∇ · u) = 0,

in Ω × (0, T ), in Ω × (0, T ),

as well as the inhomogeneous incompressible Navier-Stokes equation   in Ω × (0, T ), ∂t ̺ + ∇ · (̺u) = 0, ∂t (̺u) + ∇ · (̺u ⊗ u) + ∇P − 2ν∇ · (̺Du) = 0, in Ω × (0, T ),  ∇ · u = 0, in Ω × (0, T ),

with initial data



(̺u)(x, 0) = ̺0 (x)u0 (x),

̺(x, 0) = ̺0 (x), and homogeneous Dirichlet boundary condition

x ∈ Ω,

x ∈ Ω,

u = 0 on ∂Ω × (0, T ).

(cNS)

(icNS)

(1)

(2)

Here T > 0, γ > 1, ν, λ > 0, ̺ denotes the density of the fluid, u is the velocity of fluid and D = 21 [∇u + ∇⊤ u] stands for the strain tensor, and in the case of (icNS), P is the scalar pressure. Naturally, the boundary condition u = 0 on ∂Ω × (0, T ) is only imposed in the case of bounded domain. Energy conservation is one aspect of the Onsager’s conjecture which was announced in his celebrated paper on statistical hydrodynamics [Ons49]. More precisely, Onsager [Ons49] conjectured that, in the context of homogeneous incompressible Euler equations, kinetic energy is globally conserved for H¨older continuous solutions with the exponent greater than 1/3, while energy dissipation phenomenon occurs for H¨older continuous solutions with the exponent less than 1/3. The ‘positive’ part of the conjecture was first proved by by Eyink [Eyi94] and Constantin-E-Titi [CET94] for the whole space Rd or periodic boundary conditions, i.e. Td . Significant contributions in the case of domains with boundary were recently achieved in [BT18, BTW, DN18, NN18] for homogeneous incompressible Euler equations and in [NNT18] for inhomogeneous incompressible and compressible isentropic Euler equations. The other direction of the Onsager’s conjecture was initiated in the groundbreaking paper of Scheffer [Sch93], then has reached to its full flowering with a work of Shnirelman [Shn97], a series of celebrated works of De Lellis and Sz´ekelyhidi [DS12, DS13, DS14, BDIS15], and recently was settled by Isett in [Ise18a, Ise18] and by Buckmaster et al. in [BDSV18]. Roughly speaking, Onsager’s conjecture (for Euler equations) addresses the question how much regularities needed for a weak solution to conserve energy. In the context of classical incompressible Navier-Stokes equations, since global regularity in three dimensions has been a long standing open problem, it is natural to ask how much regularities needed for a weak solution to satisfy the energy equality (rather than energy conservation in the context of Euler equations). This topic of research has been studied already in the sixties starting with the work of Serrin [Ser62] where he asserted that the energy equality must hold if u ∈ Lp (0, T ; Lq (Td )) with 2p + dq ≤ 1. Later Shinbrot [Shi74] assumed a condition which is independent of dimensions, more precisely u ∈ Lp (0, T ; Lq (Td )) with 1 + 1q ≤ 21 and q ≥ 4. By using a new approach based on a lemma introduced by Lions, C. Yu p [Yu16] obtained the same result as in [Shi74]. All these results deal with either Ω = Rd or Ω = Td . Recently, C. Yu [Yu17a] considered the case where Ω is a bounded domain with C 2 boundary and proved the energy equality under additional condition u ∈ Ls (0, T ; Bsα,∞ (Ω)) (with s > 2 and 1 + 1s < α < 1) which enables to deal with the boundary effects. Here Bsα,∞ (Ω) denotes the Besov 2 space.

ENERGY EQUALITIES FOR COMPRESSIBLE NAVIER-STOKES EQUATIONS

3

Much less is known concerning energy equalities for compressible or inhomogeneous incompressible Navier-Stokes equations. Recent results for (cNSd) and (cNS) in Td were carried out by C. Yu [Yu17]. More precisely, he gave sufficient conditions on the regularity of the density ̺ and the velocity u for the validity of the energy equality. The framework employed in [Yu17] is remarkably different from that in the incompressible case (see [Yu16, Yu17a]). In particular, solutions defined in [Yu17, Definition 1] are required to satisfy a set of regularities which allow to deduce √ the continuity of ( ̺u)(t) in the strong topology at t = 0. The existence of this kind of solutions is guaranteed by [VY16] and [LV18]. Motivated by the aforementioned works, we present in this paper a unified approach to show energy equalities for the compressible and inhomogeneous incompressible Navier-Stokes equations (cNSd), (cNS) and (icNS). The main idea is inspired by our recent work for compressible Euler equations [NNT18], which is different from the method employed in e.g. [Yu16, Yu17, Yu17a]. In particular, we use the test function (̺ε )−1 (̺u)ε instead of uε , where the convolution is only taken in spatial variables. The choice of this test function allows to obtain mild conditions on the density ̺ and the velocity u which are in fact weaker than previous works1. For instance, when d = 3 5 we assume only u ∈ L4 (Td × (0, T )) while [Yu17] required u ∈ Lp (0, T ; Lq (Td )) where 1p + 1q ≤ 12 and q ≥ 6 (see more discussion after Theorem 1.2). Moreover, by carefully using the Dirichlet boundary conditions (in the case of bounded domains), we show that no additional regularities near the boundary on the velocity u are required. Before stating the main results, we introduce the definition of weak solutions. Definition 1.1. A couple (̺, u) is called a weak solution to (cNSd) with initial data (1) if (i) ˆ Tˆ (̺∂t ϕ + ̺u · ∇ϕ)dxdt = 0 0

(ii)

for every test function ϕ ∈ ˆ

0

T

ˆ



Ω ∞ C0 (Ω

× (0, T )).

(̺u · ∂t ψ + ̺u ⊗ u : ∇ψ + ̺γ ∇ · ψ − 2ν̺Du : ∇ψ)dxdt = 0

for every test vector field ψ ∈ C0∞ (Ω × (0, T ))d. (iii) ̺(·, t) ⇀ ̺0 in D ′ (Ω) as t → 0, i.e. ˆ ˆ lim ̺(x, t)ϕ(x)dx = ̺0 (x)ϕ(x)dx t→0





(4)

(5)



for every test function ϕ ∈ C0∞ (Ω). (iv) (̺u)(·, t) ⇀ ̺0 u0 in D ′ (Ω) as t → 0, i.e. ˆ ˆ lim (̺u)(x, t)ψ(x)dx = (̺0 u0 )(x)ψ(x)dx t→0

(3)

(6)



for every test vector field ψ ∈ C0∞ (Ω)d . Weak solutions to (cNS) and (icNS) can be defined similarly. Remark 1.1. The notion of weak solutions defined in Definition 1.1 requires only modest regularities of the density and the velocity, e.g. for (3) and (4) to make sense, one only needs ̺, ̺u, ̺u ⊗ u, ̺γ , ̺Du ∈ L1loc (Ω × (0, T )). If one were to show energy equality for weak solutions 1Admittedly,

our approach seems not directly applicable to the case with vacuum.

4

Q-H. NGUYEN, P-T. NGUYEN, AND B. Q. TANG

defined in [VY16] (or [LV18]), as it was done in [Yu17], more regularities on these solutions are already given, see [VY16, Theorem 1.2] or [Yu17, Definition 1.1]. Theorem 1.2. Let Ω = Td and (̺, u) be a weak solution of (cNSd) with initial data (1). Assume that 0 < c1 ≤ ̺(t, x) ≤ c2 < ∞ and u ∈ L∞ (0, T ; L2 (Td )) ∩ L2 (0, T ; H 1(Td )). Define α := min{γ; 2}.

(7)

sup sup |h|− α k̺(· + h, t) − ̺(·, t)kLα (T2 ) < ∞,

(8)

We assume in the case d = 2 that 1

t∈(0,T ) |h| δ}. Since Ω is a bounded, connected domain with C 2 boundary, we find r0 > 0 and a unique Cb1 -vector function n : Ω\Ωr0 → S d−1 such that the following holds true: for any r ∈ [0, r0 ), x ∈ Ωr \Ωr0 there exists a unique xr ∈ ∂Ωr such that d(x, ∂Ωr ) = |x − xr | and n(x) is the outward unit normal vector field to the boundary ∂Ωr at xr . The energy equality for (cNSd) in a domain with boundary is stated in the following theorem. Theorem 1.5. Let Ω be a bounded domain with C 2 boundary ∂Ω and (̺, u) be a weak solution of (cNSd) with initial data (1) and Dirichlet boundary condition (2). Assume that 0 < c1 ≤ ̺(x, t) ≤ c2 < ∞,

and

u ∈ L∞ (0, T ; L2 (Ω)) ∩ L2 (0, T ; H 1(Ω)).

We assume in the case d = 2 that, for each 0 < δ < 1, 1

sup sup |h|− α k̺(· + h, t) − ̺(·, t)kLα (Ωδ ) < ∞,

(12)

t∈(0,T ) |h|≤δ

and in case d = 3 that, for each 0 < δ < 1, u ∈ L4 (Ω × (0, T ))

and

1

sup sup |h|− α k̺(· + h, t) − ̺(·, t)kL 65 α (Ω ) < ∞.

t∈(0,T ) |h|≤δ

δ

Then the energy equality holds, i.e.  ˆ tˆ ˆ  ̺(x, t)γ 1 2 dx + ν (̺|u| )(x, t) + ̺|Du|2 dxdt 2 γ − 1 0 Ω Ω  ˆ  1 ̺γ0 2 = dx ∀t ∈ (0, T ). ̺0 |u0 | + 2 γ−1 Ω

(13)

(14)

Remark 1.6. In contrast to Euler equations, we do not need any additional regularities near the boundary to establish energy equality for Navier-Stokes equations in bounded domains. This is due to a careful use of homogeneous Dirichlet boundary conditions and the assumption u ∈ L2 (0, T ; H 1(Ω)) (see Lemma 2.6). As mentioned above, our method of proof is also suitable to obtain energy equalities for the compressible Navier-Stokes equation (cNS) and the inhomogeneous incompressible Navier-Stokes equation (icNS). Theorem 1.7. Let either Ω = Td or Ω ⊂ Rd be a bounded domain with C 2 boundary ∂Ω where d = 2, 3. Let (̺, u) be a weak solution to (cNS) with initial data (1) (and with Dirichlet boundary condition (2) in case Ω is a bounded domain). Assume 0 ≤ c1 ≤ ̺(x, t) ≤ c2 < +∞

and

u ∈ L∞ (0, T ; L2 (Ω)) ∩ L2 (0, T ; H 1(Ω)).

Moreover, if Ω = Td we assume (8) for d = 2 and (9) for d = 3, and if Ω is a bounded domain we assume (12) for d = 2 and (13) for d = 3. Then the energy equality holds, i.e.  ˆ tˆ ˆ tˆ ˆ  ̺(x, t)γ 1 2 2 dx + 2ν (̺|u| )(x, t) + ̺|∇u| dxdt + λ |∇ · u|2 dxds 2 γ − 1 0 Ω 0 Ω Ω ˆ  γ  ̺ 1 dx ∀t ∈ (0, T ). ̺0 |u0 |2 + 0 = 2 γ −1 Ω

6

Q-H. NGUYEN, P-T. NGUYEN, AND B. Q. TANG

Theorem 1.8. Let either Ω = Td or Ω ⊂ Rd be a bounded domain with C 2 boundary ∂Ω where d = 2, 3. Let (̺, u, P ) be a weak solution to (icNS) with initial data (1) (and with Dirichlet boundary condition (2) in case Ω is a bounded domain). Assume 0 ≤ c1 ≤ ̺(x, t) ≤ c2 < +∞,

u ∈ L∞ (0, T ; L2(Ω)) ∩L2 (0, T ; H 1(Ω)),

and P ∈ L2 (Ω×(0, T )).

Moreover, in the case d = 3 we assume additionally u ∈ L4 (Ω × (0, T )). Then the energy equality holds, i.e. ˆ ˆ tˆ ˆ 1 1 2 2 ̺(x, t)|u(x, t)| dx + ν ̺0 (x)|u0 (x)|2 dx ∀ t ∈ (0, T ). ̺|Du| dxds = 2 Ω 2 Ω 0 Ω We emphasize again that, in Theorem 1.8, the condition u ∈ L2 (0, T ; H 1(Ω)) helps to handle the boundary effects without requiring extra conditions of u near the boundary. As a consequence, we improve the results of [Yu17a] for homogeneous incompressible Navier-Stokes equations by removing the assumption u ∈ Ls (0, T ; Bsα,∞(Ω)) with 12 + 1s < α < 1. Theorem 1.9. Let Ω ⊂ R3 be a bounded domain with C 2 boundary and (u, P ) be a weak solution to the homogeneous incompressible Navier-Stokes equation  ∂t u − µ∆u + ∇ · (u ⊗ u) + ∇P = 0, in Ω × (0, T ),    ∇ · u = 0, in Ω × (0, T ), (15)  u = 0, on ∂Ω × (0, T ),    u(x, 0) = u0 (x), in Ω where µ > 0 is the viscosity. Assume u ∈ L4 (Ω × (0, T )) and there exists δ0 > 0 such that P ∈ L2 (Ω\Ωδ0 × (0, T )).

(16)

Then the energy equality holds ˆ ˆ tˆ ˆ 2 2 |u(x, t)| dx + µ |∇u(x, s)| dxds = |u0 (x)|2 dx ∀ t ∈ (0, T ). Ω

0





Remark 1.10. The assumption (16) is to deal with the boundary layer. We remark that [Yu17a] did not impose any condition on the pressure, but the author used P = 0 on the boundary (see [Yu17a, Proposition 2.3]) which is neither assumed nor implied from (15). Nevertheless, if someP how P = 0 on the boundary (in a weak sense), from the equation −∆P = di,j=1 ∂xi ∂xj (ui uj ) we obtain kP kL2 (Ω×(0,T )) ≤ Ckuk2L4 (Ω×(0,T )) and therefore (16) is automatically satisfied. The organization of this paper is as follows: In the next section, we prove some auxiliary estimates which will play important roles in our proof. The proofs of Theorems 1.2, 1.5, 1.7, 1.8 and 1.9 are presented in the last five sections respectively. Notation. Throughout the paper, C denotes generic constants which may depend on d, T , k̺kL∞ (Ω×(0,T )) , k̺−1 kL∞ (Ω×(0,T )) and other scalar parameters. We use the notation kf (s)kLp (Ω) to denote kf (·, s)kLp (Ω) . ffl ´ For any Borel set E, we denote by E f (x)dx = Ld1(E) E f (x)dx the average of f over E, where Ld (E) is the Lebesgue measure of E.

ENERGY EQUALITIES FOR COMPRESSIBLE NAVIER-STOKES EQUATIONS

7

2. Preliminaries −

1

1−|x|2 for |x| < 1 and ω(x) = 0 for Let ω : Rd → R be a standard mollifier, ´ i.e. ω(x) = c0 e |x| ≥ 1, where c0 is a constant such that Rd ω(x)dx = 1. For any ε > 0, we define the rescaled mollifier ωε (x) = ε−d ω( xε ). For any function f ∈ L1loc (Ω), its mollified version is defined as ˆ ε f (x) = (f ⋆ ωε )(x) = f (x − y)ωε(y)dy, x ∈ Ωε ,

Rd

recalling Ωε = {x ∈ Ω : d(x, ∂Ω) > ε}.

Lemma 2.1. Let 2 ≤ d ∈ N, 1 ≤ p, q ≤ ∞ and f : Td × (0, T ) → R. (i) Assume f ∈ Lp (0, T ; Lq (Td )). Then for any ε > 0, there holds d

kf ε kLp (0,T ;L∞ (Td )) ≤ Cε− q kf kLp (0,T ;Lq (Td )) ,

(17)

k∇f ε kLp (0,T ;L∞ (Td )) ≤ Cε−1− q kf kLp (0,T ;Lq (Td )) .

(18)

k∇f ε kLp (0,T ;Lq (Td )) ≤ Cε−1 kf kLp (0,T ;Lq (Td )) .

(19)

lim sup εk∇f εkLp (0,T ;Lq (Td )) = 0.

(20)

d

(ii) Assume f ∈ Lp (0, T ; Lq (Td )). Then for any ε > 0, there holds Moreover, if p, q < ∞ then

ε→0

(iii) Assume f ∈ Lp (0, T ; Lq (Td )) and g : Td × (0, T ) → R with 0 < c1 ≤ g ≤ c2 < ∞. Then for any ε > 0, there holds



∇ ≤ C(c1 , c2 )ε−1 kf kLp (0,T ;Lq (Ω)) . (21)

gε p L (0,T ;Lq (Td )) Moreover, if p, q < ∞ then





lim sup ε ∇ ε = 0. g Lp (0,T ;Lq (Td )) ε→0

(22)

(iv) Assume f ∈ L2 (0, T ; H 1(T2 )). Then for any ε > 0, there holds

k∇f ε kL2 (0,T ;L∞ (T2 )) ≤ Cε−1 kf kL2 (0,T ;H 1 (T2 )) .

(23)

lim sup εk∇f ε kLr (0,T ;L∞ (T2 )) = 0.

(24)

Moreover, for any r ∈ [1, 2], there holds ε→0

Proof. (i) By the definition of f ε and H¨older’s inequality, for a.e. x ∈ Td and s ∈ (0, T ), we have  q−1 ˆ  1q ˆ ˆ q q |f ε (x, s)| ≤ |f (x − y, s)ωε(y, s)|dy ≤ |f (x − y, s)|q dy |ωε (y)| q−1 dy Td

Td

≤ Cε

− dq

Td

kf (s)kLq (Td ) ,

where we have used ωε (x) = ε−d ω(x/ε) at the last step. This implies d

kf ε (s)kL∞ (Td ) ≤ Cε− q kf (s)kLq (Td ) ,

which in turn yields (17). Next we use H¨older’s inequality again to estimate ˆ d ε |∇f (x, s)| ≤ |f (x − y, s)||∇ωε(y)|dy ≤ Cε−1− q kf (s)kLq (Td ) . Td

8

Q-H. NGUYEN, P-T. NGUYEN, AND B. Q. TANG

This leads to (18). (ii) From the fact that and s ∈ (0, T ),

´

Rd

∇ωε (y)dy = 0 and H¨older’s inequality, we obtain, for a.e. x ∈ Td

ˆ |∇f (x, s)| ≤

[f (x − y, s) − f (x, s)]∇ωε (y)dy |y| 0 to be small, we integrate (79) with respect to ε2 on (ε1 , ε1 + ε3 ) to get ˆ ˆ ˆ ˆ ˆ ˆ (̺u)ε (̺u)ε 1 ε1 +ε3 t 1 ε1 +ε3 t ε ∂ (̺u) dxdsdε + ∇ · (̺u ⊗ u)ε dxdsdε2 t 2 ε ε ε3 ε1 ̺ ε ̺ 3 ε1 τ Ωε2 τ Ωε2 ˆ ε1 +ε3 ˆ t ˆ ˆ ε1 +ε3 ˆ t ˆ ε ε 1 2ν (̺u) ε (̺u) γ ε ∇ · (̺Du) + ∇(̺ ) dxdsdε − dxdsdε2 = 0. 2 ε3 ε1 ̺ε ε3 ε1 ̺ε τ Ωε2 τ Ωε2 (80) We denote by (E), (F ), (G) and (H) the first, second, third and forth term on the left hand side of (80) respectively. We will estimate them separately in the following subsections. 4.1. Estimate of (E). This term is estimated similarly to (A),   ˆ ˆ ˆ ˆ ˆ ˆ ε 2 1 1 ε1 +ε3 t |(̺u)ε |2 1 1 ε1 +ε3 t ε |(̺u) | (E) = ∂t ∂ ̺ dxdsdε + dxdsdε2 t 2 2 ε3 ε1 ̺ε 2 ε3 ε1 (̺ε )2 τ Ωε2 τ Ωε2   ˆ ˆ ˆ ˆ ˆ ˆ ε 2 |(̺u)ε |2 1 1 ε1 +ε3 t 1 1 ε1 +ε3 t ε |(̺u) | ∂t ∇ · (̺u) = dxdsdε − dxdsdε2 2 2 ε3 ε1 ̺ε 2 ε3 ε1 (̺ε )2 τ Ωε2 τ Ωε2 =: (E1) + (E2).

The term (E1) is desired, while (E2) will be canceled by (F 4) later. 4.2. Estimate of (F ). We estimate (F ) similarly to (B) by using integration by parts ˆ ˆ ˆ (̺u)ε 1 ε1 +ε3 t ∇ · [(̺u ⊗ u)ε − (̺u)ε ⊗ uε ]dxdsdε2 (F ) = ε ε3 ε1 ̺ τ Ωε2 ˆ ε1 +ε3 ˆ t ˆ (̺u)ε 1 ∇ · ((̺u)ε ⊗ uε )dxdsdε2 + ε ε3 ε1 ̺ τ Ωε2 ˆ ε1 +ε3 ˆ t ˆ 1 |(̺u)ε |2 =: (F 1) + ∇ · uε dxdsdε2 ε3 ε1 ̺ε τ Ωε2 ˆ ˆ ˆ 1 1 ε1 +ε3 t uε + ∇|(̺u)ε |2 dxdsdε2 ε 2 ε3 ε1 ̺ τ Ωε2 ˆ ε1 +ε3 ˆ t ˆ 11 uε =: (F 1) + |(̺u)ε |2 ε n(θ)dHd−1 (θ)dsdε2 2 ε3 ε1 ̺ τ ∂Ωε2 ˆ ε1 +ε3 ˆ t ˆ 11 |(̺u)ε |2 + ∇ · [̺ε uε − (̺u)ε ] dxdsdε2 2 ε3 ε1 (̺ε )2 τ Ωε2 ˆ ˆ ˆ |(̺u)ε |2 1 1 ε1 +ε3 t ∇ · (̺u)ε + dxdsdε2 2 ε3 ε1 (̺ε )2 τ Ωε2 =: (F 1) + (F 2)bdr + (F 3) + (F 4). The superscript “bdr” in (F 2)bdr means that this term contains a boundary layer. It’s obvious that (F 4) + (E2) = 0. The term (F 1) is estimated using integration by parts as ˆ ˆ ˆ (̺u)ε 1 ε1 +ε3 t ∇· [(̺u ⊗ u)ε − (̺u)ε ⊗ uε ]dxdsdε2 (F 1) = − ε3 ε1 ̺ε τ Ωε2 ˆ ˆ ˆ 1 ε1 +ε3 t (̺u)ε + [(̺u ⊗ u)ε − (̺u)ε ⊗ uε ] ε n(θ)dHd−1 (θ)dsdε2 ε3 ε1 ̺ τ ∂Ωε2 =: (F 11) + (F 12)bdr .

ENERGY EQUALITIES FOR COMPRESSIBLE NAVIER-STOKES EQUATIONS

17

We estimate (F 11) similarly to (B1) in (52)–(54) and therefore lim sup lim sup |(F 11)| = 0. ε→0

(81)

τ →0

The term (F 12)bdr will be treated later on, together with other boundary terms. For (F 3) it follows from integration by parts that ˆ ˆ ˆ |(̺u)ε |2 1 1 ε1 +ε3 t [̺ε uε − (̺u)ε ]∇ (F 3) = − dxdsdε2 2 ε3 ε1 (̺ε )2 τ Ωε2 ˆ ˆ ˆ 1 1 ε1 +ε3 t |(̺u)ε |2 + [̺ε uε − (̺u)ε ]n(θ) dHd−1 (θ)dsdε2 ε )2 2 ε3 ε1 (̺ τ ∂Ωε2 =: (F 31) + (F 32)bdr . The term (F 31) is estimated similarly to (B2) in (56)–(58) and therefore lim sup lim sup |(F 31)| = 0. ε→0

(82)

τ →0

It remains to estimate the boundary terms (F 12)bdr , (F 2)bdr and (F 32)bdr . As for the term (F 12)bdr , we use the coarea formula and then letting successively τ → 0 and ε → 0 to obtain lim sup lim sup |(F 12)bdr| = 0. ε→0

(83)

τ →0

Next, we deal with the term (F 2)bdr . By using the coarea formula, Holder’s inequality and Lemma 2.6, we get 1 ˆ tˆ ε ε 2u bdr |(̺u) | ε n(x)dxds lim sup lim sup |(F 2) | = lim sup lim sup ̺ ε1 ,ε→0 τ →0 ε1 ,ε→0 τ →0 2ε3 τ Ωε1 \Ωε1 +ε3 ˆ Tˆ C ≤ |u|3 dxds 2ε3 0 Ω\Ωε3 C kuk2L4 ((Ω\Ωε3 )×(0,T )) kukL2 ((Ω\Ωε3 )×(0,T )) 2ε3 ≤ Ckuk2L4 ((Ω\Ωε3 )×(0,T )) k∇ukL2 ((Ω\Ω2ε3 ×(0,T )) .



Since u ∈ L4 (Ω × (0, T )) and u ∈ L2 (0, T ; H 1(Ω)), by letting ε3 → 0, we obtain lim sup lim sup lim sup |(F 2)bdr | = 0. ε3 →0

ε1 ,ε→0

τ →0

Next we treat the term (F 32)bdr . Again, by employing the coarea formula and letting successively τ → 0 and ε → 0, we derive lim sup lim sup |(F 32)bdr | = 0. ε→0

τ →0

4.3. Estimate of (G). By integration by parts, we have ˆ ˆ ˆ ˆ ˆ ˆ (̺u)ε (̺u)ε 1 ε1 +ε3 t 1 ε1 +ε3 t γ ε ε γ ∇[(̺ ) − (̺ ) ]dxdsdε + ∇(̺ε )γ dxdsdε2 (G) = 2 ε ε ε3 ε1 ̺ ε3 ε1 ̺ τ Ωε2 τ Ωε2 ˆ ε1 +ε3 ˆ t ˆ (̺u)ε γ ε 1 ∇· [(̺ ) − (̺ε )γ ]dxdsdε2 =− ε ε3 ε1 ̺ τ Ωε2 ˆ ε1 +ε3 ˆ t ˆ 1 (̺u)ε + n(θ)[(̺γ )ε − (̺ε )γ ]dHd−1 (θ)dsdε2 ε ε3 ε1 ̺ τ ∂Ωε2

18

Q-H. NGUYEN, P-T. NGUYEN, AND B. Q. TANG

γ 1 − γ − 1 ε3

ˆ

ε1 +ε3 ε1

γ 1 γ − 1 ε3

ˆ

γ 1 (G3) = − γ − 1 ε3

ˆ

+

∇ · (̺u)ε (̺ε )γ−1 dxdsdε2

ˆ tˆ

(̺u)ε n(θ)(̺ε )γ−1 dHd−1 (θ)dsdε2

τ

ε1 +ε3

ε1

=: (G1) + (G2)

ˆ tˆ

Ωε2

τ

bdr

∂Ωε2

+ (G3) + (G4)bdr

where ε1 +ε3 ε1

ˆ tˆ τ

ε

Ωε2

ε γ−1

∇·(̺u) (̺ )

1 1 dxdsdε2 = γ − 1 ε3

ˆ

ε1 +ε3 ε1

ˆ tˆ τ

∂t (̺ε )γ dxdsdε2

Ωε2

is a desired term. The term (G1) is rewritten as ˆ ˆ ˆ 1 ε1 +ε3 t (̺u)ε − ̺ε uε (G1) = ∇[(̺γ )ε − (̺ε )γ ]dxdsdε2 ε ε3 ε1 ̺ τ Ωε2 ˆ ε1 +ε3 ˆ t ˆ 1 (∇ · uε )[(̺γ )ε − (̺ε )γ ]dxdsdε2 + ε3 ε1 τ Ωε2 ˆ ε1 +ε3 ˆ t ˆ 1 + [(̺γ )ε − (̺ε )γ ]uε n(θ)dHd−1 (θ)dsdε2 ε3 ε1 τ ∂Ωε2 =: (G12) + (G13) + (G14)bdr . We can handle (G12) and (G13) using similar arguments to estimate of (C1) in (60) and (68), and thus lim sup lim sup |(G12)| = 0 ε→0

and

τ →0

lim sup lim sup |(G13)| = 0. ε→0

τ →0

It remains to estimate the terms involving the boundary (G14)bdr , (G2)bdr and (G4)bdr . By using the coarea formula and letting successively τ → 0 and ε → 0, we assert that lim sup lim sup |(G14)bdr | = 0 and ε→0

τ →0

lim sup lim sup |(G2)bdr | = 0. ε→0

τ →0

To deal with the term (G4)bdr , we use the coarea formula (44), the fact Ld (Ω\Ωε3 ) ≈ ε3 , Holder’s inequality and Lemma 2.6 to get ˆ ˆ t 1 γ (̺u)ε n(x)(̺ε )γ−1 dxds lim sup lim sup |(G4)bdr | = lim sup lim sup γ − 1 ε ε1 ,ε→0 τ →0 ε1 ,ε→0 τ →0 3 τ Ωε1 \Ωε1 +ε3 ˆ ˆ C T ≤ |u|dxds ε3 0 Ω\Ωε3 !1/2 ˆ Tˆ C |u|2 dxds ≤ (T Ld (Ω\Ωε3 ))1/2 ε3 0 Ω\Ωε3 1

≤ Cε32 k∇ukL2 ((Ω\Ω2ε3 )×(0,T )) . Since u ∈ L2 (0, T ; H 1(Ω)), it follows that lim sup lim sup lim sup |(G4)bdr | = 0. ε3 →0

ε1 ,ε→0

τ →0

ENERGY EQUALITIES FOR COMPRESSIBLE NAVIER-STOKES EQUATIONS

19

4.4. Estimate of (H). By rewriting (H) as ˆ ˆ ˆ 2ν ε1 +ε3 t (̺u)ε − ̺ε uε (H) = − ∇ · (̺Du)ε dxdsdε2 ε3 ε1 ̺ε τ Ωε2 ˆ ˆ ˆ 1 ε1 +ε3 t ∇ · (̺Du)ε uε dxdsdε2 − 2ν ε3 ε1 τ Ωε2 ˆ ε1 +ε3 ˆ t ˆ 1 =: (H1) − 2ν (̺Du)ε uε n(θ)dHd−1 (θ)dsdε2 ε3 ε1 τ ∂Ωε2 ˆ ε1 +ε3 ˆ t ˆ 1 (̺Du)ε ∇uε dxdsdε2 + 2ν ε3 ε1 τ Ωε2 =: (H1) + (H2)bdr + (H3). Estimate (H1) using arguments similar to (73) and (75) we have lim sup lim sup |(H1)| = 0. ε→0

τ →0

The boundary term is computed using coarea formula (44) and Lemma 2.6 as follows ˆ ˆ 1 t lim sup lim sup |(H2)bdr | = lim sup lim sup 2ν (̺Du)u n(x)dxds ε ε1 ,ε→0 τ →0 ε1 ,ε→0 τ →0 3 τ Ωε1 \Ωε1 +ε3 1 k∇ukL2 ((Ω\Ωε3 )×(0,T )) kukL2 ((Ω\Ωε3 )×(0,T )) ε3 ≤ Ck∇uk2L2 ((Ω\Ω2ε )×(0,T )) . ≤C

(84)

3

2

1

Due to the assumption u ∈ L (0, T ; H (Ω)), by letting ε3 → 0, we get lim sup lim sup lim sup |(H2)bdr| = 0. ε3 →0

ε1 ,ε→0

τ →0

4.5. Conclusion of the Proof of Theorem 1.5. Collecting the estimates for (E), (F ), (G) and (H) and using similar arguments to subsection 3.5 we obtain the desired energy equality (14).  5. Proof of Theorem 1.7 The proof of Theorem 1.7 follows closely from that of Theorems 1.2 and 1.5 except we have to take extra care of the terms −2ν∆u and −λ∇(∇ · u). In the case Ω = Td , by multiplying the smoothed version of (cNS),

∂t (̺u)ε + ∇ · (̺u ⊗ u)ε + ∇(̺γ )ε − 2ν∆uε − λ∇(∇ · u)ε = 0

by (̺u)ε /̺ε we can proceed similarly to the proof of Theorem 1.2 except for the extra terms ˆ tˆ (̺u)ε (EX1) = −2ν ∆uε ε dxds (85) ̺ τ Td and (EX2) = −λ

ˆ tˆ τ

Td

∇(∇ ·

(̺u) u)ε ε ̺

ε

dxds.

By integration by parts, (EX1) = −2ν

ˆ tˆ τ

Td

ε

ε

∆u u dxds − 2ν

ˆ tˆ τ

Td

∆uε

(̺u)ε − ̺ε uε dxds ̺ε

(86)

20

Q-H. NGUYEN, P-T. NGUYEN, AND B. Q. TANG

= 2ν

ˆ tˆ τ

ε 2

Td

|∇u | dxds − 2ν

ˆ tˆ

Td

τ

∇ · (∇uε )

(̺u)ε − ̺ε uε dxds. ̺ε

By using an argument similar to that leading to (75) and the assumption u ∈ L2 (0, T ; H 1(Td )), we obtain ˆ tˆ 2 lim sup lim sup (EX1) − 2ν |∇u| dxds = 0. d ε→0 τ →0 0

T

For (EX2) we write similarly ˆ tˆ ˆ tˆ (̺u)ε − ̺ε uε ε ε (EX2) = −λ ∇(∇ · u) u dxds − λ dxds ∇(∇ · u)ε ̺ε τ Td τ Td ˆ tˆ ˆ tˆ (̺u)ε − ̺ε uε ε ε =λ (∇ · u) (∇ · u )dxds − λ ∇(∇ · u)ε dxds ̺ε τ Td τ Td and therefore obtain ˆ tˆ lim sup lim sup (EX2) − λ ε→0

τ →0

0

|∇ · u| dxds = 0. d 2

T

In the case Ω ⊂ Rd is a bounded domain with C 2 boundary, by proceeding as in section 4 and using the arguments dealing with (85) and (86), we are left to estimate only the extra boundary terms ˆ ˆ ˆ 1 ε1 +ε3 t (BEX1) = −2ν ∇uε uε n(θ)dHd−1 (θ)dsdε2, ε3 ε1 τ ∂Ωε2 and 1 (BEX2) = −λ ε3

ˆ

ε1 +ε3

ε1

ˆ tˆ τ

∂Ωε2

(∇ · u)ε uε n(θ)dHd−1 (θ)dsdε2 .

Both of these terms can be estimated by using an argument similar to the one leading to (84). Therefore we obtain lim sup lim sup lim sup |(BEX1)| = 0 ε3 →0

ε1 ,ε→0

τ →0

and

lim sup lim sup lim sup |(BEX2)| = 0. ε3 →0

ε1 ,ε→0

τ →0

Thus the proof of Theorem 1.7 is complete.



6. Proof of Theorem 1.8 Thanks to the proof of Theorems 1.2 and 1.5 we only need to take care of the terms concerning the scalar pressure. In case Ω = Td , we multiply ∂t (̺u)ε + ∇ · (̺u ⊗ u)ε + ∇P ε − 2ν∇ · (̺Du)ε = 0

by (̺u)ε /̺ε . Then we only need to take care of the following term (since other terms can be estimated similarly as in the previous Theorems) ˆ tˆ ˆ tˆ ˆ tˆ ε ε ε ε ε (̺u) − ̺ u ε (̺u) dxds = dxds + ∇P ∇P ε uε dxds ∇P ε ε ̺ ̺ d d d τ T τ T τ T ˆ tˆ ε ε ε (̺u) − ̺ u = ∇P ε dxds. ̺ε τ Td

ENERGY EQUALITIES FOR COMPRESSIBLE NAVIER-STOKES EQUATIONS

21

At the last step we have made used of the free divergence condition. Using H¨older’s inequality and Lemma 2.1 (ii), we have ˆ t ˆ ε ε ε ε ε ε (̺u) − ̺ u ε ≤ k∇P ε kL2 (Td ×(0,T )) k (̺u) − ̺ u kL2 (Td ×(0,T )) ∇P dxds ̺ε ̺ε (87) τ Td −1 ε ε ε ≤ CkP kL2 (Td ×(0,T )) ε k(̺u) − ̺ u kL2 (Td ×(0,T )) .

By assumption P ∈ L2 (Td × (0, T )) and Lemma 2.3 we obtain the desired limit ˆ t ˆ ε ε (̺u) =0 lim sup lim sup ∇P dxds ε ̺ d ε→0 τ →0 τ T and therefore finish the proof in the case Ω = Td .

In case Ω is a bounded domain, we need to take care of the boundary. Similar to the case of a torus, we only need to deal with the term ˆ ˆ ˆ ε 1 ε1 +ε3 t ε (̺u) ∇P dxds ε3 ε1 ̺ε τ Ωε2 ˆ ˆ ˆ ˆ ˆ ˆ ε ε ε 1 ε1 +ε3 t 1 ε1 +ε3 t ε (̺u) − ̺ u ∇P ∇P ε uε dxds = dxds + ε ε3 ε1 ̺ ε3 ε1 τ Ωε2 τ Ωε2 ˆ ˆ ε1 +ε3 ˆ t ˆ ˆ ˆ t ε +ε 1 3 (̺u)ε − ̺ε uε 1 1 ∇P ε P ε uε n(θ)dHd−1 (θ)dsdε2 dxds + = ε ε3 ε1 ̺ ε3 ε1 τ Ωε2 τ ∂Ωε2 =: (K1) + (K2)bdr . The term (K1) is estimated exactly as in (87) and therefore, lim sup lim sup |(K1)| = 0. ε→0

τ →0

For (K2)bdr we use the coarea formula, H¨older’s inequality and Lemma 2.6 to obtain ˆ ˆ 1 t lim sup lim sup |(K2)bdr | = lim sup lim sup P ε uε n(x)dxds ε ε1 ,ε→0 τ →0 ε1 ,ε→0 τ →0 3 τ Ωε1 \Ωε1 +ε3 1 kP kL2 ((Ω\Ωε3 )×(0,T )) kukL2 ((Ω\Ωε3 )×(0,T )) ε3 ≤ CkP kL2 ((Ω\Ωε3 )×(0,T )) k∇ukL2 ((Ω\Ω2ε3 )×(0,T )) . ≤

Since P ∈ L2 (Ω × (0, T )) and u ∈ L2 (0, T ; H 1(Ω)), by letting ε3 → 0, we derive lim sup lim sup lim sup |(K2)bdr | = 0. ε3 →0

ε1 ,ε→0

τ →0

Thus the proof of Theorem 1.8 is complete.



7. Proof of Theorem 1.9 The proof of Theorem 1.9 follows exactly from that of Theorem 1.8 except for the term relating to the pressure. However, since in this case we can take ̺ ≡ 1 and therefore (K1) = 0 trivially. From the estimate of (K2)bdr lim sup lim sup |(K2)bdr | ≤ CkP kL2(Ω\Ωε3 ×(0,T )) k∇ukL2 (Ω\Ω2ε3 ×(0,T )) ε1 ,ε→0

τ →0

22

Q-H. NGUYEN, P-T. NGUYEN, AND B. Q. TANG

we use u ∈ L2 (0, T ; H 1(Ω)) and the assumption (16) to conclude lim sup lim sup lim sup |(K2)bdr | = 0 ε3 →0

ε1 ,ε→0

τ →0

and therefore finish the proof of Theorem 1.9.



Acknowledgements: The third author is supported by the International Training Program IGDK 1754 and NAWI Graz. References [BT18] C. Bardos, E.S. Titi. Onsager’s conjecture for the incompressible Euler equations in bounded domains. Arch. Ration. Mech. Anal. 228 (2018), no. 1, 197-207. [BTW] C. Bardos, E.S. Titi, E. Wiedemann. Onsager’s Conjecture with Physical Boundaries and an Application to the Vanishing Viscosity Limit. arXiv:1803.04939. [BDIS15] T. Buckmaster, C. De Lellis, P. Isett, L.J. Sz´ekelyhidi. Anomalous dissipation for 1/5-H¨ older Euler flows. Ann. of Math. (2) 182 (2015), no. 1, 127-172. [BDSV18] T. Buckmaster, C. De Lellis, L. Sz´ekelyhidi Jr., V. Vicol. Onsager’s conjecture for admissible weak solutions, to appear in Comm. Pure Appl. Math. [CET94] P. Constantin, W. E, E.S. Titi. Onsager’s conjecture on the energy conservation for solutions of Euler’s equation. Comm. Math. Phys. 165 (1994), no. 1, 207-209. [CY] R.M. Chen, C. Yu. Onsager’s energy conservation for inhomogeneous Euler equations. arXiv:1706.08506. [DS10] C. De Lellis, L.J. Sz´ekelyhidi, On Admissibility Criteria for Weak Solutions of the Euler Equations. Arch. Rational Mech. Anal. 195 (2010), 225-260. [DS12] C. De Lellis, L.J. Sz´ekelyhidi, The h-principle and the equations of fluid dynamics. Bull. Amer. Math. Soc. (N.S.) 49 (2012), no. 3, 347-375. [DS13] C. De Lellis, L.J. Sz´ekelyhidi. Dissipative continuous Euler flows. Invent. Math. 193 (2013), no. 2, 377-407. [DS14] C. De Lellis, L.J. Sz´ekelyhidi. Dissipative Euler flows and Onsager’s conjecture. J. Eur. Math. Soc. 16 (2014), no. 7, 1467-1505. [DN18] T.D. Drivas, H.Q. Nguyen. Onsagers conjecture and anomalous dissipation on domains with boundary. to appear in SIAM J. Math. Anal. (2018) arXiv:1803.05416v2. [Eyi94] G. L. Eyink. Energy dissipation without viscosity in ideal hydrodynamics: I. Fourier analysis and local energy transfer. Phys. D 78 (1994), 222-240. [Ise18] P. Isett, On the endpoint regularity in Onsager’s conjecture. Preprint, 2017. arXiv:1706.01549. [Ise18a] P. Isett. A proof of Onsgar’s conjecture. to appear in Ann. Math. [BDIS15] T. Buckmaster, C. De Lellis, P. Isett, L.J. Sz´ekelyhidi. Anomalous dissipation for 1/5-H¨ older Euler flows. Ann. of Math. (2) 182 (2015), no. 1, 127-172. [LV18] I. Lacroix-Violet, A. Vasseur. Global weak solutions to the compressible quantum NavierStokes equation and its semi-classical limit. Journal de Math´ematiques Pures et Appliqu´ees 114 (2018): 191-210. [NN18] Q.-H. Nguyen, P.-T. Nguyen. Onsager’s conjecture on the energy conservation for solutions of Euler equations in bounded domains. to appear in Journal of Nonlinear Science (2018). [NNT18] Q.-H. Nguyen, P.-T. Nguyen, B. Q. Tang. Energy conservation for inhomogeneous incompressible and compressible Euler equations. Preprint, 2018. [Ons49] L. Onsager. Statistical Hydrodynamics. Nuovo Cimento (Supplemento), 6 (1949), 279-287. [Sch93] V. Scheffer. An inviscid flow with compact support in space-time. J. Geom. Anal. 3 (1993), 343-401 [Shn97] A. Shnirelman. On the nonuniqueness of weak solution of the Euler equation, Comm. Pure Appl. Math. 50.12 (1997) 1261-1286. [Ser62] J. Serrin. The initial value problem for the NavierStokes equations. Nonlinear Problems. Proceedings of the Symposium, Madison,Wisconsin, 1962. University of Wisconsin Press, Madison, Wisconsin, 6998, 1963. [Shi74] M. Shinbrot. The energy equation for the NavierStokes system. SIAM J. Math. Anal. 5 (1974) 948–954. [VY16] A. Vasseur. C. Yu, Existence of global weak solutions for 3D degenerate compressible Navier-Stokes equations. Invent. Math. (2016) 206:935–974. [Yu17] C. Yu. Energy conservation for weak solutions of the compressible Navier-Stokes equations. Arch. Rational Mech. Anal. 225 (2017) 1073–1087. [Yu17a] C. Yu. The energy conservation for the Navier-Stokes equations in bounded domains. Preprint 2017. [Yu16] C. Yu. A new proof of the energy conservation for the Navier-Stokes equations. Preprint, 2016.

ENERGY EQUALITIES FOR COMPRESSIBLE NAVIER-STOKES EQUATIONS

23

Quoc-Hung Nguyen Department of Mathematics, New York University Abu Dhabi, Abu Dhabi, United Arab Emirates E-mail address: [email protected] Phuoc-Tai Nguyen Faculty of Science, Department of Mathematics and Statistics, Masaryk University, 61137 Brno, Czech Republic E-mail address: [email protected], [email protected] Bao Quoc Tang Institute of Mathematics and Scientific Computing, University of Graz, Heinrichstrasse 36, 8010 Graz, Austria E-mail address: [email protected]