arxiv: v1 [math.na] 12 Mar 2018

Transcription

1 ON PRESSURE ESTIMATES FOR THE NAVIER-STOKES EQUATIONS J A FIORDILINO arxiv: v1 [matna 12 Mar 2018 Abstract Tis paper presents a simple, general tecnique to prove finite element metod (FEM) pressure stability and convergence Typically, pressure estimates are ignored in te literature However, full reliability of a numerical metod is not establised unless te solution pair (u,p) is treated Te simplicity of te proposed tecnique puts pressure estimates witin reac of many existing and future numerical metods and lends itself to te numerical analyst s toolbox 1 Introduction In tis paper, we illustrate a simple, general tecnique for proving stability and convergence of FEM pressure approximations for te incompressible Navier-Stokes equations (NSE) Pressure estimates are acievable but typically ignored due to te difficulty in obtaining tem; see, eg, Heywood and Rannacer [7, Ingram [8, and Labovscii [9 However, a complete numerical analysis establises full reliability of a numerical metod Terefore, a complete treatment sould be given to te velocity-pressure solution pair (u, p) wenever possible Consequently, scientists and engineers will be more confident in an algoritm s capabilities and limitations Te proposed metod utilizes te discrete inf-sup condition, (6), and equivalence of te dual norms X and V, Lemma 3 in Section 21 Regarding te latter, tis equivalence was establised by Galvin in [5 Galvin used tis result to prove pressurestability in L 2 (0,t ;L 2 (Ω)) by bounding te discrete time derivative of te velocity approximation Later, Zang, Hou, and Zao [14 followed te stability proof to produce an error estimate for te pressure in L 1 (0,t ;L 2 (Ω)) However, as is common in intricate, tecnical analyses, teir result was correct but a gap was present in teir analysis Herein, we illustrate te tecnique for linearly implicit Backward Euler (LIBE), (7) It will be easy to see tat our analysis can be extended to many timestepping metods including linearly implicit BDF-k and Crank-Nicolson variants Alternatively, for penalty and artificial compressibility metods, were te velocity approximation is not weakly divergence free, weaker forms of stability and convergence can be proven in te fully discrete setting [4; see, eg, Sen [12 for semi-discrete analyses Recall te NSE Let Ω R d (d = 2,3) be a convex polyedral domain wit piecewise smoot boundary Ω Given te fluid viscosity ν, u(x,0) = u 0 (x), and te body force f(x,t), te velocity u(x,t) : Ω (0,t R d and pressure p(x,t) : Ω (0,t R satisfy (1) u t +u u ν u+ p = f in Ω, u = 0 in Ω, u = 0 on Ω In Section 2, we collect necessary matematical tools In Section 3, we present te fully discrete numerical sceme, (7), tat will be analyzed A complete stability and error analysis are presented in Sections 4 and 5 In particular, te velocity and pressure approximations are proven to be unconditionally, nonlinearly, energy stable in Teorem 4 and Corollary 5 First-order, optimal convergence is proven in Teorem 6 and Corollary 7 Lastly, we state conclusions in Section 6 Supported by te DoD SMART Scolarsip Partially supported by NSF grants CBET and DMS

2 2 Matematical Preliminaries Te L 2 (Ω) inner product is (, ) and te induced norm is Te L (Ω) norm is denoted Furter, H s (Ω) denotes te Hilbert spaces of L 2 (Ω) functions wit distributional derivatives of order s 0 in L 2 (Ω) Te corresponding norms and seminorms are s and s ; note tat 0 = 0 = Define te Hilbert spaces, X := H 1 0 (Ω)d = {v H 1 (Ω) d : v = 0 on Ω}, Q := L 2 0 (Ω) = {q L2 (Ω) : (1,q) = 0}, V := {v X : (q, v) = 0 q Q} For functions v X, te Poincaré-Friedrics inequality olds, v C p v Te explicitly skew-symmetric trilinear form is denoted: It enjoys te following properties b(u,v,w) = 1 2 (u v,w) 1 (u w,v) u,v,w X 2 Lemma 1 Tere exists C 1 and C 2 suc tat for all u,v,w X, b(u,v,w) satisfies b(u,v,w) = (u v,w)+ 1 2 (( u)v,w), b(u,v,w) C 1 u v w, b(u,v,w) C 2 u u v w Proof See Lemma 21 on p 12 of [13 Te discrete time analysis will utilize te following norms, for 1 k < : v,k := max 0 n N vn k, v p,k := ( N v n p 1/p k) Te weak formulation of system (1) is: Find u : [0,t X and p : [0,t Q for ae t (0,t satisfying (2) (3) (u t,v)+b(u, u,v)+ν( u, v) (p, v) = (f,v) v X, (q, u) = 0 q Q 21 Finite Element Preliminaries Consider a quasi-uniform mes Ω = {K} of Ω wit maximum triangle diameter lengt Let X X and Q Q be conforming finite element spaces consisting of continuous piecewise polynomials of degrees j and l, respectively Moreover, tey satisfy te following approximation properties [3, 1 j,l k,m: { } (4) inf u v + (u v ) C k+1 u k+1, v X (5) inf q Q p q C m p m, 2

3 for all u X H k+1 (Ω) d and p Q H m (Ω) Furtermore, we consider tose spaces for wic te discrete inf-sup condition is satisfied, (6) inf q Q sup v X (q, v ) q v α > 0, were α is independent of Examples include te MINI-element and Taylor-Hood family of elements [10 Te space of discretely divergence free functions is defined by V := {v X : (q, v ) = 0, q Q } Te discrete inf-sup condition implies tat we may approximate functions in V well by functions in V, Lemma 2 Suppose te discrete inf-sup condition (6) olds, ten for any v V inf (v v ) C(α) inf (v v ) v V v X Proof See Capter 2, Teorem 11 on p 59 of [6 Te spaces X and V, dual to X and V, are endowed wit te following dual norms w X := sup v X (w,v ) v, w V Interestingly, tese norms are equivalent for functions in V := sup v V (w,v ) v Lemma 3 Let w V Ten, tere exists C > 0, independent of, suc tat Proof See Lemma 1 on p 243 of [5 C w X w V w X Remark: Te proof of te above lemma utilizes te H 1 (Ω) d stability of te L 2 (Ω) d ortogonal projection onto V Tis result can be extended to regular meses under certain conditions on te mes; see [1, 2 and references terein 3 Numerical Sceme Denote te fully discrete solutions by u n and pn at time levels tn = n, n = 0,1,,N, and t = N Te fully discrete approximations of (1) are LIBE: Given u n X, find (u n+1 ) (X,Q ) satisfying,p n+1 (7) ( un+1 u n,v )+b(u n,un+1,v )+ν(u n+1,v ) (p n+1, v ) = (f n+1,v ) v X, ( u n+1,q ) q Q 4 Stability Analysis In tis section, we prove tat te pressure approximation is stable in L 1 (0,t ;L 2 (Ω)) and L 2 (0,t ;L 2 (Ω)) We first state, witout proof, Teorem 4 regarding te stability of te velocity approximation Stability of te pressure approximation follows from tis teorem as a proven corollary 3

4 Teorem 4 Consider te numerical sceme (7) Suppose f L 2 (0,t ;H 1 (Ω) d ) Ten, te velocity approximation satisfies for all N 1 u N 2 + u n+1 u n 2 +ν u 2 2,0 ν 1 f 2 2, 1 + u0 2 As a consequence of te above, te following olds Corollary 5 Suppose Teorem 4 olds Ten, te pressure approximation satisfies for all N 1 α [ (C1 p n+1 (1+C 1 ) ν 2 f 2, 1 +2 t ) f 2, 1 + ( C 1 ν 1 u 0 + νt ) u 0, and α p 2,0 [ 3(1+C 1 ) C1 2C2 ν 2 +1 f 2, 1 + C1 2C2 ν 1 +ν u 0 Proof Consider equation (7) Let (v,q ) (V,Q ) and isolate te discrete time derivative Ten, (8) ( un+1 u n,v ) = b(u n,u n+1,v ) ν(u n+1,v )+(f n+1,v ) v V Te terms on te rigt-and side can be bounded using Lemma 1, te Caucy-Scwarz inequality, and duality, respectively, (9) (10) (11) b(u n,un+1,v ) C 1 u n un+1 v, ν( u n+1, v ) ν u n+1 v, (f n+1,v ) f n+1 1 v Using te above estimates in equation (8), dividing bot sides by v, and taking te supremum over v V yields (12) Lemma 3 ten implies un+1 u n V ( C 1 u n +ν ) u n+1 + f n+1 1 (13) un+1 u n [ (C1 X C 1 u n +ν) u n+1 + f n+1 1 Now, reconsider equation (7) wit v X Isolate te pressure term and use te estimates (9) - (11) Ten, (14) (p n+1, v ) ( un+1 u n,v )+ ( C 1 u n +ν) u n+1 v + f n+1 1 v 4

5 Divide bot sides by v, take te supremum over v X, and use bot te discrete inf-sup condition 6 and estimate (13) Ten, (15) [ (C1 α p n+1 (1+C 1 ) u n +ν ) u n+1 + f n+1 1 Multiplying by, summing from n = 0 to n = N 1, and using te Caucy-Scwarz inequality on te rigt-and side terms yields α [ (C1 p n+1 (1+C 1 ) u 2,0 +ν t ) u 2,0 + t f 2, 1 Lastly, application of Teorem 4 on te velocity terms gives (16) α (C1 p n+1 (1+C )[ 1 ν 2 f +2 2, 1 t ) f +( 2, 1 C 1 ν 1 u 0 + νt ) u 0 For te second result, reconsider equation (15) Square bot sides and notice tat te rigt-and side is te square of te l 1 norm of te 3-vector χ = (C 1 u n un+1,ν u n, fn+1 1 ) T Consequently, using χ 2 l 1 3 χ 2 l 2, multiplying by, and summing from n = 0 to n = N 1 we ave α 2 p 2 2,0 3(1+C 1 )2[ C1 2 u n 2 u n+1 2 +ν 2 u 2 2,0 + f 2 2, 1 Now, u n X for every 0 n N 1; tat is, max 0 n u n C < Tus, using Teorem 4 yields α 2 p 2 2,0 3(1+C 1 ) 2[ ( C 2 1 C 2 +ν 2) (17) u 2 2,0 + f 2 2, 1 3(1+C 1 )2[ (C 2 1 C2 ν 2 +1) f 2 2, 1 +( C 2 1 C2 ν 1 +ν ) u 0 2 Te result follows after taking te square root In te above, note tat growt wit respect to t is allowed in L 1 (0,t ;L 2 (Ω)) Alternatively, an unknown, bounded linear growt factor C is present wit respect to L 2 (0,t ;L 2 (Ω)) Clearly, te velocity approximation is most stable; tat is, stable uniformly in time in te appropriate norms 5 Error Analysis Herein, we prove optimal-order convergence of te pressure approximation in L 1 (0,t ;L 2 (Ω)) and L 2 (0,t ;L 2 (Ω)) Te outline of tis section is identical to te last We first state a convergence teorem for te velocity approximation, Teorem 6, below We ten prove convergence for te pressure approximation as a corollary Denote u n and p n as te true solutions at time t n = n Assume te solutions satisfy te following regularity assumptions: (18) (19) u L (0,t ;X H k+1 (Ω)), u t L 2 (0,t ;H k+1 (Ω)), u tt L 2 (0,t ;H k+1 (Ω)), p L 2 (0,t ;Q H m (Ω)) 5

6 Remark: Tese assumptions are consistent wit wat is seen in te literature; more commonly, increased regularity is demanded Te errors for te solution variables are denoted e n u = u n u n, e n p = p n p n We begin again wit a result for te velocity approximation: optimal convergence Teorem 6 For u and p satisfying (1), suppose tat u 0 X is an approximation of u 0 to witin te accuracy of te interpolant Ten tere exists constants C, C > 0 suc tat LIBE satisfies ( ) e N u 2 + e n+1 û e n u 2 +ν eû 2 2,0 Cexp(C t ) 2k + 2 +initial errors As a consequence, we ave: Corollary 7 Suppose Teorem 6 and its ypoteses old Ten tere exists a constant C > 0 suc tat LIBE satisfies ( α e n+1 p + e p 2,0 ) Cexp( C t 2 ( ) ) k ++initial errors Proof Recall te error equation for LIBE, (20) ( en+1 u e n u,v )+ν( e n+1 u, v )+b(u n+1 u n,u n+1,v )+b(e n u,u n+1,v ) +b(u n,e n+1 u,v ) (e n+1 p, v ) = ( un+1 u n u n+1 t,v ) v X Let v V and rewrite 1 (e n+1 u e n u,v ) = 1 (η n+1 η n,v ) 1 (φ n+1 φ n,v ) were e n u = (un I (u n )) (u n I (u n )) = η n φ n I (u n ) is an interpolant, suc as te Lagrange interpolant,ofuintotefiniteelementspacesuctatφ n V Moreover,notetat(e n+1 p, v ) = (p n+1 q n+1, v ), since v V, and rearrange Ten, (21) ( φn+1 φ n,v ) = ( ηn+1 η n,v )+ν( e n+1 u, v )+b(u n+1 u n,u n+1,v )+b(e n u,un+1,v ) +b(u n,en+1 u,v ) (p n+1 q n+1, v ) ( un+1 u n u n+1 t,v ) v V 6

7 Using standard estimates on te rigt-and side terms yields (22) (23) (24) (25) (26) (27) (28) ( ηn+1 η n,v ) C p 1/2 η t L 2 (t n,t n+1 ;L 2 (Ω)) v, ν( e n+1 u, v ) ν e n+1 v, u b(u n+1 u n,u n+1,v ) C 1 1/2 u t L 2 (t n,t n+1 ;L 2 (Ω)) u n+1 v, b(e n u,un+1,v ) C 1 e n u un+1 v, b(u n,en+1 u,v ) C 1 u n+1 e n+1 u v, (p n+1 q n+1, v ) d p n+1 q n+1 v, ( un+1 u n u n+1 t,v ) C 1/2 u tt L2 (t n,t n+1 ;L 2 (Ω)) v Using te above estimates in equation (21), dividing bot sides by v, taking a supremum over V and using Lemma 3 yields (29) φn+1 φ n [ X C 1 C p 1/2 η t L2 (t n,t n+1 ;L 2 (Ω)) + ( ν +C 1 u n ) e n+1 u +C 1 u n+1 ( e n u + 1/2 u t L2 (t n,t n+1 ;L 2 (Ω)) ) + d p n+1 q n+1 +C 1/2 u tt L2 (t n,t n+1 ;L 2 (Ω)) Reconsidering te error equation (20), splitting te pressure error term via (e n+1 p, v ) = (p n+1 q n+1, v ) (p n+1 q n+1, v ), and rearranging yields (30) (q n+1 φ n, v ) = ( ηn+1 η n,v ) ( φn+1,v )+ν( e n+1 u, v ) +b(u n+1 u n,u n+1,v )+b(e n u,u n+1,v )+b(u n,e n+1 u,v ) (p n+1 q n+1, v ) ( un+1 u n u n+1 t,v ) v X Using te above estimates (22) - (28), dividing by v, taking a supremum over v X, using (29), and te discrete inf-sup condition yields [ (31) α q n+1 (1+C 1 ) C p 1/2 η t L 2 (t n,t n+1 ;L 2 (Ω)) + ( ν +C 1 u n ) e n+1 u +C 1 u n+1 ( e n u + 1/2 ) u t L2 (t n,t n+1 ;L 2 (Ω)) + d p n+1 q n+1 +C 1/2 u tt L2 (t n,t n+1 ;L 2 (Ω)) Multiplying by, summing from n = 0 to n = N 1, and using te Caucy-Scwarz inequality yields (32) α q n+1 [ (1+C 1 C p 1/2 η t L 2 (0,t ;L 2 (Ω)) ) + ( ν ) t +C 1 u 2,0 eu 2,0 +C 1 u 2,0 ( ) e u 2,0 + u t L2 (0,t ;L 2 (Ω)) + dt p n+1 q n+1 2,0 +C 3/2 u tt L 2 (0,t ;L 2 (Ω)) 7

8 By te triangle inequality and l1 N l2, for an N-vector, Consequently, α e n+1 p α p n+1 q n+1 +α α t p n+1 q n+1 2,0 +α q n+1 q n+1 (33) α e n+1 p (1+C [C 1 ) p 1/2 η t L2 (0,t ;L 2 (Ω)) + ( ν ) t +C 1 u 2,0 eu 2,0 +C 1 u 2,0 ( ) e u 2,0 + u t L2 (0,t ;L 2 (Ω)) +(α+ d) t p n+1 q n+1 2,0 +C 3/2 u tt L 2 (0,t ;L 2 (Ω)) Now, consider (31), square bot sides, and use 2 l l 2, for a 7-vector Multiplying by and summing over n from n = 0 to n = N 1 yields α 2 p q 2 2,0 7(1+C 1 ) 2[ C 2 p η t 2 L 2 (0,t ;L 2 (Ω)) +( ν 2 +C 2 1C 2 ) eu 2 2,0 +C 2 1 u 2,0( eu 2 2,0 +2 u t 2 L 2 (0,t ;L 2 (Ω))) +d p q 2 2,0 +C2 2 u tt 2 L 2 (0,t ;L 2 (Ω)) Te triangle inequality gives (34) α 2 e p 2 2,0 α2 p q 2 2,0 +α2 p q 2 2,0 7(1+C 1 )2[ C 2 p η t 2 L 2 (0,t ;L 2 (Ω)) +( ν 2 +C 2 1 C2 ) eu 2 2,0 +C1 2,0( u 2 eu 2 2,0 +2 u t 2 L 2 (0,t ;L (Ω))) +(d+α 2 ) p q 2 2 2,0 +C 2 2 u tt 2 L 2 (0,t ;L 2 (Ω)) Square root (34) and add it to (33) Taking infimums over V and Q and applying Lemma 2, te approximation properties (4) - (5), and Teorem 6 yield te result 6 Conclusion We presented a simple, general tecnique to prove stability and convergence of FEM pressure approximations to te NSE Te tecnique required tat: te discrete inf-sup condition olds and te equivalence of certain dual norms Suc requirements are satisfied by te popular MINI-element and Taylor-Hood family of elements Te tecnique was illustrated on linearly implicit Backward Euler Consequently, te metod is seen to be applicable to many oter numerical scemes Acknowledgements Te autor would like to tank Professor Layton for is encouragement and suggestions Appendix For completeness, we provide a proof of te H 1 (Ω) d stability of te L 2 (Ω) d ortogonal projection onto V for quasi-uniform meses 8

9 Teorem 8 Let Ω Ω be a quasi-uniform mes and X X a conforming finite element space consisting of continuous piecewise polynomials of degree j 1 Denote P : L 2 (Ω) d V as te L 2 -ortogonal projection onto V X satisfying for all u L 2 (Ω) d : Ten, for all u H 1 (Ω) d : (Pu u,v ) = 0 v V Pu C u Proof Fix u H 1 (Ω) d and let I SZ : H 1 (Ω) d X denote te Scott-Zang interpolant [11 Consider Pu, add and subtract I SZ u, and apply te triangle inequality Ten, Pu (Pu I SZ u) + I SZ u ( ) SincePu I SZ u X, teinverseestimate v C inv 1 v olds Consequently, applying te inverse estimate to te first term and H 1 (Ω) d stability of te Scott-Zang interpolant [3 to te second yields (Pu I SZ u) + I SZ u C inv 1 Pu I SZ u +C SZ u Te triangle inequality and interpolation estimates give C inv 1 Pu I SZ u C inv 1( ) Pu u + u I SZ u C 1 u = C u Collecting constants yields te result REFERENCES [1 R E Bank and H Yserentant, On te H 1 -stability of te L 2 -projection onto finite element spaces, Numer Mat, 126 (2014), pp [2 J H Bramble, J E Pasciak, and O Steinbac, On te stability of te L 2 projection in H 1 (Ω), Matematics of Computation, 71 (2001), pp [3 A Ern and J-L Guermond, Teory and Practice of Finite Elements, Springer-Verlag, New York, 2004 [4 J A Fiordilino and M McLauglin, An Artificial Compressibility Ensemble Timestepping Algoritm for Flow Problems, arxiv preprint arxiv: (2017) [5 K J Galvin, New subgrid artificial viscosity Galerkin metods for te Navier-Stokes equations, Comput Metods Appl Mec Engrg, 200 (2011), pp [6 V Girault and P A Raviart, Finite Element Approximation of te Navier-Stokes Equations, Springer, Berlin, 1979 [7 J G Heywood and R Rannacer, Finite-Element Approximation of te Nonstationary Navier-Stokes Problem Part IV: Error Analysis for Second-Order Time Discretization, SIAM J Numer Anal, 27 (1990), pp [8 R Ingram, Approximating fast, viscous fluid flow in complicated domains, PD tesis, University of Pittsburg, Pittsburg, PA, 2011 [9 A Labovscii, A Defect Correction Metod for te Time-Dependent Navier-Stokes Equations, Tecnical report, TR-MATH 07-14, University of Pittsburg, 2007 [10 W Layton, Introduction to te Numerical Analysis of Incompressible, Viscous Flows, SIAM, Piladelpia, 2008 [11 L R Scott and S Zang, Finite element interpolation of nonsmoot functions satisfying boundary conditions, Matematics of Computation, 54 (1990), pp [12 J Sen, On a new pseudocompressibility metod for te incompressible Navier-Stokes equations, Applied Numerical Matematics, 21 (1996), pp [13 R Temam, Navier-Stokes Equations and Nonlinear Functional Analysis, SIAM, Piladelpia, 1995 [14 Y Zang, Y Hou, and J Zao, Error analysis of a fully discrete finite element variational multiscale metod for te natural convection problem, Computers & Matematics wit Applications, 68 (2014), pp