pith. sign in

arxiv: 2605.21144 · v1 · pith:JXBYOEXKnew · submitted 2026-05-20 · 🧮 math.NA · cs.NA

A Bernoulli phase-fitted finite difference method and wavenumber-explicit analysis for the one-dimensional Helmholtz equation

Ansgar J\"ungel , Panchi Li , Zhiwei Sun , Zhiwen Zhang This is my paper

Pith reviewed 2026-05-21 01:52 UTC · model grok-4.3

classification 🧮 math.NA cs.NA MSC 65N06
keywords finite difference methodHelmholtz equationphase-fitted discretizationnumerical dispersionwavenumber-explicit analysispollution-free convergenceimpedance boundary conditions
0
0 comments X

The pith

A new finite difference scheme for the one-dimensional Helmholtz equation reproduces plane waves exactly and keeps error constants uniform in wavenumber under fixed resolution.

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

The paper introduces a Bernoulli phase-fitted finite difference method for the Helmholtz equation on an interval with impedance boundary conditions. The scheme comes from a complexified Scharfetter-Gummel discretization of the one-way factorization of the Helmholtz operator. For homogeneous problems the discrete solution matches exact plane waves, so there is no numerical dispersion inside the domain and no artificial reflection at the boundaries. For inhomogeneous problems the authors prove second-order consistency and convergence with explicit dependence on the wavenumber, and they show that the constants in the error bounds stay bounded independently of wavenumber when kh is bounded away from multiples of pi and the interval length satisfies kL at least pi.

Core claim

The Bernoulli phase-fitted finite difference method, obtained by applying a complexified Scharfetter-Gummel discretization to the one-way factorization of the Helmholtz operator, produces an interior stencil and boundary closures that are exact for plane-wave solutions of the homogeneous equation. The scheme therefore introduces neither numerical dispersion nor artificial reflections. For the inhomogeneous problem the method is well-posed, second-order consistent, and convergent for all kh not equal to integer multiples of pi; the constants in the error estimates remain independent of the wavenumber under the fixed-resolution condition kh less than or equal to some s0 less than pi together,

What carries the argument

Complexified Scharfetter-Gummel discretization of the one-way factorization of the Helmholtz operator, which generates phase-fitted interior stencils and exact discrete impedance boundary conditions that match plane waves exactly.

If this is right

  • The homogeneous Helmholtz problem is solved exactly on any grid for any wavenumber.
  • Error bounds for the inhomogeneous problem remain uniform in wavenumber when kh is bounded away from multiples of pi and kL is at least pi.
  • The method is free of the pollution effect in the principal Nyquist regime.
  • Numerical tests show the scheme performs better than standard and dispersion-corrected finite difference methods.

Where Pith is reading between the lines

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

  • The exactness for plane waves could serve as a reference to benchmark other Helmholtz discretizations.
  • Similar factorizations might allow extension to variable-coefficient or higher-dimensional cases.
  • Uniform constants suggest the method can handle arbitrarily high frequencies without mesh refinement that scales with wavenumber.

Load-bearing premise

The analysis requires that the product of wavenumber and grid size avoids integer multiples of pi.

What would settle it

Solve a homogeneous Helmholtz problem whose exact solution is a known plane wave, apply the scheme on a grid with kh not a multiple of pi, and check whether the discrete solution equals the exact solution up to round-off error.

Figures

Figures reproduced from arXiv: 2605.21144 by Ansgar J\"ungel, Panchi Li, Zhiwei Sun, Zhiwen Zhang.

