Research Article A Stabilized Low Order Finite-Volume Method for the Three-Dimensional Stationary Navier-Stokes Equations

Transcription

1 Mathematical Problems in Engineering Volume 212, Article ID , 14 pages doi:1.1155/212/ Research Article A Stabilized Low Order Finite-Volume Method for the Three-Dimensional Stationary Navier-Stokes Equations Jian Li, 1, 2 Xin Zhao, 2 Jianhua Wu, 1 and Jianhong Yang 2 1 College of Mathematics and Information Science, Shaanxi Normal University, Xi an 7162, China 2 Department of Mathematics, Baoji University of Arts and Sciences, Baoji 7217, China Correspondence should be addressed to Jian Li, jiaaanli@gmail.com Received 13 July 211; Accepted 28 October 211 Academic Editor: Rafael Martinez-Guerra Copyright q 212 Jian Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. This paper proposes and analyzes a stabilized finite-volume method FVM for the threedimensional stationary Navier-Stokes equations approximated by the lowest order finite element pairs. The method studies the new stabilized FVM with the relationship between the stabilized FEM FEM and the stabilized FVM under the assumption of the uniqueness condition. The results have three prominent features in this paper. Firstly, the error analysis shows that the stabilized FVM provides an approximate solution with the optimal convergence rate of the same order as the usual stabilized FEM solution solving the stationary Navier-Stokes equations. Secondly, superconvergence results on the solutions of the stabilized FEM and stabilized FVM are derived on the H 1 -norm and the L 2 -norm for the velocity and pressure. Thirdly, residual technique is applied to obtain the L 2 -norm error for the velocity without additional regular assumption on the exact solution. 1. Introduction Recently, the development of stable mixed FEMs is a fundamental component in the search for the efficient numerical methods for solving the Navier-Stokes equations governing the flow of an incompressible fluid by using a primitive variable formulation. The object of this work is to analyze the stabilized finite volume method for solving the three-dimensional stationary Navier-Stokes equations. The importance of ensuring the compatibility of the component approximations of velocity and pressure by satisfying the so-called inf-sup condition is widely understood. The numerous mixed finite elements satisfying the inf-sup condition have been proposed over the years. However, elements not satisfying the inf-sup condition may also work well. So far, the most convenient choice of the finite element space from an implementational point of view

2 2 Mathematical Problems in Engineering would be the elements of the low polynomial order in the velocity and the pressure with an identical degree distribution for both the velocity and the pressure. This paper focuses on the stabilized method called local polynomial pressure projection for the three-dimensional Navier-Stokes equations 1 5. The proposed method is characterized by the following features. First, the method does not require approximation of derivatives, specification of mesh-dependent parameters, edge-based data structures, and a nonstandard assembly procedure. Second, this method is completely local at the element level. On the other hand, FVM has become an active area in numerical analysis. The most attractive things are that FVM can keep local conservation and have the advantages of FVM and finite difference methods. The FVM is also termed the control volume method, the covolume method, or the first-order generalized difference method. Nowadays, it is difficult in analyzing FVM to obtain L 2 -norm error estimates because trail functions and test functions are derived from different spaces. Many papers were devoted to its error analysis for secondorder elliptic and parabolic partial differential problems 6 1. Error estimates of optimal orderintheh 1 -norm are the same as those for the linear FEM 9, 11. Error estimates of optimal order in the L 2 -norm can be obtained as well 8, 9. Moreover, the FVM for generalized Stokes problems was studied by many people They analyzed this method by using a relationship between it and the FEM and obtained its error estimates through those known for the latter method. Also, it still requires H 3 smoothness assumption of the exact solution to obtain O h 2 error bound in most previous literatures. However, for the Stokes problems only the finite elements that satisfy the discrete inf-sup condition have been studied. The work of 14, 15 for the two-dimensional stationary Stokes equations is extended in this paper for the three-dimensional stationary Navier-Stokes equations approximated by lowest equal-order finite elements. Following the abstract framework of the relationship between the stabilized FEM and stabilized FVM 14, 15, the stabilized FVM is studied, and the optimal error estimate of the stabilized FVM is obtained for the three-dimensional stationary Navier-Stokes equations relying on the uniqueness condition. As far as known, there still requires much research on FVM results 16 about the velocity in L 2 -norm and superconvergence result between FEM solution and FVM solution of the three-dimensional Navier-Stokes equations. The remainder of the paper is organized as follows. In Section 2, an abstract functional setting of the three-dimensional Navier-Stokes problem is given with some basic assumptions. In Section 3, the stability of the stabilized FVM is analyzed and provided by Brouwer s fixed-point theorem. In Section 4, the optimal error estimates of the stabilized finite volume approximation for the three-dimensional stationary Navier-Stokes equations are obtained. 2. FVM Formulation Let Ω be a bounded domain in R 3, assumed to have a Lipschitz-continuous boundary Γ and to satisfy a further condition stated in A1 below. The three-dimensional stationary Navier- Stokes equations are considered as follows: νδu p u u f, in Ω, 2.1a div u, in Ω, 2.1b u Ω, on Ω, 2.1c

