Stable approximation of Helmholtz solutions by evanescent plane waves

Solutions of the Helmholtz equation are known to be well approximated by superpositions of propagative plane waves. This observation is the foundation of successful Trefftz methods. However, when too many plane waves are used, the computation of the expansion is known to be numerically unstable. We explain how this effect is due to the presence of exponentially large coefficients in the expansion and can drastically limit the efficiency of the approach. In this work, we show that the Helmholtz solutions on a disk can be exactly represented by a continuous superposition of evanescent plane waves, generalizing the standard Herglotz representation. Here, by evanescent plane waves, we mean exponential plane waves with complex-valued propagation vector, whose absolute value decays exponentially in one direction. In addition, the density in this representation is proved to be uniformly bounded in a suitable weighted Lebesgue norm, hence overcoming the instability observed with propagative plane waves and paving the way for stable discrete expansions. In view of practical implementations, discretization strategies are investigated. We construct suitable finite-dimensional sets of evanescent plane waves using sampling strategies in a parametric domain. Provided one uses sufficient oversampling and regularization, the resulting approximations are shown to be both controllably accurate and numerically stable, as supported by numerical evidence.