MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 747391 10.1155/2012/747391 747391 Research Article A Finite Element Variational Multiscale Method Based on Two Local Gauss Integrations for Stationary Conduction-Convection Problems Jiang Yu 1 Mei Liquan 1,2 Wei Huiming 3 Tian Weijun 1,4 Ge Jiatai 1 Nguyen-Xuan Hung 1 School of Mathematics and Statistics Xi'an Jiaotong University Xi'an 710049 China xjtu.edu.cn 2 Center for Computational Geosciences Xi'an Jiaotong University Xi'an 710049 China xjtu.edu.cn 3 China Nuclear Power Simulation Technology Company Limited Shenzhen 518115 China cgnpc.com.cn 4 College of Mathematics and Information Science Xianyang Normal University Xianyang 712000 China xysfxy.cn 2012 8 11 2012 2012 06 07 2012 23 10 2012 24 10 2012 2012 Copyright © 2012 Yu Jiang 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.

A new finite element variational multiscale (VMS) method based on two local Gauss integrations is proposed and analyzed for the stationary conduction-convection problems. The valuable feature of our method is that the action of stabilization operators can be performed locally at the element level with minimal additional cost. The theory analysis shows that our method is stable and has a good precision. Finally, the numerical test agrees completely with the theoretical expectations and the “ exact solution,” which show that our method is highly efficient for the stationary conduction-convection problems.

1. Introduction

The conduction-convection problems constitute an important system of equations in atmospheric dynamics and dissipative nonlinear system of equations. Many authors have worked on these problems . The governing equations couple viscous incompressible flow and heat transfer process , where the incompressible fluid is the Boussinesq approximation to the nonstationary Navier-Stokes equations. Christon et al.  summarized some relevant results for the fluid dynamics of thermally driven cavity. A multigrid (MG) technique was applied for the conduction-convection problems [11, 12]. Luo et al.  combined proper orthogonal decomposition (POD) with the Petrov-Galerkin least squares mixed finite element (PLSMFE) method for the problems. In , a Newton iterative mixed finite element method for the stationary conduction-convection problems was shown by Si et al. In , Si and He gave a defect-correction mixed finite element method for the stationary conduction-convection problems. In , an analysis of conduction natural convection conjugate heat transfer in the gap between concentric cylinders under solar irradiation was carried out. In , Boland and Layton gave an error analysis for finite element methods for steady natural convection problems. Variational multiscale (VMS) method which defines the large scales in a different way, namely, by a projection into appropriate subspaces, see Guermond , Hughes et al.  and Layton , and other literatures on VMS methods . The new finite element VMS strategy requires edge-based data structure and a subdivision of grids into patches. It does not require a specification of mesh-dependent parameters and edge-based data structure, and it is completely local at the element level. Consequently, the new VMS method under consideration can be integrated in existing codes with very little additional coding effort.

For the conduction-convection problems, we establish such system that Ω be a bounded domain in Rd(d=2  or  3), with Lipschitz-continuous boundary Ω. In this paper, we consider the stationary conduction-convection problem as follows: (1.1)-2ν·D(u)+(u·u)+p=λjT,xΩ,·u=0,xΩ,-ΔT+λu·T=0,xΩ,u=0,T=T0,xΩ, where D(u)=(u+uT)/2 is the velocity deformation tensor, (u,p,T)X×M×W,ΩRd is a bounded convex domain. u=(u1(x),u2(x))T represents the velocity vector, p(x) the pressure, T(x) the temperature, λ>0 the Grashoff number, j=(0,1)T the two-dimensional vector and ν>0 the viscosity.

The study is organized as follows. In the next section, the finite element VMS method is given. In Section 3, we give the stability. The error analysis is given in Section 4. In Section 5, we show some numerical test. The last but not least is the conclusion given in Section 6.

2. Finite Element VMS Method

