Search NASA⌕ Search

SEARCH · Search NASA

Results for “MHD modeling”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Validation of the imaging neutral particle analyzer in nearly MHD quiescent plasmas using injected beam ions on DIII-D

The imaging neutral particle analyzer (INPA) is a scintillator-based diagnostic that provides energy and radially resolved measurements of confined fast ions. To verify operation, understanding, and modeling of the diagnostic performance, distributions of deuterium beam ions in low density, nearly MHD-quiescent plasmas are measured on the DIII-D tokamak using the INPA. Here, the distribution is altered by changing the pitch-angle scattering rate with electron cyclotron heating (ECH). Experimental images are mostly reproduced with simulation of core fast ions charge exchanging with injected beam and edge cold neutrals in classical cases. Image comparisons between the experiment and simulation show agreement with less than 25% discrepancies near the beam injection energy of 50 keV. Other strong passive signals, observed in a localized region on the image, are suggested to come from charge exchange between cold neutrals and fast ions around the plasma boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Minimizing the Electromechanical Stresses in Poloidal Field Coils by Optimizing their Numbers and Locations using FREDA Framework

Poloidal field (PF) and central solenoid (CS) coils play a crucial role in sustaining the equilibrium and preserving the shape of highly confined tokamak plasmas. Ensuring that PF coil current and mechanical stress stay within superconducting and structural limitations is an important check in the design assessment. Minimizing the PF coil currents and mechanical stresses influences reliability, cost, and performance. A free-boundary MHD equilibrium code—FreeGS is employed within the fusion reactor design and assessment (FREDA) whole facility modeling (WFM) framework to construct the plasma equilibrium based on the configuration and currents in the PF coils. Here, we present the capability of the FreeGS code to minimize the currents, forces, and electromagnetic stresses on the PF coils by optimizing their number, sizes, structures, and locations while maintaining an MHD stable plasma configuration with a large confinement factor. The workflow is initialized with a configuration of plasma parameters and coils’ locations from the 0-D tokamak build systems code in the FREDA framework. Then, FreeGS is called to calculate the initial equilibrium at the minimum total current in PF coils. Thereafter, FreeGS’s internal optimizer minimizes the currents and hoop and central forces on the PF coils while maintaining the reference equilibrium. Finally, the input configuration is updated with the optimized parameters for equilibria over the ramp-up phase of a burning-plasma operation. FREDA’s whole facility optimization capability, which includes all magnetic field coil systems, blanket, vacuum vessel (VV), first wall, divertor, etc., is under development and out of the scope for this study.

Hassan, Ehab [ORNL] (ORCID:0000000181060301)↗
DOE OSTI2022

A finite-volume scheme for modeling compressible magnetohydrodynamic flows at low Mach numbers in stellar interiors

Fully compressible magnetohydrodynamic (MHD) simulations are a fundamental tool for investigating the role of dynamo amplification in the generation of magnetic fields in deep convective layers of stars. The flows that arise in such environments are characterized by low (sonic) Mach numbers ($\mathscr{M}$ son ≲ 10 –2 ). In these regimes, conventional MHD codes typically show excessive dissipation and tend to be inefficient as the Courant–Friedrichs–Lewy (CFL) constraint on the time step becomes too strict. In this work we present a new method for efficiently simulating MHD flows at low Mach numbers in a space-dependent gravitational potential while still retaining all effects of compressibility. The proposed scheme is implemented in the finite-volume SEVEN-LEAGUE HYDRO (SLH) code, and it makes use of a low-Mach version of the five-wave Harten–Lax–van Leer discontinuities (HLLD) solver to reduce numerical dissipation, an implicit–explicit time discretization technique based on Strang splitting to overcome the overly strict CFL constraint, and a well-balancing method that dramatically reduces the magnitude of spatial discretization errors in strongly stratified setups. The solenoidal constraint on the magnetic field is enforced by using a constrained transport method on a staggered grid. We carry out five verification tests, including the simulation of a small-scale dynamo in a star-like environment at $\mathscr{M}$ son ~ 10 –3 . We demonstrate that the proposed scheme can be used to accurately simulate compressible MHD flows in regimes of low Mach numbers and strongly stratified setups even with moderately coarse grids.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

On hydromagnetic wave interactions in collisionless, high- β plasmas

