A hybrid-Vlasov model based on the current advance method for the simulation of collisionless magnetized plasma
- Dipartimento di Fisica and CNISM, Universita degli Studi della Calabria, Ponte P. Bucci, Cubo 31C, 87036 Arcavacata di Rende (Serbia and Montenegro) (Italy)
- Institute of Atmospheric Physics, AS CR, Prague (Czech Republic)
- Dipartimento di Fisica and CNISM, Universita di Pisa, Pisa (Italy)
- LESIA - Observatoire de Paris, Section de Meudon 5, place Jules Janssen, 92195 Meudon Cedex (France)
We present a numerical scheme for the integration of the Vlasov-Maxwell system of equations for a non-relativistic plasma, in the hybrid approximation, where the Vlasov equation is solved for the ion distribution function and the electrons are treated as a fluid. In the Ohm equation for the electric field, effects of electron inertia have been retained, in order to include the small scale dynamics up to characteristic lengths of the order of the electron skin depth. The low frequency approximation is used by neglecting the time derivative of the electric field, i.e. the displacement current in the Ampere equation. The numerical algorithm consists in coupling the splitting method proposed by Cheng and Knorr in 1976 [C.Z. Cheng, G. Knorr, J. Comput. Phys. 22 (1976) 330-351.] and the current advance method (CAM) introduced by Matthews in 1994 [A.P. Matthews, J. Comput. Phys. 112 (1994) 102-116.] In its present version, the code solves the Vlasov-Maxwell equations in a five-dimensional phase space (2-D in the physical space and 3-D in the velocity space) and it is implemented in a parallel version to exploit the computational power of the modern massively parallel supercomputers. The structure of the algorithm and the coupling between the splitting method and the CAM method (extended to the hybrid case) is discussed in detail. Furthermore, in order to test the hybrid-Vlasov code, the numerical results on propagation and damping of linear ion-acoustic modes and time evolution of linear elliptically polarized Alfven waves (including the so-called whistler regime) are compared to the analytical solutions. Finally, the numerical results of the hybrid-Vlasov code on the parametric instability of Alfven waves are compared with those obtained using a two-fluid approach.
- OSTI ID:
- 20991601
- Journal Information:
- Journal of Computational Physics, Vol. 225, Issue 1; Other Information: DOI: 10.1016/j.jcp.2007.01.001; PII: S0021-9991(07)00002-2; Copyright (c) 2007 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
Numerical study of a Vlasov equation for systems with interacting particles
Asymmetric kinetic equilibria: Generalization of the BAS model for rotating magnetic profile and non-zero electric field
Related Subjects
GENERAL PHYSICS
ALFVEN WAVES
ALGORITHMS
ANALYTICAL SOLUTION
APPROXIMATIONS
BOLTZMANN-VLASOV EQUATION
COMPUTERIZED SIMULATION
DISTRIBUTION FUNCTIONS
ELECTRIC FIELDS
ELECTRONS
MATHEMATICAL EVOLUTION
MOMENT OF INERTIA
PARAMETRIC INSTABILITIES
PHASE SPACE
RELATIVISTIC PLASMA
SUPERCOMPUTERS