Optimization of calculations in modeling the breaking of plasma oscillations

The proposed method of optimization of calculations is of a combined nature. First, using a piecewise uniform grid, a difference scheme of the second order of accuracy of the McCormack type is constructed, and then — to weaken the stability condition — an implicit version of the scheme is used, which even in the nonlinear case allows an iteration-free implementation. When constructing a grid, a preliminary study is required, which is based on the analytical properties of the solution of the differential problem. Based on the calculations carried out, it can be concluded that the potential reduction in the amount of calculations is dozens of times without a noticeable loss of accuracy.