arxiv: v1 [math.ap] 6 Feb 2013

Transcription

1 AN ASYMPTOTIC LIMIT OF A NAVIER-STOKES SYSTEM WITH CAPILLARY EFFECTS ANSGAR JÜNGEL, CHI-KUN LIN, AND KUNG-CHIEN WU arxiv:13.199v1 [math.ap] 6 Feb 13 Abstract. A combined incompressible and vanishing capillarity limit in the barotropic compressible Navier-Stokes equations for smooth solutions is proved. The equations are considered on the two-dimensional torus with well prepared initial data. The momentum equation contains a rotational term originating from a Coriolis force, a general Kortewegtype tensor modeling capillary effects, and a density-dependent viscosity. The limiting model is the viscous quasi-geostrophic equation for the rotated velocity potential. The proof of the singular limit is based on the modulated energy method with a careful choice of the correction terms. 1. Introduction The aim of this paper is to prove a combined incompressible and vanishing capillarity limit for a two-dimensional Navier-Stokes-Korteweg system, leading to the viscous quasigeostrophic equation. We consider the (dimensionless) mass and momentum equations for the particle density ρ(x,t) and the mean velocity u(x,t) = (u 1 (x,t),u (x,t)) of a fluid in the two-dimensional torus : (1) () with initial conditions t ρ+div(ρu) = in, t >, t (ρu)+div(ρu u)+ρu + p(ρ) = div(k +S), ρ(,) = ρ, u(,) = u in. Date: February 7, 13. Mathematics Subject Classification. 35B5, 35Q35, 35Q4. Key words and phrases. Incompressible limit, vanishing capillarity limit, modulated energy, Navier- Stokes-Korteweg model, quasi-geostrophic equation. The first author acknowledges partial support from the Austrian Science Fund (FWF), grants P18, P434, and I395 and from the National Science Council of Taiwan (NSC) through the Tsungming Tu Award. The second author was supported by the NSC, grant NSC M-9-8-MY. The third author acknowledges partial support from the Tsz-Tza Foundation of the Institute of Mathematics, Academia Sinica (Taipei, Taiwan). He thanks Clément Mouhot for his kind invitation to visit the University of Cambridge during the academic years Part of this work was written during the stay of the first author at the Center of Mathematical Modeling and Scientific Computing, National Chiao Tung University (Hsinchu, Taiwan), and at the Institute of Mathematics, Academia Sinica (Taipei, Taiwan); he thanks Chi-Kun Lin and Tai-Ping Liu for their kind hospitality. 1

2 A. JÜNGEL, C.-K. LIN, AND K.-C. WU Here, ρu describes the Coriolis force, u = ( u,u 1 ), the function p(ρ) = ρ γ /γ with γ > 1 denotes the pressure of an ideal gas obeying Boyle s law, K is the Korteweg-type tension tensor and S the viscous stress tensor. More precisely, the free surface tension tensor is given by divk = κ ρ (σ (ρ) σ(ρ)), where κ >, which can be written in conservative form as (3) divk = κ div (( S(ρ) 1 ) ) S (ρ) ρ I σ(ρ) σ(ρ), where S (ρ) = ρσ (ρ), σ(ρ) is a (nonlinear) function, and I denotes the unit matrix in R. For a general introduction and the physical background of Navier-Stokes-Korteweg systems, we refer to [7, 11, 19]. In standard Korteweg models, κ(ρ) = σ (ρ) defines the capillarity coefficient [11, Formula (1.9)]. In the shallow-water equation, often σ(ρ) = ρ is used such that divk = ρ ρ (see, e.g., [5, 5]). Bresch and Desjardins [6] employed general functions σ(ρ) and suitable viscosities allowing for additional energy estimates (also see []). If σ(ρ) = ρ, the third-order term can be interpreted as a quantum correction, and system (1)-() (without the rotational term) corresponds to the so-called quantum Navier-Stokes model, derived in [8] and analyzed in [19]. The viscous stress tensor is defined by divs = div(µ(ρ)d(u)), where D(u) = 1 ( u + u ) and µ(ρ) denotes the density-dependent viscosity. Often, the viscosity in the Navier-Stokes model is assumed to be constant for the mathematical analysis [13]. Density-dependent viscosities of the form µ(ρ) = ρ were chosen in [5] and were derived, in the context of the quantum Navier-Stokes model, in [8]. The choice µ(ρ) = σ(ρ) allows one to exploit a certain entropy structure of the system [6]. Without capillary effects, system (1)-() reduces to the viscous shallow-water or viscous Saint-Venant equations, whose inviscid version was introduced in [7]. The viscous model was formally derived from the three-dimensional Navier-Stokes equations with a free moving boundary condition [14]. This derivation was generalized later to varying river topologies [5]. The existence of global weak or strong solutions to the Kortewegtype shallow-water equations was proved in [6, 7, 15, 16, 18] under various assumptions on the nonlinear functions. In [7], the authors obtained several existence results of weak solutions under various assumptions concerning the density dependency of the coefficients. The notion of weak solution involves test functions depending on the density; this allows one to circumvent the vacuum problem. Duan et al. [1] showed the existence of local classical solutions to the shallow-water model without capillary effects. For more details and references on the shallow-water system, we refer to the review [4]. The combined incompressible and vanishing capillarity limit studied in this work is based on the scaling t εt, u εu, µ(ρ) εµ(ρ) and on the choice κ = ε α ( < α,ε < 1), which gives (4) t ρ ε +div(ρ ε u ε ) = in, t >,