Figure 1
Figure 1. Figure 1: Comparison between the reference solution and the numerical solution for the plane-wave test on a coarse mesh. 5.2. A smooth manufactured-solution test. To directly validate the convergence es￾timate of Theorem 4.3, we consider the smooth manufactured solution u(x) = e ikx + r(x), r(x) := x 4 (1 − x) 4 for x ∈ (0, 1). The source term is defined by f(x) = r ′′(x) + k 2 r(x) or, more explicitly, f(x) = 12x 2… view at source ↗
Figure 2
Figure 2. Figure 2: Results for the smooth-source test: convergence with respect to h in the relative V -norm (left) and the relative L ∞-norm (right). 5.3. A nonsmooth-source test. We consider a nonsmooth source term to illustrate the robustness of the BPF scheme beyond the regularity assumptions required in Theorem 4.3. We take fx) = 50 for |x−0.5| ≤ 1/9 and f(x) = 0 else, together with the nonhomogeneous impedance boundary… view at source ↗
Figure 3
Figure 3. Figure 3: Results for the nonsmooth-source test: convergence with respect to h in the relative V -norm (left) and the relative L ∞-norm (right) 4.3. The reference solution is computed on the fine mesh h = 2−18. For the present benchmark, the numerical data indicate that ∥u∥L∞ and ∥u∥L2 scale like k −1 , whereas the discrete energy norm ∥u∥V is nearly independent of k. We therefore use the relative V -norm error as t… view at source ↗
Figure 4
Figure 4. Figure 4: Relative L ∞-error versus k for the classical FD method (left), dispersion-corrected FD method (middle), and BPF scheme evaluated at fixed values of kh (right). 6. Conclusion and outlook In this paper, we introduced a Bernoulli phase-fitted scheme for the one-dimensional Helmholtz equation with impedance boundary conditions. The scheme is derived from a complexified Scharfetter–Gummel discretization of the… view at source ↗
read the original abstract

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\pi\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<\pi$ together with $kL\ge\pi$, 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.

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 / 2 minor

Summary. The manuscript proposes a Bernoulli phase-fitted (BPF) finite difference method for the one-dimensional Helmholtz equation on (0, L) subject to impedance boundary conditions. The scheme is obtained via a complexified Scharfetter-Gummel discretization of the one-way factorization of the Helmholtz operator. For the homogeneous problem the method is exact for plane waves, eliminating both interior numerical dispersion and artificial boundary reflections. For the inhomogeneous problem the authors establish well-posedness, wavenumber-explicit stability bounds, second-order consistency, and convergence; under the fixed-resolution regime kh ≤ s0 < π together with kL ≥ π the error constants remain uniform in the wavenumber, provided kh ∉ πℤ.

Significance. If the stated results hold, the work supplies a concrete pollution-free second-order theory for the 1D Helmholtz equation in the principal Nyquist regime, a result that remains rare in the literature. The exactness property for plane waves is a strong structural feature that directly removes the usual sources of phase error and reflection. The explicit wavenumber dependence in the stability and convergence estimates, together with the uniform constants under the stated mesh and domain conditions, constitutes a useful addition to the analysis of high-frequency discretizations.

