Papers › Reduced order modeling for elliptic problems with high contrast diffusion coefficients

Reduced order modeling for elliptic problems with high contrast diffusion coefficients

21 Apr 2023arXiv:2304.10971links table onlyarchive 2025-07-28

Albert Cohen, Matthieu Dolbeault, Agustin Somacal, Wolfgang Dahmen

The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.

We consider the parametric elliptic PDE -div (a(y)∇u)=f on a spatial domain Ω, with a(y) a scalar piecewise constant diffusion coefficient taking any positive values y=(y₁, …, y_d)∈]0,∞[ᵈ on fixed subdomains Ω₁,…,Ω_d. This problem is not uniformly elliptic as the contrast κ(y)=maxyⱼ/minyⱼ can be arbitrarily high, contrarily to the Uniform Ellipticity Assumption (UEA) that is commonly made on parametric elliptic PDEs. Based on local polynomial approximations in the y variable, we construct local and global reduced model spaces Vₙ of moderate dimension n that approximate uniformly well all solutions u(y). Since the solution u(y) blows as y→0, the solution manifold is not a compact set and does not have finite n-width. Therefore, our results for approximation by such spaces are formulated in terms of relative H¹₀-projection error, that is, after normalization by u(y)_(H¹₀). We prove that this relative error decays exponentially with n, yet exhibiting the curse of dimensionality as the number d of subdomains grows. We also show similar rates for the Galerkin projection despite the fact that high contrast is well-known to deteriorate the multiplicative constant when applying Cea's lemma. We finally establish uniform estimates in relative error for the state estimation and parameter estimation inverse problems, when y is unknown and a limited number of linear measurements ℓᵢ(u) are observed. A key ingredient in our construction and analysis is the study of the convergence of u(y) to limit solutions when some of the parameters yⱼ tend to infinity.

PaperPDFCode

Code

agussomacal/romhighcontrast officialmentioned in paper report

Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.

Code Syntology ran Syntology

Not run by Syntology. Nothing on this page verifies that the listed code works.

Results from the paper archive 2025-07-28

No leaderboard rows for this paper in the archive.

Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections