Research Article  Open Access
Yamiao Zhang, Biwu Huang, Jiazhong Zhang, Zexia Zhang, "A Multilevel Finite Element Variational Multiscale Method for Incompressible NavierStokes Equations Based on Two Local Gauss Integrations", Mathematical Problems in Engineering, vol. 2017, Article ID 4917054, 13 pages, 2017. https://doi.org/10.1155/2017/4917054
A Multilevel Finite Element Variational Multiscale Method for Incompressible NavierStokes Equations Based on Two Local Gauss Integrations
Abstract
A multilevel finite element variational multiscale method is proposed and applied to the numerical simulation of incompressible NavierStokes equations. This method combines the finite element variational multiscale method based on two local Gauss integrations with the multilevel discretization using Newton correction on each step. The main idea of the multilevel finite element variational multiscale method is that the equations are first solved on a single coarse grid by finite element variational multiscale method; then finite element variational multiscale approximations are generated on a succession of refined grids by solving a linearized problem. Moreover, the stability analysis and error estimate of the multilevel finite element variational multiscale method are given. Finally, some numerical examples are presented to support the theoretical analysis and to check the efficiency of the proposed method. The results show that the multilevel finite element variational multiscale method is more efficient than the onelevel finite element variational multiscale method, and for an appropriate choice of meshes, the multilevel finite element variational multiscale method is not only timesaving but also highly accurate.
1. Introduction
The finite element method is one of the most general techniques for incompressible NavierStokes equations. However, in numerical simulations of the incompressible NavierStokes equations at higher Reynolds number, the standard Galerkin finite element method is often failed due to the domination of convection term [1–3]. This defect can be overcome to some extent by a variety of stabilizing techniques, including Galerkin least square (GLS) method [4, 5], streamlineupwind PetrovGalerkin (SUPG) method [6], defectcorrection method [7, 8], local projection stabilization (LPS) method [9], discontinuous Galerkin method [10–12], and variational multiscale (VMS) method. Among them, the VMS method is one of some successful stabilized methods that was originally introduced and developed by Hughes et al. [13–16] as a technique in the construction of subgridscale model in large eddy simulation for incompressible and compressible flows. Subsequently, it was also used as a subgrid stabilized technique to improve the accuracy of the solution which was obtained by the Galerkin finite element approximation for the convection dominated problems [17–19]. In the VMS method, the governing equations are rewritten as a coupled system with two types of scales: large scales and small scales. Furthermore, according to Collis [20], the governing equations can also be rewritten as a coupled system with three types of scales: resolved large scales, resolved small scales, and unresolved scales. As solving the above coupled system by finite element method, the traditional Galerkin approximation can be improved by adding an asymptotically consistent artificial diffusion term on the subgrid scales. Under the circumstance, it is necessary to define an appropriate operator, function spaces, and some additional dependent variables. Recently, a finite element VMS method based on two local Gauss integrations was proposed and developed by the authors of [21–25]. This finite element VMS method is simpler than the common finite element VMS method, which does not need additional dependent variables and the projection operator but keeps the same efficiency.
However, when solving nonlinear partial differential equation directly using this finite element VMS method based on two local Gauss integrations, it needs too much computing time, especially for conditions requiring a relatively finegrid discretization. A treatment for this problem is the multilevel discretization. The multilevel method is a wellestablished and efficient method for solving nonlinear partial differential equations. Some classical multilevel methods can be found in [26–32]. The main idea of the multilevel method is to solve the nonlinear equations first on a coarse grid and then solve a linearized problem corresponding to this nonlinear governing equation on a succession of refined grids. Hence, the multilevel method can save much computing time compared with the onelevel method. Nevertheless, there is a shortcoming for these classical multilevel methods in discretization of the convection dominated NavierStokes equations; that is, the nonlinear NavierStokes equations need to be first solved on the coarsest grid by standard Galerkin finite element method. This may lead to numerical oscillations in the coarse grid approximation and transmit remarkable error to the subsequent refined grid numerical approximation and hence cannot yield a correct approximation.
In this paper, by combining the multilevel discretization strategy with the finite element VMS method based on two local Gauss integrations, we propose a multilevel finite element VMS method for the convection dominated NavierStokes equations. In this method, the nonlinear governing equations are first solved on a single coarse grid by finite element VMS method based on two local Gauss integrations, and then finite element VMS approximations are generated on a succession of refined grids using Newton correction on each step by solving a linearized problem corresponding to the nonlinear governing equations. Compared with the classical multilevel methods, the proposed method can solve the convection dominated NavierStokes equations more accurately. While compared with the onelevel finite element VMS method, the proposed method can save plenty of computing time in finegrid discretization.
It is noted that a multilevel VMS method has been reported in [33]. However, there are two essential differences between the proposed method and the one studied in [33]. First, the proposed VMS method in this paper is based on two local Gauss integrations while the one in [33] is the residualbased VMS method. Second, the purpose and the refinement strategy are different. The meshes in the proposed method are global refinement in domains with a Newton correction at each level to increase the efficiency of numerical scheme based on the finite element discretization, while the meshes in [33] are local adaptive refinement in domains via a clear hierarchy, multilevel formulation for further developing the finite cell methodology. Moreover, the twolevel case of the proposed method is also different from those in [34, 35]. Compared with the twolevel projectionbased VMS method in [34], the proposed method in this paper does not need additional dependent variables and the projection operator, so it is simpler. Compared with the twolevel VMS method in [35], the proposed method uses Newton correction to obtain a better approximation than the Stokes correction used in [35]. And what is more, the theoretical analysis and results of [35] are limited to the small data condition, while the condition in this paper is no longer needed, so that the proposed method has a wider effective range than the method in [35].
This paper is organized as follows. In Section 2, the NavierStokes equations, notations, and wellknown results used throughout this paper are given. In Section 3, the finite element VMS method is introduced. Then, the multilevel finite element VMS method is presented and the corresponding stability and convergence error are analyzed in Section 4. In Section 5, three numerical examples are given to verify the theoretical analysis. Finally, conclusions are drawn in Section 6.
2. Preliminaries
In this paper, the following steady NavierStokes equations will be considered:where is a bounded domain in with a Lipschitz continuous boundary , the fluid velocity, the pressure, the prescribed body force, and the kinematic viscosity which represents the inverse of Reynolds number . For the mathematical setting of problem (1), the following Hilbert spaces are needed: The norm and the inner product in are denoted by and , respectively. is the norm of the Sobolev space or , . Using the Poincaré inequality, the norm in space equipped with and can be considered to be equivalent. We define the continuous bilinear forms and , and the trilinear term byrespectively. The trilinear term satisfies the following estimates [36]:where is a positive constant depending only on . Subsequently, and (with a subscript) will denote a positive constant which may stand for different values at its different occurrences, respectively.
With the above notations, the variational formulation of problem (1) reads as follows: find for all satisfying
The nonsingular solution of (1) is defined as follows [37]: is a nonsingular solution of (1) if there is a constant such thatHere . If is a nonsingular solution of the NavierStokes equations, satisfies the estimate [37, 38]
Lemma 1 is also needed; more details of it can be found in [39, 40].
Lemma 1. Suppose that and that the bilinear form satisfies Then the continuous bilinear form defined by satisfies the infsup condition on the product space , which means that there exists a constant such that
For the finite element discretization, let be the regular triangulations of the domain . The mesh defined by is a realpositive parameter tending to 0. We choose the conforming velocitypressure finite element space and consider the TaylorHood elements [41] where is the space of th order polynomials on . The TaylorHood elements on both triangles and rectangles are proved to satisfy the following discrete infsup condition for compatibility of the velocitypressure spaces:
3. Finite Element Variational Multiscale Method
The finite element VMS method proposed by John and Kaya [19, 42] for the steady NavierStokes equations reads find and for all satisfying This system is determined by the choices of and . The stabilization parameter acts only on the small scales which can be chosen as the scale of in order to stabilize the convective term appropriately.
An orthogonal projection operator is defined as which is orthogonal projection with the following properties: where is a positive constant, , and (for details, refer to [18]).
Based on the properties of the orthogonal projection operator , the finite element variational multiscale method proposed by Zheng et al. [21] is as follows: find for all satisfying Here, where denotes an appropriate Gauss integral over which is exact for polynomials of degree .
From [24, 35], problem (16) has the following stability and error estimate.
Theorem 2. Suppose the finite element space satisfies the LBB condition (14). Let be a nonsingular solution of (6). Then the approximate solution given by (16) satisfiesand the following error estimate
The natural energy norm given by will be used in the next section for the theoretical analysis of the proposed multilevel finite element VMS method.
4. Multilevel Finite Element VMS Method
In this section, the multilevel finite element VMS method for (16) will be presented. The stability and error estimate of it will be given. With mesh widths , the finite element space pairs satisfy that And is the solution obtained by the proposed multilevel finite element VMS method. The general steps of the proposed multilevel finite element VMS method for approximating the solution of (16) can be summarized as follows.
Step 1. Solve the following nonlinear system on the coarsest mesh . Find for all satisfying
Step 2. Update on fine mesh with Newton correction at each level.
Give and find such thatfor all with .
Obviously, the solution after first step of the multilevel finite element VMS method satisfies the results in Theorem 2. Before giving the stability analysis and error estimates, the continuous bilinear form is given byHere we assume that is a nonsingular solution of (1) and is close enough to . Thus, by Lemma 1 and (7), the continuous bilinear form satisfies the infsup conditions; namely, there is a constant such that [40]
The stability analysis and error estimate of the proposed multilevel finite element method will be discussed in Theorems 3 and 4.
Theorem 3. Suppose , , and each pair satisfies the LBB condition (14). Let be a nonsingular solution of (6). Suppose and is sufficiently small. Then the approximate solution given by the multilevel finite element method (23) with satisfiesMoreover, if there holds the scaled relation , it has the following error estimate:
Proof. Using (5), (18), (23), (25), and gives Choosing , for sufficiently small , . Then we have Letting , then (26) holds.
Subtracting (23) with from (6) yieldsThe right terms in (30) are bounded according to Taking as interpolation of into , we have Then, together with the relation , denoting the constant as yields (27).
Theorem 4. Suppose , , and each pair satisfies the LBB condition. Let be a nonsingular solution of (6). If is sufficient small, then the approximate solution given by the multilevel finite element method (23) satisfiesMoreover, if there holds the scaled relation , it has the following error estimate:with .
Proof. This theorem will be proved by the induction method. From Theorem 3, it is obvious that it holds for . Assume Theorem 4 is true for ; we need to prove that it is also true for with .
Using (5), (23), (25), and , we obtain Choosing , for sufficiently small , . Then we have Letting , then (33) holds for .
Next, the convergence error of the multilevel finite element method is analyzed. Subtracting (23) with from (6) yieldsThe right terms in (37) are bounded according to Taking as interpolation of into , we have Then, together with the relation , denoting the constant as yields (34) for . Thus, Theorem 4 follows from the induction principle.
5. Numerical Examples
In this section, three numerical examples are given: the first one is an exact solution problem, the second one is the liddriven cavity problem, and the third one is the backward facing step problem. The mentioned methods in this paper are all implemented using Freefem++ [43] and the convergence tolerance is set equal to . In the numerical examples, the triangular mesh with the TaylorHood element pair is used for the finite element discretization. For the sake of simplicity, the onelevel and the proposed multilevel finite element VMS methods are written as level finite element VMS method in the numerical examples. The nonlinear systems in these methods are all linearized by the Oseen iterative method (see [44, 45]).
5.1. Problem with an Exact Solution
This problem is a steady flow of incompressible viscous Newtonian fluid in bounded domain with on the boundaries. It has the following exact solution: where is determined by (1). The TaylorHood elements predict a convergence rate of in the energy norm (i.e., norm for the velocity and norm for the pressure), so the eddy viscosity parameter is chosen as for multilevel finite element VMS method which is the same as that used in [21].
First, is considered. The meshes for the onelevel finite element VMS method are chosen as = [; ; ; ; ], and these meshes are also the finest meshes for the multilevel finite element VMS method. The choices for the coarse meshes of the multilevel finite element VMS method are very important. If the coarse meshes are not fine enough, the optimal convergence orders of velocity and pressure can not be reached. According to Theorems 3 and 4, the coarse meshes and the fine mesh should be chosen as , which ensure optimal accuracy of the approximate solution. Here we consider the multilevel finite element VMS method for , and meshes are chosen as , respectively. The numerical results using the onelevel and the proposed multilevel finite element VMS methods are shown in Figure 1, which indicates that the proposed multilevel finite element VMS method can reach the convergence order of for velocity in norm and pressure in norm and the optimal convergence order of for velocity in the norm. These numerical results are in agreement with the theoretical predictions. Moreover, both the proposed multilevel finite element VMS method and the onelevel finite element VMS method obtain nearly the same approximation results, but the former significantly saves computing time compared with the later.
(a)
(b)
(c)
(d)
Then is considered. The mesh, CPUtime, the error of the velocity, error of the velocity, and error of the pressure for comparison of these level finite element VMS methods are tabulated in Table 1. Here refers to the approximation solution obtained by level finite element VMS method. Table 1 indicates that level finite element VMS method saves much computing time compared with level finite element VMS method at the same fine mesh, and the finer the mesh becomes, the smaller the error achieves. Furthermore, when scrutinizing the error estimates of velocity and the pressure in Table 1, it is easily observed that multilevel finite element VMS method for the velocity and the pressure is better than onelevel finite element VMS method.

