Research Article  Open Access
Aiwen Wang, Jian Li, Zhenhua Di, Xiangjun Tian, Dongxiu Xie, "A New Reduced Stabilized Mixed FiniteElement Method Based on Proper Orthogonal Decomposition for the Transient NavierStokes Equations", Mathematical Problems in Engineering, vol. 2011, Article ID 895386, 19 pages, 2011. https://doi.org/10.1155/2011/895386
A New Reduced Stabilized Mixed FiniteElement Method Based on Proper Orthogonal Decomposition for the Transient NavierStokes Equations
Abstract
A reduced stabilized mixed finiteelement (RSMFE) formulation based on proper orthogonal decomposition (POD) for the transient NavierStokes equations is presented. An ensemble of snapshots is compiled from the transient solutions derived from a stabilized mixed finiteelement (SMFE) method based on two local Gauss integrations for the twodimensional transient NavierStokes equations by using the lowest equalorder pair of finite elements. Then, the optimal orthogonal bases are reconstructed by implementing POD techniques for the ensemble snapshots. Combining POD with the SMFE formulation, a new lowdimensional and highly accurate SMFE method for the transient NavierStokes equations is obtained. The RSMFE formulation could not only greatly reduce its degrees of freedom but also circumvent the constraint of infsup stability condition. Error estimates between the SMFE solutions and the RSMFE solutions are derived. Numerical tests confirm that the errors between the RSMFE solutions and the SMFE solutions are consistent with the the theoretical results. Conclusion can be drawn that RSMFE method is feasible and efficient for solving the transient NavierStokes equations.
1. Introduction
Mixed finiteelement (MFE) methods are one of the most important approaches for solving the nonstationary NavierStokes equations [1–3]. However, some fully discrete MFE formulations for the nonstationary NavierStokes equations involve generally many degrees of freedom. In addition, the importance of ensuring the compatibility of the approximations for the velocity and pressure by satisfying the socalled infsup condition (LBB condition) is widely understood. Thus, an important problem is how to avoid the lack of the LBB stability and simplify the computational load by saving timeconsuming calculations and resource demands in the actual computational process in a way that guarantees a sufficiently accurate numerical solution.
Proper orthogonal decomposition (POD) is an effective method for approximating a large amount of data. The method essentially finds a group of orthogonal bases from the given data to approximately represent them in a least squares optimal sense. In addition, as the POD is optimal in the least squares sense, it has the property that the model decomposition is completely dependent on the given data and does not require assuming any prior knowledge of the process. Combined with a Galerkin projection procedure, POD provides a powerful method for deriving lower dimensional models of dynamical systems from a high or even infinite dimensional space. A dynamic system is generally governed by related structures or the ensemble formed by a large number of different instantaneous solutions, and the POD method can capture the temporal and spatial structures of dynamic system by applying a statistical analysis to the ensemble of data. POD provides an adequate approximation for a large amount of data with a reduced number of degrees of freedom; it alleviates the computational load and provides substantial savings in memory requirements. POD has found widespread application in a variety of fields such as signal analysis and pattern recognition [4, 5], fluid dynamics and coherent structures [6–8], optimal flow control problems [9, 10], and land surface soil moisture data assimilation [11]. In fluid dynamics, Lumley first applied the POD method to capture the large eddy coherent structures in a turbulent boundary layer [12]. This method was further applied to study the relation between the turbulent structure and a chaotic dynamic system [13]. Sirovich introduced the method of snapshots and applied it to reduce the order of POD eigenvalue problem [14]. Kunisch and Volkwein presented Galerkin POD methods for parabolic problems and a general equation in fluid dynamics [15, 16]. More recently, a finite difference scheme (FDS) and a MFE formulation for the nonstationary NavierStokes equation based on POD were derived [17, 18], respectively. Finiteelement formulation based on POD was also applied for parabolic equations and the Burgers equation [19, 20]. In other physical applications, an effective use of POD for a chemical vapor deposition reactor was demonstrated and some reducedorder FDS and MFE for the upper tropical Pacific Ocean model based on POD were presented [21–25]. An optimizing reduced FDS based on POD for the chemical vapor deposit (CVD) equations was also presented in [26]. Except for POD, the empirical orthogonal function (EOF) analysis is another effective method to extract information from large datasets in time and space [27, 28].
In order to avoid the lack of LBB stability, some kinds of stabilized techniques for the lowestorder finite elements appear in [29–44]. Luo et al.[45] has combined the POD method with a stabilize method [40] to deal with the nonstationary NavierStokes equations and obtained good results. But the stabilized mixed methods in [40, 45] are often developed using residuals of the momentum equation. These residual terms must be formulated using meshdependent parameters, whose optimal values are usually unknown. Particularly, for the lowest equalorder pairs of mixed elements such as and , pressure and velocity derivatives in the residual either vanish or are poorly approximated, causing difficulties in the application of consistent stabilization.
In this paper, we mainly consider the twodimensional transient NavierStokes equations by combining a new stabilized finiteelement method [29–31] based on two local Gauss’ integrations with POD method. This new stabilized finite element method has some prominent features: parameterfree, avoiding higherorder derivatives or edgebased data structures, and stabilization being completed locally at the element level. In this manner, we could not only ensure the stabilization of solutions of fully discrete stabilized mixed finiteelement system but also greatly reduce degrees of freedom and save timeconsuming calculations and resource demands in the actual computational process in a way that guarantees a sufficiently accurate numerical solution. we also derive the error estimates between the original SMFE solutions and the RSMFE solutions based on the POD technique. Numerical experiments show the errors between the original SMFE method and the RSMFE solutions are consistent with theoretical results.
The remainder of this paper is organized as follows. In Section 2, an abstract functional setting for the twodimensional NavierStokes equations is given, together with some basic notations. Section 3 is to state the fully discrete stabilized finiteelement method and to generate snapshots from transient solutions computed from the equation system derived by the classical SMFE formulation. In Section 4, the optimal orthogonal bases are reconstructed from the elements of the snapshots with POD method and a reduced SMFE formulation with lowerdimensional number based on POD method for the transient NavierStokes equations is developed. In Section 5, error estimates between the classical SMFE solutions and the RSMFE solutions based on the POD method are derived. In Section 6, a series of numerical experiments are given to illustrate the theoretical results. We conclude with a few remarks in the final section.
2. Functional Setting of the NavierStokes Equations
Let be a bounded domain in , assumed to have a Lipschitz continuous boundary and to satisfy further assumptions below. The transient NavierStokes equations are considered as follows: Here and are the velocity and pressure, is the viscosity, and represents the body forces, the final time, and .
For the mathematical setting of problems (2.1)(2.2), we introduce the following Sobolev spaces: Furthermore, we make a regularity assumption on the Stokes problem as follows.
Assumption H1. For a given and the Stokes problem, satisfying the following regularity result: where is the norm of the Sobolev space or , , as appropriate, and is a positive constant depending only on , which may stand for different value at its different occurrences. Subsequently, the positive constants and (with or without a subscript) will depend only on the data . Because the norm and seminorm are equivalent on , we use the same notation for them. It is well known that for each there hold the following inequalities:
Assumption H2. The initial velocity and the body force are assumed to satisfy
Now, the bilinear forms and , on and , are defined, respectively, by
Also, a generalized bilinear form on is defined by
Moreover, we define the trilinear form
By the above notations and the Hölder inequality, there hold the following estimates:
Also, the Poincare inequality holds:
where , and are positive constants depending only on .
For a given , the variational formulation of problem (2.1)(2.2) reads as follows: find such that
For convenience, we recall the discrete Gronwall Lemma that will be frequently used.
Lemma 2.1 (see [1, 45, 46]). Let , and be three positive sequences, and let be monotone and satisfy then
The following existence and uniqueness result is classical (see [1, 46]).
Theorem 2.2. Assume that (H1) and (H2) hold. Then, for any given , there exists a unique solution satisfying the following regularities: where .
3. Fully Discrete SMFE Method and Generation of Snapshots
In this section, we focus on the stabilized method proposed by [29] for the Stokes equations. Let be a real positive parameter. Finiteelement subspace of is characterized by , a partitioning of into triangles or quadrilaterals , assumed to be regular in the usual sense; that is, for some and with and , where is the diameter of element , is the diameter of the inscribed circle of element , and are the angles of in the case of a quadrilateral partitioning. The mesh parameter is given by . The finiteelement subspaces of this paper are defined by setting Then, the finiteelement pairs are coupled as follows: It is well known that this lowest equalorder finiteelement pair does not satisfy the infsup condition. We define the following local difference between a consistent and an underintegrated mass matrices the stabilized formulation [29–31]: Here, we set where is the basis function of the pressure on the domain such that its value is one at node and zero at other nodes; the symmetric and positive and are pressure mass matrix computed by using the korder and 1order Gauss integrations in each direction; respectively, also, and , are the value of and at the node . is the transpose of the matrix .
Let be the standard projection with the following properties [29–32]: where . Then we can rewrite the bilinear form by
Remark 3.1. The bilinear form in (3.7) is a symmetric, semipositive definite form generated on each local set . The term can alleviate and offset the infsup condition [29]. It differs from the stabilized term in [45]. It does not require a selection of meshdependent stabilization parameters or a calculation of higherorder derivatives. Its another valuable feature is that the action of stabilization operators can be performed locally at the element level with minimal additional cost. With the above notation, we begin by choosing an integer and defining the time step and discrete times , . We obtain the fully discrete scheme as follows: find functions and as solutions of the recursive linear algebraic equations, for all , where and is the approximation of . The solutions and to (3.8)(3.9) are expected to the approximations of and with
Theorem 3.2 (see [29–32]). Let be defined as above, then there exists a positive constant , independent of , such that
By using the same approaches as those in [45], we can prove the following result.
Theorem 3.3. Under the assumptions of Theorems 2.2 and 3.2, if and are sufficiently small and , then problem (3.8)(3.9) has a unique solution such that
Theorem 3.4 (see [32]). Under the assumptions of Theorem 3.3, the error satisfies the following bound:
for all .
If , the time step increment , and the righthand side are given, by solving problem(3.8)(3.9), we can obtain solution ensemble . Then we choose (in general, , e.g., ) instantaneous solutions from the group of solutions for problems (3.8), which are known as snapshots.
Remark 3.5. When one computes actual problems, one may obtain the ensemble of snapshots from physical system trajectories by drawing samples from experiments and interpolation (or data assimilation). For example, when one finds numerical solutions to PDES representing weather forecast, one can use the previous weather prediction results to construct the ensemble of snapshots, then restructure the POD optimal basis for the ensemble of snapshots by the following POD method, next replace finite element space with the subspace spanned by the optimal POD basis. Numerical weather forecast equation is reduced to a fully discrete algebra equation with fewer degrees of freedom. Thus, the forecast of future weather change can be quickly simulated, which is a result of major importance for reallife applications.
4. Reduced SMFE Formulation Based on POD Method
The POD method has received much attention in recent years as a tool to analyze complex physical systems. In this section, we use POD technique to deal with the snapshots in Section 3 and then use the POD basis to develop an RSMFE formulation for the transient NavierStokes equations.
Let , and let , see Section 3). Set where is the ensemble consisting of the snapshots , at least one of which is supposed to nonzero. Let denote an orthogonal basis of with . Then each member of the ensemble is expressed as where , and are the orthogonal basis corresponding to and , respectively.
The method of POD consists in finding the orthogonal basis such that, for every , the mean square error between the elements and corresponding th partial sum of (4.2) is minimized on average subject to where . A solution sequence of (4.3) and (4.4) is known as a POD basis of rank .
We introduce the correlation matrix corresponding to the snapshots by The matrix is positive semidefinite and has rank . The solution of (4.3) and (4.4) can be found in [45].
Proposition 4.1. Let denote the positive eigenvalues of , and let be the associated eigenvectors. Then a POD basis of rank is given by
where denotes the th component of the eigenvector . Furthermore, the following error formula holds:
Let , and let with , and let . Set the Ritzprojection (if is restricted to the Ritzprojection from to , it is written as ) such that and and projection denoted by, respectively,
where and . Owing to (4.8)(4.9) the linear operators and are well defined and bounded:
Lemma 4.2 (see [45]). For every , the projection operators and satisfy, respectively,
where and .
Thus, using , we can obtain the reduced SMFE formulation for problems (3.8) as follows. Find such that
Remark 4.3. Problem (3.8) includes ( is the number of triangles or quadrilaterals vertex in ) freedom degrees, while problem (4.14) includes freedom degrees. For actual science and engineering problems, the number of the vertex in are tens of thousands, even hundreds of millions, but is the number of the largest eigenvalues of snapshots from transient solutions; it is very small. For numerical example in Section 6, , but . Thus, problem (4.14) is a simplified stabilized finiteelement scheme. In addition, the future development of many natural phenomena is affected by previous information, such as biological evolution and weather change. Here, we use the existing data to construct the POD basis, which contains the information on past data. Thus, this method can not only save computational load, but also make better use of the existing information to capture the law of the future development of natural phenomena.
5. Existence and Error Analysis of Solution to the Optimizing RSMFE Formulation
This section is devoted to discussing the existence and error estimates of solutions to problem (4.14). We see from (4.6) that . Using the same approaches as proving Theorem 3.3, we could prove the following existence result for solutions of problem (4.14).
Theorem 5.1. Under the assumptions of Theorems 2.2 and 3.3, Problem (4.14) has a unique solution sequence and satisfies, for ,
In the following theorem, the errors between the solution to Problem (4.14)(4.15) and the solution to Problem (3.8) are derived.
Theorem 5.2. Under the assumptions of Theorem 5.1, if and are sufficiently small, , and , then the errors between the solutions to Problem (4.14), and the solutions to Problem (3.8) have the following error estimates, for :
Proof. Let , . Subtracting (3.8) from (4.14) yields that Taking in (5.4), since , , we deduce Using (2.12)(2.13), the Hölder inequality and the Young inequality, we see that Noting that (for ( and )), owing to (5.6), we obtain that First, we consider the case of . Summing (5.7) from to , , and noting that , using Lemma 4.2, we can derive that By using the discrete Gronwall inequality, we obtain that If and are sufficiently small, , and noting that , we find that Thanks to , we obtain By combining (5.10)(5.11) with Lemma 4.2, we obtain the error estimate result (5.2).
Next, we consider the case of ; we assume that and is the nearest point to . and are expanded into the Taylor series expansion at point . where is the number of time steps from to . If the snapshots are equably taken, then . Summing (5.7) from to , n, and using (5.12), if and are bounded, by discrete Gronwall inequality and Lemma 4.2 and (3.12), we obtain that
If , by (5.14) we obtain that Hence, combining (5.11), (5.13), and (5.15) with Lemma 4.2 yields (5.3).
Theorem 5.3. Under hypotheses of Theorems 3.4 and 5.2, the error estimates between the solution to Problem (2.1)(2.2) and the solutions to Problem (4.14) are as follows:
6. Numerical Examples
In order to illustrate and verify the theoretical results of Theorem 5.3, we present the results obtained in a simple test case. We set is the unit square and viscosity . The velocity and pressure are designed on the same uniform triangulation of . The exact solution is given by and is determined by (2.1).
All the numerical experiments have been performed using the conforming finite element for both velocity and pressure. The implicit (backward) Euler’s scheme is used for the time discretization. For simplicity, the unit square is divided into small squares with side length . In order to make , we take time step increment as .
We obtain 20 values (i.e., snapshots) outputting at time by solving the SMFE formulation. We use 7 optimal POD bases to obtain the solutions of the reduced formulation problem (4.14) as . In Table 1, we present the velocity and pressure relative error estimates, convergence rates, and CPU times using the SMFE method, and, in Table 2, we give the corresponding results obtained using the RSMFE method. In particular, as , the SMFE solutions (, )(c), the exact solutions (,) (b), and the RSMFE solutions () (a) are depicted, respectively, in Figures 1 and 3. Moreover, the difference (, ) (a) between the exact solutions and the RSMFE solutions and the difference (,) (b) between the SMFE solutions and RSMFE solutions are depicted in Figures 2 and 4, respectively. From Tables 1, 2, and Figures 1–4, we can find that the RSMFE solutions has the same accuracy as the reduced SMFE solutions and the exact solutions. As , for the SMFE Problem (3.8)(3.9), there are freedom degrees; the performing time required is 2906 seconds, while the reduced SMFE Problem (4.14) with 7 POD bases only has 7 freedom degrees and the corresponding time is only 14 seconds, that is, the required implementing time to solve the usual SMFE Problem (3.8) is as 207 times as that to do the reduced SMFE problem (4.14) with 7 POD bases, while the errors between their respective solutions do not exceed . As , Figure 5 shows the velocity relative errors between solutions with different number of optimal POD bases and solutions obtained with full bases at and . It is shown that the reduced SMFE problem (4.14) is very effective and feasible. In addition, the results obtained for the numerical examples are consistent with the theoretical ones.


