Abstract
A variant of the generalized Ohm's law that is suited for a weakly-ionized multicomponent plasma in a magnetic field is here derived. The latter takes into consideration the current due to the non-neutrality of the plasma, the current due to the Hall effect, and the currents due to the ion slip associated with each type of ion. An equation for the electric field potential applicable to a non-neutral multicomponent plasma in the presence of a magnetic field is then presented. Despite some similarities between the potential equation and the Poisson equation, it is argued that the discretization of the potential equation cannot be accomplished in the same manner by using only central differences. It is here proven (and subsequently verified through a test case) that when the plasma exhibits conjunctly a high Hall parameter and a high electrical conductivity gradient, the centered stencils introduce spurious oscillations which can lead to severe numerical error. A novel discretization of the potential equation consisting of a blend of central and upwind differences is then presented. The proposed scheme is consistently monotonic for any value of the Hall parameter and is second-order accurate except in the vicinity of discontinuities.
Original language | English (US) |
---|---|
Pages (from-to) | 1439-1453 |
Number of pages | 15 |
Journal | Journal of Computational Physics |
Volume | 230 |
Issue number | 4 |
DOIs | |
State | Published - Feb 20 2011 |
All Science Journal Classification (ASJC) codes
- Numerical Analysis
- Modeling and Simulation
- Physics and Astronomy (miscellaneous)
- General Physics and Astronomy
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics
Keywords
- Electric field potential equation
- Flux limiters
- Generalized Ohm's law
- Hall parameter
- Ion slip
- MHD
- Monotonicity
- Multicomponent plasma
- Weakly-ionized plasma