pith. sign in

arxiv: 2605.22159 · v1 · pith:3B5EQKW6new · submitted 2026-05-21 · 🧮 math.NA · cs.NA

BEM for variable coefficient second-order problems

Benedikt Gr\"a{\ss}le , Stefan A. Sauter This is my paper

Pith reviewed 2026-05-22 03:59 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords boundary element methodvariable coefficientselliptic operatorsGalerkin discretizationdata-sparse compressionfundamental solutionssecond-order problems
0
0 comments X

The pith

A boundary element method approximates the boundary operator from a volume discretization to handle variable-coefficient elliptic problems without explicit fundamental solutions.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper introduces a boundary element method for strongly elliptic operators with variable coefficients. It builds a computable approximation to the boundary integral operator through a one-time Galerkin discretization of the underlying differential operator, such as by finite elements. The resulting method keeps the dimension reduction of boundary integral approaches and supports standard data-sparse compression. A reader would care because this removes the classical requirement for explicit fundamental solutions that often do not exist or are difficult to find for variable coefficients.

Core claim

The central claim is that a computable approximation of the boundary operator can be constructed from a Galerkin discretisation of the underlying elliptic differential operator in a one-time preprocessing step. This yields an algebraic formulation that retains the dimension reduction of boundary integral methods and remains compatible with data-sparse matrix compression while extending quasi-optimal BEM discretisations to variable-coefficient problems.

What carries the argument

The central mechanism is the Galerkin discretization of the volume elliptic operator that produces the approximation to the boundary integral operator.

Load-bearing premise

The Galerkin discretization of the volume operator must produce an approximation to the boundary operator accurate enough to retain quasi-optimal convergence rates while staying compatible with data-sparse compression.

What would settle it

Numerical tests in which the convergence rate falls below the expected quasi-optimal order when the volume discretization mesh is coarsened relative to the boundary mesh would show the approximation is not accurate enough.

Figures

Figures reproduced from arXiv: 2605.22159 by Benedikt Gr\"a{\ss}le, Stefan A. Sauter.

Figure 1
Figure 1. Figure 1: Isomorphism between Bs , Xs , and Ws ⊂ V for all 0 ≤ s ≤ σreg. Theorems 3.1 and 3.2 imply that S : B → W and V : B → X are isomorphisms between the spaces B, X, and the set of solutions to (3.3) given by W := {v ∈ V : Lv = 0 in V ′ 0 in the sense of (3.3a)} ⊂ V. Associate each of these spaces Y ∈ {B, X, W} with a scale (Y s )s≥0 of densely embedded Hilbert spaces indexed by some parameter 0 < s such that Y… view at source ↗
read the original abstract

A novel boundary element method (BEM) removes the classical dependence on explicit fundamental solutions and extends quasi-optimal BEM discretisations to strongly elliptic operators with variable coefficients. The approach constructs a computable approximation of the boundary operator from a Galerkin discretisation of the underlying elliptic differential operator in a one-time preprocessing step, for instance by conforming finite elements. The resulting algebraic formulation retains the dimension reduction intrinsic to boundary integral methods and is compatible with standard data-sparse matrix compression techniques.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper proposes a novel boundary element method for strongly elliptic second-order problems with variable coefficients. It constructs a computable approximation to the boundary integral operator via a one-time conforming Galerkin discretization (e.g., finite elements) of the underlying volume elliptic operator, thereby eliminating the need for an explicit fundamental solution while preserving the dimension reduction of BEM and compatibility with data-sparse compression.

Significance. If the operator approximation error can be rigorously controlled to preserve quasi-optimality, the method would meaningfully extend BEM to a wider class of variable-coefficient problems where classical fundamental solutions are unavailable. The one-time preprocessing step and compression compatibility are practical strengths that could enable efficient high-dimensional simulations.

major comments (2)
  1. [Abstract and method construction] The central claim that quasi-optimal convergence rates are retained depends on the volume Galerkin approximation producing a boundary operator error small enough not to pollute the Céa estimate or inf-sup constant. No explicit a priori bound relating volume mesh size, coefficient variation, and boundary discretization error is derived or stated, leaving the transfer of accuracy from volume to boundary unquantified.
  2. [Abstract and method construction] Without a consistency analysis or numerical verification showing that the preprocessing error remains of higher order than the boundary mesh size h, the assertion of compatibility with standard BEM quasi-optimality cannot be verified from the given material.
