Particle-in-cell simulation of Buneman instability beyond quasilinear saturation
arXiv:1709.07478 · doi:10.1063/1.5006463
Abstract
Spatio-temporal evolution of Buneman instability has been followed numerically till its quasilinear quenching and beyond, using an in-house developed electrostatic 1D particle-in-cell simulation code. For different initial drift velocities and for a wide range of electron to ion mass ratios (m/M), growth rate obtained from simulation agrees well with the numerical solution of the fourth order dispersion relation. Quasi-linear saturation of Buneman instability occurs when ratio of electrostatic field energy density () to initial electron drift kinetic energy density () reaches up to a constant value, which as predicted by Hirose [Plasma Physics 20, 481(1978)], is independent of initial electron drift velocity but depends on electron to ion mass ratio m/M as . This result stands verified in our simulations. Growth of the instability beyond the first saturation (quasilinear saturation ) till its final saturation [Ishihara et. al., PRL 44, 1404(1980)] follows an algebraic scaling with time. In contrast to the quasilinear saturation, the ratio of final saturated electrostatic field energy density to initial kinetic energy density, is relatively independent of electron to ion mass ratio and is found to depend only on the initial drift velocity. Beyond the final saturation, electron phase space holes coupled to large amplitude ion solitary waves, a state known as coupled hole-soliton , are seen in our simulations. The propagation characteristics ( amplitude - speed relation ) of these coherent modes is found to be consistent with the theory of Saeki et. al. [PRL 80, 1224(1998)].
References in corpus (3)
Cited by in corpus (5)
- The role of noise in PIC and Vlasov simulations of the Buneman instability
- Backward waves in the nonlinear regime of the Buneman instability
- On the validity of quasilinear theory applied to the electron bump-on-tail instability
- Excitation of Electrostatic Standing Wave in the Superposition of Two Counter Propagating Relativistic Whistler Waves
- Effect of Ion Motion on Breaking of Longitudinal Relativistically Strong Plasma Waves: Khachatryan mode revisited