Search NASA⌕ Search

SEARCH · Search NASA

Results for “ITER baseline scenario”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The root cause of disruptive NTMs and paths to stable operation in DIII-D ITER baseline scenario plasmas

Analyses of the DIII-D ITER Baseline Scenario database support that the disruptive m,n=2,1 magnetic islands are pressure gradient driven, non-linear instabilities seeded in a sequence of stochastic transient magnetic perturbations, and that the current profile relaxation does not affect the m,n=2,1 island onset rate. At low torque, these Neoclassical Tearing Modes are most commonly seeded by non-linear 3-wave coupling when the differential rotation between the q=1 & q=2 rational surfaces approaches zero. Lack of statistically significant difference between the current profiles of stable and unstable states, as well as lack of correlation between the tearing mode onset rate and the current profile relaxation both reject causality between the current profile evolution and the 2,1 magnetic island onsets in these plasmas. These support that preserving the differential rotation between the q=1 and q=2 rational surfaces is key to long pulse stable operation in the plasma scenario planned for ITER, while optimization of the current profile within the explored parameter space may lead to much weaker improvements than sustaining the differential rotation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Multi-scale transport in the DIII-D ITER baseline scenario with direct electron heating and projection to ITER

Multi-scale fluctuations measured by turbulence diagnostics spanning long and short wavelength spatial scales impact energy confinement and the scale-lengths of plasma kinetic profiles in the DIII-D ITER baseline scenario with direct electron heating. Contrasting discharge phases with ECH + neutral beam injection (NBI) and NBI only at similar rotation reveal higher energy confinement and lower fluctuations when only NBI heating is used. Modeling of the core transport with TGYRO using the TGLF turbulent transport model and NEO neoclassical transport reproduces the experimental profile changes upon application of direct electron heating and indicates that multi-scale transport mechanisms are responsible for changes in the temperature and density profiles. Intermediate and high-k fluctuations appear responsible for the enhanced electron thermal flux, and intermediate-k electron modes produce an inward particle pinch that increases the inverse density scale length. Projection to ITER is performed with TGLF and indicates a density profile that has a finite scale length due to intermediate-k electron modes at low collisionality and increases the fusion gain. For a range of E×B shear, the dominant mechanism that increases fusion performance is suppression of outward low-k particle flux and increased density peaking.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The role of ion and electron-scale turbulence in setting heat and particle transport in the DIII-D ITER baseline scenario

In this work, the heat and particle transport in a DIII-D ITER Baseline Scenario (IBS) discharge has been investigated using both linear and nonlinear gyrokinetic simulations performed with the CGYRO code [J. Candy et al: Journal Comp. Phys. (2016)].These simulations were used to investigate the role that ion-scale ($k_θp_s$ < 1:0) and electron-scale ($k_θp_s$ > 1:0) turbulence play in determining heat and particle transport in the core of conditions believed to directly extrapolate to ITER operation. This investigation spans over nearly half of the plasma minor radius, from ρ = 0:45 - 0:85. To probe the nature of the transport and turbulence in these conditions and to validate the model against experimental fluxes, scans of $a/L_{T_i}$, $a/L_{T_e}$, and $a/L_n$ were performed at all radial locations. Long wavelength turbulence is found to be dominated by ITG modes with strong non-adiabatic electron effects and appears unable to reproduce ion and electron heat fluxes and electron particle fluxes simultaneously at most radial positions when single parameter scans are performed. To investigate the potential role of short wavelength turbulence, a series of electron-scale simulations are performed that indicate that experimentally relevant levels of electron heat flux could arise at sub-ion scales. Quantitative comparison of the simulated fluxes obtained from ion and electron-scale simulations with experimental levels along with the sensitivity of simulations results to changes in inputs is presented. Through extensive sensitivity scans and comparison with multiple transport channels, this work demonstrates a clear need for self-consistent modification of multiple inputs and suggests multi-scale interactions play a significant role at many of the radial locations studied. The results of this analysis help shape our understanding of the model fidelity needed to predict turbulence and transport reactor-relevant conditions and have important implications for the prediction of future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Impact of n > 1 neoclassical tearing modes on fast ions, plasma rotation, and the onset of disruptive tearing modes in DIII-D ITER baseline scenario

