H(div)-椭圆问题的若干数值方法研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
在本文中,我们研究H(div)-椭圆问题的若干数值方法。据我们所知,H(div)-椭圆问题在固体和流体力学中是一个很普遍的问题,并且在实际问题中有很重要的应用。对二阶椭圆问题做一阶最小二乘时,在实现对于非静态不可压Navier-Stokes方程组序列正则化的过程中,在增广Lagrangian混合有限元方法中,在Stokes方程组的有限元稳定化过程中,都有可能遇到H(div)-椭圆问题的求解。因此,设计数值算法离散和求解该问题是有必要的,也是有意义的。
     第一章,我们给出文中其他章节会采用的一些预备知识。
     第二章,我们设计了内罚型间断Galerkin方法离散H(div)、椭圆问题。首先给出了在能量范数意义下的最优先验误差估计。并且,我们给出了一类残量型后验误差估计子,证明了该误差估计子的可靠性和有效性。最后,数值实验验证了我们的理论结果。
     第三章,我们设计了一类并行的Robin-Robin区域分解算法求解H(div)-椭圆问题。对偏微分方程和有限元逼近问题,都给出了收敛性分析。最后,数值实验说明了该算法的有效性。
     第四章,我们设计了最优型区域分解算法求解H(div)-椭圆问题。通过适当参数的选取,论证了算法的收敛性。数值实验验证了算法的有效性。
     第五章,我们提供了关于非标准有限元离散H(div)-椭圆问题的后验误差框架。我们将该框架应用于第二章的内罚型间断Galerkin方法。并且说明该框架也可以应用于Mortar有限元方法。
     第六章,我们设计了非拟合罚有限元方法离散H(div)-椭圆界面问题。给出了能量范数意义下的最优先验误差估计。
In this dissertation, we study some numerical methods for H(div)-elliptic prob-lem. As we know, H(div)-elliptic problem is ubiquitous in solid and fluid mechanics, and have many important real applications. It may arises from, the first-order system least-squares formulation of H1-elliptic problem, the implementation of the sequential regularization method for the nonstationary incompressible Navier-Stokes equations, the mixed methods with augmented lagrangians, or the stabilized formulations of the Stokes equations. Therefore, it is so important and necessary to design numerical methods to solve this equation.
     In Chapter1, we provide some preliminaries which may be applied to the re-maining chapters throughout the dissertation.
     In Chapter2, we propose and analyze the interior penalty discontinuous Galerkin method for H(div)-elliptic problem. An optimal a priori error estimate in the energy norm is proved. In addition, a residual-based a posteriori error estimator is obtained. The estimator is proved to be both reliable and efficient in the energy norm. The numerical results presented validate our theoretical analysis.
     In Chapter3, we present a parallel Robin-Robin domain decomposition method (DDM) for H(div)-elliptic problem. The convergence of this method is proved for the continuous problem and the finite element discrete approximation problem. Some numerical testes are presented to demonstrate the effectiveness of the method.
     In Chapter4, we propose some optimization-based domain decomposition meth-ods for H(div)-elliptic problem. Convergent properties are examined by choosing proper parameters. The effectiveness of the method is validated by some numerical tests.
     In Chapter5, we provide an unified framework for an a posteriori error analysis of non-standard finite element approximations of H(div)-elliptic problem. We apply it to the interior penalty discontinuous Galerkin method in Chapter2. Furthermore, we can apply it to the mortar finite element method.
     In Chapter6, we propose an unfitted penalty finite element method for H(div)-elliptic interface problem. An optimal a priori error estimate in energy norm is ob- tained.
