English

Finite element approximation of the Hardy constant

Numerical Analysis 2024-02-06 v2 Numerical Analysis Analysis of PDEs

Abstract

We consider finite element approximations to the optimal constant for the Hardy inequality with exponent p=2p=2 in bounded domains of dimension n=1n=1 or n3n \geq 3. For finite element spaces of piecewise linear and continuous functions on a mesh of size hh, we prove that the approximate Hardy constant converges to the optimal Hardy constant at a rate proportional to 1/logh21/| \log h |^2. This result holds in dimension n=1n=1, in any dimension n3n \geq 3 if the domain is the unit ball and the finite element discretization exploits the rotational symmetry of the problem, and in dimension n=3n=3 for general finite element discretizations of the unit ball. In the first two cases, our estimates show excellent quantitative agreement with values of the discrete Hardy constant obtained computationally.

Keywords

Cite

@article{arxiv.2308.01580,
  title  = {Finite element approximation of the Hardy constant},
  author = {Francesco Della Pietra and Giovanni Fantuzzi and Liviu I. Ignat and Alba Lia Masiello and Gloria Paoli and Enrique Zuazua},
  journal= {arXiv preprint arXiv:2308.01580},
  year   = {2024}
}

Comments

Review: Significantly improved estimates compared to the original version (23 pages, 6 figures)

R2 v1 2026-06-28T11:47:05.280Z