Search NASA⌕ Search

SEARCH · Search NASA

Results for “steady state tokamak”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Progress of Divertor Heat and Particle Flux Control in EAST for Advanced Steady-State Operation in the Last 10 Years

Active control of the excessively high heat and particle fluxes on the divertor target plates is of fundamental importance to the steady state operation of tokamaks, especially for fusion reactors. A series of experiments have been carried out on this critical issue to relieve the plasma-wall interactions in the experimental advanced superconducting tokamak (EAST) in the last ten years, not only contributing to the long pulse operation of EAST itself, but also providing physical understandings and potential techniques to the next-generation devices like ITER. Here we have characterized the power deposition pattern and broadened the divertor footprint width effectively. The plasma-wetted area is actively handled using either 3-dimentional edge magnetic topology or advanced plasma equilibrium, thereby peak heat flux around the strike point is reduced. Active control of detachment or radiation compatible with core plasma performance has progressed significantly in very recent years, with a series of active feedback control modules developed and utilized successfully, based on the divertor physics advances with both experiments and simulation. The upper divertor of EAST was upgraded from graphite to active water-cooling ITER-like tungsten in 2014, exhibiting much enhanced heat removal capability. As for the particle exhaust including both fueling and impurity particles, in addition to wall conditioning and impurity source control, the efficiency of particle flux exhaust is optimized by making full use of the divertor closure and the plasma drifts in both scrape-off layer and divertor volume. These heat and particle exhaust advances contribute greatly to a series of EAST achievements like H-mode operation over 100 s. A brief near-term plan on the integrated control of divertor plasma-wall interactions in long-time scale will also be introduced, aiming to provide favorable divertor operation solution for ITER and CFETR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Design considerations for optimizing the transient CHI injector on QUEST

Transient coaxial helicity injection (transient CHI) on the QUEST (Q-shu University Experiment with Steady-State ST) spherical tokamak (ST) has recently validated the floating biased electrode configuration for solenoid-free plasma startup. In support of a significant divertor upgrade on QUEST, the details of the transient CHI injector geometry on QUEST have been examined more closely using tokamak simulation code (TSC) simulations. QUEST uses a HIT-II-like (Helicity Injected Torus) injector configuration in which the injector region is comprised of coaxial injector electrodes located in the lower part of the machine. TSC simulations indicate high injector flux operation may benefit from an increased gap between the coaxial electrodes. This is the first study to examine the differences between an open and closed electrode configuration for solenoid-free plasma startup using transient CHI. Results show that both configurations can generate similar levels of closed flux, but the closed configuration may be easier to implement in some ST designs as a much smaller portion of the injector needs to be close to the injector flux coil. Results show that increasing the electrode gap width from the present 10.8 cm to about 15–20 cm would increase the closed flux fraction by about 40%. The results presented in this paper are generally applicable to the CHI design for other STs.

coaxial helicity injection↗
DOE OSTI2020

Simulation of divertor heat flux width on EAST by BOUT++ transport code

The BOUT++ edge plasma transport code is applied to study the effects of neutral recycling, drifts and heating scheme on edge plasma profiles and the divertor heat flux width of two discharges of the experimental advanced superconducting tokamak (EAST) steady-state H-mode plasmas heated by low hybrid wave (LHW) and neutral beam (NB), respectively. Neutral recycling seems to have played an important role in the plasma density profile. The edge plasma density drops dramatically for the case w/o neutral recycling, while it can be sustained for the case w/ neutral recycling. Drifts are found to have significant influences on edge plasma profiles. Both the amplitude and width of the divertor heat flux are found to have increased a lot due to drifts. The simulated heat flux width w/ drifts for the two discharges shows reasonable agreement with the experiments. However, the width from the simulation and experiment for the LHW heated discharge is much larger than that of the NB heated discharge. Comparison with Goldston's drift-based model and more detailed analysis on the heat flux contributions from drifts versus turbulence show that drifts are the dominant factor in edge plasma transport for both LHW and NB heated discharges, while turbulence may have played a more important role in determining the heat flux width in LHW heated discharges than that in NB heated discharges. This may account for the larger heat flux width in the LHW heated discharge. The magnetic topology and the equilibrium change by the LHW power may also be a potential reason for the larger heat flux width in the LHW heated discharges, which still needs more evidence and studies to support it. Finally, additional SOLPS simulation w/o drifts for the two discharges show reasonable agreement with the counterpart from the BOUT++ simulation, suggesting the two are both suitable codes for EAST edge plasma simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Realization of thousand-second improved confinement plasma with Super I-mode in Tokamak EAST

