High precision methods for solving a system of cold plasma equations taking into account electron–ion collisions

High precision simulation algorithms are proposed and justified for modelling cold plasma oscillations taking into account electron–ion collisions in the non-relativistic case. The specific feature of the approach is the use of Lagrangian variables for approximate solution of the problem formulated initially in Eulerian variables. High accuracy is achieved both through the use of analytical solutions on trajectories of particles and due to sufficient smoothness of the solution in numerical integration of Cauchy problems. Numerical experiments clearly illustrate the obtained theoretical results. As a practical application, a simulation of the well-known breaking effect of multi-period relativistic oscillations is carried out. It is shown that with an increase in the collision coefficient one can observe that the breaking process slows down until it is completely eliminated.