Engineering topics

Dorf, M. A.

Publications and source records attributed to Dorf, M. A..

11 records

DOE OSTI2021

Gyrokinetic and extended-MHD simulations of a flow shear stabilized Z-pinch experiment

Axisymmetric (m = 0) gyrokinetic and extended-MHD simulations of a sheared-flow Z-pinch plasma are performed with the high-order finite volume code COGENT. The present gyrokinetic model solves the long-wavelength limit of the gyrokinetic equation for both ion and electron species coupled to the electrostatic gyro-Poisson equation for the electrostatic potential. The extended-MHD model is electromagnetic and includes the effects of the gyro-viscous pressure tensor, diamagnetic electron and ion heat fluxes, and generalized Ohm's law. Additionally, a prominent feature of this work is that the radial profiles for the plasma density and temperature are taken from the fusion Z-pinch experiment (FuZE), and the magnetic field profile is obtained as a solution of the MHD force balance equation. Such an approach allows to address realistic plasma parameters and provide insights into the current and planned experiments. In particular, it is demonstrated that the radial profiles play an important role in stabilization, as the embedded guiding center (E×B ) drift has a strong radial shear, which can contribute to the Z-pinch stabilization even in the absence of the fluid flow shear. The results of simulations for the FuZE plasma parameters show a decrease in the linear growth rate with an increase in the flow shear; however, full stabilization in the linear regime is not observed even for large (comparable to the Alfvén velocity) radial variations of the axial flow. Nonlinear stability properties of the FuZE plasmas are also studied, and it is found that profile broadening can have a pronounced stabilizing effect in the nonlinear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Simulation of edge localized mode heat pulse using drift-kinetic ions and Boltzmann electrons

This article reports on 1D+2V heat pulse propagation studies using the COGENT guiding center kinetic code. The model uses magnetized kinetic ions and a simple Boltzmann electron model. Results agree with previous kinetic and fluid modeling benchmark studies that correspond to the parameters of edge localized modes (ELMs) observed on the JET tokamak. The plasma parameters for the edge pedestal and ensuing ELM dynamics are in the low collisionality regime. Hence, the dominant balance between the assumed Maxwellian ELM source and collisionless parallel advection causes the ion PDF to develop a significantly anisotropic velocity distribution. Adding nonlinear Coulomb ion-ion collisions to the model acts to smooth the sharp features of the ion distribution function, but the anisotropy remains robust due to the low collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2018

Verification of 5D continuum gyrokinetic code COGENT: Studies of kinetic drift wave instability

COGENT (Continuum Gyrokinetic Edge New Technology) is a kinetic plasma simulation code that is being developed by the edge simulation laboratory (ESL) collaboration. The original version of the code has been developed in a 4D phase space (2D configuration space and 2D velocity space) to address kinetic plasma phenomena in a complex magnetic field geometry including the core, magnetic separatrix, and scrape-off layer regions. This paper is focused on extending the original 4D phase space to a 5D phase space (3D2V) to address full kinetic turbulences in a tokamak edge region. Here, we report on the current status of 5D COGENT, which presently operates in a shear-less, simple slab geometry. As a verification study, we use the problems of collision-less drift wave instability (universal instability) including its modification in the presence of collisional effects. The simulation model includes the gyrokinetic equations for ion and electron species coupled to the long-wavelength limit of the gyro-Poisson equation. Collisional effects are represented by the Krook collision model. Linear analytical results for the drift mode growth rate and real frequency are recovered, and the non-linear stage is modelled and analysed as well. In addition, extensive 5D runs have been performed to address the effects of the drift wave instability on blob/filamentary structures characteristic of a tokamak edge. Finally, a helical-shaped potential perturbation is observed to grow exponentially in time while spinning around the filament axis with electron drift frequency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

Continuum kinetic modeling of the tokamak plasma edge

The first 4D (axisymmetric) high-order continuum gyrokinetic transport simulations that span the magnetic separatrix of a tokamak are presented. The modeling is performed with the COGENT code, which is distinguished by fourth-order finite-volume discretization combined with mapped multiblock grid technology to handle the strong anisotropy of plasma transport and the complex X-point divertor geometry with high accuracy. The calculations take into account the effects of fully nonlinear Fokker-Plank collisions, electrostatic potential variations, and anomalous radial transport. Topics discussed include: (a) ion orbit loss and the associated toroidal rotation and (b) edge plasma relaxation in the presence of anomalous radial transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

Evolution of helical perturbations in a thin-shell model of an imploding liner

A thin-shell model of the liner stability has been revisited and applied to the stability of the helical perturbations. Several stages of the implosion have been identified, starting from a long initial “latent” phase of an almost resting liner, continuing to the second stage of a rapid contraction and significant perturbation growth, and then transitioning to the third stage where perturbations become ballistic and highly non-linear. The stage of stagnation and rebound is beyond the scope of this paper. An importance of vorticity conservation during the late stages is emphasized. Nonlinear evolution of perturbations is followed up to the point of the formation of cusp structures. Effects of in-surface flows and of their enhancement due to the vorticity conservation are discussed. It is shown that the pre-machined perturbations created only on the outer surface of the liner grow much slower than one could anticipate. The limitations on the thin-shell description are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

Progress with the COGENT Edge Kinetic Code: Implementing the Fokker-Plank Collision Operator