Mastering nuclear fusion, which is an abundant, safe, and environmentally competitive energy, is a great challenge for humanity. Tokamak represents one of the most promising paths toward controlled fusion. Obtaining a high-performance, steady-state, and long-pulse plasma regime remains a critical issue. Recently, a big breakthrough in steady-state operation was made on the Experimental Advanced Superconducting Tokamak (EAST). A steady-state plasma with a world-record pulse length of 1056 s was obtained, where the density and the divertor peak heat flux were well controlled, with no core impurity accumulation, and a new high-confinement and self-organizing regime (Super I-mode = I-mode + e-ITB) was discovered and demonstrated. These achievements contribute to the integration of fusion plasma technology and physics, which is essential to operate next-step devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Sustained edge-localized-modes suppression and radiative divertor with an impurity-driven instability in tokamak plasmas

We report simultaneous control of the large edge localized modes (ELMs) and divertor heat fluxes in a metal wall environment is a critical issue for steady-state operation of a tokamak fusion reactors. Here we report a sustained ELM suppression scenario achieved in the EAST tokamak compatible with radiative divertor using different seeding impurity species over a wide range of conditions. A low-n mode appears, as manifested by the oscillations of a radiation front near the X-point. This mode appears to drive strong particle transport and tungsten exhaust, which is essential to the maintenance of the ELM-stable state. We have developed a model to explain the mode excitation, by coupling the impurity radiative condensation instability to drift waves, which could explain some characteristics of the low-n mode well. The low-n mode may offer a new ELM-stable scenario compatible with radiative divertor for future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Results and lessons learned from accelerating radio frequency modeling using machine learning [slides]

The “advanced tokamak” reactor concept is a leading candidate for a steady state fusion pilot plant. An advanced tokamak (AT) sustains a majority of the required plasma current with effects resulting from maintenance of the peaked pressure at the device center. This current is augmented by auxiliary current drive sources. These auxiliary actuators may consist of neutral particle beams and/or radio frequency (RF) systems such as lower hybrid current drive (LHCD) and high harmonic fast wave (HHFW) current drive using radio and microwaves from antennas. The primary focus of this work is to develop models of RF current profile control suitable for use in integrated modeling frameworks and for real-time control in experiments. Direct physics models of RF current drive can be computationally intensive. In order to achieve predictive times appropriate for the thousands of calls needed in real-time control of experiments and for use in integrated models, we will apply modern machine learning (ML) techniques to accelerate these models and interpolate their results. To generate the fast and accurate models for use in control level algorithms and integrated modeling we need to replace present models with high dimensional interpolation of their results. We will perform additional simulations across a broader parameter range for EAST and other tokamaks in different physics regimes (Alcator C-Mod, DIII-D, WEST, CFETR, ARC, ITER) and combine them into a larger database for training and testing of the ML models. Further testing of the control level models with experimental current profile data from EAST and C-Mod tokamaks will provide additional confirmation of the control level model before integration in a tokamak control system or integrated modeling suite. ML will be used to optimize the selection of training data consisting of RF current driven at different values of density profile, temperature profile, plasma current, and wavenumber. ML will also be used to facilitate classification of current drive from these input data. The output of this effort will be a validated classifier capable of determining the current drive profiles for HHFW CD and LHCD on a mille-second timescale. This will provide a breakthrough capability enabling real-time control of RF driven current profiles in experiments including ITER ICRF and use integrated modeling frameworks requiring thousands of current profile calculations in discharge simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

To dee or not to dee: costs and benefits of altering the triangularity of a steady-state DEMO-like reactor

