Search

NT experiments have demonstrated core performance on par with positive triangularity (PT) H-mode without edge-localized modes (ELMs), encouraging further study of an NT reactor core. In this work, we use integrated modeling to scope the operating space around two NT reactor strategies: a high-field, compact fusion pilot plant concept and a low field, high aspect ratio concept. By integrating equilibrium, core transport, and edge ballooning instability models, we establish a range of operating points with less than 50 MW scrape-off layer power and fusion power comparable to positive triangularity (PT) H-mode reactor concepts. Heating and seeded impurities are leveraged to accomplish the same fusion performance and scrape-off layer exhaust power for various pressure edge boundary conditions. Scans over these pressure edge conditions accommodate any current uncertainty of the properties of the NT edge and show that the performance of an NT reactor will be extremely dependent on the edge pressure. The high-field case is found to enable lower scrape-off layer power because it is capable of reaching high fusion powers at a relatively compact size, which allows increased separatrix density without exceeding the Greenwald density limit. An increase in fusion power density is seen at weaker NT. Infinite-n ballooning instability models indicate that an NT reactor core can reach fusion powers comparable to leading PT H-mode reactor concepts while remaining ballooning-stable. Seeded krypton is leveraged to further lower scrape-off layer power since NT does not have a requirement to remain in H-mode. We contextualize the NT reactor operating space by comparing to popular PT H-mode reactor concepts, and find that NT exhibits competitive ELM-free performance with these concepts for a variety of edge conditions while maintaining relatively low scrape-off layer power., Comment: 23 pages, 14 figures. Submitted to Plasma Physics and Controlled Fusion Special Issue on Advances in the Physics Basis of Negative Triangularity Tokamaks

The first achievement of highly radiating plasmas in negative triangularity is shown with an operational space featuring high core radiation at high Greenwald fraction obtained with the injection of reactor-relevant seeded gases. These negative triangularity (NT) shape diverted discharges reach high values of normalized plasma pressure (BetaN > 2) at high radiation fraction with no ELMs. 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, ExB shear and profiles changes due to impurity radiation cooling. The underlying physics mechanism is robust 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 rector-relevant seeding gases making this regime a potential new scenario for reactor operation.

A plasma scenario with negative triangularity (NT) shaping is achieved on MAST-U for the first time. While edge localized modes (ELMs) are eventually suppressed as the triangularity is decreased below $\delta$ < -0.06, an extended period of H-mode operation with Type-III ELMs is sustained at less negative $\delta$ even through access to the second stability region for ideal ballooning modes is closed. This documents a qualitative difference from the ELM-free access conditions documented in NT scenarios on conventional aspect ratio machines. The electron temperature at the pedestal top drops across the transition to ELM-free operation, but a steady rise in core temperature as $\delta$ is decreased allows for similar normalized beta in the ELM-free NT and H-mode positive triangularity shapes.

This work explores power handling solutions for high-field, highly-radiative negative triangularity (NT) reactors based around the MANTA concept \cite{rutherford_manta_2024}. The divertor design is kept as simple as possible, opting for a standard divertor with standard leg length. FreeGS is used to create an equilibrium for the boundary region, prioritizing a short outer leg length of only $\sim$50 cm ($\sim$40\% of the minor radius). The UEDGE code package is used for the boundary plasma solution, to track plasma temperatures and fluxes to the divertor targets. It is found that for $P_\mathrm{SOL}$ = 25 MW and $n_\mathrm{sep} = 0.96 \times 10^{20}$ m$^{-3}$, conditions consistent with initial core transport modeling, little additional power mitigation is necessary. For external impurity injection of just 0.13\% Ne, the peak heat flux density at the more heavily loaded outer targets falls to 7.8 MW/m$^{2}$, while the electron temperature $T_\mathrm{e}$ remains just under 5 eV. Scans around the parameter space reveal that even at densities lower than in the primary operating scenario, $P_\mathrm{SOL}$ can be increased up to 50 MW, so long as a slightly higher fraction of extrinsic radiator is used. With less than 1\% neon (Ne) impurity content, the divertor still experiences less than 10 MW/m$^{2}$ at the outer target. Design of the plasma-facing components includes a close-fitting vacuum vessel with a tungsten inner surface as well as FLiBe-carrying cooling channels fashioned into the VV wall directly behind the divertor targets. For the seeded heat flux profile, Ansys Fluent heat transfer simulations estimate that the outer target temperature remains at just below 1550\degree C. Initial scoping of advanced divertor designs shows that for an X-divertor, detachment of the outer target becomes much simpler, and plasma fluxes to the targets drop considerably with only 0.01\% Ne content.

