Search NASA⌕ Search

Engineering topics

Candy, J. (ORCID:0000000338846485)

Publications and source records attributed to Candy, J. (ORCID:0000000338846485).

11 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Spherical tokamak physics research in preparation for the operation of NSTX-U

The National Spherical Torus Experiment Upgrade (NSTX-U) is preparing to resume operation, representing a crucial step toward realizing compact, cost-effective fusion pilot plants. In advance of this, extensive modeling and data analysis have been conducted to advance the physics basis for low-aspect-ratio, high-performance plasma regimes, focusing on three core objectives: confinement and stability, power and particle handling, and steady-state operation. Significant progress has been made in understanding the electron temperature flattening in high-β plasmas, which is shown to be driven by a complex interplay of magnetohydrodynamic instabilities (e.g. non-resonant infernal modes), fast-ion-driven Alfvén eigenmodes, and electron and ion-scale micro-instabilities, particularly Kinetic Ballooning Modes (KBMs), whose destabilization is strongly dependent on parallel magnetic field fluctuations (δB ∥ ). Furthermore, a new gyrokinetic critical pedestal model was developed, accurately predicting pedestal structure by identifying KBMs as the primary stability limit, offering a critical constraint for future high-confinement scenarios. To address the challenge of high heat flux, novel liquid lithium plasma-facing components were modeled. The analysis confirmed that lithium vapor shielding is a self-regulating mechanism for heat mitigation, while also emphasizing that strong main ion parallel flow is essential to minimize core lithium contamination. Finally, progress toward steady-state operation was anchored by developing the required physics basis and control tools. This includes predictive modeling for reversed magnetic shear sustainment, demonstrating that magnetic island-induced bootstrap current reduction is negligible in STs, and advancing real-time control and disruption avoidance capabilities. The development of high-speed surrogate models (e.g. MMMNet) provides computationally efficient tools vital for non-inductive scenario optimization and integrated, low-disruptivity operations planned for NSTX-U.

NSTX-U↗
DOE OSTI2024

NSTX-U research advancing the physics of spherical tokamaks

Abstract The objectives of NSTX-U research are to reinforce the advantages of STs while addressing the challenges. To extend confinement physics of low- A , high beta plasmas to lower collisionality levels, understanding of the transport mechanisms that set confinement performance and pedestal profiles is being advanced through gyrokinetic simulations, reduced model development, and comparison to NSTX experiment, as well as improved simulation of RF heating. To develop stable non-inductive scenarios needed for steady-state operation, various performance-limiting modes of instability were studied, including MHD, tearing modes, and energetic particle instabilities. Predictive tools were developed, covering disruptions, runaway electrons, equilibrium reconstruction, and control tools. To develop power and particle handling techniques to optimize plasma exhaust in high performance scenarios, innovative lithium-based solutions are being developed to handle the very high heat flux levels that the increased heating power and compact geometry of NSTX-U will produce, and will be seen in future STs. Predictive capabilities accounting for plasma phenomena, like edge harmonic oscillations, ELMs, and blobs, are being tested and improved. In these ways, NSTX-U researchers are advancing the physics understanding of ST plasmas to maximize the benefit that will be gained from further NSTX-U experiments and to increase confidence in projections to future devices.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Nonlinear gyrokinetic modelling of high confinement negative triangularity plasmas

Abstract Nonlinear gyrokinetic simulations correctly predict particle as well as ion and electron energy fluxes of high confinement plasmas with a negative triangularity cross sectional shape, showing that core transport in these plasmas is well described by standard gyrokinetic models. Experimentally inferred power balance fluxes are mostly reproduced within one standard deviation across a wide portion of the minor radius. Experimental conditions are reproduced by ion scale simulations, without the need to include density and temperature profile curvature effects. The experimental case is used as baseline to predict that the non-dimensional confinement scaling in negative triangularity plasmas increases strongly with plasma current while slightly degrading at increasing normalized pressure and decreasing collisionality. Recent experiments showed that low toroidal rotation negatively impacts confinement; consistent with the experiment, simulations predict that low rotational shear significantly affects confinement unless the plasma effective charge is maintained above a minimum level. Core confinement is predicted to significantly degrade in low aspect ratio devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Enhancing predictive capabilities in fusion burning plasmas through surrogate-based optimization in core transport solvers

Abstract This work presents the PORTALS framework (Rodriguez-Fernandez et al 2022 Nucl. Fusion 62 076036), which leverages surrogate modeling and optimization techniques to enable the prediction of core plasma profiles and performance with nonlinear gyrokinetic simulations at significantly reduced cost, with no loss of accuracy. The efficiency of PORTALS is benchmarked against standard methods, and its full potential is demonstrated on a unique, simultaneous 5-channel (electron temperature, ion temperature, electron density, impurity density and angular rotation) prediction of steady-state profiles in a DIII-D ITER Similar Shape plasma with GPU-accelerated, nonlinear CGYRO (Candy et al 2016 J. Comput. Phys. 324 73–93). This paper also provides general guidelines for accurate performance predictions in burning plasmas and the impact of transport modeling in fusion pilot plants studies.