3 Mathematical Problems in Engineering 3 where ν> is the viscosity, u u 1 x,u 2 x,u 3 x represents the velocity vector, p p x the pressure, and f f 1 x,f 2 x,f 3 x the prescribed body force. In order to introduce a variational formulation, we set 17 X ] 3, [ ] 3, } [H 1 Ω Y L 2 Ω M L 2 {q Ω L 2 Ω : qdx, Ω 2.2 [ ] 3 D A H 2 Ω X. As mentioned above, a further assumption on Ω is presented. A1 Assume that Ω is regular so that the unique solution v, q X, M of the steady Stokes problem Δv q g, div v inω, v Ω, 2.3 for a prescribed g Y exists and satisfies v 2 q 1 c g, 2.4 where c>isageneral constant depending on Ω. Here and after, i and i denote the usual norm and seminorm of the Sobolev space H i Ω or H i Ω 3 for i, 1, 2. We denote by, the inner product on L 2 Ω or Y. The space H 1 Ω and X are equipped with their equivalent scalar product and norm 17 u, v u, v, u u, u 1/ It is well known 18 that for each v X there hold the following inequalities: v L 4 2 1/2 v 1/4 v 3/4, v γ v 1, 2.6 where γ is a positive constant depending only on Ω. The continuous bilinear form a, on X X and d, on X M, respectively, are defined by a u, v u, v, u, v X, d v, q ) v, q ) q, div v v X, q M. 2.7 Also, the trilinear term is defined by b u, v, w u v, w 1 div u v, w u v, w u w, v, u, v, w X

4 4 Mathematical Problems in Engineering and satisfies b u, v, w c u v w. 2.9 Then the mixed variational form of 2.1a 2.1c is to seek u, p X, M such that a u, v d v, p ) d u, q ) b u, u, v f, v v, q ) X M. 2.1 The existence and uniqueness results are classical and can be found in We introduce the finite-dimensional subspace X h,m h X, M, which is characterized by τ h with mesh scale h, a partitioning of Ω into tetrahedron or hexahedron, assumed to be regular in the usual sense see Here, the space X h,m h satisfies the following approximation properties. For each v D A, p H 1 Ω, there exist approximations I h v X h and J h q M h such that u I h u h u I h u p Jh p ) ch 2 u 2 p together with the inverse inequality v h c 1 h 1 v h, u h L c 2 h 1/2 u h The stable and accurate finite element approximational solution of 2.1 requires that X h,m h satisfies the discrete inf-sup condition d sup v h X h ) v h β qh, 2.13 where β is positive constant independent of h. The main purpose of this paper is to study a stabilized FVM for the stationary 3D Navier-Stokes equations. We follow 23, 24 to obtain the dual partition K h. We first choose an arbitrary point Q in the interior of each tetrahedron K and then connect Q with the barycenters Q ijk of its 2D faces ΔP i P j P k by straight lines see Figure 1. On each face ΔP i P j P k, we connect by straight lines Q ijk with the middle points of the segments P i P j, P j P k,andp k P i, respectively. Then the contribution of K to the control volume K of a vertex P of K is the volume surrounding P by these straight lines, for example, the contribution from one simplex to the control volume K with the interfaces γ 12 and γ 13. Then, the dual finite element space can be constructed for the FVM as follows: X h { [ 3 ) ṽ L Ω ] 2 : ṽ K P K K K h ; ṽ K }. 2.14

