Full Orbit PIC in NIMROD
POSTER
Abstract
The primary goal of the Plasma Science and Innovation Center (PSI Center) is to refine and optimize existing MHD codes to achieve improved predictability for emerging concept (EC) experiments. Kinetic effects have been shown to play a dominant role in some EC experiments, particularly in FRC stability\footnote{\normalsize{E.~Belova et.al. ``Numerical Study of tilt stability of prolate field-reversed configurations,'' PoP, {\bf 7}, 4996, 2000}}. The PSI Center will extend the hybrid kinetic-MHD implementation in NIMROD\footnote{\normalsize{C.C. Kim et.al. ``Hybrid Kinetic-MHD Simulations in General Geometry,'' CPC, {\bf 164}, 448, 2004}} from the drift kinetic model to the full kinetic model to include sufficient physics to accurately account for these effects, in particular large Larmour radius. The Boris push has been implemented for particles in NIMROD. However, this places a severe timestep restriction on the particle time step. We discuss several strategies for circumventing this restriction to simulate on the MHD time scale. We will show some of the preliminary results from the implementation of full orbit (Lorentz force) particles coupled to the NIMROD code.
Authors
-
Charlson Kim
University of Washington - PSI Center