引文
[1]A. B. Abdallah, F. B. Belgacem and Y. Maday, Mortaring the two-dimensional Nedelec finite elements for the discretization of the Maxwell equations, Math. Models Methods Appl. Sci.,14 (2004),1635-1656.
    [2]Y. Achdoua, P. L. Tallecb, F. Natafc and M. Vidrascud, A domain decompo-sition preconditioner for an advection-diffusion problem, Comput. Methods. Appl. Mech. Eng,184(2000),145-170.
    [3]E. Ahusborde, M. Azaiez, M. O. Deville and E. H. Mund, An iterative domain decomposition algorithm for the grad(div) operator, Commun. Comput. Phys., 5 (2009),391-397.
    [4]M. Ainsworth, A posteriori error estimation for discontinuous Galerkin finite element approximation, SIAM J. Numer. Anal.,45(2007),1777-1798.
    [5]M. Ainsworth and J. T. Oden, A Posteriori Error Estimation in Finite Element Analysis, Wiley-Interscience, JohnWiley & Sons, New York,2000.
    [6]D. N. Arnold, An interior penalty finite element method with discontinuous elements, SIAM J. Numer. Anal.,19(1982),742-760.
    [7]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.
    [8]D. N. Arnold, R. S. Falk and R. Winther, Multigrid preconditioning in H(div) and application, Math. Comp.,66(1997),957-984.
    [9]D. N. Arnold, R. S. Falk and R. Winther, Multigrid in H(div) and H(curl), Numer. Math.,85 (2000),197-217.
    [10]I. Babuska, The finite element method for elliptic equations with discontinu-ous coeffcients, Computing,5 (1970),207-218.
    [11]I. Babuska and T. Strouboulis, The Finite Element Method and its Reliability, Clarendon Press, Oxford University Press, New York,2001.
    [12]J. W. Barrett and C. M. Elliott, Fitted and unfitted finite-element methods for elliptic equations with smooth interfaces, IMA J. Numer. Anal.,7 (1987), 283-300.
    [13]P. Bastian and C. Engwer, An unfitted finite element method using discontin-uous Galerkin, Int. J. Numer. Meth. Engng,79 (2009),1557-1576.
    [14]R. Becker, P. Hansbo and M. Larson, Energy norm a posteriori error esti-mation for discontinuous Galerkin methods, Comput. Methods Appl. Mech. Engrg.,192(2003),723-733.
    [15]R. Beck, R. Hiptmair, R. H. W. Hoppe and B. I. Wohlmuth, Residual based a-posteriori error estimators for eddy current computation, M2AN. Math. Model. Numer. Anal.,34(2000),159-182.
    [16]F. B. Belgacem, The mortar finite element method with Lagrange multipliers, Numer. Math.,84 (1999),173-198.
    [17]F. B. Belgacem, A. Buffa and Y. Maday, The mortar finite element method for 3D Maxwell equations:first results., SIAM J. Numer. Anal.,39 (2001), 880-901.
    [18]F. B. Belgacem and Y. Maday, The mortar element method for three dimen-sional finite elements, RAIRO Anal. Numer.,31 (1997),289-302.
    [19]S. Bennethum and X. Feng, A domain decomposition method for solving a Helmholz-like problem in elasticity based on the Wilson nonconfoming ele-ment, M2AN. Math. Model. Numer. Anal.,31 (1997) 1-25.
    [20]C. Bernardi, Y. Maday and A. T. Patera, A new nonconforming approach to domain decomposition:The mortar element method, in Nonlinear Partial Dif-ferential Equations and Their Applications, College de France Seminar, Pit-man Res. Notes inMath.9, H. Brezis and J. L. Lions, eds., Longman Scientific and Technical, Harlow, UK,1994.
    [21]J. H. Bramble and J. T. King, A finite element method for interface problems in domains with smooth boundaries and interfaces, Adv. Comput. Math.,6 (1996),109-138.
    [22]S. C. Brenner and L. R. Scott, The Mathematical Theory of Finite Element Methods (3nd edn), Springer-Verlag, New York/Berlin/Heidelberg,2008.
    [23]D. Bresch and J. Koko, Oprator-splitting and lagranger multiplier domain de-compsition methods for numerical simulation of two coupled Navier-Stokes fluids, Int J. Appl. Math. Comput. Sci.,16 (2006),419-429.
    [24]F. Brezzi, J. Douglas and L. D. Marini, Two families of mixed finite elements for second order elliptic problems, Numer. Math.,47(1985),217-235.
    [25]F. Brezzi and M. Fortin, Mixed and Hybrid Finite Element Methods, Springer-Verlag, New York,1991.
    [26]Z. Cai, R. Lazarov, T. Manteuffel and S. McCormick, First-order system least-squares for partial differential equations, Part I. SIAM J. Numer. Anal., 31(1994),1785-1799.
    [27]C. Carstensen, A unifying theory of a posteriori finite element error control, Numer. Math.,100 (2005),617-637.
    [28]C. Carstensen, T. Gudi and M. Jensen, A unifying theory of a posteriori error control for discontinuous Galerkin FEM, Numer. Math.,112(2009),363-379.
    [29]C. Carstensen and R. H. W. Hoppe, Convergence analysis of an adaptive edge finite element method for the 2d eddy current equations, J. Numer. Math., 13(2005),19-32.
    [30]C. Carstensen and R. H. W. Hoppe, Unified framework for an a posteriori er-ror analysis of non-standard finite element approximations of H(curl)-elliptic problems, J. Numer. Math.,17 (2009),27-44.
    [31]C. Carstensen and J. Hu, A unifying theory of a posteriori error control for nonconforming finite element methods, Numer. Math.,107(2007),473-502.
    [32]C. Carstensen, J. Hu and A. Orlando, Framework for the a posteriori error analysis of nonconforming finite elements, SIAM J. Numer. Anal.,45(2007), 68-82.
    [33]J. M. Cascon, R. H. Nochetto and K. G. Siebert, Design and convergence of AFM in H(div), Math. Models Methods Appl. Sci.,17(2007),1849-1881.
    [34]L. Chen and C. Zhang, A MATLAB Package of Adaptive Finite Element Methods, Technical report, University of Maryland, College Park, MD,2006.
    [35]Z. Chen, Y. Xiao and L. Zhang, The adaptive immersed interface finite ele-ment method for elliptic and Maxwell interface problems, J. Comput. Phys., 228 (2009),5000-5019.
    [36]Z. Chen and J. Zou, Finite element methods and their convergence for elliptic and parabolic interface problems, Numer. Math.,79 (1998),175-202.
    [37]C. Chu, I. G. Graham and T. Hou, A new multiscale finite element method for high-contrast elliptic interface problems, Math. Comp.,79 (2010),1915-1955.
    [38]P. Ciarlet, The Finite Element Method for Elliptic Problems, North-Holland, Amsterdam,1978.
    [39]P. Clement, Approximation by finite element functions using local regulariza-tion, RAIRO Anal. Numer.,2 (1975),77-84.
    [40]B. Cockburn, G. E. Karniadakis and C. W. Shu (Eds.), Discontinuous Galerkin Methods-Theory, Computation and Applications, Lecture Notes in Computa-tional Science and Engineering 11. Springer-Verlag, New York,2000.
    [41]M. Crouzeix and A. Raviart, Conforming and nonconforming finite element methods for solving the stationary Stokes equation, RAIRO. Anal. Numer, 7(1973),33-76.
    [42]Q. Deng, An analysis for a nonoverlapping domain decomposition iterative procedure, SIAM J. Sci. Comput.,18 (1997),1517-1525.
    [43]Q. Deng, A nonoverlapping domain decomposition method for nonconform-ing finte element problems. Commun. Pure Appl. Math.,2 (2003),297-310.
    [44]B. Despres, Domain decomposition method and Helmholtz problem, in Math-ematical and Numerical Aspects of Wave Propagation Phenomena, G. Cohen, L. Halpern, and P. Jply, eds. SIAM, Philadephia, PA,1991,44-51.
    [45]Q. Du, Optimization based nonoverlapping domain decomposition algorithms and their convergence, SIAM J. Numer. Anal.,39 (2001),1056-1077.
    [46]Q. Du and M. D. Gunzburger, A gradient method approach to optimization-based multidisciplinary simulations and nonoverlapping domain decomposi-tion algorithms, SIAM J. Numer. Anal.,37 (2000),1513-1541.
    [47]X. Feng, Analysis of a domain decomposition method for the nearly elastic wave equations based on mixed finite element methods, IMA J. Numer. Anal., 18 (1998),229-250.
    [48]M. J. Gander, Optimized Schwarz methods, SIAM J. Numer. Anal.,44 (2006), 699-731.
    [49]V. Girault and P. A. Raviart, Finite element methods for Navier-Stokes equa-tions, Springer-Verlag, New York,1986.
    [50]Y. Gong, B. Li and Z. Li, Immersed-interface finite-element methods for el-liptic interface problems with nonhomogeneous jump conditions, SIAM J. Numer. Anal.,46 (2008),472-495.
    [51]J. Gopalakrishnan and J. P. Pasciak, Multigrid for the mortar finite element method, SIAM J. Numer. Anal.,37 (2000),1029-1052.
    [52]M. D. Gunzburger, M. Heikenschloss and H. K. Lee, Solution of elliptic par-tial differential equations by an optimization-based domaind decompostion method, Appl. math, comput.,113 (2000),111-139.
    [53]M. D. Gunzburger and H. K. Lee, An optimization-based domain decompo-sition method for the Navier-Stokes equations, SIAM J. Numer. Anal.,37 (2000),1455-1480.
    [54]M. D. Gunzburger, J. Peterson and H. K. Lee, An optimization based do-main decomposition method for partial differential equations, Comput. Math. Appl.,37 (1999),77-93.
    [55]W. Guo and L. S. Hou, Generalizations and accelerations of Lions'nonover-lapping domain decomposition method for linear elliptic PDE, SIAM J. Nu-mer. Anal.,41 (2003),2056-2080.
    [56]G. Guyomarc'h, C. O. Lee and K. Jeon, A discontinuous Galerkin method for elliptic interface problems with application to electroporation, Commun. Numer. Meth. Engng.,25 (2009),991-1008.
    [57]A. Hansbo and P. Hansbo, An unfitted finite element method, based on Nitsche's method, for elliptic interface problems, Comput. Methods Appl. Mech. Engrg.,191 (2002),5537-5552.
    [58]R. Hiptmair, Multigrid method for H(div) in three dimensions, Electron. Tran. Numer. Anal.,6 (1997),133-152.
    [59]R. Hiptmair, J. Li and J. Zou, Convergence analysis of finite element methods for H(div;Ω)-elliptic interface problems. J. Numer. Math.,18 (2010),187-218.
    [60]R. Hiptmair, J. Li and J. Zou, Convergence analysis of finite element methods for H(curl;Ω)-elliptic interface problems. Numer. Math., to appear.
    [61]R. Hiptmair and A. Toselli, Overlapping and multilevel Schwarz methods for vector valued elliptic problems in three dimensions, in Parallel Solution of PDEs, IMA Volumes in Mathematics and its Applications, Berlin,2000, Springer-Verlag.
    [62]R. Hiptmair and J. Xu, Nodal auxiliary space preconditioning in H(curl) and H(div) spaces. SIAM J. Numer. Anal.,45(2007),2483-2509.
    [63]R. H. W. Hoppe, Mortar edge elements in R3. East-West J. Numer. Math., 7(1999),159-173.
    [64]P. Houston, I. Perugia, A. Schneebeli and D. Schotzau, Interior penal-tymethod for the indefinite time-harmonic Maxwell equations. Numer. Math., 100(2005),485-518.
    [65]P. Houston, I. Perugia and D. Schotzau, An a posteriori error indicator for discontinuous Galerkin discretizations of H(curl)-elliptic partial differential equations, IMA J. Numer. Anal.,27 (2007),122-150.
    [66]P. Houston, D. Schotzau and T. P. Wihler, Energy norm a posteriori error esti-mation of hp-adaptive discontinuous Galerkin methods for elliptic problems, Math. Models Methods Appl. Sci.,17 (2007),33-62.
    [67]Q. Hu, A non-conforming domain decomposition for Raviart-Thomas vector fields in three dimensions, Science in China Series A:Mathematics,47(2004), 190-199.
    [68]Q. Hu, S. Shu and J. Zou, A mortar edge element method with nearly op-timal convergence for three-dimensional Maxell's equations, Math. Comp., 77(2008),1333-1353.
    [69]H. Huang and Z. Li, Convergence analysis of the immersed interface method, IMA J. Numer. Anal.,19 (1999),583-608.
    [70]J. Huang and J. Zou, A mortar element method for elliptic problems with discontinuous coefficients, IMA J. Numer. Anal.,45 (2002),549-576.
    [71]O. A. Karakashian and F. Pascal, A posteriori error estimates for a discon-tinuous Galerkin approximation of second-order elliptic problems, SIAM J. Numer. Anal.,41 (2003),2374-2399.
    [72]F. Kong, Y. Ma and J. Lu, An Optimization-based domain decomposition method for numerical simulation of the incompressible Navier-Stokes flows, Numer Methods Partial Differential Eq,27 (2011),255-276.
    [73]R. LeVeque and Z. Li, The immersed interface method for elliptic equations with discontinuous coeffcients and singular sources, SIAM J. Numer. Anal., 31 (1994),1019-104.
    [74]P. Lin, A sequential regularization method for time-dependent incompressible Navier-Stokes equations, SIAM J. Numer. Anal.,34 (1997),1051-1071.
    [75]P. L. Lions, On the Schwarz alternating method. Ⅲ:A variant for nonoverlap-ping subdomains, in Third International Symposium on Domain Decomposi-tion Mehthods for Partial Differential Equations, T.F. Chan, R.Glowinski, J. Periaux, and O.B. Widlund, eds. SIAM, Philadephia, PA,1990,202-223.
    [76]J. Li, J. M. Melenk, B. Wohlmuth and J. Zou, Optimal a priori estimates for higher order finite elements for elliptic interface problems, Appl. Numer. Math.,60 (2010),19-37.
    [77]Z. Li and K. Ito, The Immersed Interface Method:Numerical solutions of PDE Involving Interface and Irregular Domains, SIAM, Philadephia,2006.
    [78]Z. Li, T. Lin and X. Wu, New cartesian grid methods for interface problems using the finite element formulation, Numer. Math.,96 (2003),61-98.
    [79]M. N. Linnick and H. F. Fasel, A high-order immersed interface method for simulating unsteady incompressible flows on irregular domains, J. Comput. Phys.,204 (2005),157-192.
    [80]S. H. Lui, A Lions nonoverlapping domain decomposition method for do-mains with an arbitrary interface, IMA J. Numer. Anal.,29 (2009),332-349.
    [81]L. Marcinkowski, The mortar element method with locally nonconforming elements. BIT,39 (1999),716-739.
    [82]R. Massjung, An hp-error estimate for an unfitted discontinuous Galerkin method applied to elliptic interface problems, RWTH 300, IGPM Report, 2009.
    [83]P. Monk, Finite Element Methods for Maxwell's Equation, Oxford University press, New York,2003.
    [84]J. C. Nedelec, Mixed finite elements in R3, Numer. Math.,35(1980),315-341.
    [85]J. C. Nedelec, A new family of mixed finite elements in R3, Numer. Math., 50(1986),57-81.
    [86]P. Neittaanmaki and S. Repin, Reliable Methods for Mathematical Modelling. Error Control and a Posteriori Estimates, Elsevier, New York,2004.
    [87]I. Perugia and D. Schotzau, The hp-local discontinuous Galerkin method for low-frequency time-harmonic Maxwell equations. Math. Comp.,72(2003), 1179-1214.
    [88]D. Pradhan, B. Shalini, N. Nataraj and A. K. Pani, A Robin-type non-overlapping domain decomposition procedure for second order elliptic prob-lems, Adv. Comput. Math,34(2011),339-368.
    [89]L. Qin and X. Xu, On a parallel Robin-type nonoverlapping domain decom-position method, SIAM J. Numer. Anal.,44 (2006),2539-2558.
    [90]L. Qin, Z. Shi and X. Xu, On the convergence rate of a parallel nonoverlapping domain decomposition method, Science in China Series A:Mathematics,51 (2008),1461-1478.
    [91]A. Quateroni and A. Valli, Domain Decomposition Methods for Partial Equa-tions, Oxford University Press, Oxfrod,1999.
    [92]D. S. Oh, Domain Decomposition Methods for Raviart-Thomas Vector Fields, Thesis, New York University,2011.
    [93]R. Rannacher and S. Turek, A simple nonconforming quadrilateral Stokes el-ement, Numer Methods Partial Differential Eq,8(1992),97-111.
    [94]P. A. Raviart and J. M. Thomas, A mixed finite element method for second order elliptic problems, Lecture Notes in Math., Berlin-Heidelberg Newyork Springer, Vol.66,1977,292-315.
    [95]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.
    [96]R. Schneider, Y. Xu and A. Zhou, An analysis of discontinuous Galerkin meth-ods for elliptic problems, Adv. Comput. Math.,25(2006),259-286.
    [97]J. Schoberl, A posteriori error estimates for Maxwell equations, Math. Comp., 77(2008),633-649.
    [98]L. R. Scott and S. Zhang, Finite element interpolation of nonsmooth functions satisfying boundary conditions, Math. Comp.,54(1994),483-493.
    [99]B. F. Smith, P. Bjorstad, W.D. Gropp, Domain Decomposition:Parallel Multi-level Algorithms for Elliptic Partial Equations, Cambridge Univerdsity Press, New York,1996.
    [100]S. K. Tomar and S. I. Repin, Efficient computable error bounds for discon-tinuous Galerkin pproximations of elliptic problems, J. Comput. Appl. Math., 226(2009),358-369.
    [101]A. Toselli, Neuman-Neuman methods for vector field problems. Electron. Tran. Numer. Anal.11 (2000) 1-24.
    [102]A. Toselli and O. B. Wildlund, Domain Decomposition Methods-Algorithms and Theory, Springer-Verlag, Berlin,2005.
    [103]R. Verfurth, A Review of A Posteriori Estimation and Adaptive Mesh-Refinement Techniques, Wiley-Teubner, New York, Stuttgart,1996.
    [104]B. I. Wohlmuth and R. H. W. Hoppe, A comparison of a posteriori error esti-mators for mixed finite element discretizations by Raviart-Thomas elements, Math. Comp.,68 (1999),1347-1378.
    [105]B. I. Wohlmuth, A. Toselli and O. B. Widlund, An iterative substructuring method for Raviart-Thomas vector fields in three dimensions, SIAM J. Numer. Anal.,37 (2000),1657-1676.
    [106]H. Wu and Y. Xiao, An unfitted hp-interface penalty finite elemen method for elliptic interface problems. arXiv:1007.2893 [math.NA].
    [107]J. Xu and J. Zou, Some nonoverlapping domain decomposition methods, SIAM Rev.,40 (1998),857-914.
    [108]X. Xu and J. Chen, Multigrid for the mortar element method for P1 noncon-forming element, Numer. Math.,88(2001),381-398.
    [109]X. Xu and R. H. W. Hoppe, On the convergence of mortar edge element meth-ods in R3. SIAM J. Numer. Anal.,43(2005),1276-1294.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700