5 Mathematical Problems in Engineering 5 P 3 γ 13 Q 13 Q 123 Q D P 1 P 1 Q 12 P 2 γ 12 Figure 1: Control volumes in three-dimensional case. Obviously, the dimensions of X h and X h are the same. Furthermore, there exists an invertible linear mapping Γ h : X h X h such that for N ) v h x v h Pj φj x, x Ω, v h X h, 2.15 j 1 with N ) Γ h v h x v h Pj χj x, j where {φ j } indicates the basis for the finite element space X h,and{χ j } denotes the basis for the finite volume space X h that are the characteristic functions associated with the dual partition K h : 1 ifx K j K h, χ j x otherwise The above idea of connecting the trial and test spaces in the Petrov-Galerkin method through the mapping Γ h was first introduced in 25, 26 in the context of elliptic problems. Furthermore, the mapping Γ h satisfies the following properties 26.

6 6 Mathematical Problems in Engineering Lemma 2.1. Let K K h.ifv h X h and 1 r,then K v h Γ h v h dx, 2.18 Γ h v h c 3 v h, v h Γ h v h L r K c 4h K v h L r K, 2.19 where h K is the diameter of the element K. Multiplying 2.1a by Γ h v h X h and integrating over the dual elements K K h, 2.1b by q h M h and over the primal elements K K h, and applying Green s formula, we define the following bilinear forms for the FVM: N ) A u h, Γ h v h v h Pj j 1 K j u h n ds, u h,v h X h, D ) N ) Γ h v h,p h v h Pj p h n ds, p h M h, K j j b u h,v h, Γ h w h u h v h, Γ h w h 1 2 div u hv h, Γ h w h, ) N ) f, Γh v h v h Pj fdx, v h X h, K j j 1 where n is the unit normal outward to K j and these terms are well posed. As noted above, this paper forces on a class of unstable velocity-pressure pairs consisting of the lowest equal-order finite elements X h {v X : v K R 1 K 3, K τ h }, M h { q M : q K R i K, i, 1, K τ h }, 2.21 where R i K, i, 1 represent piecewise constant range and continuous range on set K, R i, i,1 are spaces of polynomials, the maximum degree of which is bounded uniformly with respect to K τ h and h. The corresponding stabilized FEM is formulated as follows 3 : a u h,v h d v h, p h ) d uh,q h ) G ph,q h ) b uh, u h,v h f, v h ) Xh,M h Also, the corresponding stabilized FVM is defined for the solution u h,p h X h,m h as follows: C h uh,p h ) ; )) b uh,u h, Γ h v h f, Γ h v h ) v h,q h ) Xh,M h, 2.23

7 Mathematical Problems in Engineering 7 where C h uh,p h ) ; )) A uh, Γ h v h D Γ h v h,p h ) d uh,q h ) G ph,q h ) Obviously, the bilinear form G, can be defined by the following symmetry form: 1 G p, q ) p Π h p, q Π h q ) Note that L 2 Ω R if i 1, Π h L 2 Ω R 1 if i Here, the operator Π h satisfies the following properties: 1, 4 ) ) p, qh Πh p, q h p M, q h R i, 2.27 Πh p c 5 p p M, 2.28 p Πh p L p c 6 h ph 1,p p H 1 Ω M In particular, the L 2 -projection operator Π h can be extended to the vector case. This section concentrates on the study of a relationship between the FEM and FVM for the Stokes equations. Lemma 2.2. It holds that A u h, Γ h v h a u h,v h u h,v h X h, 2.3 with the following properties: A u h, Γ h v h A v h, Γ h u h, A u h, Γ h v h c 7 u h v h, 2.31 A v h, Γ h v h ν v h 2. Moreover, the bilinear form D, satisfies 14 D Γ h v h,q h ) d ) v h,q h ) Xh,M h Based on detailed results on existence, uniqueness, and regularity of the solution for the FVM 2.23, the following result establishes its continuity and weak coercivity.