Shaping a tokamak plasma to have a negative triangularity may allow operation in an ELM-free L-mode regime and with a larger strike-point radius, ameliorating divertor power-handling requirements. However, the shaping has a potential drawback in the form of a lower no-wall ideal beta limit, found using the MHD codes CHEASE and DCON. Using the new fusion systems code FAROES, we construct a steady-state DEMO2 reactor model. This model is essentially zero-dimensional and neglects variations in physical mechanisms like turbulence, confinement, and radiative power limits, which could have a substantial impact on the conclusions deduced herein. Keeping its shape otherwise constant, we alter the triangularity and compute the effects on the levelized cost of energy (LCOE). If the tokamak is limited to a fixed B field, then unless other means to increase performance (such as reduced turbulence, improved current drive efficiency or higher density operation) can be leveraged, a negative-triangularity reactor is strongly disfavored in the model due to lower \beta_N limits at negative triangularity, which leads to tripling of the LCOE. However, if the reactor is constrained by divertor heat fluxes and not by magnet engineering, then a negative-triangularity reactor with higher B0 could be favorable: we find a class of solutions at negative triangularity with lower peak heat flux and lower LCOE than those of the equivalent positive triangularity reactors.

CHEASE↗
DOE OSTI2022

To dee or not to dee: costs and benefits of altering the triangularity of a steady-state DEMO-like reactor

Shaping a tokamak plasma to have a negative triangularity may allow operation in an edge-localized mode-free L-mode regime and with a larger strike-point radius, ameliorating divertor power-handling requirements. However, the shaping has a potential drawback in the form of a lower no-wall ideal beta limit, found using the MHD codes CHEASE and DCON. Using the new fusion systems code FAROES, we construct a steady-state DEMO2 reactor model. This model is essentially zero-dimensional and neglects variations in physical mechanisms like turbulence, confinement, and radiative power limits, which could have a substantial impact on the conclusions deduced herein. Keeping its shape otherwise constant, we alter the triangularity and compute the effects on the levelized cost of energy (LCOE). If the tokamak is limited to a fixed B field, then unless other means to increase performance (such as reduced turbulence, improved current drive efficiency or higher density operation) can be leveraged, a negative-triangularity reactor is strongly disfavored in the model due to lower β N limits at negative triangularity, which leads to tripling of the LCOE. However, if the reactor is constrained by divertor heat fluxes and not by magnet engineering, then a negative-triangularity reactor with higher B 0 could be favorable: we find a class of solutions at negative triangularity with lower peak heat flux and lower LCOE than those of the equivalent positive triangularity reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Reply to Comment on ‘The advanced tokamak path to a compact net electric fusion pilot plant’

Abstract The Comment by Manheimer (submitted to Nucl. Fusion with this response) on our recent paper Buttery et al (2021 Nucl. Fusion 61 046028), has mischaracterized our paper with a series of misleading statements about its content, which it then seeks to refute. It also offers a series of assertions that are unsupported by refereed publication or scientific data—and sometimes in contradiction to the published record and well-established points in the scientific community. In this response we address the main thrusts and themes of the comment, while providing a point-by-point response in the appendix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Plasma self-driven current in tokamaks with magnetic islands

Magnetic island perturbations may cause a reduction in plasma self-driven current that is needed for tokamak operation. A novel effect on tokamak self-driven current revealed by global gyrokinetic simulations is due to magnetic-island-induced 3D electric potential structures, which have the same dominant mode numbers as that of the magnetic island, whereas centered at both the inner and outer edge of the island. The non-resonant potential islands are shown to drive a current through an efficient nonlinear parallel acceleration of electrons. In large aspect ratio (large-A) tokamak devices, this new effect can result in a significant global reduction of the electron bootstrap current when the island size is sufficiently large, in addition to the local current loss across the island region due to the pressure profile flattening. It is shown that there exists a critical magnetic island width for large-A tokamaks beyond which the electron bootstrap current loss is global and increases rapidly with the island size. As such, this process may introduce a size limit for tolerable magnetic islands in large-A tokamak devices in the context of steady state operation. On the other hand, the current loss caused by magnetic islands in low-A tokamaks such as spherical tokamak (ST) NSTX/U is minor. The reduction of the axisymmetric current by magnetic islands scales with the square of island width. However, the loss of the current is mainly local to the island region, and the pace of current loss as the island size increases is substantially slower compared to large-A tokamaks. In particular, the bootstrap current reduction in STs is even smaller in the reactor-relevant high-β p regime where neoclassical tearing modes are more likely to develop.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Radiative pulsed L-mode operation in ARC-class reactors

