Analysis and simulation of a variational stabilization for the Helmholtz equation with noisy Cauchy data

This article considers a Cauchy problem of Helmholtz equations whose solution is well known to be exponentially unstable with respect to the inputs. In the framework of variational quasi-reversibility method, a Fourier truncation is applied to appropriately perturb the underlying problem, which allows us to obtain a stable approximate solution. The corresponding approximate problem is of a hyperbolic equation, which is also a crucial aspect of this approach. Error estimates between the approximate and true solutions are derived with respect to the noise level. From this analysis, the Lipschitz stability with respect to the noise level follows. Some numerical examples are provided to see how our numerical algorithm works well.