Physics↗
DOE OSTI2024

Core performance predictions in projected SPARC first-campaign plasmas with nonlinear CGYRO

This work characterizes the core transport physics of SPARC early-campaign plasmas using the PORTALS-CGYRO framework. Empirical modeling of SPARC plasmas with L-mode confinement indicates an ample window of breakeven (Q > 1) without the need of H-mode operation. Extensive modeling of multi-channel (electron energy, ion energy, and electron particle) flux-matched conditions with the nonlinear CGYRO code for turbulent transport coupled to the macroscopic plasma evolution using PORTALS reveals that the maximum fusion performance to be attained will be highly dependent on the near-edge pressure. Stiff core transport conditions are found, particularly when fusion gain approaches unity, and predicted density peaking is found to be in line with empirical databases of particle source-free H-modes. Impurity optimization is identified as a potential avenue to increase fusion performance while enabling core-edge integration. Extensive validation of the quasilinear TGLF model builds confidence in reduced-model predictions. The implications of projecting L-mode performance to high-performance and burning-plasma devices is discussed, together with the importance of predicting edge conditions.

Rodriguez-Fernandez, P. (ORCID:0000000273611131)↗
DOE OSTI2024

Simultaneous reproduction of experimental profiles, fluxes, transport coefficients, and turbulence characteristics via nonlinear gyrokinetic profile predictions in a DIII-D ITER similar shape plasma

Experimental conditions obtained on the DIII-D tokamak in the ITER Similar Shape (ISS) have been compared extensively with nonlinear gyrokinetic simulation using the CGYRO code [Candy et al., J. Comput. Phys. 324, 73–93 (2016)] with comparisons spanning ion and electron heat fluxes, electron and impurity particle transport, and turbulent fluctuation levels and characteristics. Bayesian optimization techniques [Rodriguez-Fernandez et al., Nucl. Fusion 62(7), 076036 (2022)], combined with nonlinear gyrokinetics, have been used to obtain simultaneously Qi, Qe, and Γe flux-matched profiles that are found to be in good agreement with experimental profile measurements. Synthetic diagnostics were used to compare measured beam emission spectroscopy and correlation electron cyclotron emission turbulent fluctuations with nonlinear simulation. Although some disagreements exist, nonlinear simulations are found to be in generally good agreement with measured fluctuation levels, spectral shapes, and measured radial trends in low-k δne/ne and δTe/Te. Low (Li and C) and mid-Z (Ca) impurity transport was also compared with these flux-matched simulations. Fully stripped, low-Z impurities are well reproduced by the gyrokinetic modeling while clear disagreement exists in comparisons with mid-Z impurities. Nonlinear gyrokinetic investigation into the Z dependence of impurity transport in the ISS conditions is also performed, demonstrating clear trends of impurity diffusion with impurity Z (both D∝Z and D∝1/Z) that vary with the radial location studied. These trends are shown to result from the local dominance of ion temperature gradient or ∇n driven trapped electron mode turbulence and may contribute to the disagreement between simulation and experiment in mid-Z impurity transport. The results of this work represent one of the most complete validation studies of the gyrokinetic model performed to date and provide an example of new capabilities for predicting performance in future fusion devices.

Physics↗
DOE OSTI2023

Flexible, integrated modeling of tokamak stability, transport, equilibrium, and pedestal physics

The STEP (Stability, Transport, Equilibrium, and Pedestal) integrated-modeling tool has been developed in OMFIT to predict stable, tokamak equilibria self-consistently with core-transport and pedestal calculations. STEP couples theory-based codes to integrate a variety of physics, including magnetohydrodynamic stability, transport, equilibrium, pedestal formation, and current-drive, heating, and fueling. The input/output of each code is interfaced with a centralized ITER-Integrated Modelling & Analysis Suite data structure, allowing codes to be run in any order and enabling open-loop, feedback, and optimization workflows. This paradigm simplifies the integration of new codes, making STEP highly extensible. STEP has been verified against a published benchmark of six different integrated models. Core-pedestal calculations with STEP have been successfully validated against individual DIII-D H-mode discharges and across more than 500 discharges of the H98,y2 database, with a mean error in confinement time from experiment less than 19%. STEP has also reproduced results in less conventional DIII-D scenarios, including negative-central-shear and negative-triangularity plasmas. Predictive STEP modeling has been used to assess performance in several tokamak reactors. Simulations of a high-field, large-aspect-ratio reactor show significantly lower fusion power than predicted by a zero-dimensional study, demonstrating the limitations of scaling-law extrapolations. STEP predictions have found promising scenarios for an EXhaust and Confinement Integration Tokamak Experiment, including a high-pressure, 80%-bootstrap-fraction plasma. ITER modeling with STEP has shown that pellet fueling enhances fusion gain in both the baseline and advanced-inductive scenarios. Finally, STEP predictions for the SPARC baseline scenario are in good agreement with published results from the physics basis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Gyrokinetic analysis of inter-edge localized mode transport mechanisms in a DIII-D pedestal

