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.

Analysis of a $\boldsymbol{P}_1\oplus \boldsymbol{RT}_0$ finite element method for linear elasticity with Dirichlet and mixed boundary conditions

Analysis of a finite element method for linear elasticity with Dirichlet and mixed boundary conditions

A stabilized $P_1$ immersed finite element method for the interface elasticity problems

THE UPPER BOUND PROPERTY FOR SOLID MECHANICS OF THE LINEARLY CONFORMING RADIAL POINT INTERPOLATION METHOD (LC-RPIM)

Lower Order Rectangular Nonconforming Mixed Finite Elements for Plane Elasticity

Efficient and reliable divergence-conforming methods for an elasticity-poroelasticity interface problem

New lower order mixed finite element methods for linear elasticity

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

A new mixed finite element for the linear elasticity problem in 3D

A new family of efficient conforming mixed finite elements on both rectangular and cuboid meshes for linear elasticity in the symmetric formulation

A Family of Discontinuous Galerkin Mixed Methods for Nearly and Perfectly Incompressible Elasticity

Lower order mixed elements for the linear elasticity problem in 2D and 3D

An Analysis of broken $P_1$-Nonconforming Finite Element Method For Interface Problems

Comments on 'upper Bound Solution to Elasticity Problems: A Unique Property of the Linearly Conforming Point Interpolation Method (Lc-Pim)'

Intrinsic mixed finite element methods for linear Cosserat elasticity and couple stress problem

The least squares finite element method for elasticity interface problem on unfitted mesh

Lowest-degree robust finite element schemes for inhomogeneous bi-Laplace problems

Low-degree robust Hellinger-Reissner finite element schemes for planar linear elasticity with symmetric stress tensors

Geometrical mixed finite element methods for fourth order obstacle problems in linearised elasticity

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