Abstract:This work presents an abstract framework for the design, implementation, and analysis of the multiscale spectral generalized finite element method (MS-GFEM), a particular numerical multiscale method originally proposed in [I. Babuska and R. Lipton, Multiscale Model.\;\,Simul., 9 (2011), pp.~373--406]. MS-GFEM is a partition of unity method employing optimal local approximation spaces constructed from local spectral problems. We establish a general local approximation theory demonstrating exponential convergence with respect to local degrees of freedom under certain assumptions, with explicit dependence on key problem parameters. Our framework applies to a broad class of multiscale PDEs with $L^{\infty}$-coefficients in both continuous and discrete, finite element settings, including highly indefinite problems (convection-dominated diffusion, as well as the high-frequency Helmholtz, Maxwell and elastic wave equations with impedance boundary conditions), and higher-order problems. Notably, we prove a local convergence rate of $O(e^{-cn^{1/d}})$ for MS-GFEM for all these problems, improving upon the $O(e^{-cn^{1/(d+1)}})$ rate shown by Babuska and Lipton. Moreover, based on the abstract local approximation theory for MS-GFEM, we establish a unified framework for showing low-rank approximations to multiscale PDEs. This framework applies to the aforementioned problems, proving that the associated Green's functions admit an $O(|\log\epsilon|^{d})$-term separable approximation on well-separated domains with error $\epsilon>0$. Our analysis improves and generalizes the result in [M. Bebendorf and W. Hackbusch, Numerische Mathematik, 95 (2003), pp.~1-28] where an $O(|\log\epsilon|^{d+1})$-term separable approximation was proved for Poisson-type problems.

nlfem: a Flexible 2d FEM Python Code for Nonlocal Convection-Diffusion and Mechanics

NSPFEM2D: A Lightweight 2D Node-Based Smoothed Particle Finite Element Method Code for Modeling Large Deformation

A General Finite Element Method: Extension of Variational Analysis for Nonlinear Heat Conduction with Temperature-Dependent Properties and Boundary Conditions, and Its Implementation As Local Refinement.

A general finite element method

An adaptive non-symmetric finite volume and boundary element coupling method for a fluid mechanics interface problem

FEM on nonuniform meshes for nonlocal Laplacian: Semi-analytic Implementation in One Dimension

Non-intrusive implementation of a wide variety of Multiscale Finite Element Methods

A scalable domain decomposition method for FEM discretizations of nonlocal equations of integrable and fractional type

FEniCS Mechanics: A Package for Continuum Mechanics Simulations

DOLFIN: Automated Finite Element Computing

OBMeshfree: An Optimization-Based Meshfree Solver for Nonlocal Diffusion and Peridynamics Models

FDM/FEM for nonlinear convection–diffusion–reaction equations with Neumann boundary conditions—Convergence analysis for smooth and nonsmooth solutions

A localized Fourier collocation method for 2D and 3D elastic mechanics analysis: Theory and MATLAB code

An unstructured immersed finite element method for nonlinear solid mechanics

Integrating open source codes for solid modeling and automatic mesh generation in FEM applications

A unified framework for multiscale spectral generalized FEMs and low-rank approximations to multiscale PDEs

On the implementation of linear finite element method for nonlocal diffusion model over 2D domain

AutoFreeFem: Automatic code generation with FreeFEM++ and LaTex output for shape and topology optimization of non-linear multi-physics problems

MFEM: A modular finite element methods library

Local Discontinuous Galerkin for the Functional Renormalisation Group

A splitting, discontinuous Galerkin solver for the cell-by-cell electroneutral Nernst-Planck framework