We describe the interaction of parallel-propagating Alfvén waves with ion-acoustic waves and other Alfvén waves, in magnetized, high- β collisionless plasmas. This is accomplished through a combination of analytical theory and numerical fluid simulations of the Chew–Goldberger–Low (CGL) magnetohydrodynamic (MHD) equations closed by Landau-fluid heat fluxes. An asymptotic ordering is employed to simplify the CGL-MHD equations and derive solutions for the deformation of an Alfvén wave that results from its interaction with the pressure anisotropy generated either by an ion-acoustic wave or another, larger-amplitude Alfvén wave. The difference in time scales of acoustic and Alfvénic fluctuations at high- β means that interactions that are local in wavenumber space yield little modification to either mode within the time it takes the acoustic wave to Landau damp away. Instead, order-unity changes in the amplitude of Alfvénic fluctuations can result after interacting with frequency-matched acoustic waves. Additionally, we show that the propagation speed of an Alfvén-wave packet in an otherwise homogeneous background is a function of its self-generated pressure anisotropy. This allows for the eventual interaction of separate co-propagating Alfvén-wave packets of differing amplitudes. The results of the CGL-MHD simulations agree well with these predictions, suggesting that theoretical models relying on the interaction of these modes should be reconsidered in certain astrophysical environments. Applications of these results to weak Alfvénic turbulence and to the interaction between the compressive and Alfvénic cascades in strong, collisionless turbulence are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Ponderomotive electron physics captured in a single-fluid extended magnetohydrodynamics model

The ponderomotive force, arising from the interaction between electromagnetic waves and plasma, plays a critical role in laser fusion, astrophysical plasmas, and laser diagnostics. Traditionally, modeling this force requires multi-fluid or particle-in-cell simulations due to its strong coupling to electron-scale dynamics. In this work, we demonstrate that a one-fluid, two-temperature extended magnetohydrodynamics (XMHD) model—augmented with a generalized Ohm's law (GOL) including electron inertia—can accurately reproduce key ponderomotive effects. We derive the ponderomotive force within this framework using a phasor-based approach and then validate its nonlinear manifestations through direct numerical simulations in the PERSEUS code, where steepening and density modulation phenomena typically associated with kinetic-scale models are reproduced. These results and prior work establish XMHD as a robust and efficient alternative for modeling nonlinear laser–plasma dynamics, bridging the gap between ideal MHD and fully kinetic approaches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Magnetic Reconnection in Solar Flares [Slides]

Magnetic reconnection is the underlying mechanism in solar flares. Magnetic reconnection is an astrophysical process in an electrically conducting plasma where neighboring magnetic field lines running anti-parallel break and reconnection therefore exchanging the direction of their paths. The result of this process is the magnetic field energy converting to kinetic and thermal energy which is explosively released. The scale of a solar is extremely large compared to Earth so kinetic simulations are unable to portray an accurate model. In this project, we are using an astrophysical magnetohydrodynamic (MHD) simulation code known as Athena++ with the reconnection problem generator file created by Li et al. 2018. To further study this simulation, we began testing a heating/cooling term. The results of the original file will be compared with the modified file with an implemented heating/cooling term to better understand how solar flares function.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λ c , that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λ c -based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Investigating nonlinear magnetohydrodynamics in an optimized, reactor-scale quasi-axisymmetric stellarator

We use the new simulation capabilities of the extended-magnetohydrodynamic (MHD) code, M3D-C1, to investigate the nonlinear MHD properties of a reactor-scale quasisymmetric stellarator equilibrium. Our model captures the self-consistent evolution of the magnetic field, temperature, density, and flow profiles without imposing restrictions on the structure of the first. We include the effects of resistivity using a realistic temperature-dependent Spitzer model, along with a model for heat transport that captures the key physical characteristic, namely, strongly anisotropic diffusion in directions perpendicular and parallel to the magnetic field. We consider a quasi-axisymmetric, finite-pressure equilibrium that was optimized for self-consistent bootstrap current, quasi-symmetry, and energetic particle confinement. Our assessment finds that the equilibrium is highly unstable to interchange-like pressure-driven instabilities near the plasma edge. The initially unstable modes rapidly destabilize other modes in the direction of the N-fold rotational symmetry (toroidal, in this case). For this equilibrium, N = 2, meaning destabilization of a large number of even-numbered toroidal Fourier modes. Thus, field-periodicity is likely to be an important factor in the nonlinear MHD stability characteristics of optimized stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)↗
DOE OSTI2024

