Polyharmonic homogenization, rough polyharmonic splines and sparse super-localization

Abstract

We introduce a new variational method for the numerical homogenization of divergence form elliptic, parabolic and hyperbolic equations with arbitrary rough ($L^{\infty }$) coefficients. Our method does not rely on concepts of ergodicity or scale-separation but on compactness properties of the solution space and a new variational approach to homogenization. The approximation space is generated by an interpolation basis (over scattered points forming a mesh of resolution $H$) minimizing the $L^2$ norm of the source terms; its (pre-)computation involves minimizing $\mathcal{O}(H^{-d})$ quadratic (cell) problems on (super-)localized sub-domains of size $\mathcal{O}(H\ln (1/H))$. The resulting localized linear systems remain sparse and banded. The resulting interpolation basis functions are biharmonic for $d\le 3$, and polyharmonic for $d\ge 4$, for the operator $-\diiv(a\nabla \cdot)$ and can be seen as a generalization of polyharmonic splines to differential operators with arbitrary rough coefficients. The accuracy of the method ($\mathcal{O}(H)$ in energy norm and independent from aspect ratios of the mesh formed by the scattered points) is established via the introduction of a new class of higher-order Poincaré inequalities. The method bypasses (pre-)computations on the full domain and naturally generalizes to time dependent problems, it also provides a natural solution to the inverse problem of recovering the solution of a divergence form elliptic equation from a finite number of point measurements.

ESAIM: Mathematical Modelling and Numerical Analysis. Special issue (2013)

References in corpus (2)

• Inverse Problem for a Curved Quantum Guide
• Corrector theory for MsFEM and HMM in random media

Cited by in corpus (14)

• Constraint Energy Minimizing Generalized Multiscale Finite Element Method
• Kernel Flows: from learning kernels from data into the abyss
• Gamblets for opening the complexity-bottleneck of implicit schemes for hyperbolic and parabolic ODEs/PDEs with rough coefficients
• Fast Online Generalized Multiscale Finite Element Method using Constraint Energy Minimization
• Constraint Energy Minimizing Generalized Multiscale Finite Element Method in the Mixed Formulation
• Multiscale Partition of Unity
• Multiscale-Spectral GFEM and Optimal Oversampling
• A multiscale method for linear elasticity reducing Poisson locking
• Exponential Convergence for Multiscale Linear Elliptic PDEs via Adaptive Edge Basis Functions
• Adaptive Heterogeneous Multiscale Methods for immiscible two-phase flow in porous media
• A Super-Localized Generalized Finite Element Method
• Constraint Energy Minimizing Generalized Multiscale Finite Element Method for multi-continuum Richards equations
• Computational multiscale methods for nondivergence-form elliptic partial differential equations
• Uniform estimate of an iterative method for elliptic problems with rapidly oscillating coefficients