Self-consistent hybrid neoclassical-magnetohydrodynamic simulations of axisymmetric plasmas
- Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States); Princeton Univ., NJ (United States). Dept. of Astrophysical Sciences
Neoclassical effects (e.g., conductivity reduction and bootstrap currents) have a profound impact on many magnetohydrodynamic (MHD) instabilities in toroidally-confined plasmas, including tearing modes, edge-localized modes, and resistive wall modes. High-fidelity simulations of such phenomena require a multiphysics code that self-consistently couples the kinetic and fluid models. We review a hybrid formulation from the recent literatureAB that is appropriate for such studies. In particular, the formulation uses a set of time-dependent drift-kinetic equations (DKEs) to advance the non-Maxwellian part of the electron and ion distribution functions (fNM) with linearized Fokker-Planck-Landau collision operators. The form of the DKEs used were derived in a Chapman-Enskog-like fashion, ensuring that fNM carries no density, momentum, or temperature. Rather, these quantities are contained within the background Maxwellian and are evolved by a set of MHD equations which are closed by moments of fNM. We then present two DKE solvers based upon this formulation in axisymmetric toroidal geometries. The Neoclassical Ion-Electron Solver (NIES) solves the steady-state DKEs in the low-collisionality limit. Convergence and benchmark studies are discussed, providing a proof-of-principle that this new formulation can accurately reproduce results from the literature in the limit considered. We then present the DK4D code which evolves the finite-collisionality DKEs time-dependently. Computational methods used and successful benchmarks to other neoclassical models and codes are discussed. Furthermore, we couple DK4D to a reduced, transport-timescale MHD code. The resulting hybrid code is used to simulate the evolution of the current density in a large-aspect-ratio plasma in the presence of several different time-dependent pressure profiles. These simulations demonstrate the self-consistent, dynamic formation of the ohmic and bootstrap currents. In the slowly-evolving plasmas considered, these first-principle simulations are shown to verify a simpler, steady-state neoclassical model that is commonly used in transport codes. Future work will involve coupling DK4D to a spatially three-dimensional, extended MHD code, generalizing DK4D to nonaxisymmetric geometries, and simulating more quickly-evolving and realistic plasmas. The ultimate goal of this work is to perform self-consistent, hybrid simulations of complex tokamak instabilities and calculations of neoclassical toroidal viscosity. A J. J. Ramos, Phys. Plasmas 17, 082502 (2010). B J. J. Ramos, Phys. Plasmas 18, 102506 (2011).
- Research Organization:
- Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC), Fusion Energy Sciences (FES)
- Contributing Organization:
- Princeton University, NJ (United States). Dept. of Astrophysical Sciences. Program in Plasma Physics
- DOE Contract Number:
- AC02-09CH11466
- OSTI ID:
- 1365824
- Country of Publication:
- United States
- Language:
- English
Similar Records
Final Report for "Tech-X Corporation work for the SciDAC Center for Simulation of RF Wave Interactions with Magnetohydrodynamics (SWIM)"
Fluid equations in the presence of electron cyclotron current drive