pith. sign in

arxiv: 2605.19540 · v1 · pith:ASMAPSSXnew · submitted 2026-05-19 · 🧮 math.AP · math.FA

Integral Equation Methods for Scattering by Multifractal Obstacles

Simon N. Chandler-Wilde (UOR) , Gabriel Claret (MICS) , David P. Hewett (UCL) , Anna Rozanova-Pierrat (MICS) , Siavash Sadeghi (UOR) This is my paper

Pith reviewed 2026-05-20 04:16 UTC · model grok-4.3

classification 🧮 math.AP math.FA
keywords acoustic scatteringintegral equationsfractal obstaclesmultifractal setsGalerkin methodsRadon measurestrace operatorsd-sets
0
0 comments X

The pith

Operator equations for acoustic scattering reduce to integral equations on multifractal supports when trace continuity and singular integral finiteness hold.

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

This paper extends earlier work on sound-soft scattering by compact obstacles with uniform fractal dimension to the multifractal case, where different parts of the scatterer have different fractal dimensions. It gives an interpretation of the governing operator equation in terms of trace spaces on the compact set Γ inside the obstacle and relates the unknown density to the normal derivative of the scattered field in many cases. The central result is that the operator equation becomes equivalent to an integral equation integrated against a Radon measure μ supported on Γ, provided the trace from H¹(ℝⁿ) into L²(Γ, μ) is continuous and certain canonical singular integrals with respect to μ remain finite. Galerkin methods using finite-element subspaces of L²(Γ, μ) converge if and only if smooth compactly supported functions vanishing near Γ are dense in the kernel of the trace operator. These conditions are shown to hold for finite unions of d-sets, including attractors of iterated function systems, with explicit convergence rates in the latter setting.

Core claim

The operator equation is equivalent to an integral equation on Γ whenever Γ is the support of a Radon measure μ such that the trace operator from H¹(ℝⁿ) to L²(Γ, μ) is continuous and certain canonical singular integrals with respect to μ are finite; Galerkin methods based on finite element subspaces of L²(Γ, μ) converge if and only if C∞₀(ℝⁿ∖Γ) is dense in the kernel of the trace operator. These results apply in particular if Γ is a finite union of d-sets with different values of d.

What carries the argument

Equivalence of the operator equation to integration against a Radon measure μ on Γ, which holds once the trace operator from H¹ to L²(Γ, μ) is continuous and the canonical singular integrals with respect to μ are finite.

If this is right

  • The scattering problem for a multifractal obstacle reduces to solving an integral equation whose integration is with respect to the measure μ on Γ.
  • Galerkin methods with finite-element subspaces of L²(Γ, μ) converge for any finite union of d-sets.
  • Explicit rates of convergence hold when each component d-set is the attractor of an iterated function system of contracting similarities.
  • In many cases the unknown density equals a suitable notion of the normal derivative of the scattered field on Γ.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same measure-theoretic conditions could be checked for other wave-scattering problems, such as electromagnetic or elastic scattering by multifractal objects.
  • Numerical codes could adapt mesh refinement to local fractal dimension once the measure μ is chosen.
  • Physical experiments with wave scattering from surfaces whose roughness changes across the domain would provide direct tests of the predicted convergence rates.
  • Links to potential theory on fractals might yield closed-form solutions for special iterated-function-system attractors.

Load-bearing premise

The Radon measures coming from finite unions of d-sets make the trace operator from H¹ continuous into L² on the support and keep the canonical singular integrals finite.

What would settle it

A concrete Radon measure supported on a multifractal set for which either the trace operator from H¹(ℝⁿ) to L²(Γ, μ) fails to be continuous or one of the canonical singular integrals diverges.

Figures

Figures reproduced from arXiv: 2605.19540 by Anna Rozanova-Pierrat (MICS), David P. Hewett (UCL), Gabriel Claret (MICS), Siavash Sadeghi (UOR), Simon N. Chandler-Wilde (UOR).

Figure 3
Figure 3. Figure 3: A commutative diagram showing the key function spaces and the mappings between them, [PITH_FULL_IMAGE:figures/full_fig_p016_3.png] view at source ↗
read the original abstract