5.2. LidDriven Cavity Problem
Liddriven cavity problem often serves as a standard benchmark problem for incompressible NavierStokes equations because the benchmark data is available for comparison. In this problem, the computation is carried out in the region . The velocity is on the top side and on the rest of the boundary.
We compute the approximate solutions for and then compare the results of 1level finite element VMS method (), 2level finite element VMS method (), 3level finite element VMS method (), and 4level finite element VMS method () with that of Ghia et al. [1]. Figures 2–4 plot the horizontal component of the velocity along the vertical line passing through the geometrical center of the cavity and the vertical component of the velocity along the horizontal line passing through the geometrical center of the cavity, respectively. As shown in Figure 2, the results of multilevel finite element VMS method are indistinguishable from 1level finite element VMS method and the benchmark data in [1]. However, in Figures 3 and 4, the 4level finite element VMS method at and the 3level and 4level finite element VMS method at provide unsatisfactory results. This is because the nonlinear solution on the initial coarse mesh is not accurate enough. Actually, with an appropriate choice of initial coarse mesh, the multilevel finite element VMS method can get better results than onelevel finite element VMS method. Table 2 presents the computing time of these level finite element VMS methods and demonstrates that level finite element VMS method saves much computing time compared with level finite element VMS method.

(a)
(b)
(a)
(b)
(a)
(b)
5.3. Backward Facing Step Problem
Backward facing step problem is also a standard numerical example which involves a steady viscous incompressible flow over an isothermal twodimensional backward facing step. It is wellknown that the backward facing step problem possesses a corner singularity; thus it is often adopted by many researchers to verify their numerical methods. The geometry and the boundary conditions are shown in Figure 5 as in [21, 46].
In this example, the Reynolds number is which is based on the maximum inlet velocity and the height of the inlet. We compute an approximate solution on the following number of triangle elements: 9344 (for 1level finite element VMS method), 52989344 (for 2level finite element VMS method), 239252989344 (for 3level finite element VMS method), and 646239252989344 (for 4level finite element VMS method). Figures 6(a)–6(d) show the pressure field and the streamlines obtained by these level finite element VMS methods, respectively, and these results are in agreement with those of [21, 46]. In order to observe these results more clearly, a comparison of the computed velocity, pressure, and vorticity is given across the channel at in Figures 7(a)–7(d), it is obvious that the results of these level finite element VMS methods nearly have no differences with each other in velocity and vorticity, but the pressure computed by 4level finite element VMS method has a clear difference from others. The reason is that the initial grid chosen as 646 in 4level finite element VMS method is too coarse. Table 3 presents the computing time of these level finite element VMS methods; it shows that the level finite element VMS method saves much computing time compared with level finite element VMS method. From this example, it can be concluded that the proposed multilevel finite element VMS method is also efficient for problems with a corner singularity, and it is a highly effective and timesaving method if the coarse grid is chosen appropriately.

