PROJECTION METHOD III: SPATIAL DISCRETIZATION ON THE STAGGERED GRID
|
|
- Barrie Clarke
- 5 years ago
- Views:
Transcription
1 MATHEMATICS OF COMPUTATION Volume 71, Number 237, Pages S (01) Article electronically published on May 14, 2001 PROJECTION METHOD III: SPATIAL DISCRETIZATION ON THE STAGGERED GRID WEINAN E AND JIAN-GUO LIU Abstract. In E & Liu (SIAM J Numer. Anal., 1995), we studied convergence and the structure of the error for several projection methods when the spatial variable was kept continuous (we call this the semi-discrete case). In this paper, we address similar questions for the fully discrete case when the spatial variables are discretized using a staggered grid. We prove that the numerical solution in velocity has full accuracy up to the boundary, despite the fact that there are numerical boundary layers present in the semi-discrete solutions. 1. Introduction In this paper, we continue our study of convergence and error structure for the projection method. In previous work [3], we studied semi-discrete situations when the temporal variable is discretized, but the spatial variable is kept continuous. We proved that as far as velocity is concerned, the method of Chorin and Temam [2, 12] is uniformly first order accurate up to the boundary, and the second order method of [6, 9] is uniformly second order accurate up to the boundary. However, due to the presence of numerical boundary layers, pseudo-pressure has only half order of accuracy near the boundary. We characterized explicitly the structure of the numerical boundary layer, and proved that full accuracy for pressure can be recovered if we subtract the numerical boundary layer profiles from the numerical solution. We also showed, using normal mode analysis, that the boundary layer modes turn into oscillatory modes for the pressure-increment formulation of [13, 1]. For a summary of these results, we refer to [4]. For related convergence results for the semi-discrete case, we refer to [10, 11]. The presence of singular modes in the semi-discrete solutions raises serious doubts on the convergence of the fully discrete method. Indeed, it has been known for some time that accuracy and even convergence can be lost if care is not exercised at the projection step. This is documented carefully in [5]. It has also been known for some time that full accuracy is kept if the spatial discretization is done on a staggered grid. This paper is devoted to a proof of this empirical fact. Our work still leaves open the very important issue of characterizing the minimum condition the spatial discretization has to satisfy in order to guarantee accuracy. The working assumption seems to be that as long as the projection step is truly a projection, i.e. the numerical projection operator P h satisfies: (1) P h is Received by the editor May 19, 1997 and, in revised form, March 1, Mathematics Subject Classification. Primary 65M06, 76M20. Key words and phrases. Viscous incompressible flows, projection method, numerical boundary layer, finite difference, convergence. 27 c 2001 American Mathematical Society
2 28 WEINAN E AND JIAN-GUO LIU self-adjoint; (2) P 2 h = P h (we will refer to this as the projection condition), accuracy in velocity will be kept. It remains to be seen whether this is true in general. We refer to [16] for some progress in this direction. On the other hand, there are many situations in which enforcing the projection condition is difficult. Therefore we are motivated to seek numerical methods which do not require the projection condition. The gauge method [5] has so far proved to be a very attractive alternative in this regard. This paper is organized as follows. In the next section, we review the projection method, with emphasis on spatial discretization on the staggered grid (also known as the MAC grid [7]). In Section 3, we summarize our main results. In Sections 4 and 5 we present the proof for the first and second order methods respectively. 2. Review of the projection methods In primitive variables, the Navier-Stokes equation (NSE) takes the form { t u +(u )u + p = u, (2.1) u =0. Here u =(u, v)isthevelocity,p is the pressure. For simplicity, we will only consider the case when the no-slip boundary condition (BC) is supplemented to (2.1): (2.2) u =0 on Ω, where Ω is an open domain in RI 2 with smooth or piecewise smooth boundary Time discretization. The first order projection method of [2, 12] proceeds in two steps. Step 1: Computing the intermediate velocity field u : u u n +(u n )u n = u, (2.3) u =0, on Ω. Step 2: Projecting to the space of divergence-free vector fields to obtain u n+1 : { u = u n+1 + p n+1, (2.4) u n+1 =0. The projection step enforces the boundary condition: u n+1 n p n+1 (2.5) =0, or =0, on Ω. n Second order schemes. We will concentrate on the formulation in [6, 9], and refer to [4] for a summary of other second order methods. We have u u n +(u n+1/2 )u n+1/2 = u + u n, 2 u + u n = p n 1/2, on Ω, (2.6) u = u n+1 + p n+1/2, u n+1 =0, u n+1 n =0, on Ω.
3 PROJECTION METHOD III 29 In this formulation, the homogeneous Neumann BC for pressure, p n+1/2 =0, on Ω, n is retained. An inhomogeneous BC for u is introduced so that the slip velocity of u n+1 at the boundary is of order () 2. Remark. The nonlinear convection term (u n+1/2 )un+1/2 can be treated in many 3 ways. In Theorems 1 and 2, we use the Adams-Bashforth formula: 2 (un )un 1 2 (un 1 )u n 1. This explicit treatment of the convection term requires a CFL stability constraint. In the convergence analysis, this is realized by an apriori estimate (see (4.58)). This stability restriction and the aprioriestimate are no longer needed in an implicit treatment of the convection term. It is readily seen that the projection step enforces p n+1 (2.7) n = pn p0 = = n n =0, on Ω, for the numerical solution. In general this is not satisfied by the exact solution of (2.1). Therefore we expect that p n / n has O(1) error at the boundary. This causes u and p n to have numerical boundary layers Spatial discretization. We will consider the case when the MAC scheme is used to discretize in x. An illustration of the MAC mesh near the boundary is given in Figure 1. Here pressure is evaluated at the square points (i, j), the u velocity at the triangle points (i ± 1/2,j), and the v velocity at the circle points (i, j ± 1/2). The discrete divergence is computed at the square points: ( u) i,j = u i+1/2,j u i 1/2,j x + v i,j+1/2 v i,j 1/2 y. Γ y Γ x Figure 1. The MAC mesh
4 30 WEINAN E AND JIAN-GUO LIU Other differential operators are discretized as follows: ( u) i+1/2,j = u i+3/2,j 2u i+1/2,j + u i 1/2,j x 2 + u i+1/2,j+1 2u i+1/2,j + u i+1/2,j 1 y 2, ( v) i,j+1/2 = v i+1,j+1/2 2u i,j+1/2 + u i 1,j+1/2 x 2 + v i,j+3/2 2v i,j+1/2 + v i,j 1/2 y 2, (p x ) i+1/2,j = p i+1,j p i,j x (p y ) i,j+1/2 = p i,j+1 p i,j y,, ū i,j+1/2 = 1 4 (u i+1/2,j + u i 1/2,j + u i+1/2,j+1 + u i 1/2,j+1 ), v i+1/2,j = 1 4 (v i+1,j+1/2 + v i+1,j 1/2 + v i,j+1/2 + v i,j 1/2 ), N h (u,a) i+1/2,j = u i+1/2,j a i+3/2,j a i 1/2,j 2 x a i+1/2,j+1 a i+1/2,j 1 + v i+1/2,j, 2 y N h (u,b) i,j+1/2 =ū i,j+1/2 b i+1,j+1/2 b i 1,j+1/2 2 x b i,j+3/2 b i,j 1/2 + v i,j+1/2, 2 y whereweusethenotationn to denote the nonlinear convection term. Clearly the truncation errors of these approximations are of second order. The boundary condition u = 0 is imposed at the vertical physical boundary, whereas v = 0 is imposed at the ghost circle points which are x 2 to the left or right of the physical boundary. The ghost points are eliminated using linear interpolation of the boundary conditions. More explicitly, v 1/2,j +v 1/2,j =0. Similarly the boundary condition v = 0 is imposed at the horizontal physical boundary, but u = 0 is imposed at the ghost triangle points a distance of y 2 away from the physical boundary. One shortcoming of the MAC scheme is the serious constraint on geometry. Although slightly more general situations can be studied, in the present paper we willconcentrateonthesituationwhenω=[ 1, 1] [0, 2π] with periodic boundary condition in the y direction and no-slip boundary condition in the x-direction: u(x, 0,t)=u(x, 2π, t), u( 1,y,t)=0, u(1,y,t) = 0. We will use Ωtodenote the part of the boundary at x = ±1. We will always assume that x y and h =min( x, y). Notation. We will use C to denote generic constants which may depend on the norms of the exact solutions.
5 PROJECTION METHOD III Summary of results The main results of this paper are the following (the constants are independent of and h). Theorem 1. Let (u,p) be a solution of the Navier-Stokes equation (2.1) with smooth initial data u 0 (x) satisfying the compatibility condition (3.1) u 0 (x) =0, y p(x, 0) = 2 xyp(x, 0) = 0, on Ω. Let (u h,p h ) be the numerical solution of the projection method (2.3)-(2.5) coupled with the MAC spatial discretization. Assume that <<h. Then we have (3.2) u u h L + 1/2 p p h L C( + h 2 ), (3.3) where (3.4) p p h p c L C( + h 2 ), p c (x,t) 1/2 β eα 1/2 e α 1 e α x 1 / D+p x h (x 1/2,y,t) eα + 1/2 β e α x+1 /1/2 e α D x 1 +p h (x + 1/2,y,t), α = 1/2 x 2 x arccosh(1 + 2 ), β = x 1/2 (1 e α x/ ) 1. 1/2 Remark 1. As will be seen more clearly in the proof (4.52), at the scaling = Ch 2 the boundary layer in pressure looks like t 1/2 A(y, t)κ j,wherea(y, t) issmooth, κ < 1, κ does not depend on h, andj is the horizontal grid index relative to the lower boundary. Theorem 2. Let (u,p) be a smooth solution of the Navier-Stokes equation (2.1) with smooth initial data u 0 (x) satisfying the compatibility condition (3.5) x α1 α2 y u0 (x) =0, on Ω, for α 1 + α 2 6. Let (u h,p h ) be the numerical solution of the projection method (2.6) coupled with the MAC spatial discretization. Assume that 2 << h. Then we have (3.6) u u h L + 3/2 p p h L + p p h L (0,T,L 2 ) C( 2 + h 2 ), (3.7) where (3.8) p p h p c L C( + h 2 ), p c 1/2 β eα e α 1 e α x 1 /1/2 D x + p h(x 1/2,y,t) eα + 1/2 β e α x+1 /1/2 e α D+ x 1 p h(x + 1/2,y,t), α = 1/2 x 2 x arccosh(1 + ), β = x 1/2 (1 e α x/ ) 1. 1/2
6 32 WEINAN E AND JIAN-GUO LIU Remark 2. The constraint on the size of is only technical and is different from the standard stability condition. The issue of the compatibility conditions on the initial data for the general domain are very complicated and are usually nonlocal [8]. In the case of periodic BC in y and no-slip BC in x (as in the present work), this reduces to the boundary condition (3.5). The correction term in (3.8) can also be given by 1 2 hp, as proposed in Kim & Moin s original paper [9]. The proofs of these results follow the general strategy outlined in Section 3 of [3]. Step 3 is made simpler since we can now use the inverse inequality. 4. First order schemes with spatial discretization We will concentrate on the following version of the first order projection method with the standard MAC spatial discretization: u u n + N h (u n, u n )= h u, u =0, on Ω, (4.1) u = u n+1 + h p n, h u n+1 =0, n u n+1 =0, on Ω. For a = (a, b), c = (c, d), u = (u, v), we define the following discrete inner products on the grid: (4.2) N 1 ((a, c)) = x y + x y N 1 ((u, h p)) = y + x i=1 j=1 N i=1 j=1 N N 1 (( h u,p)) = y + x N a i+1/2,j c i+1/2,j i=1 j=1 N b i,j+1/2 d i,j+1/2, N u i+1/2,j (p i+1,j p i,j ) i=1 j=1 i=1 j=1 N N v i,j+1/2 (p i,j+1 p i,j ), N (u i+1/2,j u i 1/2,j )p i,j i=1 j=1 N (v i,j+1/2 v i,j 1/2 )p i,j, and discrete norms (4.3) u =((u, u)) 1/2, u =max u i,j. i,j Denote h =min( x, y).
7 PROJECTION METHOD III 33 Lemma 4.1. We have the following. (i) Inverse inequality: (4.4) (4.5) (ii) f 1 h f. Poincaré inequality: if f x=±1 =0,then f C h f, where C depends only on the domain. (iii) Suppose n u x=±1 =0; then we have (4.6) (iv) (4.7) (v) (4.8) ((u, h p)) = (( h u,p)). Suppose u x=±1 =0; then we have 2((u, h u)) h u 2 h u 2. Suppose a x=±1 =0and c n x=±1 =0; then we have ((a, N h (u, c))) c h a u W 1,. Proof. The proof of (i iii) is standard. We show (iv). Summation by parts gives (4.9) ((u, h u)) = h u 2 + j [v 0,j+1/2 (v 1,j+1/2 v 0,j+1/2 ) v N,j+1/2 (v N+1,j+1/2 v N,j+1/2 )]. Since v x=±1 =0,wehave (4.10) v 1,j+1/2 = v 0,j+1/2, v N,j+1/2 = v N+1,j+1/2. Hence (4.11) ((u, h u)) = h u 2 +2 j (v 2 0,j+1/2 v 2 N,j+1/2 ). But (4.12) h u 2 h u 2 + j [(v 1,j+1/2 v 0,j+1/2 ) 2 +(v N+1,j+1/2 v N,j+1/2 ) 2 ] = h u 2 +4 j [(v 0,j+1/2 ) 2 +(v N+1,j+1/2 ) 2 ]. Combination of (4.11) and (4.12) gives (4.7). To show (v), denote I =((a, N h (u, c))). We have I = x y i,j a i+1/2,j (u i+1/2,j D x 0 c i+1/2,j + v i+1/2,j D y 0 c i+1/2,j) (4.13) + x y i,j b i,j+1/2 (ū i,j+1/2 D x 0 d i,j+1/2 + v i,j+1/2 D y 0 d i+1/2,j),
8 34 WEINAN E AND JIAN-GUO LIU whereweusedthenotationsa =(a, b) andc =(c, d). Summation by parts gives (4.14) I = x y c i+1/2,j [D0 x (u i+1/2,j a i+1/2,j )+D y 0 ( v i+1/2,ja i+1/2,j )] i,j x y i,j d i,j+1/2 [D x 0(ū i,j+1/2 b i,j+1/2 )+D y 0 (v i,j+1/2b i,j+1/2 )] x y j (ū N+1,j+1/2 d N,j+1/2 ū N,j+1/2 d N+1,j+1/2 )D x b N+1,j+1/2 1 4 x y j (ū 1,j+1/2 d 0,j+1/2 ū 0,j+1/2 d 1,j+1/2 )D x + b 0,j+1/2. Here we have used the fact that (4.15) b 1,j+1/2 = b 0,j+1/2, b N,j+1/2 = b N+1,j+1/2. Now, (4.8) follows directly. This completes the proof of the lemma. We set ε = 1/2, ξ =(x +1)/ε, x i = 1 +i x, ξ i = i ξ, ξ = x/ε, t n = n, t n 1/2 =(n 1/2), n =1, 2,. Clearly, we have (4.16) D ξ + a(ξ i,y j,t) = a(ξ i+1,y j,t) a(ξ i,y j,t) ξ = ε a(x i+1/ε, y j,t) a(x i /ε, y j,t) x = εd x + a(x i/ε, y j,t). This shows that D ξ + = εdx +. We will use the notation (4.17) and Dξ 2 = Dξ Dξ +, D2 y = Dy Dy +, (4.18) ξ =(D ξ +, 0), y =(0,D y + ). Denote the solutions of (4.1) as (u h, u h,p h). Motivated by the discussions in 4 of [3], we make the following ansatz, valid at t n = n, n =1, 2, : u h (x,t)=u 0(x,t)+ ε j [u j (x,t)+a j (ξ, y, t)], j=2 (4.19) u h (x,t)=u 0 (x,t)+ j=2 ε j u j (x,t), p h (x,t)=p 0 (x,t)+ j=1 ε j [p j (x,t)+ϕ j (ξ, y, t)]. Note that the functions involved are defined only on the numerical grid. So these formulas and the following ones should be understood as being valid on the grid. We have (4.20) h u h = hu 0 + ε j ( h u j + ε 2 Dξ 2 a j + D y 2 a j ), j=2
9 PROJECTION METHOD III 35 (4.21) h u h = h u 0 + j=2 ε j h u j, (4.22) h p h = h p 0 + ε 1 ξ ϕ 0 + y ϕ 0 + j=1 ε j ( h p j + ε 1 ξ ϕ j + y ϕ j ), u n+1 h (x) = u 0 (x,t n+1 )+ j=2 ε j u j (x,t n+1 ) (4.23) = k=0 1 k! ε2k u (k) 0 (x,tn )+ ε j 1 k! ε2k u (k) j (x,t n ). j=2 k=0 Next we substitute these relations into (4.1) in order to determine the coefficients of ε j in (4.19). We get hierarchies of equations by collecting equal powers of ε. The first equation in (4.1) gives (4.24) u 2 + a 2 u 2 + N h (u 0, u 0 )= hu Dξ 2 a 2. For j 1, j (4.25) u j+2 + a j+2 u j+2 + N h (u k, u j k )= h u j + D ξ 2 a j+2 + D y 2 a j. k=0 The second equation in (4.1) implies that (4.26) u 2 + a 2 = u 2 + t u 0 + h p 0 + ξ ϕ 1 + y ϕ 0. For j =2l 1, l 1, l (4.27) u j+2 + a 1 j+2 = u j+2 + t u j + h p j + ξ ϕ j+1 + y ϕ j + k! u(k) j 2k+2. k=2 For j =2l, l 1, u j+2 + a j+2 = u j+2 + t u j + h p j + ξ ϕ j+1 + y ϕ j (4.28) 1 l 1 + (l +1)! u(l+1) 0 + k! u(k) j 2k+2. From the third equation in (4.1), we obtain (4.29) h u j =0, j =0, 1,. The boundary conditions become (4.30) u j + a j =0, Dx + p j 1 + D ξ + ϕ j =0, at x = 1, ξ =0, for j>0. Next we go through all these equations, one by one, to see if they are solvable. Since this is very similar to 4.1 in [3], we will only give a summary of results. The coefficients in the expansions (4.19) can be obtained successively in the following order: (4.31) k=2 t u 0 + h p 0 + N h (u 0, u 0 )= h u 0, h u 0 =0, u 0 =0, at x = ±1, u 0 (, 0) = u 0 ( ).
10 36 WEINAN E AND JIAN-GUO LIU Using the following lemma, we know that (4.31) has a smooth solution in the sense that the divided differences of various orders are bounded. Lemma 4.2. Let (u,p) be a solution of the Navier-Stokes equation (2.1) with smooth initial data u 0 (x) satisfying some compatibility conditions. Let (u 0,p 0 ) be a solution of (4.31). Then(u 0,p 0 ) is smooth in the sense that its discrete derivatives are bounded. Moreover, we have (4.32) u u 0 L + p p 0 L Ch 2. The proof of this lemma, as well as Lemma 4.3, can be found in [15]. We next have (4.33) u 2 = u 2 + t u 0 + h p 0, (4.34) This gives (4.35) where (4.36) (4.37) (4.38) { ϕ1 = D 2 ξ ϕ 1, D ξ + ϕ 1 ξ=0 = D x + p 0 x= 1. ϕ 1 (ξ, y, t) =βd x + p 0( 1,y,t) e αξ, α = 1 ξ arccosh (1 + ξ 2 /2), β = ξ(1 e α ξ ) 1, a 2 = D ξ + ϕ 1, b 2 =0, u 3 = u 3, (4.39) (4.40) ϕ 2 =0, a 3 =0, b 3 = Dy + ϕ 1, t u 2 + h p 2 + N h (u 0, u 2 )+N h (u 2, u 0 ) = h u 2 + h ( t u 0 + h p 0 ) 1 2 t 2 u 0, h u 2 =0, u 2 x= 1 = h p 0 x= 1 ξ ϕ 1 ξ=0, on Ω. With suitable initial data, we know from the following lemma that (4.40) has a smooth solution. Lemma 4.3. Let (u,p) be a solution of the linear ODE t u + h p + N h (u 0, u)+n h (u, u 0 )= h u + f, h u =0, (4.41) u = g, at x = ±1, u(, 0) = u 0 ( ), where f, g and u 0 smooth and satisfy some compatibility conditions. Then (u,p) is smooth in the sense that its divided differences of various order are bounded. Continuing in this fashion, we get { ϕ3 = Dξ 2ϕ 3 + Dy 2 ϕ 1, (4.42) D+ϕ ξ 3 ξ=0 = D+p x 2 x= 1. The solution for (4.42) is (4.43) ϕ 3 (ξ, y, t) =βd+p x 2 e αξ + β 1 (ξ + γ)d+d x y 2 p 0 x= 1 e αξ,
11 PROJECTION METHOD III 37 where (4.44) (4.45) (4.46) (4.47) β 1 = 1 (1 e α ξ )(e α ξ e α ξ ), γ = ξ e α ξ 1 e α ξ, b 4 =0, a 4 = D ξ + ϕ 3. t u 3 + h p 3 + N h (u 0, u 3 )+N h (u 3, u 0 )= h u 3 + h h p 2, h u 3 =0, u 3 x= 1 = y ϕ 1 ξ=0, { ϕ4 = D 2 ξ ϕ 4, D ξ + ϕ 4 ξ=0 = D x + p 3 x= 1, (4.48) a 5 = D+ϕ ξ 4, b 5 = D+ϕ y 3. Obviously this procedure can be continued, and we obtain { ϕj = Dξ 2ϕ j + Dy 2ϕ j 2, (4.49) D ξ + ϕ j ξ=0 = D+ x p j 1 x= 1, (4.50) (4.51) Now if we let (4.52) ϕ j = [j/2] k=0 F j,k (y)ξ k e αξ, a j = Dξ + ϕ j 1, b j = Dy + ϕ j 2, 2N U = u 0 + ε j (u j + a j ), U n = u 0 + P n = p 0 + j=1 2N j=1 2N j=1 ε j u j, ε j (p j + ϕ j )+ε 2N+1 ϕ 2N+1, then we have (4.53) U U n + N h (U n,u n )= h U + α f, U =0, at x = ±1, U = U n+1 + h P n + α+1 g, h U n+1 =0, D x + P n = n U n+1 =0, at x = ±1, where α = N 1/2; f and g are bounded and smooth if (u 0,p 0 ) is sufficiently smooth. It is easy to see that (4.54) max U n ( ) W 1, C. 0 t T
12 38 WEINAN E AND JIAN-GUO LIU We point out that the constructed solution in (4.52) is different from the discretized solution. The constructed solution consists of two parts: the leading term is the exact solution of the NSE on MAC grid (4.31) and hence is bounded (see [15]), and the remaining terms are the solutions of the linearized (linearized on the exact solution) NSE (4.40), (4.46), and hence are bounded. As a direct consequence, the truncation terms f and g in (4.53) are also bounded. For the initial data, we have (4.55) U 0 (x) =u 0 (x)+w 0 (x), where w 0 is a bounded function. Furthermore, under the compatibility condition (3.5), we can construct a better approximate initial data (4.57) U 0 (x) =u 0 (x)+ 2 w 0 (x). The construction of the initial w 0 in (4.55) and (4.57) in general is very complicated. This is exactly where we need the compatibility conditions. In the case of periodic BC for y and no-slip boundary condition in x (the domain we considered in this paper), the construction of w 0 is more straightforward. The details of the construction are very similar to those used in the analysis of the gauge method in [14]. In that paper the details of the construction of w 0 were given. Proof of Theorem 1. Assume apriorithat (4.58) max 0 t n T un W 1, C. We will justify this aprioriassumption later (see (4.67)). In the following estimates, the constant will sometimes depend on C and C, and we will show that it is uniformly bounded. Later on we will estimate C. Let (4.59) e n = U n u n, e = U u, q n = P n p n. Subtracting (4.53) from (4.1), we get the following error equation: e e n + N h (e n,u n )+N h (u n, e n )= h e + α f n, e =0, at x = ±1, (4.60) e n+1 e + h q n = α g n, h e n+1 =0, D + x qn = e n+1 n =0, at x = ±1, e 0 = α w 0. Taking the scalar product of the first equation of (4.60) with 2e and integrating by parts, we obtain e 2 e n 2 + e e n 2 + h e 2 (4.61) 2α+1 f n 2 + e 2 2 ((e, N h (e n,u n ))) 2 ((e, N h (u n, e n ))) 2α+1 f n 2 + C( e 2 + e n 2 )+ 1 2 he 2.
13 PROJECTION METHOD III 39 Here we used Lemma 4.1. Taking the scalar product of the second equation of (4.60) with 2e n+1 yields (4.62) e n+1 2 e 2 + e n+1 e 2 e n α+1 g n 2. Combining (4.61) and (4.62), we get e n+1 2 e n 2 + e e n 2 + e n+1 e 2 + h e 2 (4.63) C ( e n e n 2 )+ 2α+1 ( f n 2 + g n 2 ). Applying the discrete Gronwall lemma to the last inequality, we arrive at (4.64) e n + e e n + e n+1 e + 1/2 h e C 1 α. Using the second equation of (4.60), we have (4.65) e n + h q n C 1 α. Now by the inverse inequality (4.4) we have e n L + h q n L + h e n α (4.66) W 1, C 1 h. If N =3and α << h 2,andwechoose small enough, we will always have (4.67) e n+1 W 1, 1. Therefore in (4.58) we can choose (4.68) C =1+ max U n ( ) W 1,, n [ T ]+1 which depends only on the exact solution (u,p). This proves that (4.69) u 0 u h L + p 0 p h L 2 + 1/2 p 0 p h L + p 0 p h p c L C. But we also have, from Lemma 4.2, (4.70) u u 0 L + p p 0 L Ch 2. Thus (4.71) u u h L + p p h L 2 + 1/2 p p h L + p p h p c L C( + h 2 ). This completes the proof of Theorem Second order schemes with spatial discretization In this section we carry out the same program as in 4 for the second order method (2.6) with the standard MAC spatial discretization: u u n u + u n = h, 2 u + u n = h p n 1/2, at x = ±1, (5.1) u = u n+1 + h p n+1/2, h u n+1 =0, n u n+1 =0, at x = ±1. Here we leave out the nonlinear term, since it does not affect the major steps but substantially complicates the presentation.
14 40 WEINAN E AND JIAN-GUO LIU We begin with the following ansatz: u (x) =u 0 (x,t n )+ ε j [u j (x,tn )+a j (ξ, y, tn )], j=2 (5.2) u n (x) =u 0 (x,t n )+ j=4 ε j u j (x,t n ), p n 1/2 (x) =p 0 (x,t n 1/2 )+εϕ 1 (ξ, y, t n 1/2 )+ε 3 ϕ 3 (ξ, y, t n 1/2 ) + j=4 ε j [p j (x,t n 1/2 )+ϕ j (ξ, y, t n 1/2 )]. Here again we set ε = 1/2, ξ =(x +1)/ε, t n = n, t n 1/2 =(n 1/2), n =1, 2,. The formulas are to be understood as being valid at the grid points. Substituting (5.2) into (5.1) and collecting equal powers of ε, we get the following equations: From the first equation in (5.1), we get (5.3) u 2 + a 2 u 2 = 1 2 ( hu 0 + D 2 ξ a 2 + h u 0 ). For j 1, (5.4) u j+2 + a j+2 u j+2 = 1 2 ( hu j + D 2 ξ a j+2 + D 2 y a j + hu j ). From the third equation in (5.1), we get (5.5) u 2 + a 2 = u 2 + t u 0 + h p 0 + ξ ϕ 1, (5.6) For j =2l, u 3 + a 3 = u 3 + t u 1 + h p 1 + ξ ϕ 2 + y ϕ 1. (5.7) u j+2 + a j+2 = u j+2 + t u j + h p j + ξ ϕ j+1 + y ϕ j 1 l 1 + (l +1)! u(l+1) 0 + k! u(k) j 2k l l! hp (l) + l k=1 k=2 k=2 1 2 k k! ( ξϕ (k) j 2k+1 + yϕ (k) j 2k ). l k=1 1 2 k k! hp (k) j 2k For j =2l +1, u j+2 + a j+2 = u j+2 + t u j + h p j + ξ ϕ j+1 + y ϕ j (5.8) l+1 1 l + k! u(k) j 2k k k! ( hp (k) j 2k + ξϕ (k) j 2k+1 + yϕ (k) j 2k ). k=1 From the incompressibility condition, we get (5.9) h u j =0, for j 0.
15 PROJECTION METHOD III 41 The boundary conditions imply that for x = 1, ξ =0, (5.10) (5.11) (5.12) for j =2l, l 1, u 0 =0, u 2 + u 2 + a 2 = hp 0 + ξ ϕ 1, u 3 + u 3 + a 3 = h p 1 + ξ ϕ 2 + y ϕ 1 ; (5.13) u j + u j + a j = h p j 2 + ξ ϕ j 1 + y ϕ j 2 + ( 1)l l 1 (l 1)! hp (l 1) 0 l 2 + k=1 for j =2l +1,l 1, ( 1) k l 1 2 k k! hp (k) j 2k 2 + k=1 ( 1) k 2 k k! ( ξϕ (k) j 2k 1 + yϕ (k) j 2k 2 ); u j + u j + a j = h p j 2 + ξ ϕ j 1 + y ϕ j 2 (5.14) l 1 + k=1 ( 1) k 2 k k! ( hp (k) j 2k 2 + ξϕ (k) j 2k 1 + yϕ (k) j 2k 2 ); and for j 0 (5.15) D x +p j + D ξ +ϕ j+1 =0. Next we go through all these equations, one by one, to see if they are solvable. It can be checked that the coefficients in the expansions (5.2) can be obtained successively in the following order: (5.16) t u 0 + h p 0 = h u 0, h u 0 =0, u 0 =0, at x = ±1, (5.17) (5.18) (5.19) where (5.20) (5.21) u 2 = u 2 + t u 0 + h p 0, { ϕ1 = 1 2 D 2 ξ ϕ 1, D ξ +ϕ 1 ξ=0 = D x +p 0 x= 1, ϕ 1 = βd x +p 0 x= 1 e αξ, α = 1 ξ arccosh (1 + ξ 2 ), β = ξ(1 e α ξ ) 1, a 2 = Dξ + ϕ 1, b 2 =0.
16 42 WEINAN E AND JIAN-GUO LIU We next have (5.22) (5.23) u 3 = u 3, ϕ 2 =0, a 3 =0, b 3 = D y +ϕ 1, (5.24) The solution for (5.24) is (5.25) where (5.26) β 1 = { ϕ3 = 1 2 (D ξ 2ϕ 3 + Dy 2ϕ 1), D ξ + ϕ 3 ξ=0 =0. ϕ 3 (y, ξ, t) =β 1 (ξ + γ)d x +D 2 y p 0 x= 1 e αξ, 1 2(1 e α ξ )(e α ξ e α ξ ), γ = ξ e α ξ 1 e α ξ, (5.27) a 4 = 1 2 Dξ + tϕ 1 + D ξ + ϕ 3, b 4 =0, (5.28) u 4 = u t u t h p 0, (5.29) (5.30) (5.31) ϕ 4 =0, a 5 =0, b 5 = 1 2 Dy + t ϕ 1 + D y +ϕ 3, u 5 = u 5, t u 4 + h p 4 = h u h( 2 t u 0 + t h p 0 ) t u t hp 0, h u 4 =0, u 4 x= 1 = 1 2 ( t h p t ξ ϕ 1 ) x= 1,ξ=0. Now if we let (5.32) 2N U = u 0 + ε j (u j + a j ), U n = u 0 + j=1 2N j=1 ε j u j, then we have (5.33) P n 1/2 = p 0 + 2N j=1 ε j (p j + ϕ j )+ε 2N+1 ϕ 2N+1, U U n U + U n = h + α f, 2 U + U n = h P n 1/2, at x = ±1, U = U n+1 + h P n+1/2 + α+1 g, h U n+1 =0, D x +P n+1/2 = n U n+1 =0, at x = ±1,
17 PROJECTION METHOD III 43 where α = N 1/2; f and g are bounded and smooth if (u 0,p 0 ) is sufficiently smooth. It is easy to see that (5.34) max U n ( ) W 1, C. 0 t T For the initial approximation, we have (5.35) U 0 (x) =u 0 (x)+ 2 w 0 (x) without the extra compatibility condition, and (5.36) U 0 (x) =u 0 (x)+ 4 w 0 (x) with the compatibility condition (3.17). Proof of Theorem 2. Assume apriorithat (5.37) max 0 t n T un W 1, C. As in the proof of Theorem 1, we let (5.38) where (5.39) e n = U n u n, e = Û û, q n = P n 1/2 p n 1/2, 2û = u + u n h p n 1/2, 2Û = U + U n h P n 1/2. From (5.1) and (5.34), we get (5.40) 2(e e n ( ) + h q n 1 ) 2 hq n = h e N h(e n 1,U n 1 ) N h(u n 1, e n 1 ) 3 2 N h(e n,u n ) 3 2 N h(u n, e n )+ α f n, e =0, at x = ±1, e n+1 + e n 2e + h (q n+1 q n )= α g n, h e n+1 =0, D +q x n+1/2 = e n+1 n =0, at x = ±1, e 0 = α w 0. Taking the scalar product of the first equation of (5.40) with e and integrating by parts, we get (5.41) e 2 e n 2 + e e n e e 2 e (q n 1 Ω 2 qn ) dx + C 2α+1 f n 2 +C ( e n 2 + e n e 2 )+ 1 2 e 2.
18 44 WEINAN E AND JIAN-GUO LIU Taking the scalar product of the second equation of (5.42) with e n+1,weobtain (5.42) e n+1 2 e 2 + e n+1 e ( en+1 2 e n 2 ) 1 2 en+1 e n 2 C 2α+1 g n 2 + C e n+1 2. Combining these two estimates, we get (5.43) e n+1 2 e n 2 + e n+1 + e n 2e 2 + h e 2 + h e 2 2 ((e, h (q n 1 2 hq n ))) + C ( e n 2 + e n e n+1 2 ) +C 2α+1 ( f n 2 + g n 2 ). To estimate the first term on the right hand side of (5.43), we let (5.44) I 2((e, h (q n 1 2 hq n ))) = 2((e, h q n )) 2 (( h e, h q n )) I 1 + I 2. Using the second equation and integrating by parts, we can write the first term as I 1 = 2((e, h q n )) = 2 (( h (q n+1 q n ), h q n )) α+2 ((g n, h q n )) (5.45) = ( h q n+1 2 h q n 2 ) h (q n+1 q n ) 2 α+2 ((g n, h q n )). Since (5.46) h (q n+1 q n ) 2 = 1 2 en+1 + e n 2e α+1 g n 2 + α+1 ((g n, e n+1 + e n 2e )), we have (5.47) I 1 = ( h q n+1 2 h q n 2 )+ 1 2 en+1 + e n 2e α+1 g n 2 + α+1 ((g n, e n+1 + e n 2e )) α+2 ((g n, h q n )).
19 PROJECTION METHOD III 45 Next we rewrite the second term as (5.48) I 2 = 2 (( h e, h q n )) = (( h (q n+1 q n ), h q n )) 1 2 α+3 (( h g n, h q n )) = ( h q n+1 2 h q n 2 ) h (q n+1 q n ) α+3 (( h g n, h q n )) = ( h q n+1 2 h q n 2 )+ h e α+3 h g n 2 α+2 (( h g n, h e )) 1 2 α+3 (( h g n, h q n )). Combining these two terms, we arrive at (5.49) I = ( h q n+1 2 h q n 2 ) ( h q n+1 2 h q n 2 ) en+1 + e n 2e 2 + h e 2 + α+1 ((g n, e n+1 + e n 2e )) α+2 ((g n, h q n )) α+2 (( h g n, h e )) 1 2 α+3 (( h g n, h q n )) α+3 h g n α+1 g n 2. This gives I ( h q n+1 2 h q n 2 ) ( h q n+1 2 h q n 2 ) (5.50) en+1 + e n 2e 2 + h e 2 + e n+1 + e n 2e h q n h q n α+1 ( g n 2 + g n H 2 ). 1 Going back to (5.45), we obtain (5.51) e n+1 2 e n en+1 + e n 2e 2 + h e ( h q n+1 2 h q n 2 ) ( h q n+1 2 h q n 2 ) 3 h q n h q n 2 + C ( e n 2 + e n e n+1 2 ) +C 2α+1 ( f n 2 + g n H 2 ). 1 Gronwall s lemma gives (5.52) e n + e + h q n + 3/2 h q n + 1/2 h e C 1 α. Now by the inverse inequality (4.4) we have e n L + h e n W 1, + h q n α (5.53) L C 1 h.
20 46 WEINAN E AND JIAN-GUO LIU If N =5and α << h 2,andwechoose small enough, we will always have (5.54) e n+1 L 1. Therefore in (5.39) we can choose (5.55) C =1+ max n [ ] U n ( ) W 1,, T which depends only on the exact solution (u,p). This proves that (5.56) u 0 u h L + p 0 p h L 2 + 1/2 p 0 p h L + p 0 p h p c L C 2. From Lemma 4.2, we have (5.57) u u h L + p p h L 2 + 1/2 p p h L + p p h p c L C( 2 + h 2 ). This completes the proof of Theorem Conclusions The main conclusion from the analysis presented in this paper is that the numerical error in the projection method has the same structure as in the semi-discrete case (analyzed in [4]) after it is spatially discretized on a staggered grid. This paper leaves open the very interesting question of what happens on a regular nonstaggered grid. Some progress in this direction has been made in [16]. Acknowledgments The research of the first author was supported in part by ONR grant N and NSF grant DMS The research of the second author was supported in part by ONR grant N and NSF grants DMS and DMS References 1. J.B. Bell, P. Colella, and H.M. Glaz, A second-order projection method for the incompressible Navier Stokes equations, J. Comput. Phys. 85 (1989), MR 90i: A.J. Chorin, Numerical solution of the Navier-Stokes equations, Math. Comp. 22 (1968), MR 39: Weinan E and J.-G. Liu, Projection method I: convergence and numerical boundary layers, SIAM J. Numer. Anal. 32 (1995), ; Projection method II: Godunov-Ryabenki analysis, SIAM J. Numer. Anal. 33 (1996), MR 96e:65061; MR 97i: Weinan E and J.-G. Liu, Projection method for viscous incompressible flows, in Numerical Methods in Applied Sciences, edited by W. Cai and C.-W. Shu, Science Press New York, Ltd, New York, (1996), Weinan E and J.-G. Liu, Gauge method for viscous incompressible flows, preprint, (1996) 6. M. Fortin, R. Peyret, and R. Temam, Résolution numérique des équations de Navier-Stokes pour un fluide incompressible, J. Mecanique, 10 (1971) MR 54: F.H. Harlow and J.E. Welch, Numerical calculation of time-dependent viscous incompressible flowoffluidwithfreesurface,phys. Fluids, 8 (1965) J. Heywood and R. Rannacher, Finite element approximation of the nonstationary Navier- Stokes problem, part III. Smoothing property and higher order error estimates for spatial discretization, SIAM J. Numer. Anal., 25, (1988) MR 89k: J. Kim and P. Moin, Application of a fractional-step method to incompressible Navier-Stokes equations, J. Comp. Phys. 59 (1985), MR 87a: J. Shen, On error estimates of projection methods for Navier-Stokes equations: first order schemes, SIAM J. Numer. Anal. 29 (1992), MR 92m:35213
21 PROJECTION METHOD III J. Shen, On error estimates of some higher order projection and penalty-projection methods for Navier-Stokes equations, Numer. Math. 62 (1992), MR 93a: R. Temam, Sur l approximation de la solution des equations de Navier-Stokes par la méthode des fractionnarires II, Arch. Rational Mech. Anal. 33 (1969), MR 39: J. van Kan, A second-order accurate pressure-correction scheme for viscous incompressible flow, SIAM J. Sci. Stat. Comp. 7 (1986), MR 87h: C. Wang and J.-G. Liu, Convergence of gauge method for incompressible flow, Math.Comp., 69 (2000), MR 20001a: B.R. Wetton, Ph.D Dissertation, New York University, B.R. Wetton, Error analysis for Chorin s original fully discrete projection method and regularizations in space and time, SIAM J. Numer. Anal. 34 (1997), MR 99c:65161 Courant Institute of Mathematical Sciences, New York, New York address: weinan@cims.nyu.edu Institute for Physical Science and Technology and Department of Mathematics, University of Maryland, College Park, Maryland address: jliu@math.umd.edu
Projection Method with Spatial Discretizations
Projection Method with Spatial Discretizations Weinan E School of Mathematics Institute for Advanced Study Princeton, NJ 08540 and Jian-Guo Liu Courant Institute of Mathematical Sciences New York, NY 10012
More information(1:1) 1. The gauge formulation of the Navier-Stokes equation We start with the incompressible Navier-Stokes equation 8 >< >: u t +(u r)u + rp = 1 Re 4
Gauge Finite Element Method for Incompressible Flows Weinan E 1 Courant Institute of Mathematical Sciences New York, NY 10012 Jian-Guo Liu 2 Temple University Philadelphia, PA 19122 Abstract: We present
More informationGauge finite element method for incompressible flows
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 34: 701 710 Gauge finite element method for incompressible flows Weinan E a, *,1 and Jian-Guo Liu b,2 a Courant Institute
More informationSECOND-ORDER FULLY DISCRETIZED PROJECTION METHOD FOR INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
Tenth MSU Conference on Differential Equations and Computational Simulations. Electronic Journal of Differential Equations, Conference 3 (06), pp. 9 0. ISSN: 07-669. URL: http://ejde.math.txstate.edu or
More informationCONVERGENCE OF GAUGE METHOD FOR INCOMPRESSIBLE FLOW CHENG WANG AND JIAN-GUO LIU
MATHEMATICS OF COMPUTATION Volume 69, Number 232, Pages 135{1407 S 0025-571(00)0124-5 Article electronically published on March 24, 2000 CONVERGENCE OF GAUGE METHOD FOR INCOMPRESSIBLE FLOW CHENG WANG AND
More informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationFourth Order Convergence of Compact Finite Difference Solver for 2D Incompressible Flow
Fourth Order Convergence of Compact Finite Difference Solver for D Incompressible Flow Cheng Wang 1 Institute for Scientific Computing and Applied Mathematics and Department of Mathematics, Indiana University,
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 informationFast Solvers for Unsteady Incompressible Flow
ICFD25 Workshop p. 1/39 Fast Solvers for Unsteady Incompressible Flow David Silvester University of Manchester http://www.maths.manchester.ac.uk/~djs/ ICFD25 Workshop p. 2/39 Outline Semi-Discrete Navier-Stokes
More information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
More informationFinite Difference Methods for 3D Viscous Incompressible Flows in the Vorticity Vector Potential Formulation on Nonstaggered Grids
JOURNAL OF COMPUTATIONAL PHYSICS 138, 57 82 (1997) ARTICLE NO. CP975815 Finite Difference Methods for 3D Viscous Incompressible Flows in the Vorticity Vector Potential Formulation on Nonstaggered Grids
More informationERROR ESTIMATES FOR SEMI-DISCRETE GAUGE METHODS FOR THE NAVIER-STOKES EQUATIONS : FIRST-ORDER SCHEMES
MATHEMATICS OF COMPUTATION Volume, Number, Pages S 5-578XX- ERROR ESTIMATES FOR SEMI-DISCRETE GAUGE METHODS FOR THE NAVIER-STOKES EQUATIONS : FIRST-ORDER SCHEMES RICARDO H. NOCHETTO AND JAE-HONG PYO Abstract.
More informationA Review of Solution Techniques for Unsteady Incompressible Flow
Zeist 2009 p. 1/57 A Review of Solution Techniques for Unsteady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/57 PDEs Review : 1966 1999 Update : 2000
More informationProjection Method II : Godunov-Ryabenki Analysis
Projection Method II : Godunov-Ryabenki Analysis Weinan E 1 and Jian-Guo Liu 2 Abstract: This is the second of a series of papers on the subject of projection methods for viscous incompressible flow calculations.
More informationON THE ERROR ESTIMATES FOR THE ROTATIONAL PRESSURE-CORRECTION PROJECTION METHODS
Submitted to Math. Comp. ON THE ERROR ESTIMATES FOR THE ROTATIONAL PRESSURE-CORRECTION PROJECTION METHODS J.L. GUERMOND 1 AND JIE SHEN 2 Abstract. In this paper we study the rotational form of the pressure-correction
More informationA Second Order Energy Stable Scheme for the Cahn-Hilliard-Hele-Shaw Equations
A Second Order Energy Stable Scheme for the Cahn-Hilliard-Hele-Shaw Equations Wenbin Chen Wenqiang Feng Yuan Liu Cheng Wang Steven M. Wise August 9, 07 Abstract We present a second-order-in-time finite
More informationNONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS
NONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS J.-L. GUERMOND 1, Abstract. This paper analyzes a nonstandard form of the Stokes problem where the mass conservation
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationNumerische Mathematik
Numer. Math. 2007 07:669 705 DOI 0.007/s002-007-004-z Numerische Mathematik A fourth-order numerical method for the planetary geostrophic equations with inviscid geostrophic balance Roger Samelson Roger
More informationASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE
MATEMATICS OF COMPUTATION Volume 73, Number 246, Pages 517 523 S 0025-5718(0301570-9 Article electronically published on June 17, 2003 ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR TE POINTWISE GRADIENT
More informationRunge-Kutta-Chebyshev Projection Method
Runge-Kutta-Chebyshev Projection Method Zheming Zheng Linda Petzold Department of Mechanical Engineering, University of California Santa Barbara, Santa Barbara, CA 93106, USA June 8, 2006 This work was
More informationFinite Difference Schemes for Incompressible Flows in the Velocity Impulse Density Formulation
JOURNAL OF COMPUTATIONAL PHYSICS 130, 67 76 (1997) ARTICLE NO. CP965537 Finite Difference Schemes for Incompressible Flows in the Velocity Impulse Density Formulation Weinan E 1 Courant Institute of Mathematical
More informationON THE ERROR ESTIMATES FOR THE ROTATIONAL PRESSURE-CORRECTION PROJECTION METHODS
MATHEMATICS OF COMPUTATION Volume 73, Number 248, Pages 1719 1737 S 0025-5718(03)01621-1 Article electronically published on December 19, 2003 ON THE ERROR ESTIMATES FOR THE ROTATIONAL PRESSURE-CORRECTION
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationPointwise convergence rate for nonlinear conservation. Eitan Tadmor and Tao Tang
Pointwise convergence rate for nonlinear conservation laws Eitan Tadmor and Tao Tang Abstract. We introduce a new method to obtain pointwise error estimates for vanishing viscosity and nite dierence approximations
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationGlobal well-posedness of the primitive equations of oceanic and atmospheric dynamics
Global well-posedness of the primitive equations of oceanic and atmospheric dynamics Jinkai Li Department of Mathematics The Chinese University of Hong Kong Dynamics of Small Scales in Fluids ICERM, Feb
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationStability and Convergence Analysis of Fully Discrete Fourier Collocation Spectral Method for 3-D Viscous Burgers Equation
J Sci Comput (01) 53:10 18 DOI 10.1007/s10915-01-961-8 Stability and Convergence Analysis of Fully Discrete Fourier Collocation Spectral Method for 3-D Viscous Burgers Equation Sigal Gottlieb Cheng Wang
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
More informationON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS
ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS Abdelhafid Younsi To cite this version: Abdelhafid Younsi. ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS. 4 pages. 212. HAL Id:
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 informationGlobal Weak Solution of Planetary Geostrophic Equations with Inviscid Geostrophic Balance
Global Weak Solution o Planetary Geostrophic Equations with Inviscid Geostrophic Balance Jian-Guo Liu 1, Roger Samelson 2, Cheng Wang 3 Communicated by R. Temam) Abstract. A reormulation o the planetary
More informationIn honor of Prof. Hokee Minn
In honor of Prof. Hokee Minn It is my greatest honor to be here today to deliver the nd Minn Hokee memorial lecture. I would like to thank the Prof. Ha for the nomination and the selecting committee for
More informationAccurate Projection Methods for the Incompressible Navier-Stokes Equations
Accurate Projection Methods for the Incompressible Navier-Stokes Equations David L. Brown 1 Center for Applied Scientific Computing Lawrence Livermore National Laboratory, Livermore, CA 94551, USA E-mail:
More informationA STABILIZED NONCONFORMING QUADRILATERAL FINITE ELEMENT METHOD FOR THE GENERALIZED STOKES EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 9, Number 2, Pages 449 457 c 2012 Institute for Scientific Computing and Information A STABILIZED NONCONFORMING QUADRILATERAL FINITE ELEMENT
More informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationA Simple Compact Fourth-Order Poisson Solver on Polar Geometry
Journal of Computational Physics 182, 337 345 (2002) doi:10.1006/jcph.2002.7172 A Simple Compact Fourth-Order Poisson Solver on Polar Geometry Ming-Chih Lai Department of Applied Mathematics, National
More informationhal , version 6-26 Dec 2012
ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS ABDEHAFID YOUNSI Abstract. In this paper, we give a new regularity criterion on the uniqueness results of weak solutions for the 3D Navier-Stokes equations
More informationON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS. Chun Liu and Jie Shen
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 7, Number2, April2001 pp. 307 318 ON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS Chun Liu and Jie Shen Department
More informationA posteriori regularity of the three-dimensional Navier-Stokes equations from numerical computations
A posteriori regularity of the three-dimensional Navier-Stokes equations from numerical computations Sergei I. Chernyshenko, Aeronautics and Astronautics, School of Engineering Sciences, University of
More informationSPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM. Mirjana Stojanović University of Novi Sad, Yugoslavia
GLASNIK MATEMATIČKI Vol. 37(57(2002, 393 403 SPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM Mirjana Stojanović University of Novi Sad, Yugoslavia Abstract. We introduce piecewise interpolating
More information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More information1. Introduction. Consider nonstationary Navier Stokes equations with a nonslip boundary condition
SIAM J. SCI. COMPUT. Vol. 3, No., pp. 398 49 c 8 Society for Industrial and Applied Mathematics LONG TIME NUMERICAL SOLUTION OF THE NAVIER STOKES EQUATIONS BASED ON A SEQUENTIAL REGULARIZATION FORMULATION
More informationLow Froude Number Limit of the Rotating Shallow Water and Euler Equations
Low Froude Number Limit of the Rotating Shallow Water and Euler Equations Kung-Chien Wu Department of Pure Mathematics and Mathematical Statistics University of Cambridge, Wilberforce Road Cambridge, CB3
More informationTOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES
MATHEMATICS OF COMPUTATION Volume 67 Number 221 January 1998 Pages 73 85 S 0025-5718(98)00913-2 TOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES SIGAL GOTTLIEB AND CHI-WANG SHU Abstract. In this paper we
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationMathematical analysis of the stationary Navier-Stokes equations
Mathematical analysis of the Department of Mathematics, Sogang University, Republic of Korea The 3rd GCOE International Symposium Weaving Science Web beyond Particle Matter Hierarchy February 17-19, 2011,
More informationA Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations
A Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations Songul Kaya and Béatrice Rivière Abstract We formulate a subgrid eddy viscosity method for solving the steady-state
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 informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
More informationA RECURRENCE THEOREM ON THE SOLUTIONS TO THE 2D EULER EQUATION
ASIAN J. MATH. c 2009 International Press Vol. 13, No. 1, pp. 001 006, March 2009 001 A RECURRENCE THEOREM ON THE SOLUTIONS TO THE 2D EULER EQUATION Y. CHARLES LI Abstract. In this article, I will prove
More informationChapter 4. MAC Scheme. The MAC scheme is a numerical method used to solve the incompressible Navier-Stokes. u = 0
Chapter 4. MAC Scheme 4.1. MAC Scheme and the Staggered Grid. The MAC scheme is a numerical method used to solve the incompressible Navier-Stokes equation in the velocity-pressure formulation: (4.1.1)
More informationGlobal regularity of a modified Navier-Stokes equation
Global regularity of a modified Navier-Stokes equation Tobias Grafke, Rainer Grauer and Thomas C. Sideris Institut für Theoretische Physik I, Ruhr-Universität Bochum, Germany Department of Mathematics,
More informationVANISHING VISCOSITY IN THE PLANE FOR VORTICITY IN BORDERLINE SPACES OF BESOV TYPE
VANISHING VISCOSITY IN THE PLANE FOR VORTICITY IN BORDERLINE SPACES OF BESOV TYPE ELAINE COZZI AND JAMES P. KELLIHER Abstract. The existence and uniqueness of solutions to the Euler equations for initial
More informationVector penalty-projection methods for outflow boundary conditions with optimal second-order accuracy
Vector penalty-projection methods for outflow boundary conditions with optimal second-order accuracy Philippe Angot, Rima Cheaytou To cite this version: Philippe Angot, Rima Cheaytou. Vector penalty-projection
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationA Pressure-Correction Scheme for Rotational Navier-Stokes Equations and Its Application to Rotating Turbulent Flows
Commun. Comput. Phys. doi: 10.4208/cicp.301109.040310s Vol. 9, No. 3, pp. 740-755 March 2011 A Pressure-Correction Scheme for Rotational Navier-Stoes Equations and Its Application to Rotating Turbulent
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationA stencil of the finite-difference method for the 2D convection diffusion equation and its new iterative scheme
International Journal of Computer Mathematics Vol. 87, No. 11, September 2010, 2588 2600 A stencil of the finite-difference method for the 2D convection diffusion equation and its new iterative scheme
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationStable and high-order accurate finite difference schemes on singular grids
Center for Turbulence Research Annual Research Briefs 006 197 Stable and high-order accurate finite difference schemes on singular grids By M. Svärd AND E. van der Weide 1. Motivation and objectives The
More informationAn Equal-order DG Method for the Incompressible Navier-Stokes Equations
An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce
More informationNew Discretizations of Turbulent Flow Problems
New Discretizations of Turbulent Flow Problems Carolina Cardoso Manica and Songul Kaya Merdan Abstract A suitable discretization for the Zeroth Order Model in Large Eddy Simulation of turbulent flows is
More informationERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION. Xiaofeng Yang. (Communicated by Jie Shen)
DISCRETE AND CONTINUOUS doi:1.3934/dcdsb.29.11.157 DYNAMICAL SYSTEMS SERIES B Volume 11, Number 4, June 29 pp. 157 17 ERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION Xiaofeng Yang
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationConvergence and Error Bound Analysis for the Space-Time CESE Method
Convergence and Error Bound Analysis for the Space-Time CESE Method Daoqi Yang, 1 Shengtao Yu, Jennifer Zhao 3 1 Department of Mathematics Wayne State University Detroit, MI 480 Department of Mechanics
More informationConvergence of A Galerkin Method for 2-D Discontinuous Euler Flows Jian-Guo Liu 1 Institute for Physical Science andtechnology and Department of Mathe
Convergence of A Galerkin Method for 2-D Discontinuous Euler Flows Jian-Guo Liu 1 Institute for Physical Science andtechnology and Department of Mathematics University of Maryland College Park, MD 2742
More informationAN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION
AN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION Fareed Hussain Mangi*, Umair Ali Khan**, Intesab Hussain Sadhayo**, Rameez Akbar Talani***, Asif Ali Memon* ABSTRACT High order
More informationDecay in Time of Incompressible Flows
J. math. fluid mech. 5 (23) 231 244 1422-6928/3/3231-14 c 23 Birkhäuser Verlag, Basel DOI 1.17/s21-3-79-1 Journal of Mathematical Fluid Mechanics Decay in Time of Incompressible Flows Heinz-Otto Kreiss,
More informationVORTICITY LAYERS OF THE 2D NAVIER-STOKES EQUATIONS WITH A SLIP TYPE BOUNDARY CONDITION
VORTICITY LAYERS OF THE 2D NAVIER-STOKES EQUATIONS WITH A SLIP TYPE BOUNDARY CONDITION GUNG-MIN GIE 1 AND CHANG-YEOL JUNG 2 Abstract. We study the asymptotic behavior, at small viscosity ε, of the Navier-
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationThe Navier-Stokes Equations with Time Delay. Werner Varnhorn. Faculty of Mathematics University of Kassel, Germany
The Navier-Stokes Equations with Time Delay Werner Varnhorn Faculty of Mathematics University of Kassel, Germany AMS: 35 (A 35, D 5, K 55, Q 1), 65 M 1, 76 D 5 Abstract In the present paper we use a time
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationAvailable online at J. Math. Comput. Sci. 4 (2014), No. 3, ISSN:
Available online at http://scik.org J. Math. Comput. Sci. 4 (2014), No. 3, 587-593 ISSN: 1927-5307 A SMALLNESS REGULARITY CRITERION FOR THE 3D NAVIER-STOKES EQUATIONS IN THE LARGEST CLASS ZUJIN ZHANG School
More informationExistence, uniqueness and non-regularity of the solution to the Neumann problem for the mono-dimensional compressible Reynolds equation
Mathematical Communications 12(27), 95-1 95 Existence, uniqueness and non-regularity of the solution to the Neumann problem for the mono-dimensional compressible Reynolds equation Sanja Marušić Abstract.
More informationDUE: WEDS MARCH 26TH 2018
HOMEWORK # 2: FINITE DIFFERENCES MAPPING AND TWO DIMENSIONAL PROBLEMS DUE: WEDS MARCH 26TH 2018 NOTE: In this homework, you will choose ONE of the following three questions to perform and hand-in. Each
More informationLECTURE 5: THE METHOD OF STATIONARY PHASE
LECTURE 5: THE METHOD OF STATIONARY PHASE Some notions.. A crash course on Fourier transform For j =,, n, j = x j. D j = i j. For any multi-index α = (α,, α n ) N n. α = α + + α n. α! = α! α n!. x α =
More informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationAnalysis and Computation of Navier-Stokes equation in bounded domains
1 Analysis and Computation of Navier-Stokes equation in bounded domains Jian-Guo Liu Department of Physics and Department of Mathematics Duke University with Jie Liu (National University of Singapore)
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationISSUES RELATED TO LEAST-SQUARES FINITE ELEMENT METHODS FOR THE STOKES EQUATIONS
SIAM J. SCI. COMPUT. Vol. 20, No. 3, pp. 878 906 c 1998 Society for Industrial and Applied Mathematics ISSUES RELATED TO LEAST-SQUARES FINITE ELEMENT METHODS FOR THE STOKES EQUATIONS JENNIFER M. DEANG
More informationVector Penalty-Projection Methods for the Solution of Unsteady Incompressible Flows
Author manuscript, published in "5th International Symposium on Finite Volumes for Complex Applications, Aussois : France 2008)" Vector Penalty-Projection Methods for the Solution of Unsteady Incompressible
More informationControl of Interface Evolution in Multi-Phase Fluid Flows
Control of Interface Evolution in Multi-Phase Fluid Flows Markus Klein Department of Mathematics University of Tübingen Workshop on Numerical Methods for Optimal Control and Inverse Problems Garching,
More informationDomain decomposition schemes with high-order accuracy and unconditional stability
Domain decomposition schemes with high-order accuracy and unconditional stability Wenrui Hao Shaohong Zhu March 7, 0 Abstract Parallel finite difference schemes with high-order accuracy and unconditional
More informationOn the positivity of linear weights in WENO approximations. Abstract
On the positivity of linear weights in WENO approximations Yuanyuan Liu, Chi-Wang Shu and Mengping Zhang 3 Abstract High order accurate weighted essentially non-oscillatory (WENO) schemes have been used
More informationarxiv: v1 [math.ap] 14 Apr 2009
ILL-POSEDNESS OF BASIC EQUATIONS OF FLUID DYNAMICS IN BESOV SPACES arxiv:94.2196v1 [math.ap] 14 Apr 29 A. CHESKIDOV AND R. SHVYDKOY ABSTRACT. We give a construction of a divergence-free vector field u
More informationFourier analysis for discontinuous Galerkin and related methods. Abstract
Fourier analysis for discontinuous Galerkin and related methods Mengping Zhang and Chi-Wang Shu Abstract In this paper we review a series of recent work on using a Fourier analysis technique to study the
More information