The MANTA (Modular Adjustable Negative Triangularity ARC-class) design study investigated how negative-triangularity (NT) may be leveraged in a compact, fusion pilot plant (FPP) to take a ``power-handling first" approach. The result is a pulsed, radiative, ELM-free tokamak that satisfies and exceeds the FPP requirements described in the 2021 National Academies of Sciences, Engineering, and Medicine report ``Bringing Fusion to the U.S. Grid". A self-consistent integrated modeling workflow predicts a fusion power of 450 MW and a plasma gain of 11.5 with only 23.5 MW of power to the scrape-off layer (SOL). This low $P_\text{SOL}$ together with impurity seeding and high density at the separatrix results in a peak heat flux of just 2.8 MW/m$^{2}$. MANTA's high aspect ratio provides space for a large central solenoid (CS), resulting in ${\sim}$15 minute inductive pulses. In spite of the high B fields on the CS and the other REBCO-based magnets, the electromagnetic stresses remain below structural and critical current density limits. Iterative optimization of neutron shielding and tritium breeding blanket yield tritium self-sufficiency with a breeding ratio of 1.15, a blanket power multiplication factor of 1.11, toroidal field coil lifetimes of $3100 \pm 400$ MW-yr, and poloidal field coil lifetimes of at least $890 \pm 40$ MW-yr. Following balance of plant modeling, MANTA is projected to generate 90 MW of net electricity at an electricity gain factor of ${\sim}2.4$. Systems-level economic analysis estimates an overnight cost of US\$3.4 billion, meeting the NASEM FPP requirement that this first-of-a-kind be less than US\$5 billion. The toroidal field coil cost and replacement time are the most critical upfront and lifetime cost drivers, respectively.

Tokamak plasmas with strong negative triangularity (NT) shaping typically exhibit fundamentally different edge behavior than conventional L-mode or H-mode plasmas. Over the entire DIII-D database, plasmas with sufficiently negative triangularity are found to be inherently free of edge localized modes (ELMs), even at injected powers well above the predicted L-H power threshold. A critical triangularly ($\delta_\mathrm{crit}\simeq-0.15$), consistent with inherently ELM-free operation is identified, beyond which access to the second stability region for infinite-$n$ ballooning modes closes on DIII-D. It is also possible to close access to this region, and thereby prevent an H-mode transition, at weaker average triangularities ($\delta\lesssim\delta_\mathrm{crit}$) provided that at least one of the two x-points is still sufficiently negative. Enhanced low field side magnetic fluctuations during ELM-free operation are consistent with additional turbulence limiting the NT edge gradient. Despite the reduced upper limit on the pressure gradient imposed by ballooning stability, NT plasmas are able to support small pedestals and are typically characterized by an enhancement of edge pressure gradients beyond those found in traditional L-mode plasmas. Further, the pressure gradient inside of this small pedestal is unusually steep, allowing access to high core performance that is competitive with other ELM-free regimes previously achieved on DIII-D. Since ELM-free operation in NT is linked directly to the magnetic geometry, NT fusion pilot plants are predicted to maintain advantageous edge conditions even in burning plasma regimes, potentially eliminating reactor core-integration issues caused by ELMs.

A new burn control scheme for tokamaks is described where the total fusion power is controlled by adjusting the plasma volume fraction that is packed into power dense regions. In an example spherical tokamak burning plasma, by modifying the plasma edge squareness the total fusion power is doubled at almost constant total plasma volume and fusion power density. Therefore, increased plasma squareness could be extremely beneficial to a fusion reactor and squareness control could be desirable for power load balancing. Experiments have observed the impact of increased edge squareness on modified core plasma volume, highlighting the practical relevance of this approach., Comment: 7 pages, 7 figures

Negative triangularity (NT) tokamak configurations may be more susceptible to magneto-hydrodynamic instability, posing challenges for recent reactor designs centered around their favorable properties, such as improved confinement and operation free of edge-localized modes. In this work, we assess the vertical stability of plasmas with NT shaping and develop potential reactor solutions. When coupled with a conformal wall, NT equilibria are confirmed to be less vertically stable than equivalent positive triangularity (PT) configurations. Unlike PT, their vertical stability is degraded at higher poloidal beta. Furthermore, improvements in vertical stability at low aspect ratio do not translate to the NT geometry. NT equilibria are stabilized in PT vacuum vessels due to the increased proximity of the plasma and the wall on the outboard side, but this scenario is found to be undesirable due to reduced vertical gaps which give less spatial margin for control recovery. Instead, we demonstrate that informed positioning of passively conducting plates can lead to improved vertical stability in NT configurations on par with stability metrics expected in PT scenarios. An optimal setup for passive plates in highly elongated NT devices is presented, where plates on the outboard side of the device reduce vertical instability growth rates to 16% of their baseline value. For lower target elongations, integration of passive stabilizers with divertor concepts can lead to significant improvements in vertical stability. Plates on the inboard side of the device are also uniquely enabled in NT geometries, providing opportunity for spatial separation of vertical stability coils and passive stabilizers.

