Engineering topics

Holland, C. (ORCID:0000000160292306)

Publications and source records attributed to Holland, C. (ORCID:0000000160292306).

5 records

DOE OSTI2024

Understanding the L-H isotope effect at the DIII-D tokamak and advancements in synthetic turbulence diagnostics

Abstract It is determined that while heat flux differences between hydrogen and deuterium isotope experiments result from natural differences in carbon impurity content at DIII-D, it is not the origin of the low to high confinement mode (L-H) transition isotope effect. More specifically, a two times larger edge radial electric field in hydrogen compared to deuterium is uncovered and believed to play an important role. The origin of this radial electric field difference is determined to have two possible origins: differences in poloidal rotation and turbulent Reynolds stress in the closed field line region, and increased outer strike point temperatures and space potentials on open field lines. Experimental observations from both profile and turbulence diagnostics are supported by nonlinear gyrokinetic simulations using the code CGYRO. Simulations illustrated heat transport isotope effects in the plasma edge and shear layer resulting from differences in impurity content, electron non-adiabaticity, and main ion mass dependent E × B shear stabilization. Turbulence prediction comparisons from flux-matched CGYRO simulations to experimental measurements including electron temperature, density and velocity fluctuations are found to be in good agreement with available data. A dedicated DIII-D experiment in hydrogen was performed to seed more carbon than naturally occurring, to match deuterium experiments, and possibly reduce the L-H power threshold based on gyro-kinetic predictions. To our surprise, while ion temperature gradient (ITG) turbulence was stabilized, nodiscernible change in L-H power threshold were observed in these special hydrogen experiments. In particular, it is noticed that the edge radial electric field and Reynolds stress were observed as nearly unchanging in the presence of ITG stabilization. These experimental data have enabled a more comprehensive picture of the multitude of isotope effects at play in fusion experiments, and the important potential connection between the confined and unconfined plasma regions in regulating L-H transition dynamics.

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

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

On the origin of the DIII-D L-H power threshold isotope effect

Abstract The increased low to high confinement mode (L to H-mode) power threshold P L H in DIII-D low collisionality hydrogen plasmas (compared to deuterium) is shown to result from lower impurity (carbon) content, consistent with reduced (mass-dependent) physical and chemical sputtering of graphite. Trapped gyro-Landau fluid (TGLF) quasilinear calculations and local non-linear gyrokinetic CGYRO simulations confirm stabilization of ion temperature gradient (ITG) driven turbulence by increased carbon ion dilution as the most important isotope effect. In the plasma edge, electron non-adiabaticity is also predicted to contribute to the isotope dependence of thermal transport and P L H , however its effect is subdominant compared to changes from impurity isotopic behavior. This L-H power threshold reduction with increasing carbon content at low collisionality is in stark contrast to high collisionality results, where additional impurity content appears to increase the power necessary for H-mode access.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Development of compact tokamak fusion reactor use cases to inform future transport studies

The OMFIT STEP (Meneghini et al. , Nucl. Fusion , vol. 10, 2020, p. 1088) workflow has been used to develop inductive and steady-state H-mode core plasma scenario use cases for a $B_0 = 8 \, {\rm T}$ , $R_0 = 4 \, {\rm m}$ machine to help guide and inform future higher-fidelity studies of core transport and confinement in compact tokamak reactors. Both use cases are designed to produce 200 MW or more of net electric power in an up-down symmetric plasma with minor radius $a = 1.4 \, {\rm m}$ , elongation $\kappa = 2.0$ , triangularity $\delta = 0.5$ and effective charge $Z_{{\rm eff}} \simeq 2$ . Additional considerations based on the need for compatibility of the core with reactor-relevant power exhaust solutions and external actuators were used to guide and constrain the use case development. An extensive characterization of core transport in both scenarios is presented, the most important feature of which is the extreme sensitivity of the results to the quantitative stiffness level of the transport model used as well as the predicted critical gradients. This sensitivity is shown to arise from different levels of transport stiffness exhibited by the models, combined with the gyroBohm-normalized fluxes of the predictions being an order of magnitude larger than other H-mode plasmas. Additionally, it is shown that although heating in both plasmas is predominantly to the electrons and collisionality is low, the plasmas remain sufficiently well coupled for the ions to carry a significant fraction of the thermal transport. As neoclassical transport is negligible in these conditions, this situation inherently requires long-wavelength ion gyroradius-scale turbulence to be the dominant transport mechanism in both plasmas. These results are combined with other basic considerations to propose a simple heuristic model of transport in reactor-relevant plasmas, along with simple metrics to quantify coupling and core transport properties across burning and non-burning plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