Abstract:A variational approach is developed with a meshless discretization 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. We propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality, which we call discontinuous piecewise polynomial generalized moving least squares (DPP-GMLS). 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 5 improvement in accuracy and allows for taking an 8 larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method.

An optimization-based approach for high-order accurate discretization of conservation laws with discontinuous solutions

An Adjoint Method for a High-Order Discretization of Deforming Domain Conservation Laws for Optimization of Flow Problems

Conservative Front Tracking and Level Set Algorithms

Continuously bounds-preserving discontinuous Galerkin methods for hyperbolic conservation laws

A Time-Accurate Dynamic Mesh Approach for High-Order Shock Computation

Optimal error estimates for discontinuous Galerkin methods based on upwind-biased fluxes for linear hyperbolic equations.

An extended discontinuous Galerkin shock tracking method

High order exactly divergence-free Hybrid Discontinuous Galerkin Methods for unsteady incompressible flows

A discontinuous piecewise polynomial generalized moving least squares scheme for robust finite element analysis on arbitrary grids

A data-driven learned discretization approach in finite volume schemes for hyperbolic conservation laws and varying boundary conditions

Dissipative WENO stabilization of high-order discontinuous Galerkin methods for hyperbolic problems

An SBP-SAT Continuous Galerkin Finite Element Formulation for Smooth and Discontinuous Fields

Multi-Level Solution Strategies for Implicit Time Discontinuous Galerkin Discretizations

High Resolution Optimized High-Order Schemes for Discretization of Non-Linear Straight and Mixed Second Derivative Terms

A bound preserving cut discontinuous Galerkin method for one dimensional hyperbolic conservation laws

Discontinuous Galerkin schemes for hyperbolic systems in non-conservative variables: quasi-conservative formulation with subcell finite volume corrections

Stability and Error Analysis of an Efficient Numerical Method for Convection Dominated Parabolic PDEs with Jump Discontinuity in Source Function on Modified Layer-Adapted Mesh

Development of discontinuous Galerkin methods for hyperbolic systems that preserve a curl or a divergence constraint

High-order Discontinuous Galerkin Solver on Hybrid Anisotropic Meshes for Laminar and Turbulent Simulations