In this study, gyrokinetic simulations are used to study pedestal fluctuations for DIII-D discharge 174082 using the GENE code. Nonlinear local simulations indicate that electron heat flux has contributions from electron temperature gradient-driven transport but at levels insufficient to satisfy power balance. We show that microtearing modes (MTM) and neoclassical transport are likely to account for the remaining observed energy losses in the electron and ion channels, respectively. The MTM instabilities found in the simulations are consistent with the high-frequency fluctuations identified in the magnetic fluctuation data from Mirnov coils. The fluctuation data in this discharge also exhibit a low-frequency band of fluctuations. By modifying the equilibrium profiles and plasma β, simulations produce MHD modes, which may be responsible for these observed low-frequency fluctuations. We compare several metrics involving ratios of fluctuation amplitudes and transport quantities for both MTMs and MHD modes. This analysis suggests that the available data are consistent with the simultaneous activity of both MHD modes and MTMs provided that the former is limited largely to the particle transport channel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Equilibrium reconstruction of DIII-D plasmas using predictive modeling of the pressure profile

New workflows have been developed for predictive modeling of magnetohydrodynamic (MHD) equilibrium in tokamak plasmas. The goal of this work is to predict the MHD equilibrium in tokamak discharges without having measurements of the kinetic profiles. The workflows include a cold start tool, which constructs all the profiles and power flows needed by transport codes; a Grad–Shafranov equilibrium solver; and various codes for the sources and sinks. For validation purposes, a database of DIII-D tokamak discharges has been constructed that is comprised of scans in the plasma current, toroidal magnetic field, and triangularity. Initial efforts focused on developing a workflow utilizing an empirically derived pressure model tuned to DIII-D discharges with monotonic safety factor profiles. This workflow shows good agreement with experimental kinetic equilibrium calculations, but is limited in that it is a single fluid (equal ion and electron temperatures) model and lacks H-mode pedestal predictions. The best agreement with the H-mode database is obtained using a theory-based workflow utilizing pressure profile predictions from a coupled TGLF turbulent transport and EPED pedestal models together with external magnetics and Motional Stark Effect (MSE) data to construct the equilibrium. Here, we obtain an average root mean square error of 5.1% in the safety factor profile when comparing the predicted and experimental kinetic equilibrium. We also find good agreement with the plasma stored energy, internal inductance, and pressure profiles. Including MSE data in the theory-based workflow results in noticeably improved agreement with the q-profiles in high triangularity discharges in comparison with the results obtained with magnetic data only. The predictive equilibrium workflow is expected to have wide applications in experimental planning, between-shot analysis, and reactor studies.

Kinsey, J. E. (ORCID:0000000193347473)↗
DOE OSTI2021

Gyrokinetic simulation of turbulence and transport in the SPARC tokamak

The turbulence and transport expected in the SPARC tokamak Primary Reference Discharge (PRD) [P. Rodriguez-Fernandez et al., J. Plasma Phys. 86, 865860503 (2020)] have been investigated with the gyrokinetic code CGYRO [J. Candy et al., J. Comput. Phys. 324, 73–93 (2016)]. Linear and nonlinear simulations that focus on ion (kθρs<1.0) and electron-scale (kθρs>1.0) turbulence were used to probe the nature of the turbulence and the resulting transport in the fusion core. It is found that in the SPARC PRD, ion temperature gradient (ITG) turbulence is expected to dominate transport over most of the profile with some potential trapped electron mode impact in the near edge. Stiff turbulence is observed over a part of the plasma core such that SPARC's ion temperature profile will likely be pinned to just above the critical gradient for ITG. The role of electromagnetic turbulence, rotation, and electron-scale turbulence was investigated to provide some insight into the physics required to accurately predict SPARC performance via gyrokinetics. Additionally, predictions of impurity peaking for potential low- and high-Z SPARC first-wall materials are probed using ion-scale simulation. The dominance of low-k turbulence in SPARC provides a potential opportunity for more tractable prediction of plasma profiles using nonlinear gyrokinetics. This work is the first step toward full gyrokinetic profile prediction of SPARC kinetic profiles and the resulting fusion power and plasma gain.

Howard, N. T. (ORCID:0000000287876309)↗