Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA DISCHARGE”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Generation of shear flows induced by AE / EPM in LHD plasma

The generation of shear flows (SFs) by Alfven Eigenmodes (AEs) and energetic particle modes (EPMs) have important effects on the operation of future nuclear fusion reactors, because SFs regulate the saturation of the AEs/EPMs, the transport of EPs and thermal plasma, as well as the formation of transport barriers among other consequences. The aim of this study is the analysis of SFs generation during the saturation phase of AEs and EPMs in LHD plasma. Experiments performed in the 23rd and 24th LHD experimental campaigns are dedicated to explore the destabilization of AEs/EPMs in discharges with different heating patterns, thermal plasma and magnetic field configurations. In particular, the shots 176490 and 179697 show the destabilization of MHD bursts and energetic-ion-driven resistive interchange modes (EIC), respectively. Charge exchange spectroscopy measurements in both discharges indicate that the generation of SFs by AE/EPM is uncorrelated with the perturbation induced by the neutral beam injector (NBI). Nonlinear simulations performed using the gyro-fluid code FAR3d show the generation of zonal structures, especially SFs, induced during the saturation phase of Toroidal Alfven Eigenmodes (TAEs) triggered in the MHD burst as well as by the 1/1 EIC in the bursting phase. The simulations indicate that SFs are caused by the radial electric fields powered by energy transfers from the unstable AE/EPM towards the thermal plasma. The strongest SFs are measured during the EIC bursting phase once the 1/1 EPM overlaps with nearby resonances at the plasma periphery. Likewise, the largest SFs during the MHD burst are observed once TAEs radially overlap in the inner-middle plasma region.

AE↗
DOE OSTI2026

A Fast and Accurate Reduced-Order Model for High-Intensity Transferred Arc Discharges

Arc discharges are widely used in welding, plasma smelting, and other industrial processes, where variations in operating conditions strongly affect arc stability, temperature distribution, and energy transfer. Accurate modeling of these phenomena typically requires computationally expensive high-fidelity simulations. This study presents a hierarchy of three arc discharge models with progressively reduced physical fidelity. The high-fidelity model provides a fully physics-resolved reference, the reduced-order model reproduces the dominant thermal and flow characteristics of the arc with reasonable accuracy, and the Elenbaas-Heller model captures key trends in a simplified, rapid formulation. This hierarchy demonstrates that reduced-order model can effectively balance predictive fidelity and computational efficiency, providing practical tools for arc simulation and parametric studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Improved liquid lithium surfaces in the Lithium Tokamak Experiment-β

Advances in vacuum, surface, and lithium conditioning techniques throughout five years of continuous operations in LTX-β have produced mirror-like liquid lithium surfaces and demonstrated the feasibility of high-performance tokamak discharges fully surrounded by liquid metal without significant operational problems. Improvements in conditioning techniques and procedures, including many weeks of baking and accumulation of 70 g of Li, led to reduced residual gasses and clean Li surfaces - all while still maintaining enough operational flexibility for multiple in-vacuum diagnostic upgrades and calibrations. Coatings had a visibly clean appearance, with reflective liquid metal demonstrating good wetting and surface adhesion with films that were now macroscopically thick. Solidified Li showed large crystal grains, while surface science measurements observed reduced impurities in the lithium. Steadily improved plasma performance was achieved with liquid lithium, with discharges able to match solid Li in terms of evolution of I p and n e , including rapid density pumping indicating low recycling. There were indications of moderately increased Li impurity influx, though few significant disturbances by the large liquid surfaces on tokamak operations over hundreds of discharges. Liquid metal plasma facing components are a potential solution to the extreme heat and particle fluxes that could cause unacceptable damage to solid materials, while liquid lithium also has the potential for greatly increased confinement in the low-recycling regime. While many liquid metal approaches are possible, and numerous experiments have been conducted in test stands and small modules in fusion devices, LTX-β is the only tokamak operated while fully surrounded by liquid metal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Unified ELM Suppression on KSTAR and DIII-D via Adaptive Feedback Control Strategies