Here, COGENT is a continuum gyrokinetic code for edge plasma simulations being developed by the Edge Simulation Laboratory collaboration. The code is distinguished by application of a fourth-order finite-volume (conservative) discretization, and mapped multiblock grid technology to handle the geometric complexity of the tokamak edge. The distribution function F is discretized in v∥ – μ (parallel velocity – magnetic moment) velocity coordinates, and the code presently solves an axisymmetric full-f gyro-kinetic equation coupled to the long-wavelength limit of the gyro-Poisson equation. COGENT capabilities are extended by implementing the fully nonlinear Fokker-Plank operator to model Coulomb collisions in magnetized edge plasmas. The corresponding Rosenbluth potentials are computed by making use of a finite-difference scheme and multipole-expansion boundary conditions. Details of the numerical algorithms and results of the initial verification studies are discussed. (© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

On the applicability of the standard approaches for evaluating a neoclassical radial electric field in a tokamak edge region

The use of the standard approaches for evaluating a neoclassical radial electric field E r , i.e., the Ampere (or gyro-Poisson) equation, requires accurate calculation of the difference between the gyroaveraged electron and ion particle fluxes (or densities). In the core of a tokamak, the nontrivial difference appears only in high-order corrections to a local Maxwellian distribution due to the intrinsic ambipolarity of particle transport. The evaluation of such high-order corrections may be inconsistent with the accuracy of the standard long wavelength gyrokinetic equation (GKE), thus imposing limitations on the applicability of the standard approaches. However, in the edge of a tokamak, charge-exchange collisions with neutrals and prompt ion orbit losses can drive non-intrinsically ambipolar particle fluxes for which a nontrivial (E r -dependent) difference between the electron and ion fluxes appears already in a low order and can be accurately predicted by the long wavelength GKE. As a result, the parameter regimes where the radial electric field dynamics in the tokamak edge region is dominated by the non-intrinsically ambipolar processes, thus allowing for the use of the standard approaches, are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

Generation of multi-charged high current ion beams using the SMIS 37 gas-dynamic electron cyclotron resonance (ECR) ion source

A gas-dynamic ECR ion source (GaDIS) is distinguished by its ability to produce high current and high brightness beams of moderately charged ions. Contrary to a classical ECR ion source where the plasma confinement is determined by the slow electron scattering into an empty loss-cone, the higher density and lower electron temperature in a GaDIS plasma lead to an isotropic electron distribution with the confinement time determined by the prompt gas-dynamic flow losses. As a result, much higher ion fluxes are available, however a decrease in the confinement time of the GaDIS plasma lowers the ion charge state. The gas-dynamic ECR ion source concept has been successfully realized in the SMIS 37 experimental facility operated at the Institute of Applied Physics, Russia. The use of high-power (~100 kW) microwave (37.5 GHz) radiation provides a dense plasma (~10 13 cm -3 ) with a relatively low electron temperature (~50- 100 eV) and allows for the generation of high current (~1 A/cm 2 ) beams of multi-charged ions. In this work we report on the present status of the SMIS 37 ion source and discuss the advanced numerical modeling of ion beam extraction using the particle-in-cell code WARP

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

Numerical modelling of geodesic acoustic mode relaxation in a tokamak edge

Here, the edge of a tokamak in a high confinement (H mode) regime is characterized by steep density gradients and a large radial electric field. Recent analytical studies demonstrated that the presence of a strong radial electric field consistent with a subsonic pedestal equilibrium modifies the conventional results of the neoclassical formalism developed for the core region. In the present work we make use of the recently developed gyrokinetic code COGENT to numerically investigate neoclassical transport in a tokamak edge including the effects of a strong radial electric field. The results of numerical simulations are found to be in good qualitative agreement with the theoretical predictions and the quantitative discrepancy is discussed. In addition, the present work investigates the effects of a strong radial electric field on the relaxation of geodesic acoustic modes (GAMs) in a tokamak edge. Numerical simulations demonstrate that the presence of a strong radial electric field characteristic of a tokamak pedestal can enhance the GAM decay rate, and heuristic arguments elucidating this finding are provided.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

Simulation of neoclassical transport with the continuum gyrokinetic code COGENT

The development of the continuum gyrokinetic code COGENT for edge plasma simulations is reported. The present version of the code models a nonlinear axisymmetric 4D (R, v∥, μ) gyrokinetic equation coupled to the long-wavelength limit of the gyro-Poisson equation. Here, R is the particle gyrocenter coordinate in the poloidal plane, and v∥ and μ are the guiding center velocity parallel to the magnetic field and the magnetic moment, respectively. The COGENT code utilizes a fourth-order finite-volume (conservative) discretization combined with arbitrary mapped multiblock grid technology (nearly field-aligned on blocks) to handle the complexity of tokamak divertor geometry with high accuracy. Furthermore, topics presented are the implementation of increasingly detailed model collision operators, and the results of neoclassical transport simulations including the effects of a strong radial electric field characteristic of a tokamak pedestal under H-mode conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2012

Progress with the COGENT Edge Kinetic Code: Collision operator options

In this study, COGENT is a continuum gyrokinetic code for edge plasmas being developed by the Edge Simulation Laboratory collaboration. The code is distinguished by application of the fourth order conservative discretization, and mapped multiblock grid technology to handle the geometric complexity of the tokamak edge. It is written in v∥-μ (parallel velocity – magnetic moment) velocity coordinates, and making use of the gyrokinetic Poisson equation for the calculation of a self-consistent electric potential. In the present manuscript we report on the implementation and initial testing of a succession of increasingly detailed collision operator options, including a simple drag-diffusion operator in the parallel velocity space, Lorentz collisions, and a linearized model Fokker-Planck collision operator conserving momentum and energy (© 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