We present integrated TRANSP-kick analyses quantitatively characterizing how n > 1 magnetic islands redistribute fast ions (FIs) and modify momentum balance in low-torque DIII-D H-mode scenarios operated with the ITER normalized parameter set and shape (ITER baseline scenario, or IBS). In this plasma scenario, disruptive neoclassical tearing modes (NTMs) are most commonly seeded by nonlinear 3-wave coupling when the differential rotation between the q = 1 and q = 2 surfaces ($δf_{1, 2}$) approaches zero. Analysis of a DIII-D IBS database of shots unstable to the 2/1 NTM demonstrates that the flattening rotation profile is correlated with n > 1 NTM amplitude growth, with a strong drop in $δf_{1, 2}$ occurring when the radial magnetic perturbation amplitude at the rational surface ($δB$ r s ) is estimated to be about 20 G. The interpretive TRANSP-kick simulations show that the experimentally observed reduction in core rotation is due to the fast-ion kicks. Nonlinearity in momentum losses versus NTM amplitude is found to be correlated with overlap of FI island chains, occurring when $δB$ r s is about 35 G, in rough agreement with the experimental results. This agreement suggests the primary mechanism for the conductive momentum loss in the core is due to chaotic FI orbits and FI losses developing from overlapping resonant FI island chains.

ITER baseline scenario↗
DOE OSTI2024

Power and isotope effects in the ITER baseline scenario with tungsten and tungsten-equivalent radiators in DIII-D

Abstract Experiments in DIII-D document the ITER Baseline Scenario (IBS) at q 95 ∼ 3 and P IN / P LH ∼ 1–2, in both deuterium and hydrogen utilizing Kr and Xe as Tungsten-equivalent radiators. The power threshold for H-mode operation ( P LH ) was determined experimentally without added impurities and found to be about a factor of two higher than the scaling law. In recent IBS experiments in deuterium, intrinsic levels of metals such as Tungsten (W) or molybdenum and inconel are present that reduce the pedestal pressure by 20%–25%. A complete radiative collapse of deuterium IBS plasmas occurs at W core concentrations C W = 10 −5 . Simulations show that for core temperatures expected for ITER, the plasmas would not have a radiative collapse at C W = 1 × 10 −5 , moreover Q = 8–10 would still be achieved for C W up to 3 × 10 −5 . In contrast to deuterium, the IBS in hydrogen is not affected by intrinsic high-Z impurities, indicating that hydrogen H-modes in ITER may not inform the D-T phase with respect to W accumulation and discharge survival. Compared to deuterium, the pedestal pressure in hydrogen is ∼25% lower, with much higher ELM frequency of 150 Hz, decreasing with input power. Krypton was injected in a matrix scan of input power and impurity flow in IBS hydrogen discharges. Krypton impurity density profiles in hydrogen are similar to deuterium plasmas, but at Kr flows that are 2–3 times higher for the same input power. Krypton is transported into the core and affects the whole radius; at the highest injection rates a radiative collapse occurs at core radiation fractions of 0.3–0.35, consistent with the expected maximum W radiation fraction for ITER core plasmas. Comparing the results with previous International Tokamak Physics Activity database studies of the IBS confirms that at higher radiation fraction due to high-Z impurities, a drop in H 98 of >10% is observed. On the other hand, the results using Kr as a W-equivalent radiator indicate that metal (W) devices at lower core temperatures than ITER may provide overly pessimistic performance extrapolations to ITER for deuterium-tritium operation. The new DIII-D results support a more attractive option for the ITER Research Plan with a short hydrogen phase for system commissioning, transitioning to deuterium operations as soon as possible to provide relevant conditions for deuterium-tritium operations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modelling and experiment to stabilize disruptive tearing modes in the ITER baseline scenario in DIII-D

Abstract The achievement of high gain, stationary conditions in a tokamak scenario aimed at producing fusion energy in the ITER Project is crucial to the demonstration that this form of energy can be used in future reactors to provide cheap and clean energy globally. Disruptions are a challenge for the fusion energy field, in particular for the ‘ITER Baseline Scenario’ (IBS), as reproduced in the DIII-D tokamak. This work shows that a solution has been found for the m = 2/ n = 1 tearing modes that have consistently caused disruptions in the IBS: stable operation down to zero input torque was achieved by modifying the current density profile at the beginning of the pressure flattop and the ELM character later in the discharges, guided by previous results showing that the most likely cause of these instabilities is the current density profile. The coupling between sawteeth, n>2 modes and the 2/1 TMs is shown to not be statistically significant, nor the leading origin for the evolution towards instability. Ideal and resistive MHD modeling provide positive verification that a steeper ‘well’ in the region of the q = 2 rational surface leads to worse ideal stability, higher tearing index Δ’ and lower threshold Δ’ c for resistive instabilities, consistent with the experimental results. This provides confidence that the methods used in this work can be extrapolated to other devices and applied to avoid disruptions in ITER and pulsed fusion devices worldwide.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Radiation induced non-linear oscillations in ITER baseline scenario plasmas in DIII-D