Abstract A new ARC-class, highly-radiative, pulsed, L-mode, burning plasma scenario is developed and evaluated as a candidate for future tokamak reactors. Pulsed inductive operation alleviates the stringent current drive requirements of steady-state reactors, and operation in L-mode affords ELM-free access to ∼ 90 % core radiation fractions, significantly reducing the divertor power handling requirements. In this configuration the fusion power density can be maximized despite L-mode confinement by utilizing high-field to increase plasma densities and current. This allows us to obtain high gain in robust scenarios in compact devices with P fus > 1000 MW despite low confinement. We demonstrate the feasibility of such scenarios here; first by showing that they avoid violating 0D tokamak limits, and then by performing self-consistent integrated simulations of flattop operation including neoclassical and turbulent transport, magnetic equilibrium, and radiofrequency current drive models. Finally we examine the potential effect of introducing negative triangularity with a 0D model. Our results show high-field radiative pulsed L-mode scenarios are a promising alternative to the typical steady state advanced tokamak scenarios which have dominated tokamak reactor development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A model investigation of the impact of lower hybrid wave scattering angle on current drive profile in EAST and Alcator C-Mod

Lower hybrid current drive (LHCD) is beneficial for developing a steady-state operation scenario in a tokamak. This paper conducts a modeling investigation to identify an optimum rotation angle of the initial lower hybrid perpendicular (to the background magnetic field $\overrightarrow{B}$) wavevector for best matching the experimental RF current profile. It is hypothesized that central RF power deposition widely observed in the present-day LHCD experiments arises from wave scattering by turbulence. In a standard model without considering such interactions, the predicted power deposition profile is generally broad with off-axis peaking, not in agreement with experimental observations. A heuristic approach is adopted by introducing a spectral broadening mechanism by modifying the initial orientation of the perpendicular wavevector. The ray-tracing/Fokker–Planck solver GENRAY/CQL3D is utilized within the python-based π-scope framework. A focus is given to identify the perpendicular wavenumber orientation angle with respect to the magnetic surface normal vector at the initial ray location. Our modeling study shows that rotating the perpendicular wavevector in such a way as to increase the initial poloidal component is effective in reproducing the centrally peaked current profile observed in normal shear plasmas on both EAST and C-Mod. These waves can readily be absorbed to the central plasma, which reduces the sensitivity of the power deposition profile to a slight change of the plasma condition. The same approach is also found to help broaden the off-axis power deposition profile in a reverse-shear EAST plasma, leading to a better agreement with the experiment. The results presented here suggest that spectral modification arising from edge density fluctuations in a tokamak may need to be considered in understanding wave propagation and absorption. A further experimental and theoretical/modeling study is vital as a reverse approach is adopted in this study. Finally, our work suggests that mitigation or control measures are critical for parasitic effects occurring on the first pass in a reactor regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Observation of a new type of self-generated current in magnetized plasmas

A tokamak, a torus-shaped nuclear fusion device, needs an electric current in the plasma to produce magnetic field in the poloidal direction for confining fusion plasmas. Plasma current is conventionally generated by electromagnetic induction. However, for a steady-state fusion reactor, minimizing the inductive current is essential to extend the tokamak operating duration. Several non-inductive current drive schemes have been developed for steady-state operations such as radio-frequency waves and neutral beams. However, commercial reactors require minimal use of these external sources to maximize the fusion gain, Q, the ratio of the fusion power to the external power. Apart from these external current drives, a self-generated current, so-called bootstrap current, was predicted theoretically and demonstrated experimentally. Here, we reveal another self-generated current that can exist in a tokamak and this has not yet been discussed by present theories. We report conclusive experimental evidence of this self-generated current observed in the KSTAR tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Advances in physics understanding of high poloidal beta regime toward steady-state operation of CFETR