3 AN ASYMPTOTIC LIMIT OF A NAVIER-STOKES SYSTEM 3 (5) t (ρ ε u ε )+div(ρ ε u ε u ε )+ 1 ε ρ εu ε + 1 ε γ (ργ ε ) ε(α 1) ρ ε (σ (ρ ε ) σ(ρ ε )) = div(µ(ρ ε )D(u ε )), with the initial conditions (6) ρ ε (,) = ρ ε, u ε(,) = u ε in. The condition α < 1 is needed to control the capillary energy; see the energy identity in Lemma below. The local existence of smooth solutions to(4)-(5) is discussed in Appendix A. When letting ε, it holds ρ ε 1 and ρ ε u ε φ = ( φ/ x, φ/ x 1 ) in appropriate function spaces, where φ solves the viscous quasi-geostrophic equation [6, Chapter 6] (see Section for details) (7) (8) t ( φ φ)+( φ )( φ) = µ(1) φ in, t >, φ(,) = φ in. The objective of this paper is to make this limit rigorous. Our proof requires the (local) existence of a smooth solution to(7)-(8), which is shown in Appendix A. Several derivations of inviscid quasi-geostrophic equations have been published; see, e.g., [9, 1, 8]. The reader is also referred to the monograph [4] for a more complete discussion of this model. The viscous equation was derived rigorously for weak solutions from the shallow-water system in [5]. The proof is essentially based on the presence of the additional viscous part div(ρ u) and a friction term in the momentum equation. The novelty of the present paper is that these expressions are not needed and that more general expressions can be considered. In particular, we allow for viscous terms of the type div(µ(ρ)d(u)), and no friction is prescribed. In the following, we describe our main result. In order to simplify the presentation, we assume that the nonlinearities are given by power-law functions: σ(ρ) = ρ s, µ(ρ) = ρ m for ρ, where s > and m >. The exponents s and m cannot be chosen freely; we need to suppose that (9) < s 1, m = s+ 1 γ +1. This assumption includes the quantum Navier-Stokes model s = 1/, m = 1 and the shallow-water model with s = 1, m = 3/. Furthermore, we assume that the initial data are sufficiently regular (ensuring the local-in-time existence of smooth solutions) ρ ε H k ( ), u ε H k 1 ( ), φ H k+1 ( ), where k >, and that they are well prepared: (1) G ε (φ ε ) φ, ε 1 (ρ ε 1) φ, ρ ε u ε φ, ε α 1 ρ ε

4 4 A. JÜNGEL, C.-K. LIN, AND K.-C. WU in L ( ) as ε, where ρ ε = 1+εφ ε (this defines φ ε), (11) G ε (φ ε ) = ε sign(φ ε) h(1+εφ ε ), ρ ε = 1+εφ ε, and the internal energy h(ρ) is defined by h (ρ) = p (ρ)/ρ = ρ γ and h(1) = h (1) = (see (13) for an explicit expression). Note that the convergence ε 1 (ρ ε 1) φ in L ( ) implies that G ε (φ ε) φ in L 1 ( ) if ρ ε is bounded in L ( ) (see (17)). Theorem 1. Let < α < 1 and γ > 1. We suppose that (9) holds and that the initial data satisfy (1). Furthermore, let (ρ ε,u ε ) be the classical solution to (4)-(6) and let φ be the classical solution to (7)-(8), both on the time interval (,T). Then, as ε, Furthermore, if s < 1 for all 1 p < and 1 q <. ρ ε 1 in L (,T;L γ ( )), ρ ε u ε φ in L (,T;L γ/(γ+1) ( )). and γ (1 s) or if s = 1 and γ, ρ ε 1 in L (,T;L p ( )), ρ ε u ε φ in L (,T;L q ( )), The proof is based on the modulated energy method, first introduced by Brenier in a kinetic context [] and later extended to various models, e.g. [1, 3, ]. The idea of the method is to estimate, through its time derivative, a suitable modification of the energy by introducing in the energy the solution of the limit equation. We suggest the following form of the modulated energy: ( ρε H ε (t) = u ε φ + 1 ) G ε(φ ε ) φ +ε (α 1) σ(ρ ε ) dx t (1) + µ(ρ ε ) D(u ε ) D( φ) dx, These terms express the differences of the kinetic, internal, and Korteweg energies as well as the viscosity. Differentiating the modulated energy with respect to time and employing the evolution equations, elaborated computations lead to the inequality H ε (t) C t H ε (s)ds+o(1), t >, where o(1) denotes terms vanishing in the limit ε, uniformly in time. The Gronwall lemma then implies the result. The paper is organized as follows. In Section, we derive the energy identities for the shallow-water system and the quasi-geostrophic equation and give a formal derivation of the latter model from the former one. Theorem 1 is proved in Section 3. In the appendix, we discuss the existence of local smooth solutions to (4)-(5) and give an existence proof for local smooth solutions to (7)-(8).