Caetano et al. (Proc. R. Soc. A. 481:20230650, 2025) have proposed a formulation for sound-soft acoustic scattering by a compact scatterer O $\subset$ Rn, in which the scattered field is represented as an acoustic Newtonian potential whose density is the solution of an operator equation on a compact set $\Gamma$ $\subset$ O. In the case that $\Gamma$ is Ahlfors-David d-regular (a d-set), for some d $\in$ (n--2, n], they show, moreover, that the operator equation can be interpreted as an integral equation, the integration with respect to d-dimensional Hausdorff measure, and present a convergent Galerkin scheme for numerical computation. In this paper we make a substantial extension of these results so that they apply to more realistic fractal scatterers that are multifractal, in the sense that they have spatially varying fractal dimension. Firstly, we provide, inspired by Claret et al. (J. Math. Pures Appl. 212:103888, 2026), an interpretation of this operator equation as an equation between a trace space on $\Gamma$ and its dual, and, in many cases, relate the density to a notion of the normal derivative of the scattered field on $\Gamma$. Secondly, we show that the operator equation is equivalent to an integral equation on $\Gamma$ whenever $\Gamma$ is the support of a Radon measure $\mu$ such that: (i) the trace operator from H1(Rn) to L2($\Gamma$, $\mu$) is continuous and; (ii) certain canonical singular integrals with respect to $\mu$ are finite; and we characterise a large class of measures for which (i) and (ii) hold. Finally, we show that Galerkin methods based on finite element subspaces of L2($\Gamma$, $\mu$) are convergent if and only if, additionally, C$\infty$\_0 (Rn\$\Gamma$) is dense in the kernel of the trace operator. These results apply, in particular, if $\Gamma$ is a finite union of d-sets with different values of d. In the case that each d-set is the attractor of an iterated function system of contracting similarities, we establish rates of convergence for the Galerkin method.

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

0 major / 3 minor

Summary. The manuscript extends the integral equation formulation for sound-soft acoustic scattering by compact scatterers O ⊂ ℝⁿ to multifractal obstacles Γ that are finite unions of d-sets with varying dimensions d ∈ (n-2, n]. It recasts the operator equation in trace-space/dual form (inspired by Claret et al.), shows equivalence to an integral equation w.r.t. a Radon measure μ on Γ whenever the trace H¹(ℝⁿ) → L²(Γ, μ) is continuous and canonical singular integrals w.r.t. μ are finite, characterizes a large class of such measures (including the unions), and proves that Galerkin methods on finite-element subspaces of L²(Γ, μ) converge if and only if C∞₀(ℝⁿ ∖ Γ) is dense in the kernel of the trace operator. Convergence rates are established when each d-set is an IFS attractor.

Significance. If the central claims hold, the work meaningfully broadens the applicability of the Caetano et al. formulation from single d-regular sets to more realistic multifractal scatterers with spatially varying dimension. The explicit characterization of admissible measures, the if-and-only-if Galerkin criterion, and the control of cross terms via dimension gaps in the union case are clear strengths; the paper uses standard functional-analysis tools and prior literature without circularity or unverified embeddings. These results supply a rigorous foundation for numerical computation on complex fractals.

minor comments (3)
  1. Abstract, paragraph beginning 'Secondly': the phrasing 'certain canonical singular integrals with respect to μ are finite' is slightly vague; a forward reference to the precise definition (e.g., the principal-value integrals appearing in the single-layer or hypersingular operators) would improve clarity.
  2. Section on the union-of-d-sets case: while the dimension-gap argument for controlling cross terms is sketched, an explicit constant or inequality bounding the interaction between two d-sets with d1 ≠ d2 would make the verification easier to follow.
  3. Galerkin convergence theorem: the statement is clean, but the manuscript should note whether the finite-element subspaces are required to satisfy a uniform inverse inequality or quasi-uniformity assumption, as is standard in such analyses.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the careful and accurate summary of our manuscript, for highlighting its strengths in extending the Caetano et al. formulation to multifractal scatterers, and for the positive overall assessment. We appreciate the recommendation of minor revision and will incorporate any editorial or minor clarifications in the revised version.

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained via standard analysis

full rationale

The paper extends single d-set results from cited prior work to multifractal cases by recasting the operator equation in trace-space/dual form, then proving equivalence to an integral equation precisely when the trace operator H¹(ℝⁿ) → L²(Γ, μ) is continuous and canonical singular integrals w.r.t. μ are finite. It supplies an explicit characterization of such Radon measures and verifies the conditions hold for finite unions of d-sets (with cross terms controlled by dimension gap), deriving Galerkin convergence as an if-and-only-if statement from standard density arguments in the trace kernel. All steps rely on functional-analytic tools and independent verification rather than self-referential definitions, fitted parameters renamed as predictions, or load-bearing self-citations that reduce the central claims to unverified inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 3 axioms · 0 invented entities

The paper rests on standard results from functional analysis and geometric measure theory without introducing new free parameters or postulated entities.

axioms (3)
  • domain assumption Continuity of the trace operator from H¹(ℝⁿ) to L²(Γ, μ) for admissible measures μ
    Invoked as condition (i) for equivalence to an integral equation (abstract, 'Secondly, we show...').
  • domain assumption Finiteness of canonical singular integrals with respect to μ
    Invoked as condition (ii) for the integral-equation interpretation.
  • domain assumption Density of C∞₀(ℝⁿ∖Γ) in the kernel of the trace operator implies Galerkin convergence
    Stated as the additional condition for convergence of the numerical scheme.

