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.

Robust Weak Galerkin Finite Element Methods for Linear Elasticity with Continuous Displacement Trace Approximation

Simplified Weak Galerkin Methods for Linear Elasticity on Nonconvex Domains

New Discontinuous Galerkin Algorithms and Analysis for Linear Elasticity with Symmetric Stress Tensor

A Mixed Discontinuous Galerkin Method for Linear Elasticity with Strongly Imposed Symmetry

A Locking-Free Weak Galerkin Finite Element Method for Linear Elasticity Problems

An arbitrary order locking-free weak Galerkin method for linear elasticity problems based on a reconstruction operator

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

Generalized weak Galerkin finite element methods for second order elliptic problems

Full weak Galerkin finite element discretizations for poroelasticity problems in the primal formulation

The interpolating element-free Galerkin method for elastic large deformation problems

Weight-adjusted discontinuous Galerkin methods: matrix-valued weights and elastic wave propagation in heterogeneous media

On the enriched mixed Galerkin method with θ scheme for the elastic wave equation

A locking-free mixed enriched Galerkin method of arbitrary order for linear elasticity using the stress-displacement formulation

Robust globally divergence-free Weak Galerkin finite element method for incompressible Magnetohydrodynamics flow

Convergence and optimality of an adaptive modified weak Galerkin finite element method

Auto-Stabilized Weak Galerkin Finite Element Methods for Biharmonic Equations on Polytopal Meshes without Convexity Assumptions

$C^0$ discontinuous Galerkin finite element methods for second order linear elliptic partial differential equations in non-divergence form

A new discontinuous Galerkin spectral element method for elastic waves with physically motivated numerical fluxes

Auto-Stabilized Weak Galerkin Finite Element Methods on Polytopal Meshes without Convexity Constraints

A robust lower order mixed finite element method for a strain gradient elasticity model