5 AN ASYMPTOTIC LIMIT OF A NAVIER-STOKES SYSTEM 5. Auxiliary results In this section, we derive the energy estimates for (4)-(5) and derive formally the quasigeostrophic equation (7). Based on the definition h (ρ) = p (ρ)/ρ, h(1) = h (1) =, we can give an explicit formula for this function: 1 ( (13) h(ρ) = ρ γ 1 γ(ρ 1) ), ρ. γ(γ 1) The energy identity for (4)-(5) is given as follows. Lemma. Let (ρ ε,u ε ) be a smooth solution to (4)-(6) on (,T). Then the energy identity de ε dt +D ε =, t (,T), holds, where the energy E ε and energy dissipation D ε are defined by, respectively, ( 1 E ε = ε h(ρ ε)+ 1 ) ρ ε u ε +ε (α 1) σ(ρ ε ) dx, D ε = µ(ρ ε ) D(u ε ) dx. Proof. Multiply (4) by ε h (ρ ε ) 1 u ε ε (α 1) σ (ρ ε ) σ(ρ ε ), integrate over, and then integrate by parts: ( 1 = T ε th(ρ ε ) 1 (ρ ε h ε ) ρ ε (ρ ε u ε ) 1 u ε t ρ ε +ρ ε u ε u ε u ε ) +4ε (α 1) σ(ρ ε ) t σ(ρ ε ) ε (α 1) div(ρ ε u ε )σ (ρ ε ) σ(ρ ε ) dx. Multiplying (5) by u ε and integrating over gives, since u ε u ε =, ( = t (ρu ε ) u ε ρ ε (u ε u ε ) : u ε + 1 ε ργ 1 ε ρ ε u ε +ε (α 1) σ (ρ ε ) σ(ρ ε )div(ρ ε u ε ) µ(ρ ε )D(u ε ) : u ε )dx, where : means summation over both matrix indices. Observing that h satisfies h (ρ ε ) = ρε γ and using the identity D(u ε ) : u ε = D(u ε ), the sum of the above two equations becomes ( d 1 dt T ε h(ρ ε)+ 1 ) ρ ε u ε +ε (α 1) σ(ρ ε ) dx+ µ(ρ ε ) D(u ε ) dx =, which proves the lemma. A consequence of the energy identity is the following estimate. Lemma 3. Let (ρ ε,u ε ) be a smooth solution to (4)-(6) on (,T). Then there exists C > such that for all < ε < 1, (14) (15) ρ ε 1 L (,T;L γ ( )) Cε min{1,/γ} if γ > 1, ρ ε 1 L (,T;L ( )) Cε if γ.

6 6 A. JÜNGEL, C.-K. LIN, AND K.-C. WU Proof. If γ =, h(ρ) = 1 (ρ 1), and the result follows immediately from Lemma. Let γ >. We claim that h(ρ) ρ 1 γ /(γ(γ 1)) for ρ. Then the result follows again from the energy identity. Indeed, the function f(ρ) = ρ γ 1 γ(ρ 1) ρ 1 γ is convex in ( 1, ) and concave in (, 1 ). Since the values f() = γ and f(1 ) = γ/ 1 are positive, f on [, 1]. Furthermore, f(1) = f (1) = which implies, together with the convexity, that f in [ 1, ), proving the claim. Finally, let γ <. By [3, p. 591], h(ρ) c R ρ 1 for ρ R and h(ρ) c R ρ 1 γ for ρ > R, for some c R > and R >. Hence, using Hölder s inequality and γ <, ( γ/ ρ ε 1 γ L γ ( ) C ρ ε 1 dx) + ρ ε 1 γ dx ( C {ρ ε R} {ρ ε R} γ/ h(ρ ε )dx) +C {ρ ε>r} {ρ ε>r} h(ρ ε )dx C(ε γ +ε ) Cε γ, where here and in the following C > denotes a generic constant not depending on ε. Estimate (15) for γ follows from ρ ε 1 L ( ) T = (ρ ε 1) dx C h(ρ ε )dx Cε, which finishes the proof. We perform the formal limit ε in (4)-(5). For this, we observe that (4) can be written in terms of φ ε = (ρ ε 1)/ε as follows: t φ ε +div(φ ε u ε )+ 1 ε divu ε =. We apply the operator div (defined by div (v 1,v ) = v 1 / x + v / x 1 ) to (5) and observe that div (ρ ε u ε )/ε = divu ε/ε+div(φ ε u ε ) = t φ ε, by the above equation. Then we find that (16) t div (ρ ε u ε )+div div(ρ ε u ε u ε ) t φ ε = ε (α 1) div ( ρ ε (σ (ρ ε ) σ(ρ ε )) ) +div div(µ(ρ ε )D(u ε )). By the energy estimate, ρ ε 1 (in L (,T;L γ ( ))). Assuming that φ ε φ and u ε φ in suitable function spaces and employing the relations div div( φ φ) = ( φ )( φ), div div(d( φ)) = φ, theformallimit in(16)yields thelimit equation(7). Theinitial condition readsasφ(,) = φ, where φ = lim ε φ ε (,) in. The energy and the energy dissipation of (7) equal E = 1 ( φ +φ )dx, D = µ(1) D( φ) dx. Multiplying the limiting equation by φ and using the properties ( φ )( φ)φdx =, ( φ) dx = D( φ) dx,

7 AN ASYMPTOTIC LIMIT OF A NAVIER-STOKES SYSTEM 7 we find the energy identity of the viscous quasi-geostrophic equation: de +D =, t >. dt First, we prove the following lemma. 3. Proof of Theorem 1 Lemma 4. Let T >, γ > 1, and < α < 1. Then lim H ε(t) = uniformly in (,T), ε where H ε is defined in (1). Proof. Using the definitions of the energy and energy dissipation as well as the relation 1 G ε(φ ε ) = ε h(ρ ε ), we write t H ε (t) = (E ε +E)(t)+ (D ε +D)(s)ds+ 1 (ρ ε 1) φ dx T (G ε (φ ε ) φ ε )φdx ρ ε u ε φdx φ ε φdx T t t + (µ(ρ ε ) µ(1)) D( φ) dxds 4 µ(ρ ε )D(u ε ) : D( φ)dxds = I 1 + +I 8. The aim is to estimate dh ε /dt. To this end, we treat the integrals I j or their derivatives d term by term. By the energy estimates, (I dt 1 + I ) =. The integral I 3 cancels with a contribution originating from I 5 ; see below. The estimate of I 4,...,I 8 (or their derivatives) is performed in several steps. Step 1: estimate of I 4. L Hôpital s rule shows that for γ > 1, h(1+z) lim = 1 ( z z, lim 1 h(1+z) 1 ) = γ. z z z 6 Therefore, there exists a nonnegative function f, defined on [, ), such that h(1+z) = 1 z f(z) for z, and a function g, defined on [, ), such that f(z) 1 = zg(z) for z. Furthermore, the inequalities f(z) f() = 1 and g(z) C(1 + z (γ 3)+ ) hold, where z + = max{,z}. Finally, we claim that f(z) = h(1+z)/z (1+z) γ /(γ(γ 1)) for z and γ 4. Indeed, the function w(z) = h(1+z) z (1+z) γ /(γ(γ 1)) is convex in [, ) and w() = w () =, which implies that w(z) in [, ), proving the claim. With these preparations, we can estimate the difference G ε (φ ε ) φ ε appearing in I 4 : ( ) G ε (φ ε ) φ ε = sign(φ ε) h(1+εφε ) φ ε ε = φ ε f(εφ ε ) 1 = φ ε f(εφ ε ) 1 f(εφε )+1 = φ ε εφ ε g(εφ ε ) f(εφε )+1.

