Continuum kinetic model for simulating low-collisionality regimes in plasmas

• Main Authors: Vogman, G. V., Colella, P.
• Format: Conference Proceeding
• Language: English
• Subjects: Computational modeling; Damping; Distribution functions; Kinetic theory; Mathematical model; Plasmas; Potential energy
• ISSN: 2158-4915; 2158-4923
• DOI: 10.1109/PPC.2013.6627592

Continuum kinetic models, such as Maxwell-Boltzmann, present a viable alternative to particle-in-cell (PIC) models because they can be cast in conservation form and are not susceptible to noise. By treating the associated phase space distribution function as a continuous incompressible fluid occupying a volume of position-velocity space, evolution of the distribution function is determined by solving a 6-D advection equation. In cases where collision terms are negligible, the Boltzmann model is reduced to a Vlasov model. A 4th-order accurate continuum kinetic Vlasov model has been developed in one spatial and one velocity dimension to address the challenges associated with solving a hyperbolic governing equation in multidimensional phase space. The governing equation is cast in conservation law form and solved with a finite volume representation. Adaptive mesh refinement (AMR) is used to allow for efficient use of computational resources while maintaining desired levels of resolution. Consequently, with AMR the model is able to capture the fine structures that develop in the distribution function as it evolves in time, while using low resolution in the tail of the distribution function. The model is tested on several plasma phenomena including: weak and strong Landau damping and the two-stream instability. Conservation properties of the method are investigated.