This paper reports on the extension of our amplitude-based resonant magnetic perturbation (RMP) edge localized mode (ELM) controller to support phasing control (relative toroidal phases of RMP waveforms between rows of coils), multiple toroidal mode numbers, and new ‘jump’ and ‘probing’ strategies, all deployed on KSTAR and DIII-D. By treating the control algorithm as device-independent and adjusting only the real-time interfaces to sensors and power supplies, we have confirmed that the same finite state machine—based feedback logic can be ported between machines with minor modifications. In experiments using n = 2 RMPs on KSTAR and n = 3 on DIII-D, the controller successfully modulated RMP amplitudes in real time to sustain ELM suppression while minimizing confinement degradation. Phasing control broadened the suppression window, as it permitted the system to avoid locked-mode regions and safely access ELM-free conditions. A rotating RMP phasing scheme, integrated into the same framework, distributes divertor heat loads more uniformly, making it a promising strategy for protecting plasma-facing components during long discharges. New ‘jump’ and ‘probing’ techniques demonstrate the possibility for the controller to preempt imminent ELMs and refine the minimum required RMP amplitude without returning to ELMy conditions. Taken together, these upgrades enable extended ELM-free operation while mitigating confinement degradation, providing a practical framework for real-time ELM control in future high-performance tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Oxygenates production in a microfluidic dielectric barrier discharge device sustained in Ar/CH 4 /O 2

Reforming of methane (CH 4 ) is a process to produce syngas (CO/H 2 ) and other value-added chemicals including oxygenates such as methanol (CH 3 OH). Atmospheric pressure plasmas have the potential to be more energy efficient than traditional reforming methods as value-added chemicals can be synthesized directly in the plasma without requiring an additional step. In this paper, we discuss the results from a computational investigation of the formation of oxygenates by CH 4 oxidation in the presence of Ar, including CH 3 OH and CH 2 O, in a nanosecond pulsed dielectric barrier discharge. The plasma is formed in a microfluidic channel whose small dimensions are ideal for plasma formation at atmospheric pressure. The production and consumption mechanisms of dominant radicals and long-lived species are discussed in detail for the base case conditions of Ar/CH 4 /O 2 = 50/25/25. CH 3 OH is produced primarily by CH 3 O reacting with CH 3 O and CH 3 O 2 reacting with OH, while CH 2 O formation relies on reactions involving CH 3 O and CH 3 . The most abundant oxygenate formed is CO (produced by H abstraction from CHO). However, the greenhouse gas CO 2 is also formed as a by-product. The effects of gas mixture are examined to maximize the CH 3 OH and CH 2 O densities while decreasing the CO 2 density. Increasing the Ar percentage from 0% to 95% decreased the CH 3 OH and CH 2 O densities. At low Ar percentages, this is due to an increase in consumption of CH 3 OH and CH 2 O, while at high Ar percentages (>40% Ar), the production of CH 3 OH and CH 2 O is decreased. However, both CO and CO 2 reached peak densities at 70%–90% Ar. Changing the CH 4 /O 2 ratio while keeping 50% Ar in the discharge led to increased CH 3 OH and CH 2 O production, reaching peak densities at 35%–40% CH 4 . The CO and CO 2 densities decreased beyond 20% CH 4 , indicating that a CH 4 rich discharge is ideal for forming the desired oxygenates.

03 NATURAL GAS↗
DOE OSTI2025

Global Gyrokinetic Simulations of Isotope Effects under Ambipolar Electric Fields and Advances Toward Whole-Volume Modeling

We review global gyrokinetic simulation studies on plasma transport in the Large Helical Device using XGC-S. XGC-S is an extended version of X-point Gyrokinetic Code for stellarators and has been progressively verified throughout the code development process. Verification tests of neoclassical transport successfully demonstrate the generation of an ambipolar electric field due to ripple-trapped particles. We perform quasi-linear analyses of the ion temperature gradient mode under the influence of the ambipolar electric field. The results reveal that the ambipolar electric field and the heavy hydrogen component in mixed isotope plasmas can lead to the favorable isotope effect observed in recent deuterium experiments. We also present recent efforts in code development toward whole-volume simulations, including the helical divertor region. A mesh generation scheme based on field-line tracing and the construction of curved surfaces perpendicular to the magnetic field would be promising for global field calculations in the whole-volume simulations.

