arxiv: v1 [math.ap] 28 Dec 2018

Transcription

1 GLOBAL WELL-POSEDNESS AND TIME-DECAY ESTIMATES OF THE COMPRESSIBLE NAVIER-STOKES-KORTEWEG SYSTEM IN CRITICAL BESOV SPACES arxiv: v1 [math.ap] 8 Dec 18 NOBORU CHIKAMI AND TAKAYUKI KOBAYASHI Abstract. We consider the compressible Navier-Stokes-Korteweg system describing the dynamics of a liquid-vapor mixture with diffuse interphase. The global solutions are established under linear stability conditions in critical Besov spaces. In particular, the sound speed may be greater than or equal to zero. By fully exploiting the parabolic property of the linearized system for all frequencies, we see that there is no loss of derivative usually induced by the pressure for the standard isentropic compressible Navier-Stokes system. This enables us to apply Banach s fixed point theorem to show the existence of global solution. Furthermore, we obtain the optimal decay rates of the global solutions in the L (R d -framework. 1. Compressible Navier-Stokes-Korteweg system We consider the following barotropic compressible Navier-Stokes-Korteweg system endowed with an internal local capillarity: t ρ + div (ρu =, t (ρu + div (ρu u + (P (ρ = div (τ(ρ, u + K(ρ, (1.1 (ρ, u t= = (ρ, u, where (t, x R + R d, d, ρ = ρ(t, x R + and u = u(t, x R d are unknowns representing the fluid density and the velocity vector field, respectively. The viscous strain tensor τ, in general, is given by τ = τ(ρ, Du = µ(ρd(u + λ(ρdiv u Id, with Id denoting an identity matrix and and D(u designating the deformation tensor defined by D(u := 1 ( u + t ( u with ( u ij := i u j, where t ( u denotes the transposed matrix of u. For simplicity we assume that the viscosity coefficients λ and µ are given constants satisfying µ > and λ + µ >. We write Lu = div (µdu + (λdiv u. When µ and λ are constants, then we see that div (τ = Lu = µ u + (µ + λ div u. The pressure is assumed to be barotropic, i.e., P is a sufficiently smooth function of ρ. In the general Korteweg system, the capillary term is given by div (K(ρ, where K(ρ is the Korteweg tensor defined by K(ρ = κ(ρ ( ρ ρ Id κ(ρ ρ ρ, Date: December 31, 18. 1

2 N. CHIKAMI AND T. KOBAYASHI where ρ ρ stands for the tensor product ( j ρ k ρ jk. The coefficient κ may depend on ρ in general. Here, we assume the case of positive constant coefficient κ(ρ κ. We assume that there is no vacuum, and ρ tends to a positive constant ρ at spatial infinity. The system describes the dynamics of a liquid-vapor mixture in the setting of the Diffuse Interface approach: between the two phases lies a thin region of continuous transition and the phase changes are described through the variations of the density (with for example a Van der Waals pressure. Historically, the system is first proposed by J. D. Van der Waals (1893, and later Korteweg tensor is introduced by D. J. Korteweg (191. The system considered here is a version reformulated by J. E. Dunn and J. Serrin [13] in An important features of System (1.1 is that the pressure is non-monotone in general. Simply put, P (ρ is no longer a positive number, which inevitably leads to linear instability for some value of ρ. As such, most mathematical results (in the whole space until now deal with the condition that P (ρ is strictly positive, except for local-in-time results. There are numerous works dedicated to the study of the system with the capillarity term. For known well-posedness results for System (1.1, global smooth solutions are established by [14, 15]. Decay estimates along with the analysis of diffusive wave phenomenon are considered by T. Kobayashi and K. Tsuda [16]. The latter also obtains the existence of time-periodic solutions along with their stability in []. See also [17, 19] for the existence and decay results. In the critical framework, R. Danchin and B. Desjardins [11] first obtain the well-posedness in critical Besov spaces. Especially, they consider a general non-monotone pressure to show local existence of the solution. For results on non-local capillary terms and convergence to various models, we refer to the works by F. Charve and B. Haspot [4 6]. In [18], T. Tang and H. Gao obtain the global solution for System (1.1 in T d under a non-monotone pressure law. Two difficulties usually arise in the analysis of the compressible viscous fluids. One is the loss of derivative induced by the pressure term and the other is the appearance of the second order nonlinear term. Both of them preclude the application of standard semi-group theory and necessitate the laborious energy estimates. In contrast to the compressible viscous fluids without the capillary term, however, System (1.1 is known to have one mathematical upside, which is that there is no longer a loss of derivative due to the parabolic smoothing effect for both the velocity and the density. We may check that a scaling invariance holds for System (1.1: namely a family of rescaled functions (ρ, u (ρ ν, u ν with ρ ν (t, x = ρ(ν t, νx and u ν (t, x = νu(ν t, νx (1. leaves (1.1 invariant provided that P is changed into ν P. The pressure term, to some extent, can be regarded as a lower order term in our setting and therefore it is possible to solve the system in the scale-critical function spaces, i.e., the spaces in which the norm is kept invariant under the transforms (1.. Moreover, we clarify in this paper that it suffices to assume that the pressure is merely non-negative, i.e., P (ρ, in order to obtain global solutions in critical Besov spaces. We establish a priori estimates for the linearized system that is uniform with respect to the parameter P (ρ (See Lemma.5. We also show the decay estimates of such global small solutions under additional L 1 (R d -type regularity assumption. We remark that although we focus our analysis in the L (R d -setting, it is not difficult to extend our decay results to the L p (R d -framework in the spirit of [1]. In fact, F. Charve, R. Danchin and J. Xu have recently shown the existence of solution for System

3 WELL-POSEDNESS AND TIME-DECAY ESTIMATES OF CNSK 3 (CNSK in L p (R d -framework along with the Gevrey analyticity of the solution [7]. As a corollary of the analyticity, they obtain decay estimates of higher order norms with respect to a critical data. Note that contrary to [7] our decay estimates deal with the decay rates of critical norms with respect to a L 1 (R d data, which naturally requires more labour and cannot be deduced from their results. Furthermore, our results for the critical case P (ρ = seems completely new Notations and function spaces. Hereafter, we denote by L p (R d, 1 p, the standard Lebesgue spaces on the d-dimensional Euclidean space R d, and by l p (Z the set of sequences with summable p-th powers. Unless otherwise stated, we assume d throughout the paper. We abbreviate L p (R d = L p whenever the dimension d is not relevant, and denote the time-space Lebesgue norm by L q (,T ;L p = L q T Lp. An infinitely differentiable, complex-valued function f on R d is called a Schwartz function if for every pair of multi-indices α and β there exists a positive constant C α,β such that sup x R d x α β f = C α,β <. We denote the set of all Schwartz functions on R d by S(R d and its topological dual by S (R d. Elements in S (R d are called tempered distributions. We define the fractional derivative operator by Λ s := ( s/ = F 1 s F for s R, where Fu = û denotes the Fourier transform of u S (R d. We define a subspace of S (R d which lays a basic foundation of our analysis. For the details, we refer to []. Definition 1.1 ([] P. We denote by S h (Rd the space of tempered distributions u such that lim λ F 1 [θ(λ û] L = for any θ in S(R d. In what follows, we denote by C a generic constant that may change from line to line. We write A = B whenever there exists a positive constant c such that c 1 A B ca. We also denote A B whenever there exists arbitrary harmless constant C such that A CB. Constants depending on parameters a, b,..., c are denote by C(a, b,..., c. We agree that p 1 always denotes the Hölder conjugate of p, i.e., p + 1 p = 1, with the convention 1 =. Vectors in this article are understood as column vectors in general; however, we omit the sign of transpose t whenever it is clear from the context or simply irrelevant. We express a function u of variables (t, x as u, u(t or u(t, x depending on circumstances Besov spaces. Let {φ j } j Z be a homogeneous Littlewood-Paley dyadic decomposition. Namely, let φ S(R d be a non-negative, radially symmetric function that satisfies supp φ {ξ R d ; 1 ξ } and φ( j ξ = 1 for all ξ. Setting φ j (ξ := φ( j ξ and Φ j 1 (ξ := 1 k j φ k (ξ for j Z, we have j Z supp Φ j 1 {ξ R d ; ξ j } and supp φ j {ξ R d ; j 1 ξ j+1 } for all j Z. Given an element u in S h (Rd and a (homogeneous Littlewood-Paley decomposition {φ j } j Z, we define the (homogeneous Littlewood-Paley localization operator by j u := F 1 [ φ j û] (1.3

4 4 N. CHIKAMI AND T. KOBAYASHI and the frequency cut-off operator by Ṡju := j j j u = F 1 [ Φ j û] for j Z. Definition 1. (Homogeneous Besov space. For s R and 1 p, σ we define the homogeneous Besov spaces Ḃs p,σ(r d as follows: Ḃp,σ(R s d := {u S h (Rd ; u Ḃs p,σ (R d < }, u Ḃp,σ(R s d { := js } j u L p. For notational convenience, we abbreviate Ḃs p,σ (R d = Ḃ s p,σ. More properties of Besov spaces may be found in e.g. []. j Z l σ For some fixed j Z, we denote by u h := Ṡj u the low frequencies of u, and by u h := u u h the high frequencies of u. We also need the notation u l := js Ḃ s j u L and u h := jσ Ḃ σ j u L. (1.4 j j j j 1 The small overlap between low and high frequencies are to ensure that u l Ḃs C u l Ḃ s and u h Ḃσ C u h Ḃ σ. Lastly, we introduce the so-called Chemin-Lerner spaces. Definition 1.3 (Chemin-Lerner spaces, [, 9]. For t 1, t >, s R, 1 σ, ρ, we set u Lρ (t 1,t ;Ḃs p,σ := { js j u L ρ (t 1,t ;L p } j Z l σ. We then define the space L ρ (t 1, t ; Ḃs p,σ(r d as the set of tempered distributions u over (t 1, t R d such that lim j Ṡ j u = in S ((t 1, t R d and u Lρ (t 1,t ;Ḃs p,σ <. When t 1 = and t = T, we denote L ρ T (Ḃs p,σ(r d = L ρ (, T ; Ḃs p,σ(r d with the shortened expression of the norm Lρ T (Ḃs p,σ. In what follows, we agree that C([, T ]; Ḃs := { v C([, T ]; Ḃs ; v L T (Ḃs < } ( Momentum formulation and main results. Setting m := ρu, we recast System (1.1 in the momentum formulation, which reads as follows: t ρ + div m =, t m + div (ρ 1 m m + P = L(ρ 1 m + div (K(ρ, (CNSK (ρ, m t= = (ρ, m. The clear difference from the velocity formulation is that now the mass conservation law is linear while the capillary term is nonlinear. Also notable is that there is ρ 1 inside the second-order term, thereby requiring more regularity to the density a priori. This poses no obstacle in the actual analysis since, as is well-known, there is a parabolic smoothing effect for the density in System (CNSK. In this paper, we perform our analysis in the momentum formulation. In fact, the momentum formulation turns out to be essential for our analysis for the case γ = ; Particularly, the fact the the nonlinear terms may be written in divergence form greatly helps our analysis in d = 3, as we see in the proof of Lemma 3.4.

5 WELL-POSEDNESS AND TIME-DECAY ESTIMATES OF CNSK Global existence. Throughout this paper, we set γ := P (ρ. We first state a result of well-posedness for d and γ, which is contained in [11]. Under smallness of the initial data, the solution is globally extended when γ >. Theorem 1.4 ([11]. Let d and let P be a smooth function such that γ. For sufficiently small initial data such that ((γ + Λ(ρ ρ, m Ḃ d 1 (Rd, there exists a time T > and a unique solution (ρ, m to (CNSK satisfying ((γ + Λ(ρ ρ, m C([, T ]; Ḃ d 1 (Rd L 1 (, T ; Ḃ d +1 (Rd. Moreover, if γ > then T = and there exists a constant C depending on d, µ, λ, ρ and P such that ((γ + Λ(ρ ρ, m L (, ;Ḃ d + Λ ((γ + Λ(ρ ρ 1, m L 1 (, ;Ḃ d 1 C ((γ + Λ(ρ ρ, m d Ḃ 1. Remark 1.5. As observed in [11], as far as a local solution is concerned, the smallness of the initial data and the condition γ may be removed. The question remains whether we may obtain a global solution when the sound speed is zero, i.e., when γ =. It turns out that restricting the dimension to d 3, we have the global well-posedness for the case under an additional low-frequency assumption. Theorem 1.6. Let d 3 and let P be a smooth function such that γ =. For sufficiently small initial data such that (Λ(ρ ρ, m Ḃ d 3 (Rd Ḃ d 1 (Rd, there exists a unique solution (ρ, m to (CNSK satisfying (Λ(ρ ρ, m C([, ; Ḃ d 3 (Rd Ḃ d 1 (Rd L 1 (, ; Ḃ d 1 (Rd Ḃ d +1 (Rd. Moreover, there exists a positive constant C depending on d, µ, λ, ρ and P such that (Λ(ρ ρ, m L (, ;Ḃ d + Λ ((Λ(ρ ρ 1, m L1 (, ;Ḃ d 1 C ((Λ(ρ ρ, m d Ḃ 1. Remark 1.7. In both Theorem 1.4 and Theorem 1.6, the space C(, T ; Ḃs (Rd is defined by Definition 1.5. Let us remark that although we focus our analysis for the barotropic case, our technique can be also applied to non-isentropic cases. The results can also be naturally extended in the settings of torus and variable capillarity and viscosity coefficients Decay estimates. Next, we clarify the decay rates of the solution constructed by Theorem 1.4 and Theorem 1.6.

6 6 N. CHIKAMI AND T. KOBAYASHI We first state a result for the case γ >. Setting t := (1 + t, we introduce an auxiliary norm D(t := sup τ 1 (s+ d ((γ + Λ(ρ ρ, m l s ( d, d +1] L t (Ḃs + τ α Λ ((γ + Λ(ρ ρ, m h L t (Ḃ d 1. (1.6 Here and in what follows, we agree that the notation τ σ f, σ R, designates the function (τ, x τ σ f(τ, x. The notations l and are defined analogously L t (Ḃs h L t (Ḃ d 1 to (1.4. Under additional assumptions on the low frequencies of the initial data, one may obtain time-decay estimates that are reminiscent of those of the standard compressible Navier- Stokes equations (see [1]. Theorem 1.8. Let γ >. Let the data (a, m satisfy the assumptions of Theorem 1.4. There exist a frequency-threshold j only depending on d, µ, λ, κ and γ and a positive constant c so that if in addition D := ρ ρ l Ḃ d + m l Ḃ d c,, then the global solution (a, m given by Theorem 1.4 satisfies for all t, D(t C ( D + (Λ(ρ ρ, m h Ḃ d 1 where D(t is defined by (1.6 and α := d + 1 ε for sufficiently small ε >. We next state the result for the case γ =. We set (1.7 D(t := sup τ 1 (s+ d (Λ(ρ ρ, m l s ( d, d +1] L t (Ḃs + τ α Λ (Λ(ρ ρ, m h L. t (Ḃ d 1 (1.8 For the case γ =, we have the following. Theorem 1.9. Let γ =. Let the data (ρ, m satisfy There exist a frequency-threshold j constant c so that if in addition (Λ(ρ ρ, m Ḃ d, (Rd Ḃ d 1 (Rd. (1.9 D := (ρ ρ l Ḃ 1 d, only depending on d, µ, λ and κ and a positive + m l c Ḃ d, then the global solution ((ρ ρ, m given by Theorem 1.6 satisfies for all t, D(t C ( D + (Λ(ρ ρ, m h Ḃ d 1 where D(t is defined by (1.8 and α := d + 1 ε for sufficiently small ε >. (1.1 Remark 1.1. Note that since d 3, we have d d 3, which implies u l Ḃ d 3 u l for d 3. Thus, the regularity assumption (1.9 in Theorem 1.9 is stronger than Ḃ d,

7 WELL-POSEDNESS AND TIME-DECAY ESTIMATES OF CNSK 7 that of Theorem 1.6, and the corresponding small global solution is bound to naturally arise. The rest of the paper is organized as follows: in the second section, we study the linearized problem of (CNSK. The third section is devoted to the proof of Theorem 1.4 and Theorem 1.6. In the fourth section, we prove Theorem 1.8 and Theorem 1.9. In the Appendix, we collect some properties of Besov spaces.. Linearized problem Linearization around a constant state (ρ, gives t a + div m =, t m 1 Lm κρ a + γ a = L (Q(am + div (K(a. ρ 1 div (m m div (Q(am m + ( P (a + ρ γ a, ρ (.1 where we put γ := P a (ρ, a := ρ ρ and Q(a := ρ (a+ρ. We split m into compressible and incompressible components, i.e., m = w + P u, with w := Pm where P and P are the projectors onto divergence-free and potential vector-fields, respectively ( w := (Id + div ( 1 m. Setting v := Λ 1 div m = Λ 1 div P m, we have t a + Λv =, t v ν v γλa κρ Λ 3 a =, (. t w µ w =, where ν := ν ρ, ν := µ + λ, µ := µ ρ. (.3 We investigate the properties of the above linearized system..1. Maximal regularity estimates in Besov spaces. Let us consider the following standard heat equation with a constant diffusive coefficient ν > : { t u ν u = f, t >, x R d, (.4 u t= = u, x R d, where the outer force f : R + R d R and the initial data u : R d R are given functions or tempered distributions. Recall the following estimate for a linear heat equation, which can be found in literatures such as [] P157 and [1] P7. Lemma.1. Let T >, s R and 1 p, r, σ. Assume that u Ḃs p,σ(r d and f L r T (Ḃs + r p,σ (R d. Then (.4 has a unique solution u in L T (Ḃs p,σ(r d L r T (Ḃs+ r p,σ (R d and there exists a constant C depending only on q and r such that for all q [r, ], we have ( ν 1 q u C u Lq T (Ḃs+ q Ḃs + ν 1 r 1 f. p,σ p,σ Lr T (Ḃs + r p,σ If in addition σ is finite, then u belongs to C([, T ; Ḃs p,σ(r d. Our aim here is to prove a similar estimate for the linearized system of (CNSK.

8 8 N. CHIKAMI AND T. KOBAYASHI Two difficulties usually arise in the analysis of the compressible viscous fluids. One is the loss of derivative induced by the pressure term and the other is the appearance of the second order nonlinear term. Both of them preclude the application of standard semigroup theory and necessitate laborious energy estimates. In contrast to the compressible viscous fluids without the capillary term, System (CNSK is known to have a pleasant mathematical upside, which is that there is no longer a loss of derivative due to the parabolic smoothing effect for the density. From hereon, we consider a linear hyperbolic-parabolic equation t a + div m = f, t m 1 (.5 Lm + γ a κρ a = g ρ with initial data (a, m and derive an a priori estimate in Besov spaces using the method of Lyapunov in the frequency space. Lemma.. Let T >, s R, 1 r, σ, γ, ((γ + Λa, m Ḃs,σ(R d and ((γ + Λf, g L r (, T ; Ḃs + r,σ (R d. Suppose that (a, m is a solution for (.5 such that ((γ + Λa, m L q (, T ; Ḃs+ q,σ (Rd for q [r, ] with initial data (a, m. Then there exists some positive constant C depending only on ρ, λ, µ, κ, d, q, r and σ such that for all q [r, ], the following estimate hold: ( ((γ + Λa, m C ((γ + Λa Lq T (Ḃs+ q, m Ḃs + ((γ + Λf, g.,σ,σ Lr T (Ḃs + r,σ for all < t T. Remark.3. Note that the above constant does not depend the parameter γ and holds for both cases γ > and γ =. It is clear from the above estimates that one needs to impose a stronger control on the low frequencies at a linear level when γ >, which is in accordance with the known result for the isentropic compressible fluid. Proof of Lemma.. The argument below is inspired by [1]. However, the resulting smoothing effects, as expected, are strikingly better than those of the compressible Navier- Stokes system. After projecting m onto the compressible part by v = Λ 1 div m as before, we first treat the homogeneous system of the first two equations of System (. : { t a + Λv =, t v ν v γλa κρ Λ 3 (.6 a = where γ. Taking the Fourier transform, System (.6 is transformed into [ ] [ ] [ ] d â â = dt v Â(ξ ξ with Â(ξ := v (γ + κρ ξ ξ ν ξ. The characteristic equation for the matrix A(ξ is X + ν ξ X + (γ + κρ ξ ξ =,

9 WELL-POSEDNESS AND TIME-DECAY ESTIMATES OF CNSK 9 and Â(ξ has two distinct eigenvalues given by: λ ± = λ ± (ξ := ν (1 ξ ± 1 4(γ ξ + κρ ξ 3 ξ ν ξ 4 = ν (1 ξ ± 1 4κρ ν 4γ ν ξ Recall that γ. From the above formula, we readily gather the following: (1 If ν < 4κρ 3, then both λ ± (ξ are two complex conjugated eigenvalues λ ± (ξ = ν 4κρ (1 ξ ± is( ξ with S( ξ = ν 1 + 4γ ν ξ, regardless of the choice of γ, for all frequencies. If γ >, we have λ ± (ξ = ν (1 ξ ± is( ξ = ν ξ i γ ξ + io( ξ 3 as ξ and λ ± (ξ = ν (1 ξ ± is( ξ = ν ξ i ν as ξ.if γ =, we have λ ± (ξ = ν ( 4κρ 1 ± i ν 1 ξ. ( 1 4κρ ν 1 ξ + io(1 Thus, we may expect parabolic smoothing for both components throughout all frequencies. ( If ν = 4κρ 3 and γ >, then both λ ± (ξ are two complex conjugated eigenvalues λ ± (ξ = ν 4γ (1 ξ ± i ν ξ = ν ( ξ 1 ± i γ = ν ν ξ ξ ± i γ ξ for all frequencies. If ν = 4κρ 3 and γ =, then we have a (real double root λ ± (ξ = ν ξ. (3 If ν > 4κρ 3 and γ > then we have two regions of distinct behaviors. (i In the low frequency ( ξ < 4γρ, λ ν 4κρ 3 ± (ξ are two complex conjugated eigenvalues λ ± (ξ = ν ( 4γ (1 ξ ± is( ξ with S( ξ = ν ξ 1 4κρ ν whose real parts are negative and λ ± (ξ = ν ( ρ γ (1 ξ ± i ν ξ ν 4κρ 3 4ρ ν γ ξ + O( ξ 3 = ν ξ i γ ξ + io( ξ 3 as ξ.

10 1 N. CHIKAMI AND T. KOBAYASHI 4γρ, the two real eigen- ν 4κρ 3 (ii On the other hand, in the high frequency ( ξ > values can be expressed as λ ± (ξ = ν (1 ξ ± R( ξ with R( ξ = 1 4κρ ν 4γ ν ξ whose real parts are negative and λ ± (ξ = ν ( 1 (1 ξ ± 1 4κρ3 ν + O(1 as ξ If ν > 4κρ 3 and γ =, then both λ ± (ξ are real and purely parabolic: λ ± (ξ = ν ( 1 ± 1 4κρ ν ξ We observe that there are two striking differences from the barotropic compressible fluid: One is that the the system is linearly stable even when the sound speed is zero, i.e., γ =. Another is that the solutions in all frequency modes likely exhibit parabolic smoothing effect. These two features directly stem from the presence of the Korteweg tensor. According to the above ansatz, we seek to derive the natural regularities for (â, v by the construction of the Lyapunov function in the Fourier space. From (.7, we immediately obtain the following sets of equalities: 1 d dt â + ξ R( v â =, (.7 1 d dt v + ν ξ v (γ + κρ ξ ξ R( v â =, (.8 and d dt R(â v + ξ v (γ + κρ ξ ξ â + ν ξ R( v â =, (.9 where Rz denotes the real part of a complex number z and ( is the inner product in C (that is (z 1 z := z 1 z with z denoting the complex conjugate of z. The proof below holds for γ. We fix some real parameter η >. By adding (γ + κρ ξ (.7, (.8 and η ξ (.9, we discover that 1 d ( (γ + κρ ξ â + v η ξ R(â v dt + (ν η ξ v + η(γ + κρ ξ ξ â ην ξ 3 R(â v =. (.1 By Young s inequality, for all ε > we obtain the following lower bound for the last three terms of the above: (ν η ξ v + η(γ + κρ ξ ξ â ην ξ 3 R(â v ( ( ν η 1 + ν ξ v + η(γ + (κρ εν ξ ξ â. ε for all frequencies. Taking ε = κρ3 ν, we have κρ εν = κρ. Hence, taking η such that ( ρ η < ν + ν 1 κρ 3,

11 we obtain WELL-POSEDNESS AND TIME-DECAY ESTIMATES OF CNSK 11 (ν η ξ v + η(γ + κρ ξ ξ â ην ξ 3 R(â v C (ν, κ, ρ ξ ( (γ + ξ â + v for some constant C (ν, κ, ρ depending only on ν, κ and ρ. In the same way, there exists a positive constant C 1 (ν, κ, ρ such that C 1 (ν, κ, ρ 1 ( (γ + ξ â + v (γ + κρ ξ â + v η ξ R(â v C 1 (ν, κ, ρ ( (γ + ξ â + v. (.11 Indeed, the upper bound is trivial; In order to prove the lower bound, we use Young s inequality to show ( (γ + κρ ξ â + v η ξ R(â v γ â + (κρ ηε ξ â + 1 η ε v for all ε >. Therefore, taking ε = κρ 4η and η κρ for example, we obtain the lower bound, and thus, (.11. The above computations combined with (.1 imply 1 d dt L + ξ (γâ, ξ â, v with L := (γ + κρ ξ â + v η ξ R(â v. By (.11, there exists a constant c > independent of ξ such that 1 d dt L + c ξ L, which implies the exponential stability for all t. This combined with (.11 again gives L (t e c ξ t L (, (.1 (γâ, ξ â, v (t e c ξ t (γâ, ξ â, v (. Granted with the above estimate in the frequency space, it is now straightforward to complete the proof of the lemma : it suffices to multiply the Littlewood-Paley decomposition φ j, perform exactly the same computations as above on to obtain a bound similar to (.1: with Similarly to (.11, we have so we have obtained (γ φ j â, j φj â, φ j v L (t e cjt L (, L := (γ + κρ j φ j â + φ j v η j R( φ j â φ j v. L = (γ φj â, j φj â, φ j v, (γ φ j â, j φj â, φ j v (t e cjt (γ φ j â, j φj â, φ j v (. Then taking L (R d norms of both sides of the resulting inequality, we obtain (γ φ j â, j φj â, φ j v L e c jt (γ φ j â, j φj â, φ j v L.

12 1 N. CHIKAMI AND T. KOBAYASHI For γ >, this implies that for any integer m there exist constants depending only on m such that the action of the semigroup e tâ is given by e tâ(γ φ j â, φ j v L Ce c jt (γ φ j â, φ j v L for j m (.13 and e tâ( j φj â, φ j v L Ce c jt ( j φj â, φ j v L for m < j, (.14 where (â, v = (â, v(. Thanks to (.13 and (.14, we deduce e tâ(γ φ j â, j φj â, φ j v L Ce c jt (γ φ j â, j φj â, φ j v L for all j Z. (.15 To handle the source term, splitting the frequency regions into low and high and resorting to Duhamel s formula with the aid of (.13 and (.14, we obtain (γ φ j â, j φj â, φ j v L p e c jt (γ φ j â, j φj â, φ j v L + e c j (t τ (γ φ j f, j φj f, φj ĝ L dτ. Then, time-integration along with the convolution inequality leads to a frequency-localized time-space estimate. Summing up the resulting inequality over j Z with a weight js, we obtain the estimate for the compressible part. Since the incompressible component of the velocity w := (Id Λ 1 div u satisfies a mere heat equation t w µ w = (Id Λ 1 div g, the standard parabolic estimate (Lemma.1 for w yields the desired estimates for w as well. As a corollary of the proof of Lemma., we may deduce the following estimate which is used to establish decay estimates. Lemma.4. Let d, s 1, s R and s 1 > s. Let e tk be the semigroup associated to System (.5 with f g and initial data (a, m Ḃs, (Rd. Then there exists a constant C depending on d, s 1 and s such that the following estimate holds: e tk ((γ + Λa, m Ḃs1 Ct s 1 s ((γ + Λa, m Ḃs, for all t >. Proof. From (.15 and the corresponding estimate for w, we readily deduce e tk ( j (γ + Λa, j m L Ce c jt ( j (γ + Λa, j m L for all j Z. By Lemma 5. with r 1, we readily have { js 1 e tk } ( j (γ + Λa, j m L j Z l 1 = Ct s 1 s { } (t 1 j (s 1 s e c j t ((γ + Λa, m Ḃs j Z, l 1 Ct s 1 s ((γ + Λa, m Ḃs., 3. Nonlinear problem In this section, we show the existence and uniqueness of global solution to the nonlinear problem (.1.

13 WELL-POSEDNESS AND TIME-DECAY ESTIMATES OF CNSK Banach s fixed point theorem. We first prove the following auxiliary lemma, which is a generalization of Lemma.1 in [8] to bi- and-trilinear operators. Lemma 3.1. Let (X, X be a Banach space. Let B : X X X be a bilinear continuous operator with norm K and T : X X X X be a trilinear continuous operators with norm K 3. Let further L : X X be a continuous linear operator with norm N < 1. Then for all y X such that { 1 N (1 N } y X < min,, (K + 3K 3 Equation x = y + L(x + B(x, x + T (x, x, x has a unique solution x in the ball B X R ( of center and radius R { } 1 N := min 1, K +3K 3. In addition, x satisfies x X 1 N y X. (3.1 Proof. Although this is standard, we prove the lemma for completeness. Letting Φ : x y + L(x + B(x, x + T (x, x, x, we prove that Φ is a contraction map in some ball. Let I := y X and let x be in some closed ball BR X ( of X. Assuming R 1, we have Let R be defined by Φ(x X y X + L(x X + B(x, x X + T (x, x, x X I + NR + (K + K 3 R. Then we have from (3. ( 1 N Φ(x X Thus, suffices to take R such that R min for Φ to be a self-map on the closed ball B X R (. (3. R := I. (3.3 1 N + N + (K + K 3 R R. { } 1 N 1,. (3.4 (K + K 3 On the other hand, we have for any x 1, x BR X ( and for R 1, Φ(x 1 Φ(x X L(x 1 x X + B(x 1 x, x 1 X + B(x 1, x 1 x X + T (x 1 x, x 1, x 1 X + T (x 1, x 1 x, x 1 X + T (x, x, x 1 x X (N + K R + 3K 3 R x 1 x X (N + (K + 3K 3 R x 1 x X. Hence, in order for Φ to be a contraction, N + (K + 3K 3 R has to hold; this amounts to taking R such that { } 1 N R < min 1,. (3.5 K + 3K 3 Thus, from (3.3, (3.4 and (3.5, Banach s fixed point theorem is applicable for all y X such that { } { 1 N y 1 N 1 N (1 N } X < min 1,, i.e., y X < min,, K + 3K 3 (K + 3K 3

14 14 N. CHIKAMI AND T. KOBAYASHI The uniqueness of the fixed point in the ball naturally follows. Routine computations as in the stability estimate lead to ( Global results in L (R d -based spaces when γ >. Lemma. yields the following. Corollary 3.. Let d, T >, s R, γ and ((γ + Λf, g L 1 (, T ; Ḃs,σ(R d. Assume that (a, m is a solution for (.5 such that ((γ + Λa, m L (, T ; Ḃs,σ(R d L 1 (, T ; Ḃs+,σ (Rd. Then there exists some positive constant C depending only on ρ, P, λ, µ, κ and d such that the following estimate hold : if γ, then for all < t T. ((γ + Λa, m L t (Ḃs,σ + Λ ((γ + Λa, m L1 t (Ḃs,σ C ( ((γ + Λa, m Ḃs,σ + ((γ + Λf, g L1 t (Ḃs,σ Given initial data (a, m such that ((γ + Λa, m (.1, it suffices to show that the map Ḃ d 1 (Rd, in order to solve Φ : (b, n (a, m (3.6 with (a, m the solution to t a + div m =, t m 1 Lm κρ a + γ a = N(b, n, ρ (a, m t= = (a, m, where N(b, n := L (Q(bn + κb b 1 ρ div (n n div (Q(bn n + ( P (b + ρ γ b, has a fixed point in a ball of a suitable metric space. We define a Banach space CL T by equipped with a norm (a, m CLT CL T := L (, T ; Ḃ d 1 (Rd L 1 (, T ; Ḃ d +1 (Rd := ((γ + Λa, m L T (Ḃ d 1 + Λ ((γ + Λa, m L1 T (Ḃ d 1. By interpolation, we may readily see that (a, m CL T further satisfies (3.7 (3.8 ((γ + Λa, m L d q (, T ; Ḃ 1+ q (R d for q [1, ]. (3.9 Recall that N defined in (3.8 is a combination of bi-and-trilinear terms. We have the following nonlinear estimates.

15 WELL-POSEDNESS AND TIME-DECAY ESTIMATES OF CNSK 15 Lemma 3.3. Let d and T >. Let N be as defined in (3.8 and let P be a smooth function such that γ. Then there exists a positive constant C depending only on ρ, P, λ, µ, κ and d such that N(b, n L1 t (Ḃ d 1 C( (b, n CL t + (b, n 3 CL t (3.1 for all (b, n CL T, t T, provided further that γ >. Furthermore, there exists a positive constant C depending only on ρ, P, λ, µ, κ and d such that N(b, n L1 t (Ḃ d 1 C((1 + t (b, n CL t + (b, n 3 CL t (3.11 for all (b, n CL T, t T, provided further that γ =. Proof. Taking q = 1,, in (3.9, we in particular have (γb, Λb, n L (, T ; Ḃ d 1 (Rd L (, T ; Ḃ d (R d L 1 (, T ; Ḃ d +1 (Rd (3.1 for all (b, n CL T. Note that Q satisfies Q 1 (b = (b+ρ and Q ( = 1 ρ thanks to Lemma 5.7 and Lemma 5.5, we have >. Thus, L (Q(bn L1 t (Ḃ d C (Q(bn + Q(bdiv 1 nid L1 t (Ḃ d C (Q (b Q (( bn L1 t (Ḃ d + Q (( bn L1 t (Ḃ d + Q(bdiv n L1 t (Ḃ d ( C b L t (Ḃ d b L t (Ḃ d n L t (Ḃ d + b L t (Ḃ d n L t (Ḃ d + b L t (Ḃ d div n L1 t (Ḃ d. (3.13 Now, a standard interpolation argument ensures that the right hand-side of (3.13 is bounded by Similarly, using (3.1, we have (b, n CL T + (b, n 3 CL T. κb b L1 t (Ḃ d 1 C b L t (Ḃ d b L1 t (Ḃ d 1 C b L t (Ḃ d b L1 t (Ḃ d ; ( div (n n C n ρ L1 t (Ḃ d 1 n L1 t (Ḃ d C n L ; t (Ḃ d (3.15 div (Q(bn n L1 t (Ḃ d C Q(bn 1 n L1 t (Ḃ d C b L t (Ḃ d n L. t (Ḃ d For the pressure term when γ >, (3.16 ( P (b + ρ P (ρ b L1 t (Ḃ d C P (b + ρ 1 P (ρ L t (Ḃ d 1 b L1 t (Ḃ d C b L t (Ḃ d 1 b L1 t (Ḃ d, (3.17

16 16 N. CHIKAMI AND T. KOBAYASHI where we have used the regularity for γb in (3.1. Thus, combining (3.13, (3.14, (3.15, (3.16, and (3.17, we have N(b, n L1 t (Ḃ d 1 ( b L t (Ḃ d + 1 b L t (Ḃ d n L t (Ḃ d ( + b L t (Ḃ d div n L1 t (Ḃ d + b L1 t (Ḃ d ( + n L t (Ḃ d 1 + b L t (Ḃ d provided that γ >. Thus, we have (3.1 if γ >. If, on the other hand, γ =, we only have + b L t (Ḃ d 1 b L1 t (Ḃ d, Λb L (, T ; Ḃ d 1 (Rd L (, T ; Ḃ d (R d L 1 (, T ; Ḃ d +1 (Rd. We may write the term P (b + ρ b = (P (b + ρ and (P (b + ρ L1 t (Ḃ d C P (b + ρ 1 L1 t (Ḃ d = C P (b + ρ P (ρ L1 t (Ḃ d C b L1 t (Ḃ d, thanks to Lemma 5.7. However, we have no control on b L1 t (Ḃ d, so the above only implies Alternatively, we have (P (b + ρ L1 t (Ḃ d 1 Ct b L t (Ḃ d. P (b + ρ b L1 t (Ḃ d C P (b + ρ 1 L1 t (Ḃ d b L t C b L1 t (Ḃ d b L t (Ḃ d 1 (Ḃ d 1 Ct b L t (Ḃ d b L t (Ḃ d Ct b L. 1 t (Ḃ d 1 (3.18 Either way, we may not close the estimates globally-in-time if γ = with our given regularity in CL T : Combining (3.13, (3.14, (3.15, (3.16, and (3.18, we have N(b, n L1 t (Ḃ d 1 ( b L t (Ḃ d + 1 b L t (Ḃ d n L t (Ḃ d ( + b L t (Ḃ d div n L1 t (Ḃ d + b L1 t (Ḃ d ( + n L t (Ḃ d 1 + b L t (Ḃ d + t b L t (Ḃ d in this case. Thus, we have (3.11 if γ =. This completes the proof of Lemma 3.3. Proof of Theorem 1.4. The time-continuity follows from the standard density argument as in the corresponding proof for the linear heat equation. The proof of Theorem 1.4 readily

17 WELL-POSEDNESS AND TIME-DECAY ESTIMATES OF CNSK 17 follows from Corollary 3., Lemma 3.3 and Lemma 3.1. As the proof is standard, we only sketch the steps. We let (a L, m L be the solution to the linear equation: t a L + div m L =, t m L 1 Lm L κρ a L + γ a L =, ρ (a, m t= = (a, m and let (ã, m := (a a L, m m L. Then Corollary 3. with s d 1 and σ 1 gives (a, m CLT (a L, m L CLT + (ã, m CLT C ( ((γ + Λa, m Ḃ d 1 + N(b, n L1 t (Ḃ d 1, for all (b, n CL T. Lemma 3.3 ensures that N(b, n can be regarded as a combination of bi-and-trilinear continuous operators in CL T. We define a ball in CL T centered at the origin by B CL T R ( := {(a, m CL T ; (a, m CLT R}, (3.19 where R >. Lemma 3.1 now ensures the existence of a unique solution in the ball B CL T R ( for a sufficiently small data. Note that T = is allowed for in the case γ >, thanks to the global estimate in Lemma 3.3. By (3.1, (a, m CLT C ((γ + Λa, m d Ḃ 1 is also ensured Global solution when γ =. From hereon, we focus on the case γ =. When γ >, the term γ ρ induces a damping effect, which provides the control of b L1 t (Ḃ d. The lack of this time-integrability in the case γ = is troublesome, especially when we control the pressure term. As shown by (3.18, we need the control of b L1 t (Ḃ d, whereas we only have b L t (Ḃ d in CL T. In order to compensate this lack of global control, we impose stronger regularity on the low-frequency, which is of lower order in terms of the scaling, to gain the necessary global control of the density. To handle this case, we have to change the expressions of the capillary and pressure terms. We linearize the capillary term while keeping it in the divergence form: div K(ρ = κρ a + div (K(a, where as before we denoted a := ρ ρ. Similarly, we write the pressure as follows: By Taylor expansion, we have P (a + ρ = P (ρ + P (ρ a + a P (a, where P : R + R is smooth and vanishes at. Using the assumption that P (ρ =, we have P (a + ρ = (a P (a. With this expression, we obtain the same linear structure with nonlinear terms all written in divergence form: N(b, n = L (Q(bn + div (K(b 1 ρ div (n n div (Q(bn n + (b P (b. (3. This expression enables us to close necessary product estimates in the Besov space of regularity index s = d 3 by a sole restriction d 3.

18 18 N. CHIKAMI AND T. KOBAYASHI Given initial data (a, m such that (Λa, m Ḃ d 1 (Rd, we show that the map Φ : (b, n (a, m defined in (3.6 with (a, m the solution to (3.7 with γ =, i.e., t a + div m =, t m 1 Lm κρ a = N(b, n, ρ (a, m t= = (a, m, where N = N(b, n is as defined in (3., has a fixed point in a ball of a suitable metric space. To define such a space, for d 3, we let CL T be given by CL T := L (, T ; Ḃ d 3 (Rd Ḃ d 1 (Rd L 1 (, T ; Ḃ d 1 (Rd Ḃ d +1 (Rd equipped with a norm (a, m CLT := (Λa, m L T (Ḃ d 3 Ḃ d 1 + Λ (Λa, m L1 T (Ḃ d 3 Ḃ d 1. By interpolation, we may readily see that (a, m CL T satisfies (3.9 and, additionally, (Λa, m L d q (, T ; Ḃ 3+ q (R d for q [1, ]. (3.1 Then we have the following global nonlinear estimates for the case γ =. Lemma 3.4. Let d 3 and T >. Let N be as defined in (3.8 and let P be a smooth function such that γ =. Then there exists a positive constant C depending only on ρ, P, λ, µ, κ and d such that N(b, n L1 t (Ḃ d 3 Ḃ d 1 C( (b, n CLt + (b, n 3 CLt (3. for all (b, n CL T and t T. Remark 3.5. The restriction d 3 for (3. is used to close product estimates in Ḃ d 3 (Rd. Proof. Let d 3. Note that by (3.1, we still have (Λb, n L (, T ; Ḃ d 1 (Rd L (, T ; Ḃ d (R d L 1 (, T ; Ḃ d +1 (Rd. From this, it is clear that we have (3.13, (3.14, (3.15 and (3.16 as before. Furthermore, (b, n CL T implies (3.1: (Λb, n L (, T ; Ḃ d 3 (Rd L (, T ; Ḃ d (Rd L 1 (, T ; Ḃ d 1 (Rd. (3.3 To handle the last term (P (b+ρ L1 t (Ḃ d 1, we use the above additional low-frequency assumption on b. We have P (b + ρ b L1 t (Ḃ d C P (b + ρ 1 L1 t (Ḃ d b L t C b L t (Ḃ d b L1 t (Ḃ d, (Ḃ d 1

19 WELL-POSEDNESS AND TIME-DECAY ESTIMATES OF CNSK 19 thanks to (3.1 and (3.3. Notice that, here, we do not need to use the expression (3.. Combining all information, we have N(b, n L1 t (Ḃ d 1 ( b L t (Ḃ d + 1 b L t (Ḃ d n L t (Ḃ d ( + b L t (Ḃ d div n L1 t (Ḃ d + b L1 t (Ḃ d + b L1 t (Ḃ d ( + n L 1 + t (Ḃ d b L t (Ḃ d. In order to complete the proof of (3., we still have to bound N(b, n L1 t (Ḃ d 3 under the assumption that d 3. Thanks to (3.1, (3.3, Lemma 5.7, Lemma 5.5, we have L (Q(bn L1 t (Ḃ d 3 C (Q(bn + Q(bdiv nid L1 t (Ḃ d C (Q (b Q (( bn L1 t (Ḃ d ( C b L t (Ḃ d b L t (Ḃ d n L t (Ḃ d + b L t (Ḃ d n L1 t (Ḃ d 1 + Q (( bn L1 t (Ḃ d + Q(bdiv n L1 t (Ḃ d ; + b L t (Ḃ d n L t (Ḃ d div (K(b L1 t (Ḃ d 3 C ( b b Id κ a a L1 t (Ḃ d 1 div (n n ρ L1 t (Ḃ d 3 b L1 t (Ḃ d + b Id κ b b L1 t (Ḃ d b L t (Ḃ d + b L t (Ḃ d b L t (Ḃ d Λb L + t (Ḃ d 1 Λb L t (Ḃ d Λb L t (Ḃ d ; (3.4 (3.5 C n n L1 t (Ḃ d C n L t (Ḃ d n L t (Ḃ d ; (3.6 div (Q(bn n L1 t (Ḃ d C Q(bn 3 n L1 t (Ḃ d C Q(b L t (Ḃ d n L t (Ḃ d n L t (Ḃ d C b L t (Ḃ d n L t (Ḃ d n L t (Ḃ d ; (3.7 (b P (b L1 t (Ḃ d b P 3 (b L1 t (Ḃ d C b L1 t (Ḃ d P (b L t (Ḃ d C Λb L t (Ḃ d 3 Λb L1 t (Ḃ d. 1 (3.8 Note here that if we use the expression (3.8, then instead of (3.5 and (3.8, we end up with using Lemma 5.5 with s = d 3, which requires d 4.

20 N. CHIKAMI AND T. KOBAYASHI Summing up (3.4, (3.5, (3.6, (3.7 and (3.8, we have N(b, n L1 t (Ḃ d 3 b L t (Ḃ d ( b L t (Ḃ d n L t (Ḃ d + n L t (Ḃ d n L t (Ḃ d + b L t (Ḃ d ( n L1 t (Ḃ d + 1 b L1 t (Ḃ d + n L t (Ḃ d ( b L t (Ḃ d + n L t (Ḃ d + b L t (Ḃ d 3 b L1 t (Ḃ d, which completes the proof of (3.. Proof of Theorem 1.6. Similarly to Theorem 1.4, the proof of Theorem 1.6 readily follows from Corollary 3., Lemma 3.4 and Lemma 3.1. The details are left to the readers. 4. Proof of decay estimates In this section, we prove the decay estimates for the global small solutions that have been obtained in the previous section. Proof of Theorem 1.8. The proof below is inspired by [1]. For notational simplicity, in the following, we denote X(t := (a, m CLT = ((γ + Λa, m L t (Ḃ d 1 + Λ ((γ + Λa, m L1 t (Ḃ d 1. and X := (γ + Λa, m d Ḃ 1. Firstly, note that for a solution to (CNSK constructed by Theorem 1.4 with D(t <, we also have τ α ((γ + Λa, m h L ((γ + Λa, m t (Ḃ d h L + τ α ((γ + Λa, m 1 p,1 t (Ḃ d h L 1 p,1 t (Ḃ d +1 p,1 X(t + D(t <. (4.1 Step 1: Bounds for the low frequencies. We denote the semigroup associated to the left hand-side of System (.5 by e tk and let [ ] [ ] [ ] (γ + Λa (γ + Λa U =, U m = and F = F (a, m =, m N(a, m where N = N(a, m is as defined in (3.. By Duhamel s formula, we have the expression U(t = e tk U + e (t τk F (a, m(τdτ. Since ((γ + Λa, m Ḃ d, (Rd, the action of the semigroup on Ḃs (Rd given by Lemma.4 implies In addition, we have e tk U l Ḃ s Ct 1 (s+ d ((γ + Λa, m Ḃ d., e tk U l Ḃ s j j js j e tk U L j j j(s+ d j d j U L C U l Ḃ d, provided that s + d >, where C depends only on s, d and j. Thus, we have e tk U l Ḃ s C t 1 (s+ d U l. (4. Ḃ d,

21 WELL-POSEDNESS AND TIME-DECAY ESTIMATES OF CNSK 1 By (4. and the integral representation, we have U(t l = ((γ + Λa, m(t l e tk U Ḃ s Ḃ s l + Ḃ s ( C t 1 (s+ d U l + Ḃ d, e (t τk N(a, m(τ l Ḃdτ s t τ 1 (s+ d N(a, m(τ l Ḃ d, (4.3 where the constant C depends on j (but stays bounded as long as < j. It now comes down to bounding the nonlinear terms with our prescribed regularity by X(t and D(t defined in (1.6. We suppose now that α in D(t is some real number to be decided later. In what follows, we are going to be frequently using the fact that Indeed, we have, thanks to (4.1, sup τ d 4 ((γ + Λa, m(τ Ḃ X(t + D(t. (4.4 τ t dτ, u(τ Ḃ u l (τ Ḃ + u h (τ Ḃ u l (τ Ḃ + u h (τ Ḃ d 1 τ d 4 τ d 4 ul (τ Ḃ + τ α τ α u h (τ Ḃ d 1 τ d 4 D(t + τ α (X(t + D(t, for u {(γ + Λa, m} if d 1, i.e., d. Inequality (4.4 follows if Estimate of L(Q(am: We have t τ 1 (s+ d L(Q(am(τ l dτ Ḃ d, α d 4. (4.5 t τ 1 (s+ d Q(a(τm(τ Ḃ d, dτ t τ 1 (s+ d a(τ Ḃ m(τ Ḃ dτ thanks to Lemma 5.6. Splitting the frequencies of m, we have, for the integral pertaining m l (τ Ḃ, t τ 1 (s+ d a(τ Ḃ m l (τ Ḃ dτ sup τ t ( τ d 4 a(τ Ḃ t τ 1 (s+ d τ (1+ d dτ sup τ t t 1 (s+ d (X(t + D(tD(t, ( τ 1 (+ d m(τ l Ḃ thanks to (4.4, Lemma 5.1 and 1 (s + d < 1 + d, i.e., s < + d, which is automatically satisfied since s d + 1.

22 N. CHIKAMI AND T. KOBAYASHI As for the high frequencies, we first let t. Then we have t τ τ t 1 for all τ t, which gives t τ 1 (s+ d a(τ Ḃ m h (τ Ḃ dτ τ d 4 mh (τ Ḃ dτ D(t d m(τ Ḃ +1dτ D(t t 1 (s+ d X(tD(t. We next treat the case t. We split the integrals into I 1 and I as follows: t τ 1 (s+ d a(τ Ḃ m h (τ Ḃ dτ ( 1 = + t τ 1 (s+ d a(τ Ḃ m h (τ Ḃ dτ =: I 1 + I. 1 Note that if τ 1 and t, we have and t t t τ. As such, we have I 1 1 t τ 1 (s+ d τ d 4 m(τ h t 1 (s+ d 1 t 1 (s+ d X(tD(t. Ḃ d +1 d m(τ Ḃ +1dτ D(t On the other hand, since τ α α τ α if 1 τ t, we have I 1 t τ 1 (s+ d a(τ Ḃ m(τ h Ḃ dτ 1 1 dτ D(t t τ 1 (s+ d τ d 4 α ( τ α m(τ h Ḃ d +1 dτ D(t t τ 1 (s+ d τ d 4 α dτ D(t t 1 (s+ d D(t (4.6 (4.7 by Lemma 5.1 and 1 (s + d d 4 + α, i.e., s α. (4.8 Thus, we obtain t τ 1 (s+ d L(Q(am(τ l dτ t 1 Ḃ d (s+ d (X(t + D(tD(t, (4.9, provided that α satisfies (4.5 and (4.8. Estimate of div (K(a: Since div (K(a = a a, we have t τ 1 (s+ d div K(a(τ l dτ Ḃ d, t τ 1 (s+ d a a Ḃ d dτ, t τ 1 (s+ d a(τ Ḃ1 a(τ Ḃ 1 dτ

23 WELL-POSEDNESS AND TIME-DECAY ESTIMATES OF CNSK 3 thanks to the limiting case of the product rule (Lemma 5.6. We split the frequencies of the higher order term as before: a(τ Ḃ a l (τ Ḃ + a h (τ Ḃ. Note that we have 1 (s + d d ( + d, i.e., s + d, which is satisfied by assumption, so Lemma 5.1 implies t τ 1 (s+ d a(τ Ḃ1 a l (τ Ḃ dτ t τ 1 (s+ d τ d 4 1 (+ d dτ (X(t + D(tD(t t 1 (s+ d (X(t + D(tD(t, thanks to (4.4 and the definition of D(t. The bounds pertaining to the term a h (τ Ḃ can be carried out similarly to the bounds for L(Q(am. Let t first. Since t τ τ t 1 for all τ t, we have t τ 1 (s+ d a(τ Ḃ1 a h (τ Ḃ dτ τ d 4 a(τ hḃ dτ D(t a(τ h dτ D(t Ḃ d + t 1 (s+ d X(tD(t by the definition of X(t and D(t. We next treat the case t. Splitting the timeintegral, we have We have t τ 1 (s+ d a(τ Ḃ1 a h (τ Ḃ dτ ( 1 = + t τ 1 (s+ d a(τ Ḃ1 a h (τ Ḃ dτ =: I 3 + I I 3 t 1 (s+ d d a(τ Ḃ +dτ D(t t 1 (s+ d X(tD(t by a similar calculation to (4.6. On the other hand, for I 4, we have I 4 t 1 (s+ d D(t similarly to (4.7, thanks to Lemma 5.1 and 1 (s + d d 4 + α, i.e., (4.8. Thus, we obtain t τ 1 (s+ d div K(a(τ l dτ t 1 Ḃ d (s+ d (X(t + D(tD(t, (4.1, provided that α satisfies (4.5 and (4.8. Estimate of 1 ρ div (m m: We have t τ 1 (s+ d div (m m(τ l dτ Ḃ d, t τ 1 (s+ d (m m(τ Ḃ1 d, dτ t τ 1 (s+ d m(τ Ḃ m(τ Ḃ1 dτ

24 4 N. CHIKAMI AND T. KOBAYASHI by Lemma 5.6. Splitting the frequencies of the higher order term by m(τ Ḃ1 m l (τ Ḃ1 + m h (τ Ḃ1, we may carry out the estimates similarly to previous computations. Lemma 5.1 gives t τ 1 (s+ d m(τ Ḃ m l (τ Ḃ1 dτ t τ 1 (s+ d τ d 4 ml (τ Ḃ1 dτ D(t t τ 1 (s+ d τ d 4 1 (1+ d dτ (X(t + D(tD(t t 1 (s+ d (X(t + D(tD(t, thanks to (4.4 and the definition of D(t. For the high frequencies m h (τ Ḃ1, let t, so that we have t τ 1 (s+ d m(τ Ḃ m h (τ Ḃ1 dτ τ d 4 m(τ hḃ1 dτ D(t m(τ h dτ D(t Ḃ d +1 t 1 (s+ d X(tD(t by the definition of X(t and D(t. We next treat the case t. Splitting the timeintegral, we have t τ 1 (s+ d m(τ Ḃ m h (τ Ḃ1 dτ ( 1 = + t τ 1 (s+ d m(τ Ḃ m h (τ Ḃ1 dτ =: I 5 + I 6. 1 We have I 5 t 1 (s+ d 1 d m(τ Ḃ +1dτ D(t t 1 (s+ d X(tD(t by a similar calculation to (4.6. For I 6, we have I 6 t 1 (s+ d D(t similarly to (4.7, thanks to Lemma 5.1 and 1 (s + d d 4 + α, i.e., (4.8. Thus, we have in the end t τ 1 (s+ d div (m m(τ l Ḃ d, provided that α satisfies (4.5 and (4.8. dτ t 1 (s+ d (X(t + D(tD(t, (4.11

25 WELL-POSEDNESS AND TIME-DECAY ESTIMATES OF CNSK 5 Estimate of div (Q(am m: Lemma 5.6 and the composition estimate (Lemma 5.7 give t τ 1 (s+ d div (Q(am m (τ l dτ Ḃ d At this point, we notice that a(τ Ḃ d X(t and, t τ 1 (s+ d (Q(am m(τ Ḃ1 d dτ, t τ 1 (s+ d a(τ Ḃ d (m m(τ Ḃ1 d dτ. m m Ḃ1 d = ij m i m j Ḃ1 d ij (m i m j Ḃ d by Bernstein inequality, so the situation coincides with that of div (m m. By the same computation as that for (4.11, we readily obtain t τ 1 (s+ d div (Q(am m(τ l Ḃ d, provided that α satisfies (4.5 and (4.8. Estimate of (P (a P ( a: Lemma 5.6 and Lemma 5.7 give t τ 1 (s+ d (P (a(τ P ( a(τ l dτ Ḃ d dτ t 1 (s+ d X(t(X(t + D(tD(t,, t τ 1 (s+ d P (a(τ P ( Ḃ a(τ Ḃ1 dτ t τ 1 (s+ d a(τ Ḃ a(τ Ḃ1 dτ. By the same computation as that for (4.9, we obtain t τ 1 (s+ d (P (a(τ P ( a(τ l Ḃ d, provided that α satisfies (4.5 and (4.8. (4.1 dτ t 1 (s+ d (X(t + D(tD(t, (4.13 In the end, inserting (4.9, (4.1, (4.11, (4.1 and (4.13 into (4.3, we deduce that ( sup τ t τ 1 (s+ d U(τ l C Ḃ s U l + X(t(X(t + D(tD(t Ḃ d,, (4.14 given that (4.5 and (4.8 hold, where the constant C depends on d, µ, λ, κ and γ and the frequency-threshold j can be arbitrary.

26 6 N. CHIKAMI AND T. KOBAYASHI Step : Bounds for the high frequencies. Rewriting System (CNSK in terms of the weighted unknowns (t α a, t α m, we have t (t α a + div (t α m = αt α 1 a, t (t α m 1 L(t α m κρ (t α a + γ (t α a = αt α 1 m + t α F (a, m, (4.15 ρ (t α a, t α m t= = (,, Using Lemma. with a t α a, m t α m, f αt α 1 a, g αt α 1 v+t α F (a, m, s d +1 and q r, we have τ α ((γ+λa, m h L (γ + Λ(ατ α 1 a + ατ α 1 m + τ α F (a, m t (Ḃ d h L, +1 t (Ḃ d 1 where the constant appearing on the right hand-side does not depend on j. We first handle the lower order linear terms on the right hand-side of the above: For u {(γ + Λa, m}, we have when τ t and (γ + Λ(ατ α 1 a + ατ α 1 m h L t ατ α 1 u h L t (Ḃ d 1. (Ḃ d 1 u L t (Ḃ d 1 X(t ατ α 1 u h L (,1;Ḃ d 1 u L t (Ḃ d 1 X(t when t and τ 1. Moreover, when t and 1 τ t, by the frequency cut-off of the norm, we have ατ α 1 u h L = α j( d (1,t;Ḃ d 1 +1 j τ α 1 u L (1,t;L j j α j d +1 τ α u L (,t;l j j j( α j τ α u h L t (Ḃ d +1. Taking j Z, such that α j < 1 4, for instance, we may absorb the contributions of these terms into the left hand-side. As α only depends on the dimension and the parameter ε, so does j. Thus, we arrive at τ α ((γ+λa, m h L t (Ḃ d +1 X(t + τ α F (a, m h L t (Ḃ d 1, (4.16 where the frequency-threshold j depends only on the dimension d and the given constants such as µ, λ, κ and γ. It now comes down to estimating the nonlinear terms. Estimate of L(Q(am: τ α L(Q(am h L t ij (Ḃ d 1 τ α i j (Q(am h L t ij (Ḃ d 1 τ α ( i j (Q(am + i (Q(a j m + Q(a i j m h L t (Ḃ d 1.

27 WELL-POSEDNESS AND TIME-DECAY ESTIMATES OF CNSK 7 We have τ α i j (Q(am L t (Ḃ d τ α 1 i j (Q(a L t (Ḃ d m L t (Ḃ d 1 τ α a h L t (Ḃ d + m L t (Ḃ d X(t τ α 1 a L t (Ḃ d ; + τ α i (Q(a j m L t (Ḃ d 1 i (Q(a L t (Ḃ d τ α 1 j m L t (Ḃ d a L t (Ḃ d τ α m L t (Ḃ d X(t τ α +1 m L t (Ḃ d ; +1 τ α Q(a i j m L t Note that we have (Ḃ d 1 Q(a L t (Ḃ d τ α i j m L t a L t (Ḃ d 1 (Ḃ d 1 τ α m L t (Ḃ d +1 X(t τ α m L t (Ḃ d +1. τ α a L t (Ḃ d τ α a l + τ α a + L t (Ḃ d h L D(t, + ε t (Ḃ d + provided that α 1 ( d + ε + d, i.e., We also have α 1 (d + ε = d + 1 ε. τ α m L t (Ḃ d τ α m l + τ α m +1 L t (Ḃ d h L D(t, +1 ε t (Ḃ d +1 provided that α 1 ( d + 1 ε + d, i.e., Thus, we have provided that α satisfies (4.17. τ α L(Q(am h L t Estimate of div (K(a: Since div (K(a = a a, we have α 1 (d + 1 ε. (4.17 (Ḃ d 1 X(tD(t, (4.18 τ α div (K(a h L t (Ḃ d 1 = τ α a a h L t (Ḃ d 1 a L t (Ḃ d τ α a L t (Ḃ d X(t τ α a L t (Ḃ d. We have τ α a L t (Ḃ d τ α a l + τ α a L t (Ḃ d h L, + ε t (Ḃ d provided that α 1 ( d + ε + d, i.e., α 1 (d + ε, which is satisfied if (4.17 holds. Thus τ α div (K(a h L X(tD(t (4.19 t (Ḃ d 1 holds, provided that α satisfies (4.17.

28 8 N. CHIKAMI AND T. KOBAYASHI Estimate of 1 ρ div (m m: We have τ α div (m m h L t (Ḃ d 1 = τ α (m m + mdiv m h L t (Ḃ d 1 m L t (Ḃ d τ α 1 m L t (Ḃ d X(t τ α m L t (Ḃ d. As for the decay functional, we deduce τ α m L t (Ḃ d τ α m l + τ α m L t (Ḃ d h L, +1 ε t (Ḃ d +1 α 1 ( d + 1 ε + d, i.e., (4.17. Thus, we have provided that (4.17 holds. τ α div (m m h L t Estimate of div (Q(am m: We have (Ḃ d 1 X(tD(t, (4. τ α div (Q(am m h L t τ α Q(am m (Ḃ d h L 1 t (Ḃ d a L t (Ḃ d τ α m m L t (Ḃ d X(t j τ α j (m m L t (Ḃ d. By (4., we obtain provided that (4.17 holds. τ α div (Q(am m h L t Estimate of (P (a P ( a: We have (Ḃ d 1 X(tD(t, (4.1 τ α (P (a P ( a h L t (Ḃ d 1 a L t (Ḃ d 1 τ α a L t (Ḃ d X(t τ α a L t (Ḃ d, in which we deduce for the decay functional if (4.17 holds. Thus, this implies τ α a L t (Ḃ d τ α a l + τ α a L t (Ḃ d h L, +1 ε t (Ḃ d +1 provided that (4.17 holds. τ α div (Q(am m h L t (Ḃ d 1 X(tD(t, (4.

29 WELL-POSEDNESS AND TIME-DECAY ESTIMATES OF CNSK 9 Inserting (4.18, (4.19, (4., (4.1 and (4. into (4.16, we deduce that there exists some frequency-threshold j and a constant C depending only on d, µ, λ, κ and γ such that the following inequality holds: τ α ((γ+λa, m h L t (Ḃ d +1 C(X(t + (X(t + D(tD(t. (4.3 Conclusion. By (4.5, (4.8 and (4.17, we now know that α must satisfy { d max 4, s } α 1 (d + 1 ε for some arbitrary small ε >. Thus, the choice of α = 1 (d + 1 ε is justified under our assumption on s. Putting together (4.14 and (4.3, we finally arrive at D(t D + X(t + X(t + D(t. As Theorem 1.4 ensures that X(t X with X being small, we may now conclude that (1.7 is fulfilled for all time provided that D and X are small enough. This concludes Theorem 1.8. Proof of Theorem 1.9. As much of the argument is similar to the proof of Theorem 1.8, we only sketch the proof here. As before, for notational simplicity, we denote X(t := (a, m CLT = (Λa, m L t (Ḃ d + Λ (Λa, 1 m L1 t (Ḃ d. 1 and X := (Λa, m d Ḃ 1. Note that, for any solution to (CNSK with γ = constructed by Theorem 1.6, if in addition D(t <, we have τ α (Λa, m h L t In addition, similarly to (4.4, we have (Ḃ d 3 Ḃ d 1 p,1 (Λa, m h L t (Ḃ d 3 Ḃ d 1 p,1 + τ α (Λa, m h L t (Ḃ d +1 p,1 Starting from (4.3, our task is to bound X(t + D(t <. sup τ d 4 (Λa, m(τ Ḃ X(t + D(t. τ t t τ 1 (s+ d N(a, m(τ l dτ, Ḃ d, provided that (4.5 holds. Identical calculations leading to (4.9, (4.1, (4.11, (4.1 and (4.13 yield t τ 1 (s+ d L(Q(am(τ l Ḃ d, t τ 1 (s+ d div K(a(τ l Ḃ d, (4.4 dτ t 1 (s+ d ( X(t + D(t D(t, (4.5 dτ t 1 (s+ d ( X(t + D(t D(t, (4.6

30 3 N. CHIKAMI AND T. KOBAYASHI and t τ 1 (s+ d div (m m(τ l Ḃ d, t τ 1 (s+ d div (Q(am m(τ l Ḃ d, t τ 1 (s+ d (a P (a(τ l Ḃ d, dτ t 1 (s+ d ( X(t + D(t D(t, (4.7 dτ t 1 (s+ d X(t( X(t + D(t D(t (4.8 dτ t 1 (s+ d ( X(t + D(t, (4.9 respectively, provided that α satisfies (4.5 and (4.8. In the end, by (4.5, (4.6, (4.7, (4.8 and (4.9, we obtain ( sup τ t τ 1 (s+ d U(τ l C Ḃ s U l + X(t( X(t + Ḃ d, D(t D(t, (4.3 where the constant C depends on d, µ, λ and κ and the frequency-threshold j can be arbitrary. The bounds for the high frequencies are also carried out in the same way as the proof for Theorem 1.8. Applying Lemma. to (4.15 with γ and employing the same frequency cut-off argument leading to (4.16, we obtain τ α (Λa, m h L X(t + τ α F (a, m t (Ḃ d h L, +1 t (Ḃ d 1 where the frequency-threshold j depends only on the dimension d and the given constants such as µ, λ, κ. We begin estimating the nonlinear terms in the high-frequency regime. This too is carried out in the same way as before. By the same computations leading to (4.18, (4.19, (4., (4.1 we may easily see τ α L(Q(am h L t τ α div (K(a h L t τ α div (m m h L t (Ḃ d 1 X(t D(t, (4.31 (Ḃ d 1 X(t D(t, (4.3 (Ḃ d 1 X(t D(t, (4.33 τ α div (Q(am m h L t (Ḃ d 1 X(t D(t, (4.34 respectively, provided that α satisfies (4.17. Lastly, to estimate the pressure term, we have to slightly modify the calculation of (4.. We have τ α (P (a P ( a h L t (Ḃ d 1 a L t (Ḃ d τ α a L t X(t τ α a L t X(t( X(t + D(t, (Ḃ d 1 (Ḃ d 1 (4.35