8 8 A. JÜNGEL, C.-K. LIN, AND K.-C. WU In view of the bounds for f and g as well as the relation εφ ε = ρ ε 1, we infer that (17) G ε (φ ε ) φ ε C ε ρ ε 1 1+ρ(γ 3)+ ε f(εφε )+1. This bound allows us to estimate I 4. Indeed, if 1 < γ < 4, by (14), I 4 (t) C ε φ L (,T;L ( )) ρ ε 1 L (,T;L 1 ( )) Cεmin{1,/γ} 1 = o(1) uniformly in (,T). Here and in the following, the constant C > depends on φ and its derivatives but not on ε. If γ 4, we have, using the upper bound of f(z) for γ 4, (17), and 1+εφ ε = ρ ε, G ε (φ ε ) φ ε C ε ρ ε 1 1+ρ γ 3 ε Cρ (γ )/ ε +1 C ε ρ ε 1 ( 1+ρ ε (γ 3) (γ )/ We employ estimates (14)-(15) and Hölder s inequality to conclude that I 4 (t) C φ L (,T;L ( ))ε 1 ρ ε 1 L (,T;L ( )) ρ ε 1 L (,T;L γ ( )) ( 1+ ρ ε (γ 4)/ ) L (,T;L γ ( )) Cε /γ φ L (,T;L ( )) = o(1). Step : estimate of di 5 /dt. Inserting the momentum equation (5) and integrating by parts, it follows that di 5 dt = t (ρ ε u ε ) φdx ρ ε u ε t φdx T = ρ ε (u ε u ε ) : φdx+ 1 ρ ε u ε φdx T ε + 1 ρ γ ε ε φdxd ε (α 1) ρ ε ( σ (ρ ε ) σ(ρ ε ) ) φdx γ T + µ(ρ ε )D(u ε ) : φdx ρ ε u ε t φdx = J 1 + +J 6. We treat the integrals J 1,...,J 6 term by term. The integral J can be written as J = 1 ρ ε u ε φdx. ε The third integral vanishes since div = : J 3 = 1 ρ γ ε ε γ div( φ)dx =. Using the identity (3) and div =, we compute J 4 = ε (( S(ρ (α 1) ε ) 1 ) S (ρ ε ) ρ ε div( φ) ).

9 AN ASYMPTOTIC LIMIT OF A NAVIER-STOKES SYSTEM 9 ) ( σ(ρ ε ) σ(ρ ε )) : φ dx CH ε. Integration by parts and using div = again yields J 5 = µ(ρ ε )u ε ( φ+ div( φ))dx T µ (ρ ε )( ρ ε u ε +u ε ρ ε ) : φdx T = µ(ρ ε )u ε φdx T µ (ρ ε ) σ(ρε ) ( ρε σ (ρ ε )( ρ ε u ε )+( ρ ε u ε ) σ(ρ ε ) ) : φdx. Theassumptionsonµandσ (see(9))yieldµ (ρ ε )/( ρ ε σ (ρ ε )) = ρ m s 1/ ε. Hence, applying the Cauchy-Schwarz inequality, the last integral is bounded from above by C σ(ρ ε ) L ( ) ρ ε u ε L ( ) Cε (1 α) = o(1). We conclude that J 5 µ(ρ ε )u ε φdx+o(1). The integral J 6 remains unchanged. Finally, we estimate J 1. To this end, we add and substract the expression φ such that J 1 = K 1 + +K 4, where K 1 = ρ ε (u ε φ) (u ε φ) : φdx, T K = ρ ε φ u ε : φdx, T K 3 = ρ ε u ε φ : φdx, T K 4 = ρ ε φ φ : φdx. The first integral can be bounded by the modulated energy: K 1 C ρ ε u ε φ dx CH ε. A reformulation yields K = ρ ε u ε (( φ ) φ ) dx. We employ the continuity equation (4) to find K 3 = 1 ρ ε u ε φ dxd = 1 div(ρ ε u ε ) φ dx

10 1 A. JÜNGEL, C.-K. LIN, AND K.-C. WU = 1 t (ρ ε 1) φ dx = 1 d (ρ ε 1) φ dx+ 1 (ρ ε 1) t φ dx dt = di 3 dt +o(1). Finally, using again div =, ( K 4 = ρ ε ( φ ) φ ) φdx T = (ρ ε 1) ( ( φ ) φ ) ( φdx ( φ ) φ ) φdx T = (ρ ε 1) ( ( φ ) φ ) φdx 1 φ ( φ )dx T = (ρ ε 1) ( ( φ ) φ ) φdx+ 1 div( φ) φ dx = o(1). In the last step, we have employed estimate (14) for ρ ε 1. Summarizing the estimates for K 1,...,K 4, we have shown that J 1 CH ε di 3 dt ( ( φ ) φ ) (ρ ε u ε )dx+o(1). Then, summarizing the estimates for J 1,...,J 6, we obtain di 5 dt CH ε di 3 dt + 1 ρ ε u ε φdx ε T ( ( t + φ ) φ+µ(1) φ ) (ρ ε u ε )dx T ( ) µ(ρε ) µ(1)ρ ε uε φdx+o(1). The last integral can be estimated by employing the assumptions on µ and Hölder s inequality: µ(ρ ε ) µ(1)ρ ε ρε u ε φdx C ρ m 1/ ε ρ 1/ ε L ( ) ρ ε u ε L ( ). ρε We claim that the first factor on the right-hand side is of order o(1). To prove this statement, we consider first 1 < m < 1: ρ m 1/ ε ρ 1/ ε L ( ) ρ m 1 ε ρ ε 1 (1 m) dx ρ ε m 1 L γ ( ) ρ ε 1 (1 m) L p ( ), where p = γ(1 m)/(γ m+1). The inequality p γ is equivalent to γ 1. Note that the Hölder inequality can be applied since we supposed that m 1 γ; see (9). Second,

