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.

Piecewise equidistant meshes for quasilinear turning point problems: Technical report

A Meshless Method for Some Inverse Problems Associated with the Helmholtz Equation

Positive Solutions to a Singular Second Order Three-Point Boundary Value Problem

Analysis of Galerkin and streamline-diffusion FEMs on piecewise equidistant meshes for turning point problems exhibiting an interior layer

A Uniform Convergent Petrov-Galerkin Method for a Class of Turning Point Problems

Improvement of boundary aberration and computation efficiency in meshless method

Optimally accurate second-order numerical scheme on layer-resolving graded mesh for singularly-perturbed boundary turning point problems exhibiting power-of-type-1 and hybrid layers

A Novel Volume Constrained Smoothing Method for Meshes

A parameter-uniform hybrid scheme designed for multi-point boundary value problems that are perturbed

A fully discrete scheme on piecewise-equidistant mesh for singularly perturbed delay integro-differential equations

Simple Positivity Preserving Nonlinear Finite Volume Scheme for Subdiffusion Equations on General Non-conforming Distorted Meshes

Error analysis of discontinuous Galerkin methods on layer adapted meshes for two dimensional turning point problem

A hybrid approach for singularly perturbed parabolic problem with discontinuous data

The Helmholtz Dirichlet and Neumann problems on piecewise smooth open curves

Theoretical and numerical analysis of problems with a boundary turning point and a variable diffusion coefficient

An efficient fourth-order convergent scheme based on half-step spline function for two-point mixed boundary value problems

Linear regularized finite difference scheme for the quasilinear subdiffusion equation

A parameter uniform hybrid approach for singularly perturbed two-parameter parabolic problem with discontinuous data

Quasi-Optimal Meshes for Gradient Nonconforming Approximations

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

Numerical Solution of Quenching Problems Using Mesh-Dependent Variable Temporal Steps