Numerical simulation of oscillations in a cold but viscous plasma

The effect of viscosity on non-relativistic oscillations of cold plasma is numerically analyzed. For this purpose, an implicit difference scheme of the MacCormack type is constructed, which has a weaker restriction on stability in comparison with the explicit scheme. The scheme is implemented without iterations, which increases its computational efficiency tenfold. It is shown that taking into account the plasma viscosity can cause not only attenuation of the amplitude of plasma oscillations, but also a change in the shape of the solution. With an increase in the viscosity coefficient, a saddle point is observed in the solution, which is preserved in time.