(a) 1level finite element VMS method
(b) 2level finite element VMS method
(c) 3level finite element VMS method
(d) 4level finite element VMS method
(a) velocity at
(b) velocity at
(c) Pressure at
(d) Vorticity at
6. Conclusions
In this paper, a multilevel finite element variational multiscale method based on two local Gauss integrations is proposed and applied to the numerical simulation of steady incompressible NavierStokes equations. Theoretically, the stability and error estimate of the multilevel finite element variational multiscale method are analyzed. Numerically, three numerical examples are given. The numerical results demonstrate that, under an appropriate initial coarse mesh, the multilevel finite element variational multiscale method is more efficient than standard onelevel finite element variational multiscale method. In addition, this multilevel finite element variational multiscale method can be extended to solve other more complex fluid dynamical models which will be discussed in further work.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this article.
Acknowledgments
The project is supported by 973 Program (no. 2012CB026002), the NSF of China (no. 51305355), and the National Key Technology R&D Program of China (no. 2013BAF01B02).
References
 U. Ghia, K. N. Ghia, and C. T. Shin, “HighRe solutions for incompressible flow using the NavierStokes equations and a multigrid method,” Journal of Computational Physics, vol. 48, no. 3, pp. 387–411, 1982. View at: Publisher Site  Google Scholar
 E. Erturk, T. C. Corke, and C. Gökçöl, “Numerical solutions of 2D steady incompressible driven cavity flow at high Reynolds numbers,” International Journal for Numerical Methods in Fluids, vol. 48, no. 7, pp. 747–774, 2005. View at: Publisher Site  Google Scholar
 Y. Shang, “A twolevel subgrid stabilized Oseen iterative method for the steady NavierStokes equations,” Journal of Computational Physics, vol. 233, pp. 210–226, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 L. P. Franca and S. L. Frey, “Stabilized finite element methods. II. The incompressible NavierStokes equations,” Computer Methods in Applied Mechanics and Engineering, vol. 99, no. 23, pp. 209–233, 1992. View at: Publisher Site  Google Scholar  MathSciNet
 L. P. Franca and T. J. Hughes, “Convergence analyses of Galerkin leastsquares methods for symmetric advectivediffusive forms of the STOkes and incompressible NavierSTOkes equations,” Computer Methods in Applied Mechanics and Engineering, vol. 105, no. 2, pp. 285–298, 1993. View at: Publisher Site  Google Scholar  MathSciNet
 A. N. Brooks and T. J. Hughes, “Streamline upwind/PetrovGalerkin formulations for convection dominated flows with particular emphasis on the incompressible NavierStokes equations,” Computer Methods in Applied Mechanics and Engineering, vol. 32, no. 13, pp. 199–259, 1982. View at: Publisher Site  Google Scholar  MathSciNet
 W. Layton, H. K. Lee, and J. Peterson, “A defectcorrection method for the incompressible NavierStokes equations,” Applied Mathematics and Computation, vol. 129, no. 1, pp. 1–19, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 Z. Si and Y. He, “A defectcorrection mixed finite element method for stationary conductionconvection problems,” Mathematical Problems in Engineering, vol. 2011, Article ID 370192, 28 pages, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 D. Arndt, H. Dallmann, and G. Lube, “Local projection FEM stabilization for the timedependent incompressible NavierStokes problem,” Numerical Methods for Partial Differential Equations. An International Journal, vol. 31, no. 4, pp. 1224–1250, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 A. Montlaur, S. FernandezMendez, J. Peraire, and A. Huerta, “Discontinuous Galerkin methods for the NavierStokes equations using solenoidal approximations,” International Journal for Numerical Methods in Fluids, vol. 64, no. 5, pp. 549–564, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 N. C. Nguyen, J. Peraire, and B. Cockburn, “An implicit highorder hybridizable discontinuous Galerkin method for the incompressible NavierStokes equations,” Journal of Computational Physics, vol. 230, no. 4, pp. 1147–1170, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 B. Riviere and S. Sardar, “Penaltyfree discontinuous Galerkin methods for incompressible NavierStokes equations,” Mathematical Models and Methods in Applied Sciences, vol. 24, no. 6, pp. 1217–1236, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 T. J. R. Hughes, L. Mazzei, and K. E. Jansen, “Large Eddy Simulation and the variational multiscale method,” Computing and Visualization in Science, vol. 3, no. 12, pp. 47–59, 2000. View at: Publisher Site  Google Scholar
 T. J. R. Hughes, L. Mazzei, A. A. Oberai, and A. A. Wray, “The multiscale formulation of large eddy simulation: Decay of homogeneous isotropic turbulence,” Physics of Fluids, vol. 13, no. 2, pp. 505–512, 2001. View at: Publisher Site  Google Scholar
 T. J. R. Hughes, A. A. Oberai, and L. Mazzei, “Large eddy simulation of turbulent channel flows by the variational multiscale method,” Physics of Fluids, vol. 13, no. 6, pp. 1784–1799, 2001. View at: Publisher Site  Google Scholar
 T. J. R. Hughes, V. M. Calo, and G. Scovazzi, Variational and Multiscale Methods in Turbulence, Mechanics of the 21st Century, 2005, 153–163. View at: Publisher Site
 J.L. Guermond, “Stabilization of Galerkin approximations of transport equations by subgrid modeling,” M2AN. Mathematical Modelling and Numerical Analysis, vol. 33, no. 6, pp. 1293–1316, 1999. View at: Publisher Site  Google Scholar  MathSciNet
 W. Layton, “A connection between subgrid scale eddy viscosity and mixed methods,” Applied Mathematics and Computation, vol. 133, no. 1, pp. 147–157, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 V. John and S. Kaya, “A finite element variational multiscale method for the NavierStokes equations,” SIAM Journal on Scientific Computing, vol. 26, no. 5, pp. 1485–1503, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 S. S. Collis, “Monitoring unresolved scales in multiscale turbulence modeling,” Physics of Fluids, vol. 13, no. 6, pp. 1800–1806, 2001. View at: Publisher Site  Google Scholar
 H. Zheng, Y. Hou, F. Shi, and L. Song, “A finite element variational multiscale method for incompressible flows based on two local Gauss integrations,” Journal of Computational Physics, vol. 228, no. 16, pp. 5961–5977, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 F. Shi, H. Zheng, J. Yu, and Y. Li, “On the convergence of variational multiscale methods based on Newton's iteration for the incompressible flows,” Applied Mathematical Modelling. Simulation and Computation for Engineering and Environmental Systems, vol. 38, no. 23, pp. 5726–5742, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 J. Yu, H. Zheng, and F. Shi, “A finite element variational multiscale method for incompressible flows based on the construction of the projection basis functions,” International Journal for Numerical Methods in Fluids, vol. 70, no. 6, pp. 793–804, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Shang and S. Huang, “A parallel subgrid stabilized finite element method based on twogrid discretization for simulation of 2D/3D Steady incompressible flows,” Journal of Scientific Computing, vol. 60, no. 3, pp. 564–583, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Jiang, L. Mei, H. Wei, W. Tian, and J. Ge, “A finite element variational multiscale method based on two local gauss integrations for stationary conductionconvection problems,” Mathematical Problems in Engineering, vol. 2012, Article ID 747391, 14 pages, 2012. View at: Publisher Site  Google Scholar
 W. Layton and H. W. Lenferink, “A multilevel mesh independence principle for the NavierStokes equations,” SIAM Journal on Numerical Analysis, vol. 33, no. 1, pp. 17–30, 1996. View at: Publisher Site  Google Scholar  MathSciNet
 C. Calgaro, A. Debussche, and J. Laminie, “On a multilevel approach for the twodimensional NavierStokes equations with finite elements,” International Journal for Numerical Methods in Fluids, vol. 27, no. 14, Special Issue, pp. 241–258, 1998. View at: Publisher Site  Google Scholar  MathSciNet
 W. Layton, H. K. Lee, and J. Peterson, “Numerical solution of the stationary NavierStokes equations using a multilevel finite element method,” SIAM Journal on Scientific Computing, vol. 20, no. 1, pp. 1–12, 1998. View at: Publisher Site  Google Scholar  MathSciNet
 Y. He and K.M. Liu, “A multilevel finite element method in spacetime for the NavierStokes problem,” Numerical Methods for Partial Differential Equations. An International Journal, vol. 21, no. 6, pp. 1052–1078, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 J. Li, Y. He, and H. Xu, “A multilevel stabilized finite element method for the stationary NavierStokes equations,” Computer Methods in Applied Mechanics and Engineering, vol. 196, no. 2930, pp. 2852–2862, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 J. Li, Z. Chen, and Y. He, “A stabilized multilevel method for nonsingular finite volume solutions of the stationary 3D NavierStokes equations,” Numerische Mathematik, vol. 122, no. 2, pp. 279–304, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 X. Zhao, J. Li, J. Su, and G. Lei, “Analysis of newton multilevel stabilized finite volume method for the threedimensional stationary NavierStokes equations,” Numerical Methods for Partial Differential Equations, vol. 29, no. 6, pp. 2146–2160, 2013. View at: Publisher Site  Google Scholar
 P. Kopp, Multilevel hpFEM and the Finite Cell Method for the NavierStokes equations using a Variational Multiscale Formulation [M.S. thesis], Technische Universitat Munchen, 2017.
 V. John, S. Kaya, and W. Layton, “A twolevel variational multiscale method for convectiondominated convectiondiffusion equations,” Computer Methods in Applied Mechanics and Engineering, vol. 195, no. 3336, pp. 4594–4603, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Li, L. Mei, Y. Li, and K. Zhao, “A twolevel variational multiscale method for incompressible flows based on two local Gauss integrations,” Numerical Methods for Partial Differential Equations, vol. 29, no. 6, pp. 1986–2003, 2013. View at: Publisher Site  Google Scholar
 R. Temam, NavierStokes Equations, NorthHolland, Amsterdam, The Netherlands, 1984. View at: MathSciNet
 S. Kaya, W. Layton, and B. Rivire, “Subgrid stabilized defect correction methods for the NavierStokes equations,” SIAM Journal on Numerical Analysis, vol. 44, no. 4, pp. 1639–1654, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 Y. He, “Stability and convergence of iterative methods related to viscosities for the 2D/3D Steady NavierStokes equations,” Journal of Mathematical Analysis and Applications, vol. 423, no. 2, pp. 1129–1149, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 F. Brezzi and M. Fortin, Mixed and Hybrid Finite Element Methods, Springer, Berlin, Germany, 1991. View at: Publisher Site  MathSciNet
 W. Layton and L. Tobiska, “A twolevel method with backtracking for the NavierStokes equations,” SIAM Journal on Numerical Analysis, vol. 35, no. 5, pp. 2035–2054, 1998. View at: Publisher Site  Google Scholar  MathSciNet
 C. Taylor and P. Hood, “A numerical solution of the NavierStokes equations using the finite element technique,” Computers and Fluids, vol. 1, no. 1, pp. 73–100, 1973. View at: Publisher Site  Google Scholar  MathSciNet
 V. John and S. Kaya, “Finite element error analysis of a variational multiscale method for the NavierStokes equations,” Advances in Computational Mathematics, vol. 28, no. 1, pp. 43–61, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 “Freefem++,” version 3.38, http://www.freefem.org/. View at: Google Scholar
 Y. He and J. Li, “Convergence of three iterative methods based on the finite element discretization for the stationary NavierStokes equations,” Computer Methods in Applied Mechanics and Engineering, vol. 198, no. 1516, pp. 1351–1359, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 R. An, “Comparisons of STOkes/Oseen/Newton iteration methods for NAVierSTOkes equations with friction boundary conditions,” Applied Mathematical Modelling. Simulation and Computation for Engineering and Environmental Systems, vol. 38, no. 23, pp. 5535–5544, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 A. Masud and R. A. Khurram, “A multiscale finite element method for the incompressible NavierStokes equations,” Computer Methods in Applied Mechanics and Engineering, vol. 195, no. 1316, pp. 1750–1777, 2006. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2017 Yamiao Zhang 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.