Abstract This work shows how the radiation brought about by metals or metal-equivalent radiators such as Kr and Xe produces non-linear dynamics on otherwise stationary β N flattops of DIII-D ITER Baseline Scenario demonstration discharges. The Kr and Xe gases are used to reproduce the radiative loss rates of W in present machines that operate at core temperatures much lower than the expected ITER temperature. Experiments on DIII-D with injection of Kr and Xe, as well as with sources of intrinsic metals reach the range of radiated fraction values expected in the ITER core and experience slow oscillations in temperature and radiated power. In many cases of high radiated fraction, the core temperature decreases enough for the safety factor profile to rise above the 1/1 rational surface, naturally eliminating sawteeth and occasionally producing a persistent helical core. The oscillations can be reproduced by a modified Lotka–Volterra system for temperature and radiated fraction if diffusion and noise are included, which indicates that the interplay between temperature and radiation can be the main cause of the cyclic nature of the system. A new physics based model which includes equations for temperature, density and input power can also reproduce the oscillations observed in the experiments. The present results suggest that the non-linearity of the system can be increased by the inclusion of the inherently non-linear alpha heating term, which is proportional to ∼ n e 2 T i 2 , and obtains oscillations in the model when added to an otherwise more stationary system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Direct preemptive stabilization of $m,n$ = 2,1 neoclassical tearing modes by electron cyclotron current drive in the DIII-D low-torque ITER baseline scenario

We report the first preemptive direct stabilization of $m,n$ = 2,1 neoclassical tearing modes (NTMs) in DIII-D low-torque ITER baseline scenario plasmas by electron cyclotron waves. These experiments strongly support that the observed stabilization is achieved by replacing the missing bootstrap current in the island by electron cyclotron current drive (ECCD). These plasmas (with and without ECCD) are stable to 2,1 NTMs when the differential rotation between the q = 1 and q = 2 surfaces (Δf 1,2 ) is sustained above a critical value (Δf CRIT ≈ 1 kHz), but those evolving to low differential rotation (Δf 1,2 < Δf CRIT ) routinely develop a disruptive 2,1 island at the time and frequency of a sawtooth precursor. Preemptive, local and direct stabilization by ECCD is tested by scanning the current drive amplitude and location inside and outside of the q = 2 rational surface shot by shot. Analysis of only low differential rotation time windows, i.e. when stabilization from Δf 1,2 is absent, shows focused ECCD at q = 2 prevents the onset of 2,1 NTMs. These results are important, as they give the first demonstration of disruptive NTM control by ECCD in the plasma scenario planned for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Chaos as the cause of randomness in tearing mode onset times in DIII-D ITER baseline scenario plasmas

The first evidence of widespread chaotic plasma dynamics relevant to the onset of disruptive tearing modes is found in DIII-D tokamak discharges. This evidence is obtained by calculating positive maximal Lyapunov exponents from Mirnov signals using the Rosenstein algorithm, implying the chaotic growth and/or rotation of magnetohydrodynamic modes in DIII-D [Luxon, Nucl. Fusion 42, 614 (2002)] plasmas. Such chaotic dynamics offer an explanation for the random tearing mode onset times previously observed in ITER baseline scenario discharges. Further, the Lyapunov timescales associated with this chaos could inhibit reliable prediction of tearing mode onset beyond the angular momentum confinement timescale; this constraint should inform the design of plasma control systems on ITER.

Chaotic dynamics↗
DOE OSTI2024

Overview of the EUROfusion Tokamak Exploitation programme in support of ITER and DEMO