11 AN ASYMPTOTIC LIMIT OF A NAVIER-STOKES SYSTEM 11 let 1 < m (the case m = 1 being trivial). We compute ρ m 1/ ε ρε 1/ L ( ) ρ ε ρ ε 1 (m 1) dx ρ ε L γ ( ) ρ ε 1 (m 1) L q ( ), where q = γ(m 1)/(γ 1), and q γ if and only if m (γ + 1)/. Finally, if m (γ +1)/, we find that ρ m 1/ ε ρ 1/ ε L ( ) C ρ ε (1+ρ m ε ) ρ ε 1 dx C(1+ ρ ε m 3 L γ ( ) ) ρ ε 1 L r ( ), with r = γ/(γ m+3) satisfying r γ if and only if m (γ+1)/. We conclude that µ(ρ ε ) µ(1)ρ ε ρε u ε φdx C ρ ε 1 β L γ ( ) ρε for some β >, and together with (14), this shows that the integral is of order o(1). Therefore, di 5 dt CH ε di 3 dt + 1 ρ ε u ε φdx ε T ( (18) ( t + φ ) φ+µ(1) φ ) (ρ ε u ε )dx+o(1). Step 3: estimate of di 6 /dt. Employing (4) and (7), we can write di 6 dt = t φ ε φdx φ ε t φdx = 1 ( div(ρ ε u ε )φdx ( t + φ )( φ) µ(1) φ ) φ ε dx ε = 1 ( (19) ρ ε u ε φdx+ ( t + φ )( φ) µ(1) φ ) φ ε dx. ε We observe that the first integral on the right-hand side cancels with the corresponding integral in(18). To deal with the second integral, we employ again the momentum equation (5). We write 1 γ ργ ε = (γ 1) h(ρ ε)+ (ρ ε 1) = (γ 1) h(ρ ε )+ε φ ε. Then, because of (u ε ) = u ε, (5) is equivalent to where φ ε = ρ ε u ε εf ε, F ε = t (ρ ε u ε )+div(ρ ε u ε u ε )+ γ 1 ε h(ρ ε ) div(µ(ρ ε )D(u ε )) ε (α 1) ( S(ρ ε ) 1 (S (ρ ε ) ρ ε ) div ( σ(ρ ε ) σ(ρ ε ) )).

12 1 A. JÜNGEL, C.-K. LIN, AND K.-C. WU Replacing φ ε in the second integral in (19) by the above expression gives ( ( t + φ )( φ) µ(1) φ ) φ ε dx ( = ( t + φ )( φ) µ(1) φ ) (ρ ε u ε εfε )dx. We claim that the integral containing Fε is bounded in an appropriate space. Indeed, let ψ be a smooth (vector-valued) test function. The first term of F ε is written in weak form as follows: T T t (ρ ε u ε ) ψdxds = ρ ε u ε t ψdxds+ (ρ ε u ε )(t) ψ(t)dx T ρ ε u ε ψ()dx. These integrals are bounded if ρ ε u ε is bounded in L (,T;L 1 ( )). This is the case, since mass conservation and the energy estimate show that ρ ε u ε dx 1 ρ ε dx+ 1 ρ ε u ε dx is uniformly bounded in (,T). An integration by parts gives T T div(ρ ε u ε u ε ) ψdxds = ρ ε u ε u ε : ψdxds, and this integral is uniformly bounded, by the energy estimate. Furthermore, again integrating by parts, T ( ) γ 1 h(ρ T ε ε ) div(µ(ρ ε )D(u ε )) ψdxds T ( ) γ 1 = h(ρ ε ε )I µ(ρ ε )D(u ε ) : ψdxds, which is uniformly bounded since we can estimate T µ(ρ ε )D(u ε ) dxds 1 T µ(ρ ε )dxds+ 1 T µ(ρ ε ) D(u ε ) dxds and µ(ρ ε ) C(1+ρ γ ε ). Also the remaining terms are bounded since T ( ε (α 1) (S(ρ ε ) S(1)) 1 ) (S (ρ ε ) ρ ε ) div( σ(ρ ε ) σ(ρ ε )) ψdxds T = ε (α 1) ( (S(ρε ) S(1)) divψ + 1 S (ρ ε ) ρ ε divψ ( σ(ρ ε ) σ(ρ ε )) : ψ ) dxds.

