Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “negative triangularity”

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

MHD stability analysis against pressure and current-driven modes in the SMall Aspect Ratio Tokamak

Linear magnetohydrodynamic (MHD) simulations for the SMall Aspect Ratio Tokamak (SMART) have been carried out for the first time, for both positive (PT) and negative triangularity (NT) shaped plasmas using the MARS-F code. The MHD stability of projected SMART plasmas against internal kinks, infernal modes and edge peeling-ballooning modes have been analyzed for a wide range of realistic equilibria. A stabilization of internal kinks and infernal modes is observed when increasing the safety factor profile and reducing plasma beta. PT shaped plasmas are more stable against both internal kinks and infernal modes than their counterpart NT shaped plasmas. Toroidal flows have little impact on the MHD stability of the internal kinks, but they have a strong stabilizing effect on infernal modes, which can be further mitigated in NT shaped plasmas. The MHD stability of peeling-ballooning modes is reduced in NT shaped plasmas, as observed in conventional tokamaks.

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 OSTI2022

Effect of anisotropic fast ions on internal kink stability in DIII-D negative and positive triangularity plasmas

Recent DIII-D experiments show that sawtooth stability is strongly affected by anisotropic fast ions from neutral beam injection (NBI) in both negative and positive triangularity plasmas. Fast ions from co-current NBI are stabilizing for the sawtooth stability, resulting in longer sawtooth periods. On the other hand, fast ions from counter-current NBI are destabilizing, leading to small and frequent sawteeth. The relative change of sawtooth period and amplitude is more than a factor of two. These observations appear to hold in both plasma shapes. Non-perturbative toroidal modeling, utilizing the magnetohydrodynamic-kinetic hybrid stability code MARS-K (Liu et al 2008 Phys. Plasmas 15 112503), reveals an asymmetric dependence of the stability of the n = 1 (n is the toroidal mode number) internal kink mode on the injection direction of NBI, being qualitatively consistent with the experimentally observed sawtooth behavior. The MARS-K modeling results suggest that anisotropic fast ions affect the mode growth rate and frequency through both adiabatic and non-adiabatic contributions. Here, the asymmetry of the internal kink mode instability relative to the NBI direction is mainly due to the non-adiabatic contribution of passing fast ions, which stabilize (destabilize) the internal kink with the co-(counter-) current NBI as compared to the fluid counterpart. However, finite orbit width (FOW) correction to passing particles partially cancels the asymmetry. Trapped particles are always stabilizing due to precessional drift resonance. Modeling also shows that fast ions affect the internal kink in a similar manner in both negative and positive triangularity plasmas, although being slightly more unstable in the negative triangularity configuration already in the fluid limit. The similarity is mainly attributed to the fact that the mode is localized in the plasma core region, with very similar eigenmode structures in both negative and positive configurations. Furthermore, MARS-K modeling indicates that other factors, such as the plasma rotation and the drift kinetic effects of thermal plasmas, weakly modify the mode stability as compared to the drift kinetic resonance effects and FOW correction of fast ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Towards fully predictive gyrokinetic full- f simulations: validation and triangularity studies in TCV

Designing economical magnetic confinement fusion power plants motivates computational tools that can estimate plasma behavior from engineering parameters without direct reliance on experimental measurement of the plasma profiles. In this work, we present full-f global long-wavelength gyrokinetic simulations of edge and scrape-off layer turbulence in tokamaks that use only magnetic geometry, heating power, and particle inventory as inputs. Unlike many modeling approaches that employ free parameters fitted to experimental data, raising uncertainties when extrapolating to reactor scales. This approach directly simulates turbulence and resulting profiles through gyrokinetics without such empirical adjustments. This is achieved via an adaptive sourcing algorithm in Gkeyll that strictly controls energy injection and emulates particle sourcing due to neutral recycling. We show that the simulated kinetic profiles compare reasonably well with Thomson scattering and Langmuir probe data for Tokamak á Configuration Variable (TCV) discharge #65125, and that the simulations reproduce characteristic features such as blob transport and self-organized electric fields. Applying the same framework to study triangularity effects suggests mechanisms contributing to the improved confinement reported for negative triangularity (NT). Simulations of TCV discharges #65125 and #65130 indicate that NT increases the E x B flow shear (by about 20% in these cases), which correlates with reduced turbulent losses and a modest change in the distribution of power exhaust to the vessel wall. While the physical models contain approximations that can be refined in future work, the predictive capability demonstrated here, evolving multiple profile relaxation times with kinetic electron and ion models in hundreds of GPU hours, indicates the feasibility of using Gkeyll to support design studies of fusion devices.