Within the 9th European Framework programme, since 2021 EUROfusion is operating five tokamaks under the auspices of a single Task Force called ‘Tokamak Exploitation’. The goal is to benefit from the complementary capabilities of each machine in a coordinated way and help in developing a scientific output scalable to future largre machines. The programme of this Task Force ensures that ASDEX Upgrade, MAST-U, TCV, WEST and JET (since 2022) work together to achieve the objectives of Missions 1 and 2 of the EUROfusion Roadmap: i) demonstrate plasma scenarios that increase the success margin of ITER and satisfy the requirements of DEMO and, ii) demonstrate an integrated approach that can handle the large power leaving ITER and DEMO plasmas. The Tokamak Exploitation task force has therefore organized experiments on these two missions with the goal to strengthen the physics and operational basis for the ITER baseline scenario and for exploiting the recent plasma exhaust enhancements in all four devices (PEX: Plasma EXhaust) for exploring the solution for handling heat and particle exhaust in ITER and develop the conceptual solutions for DEMO. The ITER Baseline scenario has been developed in a similar way in ASDEX Upgrade, TCV and JET. Key risks for ITER such as disruptions and run-aways have been also investigated in TCV, ASDEX Upgrade and JET. Experiments have explored successfully different divertor configurations (standard, super-X, snowflakes) in MAST-U and TCV and studied tungsten melting in WEST and ASDEX Upgrade. The input from the smaller devices to JET has also been proven successful to set-up novel control schemes on disruption avoidance and detachment.

ASDEX upgrade↗
DOE OSTI2023

Analysis of beam ion driven Alfvén eigenmode stability induced by Tungsten contamination in EAST

Alfvén eigenmodes (AE) activity is observed in the EAST high β N and low B T discharge 93910, operation scenario dedicated to explore the ITER baseline scenario. AEs are triggered after the plasma is contaminated by Tungsten that causes an abrupt variation of the thermal plasma and energetic particles (EPs) profiles. The aim of the present study is to analyze the AE stability in the 93910 discharge using the gyro-fluid code FAR3d, identifying the AE stability trends by comparing the plasma before and after the Tungsten contamination. Tungsten contamination causes the destabilization of Toroidal AEs (TAE) and Energetic particle modes (EPMs) in the same frequency range and radial location with respect to the experimental observation and M3D-K/GTAW code results. Next, a set of parametric studies are performed to analyze the effect of the thermal plasma and EP parameters on the AE stability. The analysis indicates a lower EP β threshold for the AEs destabilization if the EP energy increases, an improved AE stability of on-axis NBI configurations due to the stronger continuum damping in the inner plasma region as well as a large enhancement of the EP drive as the thermal ion density increases due to a higher ratio of the EP and Alfven velocities. Consequently, the simulations indicate the increment of the thermal ion density after the Tungsten contamination could be the main cause of the AE/EPM destabilization.

AE↗
DOE OSTI2024

First tungsten radiation studies in DIII-D’s ITER baseline demonstration discharges

Abstract ITER Baseline Scenario plasmas were studied in DIII-D using krypton and xenon gases as a proxy for the tungsten that will be present in ITER. These impurities were chosen for having the same radiative loss rate Lz as tungsten would exhibit in the hotter ITER core. Results show that the scenario with these core radiators spans the range of impurity concentration and W radiated fraction expected for ITER, and up to 50% higher values, explored at zero injected torque, as well as 1 Nm and full co-torque injection with T ∼ 3 Nm. Stationary discharges with duration >2–4 τ R are achieved with f rad ⩾ 30% leading to a reduction in confinement of ∼10%, and a comparison with real metal radiators in the same range of f rad shows that the higher Lz at the lower temperatures in these plasmas yields too pessimistic results on the survivability and performance of this scenario in ITER. Simulations of ITER power balance including W radiation show that with concentration up to three times higher than in the DIII-D plasmas the scenario can be stationary, remaining at acceptable core radiated fraction values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Novel intrinsic helical cores and MHD dynamo flux pumping evidence in DIII-D

Evidence of a MHD dynamo brought about by intrinsic helical core structures in negative triangularity (NT) and ITER baseline scenario (IBS) plasmas demonstrates paths to operating tokamaks without seeds for deleterious global MHD modes. A stationary sawtooth- and edge localized mode-free state maintains high performance for long durations in a strong NT shape in the DIII-D tokamak. Disparity between minimum safety factor values predicted by magnetic flux diffusion theory and derived from experiment suggests anomalous flux diffusion by an MHD dynamo, driven by flows associated with quasi-interchange modes. Stability calculations with the GATO code find NT and IBS experimental equilibria unstable to quasi-interchange modes, in agreement with observed saturated activity. IBS discharges with similar magnetic winding structure occupy robust helical core states following a transition from sawtooth instabilities. These states agree with prior modeling of helical core bifurcation thresholds with the VMEC equilibrium code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Error field predictability and consequences for ITER

