Projection Method with Spatial Discretizations
|
|
- Amy Lloyd
- 5 years ago
- Views:
Transcription
1 Projection Method with Spatial Discretizations Weinan E School of Mathematics Institute for Advanced Study Princeton, NJ and Jian-Guo Liu Courant Institute of Mathematical Sciences New York, NY Abstract: In [5], we studied convergence and the structure ofthe error for several projection methods when the spatial variable was kept continuous 1. In this paper, we address similar problems for the fully discrete case when the spatial variable is discretized using the staggered grid. We prove that the numerical solution has full accuracy upto 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 ofconvergence and error structure for the projection method. In previous work [5], 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 ofchorin and Temam [?,?] is uniformly first order accurate upto the boundary, and the second order method of[?,?] is uniformly second order accurate upto the boundary. However, due to the presence ofnumerical boundary layers, pressure has only halforder ofaccuracy near the boundary. We characterized explicitly the structure ofthe numerical boundary layer, and proved that the full accuracy for pressure 1 We will refer to this situation as semi-discrete
2 can be recovered ifwe subtract the numerical boundary layer modes from the solution. We also showed using normal mode analysis that the boundary layer modes turn into oscillatory modes for the pressure-increment formulation of [?, 2]. For a summary ofthese results, we refer to [?]. The presence ofsingular modes in the semi-discrete solutions raise serious doubts on the convergence of the fully discrete method. Indeed it is known for some time that accuracy and even convergence can be lost ifcare is not exercised at the projection step. This is documented carefully in [?]. It is also known for a long time that full accuracy is kept if the spatial discretization is done on a staggered grid. This paper is devoted to a proofof this empirical fact. We work still leaves open the very important issue ofcharacterizing the minimum condition that 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 self-adjoint; (2) Ph 2 = 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. 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 does not require the projection condition. The gauge method [?] 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 [?]). In Section 3, we summarize our main results. Then in Sections 4 and 5 we present the prooffor the first and second order methods respectively. 2. Review of the Projection Methods In primitive variables, Navier-Stokes equation (NSE) takes the following form t u +(u )u + p = u, (2.1) u =0. 2
3 Here u =(u, v) is the velocity, and p is the pressure. For simplicity, we will only consider the case when the no-slip boundary condition 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[?,?] proceeds in two steps. Step 1: Computing the intermediate velocity field u. u u n +(u n )u n = u, t (2.6) u =0, on Ω. Step 2: Projecting to the space ofdivergence-free vector fields to obtain u n+1. u = u n+1 + t p n+1, (2.7) u n+1 =0. The projection step enforces the boundary condition: (2.10) u n+1 n =0, or Second order schemes: p n+1 n =0, on Ω We will concentrate on the formulation in [?,?], and refer to [?] for a summary of other second order methods. u u n t +(u n+1/2 )u n+1/2 = u + u n, 2 (2.12) u + u n = t p n 1/2, on Ω, u = u n+1 + t p n+1/2, u n+1 =0, p n+1/2 n =0, on Ω. In this formulation, homogeneous Neumann BC for pressure is retained. An inhomogeneous BC for u is introduced so that the slip velocity of u n+1 at the boundary is oforder t 2. 3
4 Remark The nonlinear convection term (u n+1/2 )u n+1/2 can be treat in many ways. In the Theorems 1 and 2, we use the Adams-Bashforth formula: 3 2 (un )u n 1 2 (un 1 )u n 1. It is readily seen that the projection step enforces (2.13) p n+1 n = pn n = = p0 n =0, on Ω for the numerical solution. In general this is not satisfied by the exact solution of (2.1). Therefore we expect that pn n has O(1) error at the boundary. As will be seen below, 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. j+1 j v i,j+1/2 p i,j u i+1/2,j j 1 i 1 i i+1 Figure 1: The MAC mesh An illustration ofthe 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 Other differential operators are discretized as: + v i,j+1/2 v i,j 1/2 y. ( 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, 4
5 ( 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 + v i+1/2,j a i+1/2,j+1 a i+1/2,j 1 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 + v i,j+1/2 b i,j+3/2 b i,j 1/2 2 y Clearly the truncation errors ofthese approximations are ofsecond 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. Similarly the boundary condition v = 0 is imposed at the horizontal physical boundary, but u = 0 is imposed at the ghost triangle points with a distance of y 2 away from the physical boundary. One shortcoming ofthe MAC scheme is the serious constraint on geometry. Although slightly more general situations can be studied, in the present paper we will on the situation when Ω = [ 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 Ω to denote the part ofthe boundary at x = ±1. We will always assume that x y and h =min( x, y). Notations: We will use C to denote generic constants which may depend on the norms ofthe exact solutions. Norms will be taken over the entire domain Ω. 3. Summary of Results and Outline of Proofs The main results ofthis paper are the following (the constants are independent of t and h): 5
6 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.13) u 0 (x) =0, y p(x, 0) = 2 xy p(x, 0) = 0, on Ω, Let (u h,p h ) be the numerical solution of the projection method (2.6), (2.7) and (2.10) coupled with the MAC spatial discretization. Assume that t <<h. Then we have (3.14) u u h L + t 1/2 p p h L C( t + h 2 ), (3.15) p p h p c L C( t + h 2 ), where (3.16) p c (x,t) t 1/2 β eα e α 1 e α x 1 / t1/2 D x + p h(x t 1/2,y,t) eα + t 1/2 β e α x+1 / t1/2 e α D+ x 1 p h(x + t 1/2,y,t), α = t1/2 2 x arccosh (1 + x 2 t ), β = x t 1/2 (1 e α x/ t1/2 ) 1. 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.17) α 1 x α 2 y u 0 (x) =0, on Ω, for α 1 + α 2 6, Let (u h,p h ) be the numerical solution of the projection method (2.12) coupled with the MAC spatial discretization. Assume that t 2 << h. Then we have (3.18) u u h L + t 3/2 p p h L + t p p h L (0,T,L 2 ) C( t 2 + h 2 ). (3.19) p p h p c L C( t + h 2 ), where (3.20) p c t 1/2 β eα e α 1 e α x 1 / t1/2 D x +p h (x t 1/2,y,t) eα + t 1/2 β e α x+1 / t1/2 e α D+ x 1 p h(x + t 1/2,y,t), 6
7 α = t1/2 2 x arccosh (1 + x t ), β = x (1 e α x/ t1/2 ) 1. 1/2 t Remark. The constraints on the size of t is only technical and is different from the standard stability condition. The proofofthese results follow the general strategy outlined in Section 3 of[5]. 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, t (4.1) u =0, on Ω, u = u n+1 + t 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 N N N ((a, c)) = x y a i+1/2,j c i+1/2,j + x y b i,j+1/2 d i,j+1/2 i=1 j=1 i=1 j=1 N 1 N N N ((u, h p)) = y u i+1/2,j (p i+1,j p i,j )+ x v i,j+1/2 (p i,j+1 p i,j ) i=1 j=1 i=1 j=1 N 1 N N N (( h u,p)) = y (u i+1/2,j u i 1/2,j )p i,j + x (v i,j+1/2 v i,j 1/2 )p i,j i=1 j=1 i=1 j=1 and discrete norms (4.3) u =((u, u)) 1/2, u =max u i,j i,j Denote h =min( x, y). Lemma 4.1 We have the following 7
8 (i) Inverse inequality: (4.4) f 1 h f, (ii) Poincare inequality: suppose f x=±1 =0,then (4.5) f h f, (iii) Suppose n u x=±1 =0,then we have (4.6) ((u, h p)) = (( h u,p)) (iv) Suppose u x=±1 =0,then we have (4.7) 2((u, h u)) h u 2 h u 2 (v) Suppose a x=±1 =0and c n x=±1 =0, then we have (4.8) ((a, N h (u, c))) c h a u W 1, Proof : The proofof(i iii) is standard. We first show (iv). Summation by parts gives (4.9) ((u, h u)) = 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 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 ] (4.12) = h u 2 +4 j [(v 0,j+1/2 ) 2 +(v N+1,j+1/2 ) 2 ] 8
9 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 0c 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) Summation by parts gives I = x y i,j c i+1/2,j [D x 0 (u i+1/2,j a i+1/2,j )+D y 0 ( v i+1/2,ja i+1/2,j )] x y i,j d i,j+1/2 [D x 0 (ū i,j+1/2b i,j+1/2 )+D y 0 (v i,j+1/2b i,j+1/2 )] (4.14) 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. Again we set ε = t 1/2, ξ =(x+1)/ε, x i = 1+i x, ξ i = i ξ, ξ = x/ε, t n = n t, t n 1/2 =(n 1/2) t, 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+ ξ = εd+ x. We will use the notation (4.17) D 2 ξ = D ξ D ξ +, D 2 y = D y D y +, and (4.18) ξ =(D ξ +, 0), y =(0,D y +), 9
10 Denote the solutions of(4.1) as (u h, u h,p h). Motivated by the discussions in 4, we make the following ansatz, valid at t n = n t, 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 = h u 0 + j=2 ε j ( h u j + ε 2 D 2 ξ a j + D 2 y a j), (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 ofequations 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 )= 2 h u 0 + 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 + Dy 2 a j. k=0 10
11 The second equation in (4.1) implies (4.26) u 2 + a 2 = u 2 + t u 0 + h p 0 + ξ ϕ 1 + y ϕ 0. For j =2l 1, l 1, (4.27) u j+2 + a j+2 = u j+2 + t u j + h p j + ξ ϕ j+1 + y ϕ j + For j =2l, l 1, (4.28) u j+2 +a j+2 = u j+2+ t u j + h p j + ξ ϕ j+1 + y ϕ j + 1 (l +1)! u(l+1) 0 + From the third equation in (4.1), we obtain l k=2 1 k! u(k) j 2k+2. l k=2 1 k! u(k) j 2k+2. (4.29) h u j =0, j =0, 1,. The boundary conditions become (4.30) u j + a j =0, D x +p j 1 + D ξ +ϕ j =0, at x = 1, ξ =0, for j>0. Next we go through all these equations, order by order, to see ifthey are solvable. Since this is very similar to what we did in 4.1, we will only give a summary ofresults. The coefficients in the expansions (4.19) can be obtained successively in the following order: (4.31) 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 ( ) Using the following lemma, we know that (4.31) has a smooth solution in the sense that the divided difference ofvarious orders are bounded. The lemma itself, as well as Lemma 4.3, belongs to the folklore of classical numerical analysis. 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). 11
12 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 We next have (4.33) u 2 = u 2 + t u 0 + h p 0, (4.34) ϕ 1 = D 2 ξ ϕ 1, D ξ +ϕ 1 ξ=0 = D x +p 0 x= 1, This gives (4.35) ϕ 1 (ξ,y,t) =βd x +p 0 ( 1,y,t) e αξ. where (4.36) α = 1 ξ arccosh (1 + ξ 2 /2), β = ξ(1 e α ξ ) 1. (4.37) a 2 = D ξ +ϕ 1, b 2 =0, (4.38) u 3 = u 3, (4.39) ϕ 2 =0, a 3 =0, b 3 = Dy +ϕ 1, (4.40) t u 2 + h p 2 + N h (u 0, u 2 )+N h (u 2, u 0 ) h u 2 =0, = h u 2 + h ( t u 0 + h p 0 ) t u 0, u 2 x= 1 = h p 0 x= 1 ξ ϕ 1 ξ=0, on Ω With a suitable initial data, we know from the following lemma that (4.40) has a smooth solution. 12
13 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, (4.41) h u =0, u = g, at x = ±1, u(, 0) = u 0 ( ) where f, g and u 0 smooth and satisfies some compatibility conditions. Then (u,p) is smooth in the sense that its divided differences of various order are bounded. (4.42) Continue in this fashion, we get The solution for (4.42) is ϕ 3 = D 2 ξ ϕ 3 + D 2 y ϕ 1, D ξ +ϕ 3 ξ=0 = D x +p 2 x= 1, (4.43) ϕ 3 (ξ,y,t) =βd x +p 2 e αξ + β 1 (ξ + γ)d x +D 2 y p 0 x= 1 e αξ. where (4.44) β 1 = 1 (1 e α ξ )(e α ξ e α ξ ), γ = ξ e α ξ 1 e α ξ (4.45) b 4 =0, a 4 = D ξ +ϕ 3. (4.46) 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.47) ϕ 4 = D 2 ξ ϕ 4, D ξ +ϕ 4 ξ=0 = D x + p 3 x= 1, (4.48) a 5 = D ξ +ϕ 4, b 5 = D y +ϕ 3. 13
14 Obviously this procedure can be continued and we obtain (4.49) ϕ j = D 2 ξ ϕ j + D 2 y ϕ j 2, D ξ +ϕ j ξ=0 = D x + p j 1 x= 1, (4.50) ϕ j = [j/2] k=0 F j,k (y)ξ k e αξ, (4.51) a j = D ξ +ϕ j 1, b j = D y +ϕ j 2, Now ifwe let (4.52) U = u 0 + U n = u 0 + P n = p 0 + 2N j=1 2N j=1 2N j=1 ε j (u j + a j ), ε j u j, ε j (p j + ϕ j )+ε 2N+1 ϕ 2N+1, then we have U U n t + N h (U n,u n )= h U + t α f, (4.53) U =0, at x = ±1, U = U n+1 + t h P n + t α+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 0 t T U n ( ) W 1, C. For the initial data, we have (4.55) U 0 (x) =u 0 (x)+ tw 0 (x) 14
15 where w 0 is a bounded function. Furthermore under the compatibility condition (3.13), we can construct a better approximate initial data (4.57) U 0 (x) =u 0 (x)+ t 2 w 0 (x). Proof of Theorem 1: Assume a priori that (4.58) max 0 t n T un W 1, C. In the following estimates, the constant will sometimes depend on C and C. 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 + t α f n, t e =0, at x = ±1, (4.60) e n+1 e t + h q n = t α g n, h e n+1 =0, D x +q n = e n+1 n =0, at x = ±1, e 0 = t α w 0. Taking the scalar product ofthe first equation of(4.60) with 2e and integrating by parts, we obtain e 2 e n 2 + e e n 2 + t h e 2 (4.61) t 2α+1 f n 2 + t e 2 2 t ((e, N h (e n,u n ))) 2 t ((e, N h (u n, e n ))) t 2α+1 f n 2 + C t( e 2 + e n 2 )+ 1 2 t he 2. Here we have used Lemma 4.1. Taking the scalar product ofthe second equation of(4.60) with 2e n+1 yields (4.62) e n+1 2 e 2 + e n+1 e 2 t e n t 2α+1 g n 2. 15
16 Combining (4.61) and (4.62), we get (4.63) e n+1 2 e n 2 + e e n 2 + e n+1 e 2 + t h e 2 C t ( e n e n 2 )+ t 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 + t 1/2 h e C 1 t α. Using the second equation of(4.60) we have (4.65) e n + t h q n C 1 t α Now by inverse inequality (4.4) we have (4.66) e n L + t h q n L + h e n W 1, C 1 t α h. Chose N =3and t α << h 2,ifwechoose t 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, +1 n [ T t which depends only on the exact solution (u,p). This proves (4.69) u 0 u h L + p 0 p h L 2 + t 1/2 p 0 p h L + p 0 p h p c L C t 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 + t 1/2 p p h L + p p h p c L C( t + h 2 ) This completes the proofoftheorem 1. 16
17 5. Second Order Schemes with Spatial Discretization In this section we carry out the same program as in 4 for the second order method (2.12) with the standard MAC spatial discretization: u u n u + u n = h, t 2 (5.1) u + u n = t h p n 1/2, at x = ±1, u = u n+1 + t 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 complicates substantially the presentation. We begin with the following ansatz: u (x) =u 0 (x,t n )+ ε j [u j(x,t n )+a j(ξ,y,t n )], 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 ε = t 1/2, ξ =(x +1)/ε, t n = n t, t n 1/2 =(n 1/2) t, 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 ). 17
18 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) u 3 + a 3 = u 3 + t u 1 + h p 1 + ξ ϕ 2 + y ϕ 1. For j =2l, (5.7) u j+2 + a j+2 = u j+2 + t u j + h p j + ξ ϕ j+1 + y ϕ j + 1 (l +1)! u(l+1) l k=1 l k=2 1 2 k k! ( ξϕ (k) j 2k+1 + yϕ (k) j 2k ). 1 k! u(k) j 2k l l! hp (l) l k k! hp (k) j 2k k=1 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 k! u(k) j 2k+2 + k=2 l k=1 1 2 k k! ( hp (k) j 2k + ξϕ (k) j 2k+1 + yϕ (k) j 2k ). From the incompressibility condition, we get (5.9) h u j =0, for j 0. The boundary conditions imply that for x = 1, ξ =0, (5.10) u 0 =0, (5.11) u 2 + u 2 + a 2 = h p 0 + ξ ϕ 1, (5.12) u 3 + u 3 + a 3 = h p 1 + ξ ϕ 2 + y ϕ 1 ; 18
19 for j =2l, l 1 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 (5.13) l 2 + k=1 ( 1) k 2 k k! hp (k) l 1 j 2k 2 + k=1 ( 1) k 2 k k! ( ξϕ (k) j 2k 1 + yϕ (k) j 2k 2 ) for j =2l +1, l 1 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, order by order, to see ifthey 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) u 2 = u 2 + t u 0 + h p 0, (5.18) ϕ 1 = 1 2 D 2 ξ ϕ 1, D ξ +ϕ 1 ξ=0 = D x + p 0 x= 1, (5.19) ϕ 1 = βd x + p 0 x= 1 e αξ, where (5.20) α = 1 ξ arccosh (1 + ξ 2 ), β = ξ(1 e α ξ ) 1. 19
20 (5.21) a 2 = D ξ +ϕ 1, b 2 =0, We next have: (5.22) u 3 = u 3, (5.23) ϕ 2 =0, a 3 =0, b 3 = D y +ϕ 1, (5.24) ϕ 3 = 1 2 (D 2 ξ ϕ 3 + D 2 y ϕ 1), D ξ +ϕ 3 ξ=0 =0, The solution for (5.24) is (5.25) ϕ 3 (y,ξ,t) =β 1 (ξ + γ)d x +D 2 y p 0 x= 1 e αξ. where (5.26) β 1 = 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) ϕ 4 =0, a 5 =0, b 5 = 1 2 Dy + tϕ 1 + D y + ϕ 3, (5.30) u 5 = u 5, (5.31) 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. (5.32) Now ifwe let 2N U = u 0 + ε j (u j + a j ), U n = u 0 + j=1 2N j=1 ε j u j, P n 1/2 = p 0 + 2N j=1 ε j (p j + ϕ j )+ε 2N+1 ϕ 2N+1, 20
21 then we have U U n t = h U + U n 2 + t α f, (5.33) U + U n = t h P n 1/2, at x = ±1, U = U n+1 + t h P n+1/2 + t α+1 g, h U n+1 =0, D x +P n+1/2 = 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 (5.34) max 0 t T U n ( ) W 1, C, For the initial approximation, we have (5.35) U 0 (x) =u 0 (x)+ t 2 w 0 (x) without the extra compatibility condition, and (5.36) U 0 (x) =u 0 (x)+ t 4 w 0 (x) with the compatibility condition (3.17). Proof of Theorem 2. Assume a priori that (5.37) max 0 t n T un W 1, C. As in the proofoftheorem 1, we let (5.38) e n = U n u n, e = Û û, q n = P n 1/2 p n 1/2. where (5.39) 2û = u + u n t h p n 1/2, 2Û = U + U n t h P n 1/2. 21
22 From (5.1) and (5.34), we get 2(e e n ( ) + h q n 1 ) t 2 t 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 )+ t α f n, (5.40) e =0, at x = ±1, e n+1 + e n 2e + h (q n+1 q n )= t α g n, t h e n+1 =0, D+ x qn+1/2 = e n+1 n =0, at x = ±1, e 0 = t α w 0. Taking the scalar product ofthe first equation of(5.40) with e and integrating by parts, we get (5.41) e 2 e n 2 + e e n t e t e 2 t e (q n 1 Ω 2 t qn ) dx + C t 2α+1 f n 2 +C t ( e n 2 + e n e 2 )+ 1 2 t e 2. Taking the scalar product ofthe 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 t 2α+1 g n 2 + C t e n+1 2. Combining the these two estimates, we get e n+1 2 e n 2 + e n+1 + e n 2e 2 + t h e 2 + t h e 2 (5.43) 2 t ((e, h (q n 1 2 t hq n ))) + C t ( e n 2 + e n e n+1 2 ) +C t 2α+1 ( f n 2 + g n 2 ). To estimate the first term on the right hand of(5.45), we let (5.44) I 2 t((e, h (q n 1 2 t hq n ))) = 2 t((e, h q n )) t 2 (( h e, h q n )) I 1 + I 2. 22
23 Using the second equation and integrating by parts, we can write the first term as I 1 = 2 t((e, h q n )) (5.45) = t 2 (( h (q n+1 q n ), h q n )) t α+2 ((g n, h q n )) = 1 2 t 2 ( h q n+1 2 h q n 2 ) t 2 h (q n+1 q n ) 2 t α+2 ((g n, h q n )). Since (5.46) 1 2 t 2 h (q n+1 q n ) 2 = 1 2 en+1 + e n 2e tα+1 g n 2 + t α+1 ((g n, e n+1 + e n 2e )). We have (5.47) I 1 = 1 2 t 2 ( h q n+1 2 h q n 2 )+ 1 2 en+1 + e n 2e tα+1 g n 2 + t α+1 ((g n, e n+1 + e n 2e )) t α+2 ((g n, h q n )). Next we rewrite the second term as I 2 = t 2 (( h e, h q n )) = 1 2 t3 (( h (q n+1 q n ), h q n )) 1 2 tα+3 (( h g n, h q n )) (5.48) = 1 4 t3 ( h q n+1 2 h q n 2 )+ 1 4 t3 h (q n+1 q n ) tα+3 (( h g n, h q n )) = 1 4 t3 ( h q n+1 2 h q n 2 )+ t h e t2α+3 h g n 2 t α+2 (( h g n, h e )) 1 2 tα+3 (( h g n, h q n )). 23
24 Combining these two terms we arrive at I = 1 2 t 2 ( h q n+1 2 h q n 2 ) 1 4 t3 ( h q n+1 2 h q n 2 ) (5.49) en+1 + e n 2e 2 + t h e 2 + t α+1 ((g n, e n+1 + e n 2e )) t α+2 ((g n, h q n )) t α+2 (( h g n, h e )) 1 2 tα+3 (( h g n, h q n )) t2α+3 h g n tα+1 g n 2. This gives I 1 2 t 2 ( h q n+1 2 h q n 2 ) 1 4 t3 ( h q n+1 2 h q n 2 ) (5.50) en+1 + e n 2e 2 + t h e 2 + t e n+1 + e n 2e 2 +2 t 3 h q n 2 +2 t 4 h q n 2 +2 t 2α+1 ( g n 2 + t g n 2 H 1 ). Going back to (5.45) we obtain e n+1 2 e n en+1 + e n 2e 2 + t h e 2 (5.51) t 2 ( h q n+1 2 h q n 2 )+ 1 4 t3 ( h q n+1 2 h q n 2 ) t 3 h q n 2 + t 4 h q n 2 + C t ( e n 2 + e n e n+1 2 ) +C t 2α+1 ( f n 2 + t g n 2 H 1 ). Gronwall lemma gives (5.52) e n + e + t h q n + t 3/2 h q n + t 1/2 h e C 1 t α. Now by inverse inequality (4.4) we have (5.53) e n L + h e n W 1, + t h q n L C 1 t α h. Chose N =5and t α << h 2,ifwechoose t small enough, we will always have (5.54) e n+1 L 1. 24
25 Therefore in (5.39) we can choose (5.55) C =1+ max [ ] U n ( ) T W 1, n t which depends only on the exact solution (u,p). This proves (5.56) u 0 u h L + p 0 p h L 2 + t 1/2 p 0 p h L + p 0 p h p c L C t 2 From Lemma 4.2, we have (5.57) u u h L + p p h L 2 + t 1/2 p p h L + p p h p c L C( t 2 + h 2 ) This completes the proofoftheorem 2. References [1] C.R. Anderson, Derivation and solution of the discrete pressure equations for the incompressible Navier-Stokes equations, preprint. [2] 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), [3] A.J. Chorin, Numerical solution of the Navier-Stokes equations, Math. Comp. 22 (1968), [4] A.J. Chorin, On the convergence of discrete approximations to the Navier-Stokes equations, Math. Comp. 23 (1969), [5] Weinan E and J.-G. Liu, Projection method I: convergence and numerical boundary layers, SIAMJ. Numer. Anal. 32 (1995), ; Projection method II: Godunov- Ryabenki analysis, SIAMJ. Numer. Anal. 33 (1996), [6] P.M. Gresho, Incompressible fluid dynamics: some fundamental formulation issues, Rev. Fluid Mech. 23 (1991),
26 [7] P.M. Gresho, Some interesting issues in incompressible fluid dynamics, both in the continuum and in numerical simulation, Adv. in Appl. Mechanics 28 (1992), [8] 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, SIAMJ. Numer. Anal., 25, , [9] G. Karniadakis, M. Israeli and S. A. Orszag, High-order splitting methods for the incompressible Navier-Stokes Equations, J. Comput. Phys., 97, No. 2, , [10] J. Kim and P. Moin, Application of a fractional-step method to incompressible Navier- Stokes equations, J. Comp. Phys. 59 (1985), [11] D. Michelson, Convergence theorems for finite difference approximations for hyperbolic quasilinear initial-boundary value problems, Math. Comp. 49 (1987), [12] H. Okamoto, On the semi-discrete finite element approximation for the nonstationary Navier-Stokes equation, J. Fac. Sci. Univ. Tokyo Sect. IA Math., 29 (1982), [13] S.A. Orszag, M. Israel, and M.O. Deville, Boundary conditions for incompressible flows, J. Scientific Computing 1 (1986), [14] Jie Shen, On error estimates of projection methods for Navier-Stokes equations: first order schemes, SIAM J. Numer. Anal , [15] Jie Shen, On error estimates of some higher order projection and penalty-projection methods for Navier-Stokes equations, Numer. Math. 62 (1992), [16] G. Strang, Accurate partial differential methods II. non-linear problems, Numerische Mathematik 6 (1964), [17] 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), [18] R. Temam, Navier Stokes Equations, North Holland, New York,
27 [19] R. Temam, Remark on the pressure Boundary condition for the projection method, Theoretical and Computational Fluid Dynamics 3 (1991), [20] J. van Kan, A second-order accurate pressure-correction scheme for viscous incompressible flow, SIAM J. Sci. Stat. Comp. 7 (1986),
PROJECTION METHOD III: SPATIAL DISCRETIZATION ON THE STAGGERED GRID
MATHEMATICS OF COMPUTATION Volume 71, Number 237, Pages 27 47 S 0025-5718(01)01313-8 Article electronically published on May 14, 2001 PROJECTION METHOD III: SPATIAL DISCRETIZATION ON THE STAGGERED GRID
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
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 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 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 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 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 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 informationESTIMATES OF LOWER ORDER DERIVATIVES OF VISCOUS FLUID FLOW PAST A ROTATING OBSTACLE
REGULARITY AND OTHER ASPECTS OF THE NAVIER STOKES EQUATIONS BANACH CENTER PUBLICATIONS, VOLUME 7 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 25 ESTIMATES OF LOWER ORDER DERIVATIVES OF
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 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 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 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 informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
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 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 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 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 informationA Multiphysics Strategy for Free Surface Flows
A Multiphysics Strategy for Free Surface Flows Edie Miglio, Simona Perotto, and Fausto Saleri MOX, Modeling and Scientific Computing, Department of Mathematics, Politecnico of Milano, via Bonardi 9, I-133
More informationThe Helically Reduced Wave Equation as a Symmetric Positive System
Utah State University DigitalCommons@USU All Physics Faculty Publications Physics 2003 The Helically Reduced Wave Equation as a Symmetric Positive System Charles G. Torre Utah State University Follow this
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 informationA study of the modelling error in two operator splitting algorithms for porous media flow
Computational Geosciences 2 (1998) 23 36 23 A study of the modelling error in two operator splitting algorithms for porous media flow K. Brusdal a,h.k.dahle a, K. Hvistendahl Karlsen a and T. Mannseth
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 informationA New Regularity Criterion for the 3D Navier-Stokes Equations via Two Entries of the Velocity Gradient
Acta Appl Math (014) 19:175 181 DOI 10.1007/s10440-013-9834-3 A New Regularity Criterion for the 3D Navier-Stokes Euations via Two Entries of the Velocity Gradient Tensor Zujin Zhang Dingxing Zhong Lin
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 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 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 informationA Projection FEM for Variable Density Incompressible Flows
Journal of Computational Physics 65, 67 88 (000) doi:0.006/jcph.000.6609, available online at http://www.idealibrary.com on A Projection FEM for Variable Density Incompressible Flows J.-L. Guermond and
More informationNumerical Analysis and Methods for PDE I
Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling
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 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 informationDIRECTION OF VORTICITY AND A REFINED BLOW-UP CRITERION FOR THE NAVIER-STOKES EQUATIONS WITH FRACTIONAL LAPLACIAN
DIRECTION OF VORTICITY AND A REFINED BLOW-UP CRITERION FOR THE NAVIER-STOKES EQUATIONS WITH FRACTIONAL LAPLACIAN KENGO NAKAI Abstract. We give a refined blow-up criterion for solutions of the D Navier-
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 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 informationImplementation of 3D Incompressible N-S Equations. Mikhail Sekachev
Implementation of 3D Incompressible N-S Equations Mikhail Sekachev Navier-Stokes Equations The motion of a viscous incompressible fluid is governed by the Navier-Stokes equations u + t u = ( u ) 0 Quick
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 informationEssentially Compact Schemes for Unsteady Viscous Incompressible Flows
Essentially Compact Schemes for Unsteady Viscous Incompressible Flows Weinan E 1 School of Mathematics Institute for Advanced Study Princeton, NJ 08540 Jian-Guo Liu 2 Department of Mathematics Temple University
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 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 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 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 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 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 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 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 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 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 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 informationAn N-order Iterative Scheme for a Nonlinear Wave Equation Containing a Nonlocal Term
Filomat 3:6 7) 755 767 DOI.98/FIL76755N Published by Faculty of Sciences and Mathematics University of Niš Serbia Available at: http://www.pmf.ni.ac.rs/filomat An N-order Iterative Scheme for a Nonlinear
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 informationNumerical dispersion and Linearized Saint-Venant Equations
Numerical dispersion and Linearized Saint-Venant Equations M. Ersoy Basque Center for Applied Mathematics 11 November 2010 Outline of the talk Outline of the talk 1 Introduction 2 The Saint-Venant equations
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 informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM ERIK BURMAN AND BENJAMIN STAMM Abstract. We propose a low order discontinuous Galerkin method for incompressible flows. Stability of the
More information