8 8 Mathematical Problems in Engineering Theorem 2.3. It holds that [14] C h uh,p h )) c u h 1 p h ) vh 1 q h ) u h,p h ) Xh,M h Moreover, ) )) Ch uh,p h, sup v h,q h X h,m h v h 1 β u h qh 1 ) ph u h,p h ) Xh,M h, 2.34 where β is independent of h. 3. Stability In this section, we analyze the results of FVM for the three-dimensional stationary Navier- Stokes equations. Firstly, we are now in a position to show the well-posedness of system 2.23 h 4c 2c 3 c 4 γh 1/2 f ν Theorem 3.1 stability. For each h> such that <h 1 2, 3.2 system 2.23 admits a solution u h,p h X h,m h. Moreover, if the viscosity ν>, the body force f Y, and the mesh size h> satisfy <h 1 4, 1 2c c 5 γν 2 f >, 3.3 then the solution u h,p h X h,m h is unique. Furthermore, it satisfies u h 2c 3γ f ν, p h β 1 c 3 γ f 4c3 c 8 ν 2 γ f 1 ). 3.4 Proof. For fixed f Y, we introduce the set B M { vh,q h ) Xh,M h : u h 2c 3γ ν f, ph β 1 c 3 γ f 4c3 c 8 ν 2 γ f 1 )}. 3.5

9 Mathematical Problems in Engineering 9 Then we define the mapping T h : X h,m h X h,m h by 19 A P h w h, Γ h v h D ) ) ) Γ h v h,p h d Th w h,q h G ph,q h b wh,t h w h, Γ h v h ) f, Γ h v h v h,q h ) Xh,M h, 3.6 where T h w h,p h : T 1 w h,t 2 p h u h,p h. We will prove that T h maps B M into B M. First, taking v h,q h u h,p h X h,m h in 3.6, andusing 2.12 and 2.19, wesee that for for all v h,q h X h,m h ν u h 2 A u h, Γ h u h G p h,p h ) since f Γ h u h 2c 2 c 4 h 1/2 w h 2 u h 2 c 3 γ f u h 4c 2c 4 c 3 γh 1/2 f u h 2 ν b w h,u h, Γ h u h b w h,u h, Γ h u h u h ) 3 w h L u h 2 u h L w h Γ h u h u h 2c 2 c 4 h 1/2 w h u h Thus, we have ν 1 4c 2c 3 c 4 γh 1/2 ) f c ν 2 3 γ f u h, 3.9 which implies u h 2c 3γ f. ν 3.1 Then, using the definition of b ;, and C h,, 2.19, setting c 8 2 max{2c 2 c 4 h 1/2,c },and the same approach as above gives that ) )) C h uh,p h, v h β ph, qh b w h ; u h, Γ h v h b w h ; u h, Γ h v h v h b w h ; u h,v h c 8 w h u h v h, ) f, Γ h v h f Γ h v h 3.11 c 3 γ f v h,

