A compact splitting scheme approach on nonuniform grids (Q2869182)

A Cauchy problem for the paraxial approximation of the Helmholtz equation is considered, i.e., the partial differential equation \(2\sqrt{-1} \kappa u_z = u_{xx} + u_{yy}\) with the wavenumber \(\kappa > 0\) and initial values at a fixed \(z=z_0\). Highly oscillatory solutions occur in the case of large constants \(\kappa\), where appropriate numerical methods are required. In the previous work by \textit{Q. Sheng} et al. [J. Comput. Appl. Math. 235, No. 15, 4452--4463 (2011; Zbl 1216.78008)], an exponential transformation is applied followed by a splitting method separating the spatial dimensions \(x\) and \(y\). In the more recent paper, the authors investigate a straightforward class of splitting methods based on the variables \(x\) and \(y\), where some coupling terms represent a degree of freedom. This finite difference method is defined on an axis-parallel grid with nonconstant step sizes and all discretisations are of second order in space. It is proved that the scheme is unconditionally asymptotically stable for highly oscillatory problems if the coupling terms and boundary terms vanish. Therefore an error amplification constant converges to one in the case of \(\kappa \rightarrow \infty\) independent of the step sizes. Finally, the authors present results of a particular Cauchy problem on a bounded domain with homogeneous Dirichlet conditions, where the splitting method yields the numerical solution.NEWLINENEWLINEFor the entire collection see [Zbl 1264.65002].