Extended-MHD simulations of disruption mitigation via massive gas injection in SPARC

Recent developments to the M3D-C1 code enable higher fidelity modeling of disruptions, and can be applied in the design verification of reactor-scale tokamaks. Among these new capabilities is a method to mesh conducting vessel structures such as coils and passive plates, packing of the toroidal mesh around gas injectors, as well as anisotropic resistivity inside the vessel structures. We present extended-magnetohydrodynamic (MHD) simulations of disruption mitigation via massive gas injection (MGI) in SPARC. The goal of this study is to inform the disruption mitigation layout of SPARC and aid in the design of an effective gas injector configuration. Fully three-dimensional simulations with M3D-C1 are carried out for various injector configurations with the primary goal of determining the effect of different MGI parameters on heat loads and vessel forces. The simulations include a model for impurity ionization, recombination, advection and radiation, as well as spatially resolved conducting structures around the plasma. A localized mixture of deuterium and neon with a small toroidal and poloidal width is injected in up to six locations. We demonstrate that M3D-C1 can model a rapid shutdown via MGI using narrow and more realistic gas plumes than in previous simulations. As a result of the q = 1 surface in the SPARC baseline case a sawtooth is observed early in the simulations. Despite the sawtooth and the onset of edge MHD instabilities, the impurity distribution remains localized around the injector locations, but enables a radiative shutdown of the plasma. We find that using the maximum of six gas injectors results in a lower peaking factor and leads to a more even distribution of radiation toroidally than using two injectors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

MARS-F/K modeling of plasma response and fast ion losses due to RMP in KSTAR

Abstract The toroidal single-fluid magnetohydrodynamic (MHD) code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) and the MHD-kinetic hybrid code MARS-K (Liu et al 2008 Phys. Plasmas 15 112503) are utilized to study the plasma response to the n = 1 ( n is the toroidal mode number) resonant magnetic perturbation (RMP), applied to suppress the type-I edge localized mode (ELM) in a KSTAR discharge. Both the resistive-rotating and ideal-static plasma models identify strong screening of the resonant radial field harmonics of the applied RMP due to the plasma response, and predict a strong edge-peeling response of the plasma which is consistent with the optimal ELM control coil current configuration adopted in experiment. The RMP-induced radial displacement of the plasma, computed by the resistive-rotating plasma model, agrees reasonably well with that reconstructed from the measured data in the plasma core. Taking into account the drift kinetic response of fast ions, MARS-K hybrid modeling also finds quantitative agreement of the plasma core fluid pressure perturbation with experiment. Based on the MARS-F computed plasma response, a guiding-center orbit-tracing simulation finds about 0.3% of fast ion losses due the n = 1 RMP in the KSTAR ELM control experiment considered. Most losses are associated with counter-current fast ions located near the plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Experimental, numerical and analytical evaluation of j × B -thrust for fast-liquid-metal-flow divertor systems of nuclear fusion devices

Abstract Divertor systems of fusion devices are exposed to intense heat loads from plasmas, which degrade solid plasma-facing components. Fast liquid metal (LM) flow divertors may be more advantageous for this purpose but have risk of piling due to intense magnetohydrodynamic (MHD) drag. However, severe deceleration of the flow could be countered with the injection of currents that are transverse to external magnetic fields, allowing to thrust the flow with j × B (Lorentz) forces. Given that the injection of currents as an approach to propel LM-divertor flows has remained experimentally understudied, this article focuses on the evaluation of j × B -thrust and finding its drawbacks. j × B -thrust was experimentally tested with free-surface-LM flows, a vertical magnetic field and an externally applied current. Experiments were reviewed with a theoretical model, showing agreement in the trends of theory and experiments. Full 3D-MHD-free-surface-flow simulations were also performed with FreeMHD and confirmed the sensitivity to unstable flow behavior in LM systems when applying external currents. Furthermore, excessive power requirements are expected for the implementation of j × B -thrust at the reactor scale, making these systems inefficient for commercial devices. This paper evidences that the simple operation of a LM-flow divertor with j × B -thrust, without any of the instabilities caused from reactor plasmas or parasitic currents, already presents intrinsic challenges.

Magnetohydrodynamics↗
DOE OSTI2021

Toroidal modeling of runaway avalanche in DIII-D discharges