Here, we introduce some notations (2.1)X=H01(Ω)d,M=L02(Ω)={φL2(Ω);Ωφdx=0},W=H1(Ω). For h>0, finite-dimension subspace (Xh,Mh,Wh)(X,M,W) is introduced which is associated with Ωe, a triangulation of Ω into triangles or quadrilaterals, assumed to be regular in the usual sense. In this study, the finite-element subspaces of personal preference are defined by setting the continuous piecewise (bi)linear velocity and pressure subspace, let τh be the regular triangulations or quadrilaterals of the domain Ω and define the mesh parameter h=maxΩeτh{diam(Ωe)},(2.2)Xh={vX:v|ΩeRl(Ωe)d  Ωeτh},Mh={qM:q|ΩeRl(Ωe)  Ωeτh},Wh={ϕM:ϕ|ΩeRl(Ωe)  Ωeτh}, where W0h=WhH01, l1 is integers. Rl(Ωe)=Pl(Ωe) if Ωe is triangular and Rl(Ωe)=Ql(Ωe) if Ωe is quadrilateral. Here (Xh,Mh) does not satisfy the discrete Ladyzhenskaya-Babuška-Brezzi (LBB) condition (2.3)supvhXhd(vh,ph)vh0βph0,phMh. Now, in order to stabilize the convective term appropriately for the higher Reynolds number and avoid the extra storage, we supply finite element VMS method that the local stabilization form of the difference between a consistent and an underintegrated mass matrices based on two local Gauss integrations at element level as the stabilize term (2.4)G(ph,qh)=ϵd(ak(ph,qh)-a1(ph,qh)). Here, (2.5)ak(ph,qh)=pGTMkqG,a1(ph,qh)=pGTM1qG,pGT=[p1,p2,pN]T,qG=[q1,q2,,qN],Mij=(ϕi,ϕj),ph=i=1Npiϕi,pi=ph(xi),phMh,i=1,2,,N,Mk=(Mijk)N×N,M1=(Mij1)N×N, the stabilization parameter ϵd(ϵd=o(h)) in this scheme acts only on the small scales, ϕi is the basis function of the velocity on the domain Ω such that its value is one at node xi and zero at other nodes, and N is the dimension of Mh. The symmetric and positive matrices Mijk, k2 and Mij1 are the stiffness matrices computed by using k-order and 1-order Gauss integrations at element level, respectively. pi and qi, i=1,2,,N are the values of ph and qh at the node xi. In detail, the stabilized term can be rewritten as (2.6)G(ph,qh)=ϵdΩeτh{Ωe,kphqhdx-Ωe,1phqhdx},ph,qhMh,G(p,q)=(p-hp,q-hq).L2-projection operator h:L2(Ω)R0 with the following properties : (2.7)(p,qh)=(hp,qh)  ,pM,qhR0;hp0cp0,pM;p-hp0chp1,pH1(Ω)M.

Lemma 2.1 (see [<xref ref-type="bibr" rid="B26">26</xref>]).

Let (Xh,Mh) be defined as above, then there exists a positive constant β independent of h, such that (2.8)|B((u,p);(v,q))|c(u1+p0)(v1+q0)(u,p),(v,q)(X,M),β(uh1+vph0)sup(vh,qh)(Xh,Mh)|B((uh,ph);(vh,qh))|v1+q0,(uh,ph)(Xh,Mh),|G(p,q)|C(I-IIh)p0(I-IIh)q0,p,qM.

Using the above notations, the VMS variational formulation of problems (1.1) reads as follows.

Find (A1)(uh,ph,Th)Xh×Mh×Wh such that (2.9)a(uh,vh)-d(ph,vh)+d(qh,uh)+b(uh,uh,vh)+G(ph,qh)=λ(jTh,vh),vhXh,φhMh;a-(Th,ψh)+λb-(uh,Th,ψh)=0,ψhW0h. Given (A2)(uhn-1,Thn-1), find (uhn,phn,Thn)Xh×Mh×Wh such that (2.10)a(uhn,vh)-d(phn,vh)+d(qh,uhn)+b(uhn,uhn-1,vh)+b(uhn-1,uhn,vh)+G(phn,qh)=b(uhn-1,uhn-1,vh)+λ(jThn,vh),vhXh,φhMh;a-(Thn,ψh)+λb-(uhn-1,Thn,ψh)=0,ψhW0h, where a(u,v)=ν(u,v),a-(T,ψ)=(T,ψ)d(q,v)=(q,divv), and (2.11)b(u,v,w)=((u·v),w)+12((divu)v,w)=12((u·)v,w)-12((u·)w,v),b-(u,T,ψ)=((u·T),w)+12((divu)T,ψ)=12((u·)T,ψ)-12((u·)ψ,T).

