High order recovery of geometric interfaces from cell-average data

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Cohen, Albert, Mula, Olga, Somacal, Agustín
Natura: Preprint
Pubblicazione: 2024
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866929552746348544
author Cohen, Albert
Mula, Olga
Somacal, Agustín
author_facet Cohen, Albert
Mula, Olga
Somacal, Agustín
contents We consider the problem of recovering characteristic functions $u:=χ_Ω$ from cell-average data on a coarse grid, and where $Ω$ is a compact set of $\mathbb{R}^d$. This task arises in very different contexts such as image processing, inverse problems, and the accurate treatment of interfaces in finite volume schemes. While linear recovery methods are known to perform poorly, nonlinear strategies based on local reconstructions of the jump interface $Γ:=\partialΩ$ by geometrically simpler interfaces may offer significant improvements. We study two main families of local reconstruction schemes, the first one based on nonlinear least-squares fitting, the second one based on the explicit computation of a polynomial-shaped curve fitting the data, which yields simpler numerical computations and high order geometric fitting. For each of them, we derive a general theoretical framework which allows us to control the recovery error by the error of best approximation up to a fixed multiplicative constant. Numerical tests in 2d illustrate the expected approximation order of these strategies. Several extensions are discussed, in particular the treatment of piecewise smooth interfaces with corners.
format Preprint
id arxiv_https___arxiv_org_abs_2402_00946
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle High order recovery of geometric interfaces from cell-average data
Cohen, Albert
Mula, Olga
Somacal, Agustín
Numerical Analysis
We consider the problem of recovering characteristic functions $u:=χ_Ω$ from cell-average data on a coarse grid, and where $Ω$ is a compact set of $\mathbb{R}^d$. This task arises in very different contexts such as image processing, inverse problems, and the accurate treatment of interfaces in finite volume schemes. While linear recovery methods are known to perform poorly, nonlinear strategies based on local reconstructions of the jump interface $Γ:=\partialΩ$ by geometrically simpler interfaces may offer significant improvements. We study two main families of local reconstruction schemes, the first one based on nonlinear least-squares fitting, the second one based on the explicit computation of a polynomial-shaped curve fitting the data, which yields simpler numerical computations and high order geometric fitting. For each of them, we derive a general theoretical framework which allows us to control the recovery error by the error of best approximation up to a fixed multiplicative constant. Numerical tests in 2d illustrate the expected approximation order of these strategies. Several extensions are discussed, in particular the treatment of piecewise smooth interfaces with corners.
title High order recovery of geometric interfaces from cell-average data
topic Numerical Analysis
url https://arxiv.org/abs/2402.00946