A Bernoulli phase-fitted finite difference method and wavenumber-explicit analysis for the one-dimensional Helmholtz equation
For computational scientists solving Helmholtz problems, this method offers a pollution-free scheme with uniform error bounds, improving upon standard finite difference methods.
The paper proposes a Bernoulli phase-fitted finite difference method for the one-dimensional Helmholtz equation that is exact for plane waves, eliminating numerical dispersion and artificial reflection. It proves wavenumber-explicit stability and second-order convergence with pollution-free error bounds under fixed-resolution conditions.
We propose a Bernoulli phase-fitted (BPF) finite difference method for the Helmholtz equation on the interval $(0, L)$ with impedance boundary conditions. The scheme is derived from a complexified Scharfetter--Gummel discretization of the one-way factorization of the Helmholtz operator. It yields both a phase-fitted interior discretization and exact discrete impedance boundary closures. For the homogeneous problem, the method is exact for plane waves, so the scheme introduces neither numerical dispersion in the interior nor artificial reflection at the boundaries. For the inhomogeneous problem, we prove well-posedness, derive wavenumber-explicit stability estimates, and establish second-order consistency and convergence valid for all $kh\notinπ\mathbb Z$, where $k$ is the wavenumber and $h$ the grid size. In particular, under the fixed-resolution condition $kh\le s_0$ for some $0<s_0<π$ together with $kL\geπ$, the constants in the error bounds remain uniform with respect to the wavenumber, yielding a pollution-free convergence theory in the principal Nyquist regime. Numerical experiments confirm the theoretical analysis and show favorable performance compared with standard and dispersion-corrected finite difference methods.