A weak Galerkin finite element method with polynomial reduction

• Lin Mu ^[2] ; Junping Wang ^[3] ; Xiu Ye ^[1]
  1. [1] University of Arkansas at Little Rock

     University of Arkansas at Little Rock

     Township of Big Rock, Estados Unidos

     
  2. [2] Michigan State University

     Michigan State University

     City of East Lansing, Estados Unidos

     
  3. [3] National Science Foundation, Arlington, USA

  Mostrar afiliaciones +
• Localización: Journal of computational and applied mathematics, ISSN 0377-0427, Vol. 285, Nº 1 (September 2015), 2015, págs. 45-58
• Idioma: inglés
• Enlaces
  • Texto completo
• Resumen

  • The weak Galerkin (WG) is a novel numerical method based on variational principles for weak functions and their weak partial derivatives defined as distributions. In the implementation, the weak partial derivatives and the weak functions are approximated by polynomials with various degrees of freedom. The accuracy and the computational complexity of the corresponding WG scheme is significantly impacted by the selection of such polynomials. This paper presents an optimal combination for the polynomial spaces that minimize the number of unknowns in the numerical scheme without compromising the accuracy of the numerical approximation. For illustrative purpose, the authors use the second order elliptic equation to demonstrate the basic ideas of polynomial reduction. Consequently, a new weak Galerkin finite element method is proposed and analyzed. Error estimates of optimal order are established for the corresponding WG approximations in both a discrete H1 norm and the standard L2 norm. In addition, the paper presents some numerical results to demonstrate the power of the WG method in dealing with finite element partitions with arbitrary polygons in 2D or polyhedra in 3D. The numerical examples include various finite element partitions such as triangular mesh, quadrilateral mesh, honeycomb mesh in 2D and mesh with deformed cubes in 3D.