INSTITUTE of MATHEMATICS. ACADEMY of SCIENCES of the CZECH REPUBLIC
|
|
- Anis Allison
- 5 years ago
- Views:
Transcription
1 INSTITUTEofMATHEMATICS Academy of Sciences Czech Republic INSTITUTE of MATHEMATICS ACADEMY of SCIENCES of the CZECH REPUBLIC Error estimates for a numerical method for the compressible Navier-Stokes system on sufficiently smooth domains Eduard Feireisl Radim Hošek David Maltese Antonín Novotný Preprint No PRAHA 2015
2
3 Error estimates for a numerical method for the compressible Navier-Stokes system on sufficiently smooth domains Eduard Feireisl Radim Hošek David Maltese Antonín Novotný Institute of Mathematics of the Academy of Sciences of the Czech Republic Žitná 25, CZ Praha 1, Czech Republic Institut Mathématiques de Toulon, EA2134, University of Toulon BP 20132, La Garde, France Abstract We derive an a priori error estimate for the numerical solution obtained by time and space discretization by the finite volume/finite element method of the barotropic Navier Stokes equations. The numerical solution on a convenient polyhedral domain approximating a sufficiently smooth bounded domain is compared with an exact solution of the barotropic Navier Stokes equations with a bounded density. The result is unconditional in the sense that there are no assumed bounds on the numerical solution. It is obtained by the combination of discrete relative energy inequality derived in [17] and several recent results in the theory of compressible Navier-Stokes equations concerning blow up criterion established in [26] and weak strong uniqueness principle established in [10]. ey words: Navier-Stokes system, finite element numerical method, finite volume numerical method, error estimates AMS classification 35Q30, 65N12, 65N30, 76N10, 76N15, 76M10, 76M12 1 Introduction We consider the compressible Navier-Stokes equations in the barotropic regime in a space-time cylinder Q T = 0, T Ω, where T > 0 is arbitrarily large and Ω R 3 is a bounded domain: t ϱ + div x ϱu = 0, 1.1 t ϱu + div x ϱu u + x pϱ = div x S x u, 1.2 In equations ϱ = ϱt, x 0 and u = ut, x R 3, t [0, T, x Ω are unknown density and velocity fields, while S and p are viscous stress and pressure characterizing the fluid via the constitutive relations S x u = µ x u + t xu 2 3 div xui, µ > 0, 1.3 The research of E.F. leading to these results has received funding from the European Research Council under the European Union s Seventh Framework Programme FP7/ / ERC Grant Agreement The Institute of Mathematics of the Academy of Sciences of the Czech Republic is supported by RVO: The research of R.H. leading to these results has received funding from the European Research Council under the European Union s Seventh Framework Programme FP7/ / ERC Grant Agreement The Institute of Mathematics of the Academy of Sciences of the Czech Republic is supported by RVO: The work of D.M. and A.N. has been supported by the MODTERCOM project within the APEX programme of the Provence-Alpes-Côte d Azur region 1
4 p C 2 0, C 1 [0,, p0 = 0, p ϱ > 0 for all ϱ 0, p ϱ lim ϱ ϱ γ 1 = p > 0, 1.4 where γ 1. Assumption p 0 > 0 excludes constitutive laws behaving as ϱ γ as ϱ 0 +. Error estimate stated in Theorem 3.8 however still holds in the case lim ϱ 0 + p ϱ > 0 at the price of some additional difficulties, ϱ γ 1 see [17] for more details. Equations are completed with the no-slip boundary conditions and initial conditions We notice that under assumption 1.3, we may write u Ω = 0, 1.5 ϱ0, = ϱ 0, u0, = u 0, ϱ 0 > 0 in Ω. 1.6 div x S x u = µ u + µ 3 xdiv x u. 1.7 The results on error estimates for numerical schemes for the compressible Navier-Stokes equations are in the mathematical literature on short supply. We refer the reader to papers of Liu [24], [25], Yovanovic [31], Gallouet et al. [17]. In [17] the authors have developed a methodology of deriving unconditional error estimates for the numerical schemes to the compressible Navier-Stokes equations and applied it to the numerical scheme discretizing the system on polyhedral domains. They have obtained error estimates for the discrete solution with respect to a classical solution of the system on the same polyhedral domain. In spite of the fact that [17] provides the first and to the best of our knowledge so far the sole error estimate for discrete solutions of a finite volume/finite element approximation to a model of compressible fluids that does not need any assumed bounds on the numerical solution itself, it has two weak points: 1 The existence of classical solutions on at least a short time interval to the compressible Navier-Stokes equations is known for smooth C 3 domains see Valli, Zajaczkowski [28] or Cho, Choe, im [2] but may not be in general true on the polyhedral domains. 2 The numerical solutions are compared with the classical exact solutions as is usual in any previous existing mathematical literature. In this paper we address both points raised above and to a certain extent remove the limitations of the theory presented in [17]. More precisely, we generalize the result of Gallouet et al. [17, Theorem 3.1] in two directions: 1 The physical domain Ω filled by the fluid and the numerical domain Ω h, h > 0 approximating the physical domain do not need to coincide. 2 If the physical domain is sufficiently smooth at least of class C 3 and the C 3 initial data satisfy natural compatibility conditions, we are able to obtain the unconditional error estimates with respect to any weak exact solution with bounded density. As in [17], and in contrast with any other error estimate literature dealing with finite volume or mixed finite volume/finite element methods for compressible fluids Yovanovich [31], Cancès et al [5], Eymard et al. [9], Villa, Villedieu [30], Rohde, Yovanovich [29], Gastaldo et al. [18] and others this result does not require any assumed bounds on the discrete solution: the sole bounds needed for the result are those provided by the numerical scheme. Moreover, in contrast with [17] and with all above mentioned papers, the exact solution is solely weak solution with bounded density. This seemingly weak hypothesis is compensated by the regularity and compatibility conditions imposed on initial data that make possible a sophisticated bootstrapping argument showing that weak solutions with bounded density are in fact strong solutions in the class investigated in [17]. These results are achieved by using the following tools: 2
5 1 The technique introduced in [17] modified in order to accommodate non-zero velocity of the exact sample solution on the boundary of the numerical domain. 2 Three fundamental recent results from the theory of compressible Navier-Stokes equations, namely Local in time existence of strong solutions in class by Cho, Choe, im [2]. Weak strong uniqueness principle proved in [10] see also [14]. Blow up criterion for strong solutions in the class by Sun, Wang, Zhang [26]. The three above mentioned items allow to show that the weak solution with bounded density emanating from the sufficiently smooth initial data is in fact a strong solution defined on the large time interval [0, T. 3 Bootstrapping argument using recent results on maximal regularity for parabolic systems by Danchin [8], Denk, Pruess, Hieber [3] and rylov [21]. The last item allows to bootstrap the strong solution in the class Cho, Choe, im [2] to the class needed for the error estimates in [17], provided a certain compatibility condition for the initial data is satisfied. 2 Preliminaries 2.1 Weak and strong solutions to the Navier-Stokes system We introduce the notion of the weak solution to system : Definition 2.1 Weak solutions. Let ϱ 0 : Ω [0, + and u 0 : Ω R 3 with finite energy E 0 = Ω 1 2 ϱ 0 u Hϱ 0 dx and finite mass 0 M 0 = Ω ϱ 0 dx. We shall say that the pair ϱ, u is a weak solution to the problem emanating from the initial data ϱ 0, u 0 if: a ϱ C weak [0, T ]; L a Ω, for a certain a > 1, ϱ 0 a.e. in 0, T, and u L 2 0, T ; W 1,2 0 Ω; R 3. b the continuity equation 1.1 is satisfied in the following weak sense ϱϕ dx τ τ = ϱ t ϕ + ϱu x ϕ dx dt, τ [0, T ], ϕ Cc [0, T ] Ω Ω 0 Ω c ϱu C weak [0, T ]; L b Ω; R 3, for a certain b > 1, and the momentum equation 1.2 is satisfied in the weak sense, ϱu ϕ dx τ τ = Ω 0 0 τ 0 Ω Ω ϱu t ϕ + ϱu u : ϕ + pϱ div ϕ dx dt µ u : x ϕ dx dt +µ + λdivudivϕ dx dt, τ [0, T ], ϕ Cc [0, T ] Ω; R 3. d The following energy inequality is satisfied 1 +Hϱ 2 ϱ u 2 dx τ τ + µ u 2 +µ+λ div u 2 dx dt 0, for a.a. τ 0, T, Ω Here and hereafter the symbol Ω 0 Ω 2.2 ϱ pz with Hϱ = ϱ dz z2 g dx τ 0 is meant for gτ, x dx g 0 x dx. 3 Ω Ω
6 In the above definition, we tacitly assume that all the integrals in the formulas are defined and we recall that C weak [0, T ]; L a Ω is the space of functions of L [0, T ]; L a Ω which are continuous as functions of time in the weak topology of the space L a Ω. We notice that the function ϱ Hϱ is a solution of the ordinary differential equation ϱh ϱ Hϱ = pϱ with the constant of integration fixed such that H1 = 0. Note that the existence of weak solutions emanating from the finite energy initial data is well-known on bounded Lipschitz domains provided γ > 3/2, see Lions [23] for large values of γ, Feireisl and coauthors [13] for γ > 3/2. Proposition 2.1. Suppose the Ω R 3 is a bounded domain of class C 3. Let r, V be a weak solution to problem in 0, T Ω, originating from the initial data r 0 C 3 Ω, r 0 > 0 in Ω, 2.5 satisfying the compatibility conditions V 0 C 3 Ω; R 3, 2.6 and such that V 0 Ω = 0, x pr 0 Ω = div x S x V 0 Ω, 2.7 Then r, V is a classical solution satisfying the bounds: 0 r r a.a. in 0, T Ω /r C[0,T ] Ω + r C 1 [0,T ] Ω + t x r C[0,T ];L 6 Ω;R t,tr C[0,T ];L 6 Ω D, 2.9 V C 1 [0,T ] Ω;R 3 + V C[0,T ];C 2 Ω;R 3 + t x V C[0,T ];L 6 Ω;R t,tv L 2 0,T ;L 6 Ω D, 2.10 where D depends on Ω, T, r, and the initial data r 0, V 0 via r 0, V 0 C 3 Ω;R 4 and min x Ω r 0x. Proof: The proof will be carried over in several steps. Step 1 According to Cho, Choe, and im [2], problem admits a strong solution unique in the class r C[0, T M ; W 1,6 Ω, t r C[0, T M ; L 6 Ω, 1/r L Q T, 2.11 V C[0, T M ]; W 2,2 Ω; R 3 L 2 0, T M ; W 2,6 Ω; R 3, t V L 2 0, T M ; W 1,2 0 Ω; R defined on a time interval [0, T M, where T M > 0 is finite or infinite and depends on the initial data. Moreover, for any TM T M, there is a constant c = ctm such that r L 0,T M ;W 1,6 Ω + t r L 0,T M ;L6 Ω + 1/r L Q T V L 0,TM ;W 2,2 Ω;R 3 + V L 2 0,TM ;W 2,6 Ω;R 3 + t V L 2 0,TM ;W 1,2 Ω c r 0 W 1,6 Ω + V 0 W 2,2 Ω. Step 2 By virtue of the weak-strong uniqueness result stated in [10, Theorem 4.1] see also [14, Theorem 4.6], the weak solution r, V coincides on the time interval [0, T M with the strong solution, the existence of which is claimed in the previous step. According to Sun, Wang, Zhang [26, Theorem 1.3], if T M then lim sup rt L Ω =. t T M 4
7 Since 2.8 holds, we infer that T M = T. At this point we conclude that couple r, V possesses regularity and that that the bound 2.13 holds with c dependent solely on T. Step 3 Since the initial data enjoy the regularity and compatibility conditions stated in , a straightforward bootstrap argument gives rise to better bounds, specifically, the solution belongs to the Valli- Zajaczkowski see [28, Theorem 2.5] class r C[0, T ]; W 3,2 Ω, t r L 2 0, T ; W 2,2 Ω, 2.14 V C[0, T ]; W 3,2 Ω L 2 0, T ; W 4,2 Ω; R 3, t V L 2 0, T ; W 2,2 Ω; R 3, 2.15 where, similarly to the previous step, the norms depend only on the initial data, r, and T. Step 4 We write equation 1.2 in the form t V 1 r div xs x V = V x V + 1 r xpr, 2.16 where, by virtue of 2.15 and a simple interpolation argument, V C 1+ν [0, T ] Ω; R 3 3, and, by the same token r C 1+ν [0, T ] Ω for some ν > 0. Consequently, by means of the standard theory of parabolic equations, see for instance Ladyzhenskaya et al. [22], we may infer that r, V is a classical solution, t V, 2 xv Hölder continuous in [0, T ] Ω and, going back to 1.1, Step 5 We write t r Hölder continuous in [0, T ] Ω x t r = x V x r V 2 xr x rdiv x V r x div x V; whence, by virtue 2.14, 2.17, 2.18, and the Sobolev embedding W 1,2 L 6, t r C[0, T ]; W 1,6 Ω Next, we differentiate 2.16 with respect to t. Denoting Z = t V we therefore obtain t Z r div xs x Z + V x Z = t div x S x V t V x V + t r r xpr, 2.20 where, in view of 2.19 and the previously established estimates, the expression on the right-hand side is bounded in C[0, T ]; L 6 Ω; R 3. Thus using the L p maximal regularity see Denk, Hieber, and Pruess [3], rylov [21] or Danchin [8, Theorem 2.2], we deduce that Finally, writing 2 t,tv = t Z L 2 0, T ; L 6 Ω; R 3, t V = Z C[0, T ]; W 1,6 Ω; R t,tr = t V x r V t x r t rdiv x V r t div x V, and using 2.19, 2.21, we obtain the desired conclusion 2 t,tr C[0, T ]; L 6 Ω. Here and hereafter, we shall use notation a b and a b. the symbol a b means that there exists c = cω, T, µ, γ > 0 such that a cb; a b means a b and b a. 5
8 2.2 Extension lemma Lemma 2.1. Under the hypotheses of Proposition 2.1, the functions r and V can be extended outside Ω in such a way that: 1 The extended functions still denoted by r and V are such that V is compactly supported in [0, T ] R 3 and r r > 0. 2 V C 1 [0,T ] R 3 ;R 3 + V C[0,T ];C 2 R 3 ;R 3 + t x V C[0,T ];L 6 R 3 ;R t,tv L 2 0,T ;L 6 R V C 1 [0,T ] Ω;R 3 + V C[0,T ];C 2 Ω;R 3 + t x V C[0,T ];L 6 Ω;R t,tv L 2 0,T ;L 6 Ω; 3 r C 1 [0,T ] R 3 + t x r C[0,T ];L 6 R 3 ;R t,tr C[0,T ];L 6 R r C 1 [0,T ] Ω + t x r C[0,T ];L 6 Ω;R t,tr C[0,T ];L 6 Ω+ V C 1 [0,T ] Ω;R 3 + V C[0,T ];C 2 Ω;R 3 + t x V C[0,T ];L 6 Ω;R t,tv L 2 0,T ;L 6 Ω; 4 t r + div x rv = 0 in 0, T R Proof: We first construct the extension of the vector field V. To this end, we follow the standard construction in the flat domain, see Adams [1, Chapter 5, Theorem 5.22] and combine it with the standard procedure of flattening of the boundary and the partition of unity technique, we get 2.22 Once this is done, we solve on the whole space the transport equation It is easy to show that the unique solution r of this equation possesses regularity and estimates stated in Remark 2.1. Here and hereafter, we denote X T R 3 a subset of L 2 0, T R 3 of couples r, V, r > 0 with finite norm r, V XT R 3 r C 1 [0,T ] R 3 + t x r C[0,T ];L 6 R 3 ;R t,tr C[0,T ];L 6 R V C 1 [0,T ] R 3 ;R 3 + V C[0,T ];C 2 R 3 ;R 3 + t x V C[0,T ];L 6 R 3 ;R t,tv L 2 0,T ;L 6 R 3 We notice that the first component of the couple belonging to X T R 3 is always strictly positive on [0, T ] R 3 and set 0 r = min t,x [0,T ] R 3rt, x, r = max t,x [0,T ] R 3rt, x Physical domain, mesh approximation The physical space is represented by a bounded domain Ω R 3 of class C 3. The numerical domains Ω h are polyhedral domains, Ω h = T, 2.27 where T is a set of tetrahedra which have the following property: If L, L, then L is either a common face, or a common edge, or a common vertex. By E, we denote the set of the faces of the element T. The set of all faces of the mesh is denoted by E; the set of faces included in the boundary Ω h of Ω h is denoted by E ext and the set of internal faces i.e E \ E ext is denoted by E int. Further, we ask V h Ω h a vertex V h Ω
9 Furthermore, we suppose that each is a tetrahedron such that ξ[] diam[] h, 2.29 where ξ[] is the radius of the largest ball contained in. The properties of this mesh needed in the sequel are formulated in the following lemma, whose easy proof is left to the reader. Lemma 2.2. There exists a positive constant d Ω depending solely on the geometric properties of Ω such that dist[x, Ω] d Ω h 2, for any x Ω h. Moreover, Ω h \ Ω Ω \ Ω h h 2. We find important to emphasize that Ω h Ω, in general. 2.4 Numerical spaces We denote by Q h Ω h the space of piecewise constant functions: Q h Ω h = {q L 2 Ω h T, q R} For a function v in CΩ h, we set v = 1 v dx for T and Π Q h vx = v 1 x, x Ω T Here and in what follows, 1 is the characteristic function of. We define the Crouzeix-Raviart space with zero traces : and with general traces E int, = L, V h,0 Ω h = {v L 2 Ω h, T, v P 1, 2.32 v ds = v L ds, E ext, V h Ω h = {v L 2 Ω, T, v P 1, E int, = L, v ds = 0}, v ds = v L ds} We denote by Π V h the standard Crouzeix-Raviart projection, and ΠV h,0 the Crouzeix-Raviart projection with zero trace, specifically, Π V h : CΩ h V h Ω h, Π V h [φ] ds x = φ ds x for all E, Π V h,0 : CΩ h V h Ω h, If v W 1,1 Ω h, we set Π V h,0 [φ] ds x = φ ds x for all E int, Each element v V h Ω h can be written in the form Π V h,0 [φ] ds x = 0 whenever E ext. v = 1 vds for E vx = E v ϕ x, x Ω h,
10 where the set {ϕ } E V h Ω h is the classical Crouzeix-Raviart basis determined by, E 2 1, ϕ ds = δ, Similarly, each element v V h,0 Ω h can be written in the form vx = E int v ϕ x, x Ω h We first recall in Lemmas the standard properties of the projection Π V h. The collection of their proofs in the requested generality can be found in the Appendix of [17] with exception of Lemma 2.8 and its Corollary 2.1. We refer to the monograph of Brezzi, Fortin [4], the Crouzeix s and Raviart s paper [6], Gallouet, Herbin, Latché [16] for the original versions of some of these proofs. We present the proof of Lemma 2.8 dealing with the comparison of projections Π V h and ΠV h,0 that we did not find in the literature. Lemma 2.3. The following estimates hold true: for all T and φ C; and for all T and φ C s. Π V h [φ] L + Π V h,0 [φ] L φ Π V h [φ] L p φ L, 2.38 h s s φ L p ;R ds, s = 1, 2, 1 p, 2.39 φ Π V h [φ] L p ;R d chs 1 s φ L p ;R ds, s = 1, 2, 1 p, 2.40 Lemma 2.4. Let 1 p. Then with any v V h Ω h. h v p v p L p Ω h, 2.41 E Lemma 2.5. The following Sobolev-type inequality holds true: v 2 L 6 Ω h x v 2 dx, 2.42 T with any v V h,0 Ω h. Lemma 2.6. There holds: T for all v C 1 Ω h, R d and all q Q h Ω h. q div Π V h [v] dx = Ω q div v dx, 2.43 Lemma 2.7 Jumps over faces in the Crouzeix-Raviart space. For all v V h,0 Ω h there holds 1 h E [v] 2,n ds x v 2 dx, 2.44 T where [v],n is a jump of v with respect to a normal n to the face, { v x v x = L E int, [v],n x = L x if n = n,, v L x v x if n = n,l n, is the normal of, that is outer w.r. to element and x E ext, [v],n x = vx, with n an exterior normal to Ω. 8
11 We will need to compare the projections Π V h and ΠV h,0. Clearly they coincide on interior elements meaning T, Ω h =. We have the following lemma for the tetrahedra with non void intersection with the boundary. Lemma 2.8. We have Π V h [φ] ΠV h,0 [φ] L + h x Π V h [φ] ΠV h,0 [φ] L ;R 3 for any φ C. sup φ L if T, Ωh, Ωh 2.45 Proof: We recall the Crouzeix-Raviart basis 2.36 and the fact that Π V h basis functions corresponding to E ext. We have and ΠV h,0 differ only in Π V h [φ] Πv h,0 [φ] L 1 ϕ E E ext φ ds L c sup φ L, 2.46 E E ext and h x Π V h [φ] ΠV h,0 [φ] L h 1 x ϕ E E ext ch sup φ L x ϕ Ω h E E ext L φ ds. L The proof is completed by E E ext x ϕ L ch 1. In fact, in the derivation of the error estimates we will use the consequence of the above observations formulated in the following two corollaries. Corollary 2.1. Let φ C 1 R 3 such that φ Ω = 0. Then we have, Π V h [φ] ΠV h,0 [φ] L = 0 if T, Ωh =, 2.47 Π V h [φ] ΠV h,0 [φ] L + h x Π V h [φ] ΠV h,0 [φ] L ;R 3 if T h, Ω h, Ω. h 2 x φ L R 3 ;R 3, 2.48 Proof: Relation 2.47 follows immediately from 2.45, as there is an empty sum on the right hand side for interior elements Ω h =. For any x Ω h there exists y Ω and thus φy = 0 such that φx dist[x, y] x φ L R 3 ;R 3 h 2 x φ L R 3 ;R 3, 2.49 where we used Lemma 2.2 for the latter inequality. The proof is completed by taking supremum over T h and combining with Note that the mesh regularity property 2.29 supplies a uniform estimate of constants c from the previous lemma, which enables to write the latter inequality in Corollary 2.2. For any φ CR 3, Π V h [φ] ΠV h,0 [φ] L p h 3/p φ L Ω h, 1 p
12 Proof: Apply inverse estimates see e.g. [20, Lemma 2.9] to We will frequently use the Poincaré, Sobolev and interpolation inequalities on tetrahedra reported in the following lemma. Lemma We have, for any v W 1,p, where 1 p. v v L p E, v v L p h v L p, 2.51 h v L p, There holds v v L p v L p, 2.53 E, v v L p for any v W 1,p, 1 p d, where p = dp d p. 3 We have, for any v W 1,p, 1 p d, where 1 q = β p + 1 β p. v L p, 2.54 v v L q ch β v L p ;R d, 2.55 v v L q ch β v L p ;R d, 2.56 We finish the section of preliminaries by recalling two algebraic inequalities 1 the imbedding inequality L a i p 1/p L a i q 1/q, 2.57 i=1 for all a = a 1,..., a L R L, 1 q p and the discrete Hölder inequality for all a = a 1,..., a L R L, b = b 1,..., b L R L, 1 q + 1 p = 1. 3 Main result i=1 L L a i b i a i q 1/q L a i p 1/p, 2.58 i=1 Here and hereafter we systematically use the following abbreviated notation: For a function v C[0, T ], L 1 Ω we set i=1 i=1 ˆφ = Π Q h [φ], φ h = Π V h [φ], φ h,0 = Π V h,0 [φ]. 3.1 v n x = vt n, x, 3.2 where t 0 = 0 t 1... t n 1 t n t n+1... t N = T is a partition of the interval [0, T ]. Finally, for a function v V h Ω h we denote h vx = T x vx1 x, div h vx = T div x vx1 x
13 In order to ensure the positivity of the approximate densities, we shall use an upwinding technique for the density in the mass equation. For q Q h Ω h and u V h,0 Ω h ; R 3, the upwinding of q with respect to u is defined, for = L E int by and we denote Up q, u Ec E int q up where a + = maxa, 0, a = mina, Numerical scheme q up = u n, = { q if u n, > 0 q L if u n, 0 E E int, 3.4 q [u n, ] + + q L [u n, ], We consider a couple ϱ n, u n = ϱ n,,h, u n,,h of numerical solutions of the following algebraic system numerical scheme: ϱ n Q h Ω h, ϱ n > 0, u n V h,0 Ω h ; R 3, n = 0, 1,..., N, 3.5 T ϱn ϱn 1 φ + T T + µ 3 E ϱ n ûn ϱ n 1 ûn 1 T T pϱ n ϱ n,up u n n, φ = 0 for any φ Q h Ω h and n = 1,..., N, E v + T E v n, + µ T 3.6 ϱ n,up û n,up [u n n, ] v 3.7 u n : v dx divu n divv dx = 0, for any v V h,0 Ω; R 3 and n = 1,..., N. The numerical solutions depend on the size h of the space discretisation and on the time step. For the sake of clarity and in order to simplify notation we will always systematically write in all formulas ϱ n, u n instead of ϱ n,,h, u n,,h. Existence of a solution to problem is well known together with the fact that any solution ϱ n 1 n N Q h Ω N satisfies ϱ n > 0 provided ϱ 0 > 0 thanks to the upwind choice in 3.6 see e.g. [15, 20]. Remark 3.1. Throughout the paper, q up in Error estimates The main result of this paper is announced in the following theorem: is defined in 3.4, where u is the numerical solution constructed Theorem 3.1. Let Ω R 3 be a bounded domain of class C 3 and let the pressure satisfy 1.4 with γ 3/2. Let {ϱ n, u n } 0 n N be a family of numerical solutions resulting from the scheme Moreover, suppose there are initial data [r 0, V 0 ] belonging to the regularity class specified in Proposition 11
14 2.1 and giving rise to a weak solution [r, V] to the initial-boundary value problem in 0, T Ω satisfying 0 rt, x r a.a. in 0, T Ω. Then [r, V] is regular and there exists a positive number C = C M 0, E 0, r, r, p C 1 [r,r], t r, r, V, t V, V, 2 V L Q T ;R 45, such that sup 1 n N where 2 t r L 1 0,T ;L γ Ω, t r L 2 0,T ;L 6γ/5γ 6 Ω;R 3, 2 t V, t V L 2 0,T ;L 6/5 Ω;R 12, Ω Ω h [ ] 1 2 ϱn û n Vt n, 2 + Hϱ n H rt n, ϱ n rt n, Hrt n + C + h a + Ω Ω h 1 n N Ω Ω h h u n x Vt n, 2 [ 1 2 ϱ0 û 0 V Hϱ 0 H r 0 ϱ 0 r 0 Hr 0 a = 2γ 3 γ if 3 2 γ 2, a = 1 2 dx ] dx 3.8 dx, otherwise. 3.9 Note that for γ = 3/2 Theorem 3.1 gives only uniform bounds on the difference of exact and numerical solution, not the convergence. 4 Uniform estimates If we take φ = 1 in formula 3.6 we get immediately the conservation of mass: n = 1,...N, ϱ n dx = Ω h ϱ 0 dx. Ω h 4.1 Next Lemma reports the standard energy estimates for the numerical scheme , see again [15, 20]. Lemma 4.1. Let ϱ n, u n be a solution of the discrete problem with the pressure p satisfying 1.4. Then there exist such that T ϱ n [minϱ n, ϱ n L, maxϱ n, ϱ n L], = L E int, n = 1,..., N, ϱ n 1,n [minϱ n 1, ϱn, maxϱ n 1, ϱn ], T, n = 1,..., N, 1 2 ϱm u m 2 + Hϱ m 1 2 ϱ0 u Hϱ 0 T m + µ x u n 2 dx + µ + λ divu n 2 dx T + [D m, u time ] + [D m, ϱ ] + [Dspace m, u ] + [Dspace m, ϱ ] = 0, 4.2 time 12
15 for all m = 1,..., N, where [D m, u time ] = [D m, ϱ time ] = [D m, u space [D m, ϱ space m T m T ] = ] = ϱ n 1 H ϱ n 1,n u n un ϱn ϱn 1 2, 4.3a 2, 4.3b m ϱ n,up u n un L 2 u n n,, 2 = L E int 4.3c m H ϱ n ϱn ϱn L 2 u n n,. 2 = L E int 4.3d We have the following corollary of Lemma 4.1 see [17, Lemma 4.1, Lemma 4.2]: Corollary 4.1. Under assumptions of Lemma 4.1, we have: 1 There exists c = cm 0, E 0 > 0 independent of n, h and such that k k N x u n 2 dx c, 4.4 N u n 2 c, 4.5 L 6 Ω h ;R 3 sup n=0,...n ϱ n ˆ u n2 L 1 Ω h c sup n=0,...n ϱ n L γ Ω h c, If the pair r, U belongs to the class 2.25 there is c = cm 0, E 0, r, r, U, U L Q T ;R 12 > 0 such that for all n = 1,..., N, sup n=0,...n Eϱ n, û n ˆrt n, Ût n c, 4.8 where Eϱ, u z, v = Ω h ϱ u v 2 + Eϱ z dx, Eϱ z = Hϱ H zϱ z Hz. 4 There exists c = cm 0, E 0, r, p C [ r,r] > 0 such that m = L E int ϱ n ϱ n L 2 m [ 1 {ϱ n 1} [max{ϱ, ϱ L }] 2 γ + 1 {ϱ n 1} ] u n n, c if γ [1, 2, 4.9 = L E int ϱ n ϱ n L 2 u n n, c if γ 2 13
16 5 Discrete relative energy inequality The starting point of our error analysis is the discrete relative energy inequality derived for the numerical scheme in [17, Theorem 5.1]. Lemma 5.1. Let ϱ n, u n be a solution of the discrete problem with the pressure p satisfying 1.4. Then there holds for all m = 1,..., N, 1 2 ϱ m u m U m 2 ϱ 0 u 0 U Eϱ m r m Eϱ 0 r 0 T T m + µ x u n U n 2 dx + µ divu n U n 2 dx T i, 3 T for any 0 r n Q h Ω h, U n V h,0 Ω h ; R 3, n = 1,..., N, where m T 1 = µ x U n : x U n u n dx + µ T 3 m T 2 = ϱ n 1 U n U n 1 U n 1 + U n 2 T 3 = T 4 = T 5 = T 6 = T m T m T m T m T E = L E = L rn ϱ n E = L U ϱ n,up n + U n L 2 pϱ n [U n n, ], H r n H r n 1, ϱ n,up H r n 1 [un n, ]. i=1 divu n divu n u n dx, u n 1, û n,up U n [u n n, ], 6 Approximate discrete relative energy inequality In this section, we transform the right hand side of the relative energy inequality 5.1 to a form that is more convenient for the comparison with the strong solution. This transformation is given in the following lemma. Lemma 6.1 Approximate relative energy inequality. Let ϱ n, u n be a solution of the discrete problem , where the pressure satisfies 1.4 with γ 3/2. Then there exists c = c M 0, E 0, r, r, p C 1 [r,r], t r, r, V, t V, V L Q T ;R 18, t 2 r L 1 0,T ;L γ Ω, t r L 2 0,T ;L 6γ/5γ 6 Ω;R 3 > 0, such that for all m = 1,..., N, we have: ϱ m û m ˆV m h,0 2 + Eϱ m ˆr m dx Ω h + m T Ω h µ x u n V n h,0 2 dx + µ 3 ϱ 0 û 0 ˆV 0 h,0 2 + Eϱ 0 ˆr 0 dx 14 divu n V n h,0 2 dx 6 S i + Rh, m +Gm, i=
17 for any couple r, V belonging to the class 2.25, where S 1 = S 2 = S 3 = m T m T m T E µ x V n h,0 : xv n h,0 un dx + µ divv n h,0 3 divv n h,0 un dx, ϱ n 1 V n h,0, V n 1 h,0, m S 4 = pϱ n div V n dx, Ω h m S 5 = ˆr n ϱ n p ˆr n Ω h ˆr n [ tr] n dx, m ϱ S 6 = Ωh n ˆr n p ˆr n u n r n dx, V n h,0, un, ϱ n,up n,up ˆV h,0, û n,up V n h,0, V n n,up h,0, ˆV h,0, n,, 6.2 and G m c m Eϱ n, û n ˆr n, ˆV n, Rh, m c + h a, 6.3 with the power a defined in 3.9 and with the functional E introduced in 4.8. Proof: We take as test functions U n = V n h,0 and rn = ˆr n in the discrete relative energy inequality 5.1. We keep the left hand side and the first term term T 1 at the right hand side as they stay. The transformation of the remaining terms at the right hand side terms T 2 T 6 is performed in the following steps: Step 1: Term T 2. We have T 2 = T 2,1 + R 2,1 + R 2,2, with T 2,1 = m T ϱ n 1 V n h,0, V n 1 h,0, V n h,0, un, 6.4 and where and R 2,1 = m T R n, 2,1 = V n 2 ϱn 1 h,0, V n 1 h,0, 2 R n 2,2 = T R n, 2,1, R 2,2 = = ϱ n 1 V n h,0, V n 1 h,0, 2 ϱn 1 m R n 2,2, [V n V n 1 ] h,0, 2, u n 1 un. We may write by virtue of the first order Taylor formula applied to function t Vt, x, = 1 [V n V n 1 ] h,0, = 1 [ 1 [ t n ] t Vz, xdz t n 1 h,0 ] dx [ 1 tn ] [ t Vz t n 1 ]h,0 xdz dx [ t V] h,0 L 0,T ;L Ω;R 3 tv, L 0,T ;L Ω;R 3 15
18 where we have used the property 2.38 of the projection Π V h,0 on the space V h,0ω h. Therefore, thanks to the mass conservation 4.1, we get To treat term R n 2,2 R n, 2,1 M 0 2 tv 2 L 0,T ;L Ω;R we use the discrete Hölder inequality and identity 4.1 in order to get R n 2,2 cm 0 t V 2 L 0,T ;W 1, Ω;R 3 1/2 + cm0 T ϱ n 1 un 1 un 2 1/2 t V L 0,T ;L Ω;R 3 ; whence, by virtue of estimate 4.2 for the upwind dissipation term 4.3a, one obtains R 2,2 cm 0, E 0, t V L Q T ;R Step 2: Term T 3. Employing the definition 3.4 of upwind quantities, we easily establish that T 3 = T 3,1 + R 3,1, m with T 3,1 = T E ϱ n,up û n,up n,up ˆV h,0, V n h,0, un n,, R 3,1 = m R n, 3,1, E int and R n, 3,1 = V n ϱn h,0, V n h,0,l 2 [u n n, ] + + ϱ n V n h,0,l V n h,0, 2 L [u n n,l ] +, = L E int. 2 2 Writing V n h,0, V n h,0,l = [V n h,0 V n h ] + V n h, V n h + V n h V n h, +V n h, V n h + V n h V n h,l + [V n h V n h,0 ] L, = L E int, and employing estimates 2.47 if Ω h =, 2.48 if Ω h to evaluate the L -norm of the first term, 2.51 then 2.40 s=1 and 2.52 after 2.40 s=1 to evaluate the L -norm of the second and third terms, and performing the same tasks at the second line, we get consequently whence V n h,0, V n h,0,l L L;R 3 ch V L L;R 9 ; 6.7 R n, 3,1 h2 c V 2 L 0,T Ω;R 9 ϱn + ϱ n L u n, = L E int, R 3,1 h c V 2 L 0,T Ω;R 9 [ m T E T = L E h ϱ n + ϱ n L 6/5 5/6 h u n 6 1/3] 1/2 h cm0, E 0, V L Q T ;R 9, 6.8 provided γ 6/5, thanks to the discrete Hölder inequality, the equivalence relation 2.29, the equivalence of norms 2.41 and energy bounds listed in Corollary 4.1. Clearly, for each face = L E int, u n n, + u n n,l = 0; whence, finally T 3,1 = m T E ϱ n,up û n,up Before the next transformation of term T 3,1, we realize that n,up ˆV h,0, V n h,0, V n h,0, u n n,. 6.9 V n h,0, Vn h,0, = [Vn h,0 Vn h ] + V n h, Vn h + Vn h Vn h, + [Vn h Vn h,0 ] ; 16
19 whence by virtue of , and 2.40 s=1, similarly as in 6.7, V n h,0, V n h,0, L ;R 3 ch xv L 0,T ;L Ω:R 3, Let us now decompose the term T 3,1 as T 3,1 = T 3,2 + R 3,2, with R 3,2 = T 3,2 = R n 3,2 = T m T E E ϱ n,up m R n 3,2, ϱ n,up n,up ˆV h,0, û n,up V n h,0, V n h,0, û n,up n,, and n,up ˆV h,0, û n,up V n h,0, V n h,0, By virtue of discrete Hölder s inequality and estimate 6.10, we get R n 3,2 c V L Q T ;R 9 T E h ϱ n,up T E h ϱ n,up γ 0 1/2γ0 û n,up T E u n û n,up n,up ˆV h,0, 2 1/2 n,. h u n û n,up q 1/q, where γ q = 1, γ 0 = min{γ, 2} and γ 3/2. For the sum in the last term of the above product, we have h u n û n,up q c h u n u n q T E T E c u n u n q L q ;R 3 + u n u n q ch 2γ 0 3 q q/2, L q ;R 3 2γ 0 x u n 2 L 2 ;R 9 T T T E where we have used the definition 3.4, the discrete Minkowski inequality, interpolation inequalities and the discrete imbedding inequality Now we can go back to the estimate of R n 3,2 taking into account the upper bounds 4.4, , in order to get provided γ 3/2, where a is given in 6.3. Finally, we rewrite term T 3,2 as R 3,2 h a cm 0, E 0, V L Q T ;R 9, 6.11 T 3,2 = T 3,3 + R 3,3, with R 3,3 = T 3,3 = R n 3,3 = T m T E E ϱ n,up m R n 3,3, ϱ n,up n,up ˆV h,0, û n,up V n h,0, V n n,up h,0, ˆV h,0, n,, and n,up ˆV h,0, û n,up V n h,0, V n h,0, û n,up n,up ˆV h,0, n, ; 6.12 whence m R 3,3 c V L Q T,R 9 Eϱ n, û n ˆr n, ˆV n h,
20 Step 3: Term T 4. Integration by parts over each T gives We may write T 4 = m T pϱ n div x V n h,0 dx. div x V n 0,h V n h L ch x V L 0,T ;L Ω;R 9, 6.14 where we have used Therefore, employing identity 2.43 we obtain T 4 = T 4,1 + R 4,1, R 4,1 = T 4,1 = m T m T pϱ n div x V n h,0 V n h dx. pϱ n div x V n dx, 6.15 Due to 1.4 and 4.7, pϱ n is bounded uniformly in L L 1 Ω; employing this fact and 6.14 we immediately get R 4,1 h ce 0, M 0, V L 0,T ;L Ω;R Step 4: Term T 5. Using the Taylor formula, we get H r n H r n 1 = H r n r n r n 1 where r n [minrn 1, rn, maxrn 1, rn ]; we infer T 5 = T 5,1 + R 5,1, with T 5,1 = R n, 2 5,1 = 1 2 H r n rn rn 1 m T ϱ n r n. r n ϱ n p r n r n 1 2 H r n r n r n 1 2, r n rn 1, R 5,1 = m T R n, 5,1, and Consequently, by the first order Taylor formula applied to function t rt, x on the interval t n 1, t n and thanks to the mass conservation 4.1 R 5,1 cm 0, r, r, p C 1 [r,r], t r L Q T Let us now decompose T 5,1 as follows: T 5,1 = T 5,2 + R 5,2, with T 5,2 = R n, 5,2 = r n ϱ n p r n r n m T r n r n 1 r n ϱ n p r n r n [ t r] n dx. [ t r] n dx, R 5,2 = m T R n, 5,2, and 6.18 In accordance with 3.2, here and in the sequel, [ t r] n x = t rt n, x. We write using twice the Taylor formula in the integral form and the Fubini theorem, R n, 5,2 = 1 tn tn p r n rϱ n n r n t 2 rzdzdsdx p r n r n tn t n 1 t n 1 ϱ n r n t 2 rz dxdzds 18 s
21 p r n tn r n ϱ n ˆr n L γ t 2 rz L γ dzds. t n 1 Therefore, by virtue of Corollary 4.1, we have estimate R 5,2 cm 0, E 0, r, r, p C 1 [r,r], t 2 r L 1 0,T ;L γ Ω Step 5: Term T 6. We decompose this term as follows: m T 6 = T 6,1 + R 6,1, R 6,1 = T 6,1 = R n,, 6,1 = m T = L E ϱ n,up ϱ n T E ϱ n H r n 1 H r n 1 H r n 1 R n,, 6,1, with H r n 1 u n n,, and u n n,, for = L E int. We will now estimate the term R n,, 6,1. We shall treat separately the cases γ 2 and γ 2. The simple case γ 2 is left to the reader. The more complicated case γ 2 will be treated as follows: We first write R n,, 6,1 h H r L Q T ;R 3 ϱn,up ϱ n [ 1 {ϱ n 1}[max{ϱ, ϱ L }] 2 γ/2 + 1 {ϱ n 1} [ 1 ] {ϱ n 1} [max{ϱ, ϱ L }] 2 γ/2 + 1 {ϱ n 1} u n n, ] h u n n,, where we have employed the first order Taylor formula applied to function x H rt n 1, x. Consequently, the application of the discrete Hölder and Young inequalities yield R 6,1 h c H r L Q T ;R 3 m [ ] 1/2 h 1 {ϱ n 1}[max{ϱ, ϱ L }] 2 γ + 1 {ϱ n 1} u n n, T E T = L E hϱ n,up ϱ n 2[ 1 {ϱ n 1} [max{ϱ, ϱ L }] 2 γ + 1 {ϱ n 1} ] u n n, h c H r L Q T ;R 3 m {[ Ω h hϱ n 6 2 γ 5 ] 6 5 T = L E T T E h u n n, hϱ n,up ϱ n 2[ 1 {ϱ n 1} [max{ϱ, ϱ L }] 2 γ + 1 {ϱ n 1} ] u n n, } 1/2 { m h c H r L Q T ;R 3 [ Ω h hϱ n 6 2 γ 5 ] 6 5 h u n 6 1/6 + m [ T = L E h cm 0, E 0, r, r, p C[r,r], r L Q T ;R 3, T E 1/2 hϱ n,up ϱ n 2[ 1 {ϱ n 1} [max{ϱ, ϱ L }] 2 γ + 1 {ϱ n 1} ] u n n, where, in order to get the last line, we have used the estimate 4.9 of the numerical dissipation to evaluate the second term, and finally equivalence of norms 2.41 p=6 together with 4.5 and 4.7, under assumption γ 12/11, to evaluate the first term. ]} 19
22 Let us now decompose the term T 6,1 as T 6,1 = T 6,2 + R 6,2, with T 6,2 = R 6,2 = m R n,, 6,2 = ϱ n E H r n 1 R n,, 6,2, and H r n 1 m T = L E H r n 1 rn 1 ϱ n H r n 1 rn 1 rn 1 [u n n, ], rn 1 [u n n, ]. Therefore, by virtue of the second order Taylor formula applied to function H, the Hölder inequality, 2.41, and 4.5, 4.7 in Corollary 4.1, we have, provided γ 6/5, R 6,2 hc H C[r,r] + H C[r,r] r m L QT ;R 3 ϱ n L γ Ω h u n L 6 Ω h ;R 3 h cm 0, E 0, r, r, p C 1 [r,r], r L Q T ;R Let us now deal with the term T 6,2. Noting that may write T 6,2 = T 6,3 + R 6,3, with T 6,3 = R 6,3 = m T m T + ϱ n H r n 1 un r n 1 dx, r n 1 dx = ϱ n H r n 1 un u n r n 1 dx m T E E r n 1 ϱ n H r n 1 rn 1 rn 1 u n u n n,. r n 1 n,, we Consequently, by virtue of Hölder s inequality, interpolation inequality 2.55 to estimate u n u n L γ 0 ;R 3 by h 5γ 0 6/2γ 0 x u n L 2 ;R 9, γ 0 = min{γ, 2} in the first term, and by the Taylor formula applied to function x rt n 1, x, then Hölder s inequality and to estimate u n u n L γ 0 ;R 3 by h 5γ 0 6/2γ 0 x u n L 2 ;R 9, we get R 6,3 h b cm 0, E 0, r, r, p C 1 [r,r] r L Q T ;R 3, b = 5γ 0 6 2γ 0, 6.21 provided γ 6/5, where we have used at the end the discrete imbedding and Hölder inequalities and finally estimates 4.4 and 4.7. Finally we write T 6,3 = T 6,4 + R 6,4, with T 6,4 = R 6,4 = m T m T where by the same token as in 6.19, ϱ n ϱ n p r n r n u n r n dx, H r n r n H r n 1 rn 1 u n dx, 6.22 R 6,4 cm 0, E 0, r, r, p C 1 [r,r], r, t r L Q T ;R 4, t r L 2 0,T ;L 6γ/5γ 6 Ω;R 3,
23 provided γ 6/5. We are now in position to conclude the proof of Lemma 6.1: we obtain the inequality 6.1 by gathering the principal terms 6.4, 6.12, 6.15, 6.18, 6.22 and the residual terms estimated in 6.5, 6.6, 6.8, 6.11, 6.13, 6.17, 6.19, 6.20, 6.21, 6.23 at the right hand side 6 i=1 T i of the discrete relative energy inequality A discrete identity satisfied by the strong solution This section is devoted to the proof of a discrete identity satisfied by any strong solution of problem in the class extended eventually to R 3 according to Lemma 2.1. This identity is stated in Lemma 7.1 below. It will be used in combination with the approximate relative energy inequality stated in Lemma 6.1 to deduce the convenient form of the relative energy inequality verified by any function being a strong solution to the compressible Navier-Stokes system. This last step is performed in the next section. Lemma 7.1 A discrete identity for strong solutions. Let ϱ n, u n be a solution of the discrete problem with the pressure satisfying 1.4, where γ 3/2. There exists c = c M 0, E 0, r, r, p C 1 [r,r], t r, r, V, t V, V, 2 V L Q T ;R 45, t 2 r L 1 0,T ;L γ Ω, t r L 2 0,T ;L 6γ/5γ 6 Ω;R 3, 2 t V, t V L 2 0,T ;L 6/5 Ω;R 12 > 0, such that for all m = 1,..., N, we have: 6 S i + R m h, = 0, 7.1 i=1 where and S 1 = S 2 = S 3 = S 4 = S 5 = 0, S 6 = m T m T m T E m m µ x V n h,0 : xv n h,0 un dx + µ divv n h,0 3 divv n h,0 un dx, r n 1 V n h,0, V n 1 h,0, V n h,0, un, r n,up n,up ˆV h,0, û n,up V n h,0, V n n,up h,0, ˆV h,0, n, Ω h pˆr n div V n dx, Ω h p ˆr n u n r n dx, R m h, h c 5/6 +, for any couple r, V belonging to 2.25 and satisfying the continuity equation 1.1 on 0, T R 3 and momentum equation 1.2 with boundary conditions 1.5 on 0, T Ω in the classical sense. 21
24 Before starting the proof we recall an auxiliary algebraic inequality whose straightforward proof is left to the reader, and introduce some notations. Lemma 7.2. Let p satisfies assumptions1.4. Let 0 a b. Then there exists c = ca, b > 0 such that for all ϱ [0, and r [a, b] there holds Eϱ r ca, b 1 R+ \[a/2,2b]ϱ + ϱ γ 1 R+ \[a/2,2b]ϱ + ϱ r 2 1 [a/2,2b] ϱ, where Eϱ r is defined in 4.8. If we take in Lemma 7.2 ϱ = ϱ n x, r = ˆr n x, a = r, b = r where r is a function belonging to class 2.25 and r, r are its lower and upper bounds, respectively, we obtain Eϱ n x ˆr n x cr, r 1 R+ \[r/2,2r]ϱ n x+ϱ n γ x1 R+ \[r/2,2r]ϱ n x+ϱ n x ˆr n x 2 1 [r/2,2r] ϱ n x. 7.2 Now, for fixed numbers r and r and fixed functions ϱ n, n = 0,..., N, we introduce the residual and essential subsets of Ω relative to ϱ n as follows: and we set Ness n = {x Ω 1 2 r ϱn x 2r}, Nres n = Ω \ Ness, n 7.3 [g] ess x = gx1 N n ess x, [g] res x = gx1 N n res x, x Ω, g L 1 Ω. Integrating inequality 7.2 we deduce cr, r [ 1 [ϱ T ]res + n γ] [ϱ + n ˆr n] 2 dx Eϱ n, u n ˆr n, V n, 7.4 res ess for any pair r, V belonging to the class 2.25 and any ϱ n Q h Ω h, ϱ n 0. We are now ready to proceed to the proof of Lemma 7.1. Proof: Since r, V satisfies 1.1 on 0, T Ω and belongs to the class 2.25, Equation 1.2 can be rewritten in the form r t V + rv V + pr µ V µ/3 div V = 0 in 0, T Ω. From this fact, we deduce the identity where R 0 = T 1 = T 3 = T 5 = 0, m m m Ω h \Ω Ω h 5 T i = R 0, 7.5 i=1 r n [ t V ] n + rv n V n + pr n µ V n µ 3 div V n V n h,0 un dx, µ V n + µ 3 div V n V n h,0 un dx, r n V n V n V n h,0 un dx, Ω h m T 6 = pr n u n dx. Ω h T 4 = m T 2 = m Ω h pr n V n h,0 dx, Ω h r n [ t V ] n V n h,0 un dx, 22
25 In the steps below, we deal with each of the terms R 0 and T i. Step 0: Term R 0. By the Hölder inequality m R 0 Ω h \ Ω 5/6 cr, p C[r,r], t r, r, V, V, 2 V L Q T ;R 43 u n L 6 Ω h + V n h,0 L 6 Ω h h 5/3 cm 0, E 0, r, p C[r,r], t r, r, V, V, 2 V L Q T ;R 43, 7.6 where we have used 4.5 and , Step 1: Term T 1. Integrating by parts, we get: T 1 = T 1,1 + R 1,1, m with T 1,1 = T µ V n h,0 : V n h,0 un + µ 3 div V n h,0 divv n h,0 un dx, and R 1,1 = I 1 + I 2, with m I 1 = µ V n V n h,0 : V n h,0 un + µ 3 divv n V n h,0 divv n h,0 un dx, I 2 = = T m T E m E µn V n µn, V n V n h,0 un + µ 3 div V n V n h,0 un n, ds [V n h,0 un],n + µ 3 div V n[ V n h,0 un],n n ds, 7.7 where in the last line n is the unit normal to the face and [ ],n is the jump over sigma with respect to n defined in Lemma 2.7. To estimate I 1, we use the Cauchy-Schwartz inequality, decompose V n Vh,0 n = Vn Vh n + Vn h Vh,0 n and employ estimates 2.40 s=2, to evaluate the norms involving V n Vh,0 n, and decompose Vh,0 n = Vn h,0 Vn h +Vn h use , 2.39 s=1, 4.4, the Minkowski inequalityt to estimate the norms involving Vh,0 n un. We get I 1 h cm 0, E 0, V, 2 V L 0,T ;L Ω;R 36. Since the integral over any face E int of the jump of a function from V h,0 Ω h is zero, we may write I 2 = m E int [ ] µn V n V n u n V n h,0,n + µ [ div V n div V n u n V n h,0 n ds; 3 ],n whence by using the first order Taylor formula applied to functions x V n x to evaluate the differences V n V n, div V n [div V n ], and Hölder s inequality, I 2 h c 2 V L Q T ;R 27 h c 2 V L Q T ;R 27 m m E int E int 1 ] h h [u n V n L h,0,n 2 ;R 3 h + 1 h [ u n V n 2 h,0 ],n L 2 ;R 3. 23
26 Therefore, R 1,1 h cm 0, E 0, V, V, 2 V L Q T,R 39, 7.8 where we have employed Lemma 2.7, 4.4 and , Step 2: Term T 2. Let us now decompose the term T 2 as T 2 = T 2,1 + R 2,1, m with T 2,1 = and R n, 2,1 = T r n 1 V n V n 1 r n r n 1 [ t V] n V n h,0 un dx + The remainder R n, 2,1 can be rewritten as follows R n, 2,1 = [ t n ] t rt, dt [ t V] n V n h,0 un dx+ 1 t n 1 whence, by the Hölder inequality, V n h,0 un dx, R 2,1 = m T r n 1 [ t V ] n V n V n 1 r n 1[ t n t n 1 tn s R n, 2,1, V n h,0 un dx. ] t 2 V z, dzds V n h,0 un dx; R n, 2,1 [ r L Q T + t r L Q T t V L Q T ;R 3 5/6 u n L 6 + V n h,0 L 6 ] + t 2 V n L 6/5 Ω;R 3 un L 6 + V n h,0 L 6. Consequently, by the same token as in 6.19 or 6.23, R 2,1 c M 0, E 0, r, t r, V, t V, V L Q T ;R 16, 2 t V L 2 0,T ;L 6/5 Ω;R 3, 7.9 where we have used the discrete Hölder and Young inequalities, the estimates 2.38, and the energy bound 4.4 from Corollary 4.1. Step 2a: Term T 2,1. We decompose the term T 2,1 as therefore, T 2,1 = T 2,2 + R 2,2, m with T 2,2 = and R n, 2,2 = T R n 2,2 = T r n 1 V n V n 1 V n h,0 un dx, R 2,2 = r n 1 r n 1 n V V n 1 V n h,0 un dx; m T R n, 2,2 h c r L Q T ;R 3 tv L Q T ;R 3 un V n h,0 L 6 Ω;R 3. R n, 2,2, Consequently, by virtue of formula 4.5 for u n and estimates 2.38, , R 2,2 h cm 0, E 0, r, V, t V, V L Q T ;R
27 Step 2b: Term T 2,2. We decompose the term T 2,2 as T 2,2 = T 2,3 + R 2,3, m with T 2,3 = T r n 1 V n h,0, V n 1 h,0, and R n, V 2,3 = r n 1 n V n 1 [ V n V n 1 ] [ V n V n 1 [ V + ]h n V n 1 ] + r n 1 r n 1 We calculate carefully [ V n V n 1 [ V ]h, n V n 1 ] I 3 = 1 rn 1 1 rn 1 tn t n 1 { t n t n 1 h, V n h,0 un dx, R 2,3 = h V n h,0 un dx V n h,0 un dx h,0, m T V n h,0 u n dx= I 1 + I 2 + I 3. ] } [[ t V z] h [ t V z] h,0 V n h,0 un dz dx [[ ] t V z] h [ t V z] h,0 V n L 6/5 ;R 3 h,0 un L 6 ;R 3 dz. Summing over polyhedra T we get simply by using the discrete Sobolev inequality T I 3 1 rn 1 tn t n 1 1 rn 1 { 1/6 V n h,0 un 6 6/5 [ L 6 ;R 3 t V z] h [ t V z] h,0 tn T T V n h,0 un L 6 Ω h ;R 3 [ tv z] h [ t V z] h,0 dz t n 1 L 6/5 Ω h ;R 3 h5/6 tn t n 1 V n h,0 un L 6 Ω h ;R 3 tv z L Ω h ;R 3 dz, where we have used estimate 2.50 to obtain the last line. As far as the term I2 is concerned, we write I2 = 1 [ t n [ rn 1 t V zdz t n 1 ]h tn ] t V zdz t n 1 h rn 1 tn t n 1 ] x [ t V z h h, L 6/5 ;R 3 un V n h,0 L 6 ;R 3, R n, 2,3, L 6/5 ;R 3 u n V n h,0 dx 5/6 } dz where we have used the Fubini theorem, Hölder s inequality and 2.51, 2.40 s=1. Further, employing the Sobolev inequality on the Crouzeix-Raviart space V h Ω h 2.42, the Hölder inequality and estimate 2.40 s=1, we get T I 2 h rn 1 un V n h,0 L 6 Ω h ;R 3 tn t n 1 x t V z dz. L 6/5 Ω h ;R 3 We reserve the similar treatment to the term I1. Resuming these calculations and summing over n from 1 to m we get by using Corollary 4.1 and estimates , 2.38, R 2,3 h 5/6 cm 0, E 0, r, V, V, t V L Q T ;R 16, t V L 2 0,T ;L 6/5 Ω;R
28 Step 2c: Term T 2,3. We rewrite this term in the form T 2,3 = T 2,4 + R 2,4, R 2,4 = with T 2,4 = and R n, 2,4 = m T m T R n, 2,4, r n 1 V n h,0, V n 1 h,0, r n 1 V n h,0, V n 1 h,0, u n V n h,0, dx, u n u n V n h,0 V n h,0, dx First, we estimate the L norm of V n h,0, n 1 V h,0, as in 6.5. Next, we decompose V n h,0 V n h,0, = V n h,0 V n h + V n h V n h, + [V n h V n h,0 ], and use 2.51 p=2 to estimate u n u n, 2.51 p=, 2.40 s=1 to estimate V n h V n h, and to evaluate [V n h V n h,0 ] L ;R 3 V n h V n h,0 L ;R 3. Thanks to the Hölder inequality and 4.4 we finally deduce R 2,4 h cm 0, E 0, r, V, t V, V L Q T ;R Step 3: Term T 3. Let us first decompose T 3 as T 3 = T 3,1 + R 3,1, m with T 3,1 = and R n, 3,1 = We have T r n V n h,0, V n V n h,0, un dx, R 3,1 = m T R n, 3,1, r n rv n n V n V n h,0 un dx + rv n n V n h,0 V n V n h,0 un dx + rv n n h,0 V n h,0, V n V n h,0 un dx + rv n n h,0, V n V n h,0 V n h,0, un u n dx. by the Taylor formula, by virtue of 2.39 s=1 and , r n r n L V n V n h,0 L ;R 3 h r n L, h V n L ;R 9, V n h,0 V n h,0, L ;R 3 V n h,0 V n h, L ;R 3 + V n h V n h, L ;R 3 + [V n h V n h,0 ] L ;R 3 by virtue of 2.51, 2.39 s= s=1 and , u n u n L ;R 3 h V n L ;R 9 h u n L ;R 9. Consequently by employing several times the Hölder inequality for integrals over and the discrete Hölder inequality for the sums over T, and using estimate 4.4, we arrive at R 3,1 h cm 0, E 0, r, r, V, V L Q T ;R
29 Now we shall deal wit term T 3,1. Integrating by parts, we get: r n V n h,0, V n V n h,0, un dx = E = E thanks to the the fact that E V n h, n,ds = 0. Next we write We may write T 3,2 = and R n 3,2 = T + T E + T E m T E E ˆr n,up T 3,1 = T 3,2 + R 3,2, r n [V n h,0, n,]v n V n h,0, un r n [V n h,0, n,]v n V n h, V n h, un, R 3,2 = m R n 3,2, ˆr n,up n,up [ ˆV h,0, n, ]V n V n n,up h, ˆV h,0, û n,up, 7.15 r n ˆr n,up [V n h,0, n,]v n V n h, V n h,0, un [ V n n,up h,0, ˆV h,0, n, ]V n V n h, V n h, un ˆr n,up n,up [ ˆV h,0, n, ]V n V n h, V n n,up h,0, ˆV h,0, u n û n,up h,. V n V n h,0, = Vn V n + V n V n h + Vn h Vn h, + [Vn h Vn h,0 ], and use several times the Taylor formula along with 2.39 s=1, 2.51, 2.40 s=1, in order to estimate r n ˆrn,up, V n Vh,0, n, Vn n,up h, ˆV h, to get the bound 3 R n 3,2 h c r W 1, Ω 1 + V W 1, Q T ;R 3 h u n We have by the Hölder inequality T 2 +c r W 1, Ω 1 + V W 1, Q T ;R 3 h u n c + + T T E T E T u n 6 L 6 ;R 3 T E T h u n u n. h u n 6 1/6 [ c u n u n 6 L 6 ;R 3 1/6 ] c T T u n 2 L 2 ;R 9 1/2, [ h u n u n c u n u n 2 L 2 ;R 3 T 1/2 1/6 1/2 ] 1/2, u n u n 2 h c u L 2 ;R 3 n 2 L 2 ;R 9 27 T
30 where we have used 2.53 p=2, p=2. Consequently, we may use 4.4 to conclude R 3,2 h c M 0, E 0, r, r, V, V L Q T ;R Finally, we replace in T 3,2 V n V n h, by V n h,0, V n h,0,. We get T 3,2 = T 3,3 + R 3,3, R 3,3 = m R n 3,3, and T 3,3 = R n 3,3 = T m T E E ˆr n,up n,up [ ˆV h,0, n, ]V n h,0, V n n,up h,0, ˆV h,0, û n,up, 7.17 r n ˆr n,up V n h,0, n, [V n V h,0 ] n [V n h V n h,0 ] n,up ˆV h,0, û n,up, committing error R n 3,3 = h c M 0, E 0, r, r, V, V L Q T ;R 15, 7.18 as in the previous step. Step 4: Terms T 4 whence We write T 4 = T 4,1 + R 4,1, T 4,1 = pr n V n dx, Ω h R 4,1 = pr n V n V n h,0 dx; Ω h R 4,1 hcr, p C[r,r], r L Q T ;R 3, 7.19 by virtue of 2.39 s=1, Next, employing the integration by parts T 4,2 = T 4,2 + R 4,2, T 4,2 = pr n div V n dx, Ω h R 4,2 = T Writing E, Ω h pr n V n n, ds = T E, Ω h V n V n h,0, = Vn V n h + Vn h V n h, + [V n h V n h,0 ], we deduce by using 2.39 s=1, 2.40 s=1, 2.52 p=, 2.47, 2.48, V n V n h,0, L ;R 3 h V n L ;R 3,. pr n V n V n h,0, n, ds. Now, we employ the fact that whence T E, Ω h ds 1; R 4,2 hcr, p C[r,r], V L Q T ;R
31 Finally, T 4,2 = T 4,3 + R 4,3, T 4,3 = pˆr n div V n dx, Ω h R 4,3 = pr n pˆr n div V n dx; Ω h 7.21 whence R 4,3 hc p C[r,r], r, V L Q T ;R Step 5: Term T 6 We decompose T 6 as m T 6 = T 6,1 + R 6,1, with T 6,1 = m R 6,1 = T T p ˆr n p r n u n r n dx; p ˆr n u n r n dx, Consequently, by the Taylor formula, Hölder inequality and estimate 4.5, 7.23 R 6,1 h cm 0, E 0, r, r, p C 1 [r,r], r L Q T ;R Gathering the formulae 7.7, 7.12, 7.17, 7.21, 7.23 and estimates for the residual terms 7.8, , , 7.19, 7.20, 7.22, 7.24 concludes the proof of Lemma A Gronwall inequality In this Section we put together the relative energy inequality 6.1 and the identity 7.1 derived in the previous section. The final inequality resulting from this manipulation is formulated in the following lemma. Lemma 8.1. Let ϱ n, u n be a solution of the discrete problem with the pressure satisfying 1.4, where γ 3/2. Then there exists a positive number c = c M 0, E 0, r, r, p C 1 [r,r], t r, r, V, t V, V, 2 V L Q T ;R 45, 2 t r L 1 0,T ;L γ Ω, t r L 2 0,T ;L 6γ/5γ 6 Ω;R 3, 2 t V, t V L 2 0,T ;L 6/5 Ω;R 12 such that for all m = 1,..., N, there holds: Eϱ m, u m ˆr m, ˆV m h,0+ µ 2 m T [ c h a + + Eϱ 0, u 0 ˆr0, ˆV ] h,0 0 + c x u n V n h,0 2 dx m Eϱ n, u n ˆr n, ˆV n h,0, with any couple r, V belonging to 2.25 and satisfying the continuity equation 1.1 on 0, T R 3 and momentum equation 1.2 with boundary conditions 1.5 on 0, T Ω in the classical sense, where a is defined in 6.3 and E is given in 4.8. Proof. We observe that m S 6 S 6 = p Ωh ˆr n ˆrn ϱ n m ˆr n V n r n dx + p Ωh ˆr n ˆrn ϱ n ˆr n u n V n r n dx., 29
32 Gathering the formulae 6.1 and 6.2, one gets Eϱ m, u m ˆr m, ˆV m h,0 Eϱ 0, u 0 ˆr0, ˆV h,0 0 + µ where P 1 = P 2 = P 3 = P 4 = m T m ϱ n 1 rn 1 T = L E m Ω h m T n V h,0, V n 1 ϱ n,up m T h,0, u n V n 0,h 2 L 2 ;R 3 V n h,0, un, 4 P i + Q, 8.1 ˆr n,up n,up ˆV h,0, û n,up V n h,0, V n h,0, V n,up h,0, n,, pϱ n p ˆr n ϱ n ˆr n pˆr n divv n, Q = R m h, + Rm h, + Gm. p ˆr n ˆrn ϱ n ˆr n u n V n r n dx, Now, we estimate conveniently the terms P i, i = 1,..., 4 in four steps. n 1 V h,0, Step 1: Term P 1. We estimate the L norm of V n h,0, by L norm of t V in the same manner as in 6.5. According to Lemma 7.2, ϱ r γ 1 R+ \[r/2,2r]ϱ cpe p ϱ r, with any p 1; in particular, ϱ r 6/5 1 R+ \[r/2,2r]ϱ ceϱ r 8.2 provided γ 6/5. We get by using the Hölder inequality, [ T ϱ n 1 T ϱ n 1 rn 1 h,0, n V h,0, V n 1 1/2 rn [r/2,2r] ϱ + i=1 V n h, un c t V L Q T ;R 3 T ϱ n 1 5/6 ] rn 1 6/5 1 R+ \[r/2,2r]ϱ V n h,0, un 6 1/6 c t V L Q T ;R 3 E 1/2 ϱ n 1, û n 1 ˆr n 1, T +E 5/6 ϱ n 1, û n 1 ˆr n 1, n 1 ˆV h,0 T V n h,0, un 6 L 6 ;R 3 1/6, ˆV n 1 h,0 where we have used 8.2 and estimate 4.8 to obtain the last line. Now, we write V n h,0, un = [V n h,0 un ] V n h,0 un + V n h,0 un and use the Minkowski inequality together with formulas 2.53, 2.42 to get 1/6 V n h,0, 1/2. un 6 L 6 ;R 3 V n h,0 un 2 L 2 ;R 3 T T Finally, employing Young s inequality, and estimate 4.8, we arrive at P 1 cδ, M 0, E 0, r, r, V, V, t V L Q T,R 15 Eϱ 0, û 0 ˆr 0, ˆV m 0 h,0 + Eϱ n, û n ˆr n, ˆV n m h,0 + δ V n h,0 un 2 L 2 ;R 3, 8.3 T 30
arxiv: v1 [math.na] 26 Aug 2015
Error estimates for a numerical method for the compressible Navier-Stokes system on sufficiently smooth domains Eduard Feireisl Radim Hošek David Maltese Antonín Novotný arxiv:1508.06432v1 [math.na] 26
More informationAsymptotic preserving error estimates for numerical solutions of compressible Navier-Stokes equations in the low Mach number regime
Asymptotic preserving error estimates for numerical solutions of compressible Navier-Stokes equations in the low Mach number regime Eduard Feireisl Mária Lukáčová-Medviďová Šárka Nečasová Antonín Novotný
More informationRelative entropies, suitable weak solutions, and weak-strong uniqueness for the compressible Navier-Stokes system
Relative entropies, suitable weak solutions, and weak-strong uniqueness for the compressible Navier-Stokes system Institute of Mathematics, Academy of Sciences of the Czech Republic, Prague joint work
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES Asymptotic preserving error estimates for numerical solutions of compressible Navier-Stokes equations in the low Mach number regime Eduard Feireisl
More informationWeak-strong uniqueness for the compressible Navier Stokes equations with a hard-sphere pressure law
Weak-strong uniqueness for the compressible Navier Stokes equations with a hard-sphere pressure law Eduard Feireisl Yong Lu Antonín Novotný Institute of Mathematics, Academy of Sciences of the Czech Republic
More informationFinite difference MAC scheme for compressible Navier-Stokes equations
Finite difference MAC scheme for compressible Navier-Stokes equations Bangwei She with R. Hošek, H. Mizerova Workshop on Mathematical Fluid Dynamics, Bad Boll May. 07 May. 11, 2018 Compressible barotropic
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. Stability of strong solutions to the Navier Stokes Fourier system
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES Stability of strong solutions to the Navier Stokes Fourier system Jan Březina Eduard Feireisl Antonín Novotný Preprint No. 7-218 PRAHA 218 Stability
More informationWeak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System
Weak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System Joshua Ballew University of Maryland College Park Applied PDE RIT March 4, 2013 Outline Description of the Model Relative Entropy Weakly
More informationMathematical theory of fluids in motion
Mathematical theory of fluids in motion Eduard Feireisl March 28, 27 Institute of Mathematics of the Academy of Sciences of the Czech Republic Žitná 25, CZ-5 67 Praha, Czech Republic Abstract The goal
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES Local existence of strong solutions and weak-strong uniqueness for the compressible Navier-Stokes system on moving domains Ondřej Kreml Šárka Nečasová
More informationON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM
ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM OLEG ZUBELEVICH DEPARTMENT OF MATHEMATICS THE BUDGET AND TREASURY ACADEMY OF THE MINISTRY OF FINANCE OF THE RUSSIAN FEDERATION 7, ZLATOUSTINSKY MALIY PER.,
More informationINSTITUTE of MATHEMATICS. ACADEMY of SCIENCES of the CZECH REPUBLIC
INSTITUTEofMATHEMATICS Academy of Sciences Czech Republic INSTITUTE of MATHEMATICS ACADEMY of SCIENCES of the CZECH REPUBLIC On weak solutions to a diffuse interface model of a binary mixture of compressible
More informationRETRACTED. INSTITUTE of MATHEMATICS. ACADEMY of SCIENCES of the CZECH REPUBLIC
INSTITUTE of MATHEMATICS ACADEMY of SCIENCES of the CZECH REPULIC Academy of Sciences Czech Republic INSTITUTE of MATHEMATICS Weak solutions to the barotropic Navier-Stokes system with slip with slip boundary
More informationOn weak solution approach to problems in fluid dynamics
On weak solution approach to problems in fluid dynamics Eduard Feireisl based on joint work with J.Březina (Tokio), C.Klingenberg, and S.Markfelder (Wuerzburg), O.Kreml (Praha), M. Lukáčová (Mainz), H.Mizerová
More informationSobolev Spaces. Chapter 10
Chapter 1 Sobolev Spaces We now define spaces H 1,p (R n ), known as Sobolev spaces. For u to belong to H 1,p (R n ), we require that u L p (R n ) and that u have weak derivatives of first order in L p
More informationThe incompressible Navier-Stokes equations in vacuum
The incompressible, Université Paris-Est Créteil Piotr Bogus law Mucha, Warsaw University Journées Jeunes EDPistes 218, Institut Elie Cartan, Université de Lorraine March 23th, 218 Incompressible Navier-Stokes
More informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationA regularity criterion for the weak solutions to the Navier-Stokes-Fourier system
A regularity criterion for the weak solutions to the Navier-Stokes-Fourier system Eduard Feireisl Antonín Novotný Yongzhong Sun Charles University in Prague, Faculty of Mathematics and Physics, Mathematical
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationarxiv: v1 [math.ap] 29 May 2018
Non-uniqueness of admissible weak solution to the Riemann problem for the full Euler system in D arxiv:805.354v [math.ap] 9 May 08 Hind Al Baba Christian Klingenberg Ondřej Kreml Václav Mácha Simon Markfelder
More informationat time t, in dimension d. The index i varies in a countable set I. We call configuration the family, denoted generically by Φ: U (x i (t) x j (t))
Notations In this chapter we investigate infinite systems of interacting particles subject to Newtonian dynamics Each particle is characterized by its position an velocity x i t, v i t R d R d at time
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. Homogenization of a non homogeneous heat conducting fluid
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES Homogenization of a non homogeneous heat conducting fluid Eduard Feireisl Yong Lu Yongzhong Sun Preprint No. 16-2018 PRAHA 2018 Homogenization of
More informationHomogenization of stationary Navier-Stokes equations in domains with tiny holes
Homogenization of stationary Navier-Stokes equations in domains with tiny holes Eduard Feireisl Yong Lu Institute of Mathematics of the Academy of Sciences of the Czech Republic Žitná 25, 115 67 Praha
More informationHomogenization of the compressible Navier Stokes equations in domains with very tiny holes
Homogenization of the compressible Navier Stokes equations in domains with very tiny holes Yong Lu Sebastian Schwarzacher Abstract We consider the homogenization problem of the compressible Navier Stokes
More informationConvergence of the MAC scheme for incompressible flows
Convergence of the MAC scheme for incompressible flows T. Galloue t?, R. Herbin?, J.-C. Latche??, K. Mallem????????? Aix-Marseille Universite I.R.S.N. Cadarache Universite d Annaba Calanque de Cortiou
More informationOn Smoothness of Suitable Weak Solutions to the Navier-Stokes Equations
On Smoothness of Suitable Weak Solutions to the Navier-Stokes Equations G. Seregin, V. Šverák Dedicated to Vsevolod Alexeevich Solonnikov Abstract We prove two sufficient conditions for local regularity
More informationKonstantinos Chrysafinos 1 and L. Steven Hou Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANALYSIS AND APPROXIMATIONS OF THE EVOLUTIONARY STOKES EQUATIONS WITH INHOMOGENEOUS
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationON THE STRONG SOLUTIONS OF THE INHOMOGENEOUS INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN A THIN DOMAIN
ON THE STRONG SOLUTIONS OF THE INHOMOGENEOUS INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN A THIN DOMAIN XIAN LIAO Abstract. In this work we will show the global existence of the strong solutions of the inhomogeneous
More informationTHE L 2 -HODGE THEORY AND REPRESENTATION ON R n
THE L 2 -HODGE THEORY AND REPRESENTATION ON R n BAISHENG YAN Abstract. We present an elementary L 2 -Hodge theory on whole R n based on the minimization principle of the calculus of variations and some
More informationMaximal dissipation principle for the complete Euler system
Maximal dissipation principle for the complete Euler system Jan Březina Eduard Feireisl December 13, 217 Tokyo Institute of Technology 2-12-1 Ookayama, Meguro-ku, Tokyo, 152-855, Japan Institute of Mathematics
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
More informationSingularities and Laplacians in Boundary Value Problems for Nonlinear Ordinary Differential Equations
Singularities and Laplacians in Boundary Value Problems for Nonlinear Ordinary Differential Equations Irena Rachůnková, Svatoslav Staněk, Department of Mathematics, Palacký University, 779 OLOMOUC, Tomkova
More informationMath The Laplacian. 1 Green s Identities, Fundamental Solution
Math. 209 The Laplacian Green s Identities, Fundamental Solution Let be a bounded open set in R n, n 2, with smooth boundary. The fact that the boundary is smooth means that at each point x the external
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationMetric Spaces and Topology
Chapter 2 Metric Spaces and Topology From an engineering perspective, the most important way to construct a topology on a set is to define the topology in terms of a metric on the set. This approach underlies
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationJUHA KINNUNEN. Harmonic Analysis
JUHA KINNUNEN Harmonic Analysis Department of Mathematics and Systems Analysis, Aalto University 27 Contents Calderón-Zygmund decomposition. Dyadic subcubes of a cube.........................2 Dyadic cubes
More informationOn the domain dependence of solutions to the compressible Navier-Stokes equations of an isothermal fluid
Nečas Center for Mathematical Modeling On the domain dependence of solutions to the compressible Navier-Stokes equations of an isothermal fluid Nikola Hlaváčová Preprint no. 211-2 Research Team 1 Mathematical
More informationThe 2-d isentropic compressible Euler equations may have infinitely many solutions which conserve energy
The -d isentropic compressible Euler equations may have infinitely many solutions which conserve energy Simon Markfelder Christian Klingenberg September 15, 017 Dept. of Mathematics, Würzburg University,
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. Dissipative measure-valued solutions to the compressible Navier-Stokes system
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES Dissipative measure-valued solutions to the compressible Navier-Stokes system Eduard Feireisl Piotr Gwiazda Agnieszka Świerczewska-Gwiazda Emil Wiedemann
More informationABOUT STEADY TRANSPORT EQUATION II SCHAUDER ESTIMATES IN DOMAINS WITH SMOOTH BOUNDARIES
PORTUGALIAE MATHEMATICA Vol. 54 Fasc. 3 1997 ABOUT STEADY TRANSPORT EQUATION II SCHAUDER ESTIMATES IN DOMAINS WITH SMOOTH BOUNDARIES Antonin Novotny Presented by Hugo Beirão da Veiga Abstract: This paper
More informationLiquid crystal flows in two dimensions
Liquid crystal flows in two dimensions Fanghua Lin Junyu Lin Changyou Wang Abstract The paper is concerned with a simplified hydrodynamic equation, proposed by Ericksen and Leslie, modeling the flow of
More informationRelation between Distributional and Leray-Hopf Solutions to the Navier-Stokes Equations
Relation between Distributional and Leray-Hopf Solutions to the Navier-Stokes Equations Giovanni P. Galdi Department of Mechanical Engineering & Materials Science and Department of Mathematics University
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES Weak solutions to the full Navier-Stokes-Fourier system with slip boundary conditions in time dependent domains Ondřej Kreml Václav Mácha Šárka Nečasová
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationSobolev spaces. May 18
Sobolev spaces May 18 2015 1 Weak derivatives The purpose of these notes is to give a very basic introduction to Sobolev spaces. More extensive treatments can e.g. be found in the classical references
More informationarxiv: v1 [math.na] 16 Oct 2007
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher arxiv:070.2987v [math.na] 6 Oct 2007 AN UNCONDITIONNALLY STABLE PRESSURE CORRECTION
More informationUniform estimates for Stokes equations in domains with small holes and applications in homogenization problems
Uniform estimates for Stokes equations in domains with small holes and applications in homogenization problems Yong Lu Abstract We consider the Dirichlet problem for the Stokes equations in a domain with
More informationOn the local well-posedness of compressible viscous flows with bounded density
On the local well-posedness of compressible viscous flows with bounded density Marius Paicu University of Bordeaux joint work with Raphaël Danchin and Francesco Fanelli Mathflows 2018, Porquerolles September
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
More informationRemarks on Bronštein s root theorem
Remarks on Bronštein s root theorem Guy Métivier January 23, 2017 1 Introduction In [Br1], M.D.Bronštein proved that the roots of hyperbolic polynomials (1.1) p(t, τ) = τ m + m p k (t)τ m k. which depend
More informationApproximation of Minimal Functions by Extreme Functions
Approximation of Minimal Functions by Extreme Functions Teresa M. Lebair and Amitabh Basu August 14, 2017 Abstract In a recent paper, Basu, Hildebrand, and Molinaro established that the set of continuous
More informationTHE STOKES SYSTEM R.E. SHOWALTER
THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationIntroduction to Real Analysis Alternative Chapter 1
Christopher Heil Introduction to Real Analysis Alternative Chapter 1 A Primer on Norms and Banach Spaces Last Updated: March 10, 2018 c 2018 by Christopher Heil Chapter 1 A Primer on Norms and Banach Spaces
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationREMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID
REMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID DRAGOŞ IFTIMIE AND JAMES P. KELLIHER Abstract. In [Math. Ann. 336 (2006), 449-489] the authors consider the two dimensional
More informationA generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem
A generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem Dave McCormick joint work with James Robinson and José Rodrigo Mathematics and Statistics Centre for Doctoral Training University
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationL p Functions. Given a measure space (X, µ) and a real number p [1, ), recall that the L p -norm of a measurable function f : X R is defined by
L p Functions Given a measure space (, µ) and a real number p [, ), recall that the L p -norm of a measurable function f : R is defined by f p = ( ) /p f p dµ Note that the L p -norm of a function f may
More informationConsistency analysis of a 1D Finite Volume scheme for barotropic Euler models
Consistency analysis of a 1D Finite Volume scheme for barotropic Euler models F Berthelin 1,3, T Goudon 1,3, and S Mineaud,3 1 Inria, Sophia Antipolis Méditerranée Research Centre, Proect COFFEE Inria,
More informationOn a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability Boundary Conditions
Proceedings of the 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August -, 5 pp36-41 On a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability
More informationOn incompressible limits for the Navier-Stokes system on unbounded domains under slip boundary conditions
Nečas Center for Mathematical Modeling On incompressible limits for the Navier-Stokes system on unbounded domains under slip boundary conditions Donatella Donatelli and Eduard Feireisl and Antonín Novotný
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationOPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS
PORTUGALIAE MATHEMATICA Vol. 59 Fasc. 2 2002 Nova Série OPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS J. Saint Jean Paulin and H. Zoubairi Abstract: We study a problem of
More informationLECTURE 15: COMPLETENESS AND CONVEXITY
LECTURE 15: COMPLETENESS AND CONVEXITY 1. The Hopf-Rinow Theorem Recall that a Riemannian manifold (M, g) is called geodesically complete if the maximal defining interval of any geodesic is R. On the other
More informationESTIMATES FOR THE MONGE-AMPERE EQUATION
GLOBAL W 2,p ESTIMATES FOR THE MONGE-AMPERE EQUATION O. SAVIN Abstract. We use a localization property of boundary sections for solutions to the Monge-Ampere equation obtain global W 2,p estimates under
More informationMath212a1413 The Lebesgue integral.
Math212a1413 The Lebesgue integral. October 28, 2014 Simple functions. In what follows, (X, F, m) is a space with a σ-field of sets, and m a measure on F. The purpose of today s lecture is to develop the
More informationOn Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1
On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1 Jie Shen Department of Mathematics, Penn State University University Park, PA 1682 Abstract. We present some
More informationChapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma
Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions
More informationSufficient conditions on Liouville type theorems for the 3D steady Navier-Stokes equations
arxiv:1805.07v1 [math.ap] 6 May 018 Sufficient conditions on Liouville type theorems for the D steady Navier-Stokes euations G. Seregin, W. Wang May 8, 018 Abstract Our aim is to prove Liouville type theorems
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationFinite Element Methods for Maxwell Equations
CHAPTER 8 Finite Element Methods for Maxwell Equations The Maxwell equations comprise four first-order partial differential equations linking the fundamental electromagnetic quantities, the electric field
More informationNATIONAL BOARD FOR HIGHER MATHEMATICS. Research Scholarships Screening Test. Saturday, January 20, Time Allowed: 150 Minutes Maximum Marks: 40
NATIONAL BOARD FOR HIGHER MATHEMATICS Research Scholarships Screening Test Saturday, January 2, 218 Time Allowed: 15 Minutes Maximum Marks: 4 Please read, carefully, the instructions that follow. INSTRUCTIONS
More informationINCOMPRESSIBLE FLUIDS IN THIN DOMAINS WITH NAVIER FRICTION BOUNDARY CONDITIONS (II) Luan Thach Hoang. IMA Preprint Series #2406.
INCOMPRESSIBLE FLUIDS IN THIN DOMAINS WITH NAVIER FRICTION BOUNDARY CONDITIONS II By Luan Thach Hoang IMA Preprint Series #2406 August 2012 INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS UNIVERSITY OF
More informationA LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE. 1.
A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE THOMAS CHEN AND NATAŠA PAVLOVIĆ Abstract. We prove a Beale-Kato-Majda criterion
More informationLecture Notes: African Institute of Mathematics Senegal, January Topic Title: A short introduction to numerical methods for elliptic PDEs
Lecture Notes: African Institute of Mathematics Senegal, January 26 opic itle: A short introduction to numerical methods for elliptic PDEs Authors and Lecturers: Gerard Awanou (University of Illinois-Chicago)
More informationMathematics of complete fluid systems
Mathematics of complete fluid systems Eduard Feireisl feireisl@math.cas.cz Institute of Mathematics of the Academy of Sciences of the Czech Republic Žitná 25, 115 67 Praha 1, Czech Republic The scientist
More informationGlobal well-posedness and decay for the viscous surface wave problem without surface tension
Global well-posedness and decay for the viscous surface wave problem without surface tension Ian Tice (joint work with Yan Guo) Université Paris-Est Créteil Laboratoire d Analyse et de Mathématiques Appliquées
More informationTHE WAVE EQUATION. d = 1: D Alembert s formula We begin with the initial value problem in 1 space dimension { u = utt u xx = 0, in (0, ) R, (2)
THE WAVE EQUATION () The free wave equation takes the form u := ( t x )u = 0, u : R t R d x R In the literature, the operator := t x is called the D Alembertian on R +d. Later we shall also consider the
More informationGLOBAL LIPSCHITZ CONTINUITY FOR MINIMA OF DEGENERATE PROBLEMS
GLOBAL LIPSCHITZ CONTINUITY FOR MINIMA OF DEGENERATE PROBLEMS PIERRE BOUSQUET AND LORENZO BRASCO Abstract. We consider the problem of minimizing the Lagrangian [F ( u+f u among functions on R N with given
More informationThe oblique derivative problem for general elliptic systems in Lipschitz domains
M. MITREA The oblique derivative problem for general elliptic systems in Lipschitz domains Let M be a smooth, oriented, connected, compact, boundaryless manifold of real dimension m, and let T M and T
More informationFormulation of the problem
TOPICAL PROBLEMS OF FLUID MECHANICS DOI: https://doi.org/.43/tpfm.27. NOTE ON THE PROBLEM OF DISSIPATIVE MEASURE-VALUED SOLUTIONS TO THE COMPRESSIBLE NON-NEWTONIAN SYSTEM H. Al Baba, 2, M. Caggio, B. Ducomet
More informationINTRODUCTION TO REAL ANALYTIC GEOMETRY
INTRODUCTION TO REAL ANALYTIC GEOMETRY KRZYSZTOF KURDYKA 1. Analytic functions in several variables 1.1. Summable families. Let (E, ) be a normed space over the field R or C, dim E
More informationVariational Theory of Solitons for a Higher Order Generalized Camassa-Holm Equation
International Journal of Mathematical Analysis Vol. 11, 2017, no. 21, 1007-1018 HIKAI Ltd, www.m-hikari.com https://doi.org/10.12988/ijma.2017.710141 Variational Theory of Solitons for a Higher Order Generalized
More informationChapter One. The Calderón-Zygmund Theory I: Ellipticity
Chapter One The Calderón-Zygmund Theory I: Ellipticity Our story begins with a classical situation: convolution with homogeneous, Calderón- Zygmund ( kernels on R n. Let S n 1 R n denote the unit sphere
More informationA CONVEX-CONCAVE ELLIPTIC PROBLEM WITH A PARAMETER ON THE BOUNDARY CONDITION
A CONVEX-CONCAVE ELLIPTIC PROBLEM WITH A PARAMETER ON THE BOUNDARY CONDITION JORGE GARCÍA-MELIÁN, JULIO D. ROSSI AND JOSÉ C. SABINA DE LIS Abstract. In this paper we study existence and multiplicity of
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
More informationb i (µ, x, s) ei ϕ(x) µ s (dx) ds (2) i=1
NONLINEAR EVOLTION EQATIONS FOR MEASRES ON INFINITE DIMENSIONAL SPACES V.I. Bogachev 1, G. Da Prato 2, M. Röckner 3, S.V. Shaposhnikov 1 The goal of this work is to prove the existence of a solution to
More informationTools from Lebesgue integration
Tools from Lebesgue integration E.P. van den Ban Fall 2005 Introduction In these notes we describe some of the basic tools from the theory of Lebesgue integration. Definitions and results will be given
More informationEXISTENCE AND REGULARITY RESULTS FOR SOME NONLINEAR PARABOLIC EQUATIONS
EXISTECE AD REGULARITY RESULTS FOR SOME OLIEAR PARABOLIC EUATIOS Lucio BOCCARDO 1 Andrea DALL AGLIO 2 Thierry GALLOUËT3 Luigi ORSIA 1 Abstract We prove summability results for the solutions of nonlinear
More informationNECESSARY CONDITIONS FOR WEIGHTED POINTWISE HARDY INEQUALITIES
NECESSARY CONDITIONS FOR WEIGHTED POINTWISE HARDY INEQUALITIES JUHA LEHRBÄCK Abstract. We establish necessary conditions for domains Ω R n which admit the pointwise (p, β)-hardy inequality u(x) Cd Ω(x)
More informationDETERMINATION OF THE BLOW-UP RATE FOR THE SEMILINEAR WAVE EQUATION
DETERMINATION OF THE LOW-UP RATE FOR THE SEMILINEAR WAVE EQUATION y FRANK MERLE and HATEM ZAAG Abstract. In this paper, we find the optimal blow-up rate for the semilinear wave equation with a power nonlinearity.
More informationFriedrich symmetric systems
viii CHAPTER 8 Friedrich symmetric systems In this chapter, we describe a theory due to Friedrich [13] for positive symmetric systems, which gives the existence and uniqueness of weak solutions of boundary
More informationASYMPTOTICALLY NONEXPANSIVE MAPPINGS IN MODULAR FUNCTION SPACES ABSTRACT
ASYMPTOTICALLY NONEXPANSIVE MAPPINGS IN MODULAR FUNCTION SPACES T. DOMINGUEZ-BENAVIDES, M.A. KHAMSI AND S. SAMADI ABSTRACT In this paper, we prove that if ρ is a convex, σ-finite modular function satisfying
More informationL p -boundedness of the Hilbert transform
L p -boundedness of the Hilbert transform Kunal Narayan Chaudhury Abstract The Hilbert transform is essentially the only singular operator in one dimension. This undoubtedly makes it one of the the most
More informationPropagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R
Propagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455 Abstract We consider semilinear
More informationAnalysis Finite and Infinite Sets The Real Numbers The Cantor Set
Analysis Finite and Infinite Sets Definition. An initial segment is {n N n n 0 }. Definition. A finite set can be put into one-to-one correspondence with an initial segment. The empty set is also considered
More informationOn some nonlinear parabolic equation involving variable exponents
On some nonlinear parabolic equation involving variable exponents Goro Akagi (Kobe University, Japan) Based on a joint work with Giulio Schimperna (Pavia Univ., Italy) Workshop DIMO-2013 Diffuse Interface
More information