Abstract A toroidal modeling tool is developed to study the runaway electron (RE) avalanche production process in tokamak plasmas, by coupling the Rosenbluth–Putvinski avalanche model (Rosenbluth and Putvinski 1997 Nucl. Fusion 37 1355) with an n = 0 magneto-hydrodynamic (MHD) solver. Initial value numerical simulations are carried out for two DIII-D discharges with different plasma shapes (one near circular, and the other with high elongation). It is found that, assuming the same level of about 1% seed current level, the Rosenbluth–Putvinski model somewhat underestimates the RE plateau current for the circular-shaped plasma, as compared with that measured in DIII-D experiments. For an elongated, higher current plasma, simulations find strong runaway current avalanche production despite the lack of measured plateau RE current in experiments. A possible reason for this discrepancy is a lack of additional RE dissipation physics in the present two-dimensional model. Systematic scans of the plasma boundary shape, at fixed pre-disruption plasma current, find that the plasma elongation helps to reduce the RE avalanche production, confirming recent results obtained with an analytic model (Fülöp et al 2020 J. Plasma Phys. 86 474860101). The effect of the plasma triangularity (either positive or negative), on the other hand, has a minor effect. On the physics side, the avalanche process involves two competing mechanisms associated with the electric field. On the one hand, a stronger electric field produces a higher instantaneous avalanche growth rate. On the other hand, a fast growing RE current quickly reduces the fraction of the conduction current together with the electric field, and hence a faster dissipation of the poloidal flux. As a final result of these two dynamic processes, the runaway plateau current is not always the largest with the strongest initial electric field. These results lay the foundation for future self-consistent inclusion of the MHD dynamics and the RE amplification process.

Physics↗
DOE OSTI2024

Solving the Orszag–Tang vortex magnetohydrodynamics problem with physics-constrained convolutional neural networks

We study the 2D Orszag–Tang vortex magnetohydrodynamics (MHD) problem through the use of physics-constrained convolutional neural networks (PCNNs) for forecasting the density, ρ, and the magnetic field, B, as well as the prediction of B given the velocity field v of the fluid. In addition to translation equivariance from the convolutional architecture, other physics constraints were embedded: absence of magnetic monopoles, non-negativity of ρ, use of only relevant variables, and the periodic boundary conditions of the problem. The use of only relevant variables and the hard constraint of non-negative ρ were found to facilitate learning greatly. The divergenceless condition ∇·B=0 was implemented as a hard constraint up to machine precision through the use of a magnetic potential to define B=∇×A. Residual networks and data augmentation were also used to improve performance. This allowed for some of the residual models to function as surrogate models and provide reasonably accurate simulations. For the prediction task, the PCNNs were evaluated against a physics-informed neural network, which had the ideal MHD induction equation as a soft constraint. Several models were able to generate highly accurate fields, which are visually almost indistinguishable and have low mean squared error. Only methods with built-in hard constraints produced physical fields with ∇·B=0. The use of PCNNs for MHD has the potential to produce physically consistent real-time simulations to serve as virtual diagnostics in cases where inferences must be made with limited observables.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Kinetic modeling of hot tail runaway electron generation during plasma disruptions using the JOREK code

The generation of runaway electrons (REs) during disruptions poses a significant challenge for the operation of tokamaks. The production of these high-energy electrons can cause substantial damage, particularly when the plasma current is high, making it a critical concern for ITER. For the high-temperature plasmas anticipated in ITER, the primary generation of REs may be dominated by the hot tail mechanism, which consists of the acceleration of hot electrons from the pre-disruption population which have not yet thermalized with the bulk following the rapid cooling of the plasma. To account for the significant 3D effects on RE production, a hot tail modeling framework has been developed within the non-linear 3D extended MHD code JOREK. This paper presents the structure of this framework, which is based on test electrons evolving in MHD fields. The verification of the method shows good agreement with the reference DREAM code for 0D test cases, as well as for axisymmetric simulations of 15 MA ITER H-mode disruption scenarios. Furthermore, a proof-of-principle application to a DIII-D case demonstrates the framework’s capability to capture for the first time the hot tail generation in 3D MHD simulations in realistic geometry. Preliminary results suggest that the production of REs is significantly reduced by stochastic losses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Physics-Based Real-time Analysis and Control to Achieve Transient-Free Operations for the ITER Era (Final Report)