minor comments (1)
  1. [Algebraic formulation] Clarify how the approximated boundary operator is assembled into the final algebraic system and whether any additional stabilization is required for stability under variable coefficients.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive report and positive assessment of the method's potential significance. We address the two major comments point by point below and agree that additional analysis is warranted to fully substantiate the claims.

read point-by-point responses

  1. Referee: [Abstract and method construction] The central claim that quasi-optimal convergence rates are retained depends on the volume Galerkin approximation producing a boundary operator error small enough not to pollute the Céa estimate or inf-sup constant. No explicit a priori bound relating volume mesh size, coefficient variation, and boundary discretization error is derived or stated, leaving the transfer of accuracy from volume to boundary unquantified.

    Authors: We agree that an explicit a priori bound would make the transfer of accuracy from the volume preprocessing step to the boundary operator more transparent. The manuscript establishes well-posedness and quasi-optimality of the resulting BEM under the assumption that the volume Galerkin error is controlled sufficiently, relying on standard Céa-type estimates for the volume problem and continuity of the boundary trace operator. To address the concern directly, we will add a dedicated subsection deriving a concrete bound on the boundary operator approximation error in terms of the volume mesh size and the variation of the coefficients, assuming standard elliptic regularity. revision: yes

  2. Referee: [Abstract and method construction] Without a consistency analysis or numerical verification showing that the preprocessing error remains of higher order than the boundary mesh size h, the assertion of compatibility with standard BEM quasi-optimality cannot be verified from the given material.

    Authors: The manuscript presents the algebraic construction and its compatibility with data-sparse compression, with consistency following from the conforming Galerkin property of the volume discretization. We acknowledge that an explicit consistency argument relating the preprocessing error to the boundary mesh size h is not detailed in the current text. In the revised version we will include a short consistency analysis showing that, by selecting the volume mesh size sufficiently small relative to h (e.g., O(h^{1+δ}) for suitable δ), the preprocessing error remains of strictly higher order and does not degrade the quasi-optimal rate. We will also add a brief numerical example confirming the observed convergence rates. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation remains self-contained

full rationale

The abstract and description present a construction in which a Galerkin discretization of the volume elliptic operator is used once to approximate the boundary integral operator. No equations are supplied that equate the final quasi-optimal rates or the algebraic BEM formulation to a fitted parameter or to a self-referential definition of the same operator. The central claim therefore rests on an independent transfer of accuracy from the volume discretization to the boundary operator, rather than on any of the enumerated circular patterns. The skeptic concern addresses the size of the approximation error relative to the boundary mesh, which is a question of quantitative analysis and not a reduction of the result to its own inputs by construction. The method is compatible with external BEM theory and data-sparse techniques without invoking self-citations or ansatzes that close a loop.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The approach rests on the standard assumption that the differential operator is strongly elliptic and that a conforming Galerkin discretization exists and can be computed once. No free parameters or new entities are introduced in the abstract.

axioms (1)
  • domain assumption The differential operator is strongly elliptic.
    Explicitly stated in the abstract as the class of operators for which the method applies.

pith-pipeline@v0.9.0 · 5599 in / 1188 out tokens · 39100 ms · 2026-05-22T03:59:16.532608+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