Experimental and modeling investigations of the high βp scenarios on the DIII-D and EAST tokamaks show advantages in high energy confinement, avoidance of n = 1 MHD, and core-edge integration with reduced heat flux, making this scenario an attractive option for CFETR steady-state operation. Experiments show that plasmas with high confinement and high density can be achieved with neutral beam injection on DIII-D (β p ~ 2.2, β N ~ 3.5, f BS ~ 50%, f Gw ~ 1.0, H 98y2 ~ 1.5) and pure RF power on EAST (β P ~ 2.0, β N ~ 1.6, f BS ~ 50%, f Gw ~ 0.8, H 98y2 > 1.3). By tailoring the current density profile, a q-profile with local (off-axis) negative shear is achieved that yields improved confinement and MHD stability. Transport analysis and simulation suggest that the combination of high density gradient and high Shafranov shift allows turbulence stabilization and higher confinement. Using on-axis ECH injection, tungsten accumulation is avoided on EAST, and this is reproduced in modeling. Reduced heat flux (by >40%) and maintenance of high core confinement is achieved with active feedback control of the radiated divertor, an important result for long pulse operation in tokamaks. In conclusion, the improved physics understanding and validated modeling tools are used to design a 1GW steady-state scenario for CFETR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Multiphysics simulations of a steady-state lower hybrid current drive antenna for the FSNF

The Fusion Nuclear Science Facility (FNSF) is a proposed tokamak reactor with the mission to investigate operation of a fusion reactor in a nuclear environment. The high neutron fluence component of the FNSF mission requires steady-state operation for extremely long pulses (t_{pulse} ∼ months) at full power. Plasma sustainment and current drive will be critical components of a successful FNSF. COMSOL Multiphysics software is used for combined radiofrequency (RF) and thermal simulations of the lower hybrid current drive (LHCD) antenna system. These simu- lations consider the resistive RF losses in the antenna including realistic surface roughness and a range of potential materials. The thermal analysis adds volumetric nuclear heating, plasma heat flux on leading edges, and electromagnetic radiation from the plasma to the RF heating calculated by COMSOL. Additional neutronics calculations have been performed to determine the impact of these antenna designs on activated waste disposal for the materials considered. The simulations show that it is technically feasible to implement a fully-active multi-junction (FAM) rather than a passive-active multi-junction (PAM) style of antenna if the septum between adjacent waveguides is sufficiently wide and the thermal conductivity of the structural material is sufficiently high. The FAM has the benefit of higher achievable power density with respect to the PAM, which results in a more compact antenna with potentially lower impact on neutron shielding and tritium breeding. These considerations point to tungsten rather than steel as the preferred structural material in constructing the antenna.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An Overview of the Hybrid Illinois Device for Research and Applications Material Analysis Test-stand (HIDRA-MAT)

The Hybrid Illinois Device for Research and Applications (HIDRA) at the University of Illinois at Urbana-Champaign is a toroidal plasma device that enables fusion plasma-material interaction testing with both stellarator and tokamak plasmas. HIDRA’s long-pulse steady state stellarator plasmas provide a testbed for plasma facing component (PFC) plasma exposures. The HIDRA Material Analysis Test-stand (HIDRA-MAT) is a material characterization module attached to HIDRA that is being designed and fabricated to include thermal desorption spectroscopy and laser induced breakdown spectroscopy systems for in-vacuo PFC characterization. A specialized rotatable sample holder positions the sample for liquid metal droplet application from a liquid metal droplet injector on HIDRA-MAT. Early experiments look to investigate the effect liquid lithium has on porous tungsten samples’ retention of H, D, and He after plasma exposure. Preliminary results from a dual residual gas analyzer system show the ability to differentiate D2 and He in HIDRA-MAT. Finally, this work aims to advance the understanding of liquid metal PFCs and further the design and development of new fusion PFCs and technologies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Benchmarking NIMROD continuum kinetic formulations through the steady-state poloidal flow

In this work, continuum kinetic formulations are employed as a mechanism to include closure physics in an extended magnetohydrodynamics model. Two continuum kinetic approaches have been implemented in the plasma fluid code NIMROD [Sovinec et al., “Nonlinear magnetohydrodynamics with high-order finite elements,” J. Comput. Phys. 195, 355 (2004)] including a Chapman–Enskog-like (CEL) formulation and a more conventional df approach. Ion kinetic closure schemes are employed to describe the neoclassical flow properties in axisymmetric toroidal geometry. In particular, predictions for steady-state values of poloidal flow profiles in tokamak geometry are provided using both the df formulation and two different solution techniques for the CEL approach. These results are benchmarked against analytic theory predictions as well as results from the drift kinetic code DK4D. The continuum kinetic formulations employed here show agreement with both the analytic theory and DK4D results, and offer a novel velocity space representation involving higher-order finite elements in pitch angle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