Hoffmann, Antoine Cyril David [Princeton Plasma Ph↗
DOE OSTI2024

The effect of plasma shaping on high density H-mode SOL profiles and fluctuations in TCV

Abstract The impact of plasma shaping on the properties of high density H-mode scrape-off layer (SOL) profiles and transport at the outer midplane has been investigated on Tokamakà configuration variable. The experimental dataset has been acquired by evolving the upper triangularity while keeping the other parameters constant. The scan comprises δ u p values between 0.0 and 0.6, excluding negative triangularity scenarios. Within this study, a transition from type-I edge localised modes to the quasi-continuous exhaust regime takes place from low to high δ u p . The modification of the upstream SOL profiles has been assessed, in terms of separatrix quantities, within the α t turbulence control parameter theoretical framework (Eich et al 2020 Nucl. Fusion 60 056016). The target parallel heat load and the upstream near-SOL density profiles have been shown to broaden significantly for increasing α t . Correspondingly, in the far SOL a density shoulder formation is observed when moving from low to high δ u p . These behaviours have been correlated with an enhancement of the SOL fluctuation level, as registered by wall-mounted Langmuir probes as well as the thermal helium beam diagnostic. Specifically, both the background and the filamentary-induced fluctuating parts of the first wall ion saturation current signal are larger at higher δ u p , with filaments being ejected more frequently into the SOL. Comparison of two pulses at the extremes of the δ u p scan range, but with otherwise same input parameters, shows that the midplane neutral pressure does not change much during the H-mode phase of the discharge. This indicates that indirect effects of the change in geometry, linked to first wall recycling sources, should not play a significant role. The total core radiation increases at high δ u p , on account of a stronger plasma–wall interaction and resulting larger carbon impurity intake from the first wall. This is likely associated to the enhanced first wall fluctuations, as well as a smaller outer gap and the close-to-double-null magnetic topology at high shaping.

Stagni, A. (ORCID:0000000180841544)↗
DOE OSTI2026

Global gyro-kinetic ion temperature gradient and trapped electron mode turbulence modelling in X -point geometry in negative and positive triangularity

Comparative modelling of ion temperature gradient (ITG)/trapped electron mode (TEM) turbulence in negative (NT) and positive (PT) triangularity plasma shapes was done using the nonlinear global full-f gyrokinetic particle code JOREK-GK in the realistic X-point tokamak geometry including the Scrape Off Layer (SOL) for TCV and DIII-D parameters. A comparison of JOREK-GK code with the gyrokinetic codes GS2 and GENE was done using NT/PT triangularity TCV L-modes parameters showing good agreement between codes in linear growth rates and clear beneficial effect of NT as compared to PT. Global non-linear modelling of the ITG/TEM saturated turbulence for realistic DIII-D NT pulses was done and compared with numerically constructed PT equilibrium with the same plasma profiles. Existence of longer correlation length of density fluctuations in PT compared to NT was demonstrated. Stronger and more sheared zonal flows are generated via Reynolds stress in NT compared to PT. These factors are stabilizing for TEM/ITG turbulence in NT and lead to smaller heat fluxes and heat conductivities in NT compared to PT configuration. Weak dependence of plasma confinement on collisionality and plasma rotation was found in modelling of DIII-D NT shots similar to the experiment. The confinement scaling with normalized ion gyro-radius ρ* was estimated both for NT and PT. Bohm-like scaling was obtained in both configurations, however with better confinement for NT compared to PT which could be favourable factor for reactor size machines.

ITG↗
DOE OSTI2018

Optimisation of confinement in a fusion reactor using a nonlinear turbulence model

We report the confinement of heat in the core of a magnetic fusion reactor is optimised using a multidimensional optimisation algorithm. For the first time in such a study, the loss of heat due to turbulence is modelled at every stage using first-principles nonlinear simulations which accurately capture the turbulent cascade and large-scale zonal flows. The simulations utilise a novel approach, with gyrofluid treatment of the small-scale drift waves and gyrokinetic treatment of the large-scale zonal flows. A simple near-circular equilibrium with standard parameters is chosen as the initial condition. The figure of merit, fusion power per unit volume, is calculated, and then two control parameters, the elongation and triangularity of the outer flux surface, are varied, with the algorithm seeking to optimise the chosen figure of merit. A twofold increase in the plasma power per unit volume is achieved by moving to higher elongation and strongly negative triangularity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Plasma shape effects on the Alfvén eigenmode spectrum through Alfvén slow-magnetosonic wave coupling

In this report the effect of the plasma shaping (triangularity and elongation) on the continuous spectrum, frequency gaps, and number of eigenmodes is studied. The plasma shaping affects the coupling between Alfvén and slow magnetosonic waves via the geodesic magnetic curvature and the plasma pressure. The Alfvén-slow coupling creates a large number of new gaps in the continuous spectrum below the Alfvén frequency where discrete modes can reside. In ideal magnetohydrodynamic simulations a large number (a few hundred or more) of potential eigenmodes are found. The number of eigenmodes is correlated with the maximum geodesic curvature and a minimum number of possible discrete eigenmodes is found at a negative triangularity of -0.3. It is hypothesized that these possible eigenmodes form a low amplitude and dense discrete spectrum, which can be studied experimentally.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Validation of NSFsim as a Grad-Shafranov equilibrium solver at DIII-D

Plasma shape is a significant factor that must be considered for any Fusion Pilot Plant (FPP) as it has significant consequences for plasma stability and core confinement. A new simulator, NSFsim, has been developed based on a historically successful code, DINA [1], offering tools to simulate both transport and plasma shape. Specifically, NSFsim is a free boundary equilibrium and transport solver and has been configured to match the properties of the DIII-D tokamak. This paper is focused on validating the Grad-Shafranov (GS) solver of NSFsim by analyzing its ability to recreate the plasma shape, the poloidal flux distribution, and the measurements of the simulated diagnostic signals originating from flux loops and magnetic probes in DIII-D. Five different plasma shapes are simulated to show the robustness of NSFsim to different plasma conditions; these shapes are Lower Single Null (LSN), Upper Single Null (USN), Double Null (DN), Inner Wall Limited (IWL), and Negative Triangularity (NT). The NSFsim results are compared against real measured signals, magnetic profile fits from EFIT [2], and another plasma equilibrium simulator, GSevolve [3]. EFIT reconstructions of shots are readily available at DIII-D, but GSevolve was manually ran by us to provide simulation data to compare against.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Microstability of β ~ 1 tokamak equilibria

High-power-density tokamaks offer a potential solution to design cost-effective fusion devices. One way to achieve high power density is to operate at a high β value (the ratio of thermal to magnetic pressure), i.e. β~1 . However, a β~1 state may be unstable to various pressure- and current-driven instabilities or have unfavourable microstability properties. To explore these possibilities, we generate β~1 equilibria and investigate their stability. First, we demonstrate the generation of high-β equilibria with the computer code VMEC. We then analyse these equilibria to determine their stability against the infinite-n ideal-ballooning mode. We follow that by engaging in a detailed microstability study using the GS2 code, beginning with assessments of electrostatic ion-temperature-gradient and trapped election mode instabilities. We observe interesting behaviour for the high-β equilibria – stabilization of these modes through two distinct mechanisms – large negative local shear and reversal of electron precession drift. Finally, we perform electromagnetic gyrokinetic simulations and observe enhanced stability in the outer core of high-β equilibria and absence of kinetic ballooning modes in the negative-triangularity, high-β equilibria. The enhanced outer-core stability of high-β equilibria is different from their lower-β counterparts and offers an alternative, potentially favourable regime of tokamak operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Shaping effects on scrape-off layer plasma turbulence: A rigorous validation of three-dimensional simulations against TCV measurements

The effect of plasma shaping on scrape-off layer (SOL) plasma turbulence is investigated through a rigorous validation exercise. Two- and three-dimensional simulations of the SOL plasma dynamics in three TCV limited discharges are carried out with the GBS code. These discharges realize an almost circular magnetic equilibrium, an elongated equilibrium, and an elongated equilibrium with negative triangularity. For the three plasma discharges, three simulations are performed, considering (i) a three-dimensional model with an explicit dependence on elongation, triangularity, and the inverse aspect ratio; (ii) a circular three-dimensional model in the infinite aspect ratio limit; and (iii) a two-dimensional model, which assumes a circular magnetic equilibrium in the infinite aspect ratio limit, cold ions, and interchange driven turbulence in the sheath limited regime. Ten validation observables common to simulations and experimental measurements from a reciprocating probe located at the TCV outer midplane are identified, and the agreement between experimental and numerical results relative to each observable is evaluated. The composite metric introduced by Ricci et al. [Phys. Plasmas 18, 032109 (2011)] is then used to assess the overall agreement between simulations and experimental measurements. It is found that the shaping model implemented in GBS improves the description of SOL plasma turbulence, taking into account the impact of elongation and triangularity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2024

Kinetic-ballooning-bifurcation in tokamak pedestals across shaping and aspect-ratio

We use a new gyrokinetic threshold model to predict a bifurcation in tokamak pedestal width-height scalings that depends strongly on plasma shaping and aspect-ratio. The bifurcation arises from the first and second stability properties of kinetic-ballooning-modes that yields wide and narrow pedestal branches, expanding the space of accessible pedestal widths and heights. The wide branch offers potential for edge-localized-mode-free pedestals with high core pressure. For negative triangularity, low-aspect-ratio configurations are predicted to give steeper pedestals than conventional-aspect-ratio. Both wide and narrow branches have been attained in tokamak experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

EFIT-Prime: Probabilistic and physics-constrained reduced-order neural network model for equilibrium reconstruction in DIII-D

We introduce EFIT-Prime, a novel machine learning surrogate model for EFIT (Equilibrium FIT) that integrates probabilistic and physics-informed methodologies to overcome typical limitations associated with deterministic and ad hoc neural network architectures. EFIT-Prime utilizes a neural architecture search-based deep ensemble for robust uncertainty quantification, providing scalable and efficient neural architectures that comprehensively quantify both data and model uncertainties. Physically informed by the Grad–Shafranov equation, EFIT-Prime applies a constraint on the current density J tor and a smoothness constraint on the first derivative of the poloidal flux, ensuring physically plausible solutions. Furthermore, the spatial location of the diagnostics is explicitly incorporated in the inputs to account for their spatial correlation. Extensive evaluations demonstrate EFIT-Prime's accuracy and robustness across diverse scenarios, most notably showing good generalization on negative-triangularity discharges that were excluded from training. Timing studies indicate an ensemble inference time of 15 ms for predicting a new equilibrium, offering the possibility of plasma control in real-time, if the model is optimized for speed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A MINi Powder INjector (MINPIN) for the SMART tokamak

Here, this article presents the design and performance characterization of a MINi Powder INjector (MINPIN) developed for boron wall conditioning (boronization) in the SMall Aspect Ratio Tokamak (SMART) at the University of Seville, Spain. SMART studies positive and negative triangularity plasmas within the low aspect ratio range of 1.4 < R/a < 3, with a major radius (R) of 0.40–0.60 m and a minor radius (a) of 0.20–0.33 m. MINPIN, in design, will gravitationally deliver boron (B) powder into the plasma, which ablates and transports the injected material. Mounted on the top of the vessel, the MINPIN will drop controlled amounts of powder using a piezoelectric blade driven by a signal generator that vibrates at a set frequency. The B will then be deposited in a layer on SMART’s stainless steel vacuum vessel walls, and this boronization process reduces the influx of intrinsic wall impurities, such as oxygen, carbon, nitrogen, and iron, which radiate power and cool the plasma. In the laboratory, we calibrated the drop rate for 1–10 mg/s, which will enable controlled boronization coating while achieving efficient powder ablation. We propose MINPIN’s operational scenarios based on the calibrated drop rate and estimated drop duration, assuming the distance from MINPIN to SMART’s plasmas.

Boron powder↗
DOE OSTI2021

Plasma performance and operational space without ELMs in DIII-D

A database of DIII-D plasmas without ELMs compares the operating space and plasma performance of stationary no-ELM regimes: RMP-ELM suppression, QH-mode, I-mode, EDA H-mode, L-mode, and negative triangularity L-mode (Neg-D). Here, operational space is documented in terms of engineering and physics parameters, revealing divergent constraints. Some operational space discriminants (like pedestal collisionality) are well known, while others, such as low torque & safety factor, or high power & density, are less commonly emphasized. Normalized performance (H-factor and betan) also discriminate the no-ELM regimes and favor regimes tolerant to power in DIII-D: RMP, QH, and Neg-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Neural network based fast prediction of β N limits in HL-2M

Artificial neural networks (NNs) are trained, based on the numerical database, to predict the no-wall and ideal-wall β N limits, due to onset of the n = 1 ( n is the toroidal mode number) ideal external kink instability, for the HL-2M tokamak. The database is constructed by toroidal computations utilizing both the equilibrium code CHEASE (Lütjens et al 1992 Comput. Phys. Commun. 69 287) and the stability code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681). The stability results show that (1) the plasma elongation generally enhances both β N limits, for either positive or negative triangularity plasmas; (2) the effect is more pronounced for positive triangularity plasmas; (3) the computed no-wall β N limit linearly scales with the plasma internal inductance, with the proportionality coefficient ranging between 1 and 5 for HL-2M; (4) the no-wall limit substantially decreases with increasing pressure peaking factor. Furthermore, both the NN model and the convolutional neural network (CNN) model are trained and tested, producing consistent results. The trained NNs predict both the no-wall and ideal-wall limits with as high as 95% accuracy, compared to those directly computed by the stability code. Additional test cases, produced by the Tokamak Simulation Code (Jardin et al 1993 Nucl. Fusion 33 371), also show reasonable performance of the trained NNs, with the relative error being within 10%. The constructed database provides effective references for the future HL-2M operations. The trained NNs can be used as a real-time monitor for disruption prevention in the HL-2M experiments, or serve as part of the integrated modeling tools for ideal kink stability analysis.

