SEARCH · Search NASA
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.
A recent publication presents the theoretical basis for extending the Newcomb criterion for the ideal magnetohydrodynamic stability of a straight, circular cylindrical plasma to an axisymmetric toroidal plasma. Here, this paper extends the method to a nonaxisymmetric plasma with stellarator symmetry and toroidal periodicity.
Magnetohydrodynamic stability at the edge pedestal in several quiescent H-mode (QH-mode) plasmas in DIII-D and JT-60U experiments was analyzed by considering plasma rotation and ion diamagnetic drift effects. It was identified that a kink/peeling mode, which is a prime candidate for a trigger of edge harmonic oscillation in QH-mode, is stabilized by plasma rotation when considering the ion diamagnetic drift simultaneously in both experiments. Furthermore, the stabilizing effect by rotation becomes more effective in case using the rotation profile of the main ion species evaluated by assuming radial force balance. In addition, when inverting the rotation direction, it was found that the kink/peeling mode is more stabilized when considering the rotation of the main ion species, though the mode is less stabilized by the rotation of impurity ion species. The result implies that the kink/peeling mode stability in QH-mode plasmas is sensitive to how the rotation profile is evaluated, but it is shown that a qualitative trend stabilizing the kink/peeling mode by rotation can be reproduced with the poloidal rotation profile of an impurity predicted numerically based on the neoclassical theory.
A new optimized quasi-helically symmetric configuration is described that has the desirable properties of improved energetic particle confinement, reduced turbulent transport by three-dimensional shaping and non-resonant divertor capabilities. The configuration presented in this paper is explicitly optimized for quasi-helical symmetry, energetic particle confinement, neoclassical confinement and stability near the axis. Post optimization, the configuration was evaluated for its performance with regard to energetic particle transport, ideal magnetohydrodynamic stability at various values of plasma pressure and ion temperature gradient instability induced turbulent transport. The effects of discrete coils on various confinement figures of merit, including energetic particle confinement, are determined by generating single-filament coils for the configuration. Preliminary divertor analysis shows that coils can be created that do not interfere with expansion of the vessel volume near the regions of outgoing heat flux, thus demonstrating the possibility of operating a non-resonant divertor.
Magnetic feedback control of the resistive-wall mode has enabled DIII-D to access stable operation at safety factor q95 = 1:9 in divertor plasmas for 150 instability growth times. Magnetohydrodynamic stability sets a hard, disruptive limit on the minimum edge safety factor achievable in a tokamak, or on the maximum plasma current at given toroidal magnetic eld. In tokamaks with a divertor, the limit occurs at q95 = 2, as con rmed in DIII-D. Since the energy con cement time scales linearly with current, this also bounds the performance of a fusion reactor. DIII-D has overcome this limit, opening a whole new high-current regime not accessible before. This result brings signi cant possible bene ts in terms of fusion performance, but it also extends resistive wall mode physics and its control to conditions never explored before. In present experiments, q95 < 2 operation is eventually halted by voltage limits reached in the feedback power supplies, not by intrinsic physics issues. Improvements to power supplies and to control algorithms have the potential to further extend this regime.
Here, the advanced physics and advanced technology tokamak power plant ARIES-ACT1 has a major radius of 6.25 m at an aspect ratio of 4.0, toroidal field of 6.0 T, strong shaping with elongation of 2.2, and triangularity of 0.63. The broadest pressure cases reached wall-stabilized β N ~ 5.75, limited by n = 3 external kink mode requiring a conducting shell at b/a = 0.3, requiring plasma rotation, feedback, and/or kinetic stabilization. The medium pressure peaking case reaches β N = 5.28 with B T = 6.75, while the peaked pressure case reaches β N < 5.15. Fast particle magnetohydrodynamic stability shows that the alpha particles are unstable, but this leads to redistribution to larger minor radius rather than loss from the plasma. Edge and divertor plasma modeling shows that 75% of the power to the divertor can be radiated with an ITER-like divertor geometry, while >95% can be radiated in a stable detached mode with an orthogonal target and wide slot geometry. The bootstrap current fraction is 91% with a q95 of 4.5, requiring ~1.1 MA of external current drive. This current is supplied with 5 MW of ion cyclotron radio frequency/fast wave and 40 MW of lower hybrid current drive. Electron cyclotron is most effective for safety factor control over ρ~0.2 to 0.6 with 20 MW. The pedestal density is ~0.9×10 20 /m 3 , and the temperature is ~4.4 keV. The H98 factor is 1.65, n/n Gr = 1.0, and the ratio of net power to threshold power is 2.8 to 3.0 in the flattop.
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.
A general stability condition for plasma-vacuum systems with resistive walls is derived by using the Frieman Rotenberg Lagrangian stability formulation [Rev. Mod. Phys. 32, 898 (1960)]. It is shown that the Lyapunov stability limit for external modes does not depend upon the gyroscopic term but upon the sign of the perturbed potential energy only. In the absence of dissipation in the plasma such as viscosity, it is expected that the flow cannot stabilize the system.
High fidelity kinetic equilibria are crucial for tokamak modeling and analysis. Manual workflows for constructing kinetic equilibria are time consuming and subject to user error, motivating development of automated equilibrium reconstruction tools to provide accurate and consistent reconstructions for downstream physics analysis. These automated tools also provide access to kinetic equilibria at large database scales, which enables the quantification of general uncertainties arising from equilibrium reconstruction techniques. In this paper, we compare a large database of DIII-D kinetic equilibria generated manually by physics experts to equilibria from automated kinetic reconstruction tools, assessing the impact of reconstruction method on equilibrium parameters and resulting magnetohydrodynamic stability calculations. We find agreement among scalar parameters, whereas profile quantities, such as the bootstrap current, show larger disagreements. We analyze ideal kink and classical tearing stability with DCON and STRIDE respectively, finding that the kink stability calculation is generally more robust than the tearing index Δ' calculation. We find that in 90% of cases, both kink stability classifications are unchanged between the manual expert and automated kinetic equilibria.
A matching method using a finite-width inner region is extended for stability analysis of magnetohydrodynamic mode including diamagnetic drift effect. The inclusion of the diamagnetic drift effect is accomplished by a newly developed ordering scheme in the outer region. The ordering scheme enables us to derive a hierarchy of generalized Newcomb equations. Higher-order equations give us correction of outer solution due to the diamagnetic drift effect as well as inertia and resistivity. By this correction, the accuracy of the dispersion relation is improved. Several numerical results are presented to demonstrate good performance of the matching method. Dropping the diamagnetic drift effect in the outer region leads to less accurate results.
A method is demonstrated to optimize a stellarator's geometry to eliminate magnetic islands and achieve other desired physics properties at the same time. For many physics quantities that have been used in stellarator optimization, including quasisymmetry, neoclassical transport, and magnetohydrodynamic stability, it is convenient to use a magnetic equilibrium representation that assures the existence of magnetic surfaces. However, this representation hides the possible presence of magnetic islands, which are typically undesirable. To include both surface-based objectives and island widths in a single optimization, two fixed-boundary equilibrium calculations are run at each iteration of the optimization: one that enforces the existence of magnetic surfaces (the Variational Moments Equilibrium Code) [S. P. Hirshman and J. C. Whitson, Phys. Fluids 26, 3553 (1983)] and one that does not (the Stepped Pressure Equilibrium Code) [Hudson et al., Phys. Plasmas 19, 112502 (2012)]. By penalizing the island residues in the objective function, the two magnetic field representations are brought into agreement during the optimization. Here, an example is presented in which, particularly on the surface where quasisymmetry was targeted, quasisymmetry is achieved more accurately than in previously published examples.
Integrated control of the toroidal current density profile, or alternatively the q-profile, and plasma stored energy is essential to achieve advanced plasma scenarios characterized by high plasma confinement, magnetohydrodynamics stability, and noninductively driven plasma current. The q-profile evolution is closely related to the evolution of the poloidal magnetic flux profile, whose dynamics is modeled by a nonlinear partial differential equation (PDE) referred to as the magnetic-flux diffusion equation (MDE). The MDE prediction depends heavily on the chosen models for the electron temperature, plasma resistivity, and non-inductive current drives. To aid control synthesis, control-oriented models for these plasma quantities are necessary to make the problem tractable. However, a relatively large deviation between the predictions by these control-oriented models and experimental data is not uncommon. For this reason, the electron temperature, plasma resistivity, and non-inductive current drives are modeled for control synthesis in this work as the product of an “uncertain” reference profile and a nonlinear function of the different auxiliary heating and current-drive (H&CD) source powers and the total plasma current. The uncertainties are quantified in such a way that the family of models arising from the modeling process is able to capture the q-profile and plasma stored energy dynamics from a typical EAST shot. A control-oriented nonlinear PDE model is developed by combining the MDE with the “uncertain” models for the electron temperature, plasma resistivity, and non-inductive current drives. This model is then rewritten into a control framework to design a controller that is robust against the modeled uncertainties. The resulting controller utilizes EAST's H&CD powers and total plasma current to regulate the q profile and plasma stored energy even when mismatches between modeled and actual dynamics are present. The effectiveness of the controller is demonstrated through nonlinear simulations.
Commonwealth Fusion Systems plans to build ARC as the first fusion power plant at a site in Chesterfield County, Virginia, USA by the early 2030s. We present an overview of analysis comprising the physics basis of the ARC V3A design, a high-magnetic-field tokamak with 𝐵 0 =11.4 T, 𝐼 𝑝 =12.0 MA, 𝑅 0 =4.62 m, 𝑎 =1.18 m. ARC V3A is designed to produce 𝑃 𝑓𝑢𝑠 ≈1.13 GW DT fusion power and deliver ⩾400 MW net electric power to the grid. This overview includes quantitative analysis of fundamental issues for design of and operational plasma scenarios for a tokamak power plant, and lays out the design targets and strategic choices for ARC, including empirical fusion performance projections, assessment of H-mode access, ion cyclotron resonance heating simulations, alpha particle physics and time-dependent full-pulse simulations. This is complemented by topical papers on fusion performance and transport, disruption physics, boundary physics and magnetohydrodynamic stability. Critically, these studies identify key model uncertainties and physics risks to be retired through SPARC operation. Due to the modular nature of ARC, early results from SPARC can be incorporated into the design of the first ARC as well as subsequent replacements of the ARC vacuum vessel.
Scientific applications of artificial intelligence (AI) often remain limited by device-specific training and unexplained “black-box” approaches, creating fundamental barriers to cross-system generalization. This challenge is critical for nuclear fusion, where future reactors will have limited operational data for AI training. Here, we demonstrate that our neural network, trained solely on megahertz-scale turbulence measurements from one machine (DIII-D), forecasts Type-I edge localized mode (ELM) onsets in a different tokamak (KSTAR) through zero-shot weight transfer following physics-consistent preprocessing without device-specific retraining. Through an explainable AI framework combining gradient-weighted class activation mapping with physics validation, we reveal that our network can internalize physics relationships governing the ELM instabilities rather than memorizing device-specific patterns. The network perceives spatiotemporal features that correlate consistently with independently calculated instability growth rates, magnetohydrodynamic stability limits, and pedestal structure dynamics. Statistical analyses of dimensionally-reduced saliency features reveal the identical triangular features between the saliency representations, instability growth rates, and prediction probability across tokamaks, providing evidence that our forecasting system can show tokamak-agnostic generalization. This work contributes to a foundation for explainable scientific AI systems, where cross-system developments are essential for transcending traditional domain-specific constraints.
The problem of ideal magnetohydrodynamic stability of plasmas in a magnetosphere-atmosphere system, in which the unperturbed magnetic field is assumed to be perpendicular to the plasma-atmosphere interface (ionospheric surface), is investigated by means of an extended magnetospheric energy principle. The derivation of the principle and conditions under which it applies to a real terrestrial magnetosphere is given. In the principle, the atmosphere is considered to be a very heavy and compressible gas with finite pressure. A thin ionospheric layer is taken into account as boundary conditions, but energetics within it are neglected. The solid-earth surface is assumed to be a perfectly conducting wall for perturbations. For a perturbation that satisfies either rigid or horizontally free boundary conditions at the plasma-atmosphere interface, the self-adjointness of the force operator is satisfied and an extended magnetospheric energy principle can be developed on the basis of the extended energy principle for fusion plasmas. These two boundary conditions are shown to be realized in the magnetosphere when the ionospheric conductivity is either very large or very small. Whereas in fusion plasmas the perturbed magnetic energy in the vacuum makes a stabilizing contribution to the potential energy, in the magnetosphere the perturbed magnetic energy in the atmosphere makes no such stabilizing contribution. This is due to the difference of the assumed field configurations of the magnetospheric and fusion plasmas. The ionospheric surface makes a destabilizing negative contribution to the potential energy owing to a horizontal plasma displacement on the spherical ionospheric surface. The method is applied to magnetospheric ballooning and interchange instabilities. The existence of a new type of magnetospheric interchange instability is shown and its structure in the magnetosphere-atmosphere system is clarified. Possible consequences of the instabilities and their relevance to magnetospheric physics are discussed.
The generalized conditional symmetry (GCS) method is applied to a specific case of the Grad–Shafranov (GS) equation, in cylindrical geometry assuming the existence of an axial symmetry. We investigate the conditions that yield the GS equation admitting a special class of second-order GCSs. The determining system for the unknown arbitrary functions is solved in several special cases and new exact solutions, including solitary waves, different in form and structure from the ones obtained using other nonclassical symmetry methods, are pointed out. Several plots of the level sets or flux surfaces of the new solutions as well as surfaces with vanishing flow are displayed. The obtained solutions can be useful for studying plasma equilibrium, transport phenomena, and magnetohydrodynamic stability.
The extended magnetohydrodynamic stability of interchange modes is studied in two configurations. In slab geometry, a local dispersion relation for the gravitational interchange mode (g-mode) with three different extensions of the MHD model [P. Zhu, et al., Phys. Rev. Lett. 101, 085005 (2008)] is analyzed. Our results delineate where drifts stablize the g-mode with gyroviscosity alone and with a two-fluid Ohm’s law alone. Including the two-fluid Ohm’s law produces an ion drift wave that interacts with the g-mode. This interaction then gives rise to a second instability at finite k y . A second instability is also observed in numerical extended MHD computations of linear interchange in cylindrical screw-pinch equilibria, the second configuration. Particularly with incomplete models, this mode limits the regions of stability for physically realistic conditions. But, applying a consistent two-temperature extended MHD model that includes the diamagnetic heat flux density ($\vec{q}$ * ) makes the onset of the second mode occur at larger Hall parameter. For conditions relevant to the SSPX experiment [E.B. Hooper, Plasma Phys. Controlled Fusion 54, 113001 (2012)], significant stabilization is observed for Suydam parameters as large as unity (D s ≲1).
The use of supersonic rotation of a plasma in mirror geometry has distinct advantages for thermonuclear fusion. The device is steady state, there are no disruptions, the loss cone is almost closed, sheared rotation stabilizes magnetohydrodynamic instabilities as well as plasma turbulence, there are no runaway electrons, and the coil configuration is simple. In this work, we examine the effect of rotation on mirror confinement using a full cyclotron orbit code. The full cyclotron simulations give a much more complete description of the particle energy distribution and losses than the use of guiding center equations. Both collisionless loss as a function of rotation and the effect of collisions are investigated. Although the cross field diffusion is classical, we find that the local rotating Maxwellian is increased to higher energy, increasing the fusion rate and also enhancing the radial diffusion. We find a loss channel not envisioned with a guiding center treatment, but a design can be chosen that can satisfy the Lawson criterion for ions. Of course, the rotation has a minimal effect on the alpha particle birth distribution, so there is initially loss through the usual loss cone, just as in a mirror with no rotation. However after this loss, the alphas slow down on the electrons with little pitch angle scattering until reaching low energy, so over half of the initial alpha energy is transferred to the electrons. The important problem of energy confinement, with losses primarily through the electron channel, is not addressed in this work. In conclusion, we also discuss the use of rotating mirror geometry to produce an ion thruster.
Fast ion (FI) loss properties in the presence of tearing mode and internal kink perturbations are numerically investigated for discharges in the MAST-U spherical tokamak, utilizing the MARS-F magnetohydrodynamic stability code and the REORBIT test particle guiding-center orbit-following module. Here, tracing about 100 000 particle markers sampled from the equilibrium distribution of the neutral-beam injection induced FIs, it is found that about 10% out of the total strike the limiting surface (including the divertor surface) in MAST-U discharge 46943, assuming a maximum perturbation of 100 G inside the plasma (corresponding to ~6 G at the Mirnov probe location at the outboard mid-plane). Detailed particle tracing, assuming a uniform initial distribution in the 2D phase space (at given radial locations), reveals that initially counter-current FIs launched near the plasma edge are subject to significant prompt losses, while almost all initially co-current ions remain well confined at the assumed perturbation level. Most lost FIs strike the lower-half of the limiting surface. Finite gyro-radius effects prevent lost ions from striking the top-outer corner of the super-X divertor chamber. A scan of the perturbation level (based on discharge 45163) reveals, not surprisingly, an approximately linear scaling of the particle loss fraction (for counter-current FIs) with respect to the perturbation amplitude.