A gyrokinetic threshold model for pedestal width-height scaling prediction is applied to multiple devices and to a shaping and aspect-ratio scan giving $\Delta_{\mathrm{ped}} = 0.92 A^{1.04} \kappa^{-1.24} 0.38^{\delta} \beta_{\theta,\mathrm{ped}}^{1.05}$ for pedestal width $\Delta_{\mathrm{ped}}$, aspect-ratio $A$, elongation $\kappa$, triangularity $\delta$, and normalized pedestal height $\beta_{\theta,\mathrm{ped}}$. We also find a width-transport scaling $\Delta_{\mathrm{ped} } = 0.028 \left(q_e/\Gamma_e - 1.7 \right)^{1.5} \sim \eta_e ^{1.5}$ where $q_e$ and $\Gamma_e$ are turbulent electron heat and particle fluxes and $\eta_e = \nabla \ln T_e / \nabla \ln n_e$ for electron temperature $T_e$ and density $n_e$. Pedestals close to those limited by kinetic-ballooning-modes (KBMs) have modified turbulent transport properties compared to strongly driven KBMs. The role of flow shear is studied as a width-height scaling constraint and pedestal saturation mechanism for a standard and wide pedestal discharge., Comment: 34 pages, 16 figures

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., Comment: 7 pages, 7 figures

In this paper, we present a new static and time-dependent MagnetoHydroDynamic (MHD) equilibrium code, TokaMaker, for axisymmetric configurations of magnetized plasmas, based on the well-known Grad-Shafranov equation. This code utilizes finite element methods on an unstructured triangular grid to enable capturing accurate machine geometry and simple mesh generation from engineering-like descriptions of present and future devices. The new code is designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. A detailed description of the numerical methods of the code, including a novel formulation of the boundary conditions for free-boundary equilibria, and validation of the implementation of those methods using both analytic test cases and cross-code validation is shown. Results show expected convergence across tested polynomial orders for analytic and cross-code test cases.

Strongly-shaped diverted negative triangularity (NT) plasmas in the DIII-D tokamak demonstrate simultaneous access to high normalized current, pressure, density, and confinement. NT plasmas are shown to exist across an expansive parameter space compatible with high fusion power production, revealing surprisingly good core stability properties that compare favorably to conventional positive triangularity plasmas in DIII-D. Non-dimensionalizing the operating space, edge safety factors below 3, normalized betas above 3, Greenwald density fractions above 1, and high-confinement mode (H-mode) confinement qualities above 1 are simultaneously observed, all with a robustly stable edge free from deleterious edge-localized mode instabilities. Scaling of the confinement time with engineering parameters reveals at least a linear dependence on plasma current although with significant power degradation, both in excess of expected H-mode scalings. These results increase confidence that NT plasmas are a viable approach to realize fusion power and open directions for future detailed study.

A theoretical model is presented that for the first time matches experimental measurements of the pedestal width-height Diallo scaling in the low-aspect-ratio high-$\beta$ tokamak NSTX. Combining linear gyrokinetics with self-consistent pedestal equilibrium variation, kinetic-ballooning, rather than ideal-ballooning plasma instability, is shown to limit achievable confinement in spherical tokamak pedestals. Simulations are used to find the novel Gyrokinetic Critical Pedestal constraint, which determines the steepest pressure profile a pedestal can sustain subject to gyrokinetic instability. Gyrokinetic width-height scaling expressions for NSTX pedestals with varying density and temperature profiles are obtained. These scalings for spherical tokamaks depart significantly from that of conventional aspect ratio tokamaks., Comment: 7 pages, 5 figures. Nuclear Fusion 2024

In a series of high performance diverted discharges on DIII-D, we demonstrate that strong negative triangularity (NT) shaping robustly suppresses all edge-localized mode (ELM) activity over a wide range of plasma conditions: $\langle n\rangle=0.1-1.5\times10^{20}$m$^{-3}$, $P_\mathrm{aux}=0-15$MW and $|B_\mathrm{t}|=1-2.2$T, corresponding to $P_\mathrm{loss}/P_\mathrm{LH08}\sim8$. The full dataset is consistent with the theoretical prediction that magnetic shear in the NT edge inhibits access to ELMing H-mode regimes; all experimental pressure profiles are found to be at or below the infinite-$n$ ballooning stability limit. Importantly, we also report enhanced edge pressure gradients at strong NT that are significantly steeper than in traditional ELM-free L-mode plasmas and provide significant promise for NT reactor integration., Comment: 5 pages, 5 figures