Physics↗
DOE OSTI2019

Central ion cyclotron emission in the DIII-D tokamak

Collective ion cyclotron emission (ICE) at the ion cyclotron frequency and its harmonics is a potential passive diagnostic of the fast-ion distribution in fusion reactors. ICE is observed in most plasmas in the DIII-D tokamak and is most strongly excited by the fast ions from neutral beam injection. The conventional outboard-edge ICE is detected in H-mode plasmas. However, weaker centrally-localized ICE is measured in L-mode plasmas, including those with negative triangularity shapes. Similar ICE spectra are found with both ICE diagnostics systems, the dedicated magnetic probes and the instrumented antenna straps. Many differences in the behavior of this central ICE are generated by varying the deuterium beam injection angle into deuterium plasmas. The co-current 'near-perpendicular' beam excites the most central ICE from the co-current beams, with this emission detected from the fundamental to the fifth ICE harmonic. However, the counter-current 'near-tangential' beam destabilizes the highest amounts of centrally-localized ICE. This emission is spectrally broader than that driven by the co-current beams and is observed up to its seventh ICE harmonic. The central ICE excited by this co-current beam correlated strongly on the local electron density and related parameters (plasma current and neutron rate) and increased with deeper fast-ion loss boundaries towards the magnetic axis. This was also the case with second harmonic ICE driven by the counter-current beam but not with its stronger third harmonic emission. The central ICE harmonics destabilized by both beams are observed to have different temporal dynamics. The central ICE amplitude responds rapidly to transient MHD events; it dropped and recovered in less than a millisecond at each sawtooth event. ICE frequency splitting is triggered by both the co-current and counter-current 'near-tangential' beams. The data presented provide opportunities to test and validate models of excitation of ICE by energetic ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