UNIFORM ERROR ESTIMATES OF THE CONSERVATIVE FINITE DIFFERENCE METHOD FOR THE ZAKHAROV SYSTEM IN THE SUBSONIC LIMIT REGIME

We rigorously analyze the error estimates of the conservative finite difference method (CNFD) for the Zakharov system (ZS) with a dimensionless parameter epsilon is an element of (0, 1], which is inversely proportional to the ion acoustic speed. When epsilon -> 0(+), ZS collapses to the standard nonlinear Schrodinger equation (NLS). In the subsonic limit regime, i.e., epsilon -> 0(+), there exist highly oscillatory initial layers in the solution. The initial layers propagate with O(epsilon) wavelength in time, O(1) and O(epsilon(2)) amplitudes, for the ill-prepared initial data and well-prepared initial data, respectively. This oscillatory behavior brings significant difficulties in analyzing the errors of numerical methods for solving the Zakharov system. In this work, we show the CNFD possesses the error bounds h(2)/epsilon + tau(2)/epsilon(3) in the energy norm for the ill-prepared initial data, where h is mesh size and tau is time step. For the well-prepared initial data, CNFD is uniformly convergent for epsilon is an element of(0, 1], with second-order accuracy in space and O(tau(4/3)) accuracy in time. The main tools involved in the analysis include cut-off technique, energy methods, epsilon-dependent error estimates of the ZS, and epsilon-dependent error bounds between the numerical approximate solution of the ZS and the solution of the limit NLS. Our approach works in one, two and three dimensions, and can be easily extended to the generalized Zakharov system and nonconservative schemes. Numerical results suggest that the error bounds are sharp for the plasma densities and the error bounds of the CNFD for the electric fields are the same as those of the splitting methods.