An Interpolated Galerkin Finite Element Method for the Poisson Equation

摘要

Abstract We develop a new approach to construct finite element methods to solve the Poisson equation. The idea is to use the pointwise Laplacian as a degree of freedom followed by interpolating the solution at the degree of freedom by the given right-hand side function in the partial differential equation. The finite element solution is then the Galerkin projection in a smaller vector space. This idea is similar to that of interpolating the boundary condition in the standard finite element method. Our approach results in a smaller system of equations and of a better condition number. The number of unknowns on each element is reduced significantly from (k2+3k+2)/2documentclass12pt{minimal} usepackage{amsmath} usepackage{wasysym} usepackage{amsfonts} usepackage{amssymb} usepackage{amsbsy} usepackage{mathrsfs} usepackage{upgreek} setlength{oddsidemargin}{-69pt} begin{document}$$(k^2+3k+2)/2$$end{document} to 3k for the Pkdocumentclass12pt{minimal} usepackage{amsmath} usepackage{wasysym} usepackage{amsfonts} usepackage{amssymb} usepackage{amsbsy} usepackage{mathrsfs} usepackage{upgreek} setlength{oddsidemargin}{-69pt} begin{document}$$P_k$$end{document} (k≥3documentclass12pt{minimal} usepackage{amsmath} usepackage{wasysym} usepackage{amsfonts} usepackage{amssymb} usepackage{amsbsy} usepackage{mathrsfs} usepackage{upgreek} setlength{oddsidemargin}{-69pt} begin{document}$$kge 3$$end{document}) finite element. We construct bivariate P2documentclass12pt{minimal} usepackage{amsmath} usepackage{wasysym} usepackage{amsfonts} usepackage{amssymb} usepackage{amsbsy} usepackage{mathrsfs} usepackage{upgreek} setlength{oddsidemargin}{-69pt} begin{document}$$P_2$$end{document} conforming and nonconforming, and Pkdocumentclass12pt{minimal} usepackage{amsmath} usepackage{wasysym} usepackage{amsfonts} usepackage{amssymb} usepackage{amsbsy} usepackage{mathrsfs} usepackage{upgreek} setlength{oddsidemargin}{-69pt} begin{document}$$P_k$$end{document} (k≥3documentclass12pt{minimal} usepackage{amsmath} usepackage{wasysym} usepackage{amsfonts} usepackage{amssymb} usepackage{amsbsy} usepackage{mathrsfs} usepackage{upgreek} setlength{oddsidemargin}{-69pt} begin{document}$$kge 3$$end{document}) conforming interpolated Galerkin finite elements on triangular grids; prove their optimal order of convergence; and confirm our findings by numerical tests.