(a)
(b)
(c)
(a)
(b)
(a)
(b)
(c)
(a)
(b)
7. Conclusions
In this paper, we have combined the POD techniques with a SMFE formulation based on two local Gauss’ integrations to derive a reduced SMFE method for the transient NavierStokes equations. The discretization uses a pair of spaces of finite elements over triangles or over quadrilateral elements. This SMFE method differs from that in [45]. It has some prominent features: parameterfree, avoiding higherorder derivatives or edgebased data structures, and stabilization being completed locally at the element level. We have also analyzed the errors between the solutions of their usual SMFE formulation and the solutions of the reduced SMFE based on POD basis and discussed theoretically the relation of the number of snapshots and the number of solutions at all time instances, which have shown that our present method has improved and innovated the existing methods. We have validated the correctness of our theoretical results with numerical examples. Though snapshots and the POD basis of our numerical examples are constructed with the solutions of the usual SMFE formulation, when one computes actual problems, one may structure the snapshots and the POD basis with interpolation or data assimilation by drawing samples from experiments, then solve Problem (4.14), while it is unnecessary to solve Problem (3.8). Thus, the timeconsuming calculations and resource demands in the computational process are greatly saved, and the computational efficiency is vastly improved. Therefore, the method in this paper holds a good prospect of extensive applications.
Future research work in this area aims at addressing some practical engineering problems arising in the fluid dynamics and more complicated PDES, extending the optimizing reduced SMFE formulation, applying it to a realistic atmosphere quality forecast system, and to a set of more complicated nonlinear PDES, for instance, 3D realistic model equations coupling strongly nonlinear properties, nonhomogeneous variable flux, and boundary.
Acknowledgments
This work is supported by the Scientific Research Common Program of the Beijing Municipal Commission of Education (Grant no. KM201110772019), the Knowledge Innovation Program of the Chinese Academy of Sciences (Grant no. KZCX2EWQN207), the National Natural Science Foundation of China (Grant no. 11071193 and 41075076), the Academic Human Resources Development in the Institutions of Higher Learning under the Jurisdiction of Beijing Municipality (Grant no. PHR201006116), the Natural Science New Star of Science and Technologies Research Plan in Shaanxi Province of China (Grant no. 2011kjxx12), Research Program of the Education Department of the Shaanxi Province (Grant no. 11JK0490). The author would like to thank the referees for their valuable comments and suggestions.
References
 V. Girault and P.A. Raviart, Finite Element Methods for NavierStokes Equations: Theory and Algorithms, vol. 5 of Springer Series in Computational Mathematics, Springer, New York, NY, USA, 1986. View at: Zentralblatt MATH
 J. G. Heywood and R. Rannacher, “Finite element approximation of the nonstationary NavierStokes problem: I. Regularity of solutions and secondorder error estimates for spatial discretization,” SIAM Journal on Numerical Analysis, vol. 19, no. 2, pp. 275–311, 1982. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. H. Luo, “The third order estimate of mixed finite element for the NavierStokes problems,” Chinese Quarterly Journal of Mathematics, vol. 10, pp. 9–12, 1995. View at: Google Scholar
 K. Fukunaga, Introduction to Statistical Pattern Recognition, Computer Science and Scientific Computing, Academic Press, New York, NY, USA, 2nd edition, 1990.
 I. T. Jolliffe, Principal Component Analysis, Springer Series in Statistics, Springer, New York, NY, USA, 2nd edition, 2002.
 P. Holmes, J. L. Lumley, and G. Berkooz, Turbulence, Coherent Structures, Dynamical Systems and Symmetry, Cambridge Monographs on Mechanics, Cambridge University Press, Cambridge, UK, 1996. View at: Publisher Site
 P. Moin and R. D. Moser, “Characteristiceddy decomposition of turbulence in a channel,” Journal of Fluid Mechanics, vol. 200, pp. 471–509, 1989. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Rajaee, S. K. F. Karlsson, and L. Sirovich, “Lowdimensional description of freeshearflow coherent structures and their dynamical behaviour,” Journal of Fluid Mechanics, vol. 258, pp. 1–29, 1994. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. D. Joslin, M. D. Gunzburger, R. A. Nicolaides, G. Erlebacher, and M. Y. Hussaini, “A selfcontained automated methodology for optimal flow control validated for transition delay,” AIAA Journal, vol. 35, no. 5, pp. 816–824, 1997. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 H. V. Ly and H. T. Tran, “Proper orthogonal decomposition for flow calculations and optimal control in a horizontal CVD reactor,” Quarterly of Applied Mathematics, vol. 60, no. 4, pp. 631–656, 2002. View at: Google Scholar  Zentralblatt MATH
 X. J. Tian and Z. H. Xie, “An explicit fourdimensional variational data assimilation method based on the proper orthogonal decomposition: theoretics and evaluation,” Science in China D, vol. 52, no. 2, pp. 279–286, 2009. View at: Publisher Site  Google Scholar
 J. Lumley, “Coherent structures in turbulence,” in Transition and Turbulence, R. E. Meyer, Ed., Academic Press, New York, NY, USA, 1981. View at: Google Scholar  Zentralblatt MATH
 N. Aubry, P. Holmes, J. L. Lumley, and E. Stone, “The dynamics of coherent structures in the wall region of a turbulent boundary layer,” Journal of Fluid Mechanics, vol. 192, pp. 115–173, 1988. View at: Publisher Site  Google Scholar
 L. Sirovich, “Turbulence and the dynamics of coherent structures: part I–III,” Quarterly of Applied Mathematics, vol. 45, no. 3, pp. 561–590, 1987. View at: Google Scholar  Zentralblatt MATH
 K. Kunisch and S. Volkwein, “Galerkin proper orthogonal decomposition methods for parabolic problems,” Numerische Mathematik, vol. 90, no. 1, pp. 117–148, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 K. Kunisch and S. Volkwein, “Galerkin proper orthogonal decomposition methods for a general equation in fluid dynamics,” SIAM Journal on Numerical Analysis, vol. 40, no. 2, p. 492–515 (electronic), 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. D. Luo, R. W. Wang, and J. Zhu, “Finite difference scheme based on proper orthogonal decomposition for the nonstationary NavierStokes equations,” Science in China A, vol. 50, no. 8, pp. 1186–1196, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. D. Luo, J. Chen, I. M. Navon, and X. Z. Yang, “Mixed finite element formulation and error estimates based on proper orthogonal decomposition for the nonstationary NavierStokes equations,” SIAM Journal on Numerical Analysis, vol. 47, no. 1, pp. 1–19, 2008. View at: Publisher Site  Google Scholar
 Z. D. Luo, J. Chen, P. Sun, and X. Z. Yang, “Finite element formulation based on proper orthogonal decomposition for parabolic equations,” Science in China A, vol. 52, no. 3, pp. 585–596, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. D. Luo, X. Z. Yang, and Y. J. Zhou, “A reduced finite element formulation based on proper orthogonal decomposition for Burgers equation,” Applied Numerical Mathematics, vol. 59, no. 8, pp. 1933–1946, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Y. H. Cao, J. Zhu, Z. D. Luo, and I. M. Navon, “Reducedorder modeling of the upper tropical Pacific Ocean model using proper orthogonal decomposition,” Computers & Mathematics with Applications, vol. 52, no. 89, pp. 1373–1386, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Y. H. Cao, J. Zhu, I. M. Navon, and Z. D. Luo, “A reducedorder approach to fourdimensional variational data assimilation using proper orthogonal decomposition,” International Journal for Numerical Methods in Fluids, vol. 53, no. 10, pp. 1571–1583, 2007. View at: Publisher Site  Google Scholar
 Z. D. Luo, J. Zhu, R. W. Wang, and I. M. Navon, “Proper orthogonal decomposition approach and error estimation of mixed finite element methods for the tropical Pacific Ocean reduced gravity model,” Computer Methods in Applied Mechanics and Engineering, vol. 196, no. 41–44, pp. 4184–4195, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. D. Luo, J. Chen, J. Zhu, R. W. Wang, and I. M. Navon, “An optimizing reduced order FDS for the tropical Pacific Ocean reduced gravity model,” International Journal for Numerical Methods in Fluids, vol. 55, no. 2, pp. 143–161, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. W. Wang, J. Zhu, Z. D. Luo, and I. M. Navon, “An equationfree reduced order modeling approach to tropic Pacific simulation,” Advances in Geosciences, vol. 12, pp. 1–16, 2007. View at: Publisher Site  Google Scholar
 J. Du, J. Zhu, Z. D. Luo, and I. M. Navon, “An optimizing finite difference scheme based on proper orthogonal decomposition for CVD equations,” International Journal for Numerical Methods in Biomedical Engineering, vol. 27, no. 1, pp. 78–94, 2011. View at: Publisher Site  Google Scholar
 H. L. Yu and H. J. Chu, “Understanding spacetime patterns of groundwater system by empirical orthogonal functions: a case study in the Choshui River alluvial fan, Taiwan,” Journal of Hydrology, vol. 381, no. 34, pp. 239–247, 2010. View at: Publisher Site  Google Scholar
 T. Bernard, O. Krol, H. Linke, and T. Rauschenbach, “Optimal management of regional water supply systems using a reduced finiteelement groundwater model,” Automatisierungstechnik, vol. 57, no. 12, pp. 593–600, 2009. View at: Publisher Site  Google Scholar
 J. Li and Y. N. He, “A stabilized finite element method based on two local Gauss integrations for the Stokes equations,” Journal of Computational and Applied Mathematics, vol. 214, no. 1, pp. 58–65, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. Li, Y. N. He, and Z. X. Chen, “A new stabilized finite element method for the transient NavierStokes equations,” Computer Methods in Applied Mechanics and Engineering, vol. 197, no. 1–4, pp. 22–35, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. Li, Finite element methods for the twodimensional impressible NavierStokes equations, Ph.D. thesis, Xi'an Jiaotong University, 2007. View at: Zentralblatt MATH
 S. Li and Y. R. Hou, “A fully discrete stabilized finite element method for the timedependent NavierStokes equations,” Applied Mathematics and Computation, vol. 215, no. 1, pp. 85–99, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. B. Bochev, C. R. Dohrmann, and M. D. Gunzburger, “stabilization of loworder mixed finite elements,” SIAM Journal on Numerical Analysis, vol. 42, pp. 1189–1208, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. W. Wang and J. Li, “Stabilization of the lowestorder mixed finite elements for the steady NavierStokes equations,” Chinese Journal of Engineering Mathematics, vol. 17, pp. 249–257, 2010. View at: Google Scholar
 Y. N. He, A. W. Wang, and L. Q. Mei, “Stabilized finiteelement method for the stationary NavierStokes equations,” Journal of Engineering Mathematics, vol. 51, no. 4, pp. 367–380, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 T. J. R. Hughes, L. P. Franca, and M. Balestra, “A new finite element formulation for computational fluid dynamics: V. Circumventing the BabuskaBrezzi condition: a stable PetrovGalerkin formulation of the Stokes problem accommodating equalorder interpolations,” Computer Methods in Applied Mechanics and Engineering, vol. 59, no. 1, pp. 85–99, 1986. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 T. J. R. Hughes and L. P. Franca, “A new finite element formulation for computational fluid dynamics: VII. The Stokes problem with various wellposed boundary conditions: symmetric formulations that converge for all velocity/pressure spaces,” Computer Methods in Applied Mechanics and Engineering, vol. 65, no. 1, pp. 85–96, 1987. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. Douglas and J. P. Wang, “An absolutely stabilized finite element method for the Stokes problem,” Mathematics of Computation, vol. 52, no. 186, pp. 495–508, 1989. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 C. Johnson and J. Saranen, “Streamline diffusion methods for the incompressible Euler and NavierStokes equations,” Mathematics of Computation, vol. 47, no. 175, pp. 1–18, 1986. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 B. Kragel, Streamline diffusion POD models in optimization, Ph.D. thesis, Universitat Trier, 2005.
 P. Hansbo and A. Szepessy, “A velocitypressure streamline diffusion finite element method for the incompressible NavierStokes equations,” Computer Methods in Applied Mechanics and Engineering, vol. 84, no. 2, pp. 175–192, 1990. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. A. Behr, L. P. Franca, and T. E. Tezduyar, “Stabilized finite element methods for the velocitypressurestress formulation of incompressible flows,” Computer Methods in Applied Mechanics and Engineering, vol. 104, no. 1, pp. 31–48, 1993. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. Blasco and R. Codina, “Space and time error estimates for a first order, pressure stabilized finite element method for the incompressible NavierStokes equations,” Applied Numerical Mathematics, vol. 38, no. 4, pp. 475–497, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. Codina and J. Blasco, “Analysis of a pressurestabilized finite element approximation of the stationary NavierStokes equations,” Numerische Mathematik, vol. 87, no. 1, pp. 59–81, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. D. Luo, J. Du, Z. Xie, and Y. Guo, “A reduced stabilized mixed finite element formulation based on proper orthogonal decomposition for the nonstationary NavierStokes equations,” International Journal for Numerical Methods in Engineering, vol. 88, no. 1, pp. 31–46, 2011. View at: Publisher Site  Google Scholar
 R. Temam, NavierStokes Equations: Theory and Numerical Analysis, vol. 2 of Studies in Mathematics and its Applications, NorthHolland, Amsterdam, The Netherlands, 3rd edition, 1984.
Copyright
Copyright © 2011 Aiwen Wang 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.