SEARCH · Engineering Papers
Results for “TOKAMAK”
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.
Integration of the ITER Tokamak C-model with the Building Model
The pellet injector system is designed to accommodate the requirement of continuous plasma refueling during ITER Tokamak operations. Through pellet flight tubes/injection lines, the pellets can directly reach the interior of the tokamak vacuum vessel. The pellet injector system is located at a port cell outside of the Tokamak bio-shield. The plasma radiation source in the tokamak vacuum vessel, together with the highly activated surrounding water in port cells, present major radiation challenges for the pellet injector system during ITER operation and after shutdown. These sources produce a radiation field encompassing the entire Tokamak complex. The ITER neutronics reference model C-model is an ITER radiation transport model. It represents a 40 deg. regular torus sector out to the bio-shield, which consists of 38 individual systems. The Monte Carlo N-Particle (MCNP) input of the C-model consists of more than a hundred thousand cell definitions and surface definitions, and over three hundred material compositions. The latest tokamak complex MCNP model is significantly more detailed than any previous releases. It consists of 36,863 cell definitions, 57,105 surface definitions, and a total of 21 different types of materials. The C-model and tokamak building model share the bio-shield as the boundary between them. In general, the radiation fields in the building due to the plasma source during tokamak operation are determined by generating the radiation source in the inner faces of the bio-shield. An inherent limitation of this approach is that it can only approximate the directions in which streaming paths are produced through the bio-shield plugs. To better estimate the shutdown dose rate for the pellet injector system due to the plasma source, we have decided to integrate the tokamak model (C-model) with the complex building model. Both MCNP models are provided by the ITER International Organization (IO) to the Oak Ridge National Laboratory (ORNL) Radiation Transport (RT) group.
Long plasma duration operation analyses with an international multi-machine (tokamaks and stellarators) database
Combined high-fusion performance and long-pulse operation is one of the key integration challenges for fusion energy development in magnetic devices. Addressing these challenges requires an integrated vision of physics and engineering aspects with the purpose of simultaneously increasing time duration and fusion performance. Significant progress has been made in tokamaks and stellarators, including very recent achievement in duration and/or performance. This progress is reviewed by analyzing the experimental data (109 plasma pulses with a total of 3200 data points, i.e. on average 29 data per pulse) provided by ten tokamaks (in alphabetical order: Axially Symmetric Divertor Experiment Upgrade, DIII-D, Experimental Advanced Superconducting Tokamak, Joint European Torus, JT-60 Upgrade, Korea Superconducting Tokamak Advanced Research, tokamak à configuration variable, Tokamak Fusion Test Reactor, Tore Supra, W Environment in Steady-State Tokamak) and two stellarators (Large Helical Device and Wendelstein 7-X) expanding the pioneering work of Kikuchi (Kikuchi M. and Azumi M. 2015 Frontiers in Fusion Research II: Introduction to Modern Tokamak Physics (Springer)). Data have been gathered up to January 2022 and coordination has been provided by the recently created International Energy Agency-International Atomic Energy Agency international Coordination on International Challenges on Long duration OPeration group. By exploiting the multi-machine international database, recent progress in terms of injected energies (e.g. 1730 MJ in L-mode, 425 MJ in H-mode), durations (1056 s in L-mode, 101 s in H-mode), injected powers, and sustained performance will be reviewed. Progress has been made to sustain long-pulse operation in tokamaks and stellarators with superconducting coils, actively cooled components, and/or with metallic walls. The graph of the fusion triple products as a function of duration shows a dramatic reduction of, at least two orders of magnitude when increasing the plasma duration from less than 1 s to 100 s. Indeed, long-pulse operation is usually reached in dominant electron-heating modes at reduced density (current drive optimization) but with low ion temperatures ranging from 1 to 3 keV for discharges above 100 s. Difficulties in extending the duration may arise from coupling high heating powers over long durations and the evolving plasma-wall interaction towards an unstable operational domain. Possible causes limiting the duration and critical issues to be addressed prior to ITER operation and DEMO design are reported and analyzed.
Tokamak Energy’s pre-concept design for a fusion power plant: an overview of ST-E1
Climate change and rapidly rising energy demand, driven in part by artificial intelligence and data-centre growth, create an urgent need for stable, low-carbon, and abundant power. Fusion is a promising long-term solution, yet its commercialisation faces a fundamental paradox in today’s investment environment: pilot plants are essential to de-risk physics, engineering, and operations, but their limited lifetime energy output and high upfront costs make them difficult to finance. This paper presents Tokamak Energy’s response: ST-E1, a pre-concept design for a low-aspect-ratio tokamak power plant engineered specifically to overcome this challenge. ST-E1 is designed from the outset for phased operation—pilot and commercial phases, with an upgrade phase in between—with emphasis on commercial viability, maintainability, nuclear engineering, modularity, and upgradability. A key design principle is the deliberate separation of long-lived assets, such as the magnet cage and vacuum vessel, from replaceable in-vessel systems. This provides an attractive and credible investment approach to generate operational data and de-risk key technologies while preserving most capital-intensive assets for later commercial phases. The architecture supports continuous optimisation toward high net electric power (targeting 800–1000 MW net electric), a normalised capital expenditure of $\$$ 12–14k/kW of net electric power, and high availability (targeting > 80%). A tokamak core with a 5 m major radius, aspect ratio of 2.3, and on-plasma axis toroidal field of 5.25 T was selected to meet these objectives. This paper summarises the ST-E1 design philosophy, principal features, and development methodology. It introduces a Focus Collection of 11 papers detailing the pre-concept design of the entire tokamak and corresponding plant.
Integration of the Full Tokamak Reference Model with the Complex Model for ITER Neutronic Analysis
The ITER International Organization has developed a number of reference Monte Carlo N-Particle (MCNP) models including the tokamak machine C-model, the Tokamak Complex model, and the neutral beam injection (NBI) systems model. The Tokamak Complex model primarily describes building structures beyond the bioshield. Representation of the tokamak and its systems are not included in this model. The Oak Ridge National Laboratory Radiation Transport Group has conducted two ITER neutronic analysis model integrations: (1) integration of the tokamak C-model with the Tokamak Complex model for shutdown dose rate characterization in Port Cell 16 at level B1, and (2) integration of the NBI model with the Tokamak Complex model for estimating the spatial distribution of biological dose rate at levels L1, L2, and L3 of the Tokamak Complex. The integrated models were further extended to include models of system components that are essential to the neutronic analyses. Furthermore, this paper presents the approach and computer tools used to integrate existing reference models, describes the additional design details implemented in the integrated models, and provides representative neutronic calculations based on the extended models.
Disruption Mitigation Solutions for Long Pulse Tokamaks (Final Technical Report)
The project "Disruption Mitigation Solutions for Long Pulse Tokamaks" was a collaborative effort between ORNL, GA, and UCSD. The main goal of the project was to study the injection of shattered cryogenic pellets (SPI) as a method of reducing wall damage from global tokamak plasma instabilities called disruptions. This work was primarily carried out in the JET tokamak in the UK. The team already had extensive experience studying PSI in the medium sized DIII-D tokamak in the US, but adding data from the larger JET tokamak enabled more confident extrapolation to the even larger planned ITER tokamak in France. Each of the collaborating institutions focused on different aspects of SPI disruption mitigation. The UCSD collaboration focused on the volume recombination of post-disruption runaway electron (RE) beams by massive hydrogen injection. Disruptions can sometimes result in the formation of toroidal (ring shaped) beams of highly energetic (relativistic) electron beams, which can cause serious local wall damage, even in present-day tokamaks. A promising method to reduce wall damage from the RE beams is to inject massive amounts of H 2 or D 2 into the beams using SPI. although the precise method by which the injected H 2 reducses RE wall damage is still being investigated, it seems clear that a necessary and required step for the reduction to occur is for the H 2 injection to cause volume recombination (free electron density reduction) of the background cold plasma surrounding the fast electron beam. The UCSD collaboration focused on answering the question of if this recombination could be achievable in the future reactor-size tokamak ITER using realistic levels of H 2 injection. A model capturing the most essential physics was developed and was constrained using both DIII-D and JET data. This model was then used to extrapolate to ITER. It was predicted that RE beam recombination should be achievable in ITER, although in the worst-case scenari (a 10 MA Re beam in ITER), recombination is barely achieved within the maximum allowed levels of H 2 injection.
Plasma self-driven current in tokamaks with magnetic islands
Magnetic island perturbations may cause a reduction in plasma self-driven current that is needed for tokamak operation. A novel effect on tokamak self-driven current revealed by global gyrokinetic simulations is due to magnetic-island-induced 3D electric potential structures, which have the same dominant mode numbers as that of the magnetic island, whereas centered at both the inner and outer edge of the island. The non-resonant potential islands are shown to drive a current through an efficient nonlinear parallel acceleration of electrons. In large aspect ratio (large-A) tokamak devices, this new effect can result in a significant global reduction of the electron bootstrap current when the island size is sufficiently large, in addition to the local current loss across the island region due to the pressure profile flattening. It is shown that there exists a critical magnetic island width for large-A tokamaks beyond which the electron bootstrap current loss is global and increases rapidly with the island size. As such, this process may introduce a size limit for tolerable magnetic islands in large-A tokamak devices in the context of steady state operation. On the other hand, the current loss caused by magnetic islands in low-A tokamaks such as spherical tokamak (ST) NSTX/U is minor. The reduction of the axisymmetric current by magnetic islands scales with the square of island width. However, the loss of the current is mainly local to the island region, and the pace of current loss as the island size increases is substantially slower compared to large-A tokamaks. In particular, the bootstrap current reduction in STs is even smaller in the reactor-relevant high-β p regime where neoclassical tearing modes are more likely to develop.
Computational analysis of ion orbital loss in diverted positive- and negative-triangularity tokamaks
The mechanism of ion orbital loss in diverted negative-triangularity (NT) tokamaks is significantly different from that of positive-triangularity (PT) tokamaks. The orbital loss from the edge of NT tokamaks can be quite large, because the separatrix X-point is located on the outboard side, meaning that trapped ions crossing the separatrix are lost to the inboard divertor plate after a mirror reflection on the side of the strong magnetic field. NT tokamaks have difficulties in building up steep pressure gradients in the edge regions. However, the majority of trapped ions crossing the separatrix in PT tokamaks survive (through complete banana motions) without reaching the X-point region. For NT tokamaks, topological changes in the trapped ion trajectories can be achieved in the presence of the radial electric field, thus reducing the orbital loss, whereas for PT tokamaks, the orbital loss is relatively insensitive to the radial electric field.
Scenario adaptive disruption prediction study for next generation burning-plasma tokamaks
Next generation High Performance (HP) tokamaks risk damage from unmitigated disruptions at high current and power. Achieving reliable disruption prediction for a device's HP operation based on its Low Performance (LP) data is a key to its success. In this letter, through explorative data analysis and dedicated numerical experiments on multiple existing tokamaks, we demonstrate how the operational regimes of tokamaks can affect the power of a trained disruption predictor. First, our results suggest data-driven disruption predictors trained on abundant LP discharges work poorly on the HP regime of the same tokamak, which is a consequence of the distinct distributions of the tightly correlated signals related to disruptions in these two regimes. Second, we find that matching operational parameters among tokamaks strongly improves cross-machine accuracy which implies our model learns from the underlying scalings of dimensionless physics parameters like q 95 , β p and confirms the importance of these parameters in disruption physics and cross machine domain matching from the data-driven perspective. Finally, our results show in the absence of HP data from the target devices, the best predictivity of the HP regime for the target machine can be achieved by combining LP data from the target with HP data from other machines. Furthermore, these results provide a possible disruption predictor development strategy for next generation tokamaks, such as ITER and SPARC, and highlight the importance of developing baseline scenario discharges of future tokamaks on existing machines to collect more relevant disruptive data.
Progress in pedestal and edge physics: Chapter 3 of the special issue: on the path to tokamak burning plasma operation
This paper describes the extensive progress that has been made in the understanding of tokamak pedestal physics since the 2007 publication of ‘Progress in the ITER Physics Basis’ (Ikeda 2007 Nucl. Fusion 47 E01–S500). It serves as Chapter 3 of the 2025 Nuclear Fusion Special Issue titled ‘On the Path to Tokamak Burning Plasma Operation’ (Campbell et al 2025 Nucl. Fusion ). This review was compiled by the pedestal and edge physics (PEP) community affiliated with the International Tokamak Physics Activity organization. It attempts to collect in one place citations to the majority of published literature on the pedestal physics topics that will be most important for the operation of a future power producing burning plasma tokamak. These include citations to publications describing the physics of the pedestal plasmas in many operating tokamaks worldwide and the pedestal physics projections for several near-term future devices including ITER. Descriptions of experimental results, interpretive modeling and predictive extrapolations are integrated together and comprehensive references are provided. This review is organized around four primary technical sections, viz.: pedestal structure, edge localized mode (ELM) characteristics, ELM control and regimes without large ELMs. Key results from many of the references are described briefly and set into the tokamak burning plasma power plant context. In addition, different perspectives on pedestal physics topics that are currently under debate within the community are also described, to provide guidance on needs for future research. Finally, attempts are made to describe conclusions from all of this progress consistent with discussions by the pedestal physics community at this time. The goal of this review is to provide a useful reference document for pedestal physics researchers going forward toward operation of a burning tokamak fusion plasma.
Basic Physics of Tokamak Transport (Final Technical Report)
The goal of this grant has been to study the basic physics of various sources of anomalous transport in tokamaks. Anomalous transport in tokamaks continues to be one of the major problems in magnetic fusion research. As a tokamak is not a physics device by design, direct experimental observation and identification of the instabilities responsible for transport, as well as physics studies of the transport in tokamaks, have been difficult and of limited value. It is noted that direct experimental observation, identification and physics study of microinstabilities including ITG, ETG, and trapped electron/ion modes in tokamaks has been very difficult and nearly impossible. The primary reasons are co-existence of many instabilities, their broadband fluctuation spectra, lack of flexibility for parameter scans and absence of good local diagnostics. This has motivated us to study the suspected tokamak instabilities and their transport consequences in a simpler, steady state Columbia Linear Machine (CLM) with collisionless plasma and the flexibility of wide parameter variations. Earlier work as part of this grant was focused on both ITG turbulence, widely believed to be a primary source of ion thermal transport in tokamaks, and the effects of isotope scaling on transport levels. Prior work from our research team has produced and definitively identified both the slab and toroidal branches of this instability and determined the physics criteria for their existence. All the experimentally observed linear physics corroborate well with theoretical predictions. However, one of the large areas of research dealt with turbulent transport results that indicate some significant differences between our experimental results and most theoretical predictions.
Effect of 3D magnetic perturbations on divertor conditions and detachment in tokamak and stellarator
Enhanced perpendicular heat and momentum transport induces parallel pressure loss leading to divertor detachment, which can be produced by the increase of density in 2D tokamaks. However, in the 3D configurations such as tokamaks with 3D fields and stellarators, the fraction of perpendicular transport can be higher even in a lower density regime, which could lead to the early transition to detachment without passing through the high-recycling regime. 3D fields applied to the limiter tokamak plasmas produce edge stochastic layers close to the last closed flux surface (LCFS), which can allow for enhanced perpendicular transport and indeed the absence of high recycling regime and early detachment have been observed in TEXTOR and Tore Supra. However, in the X-point divertor tokamaks with the applied 3D fields, the parallel transport is still dominant and the detachment facilitation has not been observed yet. Rather, 3D fields affected detachment adversely under certain conditions, either by preventing detachment onset as seen in DIII-D or by re-attaching the existing detached plasma as shown in NSTX. The possible way for strong 3D effects to induce access to the early detachment in divertor tokamaks appears to be via significant perpendicular loss of parallel momentum by frictional force for the counter-streaming flows between neighboring flow channels in the divertor. In principle, the adjacent lobes in the 3D divertor tokamak may generate the counter-streaming flow channels. However, an EMC3-EIRENE simulation for ITER H-mode plasmas demonstrated that screened RMP leads to significantly reduced counter-flows near the divertor target, therefore the momentum loss effect leading to detachment facilitation is expected to be small. This is consistent with the observation in LHD, which showed screening (amplification) of RMP fields in the attachment (stable detachment) case. In conclusion, work for optimal parameter window for best divertor operation scenario is needed particularly for the 3D divertor tokamak configuration.
Current generation by helicons and LH waves in modern tokamaks and reactors FNSF-AT, ITER and DEMO. Scenarios, modeling and antennae
The Innovative concept and 3D full wave code modeling Off-axis current drive by RF waves in large scale tokamaks, reactors FNSF-AT, ITER and DEMO for steady state operation with high efficiency was proposed [1] to overcome problems well known for LH method [2]. The scheme uses the helicons radiation (fast magnetosonic waves at high (20–40) IC frequency harmonics) at frequencies of 500–1000 MHz, propagating in the outer regions of the plasmas with a rotational transform. It is expected that the current generated by Helicons will help to have regimes with negative magnetic shear and internal transport barrier to ensure stability at high normalized plasma pressure β{sub N} > 3 (the so-called Advanced scenarios) of interest for FNSF and the commercial reactor. Modeling with full wave three-dimensional codes PSTELION and STELEC2 showed flexible control of the current profile in the reactor plasmas of ITER, FNSF-AT and DEMO [2,3], using multiple frequencies, the positions of the antennae and toroidal waves slow down. Also presented are the results of simulations of current generation by helicons in tokamaks DIII-D, T-15MD and JT-60SA [3]. In DEMO and Power Plant antenna is strongly simplified, being some analoge of mirrors based ECRF launcher, as will be shown. For spherical tokamaks the Helicons excitation scheme does not provide efficient Off-axis CD profile flexibility due to strong coupling of helicons with O-mode, also through the boundary conditions in low aspect machines, and intrinsic large amount of trapped electrons, as is shown by STELION modeling for the NSTX tokamak. Brief history of Helicons experimental and modeling exploration in straight plasmas, tokamaks and tokamak based fusion Reactors projects is given, including planned joint DIII-D – Kurchatov Institute experiment on helicons CD [1].
Axisymmetric electrostatic magnetohydrodynamic oscillations in tokamaks with general cross-sections and toroidal flow
The frequency spectrum and mode structure of axisymmetric electrostatic oscillations [the zonal flow (ZF), sound waves (SW), geodesic acoustic modes (GAM), and electrostatic mean flows (EMF)] in tokamaks with general cross-sections and toroidal flows are studied analytically using the electrostatic approximation for magnetohydrodynamic modes. These modes constitute the “electrostatic continua.” Starting from the energy principle for a tokamak plasma with toroidal rotation, we showed that these modes are completely stable. The ZF, the SW, and the EMF could all be viewed as special cases of the general GAM. The Euler equations for the general GAM are obtained and are solved analytically for both the low and high range of Mach numbers. The solution consists of the usual countable infinite set of eigen-modes with discrete eigen-frequencies, and two modes with lower frequencies. The countable infinite set is identified with the regular GAM. The lower frequency mode, which is also divergence free as the plasma rotation tends to zero, is identified as the ZF. The other lower (zero) frequency mode is a pure geodesic E×B flow and not divergence free is identified as the EMF. The frequency of the EMF is shown to be exactly 0 independent of plasma cross-section or its flow Mach number. We also show that in general, sound waves with no geodesic components are (almost) completely lost in tokamaks with a general cross-sectional shape. The exception is the special case of strict up-down symmetry. In this case, half of the GAMs would have no geodesic displacements. They are identified as the SW. Present day tokamaks, although not strictly up-down symmetric, usually are only slightly up-down asymmetric. They are expected to share the property with the up-down symmetric tokamak in that half of the GAMs would be more sound wave-like, i.e., have much weaker coupling to the geodesic components than the other half of non-sound-wave-like modes with stronger coupling to the geodesic displacements. Based on the general notion that the geodesic component of the GAM is more effective in tearing up the eddies in the electrostatic turbulence, it is important to preferentially excite the GAMs that are non-sound-wave like to maximize the efficiency on turbulence suppression through external means. Finally, approximate formulae for the frequencies of the EMF, ZF, SW, and the GAM for a large aspect ratio circular tokamak rotating at low Mach numbers are also provided.
Advances in prediction of tokamak experiments with theory-based models
The successful validation of theory-based models of transport, MHD stability, heating and current drive, with tokamak measurements over the last 20 years, has laid the foundation for a new era where these models can be routinely used in a "predict first" approach to design and predict the outcomes of experiments on tokamaks today. The capability to predict the plasma confinement and core profiles with a quantified uncertainty, based on a multi-machine, international, database of experience, will provide confidence that a proposed discharge will remain within the operational limits of the tokamak. Developing this predictive capability for the first generation of burning plasma devices, beginning with ITER, and progressing to tokamak demonstration reactors, is a critical mission of fusion energy research. Major advances have been made implementing this predict first methodology on today's tokamaks. In this work, an overview of several of these recent advances will be presented, providing the integrated modeling foundations of the experimental successes. The first steps to include boundary plasmas, and tokamak control systems, have been made. A commitment to predicting experiments as part of the planning process is needed in order to collect predictive accuracy data and evolve the models and software into a robust whole discharge pulse design simulator.
Design of vacuum vessel for Aditya Upgrade Tokamak
Aditya tokamak is a limiter tokamak device with a torus vacuum vessel having square cross section. The existing Aditya Tokamak is being upgraded into a machine with divertor operation and good plasma control to support the future Indian Fusion program in a big way. The up-gradation of Aditya Tokamak to diverter tokamak needs additional magnetic field coils for diverter operation, close to the vessel. To accommodate the coils, the existing vessel of square cross section will be replaced with a circular cross section. This vacuum vessel will be made up of stainless steel 304L. It is bakeable up to 150° C. The ultimate vacuum is ∼ 10{sup -9} torr. Two numbers of identical semi torus will be assembled to torus shaped vessel. The end flanges of the semi torus are connected with sealing double Viton 'O' rings to form a torus. The major radius of torus vessel is 750 mm and inner radius is 305 mm. Its wall thickness is ∼8 mm. The vacuum vessel will have 110 various sizes of ports/openings including 4 tangential ports. The nonstandard ports like rectangular and triangular ports are blanked with wire seals. In this paper, the design of vacuum vessel for Aditya Upgrade tokamak is presented. (author)