摘要
在当代科学与工程计算领域,许多重要的实际问题是复杂的多尺度问题,比如复合材料的热(电)性质、多孔介质的流体分析、湍流现象、集成电路设计、化学反应的时间尺度等等。多尺度问题的复杂性和多样性催生了各种自有特色的多尺度计算方法。
本文利用求解局部偏微分方程获得的基函数具备良好的性态,该基函数由方程的微分算子性质所决定,可以自适应地将微观尺度的信息带入到宏观尺度。针对带多个尺度参数ε_k的强振荡变系数多尺度椭圆问题,重点研究多尺度有限元方法(MFEM)和单位分解法(PUM),大大减少了计算资源的消耗,并得到了精确的、收敛的数值结果。研究奇异摄动反应扩散问题的多尺度有限元解法,发现利用多尺度基改善了奇异摄动的边界层误差。我们给出了多尺度方法的数学理论,通过实验证明了该方法用有限的计算资源可以获得好的结果,并对算例结果做了相应的数值分析。
多尺度有限元法通过在宏观尺度上求解多尺度基函数的局部方程,将微观尺度的信息带入基函数中,使得宏观粗网格的解具有良好的计算精度,能够节约大量的计算资源.单位分解法是一种重要的无网格方法,利用不依赖于网格的单位分解函数和灵活的局部逼近空间,可以从全局上很好地逼近未知解。
本文获得了以下主要结果:
1.使用多尺度有限元法求解强振荡变系数问题,在系数可分离变量时限定基函数的扰动边界条件所得到的结果精度更好,而且L~2范数不会出现负阶收敛。
2.当系数非分离变量时,我们使网格边界以有理角或无理角穿越周期系数和非周期系数的场域,多尺度基扰动边界条件当ε_k→0时可以获得各方向上单元的完整信息因而有更好的求解精度;反之当ε_k(?)0且网格以无理角穿过场域时得到的单元信息不完全,这种情形下限定多尺度基线性边界条件才能获得更好的精度。
3.用矩形2阶再生插值的单位分解方案处理强振荡变系数问题,与经典有限元法相比证明了单位分解法的效率更高,且有2阶收敛率。
4.研究奇异摄动反应扩散问题,通过在边界层区域使用多尺度基函数、内部光滑区域使用标准线性基函数来构造多尺度有限元空间,探讨了适当厚度边界层区域的选取问题,使用多尺度有限元法在一致的粗网格获得了与ε大小无关的一致收敛结果。
文章的创新之处在于针对多尺度有限元法的多尺度基函数,精心设计了系统的实验方案研究不同的线性边界条件与扰动边界条件的效果,针对周期函数和非周期函数的场域考虑尺度参数ε_k遍历性的影响,给出了具体环境下限定何种基函数边界条件的原则。针对奇异摄动反应扩散问题,在合适的边界厚度区域使用多尺度基可以改善边界层误差,在一致粗网格上得到了ε的一致收敛。
文章的内容安排如下:
第一章简要介绍多尺度问题背景和各种多尺度计算方法,并统一文章中用到的符号。
第二章给出多尺度有限元法的理论基础,指出多尺度基的优点在于能够如实反映微分算子的性质如高阶振荡性,按尺度参数ε_k与网格步长h的大小关系分情形证明了多尺度解的H~1误差估计和L~2误差估计的收敛性定理。
第三章使用多尺度有限元法求解强振荡变系数问题,对多尺度基函数的子问题定义线性边界条件与扰动边界条件。分别考虑系数可分离变量与非分离变量的情况,通过系统地研究数值实验给出了多尺度有限元法基函数最佳边界条件的选取原则。
第四章针对强振荡变系数问题使用单位分解法,介绍该方法的理论基础,给出问题的矩形2阶再生插值的单位分解方案,通过实验验证了单位分解法处理变系数问题有良好的效果。
第五章研究奇异摄动反应扩散问题,在边界区域利用多尺度基可以改善边界层误差获得高精度,多尺度有限元法可以在一致的粗网格上得到2阶收敛的L~2范数误差,1阶收敛的能量范数误差。
In modern scientific and engineering computation,many important practical problems are complicated and troublesome multiscale problems,such as composite materials with thermal/electrical conductivity,flow through the heterogeneous porous media,turbulent transport problem,designing of the integrate circuit,and time scale of the chemical reactions,etc..There come up with many kinds of multiscale computation methods because of the intrinsic complexity of the multiscale problems.
This paper acquires the optimal multiscale base functions by solving the local partial differential equations,these base functions are adapted to the local property of the differential operator and can capture the small scales information onto the large scales.We put emphases on studying the multiscale finite element method(MFEM) and the partition of unity method(PUM) for the rapidly oscillating coefficients multiscale elliptic problem with several scale parametersε_k,the methods can greatly reduce the costing of the computer resources and obtain accurate and convergent numerical results.And we investigate the MFEM for solving the singularly perturbed reaction-diffusion problem,find that the method improves the boundary layer error with the help of the multiscale bases.We give the mathematical theories of our multiscale methods,and through numerical experiments we approve that the methods have the abilities to obtain good results with just costing a few handy computer resources,and achieve the numerical analysis for the experiment results.
The multiscale finite element method is capable of accurately capturing the large scale behaviors of the solution by resolving the local PDEs for the multiscale base functions,which extract the small scale information of the solution,as a consequence the method acquires good accuracy on coarse meshes and saves plenty of computer resources.The partition of unity method is an important kind of meshless methods,it is made up of two parts,that is,the partition of unity function which may be independent of the grid,plus the agile local approximation space, thus can approximate the global unknown more precisely.
In this paper we achieve the following main results.
1.We study the MFEM for solving the elliptic problem with rapidly oscillating coefficients,when the coefficients are variables separated,the results of the MFEM with the oscillatory bases are more precise and don't show negative convergence rate for the L~2 norm.
2.When the coefficients are not variables separated,we make the mesh grids cut across the periodic and nonperiodic coefficients' fields under rational or irrational angle,the multiscale bases with the oscillatory boundary condition are able to capture the entire period information from each direction in the case ofε_k→0, thus may get better accuracy.On the contrary whenε_k(?)0 and the grids go through the fields under irrational angle,the inadequate element information by the oscillatory bases is inferior to those by the linear ones,in this case we use the linear multiscale bases may gain preciser results.
3.We propose the rectangular second order reproducing interpolation partition of unity for the rapidly oscillating coefficients problem,comparing with the FEM we demonstrate the high efficiency of the PUM and obtain the second order convergence rate.
4.We investigate the singular perturbed reaction-diffusion problem,and construct the multiscale finite element space by enriching the multiscale bases in the boundary layers domain plus the standard linear bases in the inner smooth domain, and discuss the issue of the optimal width of boundary layers domain,we can acquire independentlyε-uniform convergence on the uniform coarse meshes by the MFEM.
The innovations in this paper are represented that we design the systematic numerical experiments elaborately on studying the linear and the oscillatory boundary conditions for the multiscale bases in the MFEM,consider the ergodicity of the scale parametersε_k in the periodic and non-periodic cocfficients' fields,and present the guideline for choosing of the optimal bases boundary condition under different cases.For the singular perturbed reaction-diffusion problem,we use the multiscale bases in the appropriate boundary domains can improve the boundary layer error and on the uniform coarse meshes we obtainε-uniform convergence.
The paper is organized as follows.
In chapter one we introduce the multiscale problem backgrounds and existing kinds of multiscale computation methods,and unify the denotations used in this paper.
In chapter two the theoretical foundation of the multiscale finite element method is provided,we point out that the multiscale bases can reflect the property of the differential operator,for instance highly oscillating property,and prove the convergence theorems of H~1 error estimate and L~2 error estimate for the multiscale solution according to different cases between the scale parametersε_k and the coarse mesh size h.
In chapter three we study the MFEM for solving the elliptic problem with rapidly oscillating coefficients,define the linear boundary condition and the oscillatory boundary condition of the multiscale bases for the local PDEs.Considering the separated variables and the non-separated variables coefficients respectively, we present the guideline for defining the multiscale bases optimal boundary of the MFEM through systematic numerical studies.
In chapter four we apply the PUM for solving the rapidly oscillating coefficients elliptic problem and provide the theories of this method,we propose the rectangular second order reproducing interpolation partition of unity,through numerical experiment we demonstrate the efficiency of the PUM.
In the last chapter five the two-dimensional singular perturbed reactiondiffusion problem is studied,in the appropriate boundary domains the multiscale bases can improve the boundary layer error and obtain high accuracy,thus on the uniform coarse meshes the MFEM acquires second order convergence rate for the L~2 norm error and first order convergence rate for the energy norm error,respectively.
引文
[1]J.Aarnes and T.Y.Hou,Multiscale domain decomposition methods for elliptic problems with high aspect ratios,Acta Math.Appl.Sinica,2002,18(1):63-76.
[2]R.A.Adams,Sobolev spaces,Academic Press,New York,1975.
[3]G.Allalre,Homogenization and two-scale convergence,SIAM J.Math.Anal.,1992,23(6):1482-1518.
[4]G.Allaire and R.Brizzi,A multiscale finite element method for numerical homogenization,Multiscale Model.Simul.,2005,4(3):790-812.
[5]R.Araya and F.Valentin,A multiscale a posteriori error estimate,Comput.Methods Appl.Mech.Engrg.,2005,194(18-20):2077-2094.
[6]I.Babuska and E.Osborn,Generalized finite element methods:Their performance and their relation to mixed methods,SIAM J.Numer.Anal.,1983,20:510-536.
[7]I.Babuska,G.Caloz and E.Osborn,Special finite element methods for a class of second order elliptic problems with rough coefficients,SIAM J.Numer.Anal.,1994,31(4):945-981.
[8]I.Babuska,B.Andersson,B.Guo,et.al,Finite element method for solving problems with singular solutions,J.Comput.Appl.Math.,1996,74(1):51-70.
[9]I.Babuska and J.M.Melenk,The partition of unity method,Int.J.Numer.Meth.Engrg.,1997,40(4):727-758.
[10]I.Babuska,U.Banerjee and J.E.Osborn,Generalized finite element methods:main ideas,results and perspective,Int.J.Comput.Meth.,2004,1(1):67-103.
[11]J.Bear,Use of models in decision making,Transport and Reactive Processes in Aquifers,Rotterdam,1994,3-9.
[12]T.Belytschko,Y.Y.Lu,and L.Gu,Element-free Galerkin methods,Int.J.Numer.Meth.Engrg.,1994,37(2):229-256.
[13]T.Belytschko,Y.Krongauz,D.Organ,et.al,Meshless methods:An overview and recent developments,Comput.Methods Appl.Mech.Engrg.,1996,139(1):3-47.
[14]A.Bensoussan,J.L.Lion and G.Papanicolaou,Asymptotic analysis for periodic structure,North-Holland,Amsterdam,1978.
[15]A.Brandt,Multi-level adaptive solutions to boundary value problems,Math.Comput.,1977,31(138):333-390.
[16]S.C.Brenner and L.R.Scott,The mathematical theory of fnite element methods,Springer-Verlag,New York,1994.
[17]J.R.Chen and J.Z.Cui,A multiscale finite element method for elliptic problems with highly oscillatory coefficients,Appl.Numer.Math.,2004,50(1):1-13.
[18]S.Q.Chen,W.N.E,Y.X.Liu and C.W.Shu,A discontinuous Galerkin implementation of a domain decomposition method for kinetic-hydrodynamic coupling multiscale problems in gas dynamics and device simulations,J.Comput.Phys.,2007,225(2):1314-1330.
[19]Y.Chen,Uniform pointwise convergence for a singularly perturbed problem using arc-length equidistribution,J.Comput.Appl.Math.,2003,159(1):25-34.
[20]Y.Chen,Uniform convergence analysis of finite difference approximations for singular perturbation problems on an adapted grid,Advan.Comput.Math.,2006,24(1-4):197-212.
[21]Z.M.Chen and T.Y.Hou,A mixed multiscale finite clement method for elliptic problems with oscillating coefficients,Math.Comput.,2002,72(242):541-576.
[22]Z.M.Chen and X.Y.Yue,Numerical homogenization of well singularities in the flow transport through heterogeneous porous media,Multiscale Model.Sim.,2003,1(2):260-303.
[23]P.G.Ciarlet,The finite element method for elliptic problems,North-Holland,Amsterdam,1978.
[24]J.Z.Cui and L.Q.Cao,Two-scale asymptotic analysis method for a class of elliptic boundary value problems with small periodic coefficient,Math.Numer.Sci.,1999,21(1):19-28.
[25]M.Dorobantu and B.Engquist,Wavelet-based numerical homogenization,SIAM J.Numer.Anal.,1998,35(2):540-559.
[26]C.A.Duarte and J.T.Oden,An h-p adaptive method using clouds,Comput.Methods Appl.Mech.Engrg.,1996,139(1):237-262.
[27]W.N.E and B.Engquist,The heterogeneous multiscale methods,Comm.Math.Sci.,2002,1(1):87-132.
[28]W.N.E,P.B.Ming and P.W.Zhang,Analysis of the heterogeneous multiscale method for elliptic homogenization problems,J.Amer.Math.Soc.,2004,18(1):121-156.
[29]W.N.E,B.Engquist,X.T.Li,et.al,Heterogeneous multiscale methods:A review,Commun.Comput.Phys.,2007,2(3):367-450.
[30]Y.Efendiev,T.Y.Hou and X.H.Wu,Convergence of a nonconforming multiscale finite element methods,SIAM J.Numer.Anal.,2000,37(3):888-910.
[31]Y.Efendiev and X.H.Wu,Multiscale finite element for problems with highly oscillatory coefficients,Numer.Math.,2002,90(3):459-486.
[32]B.Engquist and E.Luo,Convergence of a multigrid method for elliptic equations with highly oscillatory coefficients,SIAM J.Numer.Anal.,1997,34(6):2254-2273.
[33]J.Fish and Z.Yuan,Multiscale enrichment based on partition of unity,Int.J.Numer.Meth.Engrg.,2005,62(10):1341-1359.
[34]L.P.Franca,A.L.Madureira and F.Valentin,Towards multiscale functions:enriching finite element spaces with local but not bubble-like functions,Comput.Methods Appl.Mech.Engrg.,2005,194(27-29):3006-3021.
[35]X.G.He and L.Ren,Finite volume multiscale finite element method for solving the groundwater flow problems in heterogeneous porous media,Water Resour.Res.,2005,41,W10417 doi:10.1029/2004WR003934.
[36]X.G.He and L.Ren,A modified multiscale finite element method for well-driven flow problems in heterogeneous porous media,J.Hydrol.,2006,329(3-4):674-684.
[37]U.Hornung,Homogenization and porous media,Springer-Verlag,New York,1996.
[38]T.Y.Hou and X.H.Wu,A multiscale finite element method for elliptic problems in composite materials and porous media,J.Comput.Phys.,1997,134(1):168-189.
[39]T.Y.Hou,X.H.Wu and Z.Q.Cai,Convergence of a multiscale finite element method for elliptic problems with rapidly oscillating coefficients,Math.Comput.,1999,68(227):913-943.
[40]T.Y.Hou,X.H.Wu and Y.Zhang,Removing the cell resonance error in the multiscale finite element method via a Petrove-Galerkin formulation,Comm.Math.Sci.,2004,2(2):185-205.
[41]W.Z.Huang,Y.H.Ren and R.D.Russell,Moving mesh methods based on moving mesh partial differential equations,J.Comput.Phys.,1994,113(2):279-290.
[42]Y.Q.Huang,S.Shu and H.Y.Yu,A posteriori error estimate of generalized finite element method for elliptic problems with rough coefficients,Sino-Italian Symposium on Applied and Computational Mathematics,Beijing,1998,27-29.
[43]Y.Q.Huang and J.C.Xu,A partition of unity finite element method for elliptic problems with highly oscillating coefficients,in Proceedings for the Workshop on Scientific Computing,Hong Kong,1999,1-8.
[44]Y.Q.Huang and J.C.Xu,A conforming finite element method for overlapping and nonmatching grids,Math.Comput.,2002,72(243):1057-1066.
[45]T.J.R.Hughes,Multiscale phenomena:Green's functions,the Dirichlet-to-Neumann formulation,subgrid scale models,bubbles and the origins of stabilized methods,Comput.Methods Appl.Mech.Engrg.,1995,127(1):387-401.
[46]T.J.R.Hughes,G.R.Feijoo,L.Mazzei and J.B.Quincy,The variational multiscale method-a paradigm for computational mechanics,Comput.Methods Appl.Mech.Engrg.,1998,166(1):3-24.
[47]P.Jenny,S.H.Lee and H.A.Tchelepi,Multisale finite-volume method for elliptic problems in subsurface flow simulation,J.Comput.Phys.,2003,187(1):47-67.
[48]P.Jenny,S.H.Lee and H.A.Tchelepi,Adaptive multiscale finite-volume method for multiphase flow and transport in porous media,Multiscale Model.Simul.,2004,3(1):50-64.
[49]R.Krause and E.Rank,Multiscale computations with a combination of the hand p-versions finite element method,Comput.Methods Appl.Mech.Engrg.,2003,192(35):3959-3983.
[50]N.Kopteva and M.Stynes,A robust adaptive method for a quasilinear onedimensional convection-diffusion problem,SIAM J.Numer.Anal.,2001,39(4):1446-1467.
[51]T.LinB and M.Stynes,The SDFEM on Shishkin meshes for linear convectiondiffusion problems,Numer.Math.,2001,87(3):457-484.
[52]W.K.Liu,S.F.Li and T.Belytschko,Reproducing kernel particle method(I)Methodology and convergence,Comput.Methods Appl.Mech.Engrg.,1994,2:519-534.
[53]W.K.Liu,S.F.Li and T.Betytschko,Moving least-square reproducing kernel methods(I) Methodology and convergence,Comput.Methods Appl.Mech.Engrg.,1997,143(1-2):113-154.
[54]J.A.Mackenzie,Uniform convergence analysis of an upwind finite-difference approximation of a convection-diffusion boundary value problem on an adaptive grid,IMA J.Numer.Anal.,1999,19(2):233-249.
[55]S.Mehraeen and J.S.Chen,Wavelet-based multi-scale projection method in homogenization of heterogeneous media,Finite Elem.Anal.Des.,2004,40(12):1665-1679.
[56]J.M.Melenk and I.Babuska,The partition of malty finite element method:Basic theory and applications,Comput.Methods Appl.Mech.Engrg.,1996,139(1-4):289-314.
[57]P.B.Ming and X.Y.Yue,Numerical methods for multiscale elliptic problems,J.Comput.Phys.,2006,214(1):421-445.
[58]S.Moskow and M.Vogelius,First-order corrections to the homogenized eigenvalues of a periodic composite medium,SIAM J.Appl.Math.,1993,53(6):1636-1668.
[59]J.T.Oden,C.A.Duarte and O.C.Zienkiewicz,A new cloud-based hp finite element method,Comput.Methods Appl.Mech.Engrg.,1998,153(1):117-126.
[60]E.Onate and S.Idelsohn,A mesh-free finite point method for advective-diffusive transport and fluid flow problems,Comput.Mech.,1998,21(4-5):283-292.
[61]P.J.Park and T.Y.Hou,Multiscale numerical methods for singularly perturbed convection-diffusion equations,Int.J.Comput.Meth.,2004,1(1):17-65.
[62]S.Pruess and C.T.Fulton,Mathematical software for Sturm-Liouville problems,ACM Trans.Math.Soft.,1993,19(3):360-376.
[63]S.Pruess,C.T.Fulton and Y.T.Xie,An asymptotic numerical method for a class of singular Sturm-Liouville problems,SIAM J.Numer.Anal.,1995,32(5):1658-1676.
[64]A.Quarteroni and A.Valli,Numerical approximation of partial differential equations,Springer-Verlag,New York,1998.
[65]W.Q.Ren and W.N.E,Heterogeneous multiscale method for the modeling of complex fluids and micro-fluidics,J.Comput.Phys.,2005,204(1):1-26.
[66]H.G.Roos,M.Stynes and L.Tobiska,Numerical methods for singularly perturbed differential equations,The Springer Press,Berlin,1996.
[67]G.Sangalli,Capturing small scales in elliptic problems using a residual-free bubbles finite element method,Multiscale Model.Sim.,2003,1(3):485-503.
[68]M.Sarkis,Nonstandard coarse spaces and Schwarz methods for elliptic problems with discontinuous coefficients using non-conforming elements,Numer.Math.,1997,77(3):383-406.
[69]T.Stroubouiis,I.Babuska and K.Copps,The design and analysis of the generalized finite element method,Comput.Methods Appl.Mech.Engrg.,2000,181(1-3):43-69.
[70]T.Strouboulis,K.Copps and I.Babuska,The generalized finite element method,Comput.Methods Appl.Mech.Engrg.,2001,190(32-33):4081-4193.
[71]M.Stynes and E.O'Riordan,A uniformly convergent Galerkin method on a Shishkin mesh for a convection-diffusion problem,J.Math.Anal.Appl.,1997,214(1):36-54.
[72]H.Z.Tang and T.Tang,Adaptive mesh methods for one- and two-dimensional hyperbolic conservation laws,SIAM J.Numer.Anal.,2003,41(2):487-515.
[73]V.Thomee,From finite differences to finite elements:A short history of numerical analysis of partial differential equations,J.Comput.Appl.Math.,2001,128(1-2):1-54.
[74]W.X.Wang and Y.Takao,Isoparametric finite point method in computational mechanics,Comput.Mech.,2004,33(6):481-490.
[75]J.Xu,An introduction to multigrid convergence theory,Lecture notes in Winter School on Iterative Methods,1997,1-61.
[76]J.Xu and L.T.Zikatanov,On multigrid methods for generalized finite element methods,Comput.Sci.Eng.,2003,26:401-418.
[77]X.Y.Yue and W.N.E,The local microscale problem in the multiscale modeling of strongly heterogeneous media:Effects of boundary conditions and cell size,J.Comput.Phys.,2007,222(2):556-572.
[78]O.C.Zienkiewicz,The background of error estimation and adaptivity in finite element computations,Comput.Methods Appl.Mech.Engrg.,2006,195(4-6):207-213.
[79]陈传淼,黄云清,有限元高精度理论,湖南科学技术出版社,长沙,1995.
[80]都志辉,高性能计算并行编程技术-MPI并行程序设计,清华大学出版社,北京,2001.
[81]李荣华,偏微分方程数值解法,高等教育出版社,北京,2005.
[82]莫则尧,袁国兴,消息传递并行编程环境MPI,科学出版社,北京,2001.
[83]苏煜城,吴启光,奇异摄动问题数值方法引论,重庆出版社,重庆,1992.
[84]王烈衡,许学军,有限元方法的数学理论,科学出版社,北京,2004.
[85]张雄,刘岩,无网格法,清华大学出版社,北京,2004.