skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids: LDRD Final Report.

Abstract

This report summarizes the work performed under a one-year LDRD project aiming to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. In this project, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality. At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokesmore » problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5x improvement in accuracy and allows for taking an 8x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method. The report concludes with a brief summary of ongoing projects and collaborations that utilize or extend the products of this work.« less

Authors:
; ; ;
Publication Date:
Research Org.:
Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1669732
Report Number(s):
SAND2020-9980
690989
DOE Contract Number:  
AC04-94AL85000
Resource Type:
Technical Report
Country of Publication:
United States
Language:
English

Citation Formats

Kuberry, Paul Allen, Bochev, Pavel B., Koester, Jacob, and Trask, Nathaniel Albert. Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids: LDRD Final Report.. United States: N. p., 2020. Web. doi:10.2172/1669732.
Kuberry, Paul Allen, Bochev, Pavel B., Koester, Jacob, & Trask, Nathaniel Albert. Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids: LDRD Final Report.. United States. https://doi.org/10.2172/1669732
Kuberry, Paul Allen, Bochev, Pavel B., Koester, Jacob, and Trask, Nathaniel Albert. 2020. "Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids: LDRD Final Report.". United States. https://doi.org/10.2172/1669732. https://www.osti.gov/servlets/purl/1669732.
@article{osti_1669732,
title = {Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids: LDRD Final Report.},
author = {Kuberry, Paul Allen and Bochev, Pavel B. and Koester, Jacob and Trask, Nathaniel Albert},
abstractNote = {This report summarizes the work performed under a one-year LDRD project aiming to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. In this project, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality. At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5x improvement in accuracy and allows for taking an 8x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method. The report concludes with a brief summary of ongoing projects and collaborations that utilize or extend the products of this work.},
doi = {10.2172/1669732},
url = {https://www.osti.gov/biblio/1669732}, journal = {},
number = ,
volume = ,
place = {United States},
year = {2020},
month = {9}
}