A new ARC-class, highly-radiative, pulsed, L-mode, burning plasma scenario is developed and evaluated as a candidate for future tokamak reactors. Pulsed inductive operation alleviates the stringent current drive requirements of steady-state reactors, and operation in L-mode affords ELM-free access to $\sim90\%$ core radiation fractions, significantly reducing the divertor power handling requirements. In this configuration the fusion power density can be maximized despite L-mode confinement by utilizing high-field to increase plasma densities and current. This allows us to obtain high gain in robust scenarios in compact devices with $P_\mathrm{fus} > 1000\,$MW despite low confinement. We demonstrate the feasibility of such scenarios here; first by showing that they avoid violating 0-D tokamak limits, and then by performing self-consistent integrated simulations of flattop operation including neoclassical and turbulent transport, magnetic equilibrium, and RF current drive models. Finally we examine the potential effect of introducing negative triangularity with a 0-D model. Our results show high-field radiative pulsed L-mode scenarios are a promising alternative to the typical steady state advanced tokamak scenarios which have dominated tokamak reactor development.

Instability to high toroidal mode number ($n$) ballooning modes has been proposed as the primary gradient-limiting mechanism for tokamak equilibria with negative triangularity ($\delta$) shaping, preventing access to strong H-mode regimes when $\delta\ll0$. To understand how this mechanism extrapolates to reactor conditions, we model the infinite-$n$ ballooning stability as a function of internal profiles and equilibrium shape using a combination of the CHEASE and BALOO codes. While the critical $\delta$ required for avoiding $2^\mathrm{nd}$ stability to high-$n$ modes is observed to depend in a complicated way on various shaping parameters, including the equilibrium aspect ratio, elongation and squareness, equilibria with negative triangularity are robustly prohibited from accessing the $2^\mathrm{nd}$ stability region, offering the prediction that that negative triangularity reactors should maintain L-mode-like operation. In order to access high-$n$ $2^\mathrm{nd}$ stability, the local shear over the entire bad curvature region must be sufficiently negative to overcome curvature destabilization on the low field side. Scalings of the ballooning-limited pedestal height are provided as a function of plasma parameters to aid future scenario design.

This paper introduces the Scrape-Off-Layer Fast Ion (SOLFI) code, which is a new and versatile full-orbit Monte Carlo particle tracer developed to follow fast ion orbits inside and outside the separatrix in tokamaks. SOLFI is benchmarked in a simple straight mirror geometry, showing that the code conserves particle energy and magnetic moment, obtains the correct passing boundary for particles moving in the magnetic mirror field with an imposed electrostatic field, and correctly observes equal ion and electron current at the ambipolar potential predicted from analytical theory. This result has consequences for collisionless scrape-off-layers in spherical tokamaks. We then utilize SOLFI for fundamental physics studies in novel tokamak geometries, exploring the effect of shaping on the trapped particle fraction and bounce locations in tokamaks and demonstrating that negative triangularity can be used to maximize the fraction of particles bouncing in the good-curvature region, potentially leading to enhanced confinement., Comment: 14 pages, 10 figures

Recently, theoretical analysis has identified plasma microinstabilities as the primary mechanism responsible for anomalous heat transport in tokamaks. In particular, the microtearing mode (MTM) has been credited with the production of intense electron heat fluxes, most notably through a thin self-organized boundary layer called the pedestal. Here we exploit a novel, time-dependent analysis to compile explicit experimental evidence that MTMs are active in the pedestal region. The expected frequency of pedestal MTMs, calculated as a function of time from plasma profile measurements, is shown in a dedicated experiment to be in excellent agreement with observed magnetic turbulence fluctuations. Further, fast perturbations of the plasma equilibrium are introduced to decouple the instability drive and resonant location, providing a compelling validation of the analytical model. This analysis offers strong evidence of edge MTMs, validating the existing theoretical work and highlighting the important role of MTMs in regulating electron heat flow in tokamaks., Comment: 8 pages, 5 figures

Experiments searching for neutrinoless double beta decay ($0\nu\beta\beta$) require precise energy calibration and extremely low backgrounds. One of the most popular isotopes for $0\nu\beta\beta$ experiments is $^{136}$Xe. In support of these experiments, the neutron inelastic scattering properties of this isotope have been measured at the GErmanium Array for Neutron Induced Excitations (GEANIE) at the Los Alamos Neutron Science Center. Time-of-flight techniques are utilized with high-purity germanium detectors to search for inelastic scattering $\gamma$ rays for neutron energies between 0.7 and 100 MeV. Limits are set on production of yet-unobserved $\gamma$ rays in the energy range critical for $0\nu\beta\beta$ studies, and measurements are made of multiple $\gamma$-ray production cross sections. In particular, we have measured the production of the 1313 keV $\gamma$ ray which comes from the transition of the first-excited to ground state of $^{136}$Xe. This neutron-induced $\gamma$ line may be useful for a novel energy calibration technique, described in this paper., Comment: 8 pages, 8 figures, submitted to Phys Rev C