Algorithm Development for the Two-Fluid Plasma Model
Abstract
A new algorithm is developed based on the two-fluid plasma model that is more physically accurate and capable than MHD models. The algorithm uses high-order spatial and temporal accuracy to simulate time-dependent, three-dimensional plasma phenomena. High-order spatial accuracy is accomplished using a discontinuous Galerkin finite element method that has provided up to 16th order accuracy. The temporal evolution is advanced using a 3rd order Runge-Kutta method. The numerical fluxes are calculated using an approximate Riemann solver based on the two-fluid plasma model. The source terms of the two-fluid plasma model couple the electron and ion fluids to the electromagnetic fields. The simultaneous solution and evolution must be tightly coupled to prevent unstable numerical oscillations. Asymptotic approximations are individually applied to the two-fluid plasma model to approach the Hall-MHD plasma model. An improved method of plasma simulation is found by using the two-fluid plasma model with an artificially increased electron to ion mass ratio and decreased speed of light. Multiscale effects are discovered in current-carrying plasma where small-scale electron instabilities lead to ion shocks that produce large-scale disruptions on the plasma.
Document Details
- Document Type
- Technical Report
- Publication Date
- Feb 17, 2009
- Accession Number
- ADA495369
Entities
People
- Uri Shumlak
Organizations
- University of Washington