Basic Plasma Phenomena and Gas Discharges↗
DOE OSTI2024

Energy amplification in plasma crystals due to multiple torsions

Interacting torsions are examined within a two-dimensional monolayer crystal suspended in an argon complex plasma for 1–10 W discharge powers and pressures of 135–155 mTorr. Two torsions embedded in a lattice are shown to amplify the kinetic energy and range of motion of particles located between the torsions to nearly three times that observed in single torsion systems. It is also shown that multiple torsions can interact via amplified particle energy when separated by up to 14 interparticle distances (Δ). The torsion separation distance also showed a positive linear trend with power and a slightly positive correlation with the pressure. This amplification of energy is possible due to the fact that multiple torsions in a lattice increase the interparticle distance of the lattice by 16% more than single torsion systems, leading to additional freedom of motion in the lattice plane. These combined findings show that multiple torsions heat the lattice differently depending on their separation from the other torsion. The midpoint particles between torsions absorb the majority of energy from the two torsions, and energy addition at the midpoint is nonlinear. The addition of more torsions to the lattice may lead to melting of the plasma crystal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Spatial profile of argon (1s 5 ) metastables in an electron beam generated plasma

Electron beams with an applied magnetic field generate a secondary cold plasma with a selective chemical composition, featuring low-energy ions and metastable species in the discharge periphery, ideal for low-damage plasma treatment of material substrates. In this work, we studied the plasma generated by an e-beam using a 4 kV voltage in a pure argon gas environment under a magnetic field of 150 G and in the pressure range of 25–90 mTorr. We measured the absolute spatial density profile of argon (1 s 5 ) metastables in an electron beam generated plasma by laser-induced fluorescence and found it to be of the order of 10 16 m −3 . The electron temperature and the electron density measured by a Langmuir probe were of the order of 10 16 m −3 and less than an eV respectively. Electron-impact quenching was identified as a significant loss mechanism for the Ar(1s 5 ) state, leading to the saturation of the metastable density at higher pressures. Outside the primary ionization region, the spatial distribution of argon metastables followed a linear diffusion profile, indicating negligible additional production in those regions.

EEDF↗
DOE OSTI2026

Validation of the GFS model for gyrokinetic stability of NSTX pedestal data

This study presents a large database validation of the gyro fluid system (GFS) model for linear gyrokinetic stability for high-mode (H-mode) edge transport barrier conditions in the national spherical torus experiment (NSTX) tokamak. The database of linear stability calculations with the CGYRO gyrokinetic code was produced using plasma profile measurements from NSTX discharges to identify kinetic ballooning modes (KBM), trapped electron modes (TEM), and micro-tearing modes (MTM) that limit the pressure profile gradient in the H-mode barrier. A novel Bayesian optimization approach determines optimal resolution parameters for GFS specifically for spherical tokamak pedestal conditions. Our results demonstrate that GFS, with optimized resolution, can achieve accurate linear stability analysis in NSTX pedestal conditions for reduced resolution compared to CGYRO. GFS can accurately find the KBM, TEM, and MTM instability branches. Parametric analysis reveals that GFS accuracy in this extreme pedestal parameter range is degraded for low magnetic shear and near the separatrix conditions. These findings establish GFS as a fast linear eigenmode solver for spherical tokamak pedestal gyrokinetic stability and demonstrate a systematic methodology for determining the optimum resolution settings.