This is the final report for the early career award for Prof. Egemen Kolemen that aimed to develop the physics and technical basis of real-time transient – i.e. disruptions and Edge Localized Modes (ELMs) – avoidance and feedback control at the DIII-D nuclear fusion facility for the ITER era. During this project, we reconstructed high quality real-time DIII-D equilibria using kinetic information from the motional Stark effect (MSE), Thomson Scattering (TS), and Charge Exchange Recombination (CER) spectroscopy. These equilibria were analyzed in real-time to determine the low-n ideal MHD stability and the pedestal stability thresholds, and to obtain plasma response models. These calculations were used to understand which physical process evolves to lead to the onset of disruption in ITER-relevant scenarios at DIII-D, and how disruption can be predicted using stability analysis as opposed to the current methods of correlation of various signals, which are descriptive but lack predictive power. In the disruption studies, we focused on tearing modes, since these constitute the vast majority of disruptions in this regime. With this knowledge, real-time disruption avoidance algorithms were developed and tested. The enhanced situation awareness and physics insight were used in advanced real-time feedback control of pedestal and core profiles, which allowed the pursuit and sustainment of transient-free high-confinement regimes that cannot be achieved via feed-forward design based on in-between shot analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Coupled model for liquid lithium plasma facing components

Numerical analysis provides the design choice and operating window of liquid metal Plasma Facing Components (PFC) concepts. Coupled analysis of boundary plasma together with the surrounding boundary structures is required. Here, to achieve this goal, PPPL is developing a comprehensive multi-physics model for modeling of PFCs in fusion devices. The model includes the fluid-kinetic code SOLPS-ITER and the flow and heat transfer code CFX from ANSYS. SOLPS-ITER was augmented with a liquid metal boundary condition algorithm, allowing direct two-way coupling of the plasma analysis with the two-dimensional analytical slab flow model which includes heat convection in the liquid metal PFC. The target heat flux resulting from this coupled analysis is used as a boundary condition for detailed 3D Computational Fluid Dynamics (CFD) Magneto Hydro Dynamics (MHD) and heat transfer analysis. A new formulation of MHD equations is introduced in the numerical procedure ensuring current conservation of the discretized equations. Results of the 3D analysis are used for final validation of the coupled model. A PFC design where a porous wall is used to stabilize the liquid metal surface, while MHD drive is used to push the liquid metal flow inside the PFC, will be investigated in the regimes where vapor shielding is created for enhanced volumetric plasma heat dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Integrated operation scenarios: Chapter 6 of the special issue: on the path to tokamak burning plasma operation

Here we report the progress of the development and optimization of operational scenarios for ITER and beyond, focusing upon baseline, hybrid, and steady-state scenarios since 2007. This includes advancements made by the integrated operation scenarios (IOS) topical group of the international tokamak physical activity as well as contributions from the broader tokamak community. The key area of research involves developing IOSs that encompass tokamak physics, operation, and technology by utilizing integrated modeling and control strategies. This requires leveraging available actuators to simultaneously control plasma position and shape, MHD activities that could lead to disruptions, transport, plasma-wall interaction and power exhaust, fuel cycle, fusion burn, and tritium breeding. The control extends from the plasma initiation phase, through the current ramp-up, flattop, start and end of the fusion burn, and current ramp-down, to the plasma termination phase. A review of the currently developed scenarios and modeling is provided in terms of (i) optimizing plasma initiation in ITER, (ii) preparing for the low activation phase to fully commission all tokamak systems and establish and validate physics and scenario conditions in preparation for deuterim-tritium (DT) operation, (iii) developing and preparing baseline and hybrid scenarios to demonstrate the feasibility of achieving these regimes within device constraints, (iv) exploring steady-state scenarios to meet ITER’s steady-state goals, (v) evaluating and preparing actuators for ITER, (vi) developing integrated control solutions using shared actuators. The most notable achievements include; (i) the development of ITER demonstration discharges by matching various dimensionless parameters, (ii) the development of scenarios in an ITER-like tungsten environment and DT operation, and (iii) the development of scenarios in superconducting tokamaks, enabling long-pulse operations with similar coil constraints to ITER. Along with these significant achievements, outstanding issues and recommendations for further research and development are provided. Importantly, this study goes beyond simply updating the ITER Physics Basis; it carries profound implications for the broader field of burning plasma research, offering valuable insights and guidance for the next generation of fusion experiments and devices.

ITER↗