pith-pipeline@v0.9.0 · 6004 in / 1738 out tokens · 61973 ms · 2026-05-20T04:16:03.876508+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]

    D. R. Adams and L. I. Hedberg. Function Spaces and Potential Theory . Springer, 2012

    work page 2012

  2. [2]

    S. Agmon. Spectral properties of Schr¨ odinger operators and scattering theory. Ann. Sc. norm. super. Pisa - Cl. Sci. (4) , 2:151–218, 1975

    work page 1975

  3. [3]

    M. Biegert. On traces of Sobolev functions on the boundary of extension domains. Proc. Am. Math. Soc., 137:4169–4176, 2009

    work page 2009

  4. [4]

    Caetano, D

    A. Caetano, D. P. Hewett, and A. Moiola. Density results for Sobolev, Besov and Triebel- Lizorkin spaces on rough sets. J. Funct. Anal., 281:109019, 2021

    work page 2021

  5. [5]

    A. M. Caetano, S. N. Chandler-Wilde, X. Claeys, A. Gibbs, D. P. Hewett, and A. Moiola. Integral equation methods for acoustic scattering by fractals. Proc. R. Soc. A , 481:20230650, 2025

    work page 2025

  6. [6]

    A. M. Caetano, S. N. Chandler-Wilde, A. Gibbs, D. P. Hewett, and A. Moiola. A Hausdorff- measure boundary element method for acoustic scattering by fractal screens. Numer. Math. , 156:463–532, 2024

    work page 2024

  7. [7]

    A Hausdorff-measure boundary element method for acoustics scattering by fractal screens

    A M Caetano, S N Chandler-Wilde, A Gibbs, D P Hewett, and A Moiola. A Hausdorff-measure boundary element method for acoustics scattering by fractal screens. Numer. Math., 156:463– 532, 2024

    work page 2024

  8. [8]

    A. M. Caetano, S. N. Chandler-Wilde, and D. P. Hewett. Properties of IFS attractors with non- empty interiors and associated function spaces and scattering problems, 2025. arXiv preprint arXiv:2511.15213

    work page arXiv 2025

  9. [9]

    Carvalho and A

    A. Carvalho and A. Caetano. On the Hausdorff dimension of continuous functions belonging to H¨ older and Besov spaces on fractald-sets. J. Fourier Anal. Appl. , 18:386–409, 2012

    work page 2012

  10. [10]

    S. N. Chandler-Wilde, D. P. Hewett, and A. Moiola. Sobolev spaces on non-Lipschitz subsets of Rn with application to boundary integral equations on fractal screens. Integr. Equ. Oper. Theory, 87:179–224, 2017

    work page 2017

  11. [11]

    S. N. Chandler-Wilde, D. P. Hewett, A. Moiola, and J. Besson. Boundary element methods for acoustic scattering by fractal screens. Numer. Math., 147:785–837, 2021

    work page 2021

  12. [12]

    S. N. Chandler-Wilde and P. Monk. Wave-number-explicit bounds in time-harmonic scattering. SIAM J. Math. Anal. , 39:1428–1455, 2008

    work page 2008

  13. [13]

    S. N. Chandler-Wilde and S. Sadeghi. Integral equation methods for scattering by general compact obstacles: wavenumber-explicit estimates, 2026. arXiv preprint arXiv:2601.19456

    work page arXiv 2026

  14. [14]

    Claret, M

    G. Claret, M. Hinz, A. Rozanova-Pierrat, and A. Teplyaev. Layer potential operators for transmission problems on extension domains. J. Math. Pures Appl. , 212:103888, 2026

    work page 2026

  15. [15]

    Colton and R

    D. Colton and R. Kress. Integral Equation Methods in Scattering Theory . John Wiley, 1983

    work page 1983

  16. [16]

    L. C. Evans and R. E. Gariepy. Measure Theory and Fine Properties of Functions. CRC Press, 2015. 35

    work page 2015

  17. [17]

    Falconer

    K. Falconer. Fractal Geometry: Mathematical Foundations and Applications, 3rd Ed. Wiley, 2014

    work page 2014

  18. [18]

    Gibbs, D

    A. Gibbs, D. P. Hewett, and B. Major. Numerical evaluation of singular integrals on non-disjoint self-similar fractal sets. Numer. Alg., 97:311–343, 2024

    work page 2024

  19. [19]

    Gibbs, D

    A. Gibbs, D. P. Hewett, and A. Moiola. Numerical quadrature for singular integrals on fractals. Numer. Alg., 92:2071–2124, 2023

    work page 2071

  20. [20]

    Haj lasz, P

    P. Haj lasz, P. Koskela, and H. Tuominen. Sobolev embeddings, extensions and measure density condition. J. Funct. Anal., 254(5):1217–1234, 2008

    work page 2008

  21. [21]

    D. P. Hewett and A. Moiola. On the maximal Sobolev regularity of distributions supported by subsets of Euclidean space. Appl. Anal., 15:731–770, 2017

    work page 2017

  22. [22]

    M. Hinz, S. N. Chandler-Wilde, and D. P. Hewett. Kernels of trace operators via fine continuity. arXiv preprint arXiv:2507.04536v1 , 2025

    work page arXiv 2025

  23. [23]

    M. Hinz, A. Rozanova-Pierrat, and A. Teplyaev. Boundary value problems on non-Lipschitz uniform domains: stability, compactness and the existence of optimal shapes. Asymptotic Anal., 134(1-2):25–61, 2023

    work page 2023

  24. [24]

    P. Joly, M. Kachanovska, and Z. Moitier. High-order numerical integration on self-affine sets,

    work page

  25. [25]

    arXiv preprint arXiv:2410.00637

    work page arXiv

  26. [26]

    A. Jonsson. The trace of potentials on general sets. Ark. Mat., 17:1–18, 1979

    work page 1979

  27. [27]

    Jonsson and H

    A. Jonsson and H. Wallin. Function Spaces on Subsets of Rn, Math. Reports., vol. 2 . Harwood Acad. Publ., 1984

    work page 1984

  28. [28]

    Kr´ al.Integral Operators in Potential Theory

    J. Kr´ al.Integral Operators in Potential Theory . Springer, 1980

    work page 1980

  29. [29]

    Levitin, D

    M. Levitin, D. Mangoubi, and I. Polterovich. Topics in Spectral Geometry. American Mathe- matical Society, 2023

    work page 2023

  30. [30]

    P. Mattila. Geometry of Sets and Measures in Euclidean Spaces: Fractals and Rectifiability . CUP, 1995

    work page 1995

  31. [31]

    W. C. H. McLean. Strongly Elliptic Systems and Boundary Integral Equations . Cambridge University Press, 2000

    work page 2000

  32. [32]

    Mitrea and M

    M. Mitrea and M. Taylor. Boundary layer methods for Lipschitz domains in Riemannian manifolds. J. Funct. Anal., 163:181–251, 1999

    work page 1999

  33. [33]

    Digital library of mathematical functions, 2023

    NIST. Digital library of mathematical functions, 2023. See https://dlmf.nist.gov

    work page 2023

  34. [34]

    L. Riddle. Classic iterated function systems. See https://larryriddle.agnesscott.org/ ifs/ifs.htm (accessed 18 July 2023)

    work page 2023

  35. [35]

    Rozanova-Pierrat

    A. Rozanova-Pierrat. Generalization of Rellich–Kondrachov theorem and trace compactness for fractal boundaries. In M. R. Lancia and A. Rozanova-Pierrat, editors, Fractals in Engineering: Theoretical Aspects and Numerical Approximations, pages 155–173. Springer, 2021

    work page 2021

  36. [36]

    Rozanova-Pierrat

    A. Rozanova-Pierrat. On analysis of problems of mathematical physics with non-Lipschitz boundaries, 2025. arXiv preprint arXiv:2509.02176. 36

    work page arXiv 2025

  37. [37]

    W. Rudin. Real and Complex Analysis, 3rd ed. McGraw-Hill, New York, 1986

    work page 1986

  38. [38]

    W. Rudin. Functional Analysis, 2nd ed. McGraw-Hill, New York, 1991

    work page 1991

  39. [39]

    S. A. Sauter and C. Schwab. Boundary Element Methods . Springer-Verlag, Berlin, 2011

    work page 2011

  40. [40]

    H. A. Schenck. Improved integral formulation for acoustic radiation problems. J. Acoust. Soc. Am., 44:41–58, 1968

    work page 1968

  41. [41]

    H. Triebel. Fractals and Spectra. Birkh¨ auser, 1997

    work page 1997

  42. [42]

    I. Ja. ˇSne˘ ıberg. Spectral properties of linear operators in interpolation families of Banach spaces. Mat. Issled., 9:214–229, 1974

    work page 1974

  43. [43]

    G. N. Watson. Theory of Bessel Functions, 2nd Ed. Cambridge University Press, 1966

    work page 1966

  44. [44]

    T. W. Wu and A. F. Seybert. A weighted residual formulation for the CHIEF method in acoustics. J. Acoust. Soc. Am., 90:1608–1614, 1991. 37

    work page 1991