There exists a constant C which only depends on Ω, such that

u0Cu0,u0,4Cu0,for  all  uH01(Ω)d(orH01(Ω)),

u0,4Cu1,for  all  uH1(Ω)d

u0,421/2u01/2u01/2,for  all  uH01(Ω)d(orH01(Ω))

Assuming ΩCk,α(k0,α>0), then, for T0Ck,α(Ω), there exists an extension T0 in C0k,α(Rd), such that (2.12)T0k,qε,k0,1q,

where ε is an arbitrary positive constant.

b(·,·,·)andb-(·,·,·) have the following properties.

For all uX,v,wX,T·φH01(Ω), there holds that (2.13)b(u,v,w)=-b(u,w,v),b-(u,T,ψ)=-b-(u,ψ,T).

For all uX,vH1(Ω)d,TH1(Ω),for  all  wX(or  φH01(Ω)), there holds that (2.14)|b(u,v,w)|Nu0v0w0,|b-(u,T,φ)|N-u0T0φ0,

where  N=supu,v,w|b(u,v,w)|/u0v0w0,N-=supu,v,w|b-(u,T,φ)|/u0T0φ0.

3. Stability Analysis Lemma 3.1.

The trilinear form b satisfies the following estimate: (3.1)|b(uh,vh,w)|+|b(vh,uh,w)|+|b(w,uh,vh)|C|logh|1/2uh0vh0w0.

Theorem 3.2.

Suppose that (B1)-(B3) are valid and ε is a positive constant number, such that (3.2)64C2Nλε3ν2<1,16C2λ2N-ε3ν<1,T00ε4,T00Cε4. Then (uhm,Thm) defined by (A2) satisfies (3.3)uhm08C2λε3ν,Thm0ε.

Proof.

We prove this theorem by the inductive method. For m=1, (3.3) holds obviously. Assuming that (3.3) holds for m=n-1, we want to prove that it holds for m=n. We estimate Δuhn firstly. Letting vh=uhn,qh=0 in the first equation of (2.10) and using (2.13), we get (3.4)a(uhn,uhn)+b(uhn,uhn-1,uhn)=b(uhn-1,uhn-1,uhn)+λ(jThn,uhn). Setting Thn-1=khn-1+T0 and using (2.14), we have (3.5)νuhn0Nuhn0uhn-10+Nuhn-102+C2λkhn-10+CλT00. Letting Thn=khn+T0,ψ=khn in the second equation of (2.10), we can obtain (3.6)a-(khn,khn)=-λb-(uhn,T0,khn)-a-(T0,khn). Using (2.12), (2.14), and (3.2), we get (3.7)khn-10λN-uhn-10T00+T00λN-ε4uhn-10+T003ε83ε4,(ν-Nuhn-10)uhn0Nuhn-102+C2λεC2λε+64C4N9ν2λ2ε24C2λε3. Using (3.2), we have ν-Nuhn-107ν/8. Then, (3.8)uhn08C2λε3ν. Combining (2.12), (2.14), (3.2), and (3.6), we arrive at (3.9)khn0λN-uhn0T00+T003ε4,(3.10)Thn0khn0+T00ε. Therefore, we finish the proof.

4. Error Analysis

In this section, we establish the H1-bound of the error uhn-u,Thn-T and L2-bounds of the error phn-p. Setting (en,μn,ηn)=(uhn-uh,phn-ph,Thn-Th). Firstly, we give some Lemmas.

Lemma 4.1.