Yang, Minglei [Oak Ridge National Laboratory (ORNL↗
DOE OSTI2026

Multi‐Scale Model‐Informed Deep Learning for Plasma‐Nanoparticle Interaction

The Overarching Goal of this proposed research is to understand and quantitively determine the interactions between non-thermal plasma (hot electrons, reactive radicals, vibrationally excited species) and surface reactions on influencing the activity and selectivity of the desired reactions via developing multi-scale model informed deep learning algorithm. Investigating non-thermal plasma-surface interaction is feasible due to the low bulk temperature in the discharge region. To investigate the role of plasma-nanoparticle interaction on enhancing the reaction kinetics, we will focus on ammonia cracking to generate clean hydrogen over earth-abundant, non-critical metallic nanoparticles, which is of great significance for decarbonization. We hypothesize that (1) reactive radicals interacting with surface reaction species via Eley–Rideal mechanism will significantly lower the energetics of the potential rate-limiting step of nitrogen formation; (2) the surface will be charged heterogeneously under non-thermal plasma conditions and the charged site will lower the energetics of ammonia cracking through Langmuir– Hinshelwood mechanism; (3) vibrationally excited ammonia will further promote the initial N-H bond cleavage. To access the hypothesis, we will (1) reveal the surface charge effects on tunning the reaction energetics via interpretable, physics-informed deep learning accelerated density functional theory (DFT) calculations; (2) determine the reactive radicals interacting with surface reaction species on tuning the reaction energetics via DFT; (3) reveal the surface charge effects on tunning the reaction energetics via DFT and deep learning models, (4) quantify how vibrationally excited species, reactive radicals, and surface charging effects on enhancing the catalysis via developing DFT-based microkinetic modeling (MKM) and active learning. Deep and active learning of plasma-nanoparticle interactions effects on enhancing ammonia cracking to generate hydrogen represents a new paradigm for designing high performance plasma materials. The fundamental science of how plasma-nanoparticle interactions will change the plasma kinetics and will improve the energy efficiency for decarbonization and sustainability. The interpretable and physics-informed machine learning model will accelerate low temperature plasma chemistry and material discovery with physics rules and model interpretation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Initial design concepts for solid boron injection in ITER

As part of ITER’s consideration to change its first wall material from beryllium to tungsten, the ITER Organization has proposed studying the feasibility of real-time solid boron injection (SBI) into the plasma to coat the walls and divertor to supplement glow discharge boronization (GDB). Boron deposits getter oxygen and reduce sputtering of tungsten from plasma facing components (PFCs). Particularly in areas with significant plasma wall interactions, boron coatings are expected to be short-lived under high performance plasma conditions. The proposed SBI system aims to maintain boron layers in these areas to avoid excessive radiation from tungsten in the plasma as a risk mitigation to ensure ITER will be able to reach and sustain Q = 10 conditions. The system will be used sparingly, as redeposition of boron can lead to significant tritium retention, which must be minimized in ITER to comply with nuclear safety concerns. SBI is proposed to limit and precisely control the amount of boron injected in real-time during plasma operation. Here, some of the design requirements and initial concepts for an SBI system in ITER are presented based on previous results carried out with SBI systems on a number of tokamaks and stellarators around the world. Previous results using SBI systems installed by PPPL have injected boron particles from 5 µm–2 mm diameter at calibrated rates of 2–200 mg/s in real-time during plasma operation on AUG, DIII-D, EAST, KSTAR, LHD, TFTR, WEST, and W7-X, leading to improved wall conditions with reduced plasma impurity concentrations and radiated power and improved plasma performance. The boron is ionized in the plasma edge and then deposited on plasma-wetted surfaces. On AUG, EAST and WEST reduced tungsten sputtering sources were observed following several discharges with SBI. Extrapolation of these SBI results are presented to estimate the amount of boron needed for wall conditioning in ITER. Real-time SBI control requirements and plasma operation scenarios for ITER are also described.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Azimuthal and axial structures in 3D particle-in-cell simulation of Penning discharge

We report the results of a 3D particle-in-cell simulation of cylindrical Penning discharge with the axial magnetic field in the so-called reflex configuration, where the cathode and anti-cathode are biased to the same negative potential. The discharge is supported by thermal electron emission from the cathode. Electron and ion collisions, including ionization, are fully accounted for. The emphasis is on a specific regime in which the plasma potential at the center of the discharge is positive with respect to the chamber walls, serving as an anode. Spatial and temporal scales of the observed azimuthal and axial fluctuations and structures are characterized. It is suggested that azimuthal structures are caused by the dissipative gradient-drift instability. We find that the axial fluctuations related to the plasma-beam instabilities are weakly correlated with the azimuthal perturbations of the density, so that the azimuthal modes rotate as a whole and do not show any axial shear. The behavior of the electric potential is more involved, demonstrating intermittent standing wave and propagating structures in the axial direction that modulate the electron transport, producing the standing wave pattern (along the z-direction) in the radial electron flux.

3D PIC simulation of discharge↗
DOE OSTI2025

Combining physics-based and data-driven models for quantitatively accurate plasma profile prediction that extrapolates well; with application to DIII-D, AUG, and ITER tokamaks

For design, scenario planning, and control, ITER and all other envisioned tokamaks rely on a variety of statistical and physics-based models to extrapolate to unseen regimes; most notably from low plasma current to high. A 'meta-learning' methodology for combining the accuracy of data-driven models with the generalizability of physics-based models is described and tested, yielding a 5–10 percent improvement in performance beyond either alone for the task of extrapolating time-dependent plasma profile prediction from low- to high- plasma current DIII-D tokamak discharges. Meanwhile, it is shown that both machine learning models extrapolated far-distribution and state-of-the-art 'physics-based' profile predictors fare worse than merely assuming plasma profiles do not change from their initial values. Finally, a variety of other mechanisms for helping data-driven models generalize—transfer learning, adding contextual information from physics simulators, and adding data from the ASDEX Upgrade tokamak—are attempted for similar extrapolation tasks but, in the methodology used in this paper, yield no significant improvement beyond simple data-driven models. Results are summarized in figures 15 and 16.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A study of resistive peeling–ballooning modes across spherical tokamaks

We investigate how non-ideal-magnetohydrodynamics (MHD) effects, in particular plasma resistivity, impact the peeling–ballooning stability thresholds in spherical tokamaks. This analysis follows the discovery of resistive kink-peeling modes in ELMing National Spherical Torus Experiment (NSTX) discharges. In the present study we extend this modeling to ELMing pulses in the Mega Ampere Spherical Tokamak (MAST) and MAST—Upgrade (MAST-U), where we find a clear resistive scaling for peeling–ballooning modes. While in NSTX ideal-MHD predicts stability for ELMing discharges, in MAST-U we find that the plasma is slightly unstable to peeling–ballooning modes, but is fully stabilized once diamagnetic effects are considered in terms of a growth rate normalization. A resistive power law scaling is calculated for these modes on MAST-U, which lies in between that of tearing modes and resistive interchange modes. A comparison between M3D-C1 and NIMROD shows reasonable agreement for this scaling. Resistivity destabilizes the modes and the peeling–ballooning unstable domain is considerably expanded in both, MAST and MAST-U. In addition to the MAST/-U pulses we also analyze resistive PB stability in a NSTX-similarity discharge on DIII-D. While having a different aspect ratio from NSTX, this discharge uses NSTX-like shaping parameters, toroidal field and plasma current. By considering these discharges alongside NSTX cases, we identify conditions influencing the onset of resistive peeling–ballooning modes. Furthermore, our findings indicate that magnetic shear in the pedestal region is closely linked to the emergence of resistive edge modes.

MAST-U↗
DOE OSTI2025

Achievement of highly radiating plasma in negative triangularity and effect of reactor-relevant seeded impurities on confinement and transport

The first achievement of highly radiating plasmas in negative triangularity (NT) is shown with an operational space featuring high core radiation at high Greenwald fraction obtained with the injection of reactor-relevant seeded gases. These NT shape diverted discharges reach high values of normalized plasma pressure ( β N > 2) at high radiation fraction with no edge localized modes (ELMs) and a simultaneous reduction of the divertor heat flux. Here, we demonstrate that as long as the impurity level in the core is kept low to avoid excessive fuel dilution and impurity accumulation, integration of NT configuration with high radiation fraction not only is achievable but it can lead to confinement improvement with stabilization effects originating from collisionality, E × B shear and profiles changes due to impurity radiation cooling. The underlying physics mechanism is robust and reproducible and holds for a variety of impurity species. The absence of the requirement to stay in H-mode translates in a higher core radiation fraction potentially allowed in NT shape effectively mitigating the power exhaust issue. The results presented here demonstrate a path to high performance, ELM free and highly radiative regime with reactor-relevant seeding gases making this regime a potential new scenario for reactor design.

confinement improvement↗
DOE OSTI2025

Non-Equilibrium in a Dust-Forming Low-Temperature Plasma: A CARS Study

Dust-forming low-temperature plasmas are versatile systems for the production of nanoparticles with tunable functionalities. While attractive from a materials processing point of view, these systems are inherently complex, with several plasma-induced phenomena determining the properties of the produced materials. Here, we characterize a carbon nanoparticle-forming plasma using coherent anti-Stokes Raman spectroscopy (CARS), with the primary goal of measuring gas temperature. While gas temperature is typically assumed to be at or slightly above room temperature in these reactors, we measure gas temperatures exceeding 1000 K under typical process conditions. We find a correlation between the gas temperature and the nanoparticle yield, suggesting that the particle nucleation and growth process releases energy within the reaction volume, leading to significant gas heating. In addition, we find that the relaxation of vibrationally excited species at the particle surfaces is a major contributor to their heating. In conclusion, these results underscore the complexity of these systems and the need for their more in-depth characterization using advanced techniques such as CARS.

Basic Plasma Phenomena and Gas Discharges↗
DOE OSTI2026

Validation of an integrated modeling framework for investigating 3D plasma responses in tokamak plasmas

As contemporary experimental tokamaks are pushed toward reactor-relevant operation, they provide essential testbeds for demonstrating and exploring ELM control strategies for deployment in future fusion pilot plants. Accurate predictions of full plasma responses are essential to guide and optimize these demonstrations. This paper introduces and validates an integrated modeling framework over a historical range of DIII-D operational space. The integrated modeling framework uses only scalar plasma parameters and optional reference boundary to self-consistently and flexibly scan through tokamak operational space and estimate the corresponding ELM-suppression relevant plasma response over 3D coil phase space. The framework generates tightly converged equilibria that satisfy a target set of plasma parameters (I p , β N , l i ), with kinetic profiles constrained by an EPED(NN)-computed pedestal and empirical core model. The plasma response of these modeled equilibria is calculated with GPEC. Validated against 55 distinct DIII-D equilibria, the framework consistently reproduces experimentally constrained equilibrium pressure, q, and other representative profiles using the time-varying information from only evolving scalar plasma parameters; plasma response validation is performed for the n=3 perturbation. In particular, the plasma response of a DIII-D discharge scanning q 95 in search of ELM suppression windows is quantitatively reproduced. By accurately reproducing experimental equilibria and plasma responses across wide parameter variation, the framework supports potential for synthetic parameter scans in key operational, stability, and plasma response dimensions to investigate RMP ELM-suppression experiments and inform predictive RMP scenario optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Sensitivities of time-dependent temperature profile predictions for NSTX with the multi-mode model

The Multi-Mode Model (MMM) for turbulent transport was applied to a large set of well-analyzed discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate its sensitivities to a wide range of plasma conditions. MMM calculations were performed for hundreds of milliseconds in each discharge by performing time-dependent predictive simulations with the 1.5D tokamak integrated modeling code TRANSP. A closely related study (Lestz et al 2025 Plasma Phys. Control. Fusion 67 105029) concluded that MMM predicted electron and ion temperature profiles that were in reasonable agreement with NSTX observations, generally outperforming a different reduced transport model, TGLF. This finding motivates the more thorough investigation of the characteristics of the MMM predictions conducted in this work. The simulations with MMM have electron energy transport dominated by electron temperature gradient modes in the examined discharges with relatively low plasma β (ratio of kinetic plasma pressure to magnetic field pressure) and high collisionality, transitioning to a mixture of different modes for higher β and lower collisionality. The thermal ion diffusivity predicted by MMM is much smaller than the neoclassical contribution, in line with previous experimental analysis of NSTX. Nonetheless, the electron and ion temperature profiles are coupled via collisional energy exchange and thus sensitive to which transport channels are predicted. The time-dependent simulations with MMM are robust to the simulation start time, converging to remarkably similar temperature profiles later during the discharge. MMM typically overpredicts confinement relative to NSTX observations, leading to the prediction of overly steep temperature profiles. Plasmas with spatially broader temperature profiles, higher plasma β, and longer energy confinement times tend to be predicted by MMM with better agreement with the experiment. As a result, these findings provide useful context for understanding the regime-dependent tendencies of MMM in anticipation of self-consistent, time-dependent predictive simulations of NSTX-U discharges with these same modeling tools.

MMM↗