13 AN ASYMPTOTIC LIMIT OF A NAVIER-STOKES SYSTEM 13 Using the Hölder continuity of S(z) = (s/)z s, z, the first summand can be estimated by C ρ ε 1 min{1,s}. Weinfer that thecorresponding integral isof order o(1). We formulate the second summand as 1 t ε(α 1) (s 1) σ(ρ ε ) divψdxds. In view of the energy estimate, this integral as well as the third summand are uniformly bounded. This shows that ( ( t + φ )( φ) µ(1) φ ) φ ε dx ( = ( t + φ )( φ) µ(1) φ ) (ρ ε u ε )dx+o(1), and consequently, (19) becomes di 6 dt = 1 ρ ε u ε φdx ε T ( + ( t + φ )( φ) µ(1) φ ) (ρ ε u ε )dx+o(1). Step 4: estimate of di 7 /dt. The function µ satisfies µ(z) µ(1) = z m 1 z 1 m if m 1 and µ(z) µ(1) C(1+z m 1 ) z 1 if m > 1, for z. Therefore, if m 1, taking into account (14), di 7 dt ρ ε 1 m L (,T;L γ ( )) D( φ) L (,T;L γ/(γ m) ( )) Cεmmin{1,/γ}. Moreover, if 1 < m (γ +1)/, using Hölder s inequality, di 7 dt C( ) 1+ ρ ε L (,T;L (m 1)γ/(γ 1) ( )) ρε 1 L (,T;L γ ( )) Cε min{1,/γ}. The norm of ρ ε is uniformly bounded since (m 1)γ/(γ 1) γ is equivalent to m γ. Step 5: estimate of di 8 /dt. Integration by parts yields di 8 dt = µ (ρ ε ) ρ ε u ε : φdx+ µ(ρ ε )u ε φdx T µ (ρ ε ) = ρε σ (ρ ε ) σ(ρ ε) ( ρ ε u ε ) : φdx + (µ(ρ ε ) µ(1)ρ ε )u ε φdx+µ(1) ρ ε u ε φdx. By definition of µ and σ (see (9)), it follows that di 8 dt C σ(ρ ε) L ( ) ρ ε u ε L ( ) +C ρ m 1/ +µ(1) ρ ε u ε φdx. ε ρ 1/ ε L ( ) ρ ε u ε L ( )

14 14 A. JÜNGEL, C.-K. LIN, AND K.-C. WU Because of the energy estimate, the first summand is of order o(1). The second summand has been estimated in Step, and it has been found that it is also of order o(1). This shows that di 8 dt µ(1) ρ ε u ε ψdx+o(1). Step 6: conclusion. Adding the estimates for di 4 /dt,...,di 8 /dt, most of the integrals cancel, and we end up with Integrating over (, t) gives dh ε dt H ε (t) H ε ()+C CH ε + di 4 dt +o(1). t H ε (s)ds+i 4 (t) I 4 ()+o(1). By Step 1, I 4 (t) = o(1). Furthermore, I 4 () = o(1) by assumption. It holds that H ε () = o(1) since ρ ε(u ε φ ) L ( ) ρ εu ε φ L ( ) + (1 ρ ε) φ L ( ) ρ εu ε φ L ( ) + 1 ρ ε L ( ) φ L ( ) = o(1) andsince theinitialdataarewell prepared. ThentheGronwall lemma implies thath ε (t) = o(1) finishing the proof. We are now in the position to prove Theorem 1 which is a consequence of Lemma 4. We observe that by (14), ρ ε 1 in L (,T;L γ ( )) and, using the Hölder inequality and γ/(γ +1) < γ, ρ ε u ε φ L (,T;L γ/(γ+1) ( )) ρ ε L (,T;L γ ( )) ρ ε (u ε φ) L (,T;L ( )) () + ρ ε 1 L (,T;L γ/(γ+1) ( )) φ L (,T;L ( )) C ρ ε (u ε φ) L (,T;L ( )) +C ρ ε 1 L (,T;L γ ( )). We conclude that ρ ε u ε φ in L (,T;L γ/(γ+1) ( )). Next, let γ (1 s) and < s < 1/. Because of assumption (9), i.e. γ s, we have γ/(γ +(1 s)) γ, and hence, ρ ε 1 in L (,T;L γ/(γ+(1 s)) ( )) as ε. Furthermore, since α < 1, σ(ρ ε ) in L (,T;L ( )) as ε and thus, by Hölder s inequality, (1) (ρ ε 1) L (,T;L γ/(γ+(1 s)) ( )) = σ (ρ ε ) 1 σ(ρ ε ) L (,T;L γ/(γ+(1 s)) ( )) ρ ε 1 s L (,T;L γ ( )) σ(ρ ε) L (,T;L ( )).

15 AN ASYMPTOTIC LIMIT OF A NAVIER-STOKES SYSTEM 15 We infer that ρ ε 1 in L (,T;W 1,γ/(γ+(1 s)) ( )). Because of the continuous embeddingw 1,γ/(γ+(1 s)) ( ) L γ/(1 s) ( ),thisimpliesthatρ ε 1inL (,T;L γ/(1 s) ( )). Since γ/(γ + (1 s) ) γ/(1 s), this gives ρ ε 1 in L (,T;L γ/(γ+(1 s)) ( )). Applying the same procedure as in (1) again, we obtain (ρ ε 1) L (,T;L γ/(γ+(1 s)) ( )) ρ ε 1 s L (,T;L γ/(1 s) ( )) σ(ρ ε) L (,T;L ( )). Hence, ρ ε 1 in L (,T;W 1,γ/(γ+(1 s)) ( )) and in L (,T;L γ/(1 s) ( )). Repeating this argument, we conclude that ρ ε 1 in L (,T;L p ( )) for all p <. For the momentum, we obtain for p 1 ρ ε u ε φ L (,T;L p/(p+1) ( )) ρ ε L (,T;L p ( )) ρ ε (u ε φ) L (,T;L ( )) + ρ ε 1 L (,T;L p/(p+1) ( )) φ L (,T;L ( )) C ρ ε (u ε φ) L (,T;L ( )) +C ρ ε 1 L (,T;L p ( )). This shows that ρ ε u ε φ in L (,T;L q ( )) for all q <. Finally, let γ and s = 1. Then ρ ε 1 in L (,T;H 1 ( )) and, by the continuous embedding H 1 ( ) L p ( ) for all p <, also ρ ε 1 in L (,T;L p ( )) for all p <. The theorem is proved. Appendix A. Local existence of smooth solutions The local existence of smooth solutions to the Navier-Stokes-Korteweg system (4)-(5) can be shown similarly as in [1]. We only sketch the proof since it is highly technical and does not involve new ideas. First, we rewrite (4)-(5), setting ρ = ρ ε, u = u ε, and ε = 1. Taking the divergence of (5) and replacing div t (ρu) by (4), which has been differentiated with respect to time, we obtain tt ρ 1 γ ργ +ρσ (ρ) ρ = divdiv(ρu u) div(ρu ) +divdiv(µ(ρ)d(u))+f[ρ], where F[ρ] = div(ρ (σ (ρ) σ(ρ))) ρσ (ρ) ρ involves only three derivatives. This formulation allows one to treat the momentum equation as a nonlinear fourth-order wave equation for which existence and regularity results can be applied. In order to derive some regularity for the velocity, Li and Marcati [1] assumed that curlu =. Then u is reconstructed from the problem divv = 1 ρ ( tρ+ ρ u), curlv =, v(t)dx = ū(t). Theorem.1 in [1] gives the existence of a unique solution u H s+1 ( ) to this problem, provided that the right-hand side satisfies ( t ρ+ ρ u)/ρ H s ( ). Actually, Li and Marcati replace the right-hand side by ( t ρ+ ρ u)/ψ, where ψ solves the mass equation t ψ +ψdivv +u ρ =, t >, ψ() = ρ.