In , If B1-B3 hold, (u,p, T)Hm+1(Ω)×Hm(Ω)×Hm+1(Ω) and (uh,ph,Th)Xh×Mh×Wh are the solution of problem (A1) and (A2), respectively, then there holds that (4.1)(u-uh)0+p-ph0+(T-Th)0Chm(um+1+pm+Tm+1).

Lemma 4.2.

Under the assumptions of Theorem 3.2, (A2) has a unique solution (uh,ph,Th)Xh×Mh×Wh, such that T|Ω=T0 and (4.2)uh08C2λε3ν,Th0ε.

The detail proof we can see [4, 13, 14].

Theorem 4.3.

Under the assumption of Theorem 3.2, there holds (4.3)en0C2λε2n-33ν,ηn0ε2n+1,μn0β-1{νε2+4C2λε3,n=1(ν+2Nε)C2λε2n-33ν+N(C2λε2n-43ν)2+C2λε2n,n2.

Proof.

Subtracting (2.10) from (2.9), we get the following error equations, namely (en,μn,ηn) satisfies (4.4)a(en,vh)-d(μn,vh)+d(qh,en)+b(en,uhn-1,vh)+b(uhn-1,en,vh)+G(μn,qh)=b(en-1,en-1,vh)+λ(jηn,vh),(4.5)a-(ηn,ψh)+λb-(en,Thn,ψh)+λb-(uhn-1,ηn,qh)=0. Here, let ψh=ηn, in (4.5), then we have (4.6)a-(ηn,ηn)+λb-(en,Thn,ηn)=0. By using (2.14), we get (4.7)ηn0λN-εen0. In (4.4), we take vh=enXh,qh=μn, then (4.8)a(en,en)+b(en,uhn-1,en)+b(uhn-1,en,en)+G(μn,μn)=b(en-1,en-1,en)+λ(jηn,en). Using (2.13) and (2.14), we have (4.9)νen0+G(μn,μn)Nen0uhn-10+Nen-102+C2ληn0, then, we obtain (4.10)(ν-Nuhn-10)en0Nen-102+C2ληn0. By using ν-Nuhn-107ν/8. Equations (3.3) and (4.2), we get (4.11)78νen0Nen-102+C2ληn0(Nuhn-10+Nuh0+C2λ2N-ε)en-10(16NC2λε3ν+C2λ2N-ε)en-10=7ν16en-10,en012en-10. From the inductive method, we know, for n=1, subtracting (2.10) from (2.9), we can get (4.12)a(e1,vh)-d(μ1,vh)+d(qh,e1)+b(uh,uh,vh)+G(μ1,qh)=λ(jTn,vh). Letting vh=e1,qh=μ1 in (4.12) and using (2.14), we have (4.13)e10+G(μ1,μ1)ν-1Nuh02+ν-1C2λTh064C4λ2Nε29ν3+C2λεν4C2λε3ν, then (4.14)e104C2λε3ν. By (4.7), we have (4.15)η10λN-εen0λN-4C2λε23νε4. Letting qh=0 in (4.12), (2.14), and (3.9), using Lemma 2.1, we get (4.16)βμ10νe10+Nuh02+CλTh0νε2+4C2λε3ν. Assuming that (4.3) is true for n=k-1, using (4.7) and (4.11), we know that both of them are valid for n=k. Using (4.7) holds for n=k, we let qh=0 in (4.4) and using Lemma 2.1, (4.5), and (3.3), we have (4.17)βμn0(ν+2Nε)en0+Nen-102+C2ληn-10(ν+2Nε)C2λε2n-33ν+N(C2λε2n-43ν)2+C2λε2n.

Theorem 4.4.

Under the assumptions of Theorem 4.3, then there holds that (4.18)limn(uhn-uhn-10+(uhn-uhn-1)0)=0,en0+μn0+ηn0F|logh|1/2(uhn-uhn-1)0uhn-uhn-10+Hε2n+1, where F and H are two positive constants.

Proof.

By using (B1) and triangle inequality, we have (4.19)uhn-uhn-10+(uhn-uhn-1)0(C+1)(en0+en-10). Using Theorem 4.3, letting n, we obtain (4.18). Taking vh=en,qh=μn in (4.4) and using (2.14), we get (4.20)a(en,en)+b(en,uhn-1,en)+G(μn,μn)=-b(uhn-uhn-1,uhn-uhn-1,en)+λ(jηn,en). By (2.14) and Lemma 3.1, we deduce (4.21)(ν-Nuhn-10)en0+G(μn,μn)F|logh|1/2(uhn-uhn-1)0uhn-uhn-10+F2ληn0. Combining (3.3) and (4.7), we obtain (4.22)(ν-8Nε3ν)en0F|logh|1/2(uhn-uhn-1)0uhn-uhn-10+F2λ2N-εen-10. Using (3.2), we get (4.23)en0F|logh|1/2(uhn-uhn-1)0uhn-uhn-10+Hε2n+1. Combining (3.2), (4.7), and (4.17), we get (4.24)ηn0F|logh|1/2(uhn-uhn-1)0uhn-uhn-10+Hε2n+1,μn0F|logh|1/2(uhn-uhn-1)0uhn-uhn-10+Hε2n+1. Here, we complete the proof.

Theorem 4.5.

Under the assumptions of Theorem 4.3, the following inequality: (4.25)(u-uhn)0+p-phn0+(T-Thn)0F1hm(um+1+pm+Tm+1)+F|logh|1/2(uhn-uhn-1)0uhn-uhn-10+Hε2n+1, holds, where F1 and H are the positive constants.

Proof.

By Lemma 4.1, Theorem 4.4, and the triangle inequality, this theorem is obviously true.

5. Numerical Test

This section presents the numerical results that complement the theoretical analysis.

5.1. Convergence Analysis

In our experiment, Ω=[0,1]×[0,1] is the unit square in R2. Let T0=0 on left and lower boundary of the cavity, T/n=0 on upper boundary of the cavity, and T0=4y(1-y) on right boundary of the cavity (see Figure 1). Physics model of the cavity flows: t=0, that is, n=0 initial values on boundary. In general, we cannot know the exact solution of the stationary conduction-convection equations. In order to get the exact solution, we design the procedure as follows. Firstly, solving the stationary conduction-convection equations by using the P2-P1-P2 finite element pair, which holds stability, on the finer mesh, we take the solution as the exact solution. Secondly, the absolute error is obtained by comparing the exact solution and the finite element solutions with VMS methods. Finally, we can easily obtain errors and convergence rates.

From (a) to (c): physics model of the cavity flows, vertical midlines for Re = 2000, h = 1/100, horizontal midlines for Re = 2000, h = 1/100.

5.2. Driven Cavity

In this experiment, Ω=[0,1]×[0,1] is the unit square in R2. Let T0=0 on left and lower boundary of the cavity, T/n=0 on upper boundary of the cavity, and T0=4y(1-y) on right boundary of the cavity (see Figure 1). Physics model of the cavity flows: t=0, that is, n=0 initial values on boundary. Solving the stationary conduction-convection equations by using the P2-P1-P2 finite element pair, which holds stability results, on the finer mesh, we take the solution as the exact solution. From Figures 1 and 2, we know that the solution of finite element VMS using P1-P1-P1 element agree completely with the “exact solution.” In Figure 3, we choose Re=2000, divide the cavity into M×N=100×100, from left to right shows the numerical streamline, the numerical isobar, and the numerical isotherms. In Figure 4, we choose Re=3000, divide the cavity into M×N=100×100, from left to right shows the numerical streamline, the numerical isobar, and the numerical isotherms.

From (a) to (b): vertical midlines for Re = 3000, h = 1/100, horizontal midlines for Re = 3000, h = 1/100.

For Re = 2000, h = 1/100, from (a) to (c): velocity streamlines, the pressure level lines, numerical isotherms.

For Re = 3000, h = 1/100, from (a) to (c): velocity streamlines, the pressure level lines, numerical isotherms.

Remark 5.1.

Our VMS finite element method based on two local Gauss integrations and εd=0.1h is suitable for the Sobolev space. Throughout the paper, our analysis and numerical tests are all carried out for the P1-P1-P1 element (see Tables 1 and 2).

VMS: P1-P1-P1 element.

1 / h u - u h 0 u - u h 1 T - T h 0 T - T h 1 p - p h 0
10 0.000194122 0.00493006 0.00740049 0.277241 0.00506075
20 4.91998 e - 005 0.00252269 0.00208824 0.153561 0.00308312
40 1.21288 e - 005 0.00126459 0.0005746 0.0838877 0.00180539
60 5.35135 e - 006 0.000842093 0.000266991 0.0583954 0.00131444
80 2.98429 e - 006 0.000630808 0.000154418 0.0457979 0.00105007

VMS: P1-P1-P1 element.

1 / h u L 2 rate u H 1 rate T L 2 rate T H 1 rate p L 2 rate
10 / / / / /
20 1.9802 0.9666 1.8253 0.8523 0.7150
40 2.0202 0.9963 1.8617 0.8723 0.7721
60 2.0180 1.0028 1.8903 0.8934 0.7827
80 2.0300 1.0042 1.9033 0.8447 0.7806
6. Conclusion

In this paper, we studied a finite element VMS algorithm based on two local Gauss integrations to solve the stationary conduction-convection problem. From Figures 1 and 2, we see that the solution of VMS using P1-P1-P1 and εd=0.1h agrees completely with the “exact solution,” which shows that our method is highly efficient for the stationary conduction-convection problems. Numerical tests tell us that VMS finite element method based on two local Gauss integrations is very effective.

Acknowledgments

The project is supported by NSF of China (10971164) and the Research Foundation of Xianyang Normal University (06xsyk265).

Chida K. Surface temperature of a flat plate of finite thickness under conjugate laminar forced convection heat transfer condition International Journal of Heat and Mass Transfer 1999 43 4 639 642 2-s2.0-0034141094 García J. A. M. Cabeza J. M. G. Rodríguez A. C. Two-dimensional non-linear inverse heat conduction problem based on the singular value decomposition International Journal of Thermal Sciences 2009 48 6 1081 1093 2-s2.0-62849105354 10.1016/j.ijthermalsci.2008.09.002 Kim D. C. Choi Y. D. Analysis of conduction—natural convection conjugate heat transfer in the gap between concentric cylinders under solar irradiation International Journal of Thermal Sciences 2009 48 6 1247 1258 2-s2.0-62749083648 10.1016/j.ijthermalsci.2008.11.005 Luo Z. D. The Bases and Applications of Mixed Finite Element Methods 2006 Beijing, China Science Press Luo Z. D. Lu X. M. A least-squares Galerkin/Petrov mixed finite element method for stationary conduction-convection problems Mathematica Numerica Sinica 2003 25 2 231 244 2000712 Naveira C. P. Lachi M. Cotta R. M. Padet J. Hybrid formulation and solution for transient conjugated conduction-external convection International Journal of Heat and Mass Transfer 2009 52 1-2 112 123 2-s2.0-57049101468 10.1016/j.ijheatmasstransfer.2008.05.034 ZBL1156.80357 Wang Q. W. Yang M. Tao W. Q. Natural convection in a square enclosure with an internal isolated vertical plate Wärme-und Stoffübertragung 1994 29 3 161 169 2-s2.0-0347512612 10.1007/BF01548600 Yang M. Tao W. Q. Wang Q. W. Lue S. S. On identical problems of natural convection in enclosures and applications of the identity character Journal of Thermal Science 1993 2 2 116 125 2-s2.0-0027605761 10.1007/BF02718267 DiBenedetto E. Friedman A. Conduction-convection problems with change of phase Journal of Differential Equations 1986 62 2 129 185 2-s2.0-38249042743 833415 ZBL0593.35085 10.1016/0022-0396(86)90095-1 Christon M. A. Gresho P. M. Sutton S. B. Computational predictability of time-dependent natural convection flows in enclosures (including a benchmark solution) International Journal for Numerical Methods in Fluids 2002 40 8 953 980 2-s2.0-0037146118 10.1002/fld.395 ZBL1025.76042 Mesquita M. S. de Lemos M. J. S. Optimal multigrid solutions of two-dimensional convection-conduction problems Applied Mathematics and Computation 2004 152 3 725 742 2-s2.0-2342641512 10.1016/S0096-3003(03)00591-5 2062773 ZBL1077.65508 Cheikh N. B. Beya B. B. Lili T. A multigrid method for solving the Navier-Stokes/Boussinesq equations Communications in Numerical Methods in Engineering 2008 24 8 671 681 2-s2.0-49549085467 10.1002/cnm.980 2444788 ZBL1159.76355 Luo Z. Chen J. Navon I. M. Zhu J. An optimizing reduced PLSMFE formulation for non-stationary conduction-convection problems International Journal for Numerical Methods in Fluids 2009 60 4 409 436 2-s2.0-67649961644 10.1002/fld.1900 2519402 ZBL1161.76032 Si Z. Zhang T. Wang K. A Newton iterative mixed finite element method for stationary conduction-convection problems International Journal of Computational Fluid Dynamics 2010 24 3 135 141 2-s2.0-77954444721 10.1080/10618562.2010.495931 2682925 ZBL1206.65249 Si Z. He Y. A defect-correction mixed finite element method for stationary conduction-convection problems Mathematical Problems in Engineering 2011 2011 28 370192 10.1155/2011/370192 2764986 ZBL1204.76021 Boland J. Layton W. Error analysis for finite element methods for steady natural convection problems Numerical Functional Analysis and Optimization 1990 11 5-6 449 483 10.1080/01630569008816383 1079286 ZBL0714.76090 Guermond J. L. Stabilization of Galerkin approximations of transport equations by subgrid modelling Mathematical Modelling and Numerical Analysis 1999 33 6 1293 1316 2-s2.0-0033264647 10.1051/m2an:1999145 1736900 ZBL0946.65112 Hughes T. J. R. Mazzei L. Oberai A. A. Wray A. A. The multiscale formulation of large eddy simulation: decay of homogeneous isotropic turbulence Physics of Fluids 2001 13 2 505 512 2-s2.0-0035065963 10.1063/1.1332391 ZBL1184.76236 Hughes T. Mazzei L. Jansen K. Large eddy simulation and the variational multiscale method Computing and Visualization in Science 2000 3 1-2 47 59 10.1007/s007910050051 ZBL0998.76040 Hughes T. J. R. Multiscale phenomena: green's functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods Computer Methods in Applied Mechanics and Engineering 1995 127 1–4 387 401 2-s2.0-0029405472 ZBL0866.76044 Layton W. A connection between subgrid scale eddy viscosity and mixed methods Applied Mathematics and Computation 2002 133 1 147 157 2-s2.0-0037175406 10.1016/S0096-3003(01)00228-4 ZBL1024.76026 John V. Kaya S. A finite element variational multiscale method for the Navier-Stokes equations SIAM Journal on Scientific Computing 2005 26 5 1485 1503 2-s2.0-25444505626 10.1137/030601533 ZBL1073.76054 Zheng H. Hou Y. Shi F. Song L. A finite element variational multiscale method for incompressible flows based on two local gauss integrations Journal of Computational Physics 2009 228 16 5961 5977 2-s2.0-67649361339 10.1016/j.jcp.2009.05.006 ZBL1168.76028 Masud A. Khurram R. A. A multiscale finite element method for the incompressible Navier-Stokes equations Computer Methods in Applied Mechanics and Engineering 2006 195 13–16 1750 1777 2-s2.0-30944434045 10.1016/j.cma.2005.05.048 ZBL1178.76233 Li J. He Y. Chen Z. A new stabilized finite element method for the transient Navier-Stokes equations Computer Methods in Applied Mechanics and Engineering 2007 197 1–4 22 35 2-s2.0-35348842116 10.1016/j.cma.2007.06.029 ZBL1169.76392 Li J. He Y. A stabilized finite element method based on two local Gauss integrations for the Stokes equations Journal of Computational and Applied Mathematics 2008 214 1 58 65 2-s2.0-38949209539 10.1016/j.cam.2007.02.015 ZBL1132.35436