44 extracted references · 44 canonical work pages

  1. [1]

    M. S. Agranovich. Spectral problems for second-order strongly elliptic systems in smooth and non-smooth domains. Russ. Math. Surv. , 57(5):847--920, October 2002

    work page 2002

  2. [2]

    Babu s ka

    I. Babu s ka. The finite element method for elliptic equations with discontinuous coefficients. Computing , 5(3):207--213, September 1970

    work page 1970

  3. [3]

    A. Barton. Layer potentials for general linear elliptic systems. Electron J Differ. Equ. , Paper No. 309:23, 2017

    work page 2017

  4. [4]

    B \'e cache, A.-S

    \'E . B \'e cache, A.-S. Bonnet-Ben Dhia , S. Fliss, and A. Tonnoir. The half-space matching method for elastodynamic scattering problems in unbounded domains. Journal of Computational Physics , 490:112320, October 2023

    work page 2023

  5. [5]

    Bebendorf

    M. Bebendorf. Efficient inversion of the Galerkin matrix of general second-order elliptic operators with nonsmooth coefficients. Math. Comp. , 74(251):1179--1200, September 2004

    work page 2004

  6. [6]

    Bernardi

    C. Bernardi. Optimal finite-element interpolation on curved domains. SIAM J. Numer. Anal. , 26(5):1212--1240, October 1989

    work page 1989

  7. [7]

    Bebendorf and W

    M. Bebendorf and W. Hackbusch. Existence of H -matrix approximants to the inverse FE-matrix of elliptic operators with L ^ -coefficients. Numerische Mathematik , 95(1):1--28, July 2003

    work page 2003

  8. [8]

    O. P. Bruno and T. Yin. A windowed Green function method for elastic scattering problems on a half-space. Computer Methods in Applied Mechanics and Engineering , 376:113651, April 2021

    work page 2021

  9. [9]

    S. Börm. Approximation of solution operators of elliptic partial differential equations by H and H ^2 -matrices. Numer. Math. , 115(2):165--193, April 2010

    work page 2010

  10. [10]

    M. H. Cho and W. Cai. A parallel fast algorithm for computing the Helmholtz integral operator in 3- D layered media. Journal of Computational Physics , 231(17):5910--5925, July 2012

    work page 2012

  11. [11]

    Chandler-Wilde and D

    S. Chandler-Wilde and D. Hothersall. A uniformly valid far field asymptotic expansion of the green function for two-dimensional propagation above a homogeneous impedance plane. Journal of Sound and Vibration , 182(5):665--675, May 1995

    work page 1995

  12. [12]

    P. G. Ciarlet. The Finite Element Method for Elliptic Problems . Classics in Applied Mathematics . Society for Industrial and Applied Mathematics , January 2002

    work page 2002

  13. [13]

    Costabel

    M. Costabel. Boundary integral operators on Lipschitz domains: Elementary results. SIAM J. Math. Anal. , 19(3):613--626, May 1988

    work page 1988

  14. [14]

    Dahmen, B

    W. Dahmen, B. Faermann, I. G. Graham, W. Hackbusch, and S. A. Sauter. Inverse inequalities on non-quasi-uniform meshes and application to the mortar element method. Math. Comp. , 73(247):1107--1138, 2004

    work page 2004

  15. [15]

    Dur \'a n, R

    M. Dur \'a n, R. Hein, and J.-C. N \'e d \'e lec. Computing numerically the Green 's function of the half-plane Helmholtz operator with impedance boundary conditions. Numer. Math. , 107(2):295--314, August 2007

    work page 2007

  16. [16]

    Dolzmann and S

    G. Dolzmann and S. M \"u ller. Estimates for Green 's matrices of elliptic systems by L ^p theory. Manuscripta Math , 88(1):261--273, December 1995

    work page 1995

  17. [17]

    K. Feng. Finite element method and natural boundary reduction. In Proc. Int . Congr . Math . , pages 1439--1453, 1983

    work page 1983

  18. [18]

    Florian, R

    F. Florian, R. Hiptmair, and S. A. Sauter. Skeleton integral equations for acoustic transmission problems with varying coefficients. SIAM J. Math. Anal. , 56(5):6232--6267, October 2024

    work page 2024

  19. [19]

    Faustmann, J

    M. Faustmann, J. M. Melenk, and D. Praetorius. H -matrix approximability of the inverses of FEM matrices. Numer. Math. , 131(4):615--642, December 2015

    work page 2015

  20. [20]

    Gr \"a le, R

    B. Gr \"a le, R. Hiptmair, and S. A. Sauter. Stable skeleton integral equations for general coefficient Helmholtz transmission problems. arXiv:2507.00991 , 2025

    work page arXiv 2025

  21. [21]

    Gr \"a le

    B. Gr \"a le. Optimal trace norms for Helmholtz problems. arXiv:2506.11944 , 2025

    work page arXiv 2025

  22. [22]

    Grisvard

    P. Grisvard. Elliptic Problems in Nonsmooth Domains . Society for Industrial and Applied Mathematics , 1985

    work page 1985

  23. [23]

    Grisvard

    P. Grisvard. Singularities in Boundary Value Problems . Number 22 in Research notes in applied mathematics. Masson, Paris, 1992

    work page 1992

  24. [24]

    Gr \"a le and S

    B. Gr \"a le and S. A. Sauter. Dirichlet-to- Neumann operator for the Helmholtz problem with general wavenumbers on the n -sphere. arXiv:2503.18837 , 2025

    work page arXiv 2025

  25. [25]

    Gr \"u ter and K.-O

    M. Gr \"u ter and K.-O. Widman. The Green function for uniformly elliptic equations. Manuscripta Math , 37(3):303--342, October 1982

    work page 1982

  26. [26]

    Hackbusch

    W. Hackbusch. Hierarchical Matrices : Algorithms and Analysis . Springer Series in Computational Mathematics . Springer Berlin Heidelberg, Berlin, Heidelberg, 2015

    work page 2015

  27. [27]

    Hofmann and S

    S. Hofmann and S. Kim. The Green function estimates for strongly elliptic systems of second order. manuscripta math. , 124(2):139--172, September 2007

    work page 2007

  28. [28]

    Haller-Dintelmann , H

    R. Haller-Dintelmann , H. Meinlschmidt, and W. Wollner. Higher regularity for solutions to elliptic systems in divergence form subject to mixed boundary conditions. Annali di Matematica , 198(4):1227--1241, August 2019

    work page 2019

  29. [29]

    Jochmann

    F. Jochmann. An H^s -regularity result for the gradient of solutions to elliptic equations with mixed boundary conditions. Journal of Mathematical Analysis and Applications , 238(2):429--450, October 1999

    work page 1999

  30. [30]

    J. B. Keller and D. Givoli. Exact non-reflecting boundary conditions. Journal of Computational Physics , 82(1):172--192, May 1989

    work page 1989

  31. [31]

    Kim and G

    S. Kim and G. Sakellaris. Green's function for second order elliptic equations with singular lower order coefficients. Communications in Partial Differential Equations , 44(3):228--270, March 2019

    work page 2019

  32. [32]

    W. Lu, Y. Y. Lu, and J. Qian. Perfectly matched layer boundary integral equation method for wave scattering in a layered medium. SIAM J. Appl. Math. , 78(1):246--265, January 2018

    work page 2018

  33. [33]

    Lions and E

    J.-L. Lions and E. Magenes. Non-homogeneous boundary value problems and applications. Vol . I . Springer-Verlag, New York-Heidelberg, 1972

    work page 1972

  34. [34]

    C. Lin, J. M. Melenk, and S. Sauter. An explicit factorization of the Green 's function for an acoustic half-space problem with impedance boundary conditions into an oscillatory exponential and a slowly varying function. Math Methods in App Sciences , 48(13):12807--12812, September 2025

    work page 2025

  35. [35]

    J. Li, J. M. Melenk, B. Wohlmuth, and J. Zou. Optimal a priori estimates for higher order finite elements for elliptic interface problems. Applied Numerical Mathematics , 60(1-2):19--37, January 2010

    work page 2010

  36. [36]

    Littman, G

    W. Littman, G. Stampacchia, and H. F. Weinberger. Regular points for elliptic equations with discontinuous coefficients. Ann. Della Scuola Norm. Super. Pisa - Cl. Sci. , 17(1-2):43--77, 1963

    work page 1963

  37. [37]

    L. Li, J. Yang, B. Zhang, and H. Zhang. Uniform far-field asymptotics of the two-layered Green function in two dimensions and application to wave scattering in a two-layered medium. SIAM J. Math. Anal. , 56(3):4143--4184, June 2024

    work page 2024

  38. [38]

    W. McLean. Strongly elliptic systems and boundary integral equations . Cambridge University Press, Cambridge, 1. publ edition, 2000

    work page 2000

  39. [39]

    R. C. Maccamy and S. P. Marin. A finite element method for exterior interface problems. International Journal of Mathematics and Mathematical Sciences , 3(2):311--350, January 1980

    work page 1980

  40. [40]

    Mayboroda and B

    S. Mayboroda and B. Poggi. Exponential decay estimates for fundamental solutions of Schr \"o dinger-type operators. Trans. Amer. Math. Soc. , 372(6):4313--4357, April 2019

    work page 2019

  41. [41]

    J. M. Melenk and S. A. Sauter. Convergence analysis for finite element discretizations of the Helmholtz equation with Dirichlet-to-Neumann boundary conditions. Math. Comp. , 79(272):1871--1914, April 2010

    work page 1914

  42. [42]

    M. A. Olshanskii and A. Reusken. Trace Finite Element Methods for PDEs on Surfaces . In S. P. A. Bordas, E. Burman, M. G. Larson, and M. A. Olshanskii, editors, Geometrically Unfitted Finite Element Methods and Applications , volume 121, pages 211--258. Springer International Publishing, Cham, 2017

    work page 2017

  43. [43]

    Savar \'e

    G. Savar \'e . Regularity results for elliptic equations in Lipschitz domains. Journal of Functional Analysis , 152(1):176--201, January 1998

    work page 1998

  44. [44]

    S. A. Sauter and C. Schwab. Boundary element methods . Number 39 in Springer Series in Computational Mathematics . Springer Berlin Heidelberg, Berlin, Heidelberg, 2011

    work page 2011