Hybrid Modeling of Plasmas

Mats Holmström ()
Additional contact information
Mats Holmström: Swedish Institute of Space Physics

A chapter in Numerical Mathematics and Advanced Applications 2009, 2010, pp 451-458 from Springer

Abstract: Abstract Space plasmas are often modeled as a magnetohydrodynamic (MHD) fluid. However, many observed phenomena cannot be captured by fluid models, e.g., non-Maxwellian velocity distributions and finite gyro radius effects. Therefore kinetic models are used, where also the velocity space is resolved. This leads to a six-dimensional problem, making the computational demands of velocity space grids prohibitive. Particle in cell (PIC) methods discretize velocity space by representing the charge distribution as discrete particles, and the electromagnetic fields are stored on a spatial grid. For the study of global problems in space physics, such as the interaction of a planet with the solar wind, it is difficult to resolve the electron spatial and temporal scales. Often a hybrid model is then used, where ions are represented as particles, and electrons are modeled as a fluid. Then the ion motions govern the spatial and temporal scales of the model. Here we present the mathematical and numerical details of a general hybrid model for plasmas. All grid quantities are stored at cell centers on the grid. The most common discretization of the fields in PIC solvers is to have the electric and magnetic fields staggered, introduced by Yee [IEEE Transactions on Antennas and Propagation 14:302–307 1966]. This automatically ensures that ∇ ⋅ B = 0, down to round-off errors. Here we instead present a cell centered discretization of the magnetic field. That the standard cell centered second order stencil for ∇ × E in Faraday’s law will preserve ∇ ⋅ B = 0 was noted by Tóth [Journal of Computational Physics 161:605–652 2000]. The advantage of a cell centered discretization is ease of implementation, and the possibility to use available solvers that only handle cell centered variables. We also show that the proposed method has very good energy conservation for a simple test problem in one-, two-, and three dimensions, when compared to a commonly used algorithm.

Keywords: Solar Wind; Hybrid Simulation; Solar Wind Interaction; Hybrid Equation; Hybrid Code (search for similar items in EconPapers)
Date: 2010
References: Add references at CitEc
Citations:

There are no downloads for this item, see the EconPapers FAQ for hints about obtaining it.

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:spr:sprchp:978-3-642-11795-4_48

Ordering information: This item can be ordered from
http://www.springer.com/9783642117954

DOI: 10.1007/978-3-642-11795-4_48

Access Statistics for this chapter

More chapters in Springer Books from Springer
Bibliographic data for series maintained by Sonal Shukla () and Springer Nature Abstracting and Indexing ().