Skip to content
DraconiticResearch
Search Sources API
197,688 studies31 sourcesUpdated 14 min ago
PreprintarXiv licencearXiv

Quadrature by Matched Asymptotic Expansion (QBMAX): Accelerating Evaluation of Layer Potentials with Boundary Layers

Abstract

Evaluating layer potentials with rapidly decaying kernels becomes computationally challenging in the presence of boundary layers, as occurs for Helmholtz layer potentials with large imaginary wavenumbers. We introduce Quadrature by Matched Asymptotic Expansion (QBMAX), a high-order method that extends the Quadrature by Expansion (QBX) framework by incorporating the kernel's asymptotic decay behavior directly into the Taylor expansion of the potential. Numerical results in both 2D and 3D show improvements of up to four digits over QBX for the tested large-parameter cases, while the methods can perform similarly for smaller parameters or more strongly curved geometries. A symbolic weighted operation-count model indicates comparable costs for forming the two expansion kernels. The Taylor expansions have the same formal order, while QBMAX reduces the truncation-error constants in the flat-boundary model and generally exhibits less convergence degradation in the reported experiments.

Related

  1. Time-asymptotic expansion and fluid approximation of the 1D Boltzmann equation with hard potentials
    Hui-Tzu Chang, Haitao Wang, Kung-Chien Wu
    arXivPreprintarXiv licenceFull textAnalysis & PDEsarXiv
  2. AMPERe: an exactly divergence-free quadrature-by-expansion evaluator for near-surface Biot--Savart and virtual-casing fields in fusion geometries
    H. Betar, Y. Sarazin, V. Grandgirard +1
    arXivPreprintarXiv licenceFull textMathematical physicsComputational physicsarXiv
  3. Numerical evaluation of one-loop integrals using quadrature rules
    Jürgen Dölz, Claude Duhr, Bartul Kovačić +1
    arXivPreprintarXiv licenceFull textNumerical analysis & scientific computingHigh-energy & particle physicsarXiv
  4. Asymptotic expansion of smooth functions in deterministic and iid Haar unitary matrices, and application to tensor products of matrices
    Félix Parraud
    The Annals of ProbabilityPeer-reviewedLicence not statedFull text10 citedStatisticsOpenAlex
  5. Performance Evaluation of an Adaptive Quadrature and a Double Exponential Formula Using Arbitrary-Precision Floating-Point Arithmetic
    Tomonori Kouya
    arXivPreprintCC BY 4.0Full textNumerical analysis & scientific computingarXiv
BibTeX
@misc{lin2026quadrature,
  title = {{Quadrature by Matched Asymptotic Expansion (QBMAX): Accelerating Evaluation of Layer Potentials with Boundary Layers}},
  author = {Chaoqi Lin and Xiaoyu Wei and Andreas Klöckner},
  howpublished = {arXiv},
  year = {2026},
  month = sep,
  doi = {10.48550/arxiv.2609.10857},
  eprint = {2609.10857},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2609.10857},
}
RIS
TY  - UNPB
TI  - Quadrature by Matched Asymptotic Expansion (QBMAX): Accelerating Evaluation of Layer Potentials with Boundary Layers
AU  - Chaoqi Lin
AU  - Xiaoyu Wei
AU  - Andreas Klöckner
PY  - 2026
DA  - 2026/09/12
PB  - arXiv
DO  - 10.48550/arxiv.2609.10857
UR  - https://arxiv.org/abs/2609.10857
AB  - Evaluating layer potentials with rapidly decaying kernels becomes computationally challenging in the presence of boundary layers, as occurs for Helmholtz layer potentials with large imaginary wavenumbers. We introduce Quadrature by Matched Asymptotic Expansion (QBMAX), a high-order method that extends the Quadrature by Expansion (QBX) framework by incorporating the kernel's asymptotic decay behavior directly into the Taylor expansion of the potential. Numerical results in both 2D and 3D show improvements of up to four digits over QBX for the tested large-parameter cases, while the methods can perform similarly for smaller parameters or more strongly curved geometries. A symbolic weighted operation-count model indicates comparable costs for forming the two expansion kernels. The Taylor expansions have the same formal order, while QBMAX reduces the truncation-error constants in the flat-boundary model and generally exhibits less convergence degradation in the reported experiments.
ER  -