10 1 Mathematical Problems in Engineering which, together with 3.1,gives ph β 1 c 3 γ f 4c 3 c 8 ν 2 γ ) f Since the mapping T h is well defined, it follows from Brouwer s fixed-point theorem that there exists a solution to system To prove uniqueness, assume that u 1,p 1 and u 2,p 2 are two solutions to Then we see that C h u1 u 2,p 1 p 2 )) b u1 u 2 ; u 1, Γ h v h b u 2 ; u 1 u 2, Γ h v h Letting v h,q h u 1 u 2,p 1 p 2 e, η,weobtain C h e, η e, η )) ν e 2, b e; u 1, Γ h e b u 2 ; e, Γ h e b e; u 1, Γ h e e b e; u 1,e b u 2 ; e, Γ h e e 2c 2 c 4 h 1/2 u 1 u 2 e 2 c u 1 e 2 2 νh c c 5 γν 1 ) f e 2, 3.14 which together with 3.3 and 3.13,gives ν 1 2c c 5 γν 2 ) f e 2, 3.15 which shows that e by 3.15 ; thatis,u 1 u 2. Next, applying 3.3 to 3.13 and 2.34 yields that p 1 p 2. Therefore, it follows that 2.23 has a unique solution. 4. Optimal Error Estimates Theorem 4.1 optimal error and superconvergent results. Assume that h> satisfies 3.2 and f Y and ν> satisfy 3.2.Let u, p X, M and u h,p h X h,m h be the solution of 2.1 and 2.23, respectively. Then it holds u u h 1 p ph κh u 2 p 1 f ). 4.1 Also, if f H 1 Ω 3, there holds for the solution u h, p h of 2.22 that u h u h 1 ph p h κh 3/2 u 2 p 1 f 1 ). 4.2 Proof. Subtracting 2.1 from 2.23 gives that C h e, η ) ; )) b e, uh, Γ h v h b u h,e,v h b u h,u h,v h Γ h v h, 4.3

11 Mathematical Problems in Engineering 11 with e, η u h u h, p h p h.by v h,q h e, η, it follows that ν e 2 G η, η ) b e, u h,e b u h,u h, Γ h e e f, v h Γ h v h ). 4.4 Using Theorem 3.1, 2.12, 2.23, and 2.25 gives b e, u h,e c u e 2 2c 2 c 4 h 1/2 u h e 2 4c 2c 4 c 5 γh 1/2 f ν 2 e 2 νh e 2, b u h,u h, Γ h e e u h Π h u h u h 1 2 div u h u h Π h u h,e Γ h e) 4.5 u h u h Π h u h e Γ h e ch 3/2 u h e 2. Similarly, by Lemma 2.1 and 2.25, we have f, e Γ h e ) f Πh f, e Γ h e ) Ch i f i e Γ h e 4.6 ch 2 i 1 f 2 i ν 4 e 2 1, i, 1. Combining the above inequalities with 4.3 gives e 1 c h 3/2 h i 1) fi, i, In the same argument, it follows from 2.34 that η c h 3/2 h i 1) fi, i, Noting that 3 u u h 1 p ph ch u 2 p 1 f , and using a triangle inequality completes the proof of Theorem 4.1.

12 12 Mathematical Problems in Engineering As noted above, it is still difficult to achieve an optimal error estimate for the velocity in the L 2 -norm for the three-dimensional stationary Navier-Stokes equations. Here, the following dual problem is proposed and analyzed: a v, Φ d v, Ψ d Φ,q ) b u; v, Φ b v; u, Φ u u h,v. 4.1 Because of convexity of the domain Ω, this problem has a unique solution that satisfies the regularity property 18 Φ 2 Ψ 1 C u u h Below set Φ h, Ψ h I h Φ,J h Ψ X h,m h, which satisfies, by 3.2, Φ Φ h h Φ Φ h 1 Ψ Ψ h Ch 2 Φ 2 Ψ Theorem 4.2 optimal L 2 -error for the velocity. Let u, p be the solution of 2.1a 2.1c and let u h,p h be the solution of 4.3. Then, under the assumptions of Theorem 4.1, it holds u u h Ch 2 u 2 p 1 f 1 ) Proof. Multiplying 2.1a and 2.1b by Γ h Φ h X h and Ψ h M h and integrating over the dual elements K and the primary elements K, respectively, and adding the resulting equations to 2.23 with v h,q h Φ h, Ψ h,weseethat A e, Γ h Φ h D Γ h Φ h,η ) d e, Ψ h G η, Ψ h ) b e; u, Γ h Φ h b u; e, Γ h Φ h b e; e, Γ h Φ h G p, Ψ h 4.14 where e, η u u h,p p h. Subtracting 4.14 from 4.1 with v, q e, η to obtain e 2 a e, Φ Φ h d e, Ψ Ψ h d Φ Φ h,η ) G ) ) η, Ψ h G p, Ψh a e, Φh A e, Γ h Φ h d Φ h,η ) D Γ h Φ h,η ) b u; e, Φ Γ h Φ h b e; u, Φ Γ h Φ h b e; e, Γ h Φ h a e, Φ Φ h d e, Ψ Ψ h d Φ Φ h,η ) G ) ) η, Ψ h G p, Ψh b e; e, Γh Φ h b u; e, Φ Γ h Φ h b e; u, Φ Γ h Φ h f u u, Φ h Γ h Φ h ). 4.15

