deal.II Implementation of a Weak Galerkin Finite Element Solver for Darcy Flow
- Colorado State Univ., Fort Collins, CO (United States); Kitware, Inc.
- Colorado State Univ., Fort Collins, CO (United States)
- Kitware, Inc., Clifton Park, NY (United States)
This paper presents a weak Galerkin (WG) finite element solver for Darcy flow and its implementation on the deal.II platform. The solver works for quadrilateral and hexahedral meshes in a unified way. It approximates pressure by Q-type degree $k(≥ 0)$ polynomials separately defined in element interiors and on edges/faces. In this work, numerical velocity is obtained in the unmapped Raviart-Thomas space $$RT_{[k]}$$ via postprocessing based on the novel concepts of discrete weak gradients. The solver is locally mass-conservative and produces continuous normal fluxes. The implementation in deal.II allows polynomial degrees up to 5. Numerical experiments show that our new WG solver performs better than the classical mixed finite element methods.
- Research Organization:
- Kitware, Inc., Clifton Park, NY (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR) (SC-21); USDOE Office of Science (SC), Basic Energy Sciences (BES) (SC-22). Chemical Sciences, Geosciences & Biosciences Division; National Science Foundation (NSF)
- DOE Contract Number:
- SC0019609
- OSTI ID:
- 1594093
- Resource Type:
- Conference paper/presentation
- Conference Information:
- International Conference on Computational Science, Faro (Portugal), 12-14 Jun 2019; From, proceedings, part IV
- Country of Publication:
- United States
- Language:
- English
Similar Records
Weak Galerkin finite element methods for Darcy flow: Anisotropy and heterogeneity
Journal Article
·
Sat Nov 01 00:00:00 EDT 2014
· Journal of Computational Physics
·
OSTI ID:1236930