Uniform error bound of a conservative fourth-order compact finite difference scheme for the Zakharov system in the subsonic regime

We present rigorous analysis on the error bound and conservation laws of a fourth-order compact finite difference scheme for Zakharov system (ZS) with a dimensionless parameter epsilon is an element of (0,1], which is inversely proportional to the acoustic speed. In the subsonic limit regime, i.e., 0 < epsilon << 1, the solutions have highly oscillatory waves and outgoing initial layers due to the perturbation from wave operator in ZS and the incompatibility of the initial data. The solutions propagate with O(epsilon) wavelength in time, O(1/epsilon) speed in space, and O(epsilon(2)) and O(1) amplitudes for well-prepared and ill-prepared initial data respectively. The high oscillation brings noticeable difficulties in analyzing the error bounds of numerical methods to the ZS. In this work, with h the mesh size and tau the time step, we give a uniform error bound h(4)+ tau(2 alpha dagger/3) for the well- and less-ill-prepared initial data and an error bound h(4)/epsilon + tau(2)/epsilon(3) for the ill-prepared initial data with tools including energy methods and cut-off techniques. The compact scheme provides much better spatial resolution than general second-order methods and reduces the computational cost a lot. Numerical simulations are also provided to confirm our theoretical analysis.

Automated summary

Rigorous analysis was conducted on the error bound and conservation laws of a fourth-order compact finite difference scheme for the Zakharov system, providing error bounds for well-prepared and ill-prepared initial data, with high oscillation of solutions posing challenges in analyzing error bounds of numerical methods.