13 Mathematical Problems in Engineering 13 Obviously, we deduce from Theorem 3.1, , 4.11, the inverse inequality 2.12, and the Cauchy inequality that a e, Φ Φh d e, Ψ Ψ h d Φ Φ h,η ) c e 1 η ) Φ Φh 1 Ψ Ψ h ch 2 u 2 p 1 ) Φ 2 Ψ 1 ch 2 u 2 ) p1 e, ) ) G η, Ψh G p, Ψh ch p Πp ) η Ψ 1 ch 2 u 2 ) p1 e, ) [ f u u, Φ h Γ h Φ h f Πh f ] u u Π h u u, Φ h Γ h Φ h ch 2 ) f1 u u Φh 1 ch 2 f ) 1 u 1/2 u 3/2 2 u 2 1,4 e. b u; e, Φ Γ h Φ h b e; u, Φ Γ h Φ h c u 2 e 1 Φ h Γ h Φ h Φ Φ h ch 2 u 2 p 1 ) Φ 1 ch 2 u 2 p 1 ) e, b e; e, Γ h Φ h b e; e, Γ h Φ h Φ h b e; e, Φ h ) c e,4 e 1 Γ h Φ h Φ h,4 e 2 1 Φ h 1 ch e 1/4 e 7/4 1 Φ h,4 c e 2 1 Φ h 1 ch 2 u 2 p 1 ) e Combining all these inequalities with 4.15 yields In this paper, we have obtained optimal and convergent results of the stabilized mixed finite volume method for the stationary Navier-Stokes equations approximated by the low order finite elements. Furthermore, we could apply the same technique presented to develop and obtain the corresponding results of other stabilized mixed finite volume methods in two or three dimensions. Acknowledgments This research was supported in part by the NSF of China no and , Program for New Century Excellent Talents in University, Natural Science New Star of Science and Technologies Research Plan in Shaanxi Province of China no. 211kjxx12, Research Program of Education Department of Shaanxi Province no. 11JK49, the project sponsored by SRF for ROCS, SEM, and Key project of Baoji University of Arts and Science no. ZK11157.

