A new code for electrostatic simulation by numerical integration of the Vlasov and Ampere equations using MacCormack’s method

We present anewsimulationcode for electrostatic waves in one dimension which uses the Vlasovequation to integrate the distribution function and Ampère'sequation to integrate the electric field forward in time. Previous Vlasovcodes used the Vlasov and Poisson equations. Using Ampère'sequation has two advantages. First, boundary conditions do not have to be set on the electric field. Second, it forms a logical basis for an electromagnetic code since the time integration of the electric and magnetic fields is treated in a similar way. MacCormack'smethod is used to integrate the Vlasovequation, which was found to be easy to implement and reliable. A stability analysis is presented for the MacCormack scheme applied to the Vlasovequation. Conditions for stability are more stringent than the simple Courant–Friedrich's–Lewy (CFL) conditions for the spatial and velocity grids. We provide a simple linear function which when combined with the CFL conditions should ensure stability. Simulation results for Landau damping are in excellent agreement with numerical solutions of the linear dispersion relation for a wide range of wavelengths. The code is also able to retain phase memory as demonstrated by the recurrence effect and reproduce the effects of particle trapping. The use of Ampère'sequation enables standing and traveling waves to be produced depending on whether the current is zero or non-zero, respectively. In simulations where the initial current is non-zero and Maxwell'sequations are satisfied initially, additional standing waves may be set up, which could be important when the coupling of wave fields from a transmitter to a plasma is considered.

Details