Lightning Algorithm for Inhomogeneous Modified Helmholtz Equations on Cornered Domains and Their Applications
Abstract
We develop in this paper an efficient solver for the inhomogeneous modified Helmholtz equations in two-dimensional complex domains with particular attention to resolving corner singularities. We represent the solution as a summation of a particular solution and the solution of the modified Helmholtz equation. Built upon suitable domain embedding methods, we compute the particular solution using the Fourier spectral method and/or Legendre--Fourier spectral method. To treat the corner singularities, we develop a lightning algorithm by employing an expansion of the modified Bessel functions of the second kind, coupled with tapered exponentially clustering poles, to approximate the corner singularities, and using an expansion of the modified Bessel functions of the first kind to approximate the smooth part of the solution. Finally, we determine the unknown expansion coefficients by solving a least-squares problem to satisfy the boundary conditions. Numerical experiments demonstrate that root-exponential convergence can be achieved for the modified Helmholtz equation with corner singularities. To further demonstrate the applicability of the proposed algorithm, it has been successfully applied to solve some typical linear and nonlinear parabolic equations. We have also made a first attempt in extending the lightning algorithm to a three-dimensional modified Helmholtz equation.
Related
- On the Frequencies of Scattering on Traps with Micro-Inhomogeneous Structures—Helmholtz Resonators
- Inverse source problem for the stochastic fractional Helmholtz equation in inhomogeneous media
- Wavenumber-explicit stability of three-dimensional Cartesian UPMLs with piecewise-constant profiles for Helmholtz and Maxwell equations
- Freidlin-Wentzell type exit-time estimates for time-inhomogeneous diffusions and their applications
- Quasi-Helmholtz Calderón Multiplicative Preconditioning for Higher-Order Global Multi-Trace Integral Equations
BibTeX
@misc{kong2026lightning,
title = {{Lightning Algorithm for Inhomogeneous Modified Helmholtz Equations on Cornered Domains and Their Applications}},
author = {Desong Kong and Jie Shen},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.26540},
eprint = {2610.26540},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.26540},
}
RIS
TY - UNPB
TI - Lightning Algorithm for Inhomogeneous Modified Helmholtz Equations on Cornered Domains and Their Applications
AU - Desong Kong
AU - Jie Shen
PY - 2026
DA - 2026/10/19
PB - arXiv
DO - 10.48550/arxiv.2610.26540
UR - https://arxiv.org/abs/2610.26540
AB - We develop in this paper an efficient solver for the inhomogeneous modified Helmholtz equations in two-dimensional complex domains with particular attention to resolving corner singularities. We represent the solution as a summation of a particular solution and the solution of the modified Helmholtz equation. Built upon suitable domain embedding methods, we compute the particular solution using the Fourier spectral method and/or Legendre--Fourier spectral method. To treat the corner singularities, we develop a lightning algorithm by employing an expansion of the modified Bessel functions of the second kind, coupled with tapered exponentially clustering poles, to approximate the corner singularities, and using an expansion of the modified Bessel functions of the first kind to approximate the smooth part of the solution. Finally, we determine the unknown expansion coefficients by solving a least-squares problem to satisfy the boundary conditions. Numerical experiments demonstrate that root-exponential convergence can be achieved for the modified Helmholtz equation with corner singularities. To further demonstrate the applicability of the proposed algorithm, it has been successfully applied to solve some typical linear and nonlinear parabolic equations. We have also made a first attempt in extending the lightning algorithm to a three-dimensional modified Helmholtz equation.
ER -