摘要
间断有限元方法能够应用于流体力学计算领域,并受到人们的广泛关注,主要在于它具有以下优点:由于近似解的间断性假设,对网格正则性要求不高,不需要考虑像一般有限元方法中连续性的限制条件就可以对网格进行加密或减疏处理,而且不同的单元可以采用不同形式、不同次数的逼近多项式;它是局部质量守恒的;它几乎没有数值扩散,能够灵活处理数值解含有大梯度大变形和间断的问题,可以准确地捕捉到激波(或大梯度变化区域)的位置,有效地避免数值振荡;间断有限元方法比协调有限元更容易实现hp自适应;而且,在恰当的网格上,改变逼近空间多项式的次数p可以获得按指数收敛的数值解。
虽然有很多学者致力于渗流驱动问题数值方法的研究,但关于用间断有限元方法求解这类问题的较少。对于多孔介质渗流驱动问题(包括可压流和不可压流),我们采用间断有限元方法对空间进行离散求解。
对于不可压的情形,S.Sun和M.F.Wheeler等人作了大量工作,他们建立了Oden-Babu(?)ka-Baumann DG(OBB-DG),non-symmetric interior penaltyGalerkin(NIPG),symmetric interior penalty Galerkin(SIPG)和incomplete interiorpenalty Galerkin(IIPG)四种间断有限元方法半离散形式的能量模误差估计和SIPG的L~2(L~2)模和负模误差估计。但他们并没有用到对偶技巧。我们利用对偶论证,建立了上述四种半离散间断有限元方法的L~2(L~2)模一致后验误差估计式,对全离散格式也进行了分析和研究。
H.Chen和M.R.Cui对于可压缩渗流驱动问题的混合有限元/间断有限元耦合格式做了hp误差分析。但他们的分析仅只限于没有弥散的特殊情况:实际中,分子扩散和弥散的影响都是存在的,这给分析带来了一定的困难。对同时具有分子扩散和弥散效应的可压缩流方程和反应输运方程的耦合问题,我们采用间断有限元方法求解,运用归纳假设导出了先验hp误差估计,运用对偶论证导出了后验hp误差估计;也对于混合有限元/间断有限元耦合格式求解的思想作了论证,导出了先验和后验hp误差估计。在对可压缩渗流驱动问题的先验误差分析过程中,我们采用归纳假设,不象S.Sun所用的引入剪切算子,这样避免了选择该算子中那个大的正常数。在混合有限元方法的超收敛性方面国内外有大量的工作,但对间断有限元方法超收敛性的研究较少,尤其是针对渗流驱动问题的超收敛研究几乎没有。对同时具有分子扩散和弥散效应的不可压缩和可压缩渗流驱动问题进行了混合有限元/间断有限元耦合格式的超收敛分析研究。
Discontinuous Galerkin methods are very attractive for numerical simulations in computational fluid dynamics because of their physical and numerical properties. Firstly,it is flexible that allows for general non-conforming meshes with variable degrees of approximation and of higher accuracy.Secondly,it is locally mass conservative.Thirdly,it has less numerical diffusion and can also handle rough and discontinuous coefficient problems.Fourthly,it is easier for hp-adaptivity than conforming approaches because the information over cell boundaries is almost decoupled. In addition,with appropriate meshing,varying p can yield exponential convergence rates.
Incompressible and compressible miscible displacement problems in porous media are investigated.Several discontinuous Galerkin(DG) finite element methods are used for spatial discretization.For the incompressible case,many scholars such as S.Sun and M.F.Wheeler,have done a lot of work.They established the error estimates in a energy norm for four semi-discrete primal discontinuous Galerkin methods,i.e.,Oden-Babu(?)ka-Baumann DG(OBB-DG),non-symmetric interior penalty Galerkin(NIPG),symmetric interior penalty Galerkin(SIPG), and incomplete interior penalty Galerkin(IIPG) without duality assumption.The estimates in L~2(L~2) norm and in negative norm are also derived for SIPG by them. Based on the duality argument,we work out a unified a posteriori error estimate in L~2(L~2) norm for four semi-discrete and full-discrete primal DG methods above. H.Chen and M.R.Cui present hp error analysis for a combined mixed and discontinuous Galerkin method(MFE/DG) for compressible miscible displacement problems.But only the dispersion-free case is considered.In practical,the effects of molecular diffusion and dispersion are included,which makes our analysis more complicate and difficult.Based on the work of H.Chert and M.R.Cui,we consider the complete compressible case with no restrictions on the diffusion/dispersion tensor. Several DG approximations and combined MFE/DG methods are applied.A priori error estimates and a posteriori one are obtained using the induction hypothesis instead of the cut-off operator and duality technique,respectively.Finally, some superconvergence analysis of combined MFE/DG methods are investigated for incompressible and compressible miscible displacement problems.
引文
[1]S.Adjerid and A.Klauser,Superconvergence of discontinuous finite element solutions for transient convection-diffusion problems,J.Sci.Gomput.,Vol.22-23,No.1-3,(2005),5-24.
[2]T.Arbogast,S.Bryant,C.Dawson,F.Saaf,C.Wang and M.Wheeler,Computational methods for multiphase flow and reactive transport problems arising in subsurface contaminant remediation,J.Comput.Appl.Math.,74(1996),19-32.
[3]D.N.Arnold,An interior penalty finite element with discontinuous element,SIAM J.Numer.Anal.,19(1982),742-760.
[4]D.N.Arnold,F.Brezzi,B.Cockburn and L.D.Marini,Unified analysis of discontinuous Galerkin methods for elliptic problems,SIAM J.Numer.Anal.,39(2002),1749-1779.
[5]Peter Bastian and Beatrice Riviere,Superconvergence and H(div) projection for discontinuous Galerkin methods,Int.J.Numer.Meth.Fluids,42(2003),1043-1057.
[6]F.Brezzi and M.Fortin,Mixed and hybrid finite element methods,Springer-Verlag,New York,(1991).
[7]I.Babuska and M.Suri,The h-p version of the finite element method with quasiuniform meshes,RAIRO Math.Model Numer.Anal.,21(1987),199-238.
[8]I.Babuska and M.Suri,The optimal convergence rate of the p-version of the finite element method,SIAM J.Numer.Anal.,24(1987),759-776.
[9]J.H.Bramble and A.H.Schatz,Estimates for spline projections,RAIRO Anal.Numer.,10(1976),5-37.
[10]J.H.Bramble and A.H.Schatz,Higher order local accuracy by averaging in the finite element method,Math.Comp.,31(1977),94-111.
[11]F.Brezzi and B.Cockburn,L.D.Marini,E.S(u|¨)li,Stabilization mechanisms in discontinuous Galerkin finite element methods,Comput.Methods Appl.Mech.Engrg.,195(2006),3293-3310.
[12]Paul CastiUo,Performance of discontinuous Galerkin methods for elliptic PDE's,SIAM J.Sci.Comput.,Vol.24,No.2,(2002),524-547.
[13]P.Castillo,A superconvergence result for discontinuous Galerkin methods applied to elliptic problems,Comput.Methods Appl.Mech.Engrg.,192(2003),4675-4685.
[14]Y.P.Chen,Application of superconvergence to a model for compressible miscible displacement,Numer.Math.J.of Chinese University,Vol.7,No.1,(1998),25-37.
[15]Y.P.Chen,Superconvergence of the mixed finite element methods for semilinear elliptic problems,CSIAM'96,1996,517-524.(in Chinese)
[16]Y.P.Chen,Superconvergence of a fully discrete procedure for two-phase displacement problems,Nature Science J.of Shandong University,Vol.32,No.1,(1997),1-8.
[17]Y.P.Chen,Superconvergence to the full discrete FEM for compressible miscible displacement with dispersion,CSIAM'98,1998,592-597.(in Chinese)
[18]Y.P.Chen,Superconvergence in the simulation of miscible displacement,J.Sys.Sci & Math.Scis.,Vol.18,No.3,1998,328-334.(in Chinese)
[19]Y.P.Chen,Superconvergence of the full-discrete F.E.M.for compressible miscible displacement:the full tensor case,Numer.Math.J.of Chinese University,Vol.9,No.1,2000,71-82.
[20]H.S.Chen,Superconvergence properties of discontinuous Galerkin methods for two-point boundary value problems,Inter.J.Nurner.Anal.Modeling,Vol.3,No.2,(2006),163-185.
[21]陈焕祯,不可压混溶驱替问题的流线扩散-混合元数值模拟,系统科学与数学,Vol.18,No.1(1998),87-95.
[22]陈传淼,有限元方法及其提高精度的分析,湖南科技出版社,1982.
[23]A.J.Cheng,Optimal error estimate in L~∞(J;L~2(Ω)) norm of FEM for a model for compressible miscible displacement in porous media,Numer.Math.,Vol.4,No.2(1995),222-236.
[24]Z.X.Chen and Hongsen Chen,Pointwise error estimates of the discontinuous Galerkin methods with penalty for second order elliptic problems,SIAM J.Numet.Anal.,Vol.42,No.3,(2004),1146-1166.
[25]陈传淼&黄云清,有限元高精度理论,湖南科技出版社,1995.
[26]H.Chen and Y.P.Chen,A Combined Mixed Finite Element And discontinuous Galerkin method for Compressible miscible displacement problem,Natur.Sci.J.Xiangtan Univ.,Vol.26,No.2(2004),119-126.
[27]P.Castillo,B.Cockburn,I.Perugia and D.Sch(o|¨)tzau,An a priori error analysis of the local discontinuous Galerkin method for elliptic problems,SIAM J.Numer.Anal.,Vol.38,No.5,(2000),1676-1706.
[28]J.L.D.Calle,P.R.B.Devloo and S.M.Gomes,Stabilized discontinuous Galerkin method for hyperbolic equations,Comput.Methods Appl.Mech.Engrg.,194(2005),1861-1874.
[29]B.Cockburn,G.Kanschat,I.Perugia and D.Sch(o|¨)tzau,Superconvergence of the local discontinuous Galerkin method for elliptic problems on Castesian grids,SIAM J.Numer.Anal.,Vol.39,No.1,(2001),264-285.
[30]B.Cockburn,G.E.Karniadakis and C.W.shu,The development of discontinuous Galerkin methods,in:B.Cockburn,G.E.Karniadakis,C.W.shu(Eds.),Discontinuous Galerkin Finite Element Methods,Theory,Computation and Application,Lecture Notes in Computational Science.and Engineering,Vol.11,Springer,New York,2000,3-50.
[31]B.Cockburn and C.W.Shu,The local discontinuous Galerkin method for timedependent convection-diffusion systems,SIAM J.Numer.Anal.,35(1998),2440-2463.
[32]B.Cockburn and C.W.Shu,TVB Runge-Kutta local projection discontinuous Galerkin finite element method for scalar conservations laws Ⅱ:Genaral framework,Math.Comp.,52(1989),411-435.
[33]B.Cockburn,S.Hou and C.W.Shu,TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservations laws Ⅳ:The multidimensional case,Math.Comp.,54(1990),545-581.
[34]M.R.Cui,A combined mixed finite element and discontinuous Galerkin method for compressible miscible displacement problem in porous media,J.Comp.Appl.Math.,Vol.198,No.1,2007,19-34.
[35]V.Dolejsf,M.Feistauer and V.Sobotikova,Analysis of the discontinuous Galerkin method for nonlinear eonvection-diffusion problems,Comput.Methods Appl.Mech.Engrg.,194(2005),2709-2733.
[36]C.Dawson,S.Sun and M.F.Wheeler,Compatible algorithm for coupled flow and transport,Comput.Methods Appl.Mech.Engrg.,193(2004),2565-2580.
[37]J.Douglas Jr.,R.E.Ewing and M.F.wheeler,The approximation of the pressure by a mixed method in the simulation of miscible displacement,RAIRO Numer.Anal.,17(1983),17-33.
[38]J.Douglas Jr.,R.E.Ewing and M.F.wheeler,A time-discretization procedure for a mixed finite element approximation of miscible displacement in porous media,RAIRO Numer.Anal.,17(1983),249-265.
[39]C.Dawson and J.Proft,A priori error estimates for interior penalty versions of the local discontinuous Galerkin method applied to transport equations,Numer.Methods for Part.Diff.Eqns.,Vol.17,No.6(2001),545-564.
[40]J.Douglas Jr.and J.E.Roberts,Numerical methods for a model for compressible miscible displacement in porous media,Math.Comp.,41(1983),441-459.
[41]J.Jr.Douglas and J.E.Roberts,Global estimates for mixed finite element methods for second order elliptic equations,Math.Comp.,44(1985),39-52.
[42]J.Jr.Douglas,superconvergence in the pressure in the simulation of miscible displacement,SIAM J.Numer.Anal.,Vol.22,No.5(1985),962-969.
[43]J.Jr.Douglas and F.A.Milner,Interior and superconvergence estimates for mixed methods for second order elliptic equations,RAIRO Model Math.and Numer.,19(1985),397-428.
[44]J.Jr.Douglas and J.Wang,Superconvergence of mixed finite elements on rectanglar domains,Calaolo,26(1989),121-133.
[45]杜宁,多孔介质中两相完全可压缩可混溶驱动问题的修正迎风差分法,工程数学学报,Vol.22,No.4(2005),590-598.
[46]C.Dawson,S.Sun and M.F.Wheeler,Compatible algorithms for coupled flow and transport,Comput.Methods Appl.Mech.Engrg.,193(2004),2565-2580.
[47]R.E.Ewing,Mathematical modeling and simulation for fluid flow in porous media,Mathematical modeling,Russian Academy of Science(A.A.Samarskii,ed.),2001, 117-127.
[48]A.Ern and J.Proft,A posteriori discontinuous Galerkin error estimates for transient convection-diffusion equations,Applied Mathematics Letters,18(2005),833-841.
[49]R.E.Ewing,R.D.Lazarov and J.Wang,Superconvergence of the velocity along the Gauss lines in mixed finite element methods,SIAM J.Numer.Anal.,28(1991),1015-1029.
[50]R.E.Ewing,J.Shen and J.Wang,Application of superconvergence to problems in the simulation of miscible displacement,Comput.Methods appl.Mech.Engrg.,89(1991),73-84.
[51]S.F.F.Garcia,Improved error estimates for mixed finite-element approximations for nonlinear parabolic equations:the continuous-time case,Numer.Methods Partial Different.Equations,10(1994),129-147.
[52]洪敏纯,多孔介质中可压溶混流的Gslerkin方法,计算数学,10(1998),119-128.
[53]P.Houston,C.Schwab and E.S(u|)li,Discontinuous hp-finite element methods for advection-diffusion-reaction problems,SIAM J.Numer.Anal.,39(2002),2133-2163.
[54]C.Johnson,Discontinuous Galerkin finite element methods for second order hyperbolic problems,Comput.Methods Appl.Mech.Engry.,107(1993),117-129.
[55]O.A.Karakashian,F.Pascal,A posteriori error estimates for a discontinuous Galerkin approximation of second-order elliptic problems,SIAM J.Numer.Anal.,41(2003),2374-2399.
[56]M.Larson and T.Barth,A posteriori error estimation for adaptive discontinuous Galerkin approximations of hyperbolic systems,in Discontinuous Galerkin Methods,B.Cockburn,G:Karniadakis,and C.Shu,eds.,Springer-Verlag in Berlin,2000.
[57]Q.Li and S.Chou,Mixed methods for compressible miscible displacement with the effect of moleeular dispersion,Acta math.appl.Sinica,Vol.11,No.2(1995),123-138.
[58]李潜,可压溶混流的混合元逼近,系统科学与数学,Vol.18,No.3(1998),265-278.
[59]梁栋,两相渗流问题特征混合元数值模拟方法,中国科学(A辑),5(1991),478-485.
[60]梁栋,两数值模拟混溶驱动问题的迎风格式及理论分析,高校应用数学学报,9(1994),118-127.
[61]P.C.Lichtner,C.I.Steefel,and E.H.Oelkers,Reactive transport in porous media,Reviews in Mineralogy,34(1996),438.
[62]L.Machiels,A posteriori finite element bounds for output functionals of discontinuous Galerkin discretizations of parabolic problems,Comput.Methods Appl.Mech.Engrg.,190(2001),3401-3411.
[63]M.T.Nakata,A.Weiser and M.F.Wheeler,Some superconvergence results for mixed element methods for elliptic problems on rectangular domains,The Mathematics of Finite Elements and Applications V(J.R.Whiteman,ed.),Academic Press,1985,367-389.
[64]M.T.Nakata,A.Weiser and M.F.Wheeler,Some superconvergence results for mixed element methods for linear parabolic problems,Lecture notes in Mathematics,1230,Springer,1986,167-174.
[65]J.T.Oden,I.Babuska and C.E.Baumann,A discontinuous hp finite element method for diffusion problems,J.Comput.Phys.,146(1998),491-516.
[66]M.R.Ohm,H.Y.Lee and J.Y.Shin,Error estimates for discontinuous Galerkin method for nonlinear parabolic equations,J.Math.Anal.Appl.,315(2006),132-143.
[67]秦树杰,孙澈,一类非定常线性对流占优扩散问题的间断有限元分析,南开大学学报(自然科学版),Vol.34,No.1,2001,70-77.
[68]G.R.Richter,The discontinuous Galerkin method with diffusion,Math.Comp.,58(1992),631-644.
[69]B.RIVIeRE,Discontinuous Galerkin Methods for solving the miscible displacement problem in Porous Media,Ph.D.thesis,The university of Texas at Austin,2000.
[70]A.Romkes,S.Prudhomme and J.T.Oden,A priori error analysis of a stabilized discontinuous Galerkin method,Comp.Math.Appl.,46(2003),1289-1311.
[71]R.A.Raviart,J.M.Thomas,A mixed finite element method for second order elliptic problems,Mathematics Aspects of the Finite Element Method,Lecture notes in Mathematics,Vol.606,Springer,New York,(1977),292-315.
[72]B.Riviere,and M.F.Wheeler,A posteriori error estimates and mesh adaptation strategy for discontinuous Galerkin methods applied to diffusion problems,Comp.Math.Appl.,46(2003),141-163.
[73]B.Riviere,M.F.Wheeler and V.Girault,A priori error estimates for finite element methods based on discontinuous approximation spaces for elliptic problems,SIAM J.Numer.Anal.,39(2001),902-931.
[74]B.Riviere and M.F.Wheeler,Discontinuous Galerkin methods for coupled flow and transport problems,Comm.Numer.Methods Eng.,18(2002),63-68.
[75]S.Sun,Discontinuous Galerkin Methods for Reactive Transport in Porous Media,Ph.D.thesis,The university of Texas at Austin,2003.
[76]C.Sun,S.Qin,THE full discrete discontinuous finite element analysis for first-order linear hyperbolic equation,J.Comput.Math.,Vol.17,No.1,(1999),97-112.
[77]S.Sun,B.Riviere and M.F.Wheeler,A combined mixed finite element and discontinuous Galerkin method for miscible displacement problem in porous media,In Recent progress in computational and applied PDEs,conference proceedings for the international conference held in Zhangjiajie in July 2001,321-348.
[78]S.Sun,B.Riviere and M.F.Wheeler,A combined mixed finite element and discontinuous Galerkin method for miscible displacement problem in porous media,Recent Progress in Computational and Applied PDEs,Kluwer Academic Publishers,Plenum Press,Dordrecht,New York,(2002),323-351.
[79]L.B.Wahlbin,Superconvergence in Galerkin finite element methods,Lecture notes in Mathematics,Vol.1605,Springer-Verlag,1995.
[80]J.Wang,Superconvergence and extrapolation for mixed finite element methods on rectangular domains,Math.Comput.,Vol.56,No.194(1991),477-503.
[81]S.Sun and M.F.Wheeler,A posteriori error analyses for symmetric discontinuous Galerkin approximations of reactive transport problems,TICAM report 03-19,Institute for Computational Engineering and Sciences,The University of Texas at Austin,Austin,Texas,2003.
[82]S.Sun and M.F.Wheeler,A posteriori error estimation and dynamic adaptivity for symmetric discontinuous Galerkin approximations of reactive transport problems,Computer Methods in Appl.Mechanics and Engineering,Vol.195,No.7-8(2006),632-652.
[83]S.Sun,M.F.Wheeler,Discontinuous Galerldn methods for coupled flow and reactive transport problems,Appl.Numer.Math.,52(2005),273-298.
[84]S.Sun and M.F.Wheeler,Symmetric and non-symmetric discontinuous Galerkin methods for reactive transport problem in porous media,SIAM J.Numer.Anal.,Vol.43,No.1,(2005),195-219.
[85]S.Sun and M.F.Wheeler,L~2(H~1) norm a posteriori error estimation for discontinuous Galerkin approximations of reactive transport problems,Journal of Scientific Computing,22(2005),511-540.
[86]M.F.Wheeler,S.Sun,O.Eslinger,and B.Riviere,Discontinuous Galerkin method for modeling flow and reactive transport in porous media,In W.Wendland,editor,Analysis and Simulation of Multified Problem,Springer Verlag,August 2003,37-58.
[87]朱起定&林群,有限元超收敛理论,湖南科技出版社,1989.