Advanced physics calculations using a multi-fluid plasma model

Abstract The multi-fluid plasma model is derived from moments of the Boltzmann equation and typically has two fluids representing electron and ion species. Large mass differences between electrons and ions introduce disparate temporal and spatial scales and require a numerical algorithm with sufficient accuracy to capture the multiple scales. Source terms of the multi-fluid plasma model couple the fluids to themselves (interspecies interactions) and to the electromagnetic fields. The numerical algorithm must treat the inherent stiffness introduced by the multiple physical effects of the model and tightly couple the source terms of the governing equations. A discontinuous Galerkin method is implemented for the spatial representation. Time integration is investigated using explicit, implicit, semi-implicit methods. Semi-implicit treatment is accomplished using a physics-based splitting. The algorithm is applied to study drift turbulence in field reversed configuration plasmas to illustrate the physical accuracy of the model. The algorithm is also applied to plasma sheath formation which demonstrates Langmuir wave propagation.

[1]  Ammar Hakim.  Two-Fluid Physics and Field Reversed Configurations , 2006 .

[2]  Uri Shumlak,et al.  Two-fluid physics and field-reversed configurationsa) , 2007 .

[3]  Uri Shumlak,et al.  An Implicit Scheme for Nonideal Magnetohydrodynamics , 1997 .

[4]  L. L. Lao,et al.  Magnetohydrodynamic simulations of massive gas injection into Alcator C-Mod and DIII-D plasmasa) , 2008 .

[5]  Uri Shumlak,et al.  Approximate Riemann solver for the two-fluid plasma model , 2003 .

[6]  Yasushi Matsunaga,et al.  Kinetic simulation of nonlinear phenomena of an ion acoustic wave in gas discharge plasma with convective scheme , 2001 .

[7]  Boris B. Kadomtsev,et al.  Reviews of Plasma Physics , 2012 .

[8]  Uri Shumlak,et al.  A discontinuous Galerkin method for the full two-fluid plasma model , 2005, Comput. Phys. Commun..

[9]  Michael Hesse,et al.  Geospace Environmental Modeling (GEM) magnetic reconnection challenge , 2001 .

[10]  G. Knorr,et al.  The integration of the vlasov equation in configuration space , 1976 .

[11]  K. Murawski,et al.  Numerical simulations of impulsively generated vertical oscillations in a solar coronal arcade loop , 2006 .

[12]  S. Fromang,et al.  MHD simulations of the magnetorotational instability in a shearing box with zero net flux: the case Pm = 4 , 2010, 1004.2382.

[13]  van der Jjam Joost Mullen,et al.  Kinetic simulation of an extreme ultraviolet radiation driven plasma near a multilayer mirror , 2006 .

[14]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[15]  Claus-Dieter Munz,et al.  A three-dimensional finite-volume solver for the Maxwell equations with divergence cleaning on unstructured meshes , 2000 .

[16]  Chi-Wang Shu,et al.  Runge–Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems , 2001, J. Sci. Comput..

[17]  S. I. Braginskii.  Transport Processes in a Plasma , 1965 .

[18]  D. A. Knoll,et al.  An Implicit Energy-Conservative 2D Fokker—Planck Algorithm , 2000 .

[19]  Chi-Wang Shu,et al.  The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV. The multidimensional case , 1990 .

[20]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[21]  Uri Shumlak,et al.  A high resolution wave propagation scheme for ideal Two-Fluid plasma equations , 2006, J. Comput. Phys..

[22]  Subrata Roy,et al.  Modeling low pressure collisional plasma sheath with space-charge effect , 2003 .

[23]  Uri Shumlak,et al.  Nonlinear full two-fluid study of m=0 sausage instabilities in an axisymmetric Z pinch , 2006 .

[24]  Bhuvana Srinivasan.  Numerical Methods for 3-dimensional Magnetic Confinement Configurations using Two-Fluid Plasma Equations , 2010 .