minor comments (2)
  1. [§3] §3 (or the section containing the discrete Green's function estimate): the proof that the constants remain uniform when kL ≥ π relies on a lower bound for the discrete Green's function; a short remark clarifying how the condition kL ≥ π enters the argument would help readers follow the uniformity claim.
  2. [Numerical experiments] Numerical experiments section: the reported comparisons with standard and dispersion-corrected schemes would be strengthened by stating the precise values of kh and kL used in each test, together with the observed L^∞ or L^2 error magnitudes.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of the manuscript, the recognition of the exactness property for plane waves, and the recommendation for minor revision. We are pleased that the wavenumber-explicit stability and uniform convergence results under the fixed-resolution regime were viewed as a useful contribution.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The derivation begins with a complexified Scharfetter-Gummel discretization of the one-way factorization of the Helmholtz operator, which directly produces both the phase-fitted interior scheme and the exact discrete impedance boundary closures. Exactness for plane waves on the homogeneous problem is an immediate algebraic consequence of this construction rather than a fitted or self-referential claim. Well-posedness, wavenumber-explicit stability, second-order consistency, and convergence (including uniformity of constants under kh ≤ s0 < π and kL ≥ π) are established by direct proofs that control the discrete Green's function under the explicit proviso kh ∉ πℤ; these arguments do not invoke self-citations, imported uniqueness theorems, or ansatzes smuggled from prior work. The scheme is therefore self-contained against its own stated hypotheses and external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Review based on abstract only; ledger populated from elements explicitly named in the abstract.

axioms (1)
  • domain assumption The Helmholtz operator admits a one-way factorization suitable for complexified Scharfetter-Gummel discretization
    Invoked to derive the phase-fitted interior scheme and exact boundary closures.

pith-pipeline@v0.9.0 · 5752 in / 1163 out tokens · 28977 ms · 2026-05-21T01:52:27.716214+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

47 extracted references · 47 canonical work pages

  1. [1]

    and Herrera, C\'

    Cai, Zhiqiang and Doktorova, Anastassia and Falgout, Robert D. and Herrera, C\'. Efficient Shallow. SIAM Journal on Scientific Computing , volume =

    work page

  2. [2]

    Multigrid methods for

    Ernst, Oliver G and Gander, Martin J , booktitle=. Multigrid methods for

    work page

  3. [3]

    Multigrid methods for

    Ernst, Oliver G and Gander, Martin J , journal=. Multigrid methods for

    work page

  4. [4]

    Exact finite difference schemes for solving

    Wong, Yau Shu and Li, Guangrui , journal=. Exact finite difference schemes for solving

    work page

  5. [5]

    Asymptotic dispersion correction in general finite difference schemes for

    Cocquet, Pierre-Henri and Gander, Martin J , journal=. Asymptotic dispersion correction in general finite difference schemes for. 2024 , publisher=

    work page 2024

  6. [6]

    Is the pollution effect of the

    Babu. Is the pollution effect of the. SIAM J. Numer. Anal. , volume=. 1997 , publisher=

    work page 1997

  7. [7]

    Fourier Analysis of Finite Difference Schemes for the

    Gander, Martin J and Zhang, Hui , journal=. Fourier Analysis of Finite Difference Schemes for the

    work page

  8. [8]

    POLLUTION-FREE FINITE DIFFERENCE SCHEMES FOR NON-HOMOGENEOUS

    Wang, Kun and Wong, Yau Shu , journal=. POLLUTION-FREE FINITE DIFFERENCE SCHEMES FOR NON-HOMOGENEOUS

    work page

  9. [9]

    Ihlenburg and I

    F. Ihlenburg and I. Babu. Finite element solution of the. Comput. Math. Appl. , volume =. 1995 , issn =

    work page 1995

  10. [10]

    Melenk, J. M. and Sauter, S. , title =. SIAM J. Numer. Anal. , volume =

    work page

  11. [11]

    SIAM Journal on Numerical Analysis , volume =

    Zhou, Yu and Wu, Haijun , title =. SIAM Journal on Numerical Analysis , volume =

    work page

  12. [12]

    Melenk, J. M. and Parsania, A. and Sauter, S. , title=. J. Sci. Comput. , year=

    work page

  13. [13]

    and Spence, E

    Galkowski, J. and Spence, E. A. , title =. SIAM Rev. , volume =

    work page

  14. [14]

    Multiscale Model

    Hauck, Moritz and Peterseim, Daniel , title =. Multiscale Model. Simul. , volume =

    work page

  15. [15]

    ELIMINATING THE POLLUTION EFFECT IN

    Daniel Peterseim , journal =. ELIMINATING THE POLLUTION EFFECT IN

    work page

  16. [16]

    Panchi Li and Zhiwen Zhang , year=. A. 2507.16647 , archivePrefix=

    work page arXiv

  17. [17]

    Deraemaeker, Arnaud and Babuška, Ivo and Bouillard, Philippe , title =. Int. J. Numer. Meth. Engin. , volume =

    work page

  18. [18]

    Feng, Xiaobing and Wu, Haijun , title =. SIAM J. Numer. Anal. , volume =

    work page

  19. [19]

    hp-DISCONTINUOUS

    Xiaobing Feng and Haijun Wu , journal =. hp-DISCONTINUOUS

    work page

  20. [20]

    Wavenumber-explicit convergence of the. Comput. Math. Appl. , volume =. 2022 , issn =

    work page 2022

  21. [21]

    Convergence analysis for finite element discretizations of the

    Melenk, J and Sauter, Stefan , journal=. Convergence analysis for finite element discretizations of the

    work page

  22. [22]

    SIAM Journal on Numerical Analysis , volume =

    Ma, Chupeng and Alber, Christian and Scheichl, Robert , title =. SIAM Journal on Numerical Analysis , volume =

    work page

  23. [23]

    and Wang, Yixuan , title =

    Chen, Yifan and Hou, Thomas Y. and Wang, Yixuan , title =. Multiscale Modeling & Simulation , volume =

    work page

  24. [24]

    Multiscale Modeling & Simulation , volume =

    Fu, Shubin and Li, Guanglian and Craster, Richard and Guenneau, Sebastien , title =. Multiscale Modeling & Simulation , volume =

    work page

  25. [25]

    Multiscale Modeling & Simulation , volume =

    Ohlberger, Mario and Verfurth, Barbara , title =. Multiscale Modeling & Simulation , volume =

    work page

  26. [26]

    Generalized Finite Element Methods:

    Babu. Generalized Finite Element Methods:. SIAM Journal on Numerical Analysis , volume =

    work page

  27. [27]

    Is the Pollution Effect of the

    Babu. Is the Pollution Effect of the. SIAM Journal on Numerical Analysis , volume =

    work page

  28. [28]

    A dispersion minimizing finite difference scheme for the. Comput. Math. Appl. , volume =. 2017 , issn =

    work page 2017

  29. [29]

    A dispersion minimizing finite difference scheme and preconditioned solver for the. J. Comput. Phys. , volume =. 2012 , issn =

    work page 2012

  30. [30]

    A generalized optimal fourth-order finite difference scheme for a. J. Comput. Appl. Math. , volume =. 2021 , issn =

    work page 2021

  31. [31]

    An optimal compact sixth-order finite difference scheme for the. Comput. Math. Appl. , volume =. 2018 , issn =

    work page 2018

  32. [32]

    Sixth-Order Compact Finite Difference Method for. Commun. Comput. Phys. , author =. 2023 , month=

    work page 2023

  33. [33]

    Dastour, Hatef and Liao, Wenyuan , title=. Numer. Algor. , year=

    work page

  34. [34]

    A finite difference method with optimized dispersion correction for the

    Cocquet, Pierre-Henri and Gander, Martin J and Xiang, Xueshuang , booktitle=. A finite difference method with optimized dispersion correction for the. 2017 , organization=

    work page 2017

  35. [35]

    Dispersion correction for

    Cocquet, Pierre-Henri and Gander, Martin J and Xiang, Xueshuang , booktitle=. Dispersion correction for. 2018 , organization=

    work page 2018

  36. [36]

    Closed form dispersion corrections including a real shifted wavenumber for finite difference discretizations of

    Cocquet, Pierre-Henri and Gander, Martin J and Xiang, Xueshuang , journal=. Closed form dispersion corrections including a real shifted wavenumber for finite difference discretizations of. 2021 , publisher=

    work page 2021

  37. [37]

    Scharfetter and H

    D. Scharfetter and H. Gummel , title =. IEEE Trans. Electron. Devices , volume =. 1969 , pages =

    work page 1969

  38. [38]

    Farrell and E

    P. Farrell and E. Gartland , title =. Computational Methods for Boundary and Interior Layers in Several Dimensions , editor =. 1991 , address =

    work page 1991

  39. [39]

    Bessemoulin--Chatard , title =

    M. Bessemoulin--Chatard , title =. Numer. Math. , volume =. 2012 , pages =

    work page 2012

  40. [40]

    Chainais-Hillairet and J

    C. Chainais-Hillairet and J. Droniou , title =. IMA J. Numer. Anal. , volume =. 2011 , pages =

    work page 2011

  41. [41]

    Canc\'es and C

    C. Canc\'es and C. Chainais-Hillairet and J. Fuhrmann and B. Gaudeul , title =. IMA J. Numer. Anal. , volume =. 2019 , pages =

    work page 2019

  42. [42]

    Finite Element Solution of the

    Ihlenburg, Frank and Babu. Finite Element Solution of the. SIAM J. Numer. Anal. , volume =

    work page

  43. [43]

    On Accuracy Conditions for the Numerical Computation of Waves , author =. J. Comput. Phys. , volume =

    work page

  44. [44]

    Application of an Ultra Weak Variational Formulation of Elliptic

    Cessenat, Olivier and Despres, Bruno , year = 1998, month = feb, journal =. Application of an Ultra Weak Variational Formulation of Elliptic

    work page 1998

  45. [45]

    and Moiola, A

    Hiptmair, R. and Moiola, A. and Perugia, I. , year = 2011, month = jan, journal =. Plane Wave Discontinuous Galerkin Methods for the

    work page 2011

  46. [46]

    A Survey of

    Hiptmair, Ralf and Moiola, Andrea and Perugia, Ilaria , year = 2016, booktitle=. A Survey of

    work page 2016

  47. [47]

    Trefftz Discontinuous

    Hiptmair, Ralf and Moiola, Andrea and Perugia, Ilaria , year = 2014, month = may, journal =. Trefftz Discontinuous

    work page 2014