16 16 A. JÜNGEL, C.-K. LIN, AND K.-C. WU The reason is that this equation can be solved explicitly, yielding strictly positive solutions ψ. The existence proof is based on an iteration scheme: Given (ρ p,ψ p,u p,v p ), solve divv p+1 = f p (t), curlv p+1 =, v p+1 (t)dx = ū(t), t ψ p+1 +ψ p+1 divv p +u p ρ p =, t >, ψ() = ρ, tt ρ p+1 1 γ ργ p+1 +ψ pσ (ψ p ) ρ p+1 = g p (t), t >, ρ p+1 () = ρ, t ρ p+1 () = ρ divu ρ u, t u p+1 +u p+1 = h p (t), where f p (t), g p (t), and h p (t) contain the remaining terms (see [1, Section 3] for details). The existence of solutions to these linear problems follows from ODE theory and the theory of wave equations. The main effort is now to derive uniform estimates in Sobolev spaces H k ( ). This is done by multiplying the above equations by suitable test functions and assuming that T > is sufficiently small. By compactness, there exists a subsequence of (ρ p,ψ p,u p,v p ) which converges in a suitable Sobolev space to (ρ,ψ,u,v) as p. This limit allows us also to show that ρ = ψ and u = v. This shows the existence of local smooth solutions under the assumption of irrotational flow curl u =. Next, we prove the existence of local smooth solutions to the quasi-geostrophic equation (7). We set µ := µ(1) >. Theorem 5 (Local existence for the quasi-geostrophic equation). Let φ C ( ). Then there exists T > and a smooth solution φ to (7)-(8) for t T. Proof. The idea of the proof is to apply the theory of linear semigroups. Let p > and let A p : W,p ( ) R, A p (u) = µ u + u. Then A p is a sectorial operator satisfying R(λ) = 1 for all λ σ(a p ), where σ(a p ) denotes the spectrum of A p. Consequently, A p possesses the fractional powers A β p for β, defined on the domain X β,p = D(A β p). This space, endowed with its graph norm, satisfies X β,p W k,q ( ) if k /q < β /p, q p [17, Theorem 1.6.1]. Let max{1 1/p,1/+1/(p)} < β < 1 and set X := X β,p. The operator A p generates an analytical semigroup e tap (t ) [17, Theorem 1.3.4], and the following estimates hold for all t > [17, Theorem 1.4.3]: A p e tap u L p ( ) Ct β e δt u L p ( ), (e tap I)v L p ( ) Ct β A p v L p ( ) Ct β v X for < δ < 1, u L p ( ), and v X. Next, we reformulate (7). Set u = φ φ. Then (7) can be written as a system of two second-order equations: () (3) φ+φ = u in, t >, t u µ u+u = ( φ )(φ u)+µ(u φ)+u.