14 14 Mathematical Problems in Engineering References 1 P. Bochev, C. R. Dohrmann, and M. D. Gunzburger, Stabilization of low-order mixed finite elements for the Stokes equations, SIAM Journal on Numerical Analysis, vol. 44, no. 1, pp , P. Bochev, C. R. Dohrmann, and M. D. Gunzburger, A computational study of stabilized, low-order C finite element approximations of Darcy equations, Computational Mechanics, vol.38,no.4,pp , Y. He and J. Li, A stabilized finite element method based on local polynomial pressure projection for the stationary Navier-Stokes equations, Applied Numerical Mathematics, vol. 58, no. 1, pp , J. Li and Y. He, A stabilized finite element method based on two local Gauss integrations for the Stokes equations, Computational and Applied Mathematics, vol. 214, no. 1, pp , J. Li, Y. He, and Z. Chen, A new stabilized finite element method for the transient Navier-Stokes equations, Computer Methods in Applied Mechanics and Engineering, vol. 197, pp , S. H. Chou and Q. Li, Error estimates in L 2,H 1 and L and in co-volume methods for elliptic and parabolic problems: a unified approach, Mathematics of Computation, vol. 69, no. 229, pp , 2. 7 S. H. Chou and P. S. Vassilevski, A general mixed covolume framework for constructing conservative schemes for elliptic problems, Mathematics of Computation, vol. 68, no. 227, pp , Z. Chen, R. Li, and A. Zhou, A note on the optimal L 2 -estimate of the finite volume element method, Advances in Computational Mathematics, vol. 16, no. 4, pp , R. E. Ewing, T. Lin, and Y. Lin, On the accuracy of the finite volume element method based on piecewise linear polynomials, SIAM Journal on Numerical Analysis, vol. 39, no. 6, pp , R. E. Ewing, R. Lazarov, Y. Lin et al., Finite volume element approximations of nonlocal reactive flows in porous media, Numerical Methods for Partial Differential Equations, vol. 16, no. 3, pp , X. Ye, On the relationship between finite volume and finite element methods applied to the Stokes equations, Numerical Methods for Partial Differential Equations, vol. 17, no. 5, pp , S. H. Chou and D. Y. Kwak, Analysis and convergence of a MAC-like scheme for the generalized Stokes problem, Numerical Methods for Partial Differential Equations, vol. 13, no. 2, pp , S. H. Chou and D. Y. Kwak, A covolume method based on rotated bilinears for the generalized Stokes problem, SIAM Journal on Numerical Analysis, vol. 35, no. 2, pp , J. Li and Z. Chen, A new stabilized finite volume method for the stationary Stokes equations, Advances in Computational Mathematics, vol. 3, no. 2, pp , J. Li, L. Shen, and Z. Chen, Convergence and stability of a stabilized finite volume method for the stationary Navier-Stokes equations, BIT Numerical Mathematics, vol. 5, no. 4, pp , G. He and Y. He, The finite volume method based on stabilized finite element for the stationary Navier-Stokes problem, Computational and Applied Mathematics, vol. 25, no. 1, pp , R. A. Adams, Sobolev Spaces, Academic Press, New York, NY, USA, R. Temam, Navier-Stokes Equations, North-Holland, Amsterdam, The Netherlands, Y. He, A. Wang, and L. Mei, Stabilized finite-element method for the stationary Navier-Stokes equations, Engineering Mathematics, vol. 51, no. 4, pp , V. Girault and P.-A. Raviart, Finite Element Methods for Navier-Stokes Equations: Theory and Algorithms, Springer, Berlin, Germany, P. G. Ciarlet, The Finite Element Method for Elliptic Problems, North-Holland, Amsterdam, The Netherlands, Z. Chen, Finite Element Methods and Their Applications, Springer, Heidelberg, Germany, C. Carstensen, R. Lazarov, and S. Tomov, Explicit and averaging a posteriori error estimates for adaptive finite volume methods, SIAM Journal on Numerical Analysis, vol. 42, no. 6, pp , J. Xu and Q. Zou, Analysis of linear and quadratic simplicial finite volume methods for elliptic equations, Numerische Mathematik, vol. 111, no. 3, pp , R. Li and P. Zhu, Generalized difference methods for second order elliptic partial differential equations I -triangle grids, Numerical Mathematics A Chinese Universities 2, pp , R. Ewing, R. Lazarov, and Y. Lin, Finite volume element approximations of nonlocal reactive flows in porous media, Numerical Methods for Partial Differential Equations, vol. 16, no. 3, pp , 2.

15 Advances in Operations Research Advances in Decision Sciences Mathematical Problems in Engineering Algebra Probability and Statistics The Scientific World Journal International Differential Equations Submit your manuscripts at International Advances in Combinatorics Mathematical Physics Complex Analysis International Mathematics and Mathematical Sciences Stochastic Analysis Abstract and Applied Analysis International Mathematics Discrete Dynamics in Nature and Society Discrete Mathematics Applied Mathematics Function Spaces Optimization