Abstract ITER coil tolerances are re-evaluated using the modern understanding of coupling to least-stable plasma modes and an updated center-line-traced model of ITER’s coil windings. This reassessment finds the tolerances to be conservative through a statistical, linear study of n = 1 error fields (EFs) due to tilted, shifted misplacements and nominal windings of central solenoid and poloidal field coils within tolerance. We also show that a model-based correction scheme remains effective even when metrology quality is sub-optimal, and compare this to projected empirical correction schemes. We begin with an analysis of the necessity of error field correction (EFC) for daily operation in ITER using scalign laws for the EF penetration threshold. We then consider the predictability of EF dominant mode overlap across early planned ITER scenarios and, as measuring EFs in high power scenarios can pose risks to the device, the potential for extrapolation to the ITER Baseline Scenario (IBS). We find that carefully designing a scenario matching currents proportionally to those of the IBS is far more important than plasma shape or profiles in accurately measuring an optimal correction current set.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Saturation of Fishbone Instability by Self-Generated Zonal Flows in Tokamak Plasmas

Gyrokinetic simulations of the fishbone instability in DIII-D tokamak plasmas find that self-generated zonal flows can dominate the nonlinear saturation by preventing coherent structures from persisting or drifting in the energetic particle phase space when the mode frequency down-chirps. Results from the simulation with zonal flows agree quantitatively, for the first time, with experimental measurements of the fishbone saturation amplitude and energetic particle transport. Moreover, the fishbone-induced zonal flows are likely responsible for the formation of an internal transport barrier that was observed after fishbone bursts in this DIII-D experiment. Finally, gyrokinetic simulations of a related ITER baseline scenario show that the fishbone induces insignificant energetic particle redistribution and may enable high performance scenarios in ITER burning plasma experiments. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Disruptive neoclassical tearing mode seeding in DIII-D with implications for ITER

New studies identify the critical parameters and physics governing disruptive neoclassical tearing mode (NTM) onset. An m/n = 2/1 mode in DIII-D that begins to grow robustly after a seeding event (edge localized mode ELM or sawtooth precursor and crash) causes the mode rotation to drop close to the plasma's Er = 0 rest frame; this condition opens the stabilizing ion-polarization current 'gate' and destabilizes an otherwise marginally stable NTM. Our new experimental and theoretical insights and novel toroidal theory-based modeling are benchmarked and scalable to ITER and other future experiments. Here, the nominal ITER rotation at q = 2 is found to be stabilizing ('gate closed') except for MHD-induced transients that could 'open the gate'. Extrapolating from the DIII-D ITER baseline scenario (IBS) discharges, MHD transients are much more likely to destabilize problematic robustly growing 2/1 NTMs in ITER; this makes predictions of seeding and control of both ELMs and sawteeth imperative for more than just minimizing divertor pulsed-heat loading.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Prediction of performance and turbulence in ITER burning plasmas via nonlinear gyrokinetic profile prediction

Burning plasma performance, transport, and the effect of hydrogen isotope (H, D, D-T fuel mix) on confinement has been predicted for ITER baseline scenario (IBS) conditions using nonlinear gyrokinetic profile predictions. Accelerated by surrogate modeling (Rodriguez-Fernandez et al 2022 Nucl. Fusion 62 076036), high fidelity, nonlinear gyrokinetic simulations performed with the CGYRO code (Candy et al 2016 J. Comput. Phys. 324 73), were used to predict profiles of T i , T e , and n e while including the effects of alpha heating, auxiliary power (NBI + ECH), collisional energy exchange, and radiation losses inside of $r/a$ = 0.9. Predicted profiles and resulting energy confinement are found to produce fusion power and gain that are approximately consistent with mission goals ($P_\textrm{fusion} = 500$ MW at Q = 10) for the baseline scenario and exhibit energy confinement that is within 1σ of the H-mode energy confinement scaling. The power of the surrogate modeling technique is demonstrated through the prediction of alternative ITER scenarios with reduced computational cost. These scenarios include conditions with maximized fusion gain and an investigation of potential resonant magnetic perturbation (RMP) effects on performance with a minimal number of gyrokinetic profile iterations required (3–6). These predictions highlight the stiff ITG nature of the core turbulence predicted in the ITER baseline and demonstrate that $Q \gt$ 17 conditions may be accessible by reducing auxiliary input power while operating in IBS conditions. Prediction of full kinetic profiles allowed for the projection of hydrogen isotope effects around ITER baseline conditions. The gyrokinetic fuel ion species was varied from H, D, and 50/50 D-T and kinetic profiles were predicted. Results indicate that a weak or negligible isotope effect will be observed to arise from core turbulence in IBS conditions. The resulting energy confinement, turbulence, and density peaking, and the implications for ITER operations will be discussed.

gyrokinetics↗