17 AN ASYMPTOTIC LIMIT OF A NAVIER-STOKES SYSTEM 17 We employ a fixed-point argument. Let T > and R >. We introduce the spaces Y = C ([,T];X) and B R = {u Y : u u Y R}, where u = φ +φ C ( ). Given u Y C ([,T];L p ( )), let φ L (,T;W,p ( )) be the unique solution to () satisfying the elliptic estimate φ W,p ( ) C u L p ( ). Then define J(u) = e tap u + t e (t s)ap F(φ(s),u(s))ds, where F(φ,u) = ( φ )(φ u)+µ(u φ)+u. Using the continuous embedding W,p ( ) W 1,p ( ) and the elliptic estimate for φ, we infer the estimate ) F(φ,u) L (,T;L p ( )) C u L (,T;W 1,p (T ))( 1+ u L (,T;W 1,p ( )) ) C u L (,T;X)( 1+ u L (,T;X) = C u Y (1+ u Y ). The last inequality follows from the embedding X W 1,p ( ) which holds for β > 1/+1/(p). We show that J maps B R into B R and that J : B R B R is a contraction for sufficiently small T >. Let T > be such that (e tap I)u L p ( ) CT β u X R/. Then, for u B R, J(u) u Y sup (e tap I)u L p ( ) <t<t + sup <t<t t R + sup <t<t A p e tap F(φ(s),u(s)) X ds t (t s) β e δ(t s) F(φ(s),u(s)) X ds R + CT1 β 1 β u Y(1+ u Y ) R, if T > is sufficiently small, using that u B R. Thus J(u) B R. In a similar way, we show that, for given u, v B R, J(u) J(v) Y CT1 β 1 β ( u Y + v Y ) u v Y. Again, choosing T > small enough, J becomes a contraction, and the fixed-point theorem of Banach provides the existence and uniqueness of a mild solution on [,T]. It remains to prove that the mild solution is smooth. Since β > 1 1/p, we have X W,p/ ( ) and hence u L (,T;W,p/ ( )) L (,T;W 1,p (T d )). Furthermore, φ L (,T;W 1,p ( )) L (,T;L ( )) (here, we use p > ). This shows that t u + A p (u) L (,T;L p ( )). Parabolic theory implies that u L q (,T;W,p ( )) for all q <. This improves the regularity of φ to φ L q (,T;W 4,p ( )). Hence, t u+a p (u) L q (,T;L ( )), and a bootstrap procedure finishes the proof.

18 18 A. JÜNGEL, C.-K. LIN, AND K.-C. WU References [1] M. Bostan. The Vlasov-Maxwell system with strong initial magnetic field: guiding-center approximation. Multiscale Model. Simul. 6 (7), [] Y. Brenier. Convergence of the Vlasov-Poisson system to the incompressible Euler equations. Commun. Part. Diff. Eqs. 5 (), [3] Y. Brenier, R. Natalini, and M. Puel. On a relaxation approximation of the incompressible Navier- Stokes equations. Proc. Amer. Math. Soc. 13 (4), [4] D. Bresch. Shallow-water equations and related topics. In: C. Dafermos et al. (eds.), Handbook of Differential Equations: Evolutionary Equations, Vol. V, pp Elsevier, Amsterdam, 9. [5] D. Bresch and B. Desjardins. Existence of global weak solutions for a D viscous shallow water equations and convergence to the quasi-geostrophic model. Commun. Math. Phys. 38 (3), [6] D. Bresch and B. Desjardins. Some diffusive capillary models of Korteweg type. C. R. Acad. Sci. Paris, Sec. Mécanique 33 (4), [7] D. Bresch, B. Desjardins, and C.-K. Lin. On some compressible fluid models: Korteweg, lubrication, and shallow water systems. Commun. Part. Diff. Eqs. 8, (3), [8] S. Brull and F. Méhats. Derivation of viscous correction terms for the isothermal quantum Euler model. Z. Angew. Math. Mech. 9 (1), [9] B. Desjardins and E. Grenier. Derivation of quasi-geostrophic potential vorticity equations. Adv. Diff. Eqs. 3 (1998), [1] B. Duan, L. Luo, and Y. Zheng. Local existence of classical solutions to shallow water equations with Cauchy data containing vacuum. SIAM J. Math. Anal. 44 (1), [11] J. Dunn and J. Serrin. On the thermomechanics of interstitial working. Archive Ration. Mech. Anal. 88 (1985), [1] P. Embid and A. Majda. Low Froude number limiting dynamics for stably stratified flow with small or finite Rossby numbers. Geosphys. Astrophys. Fluid Dyn. 87 (1998), 1-5. [13] E. Feireisl. Dynamics of Viscous Compressible Fluids. Oxford University Press, Oxford, 4. [14] J.-F. Gerbeau and B. Perthame. Derivation of viscous Saint-Venant system for laminar shallow water; numerical validation. Discrete Contin. Dyn. Sys. B 1 (1), [15] B. Haspot. Existence of global strong solutions for the shallow-water equations with large initial data. Preprint, [16] B. Haspot. Existence of global strong solutions for the Saint-Venant system with large initial data on the irrotational part of the velocity. C. R. Acad. Sci. Paris Math. 35 (1), [17] D. Henry. Geometric Theory of Semilinear Parabolic Equations. Springer, Berlin, [18] L. Hsiao and H.-L. Li. Dissipation and dispersion approximation to hydrodynamical equations and asymptotic limit. J. Part. Diff. Eqs. 1 (8), [19] A. Jüngel. Global weak solutions to compressible Navier-Stokes equations for quantum fluids. SIAM J. Math. Anal. 4 (1), [] A. Jüngel. Effective velocity in Navier-Stokes equations with third-order derivatives. Nonlin. Anal. 74 (11), [1] H.-L. Li and P. Marcati. Existence and asymptotic behavior of multi-dimensional quantum hydrodynamic model for semiconductors. Commun. Math. Phys. 45 (4), [] C.-K. Lin and K.-C. Wu. Hydrodynamic limits of the nonlinear Klein-Gordon equation. J. Math. Pures Appl. 98 (1), [3] P.-L. Lions and N. Masmoudi. Incompressible limit for a viscous compressible fluid. J. Math. Pure Appl. 77 (1998), [4] A. Majda. Introduction to PDEs and Waves for the Atmosphere and Ocean. Courant Lecture Notes in Mathematics, Vol. 9. Amer. Math. Soc., Providence, 3.

19 AN ASYMPTOTIC LIMIT OF A NAVIER-STOKES SYSTEM 19 [5] F. Marche. Derivation of a new two-dimensional viscous shallow water model with verying topology, bottom friction and capillary effects. Europ. J. Mech. B/Fluids 6 (7), [6] J. Pedlosky. Geophysical Fluid Dynamics. Springer, New York, 199. [7] A. de Saint-Venant. Théorie du mouvement non-permanent des eaux, avec application aux crues des rivières et à l introduction des marées dans leur lit. C. R. Acad. Sci. Paris 73 (1871), [8] S. Schochet. Singular limits in bounded domains for quasilinear symmetric hyperbolic systems having a vorticity equation. J. Differ. Eqs. 68 (1987), Institute for Analysis and Scientific Computing, Vienna University of Technology, Wiedner Hauptstraße 8 1, 14 Wien, Austria address: juengel@tuwien.ac.at Department of Applied Mathematics and Center of Mathematical Modeling and Scientific Computing, National Chiao Tung University, Hsinchu 31, Taiwan address: cklin@math.nctu.edu.tw Department of Pure Mathematics and Mathematical Statistics, University of Cambridge, Wilberforce Road, Cambridge CB3 WB, United Kingdom address: kcw8@dpmms.cam.ac.uk