Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MAGNETIC CONFINEMENT”

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 703 records · Page 39

Observation of a reduced-turbulence regime with boron powder injection in a stellarator

Abstract In state-of-the-art stellarators, turbulence is a major cause of the degradation of plasma confinement. To maximize confinement, which eventually determines the amount of nuclear fusion reactions, turbulent transport needs to be reduced. Here we report the observation of a confinement regime in a stellarator plasma that is characterized by increased confinement and reduced turbulent fluctuations. The transition to this regime is driven by the injection of submillimetric boron powder grains into the plasma. With the line-averaged electron density being kept constant, we observe a substantial increase of stored energy and electron and ion temperatures. At the same time, the amplitude of the plasma turbulent fluctuations is halved. While lower frequency fluctuations are damped, higher frequency modes in the range between 100 and 200 kHz are excited. We have observed this regime for different heating schemes, namely with both electron and ion cyclotron resonant radio frequencies and neutral beams, for both directions of the magnetic field and both hydrogen and deuterium plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Diagnosing magnetized liner inertial fusion experiments on Z

Magnetized Liner Inertial Fusion experiments performed at Sandia's Z facility have demonstrated significant thermonuclear fusion neutron yields (∼10{sup 12} DD neutrons) from multi-keV deuterium plasmas inertially confined by slow (∼10 cm/μs), stable, cylindrical implosions. Effective magnetic confinement of charged fusion reactants and products is signaled by high secondary DT neutron yields above 10{sup 10}. Analysis of extensive power, imaging, and spectroscopic x-ray measurements provides a detailed picture of ∼3 keV temperatures, 0.3 g/cm{sup 3} densities, gradients, and mix in the fuel and liner over the 1–2 ns stagnation duration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the non-stiffness of edge transport in L-mode tokamak plasmas

Transport analyses using first-principle turbulence codes and 11/2 -D transport codes usually study radial transport properties between the tokamak plasma magnetic axis and a normalized minor radius around 0.8. In this region, heat transport shows significantly stiff properties resulting in temperature scalelength values (R∕L{sub T}) that are relatively independent of the level of the radial heat flux. We have studied experimentally in the tokamak à configuration variable [F. Hofmann et al., Plasma Phys. Controlled Fusion 36, B277 (1994)] the radial electron transport properties of the edge region, close to the last closed flux surface, namely, between ρ{sub V}=√(V/V{sub edge})=0.8 to 1. It is shown that electron transport is not stiff in this region and high R∕L{sub Te} values (∼20–40) can be attained even for L-mode confinement. We can define a “pedestal” location, already in L-mode regimes, where the transport characteristics change from constant logarithmic gradient, inside ρ{sub V} = 0.8, to constant gradient between 0.8 and 1.0. In particular, we demonstrate, with well resolved T{sub e} and n{sub e} profiles, that the confinement improvement with plasma current I{sub p}, with or without auxiliary heating, is due to this non-stiff edge region. This new result is used to explain the significant confinement improvement observed with negative triangularity, which could not be explained by theory to date. Preliminary local gyrokinetic simulations are now consistent with an edge, less stiff, region that is more sensitive to triangularity than further inside. We also show that increasing the electron cyclotron heating power increases the edge temperature inverse scalelength, in contrast to the value in the main plasma region. The dependence of confinement on density in ohmic plasmas is also studied and brings new insight in the understanding of the transition between linear and saturated confinement regimes, as well as of the density limit and appearance of a 2/1 tearing mode. The results presented in this paper provide an important new perspective with regards to radial transport in tokamak plasmas which goes beyond L-mode plasmas and explains some previous puzzling results. It is proposed that understanding the transport properties in this edge non-stiff region will also help in understanding the improved and high confinement edge properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Excitation of helicons by current antennas

Depending on the angle θ between the wave vector and the magnetic field, helicons are conventionally divided into two branches: proper helicons (H mode), propagating at small θ, and Trivelpiece–Gould waves (TG mode), propagating at large θ. The latter are close to potential waves and have a significant electric component along the external magnetic field. It is believed that it is these waves that provide electron heating in helicon discharges. There is also commonly believed that current antennas, widely used to ignite helicon discharges, excite essentially nonpotential Н modes, which then transform into TG modes due to plasma inhomogeneity. In this work, it is demonstrated that electromagnetic energy can also be efficiently introduced in plasma by means of TG modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multi-field characteristics and eigenmode spatial structure of geodesic acoustic modes in DIII-D L-mode plasmas

The geodesic acoustic mode (GAM), a coherent form of the zonal flow, plays a critical role in turbulence regulation and cross-magnetic-field transport. In the DIII-D tokamak, unique information on multi-field characteristics and radial structure of eigenmode GAMs has been measured. Two simultaneous and distinct, radially overlapping eigenmode GAMs (i.e., constant frequency vs. radius) have been observed in the poloidal E×B flow in L-mode plasmas. As the plasma transitions from an L-mode to an Ohmic regime, one of these eigenmode GAMs becomes a continuum GAM (frequency responds to local parameters), while the second decays below the noise level. The eigenmode GAMs occupy a radial range of ρ = 0.6–0.8 and 0.75–0.95, respectively. In addition, oscillations at the GAM frequency are observed for the first time in multiple plasma parameters, including n{sub e}, T{sub e}, and B{sub θ}. The magnitude of T(tilde sign){sub e}/T{sub e} at the GAM frequency (the magnitude is similar to that of ñ{sub e}/n{sub e}) and measured n{sub e}–T{sub e} cross-phase (∼140° at the GAM frequency) together indicate that the GAM pressure perturbation is not determined solely by ñ{sub e}. The magnetic GAM behavior, a feature only rarely reported, is significantly stronger (×18) on the high-field side of the tokamak, suggesting an anti-ballooning nature. Finally, the GAM is also observed to directly modify intermediate-wavenumber ñ{sub e} levels (kρ{sub s} ∼ 1.1). The simultaneous temperature, density, flow fluctuations, density-temperature cross-phase, and magnetic behavior present a new perspective on the underlying physics of the GAM.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Polywell Revisited

The Polywell fusion concept, originally proposed by Robert W. Bussard in 1985, has been investigated for over four decades as a potential solution for achieving net fusion energy in a compact and economically viable reactor. It combines two distinct approaches: high-beta magnetic cusp confinement of electrons using polyhedral coil configurations and electrostatic ion confinement via a potential well formed by injected electron beams. While the hybrid nature of the Polywell system offers advantages in plasma stability and engineering simplicity, previous efforts have been limited by persistent challenges in achieving sufficient plasma confinement required to generate a net energy gain. Here, in this study, we examine previous works and identify limitations of several Polywell embodiments that have historically impeded progress. We present an updated Polywell physics model incorporating experimental findings and recent first-principles particle-in-cell simulations. This updated model outlines a credible path toward overcoming confinement losses and achieving net energy gain using deuterium-tritium (D-T) fuels. Our findings provide a renewed scientific basis for the continued development of the Polywell fusion concept as a practical and scalable approach to fusion energy.

Electrostatic potential well↗

Spatially resolved electron density and electron energy distribution function in Ar magnetron plasmas used for sputter-deposition of ZnO-based thin films

Langmuir probe and trace rare gases optical emission spectroscopy were used to analyze the spatial structure of the electron density and electron energy distribution function (EEDF) in a cylindrical Ar magnetron plasma reactor used for sputter-deposition of ZnO-based thin films. While a typical Bessel (zero order) diffusion profile was observed along the radial direction for the number density of charged particles at 21 cm from the ZnO target, a significant rise of these populations with respect to the Bessel function was seen in the center of the reactor at 4 cm from the magnetron surface. As for the EEDF, it was found to transform from a more or less Maxwellian far from the target to a two-temperature Maxwellian with a depletion of high-energy electrons where magnetic field confinement effects become important. No significant change in the behavior of the electron density and EEDF across a wide range of pressures (5–100 mTorr) and self-bias voltages (115–300 V) was observed during magnetron sputtering of Zn, ZnO, and In{sub 2}O{sub 3} targets. This indicates that sputtering of Zn, In, and O atoms do not play a very significant role on the electron particle balance and electron heating dynamics, at least over the range of experimental conditions investigated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Integration of full divertor detachment with improved core confinement for tokamak fusion plasmas

Divertor detachment offers a promising solution to the challenge of plasma-wall interactions for steady-state operation of fusion reactors. Here, we demonstrate the excellent compatibility of actively controlled full divertor detachment with a high-performance (β N ~ 3, H 98 ~ 1.5) core plasma, using high-β p (poloidal beta, β p > 2) scenario characterized by a sustained core internal transport barrier (ITB) and a modest edge transport barrier (ETB) in DIII-D tokamak. The high-β p high-confinement scenario facilitates divertor detachment which, in turn, promotes the development of an even stronger ITB at large radius with a weaker ETB. This self-organized synergy between ITB and ETB, leads to a net gain in energy confinement, in contrast to the net confinement loss caused by divertor detachment in standard H-modes. These results show the potential of integrating excellent core plasma performance with an efficient divertor solution, an essential step towards steady-state operation of reactor-grade plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interplay between beam-driven chirping modes and plasma confinement transitions in spherical tokamak ST40

Experiments on the high field spherical tokamak ST40 have led to the recent observation of interplay between beam-driven modes of sweeping frequency (chirping modes) and transitions to the enhanced global confinement regime (H-mode) and back to the low confinement regime (L-mode). The H-modes of plasma confinement are identified from decreased intensity of D α signal and from clear distinctions in the edge gradients of the visible plasma boundary (observed as a sharp plasma edge in camera images). The beam-driven chirping modes, identified as ideal magnetohydrodynamics beta-induced Alfvén acoustic eigenmodes modes, are observed in Mirnov coil signals, interferometry, and soft x-ray diagnostics. A moderate amplitude ‘primer’ chirping mode usually precedes an H–L transition. This is followed by a ‘dominant’ chirping mode with higher amplitude during the L-mode. The L–H transition back to the improved confinement occurs on a longer time scale of tens of ms, consistent with the slowing down time scale of fast beam ions. A dramatic decrease in toroidal plasma rotation is systematically observed associated with chirping modes sweeping down to zero frequency. Resonance maps built for the beam-driven chirping modes with the ASCOT (accelerated simulation of charged particle orbits in torodoial devices) code show that the resonant beam ions have orbits near the trapped-passing boundary. The ASCOT modelling assesses how losses of the resonant fast ions caused by the chirping modes with high enough amplitude modify the torque, potentially affecting the plasma rotation.

magnetic confinement fusion↗

Elevating zero dimensional global scaling predictions to self-consistent theory-based simulations

In this work, we have developed an innovative workflow, Stability, Transport, Equilibrium, and Pedestal (STEP)-zero-dimensional (0D), within the OMFIT integrated modeling framework. Through systematic validation against the International Tokamak Physics Activity global H-mode confinement database, we demonstrated that STEP-0D, on average, predicts the energy confinement time with a mean relative error of less than 19%. Moreover, this workflow showed promising potential in predicting plasmas for proposed fusion reactors such as the affordable, robust, compact (ARC) reactor, the European demonstration power plant (EU-DEMO), and the China fusion engineering test reactor (CFETR) indicating moderate H-factors between 0.9 and 1.2. STEP-0D allows theory-based prediction of tokamak scenarios, beginning with 0D quantities. The workflow initiates with the PRO-create module, generating physically consistent plasma profiles and equilibrium using the same 0D quantities as the IPB98(y,2) confinement scaling. This sets the starting point for the STEP module, which further iterates between theory-based physics models of equilibrium, core transport, and pedestal to yield a self-consistent solution. Given these attributes, STEP-0D not only improves the accuracy of predicting plasma performance but also provides a path toward a novel fusion power plant design workflow. When integrated with engineering and costing models within an optimization, this new approach could eliminate the iterative reconciliation between plasma models of varying fidelity. This potential for a more efficient design process underpins STEP-0D's significant contribution to future fusion power plant development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetic-Field Effects in Unstable High-Energy-Density Plasmas

The influence of magnetic fields in high-energy-density (HED) systems is a widely unexplored field of research that can play an important role in magnetized inertial confinement fusion (ICF) and in astrophysical systems. Hydrodynamic instabilities are presently believed to be the dominant degradation mechanism in ICF implosions. Therefore, understanding the influence of an imposed B-field on RT growth is important for future mitigation strategies as well as fundamental physics of ICF implosions that utilize B-fields for enhanced performance. The Crab Nebula, one of the most observed objects in our universe, contains elongated spikes of material created by hydrodynamic instabilities, however these spikes do not break-up and become turbulent as they evolve late in time. One hypothesis for this behavior is that B-fields around these spikes prevent vortex generation typical in these types of systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Turbulence Nonlinearities Shed Light on Geometric Asymmetry in Tokamak Confinement Transitions

We report a comprehensive study of fully frequency-resolved nonlinear kinetic energy transfer has been performed for the first time in a diverted tokamak, providing new insight into the parametric dependences of edge turbulence transitions. Measurements using gas puff imaging in the turbulent L -mode state illuminate the source of the long known but as yet unexplained “favorable-unfavorable” geometric asymmetry of the power threshold for transition to the turbulence-suppressed H mode. Results from the recently discovered I mode point to a competition between zonal flow (ZF) and geodesic-acoustic modes (GAM) for turbulent energy, while showing new evidence that the I-to-H transition is still dominated by ZFs. The availability of nonlinear drive for the GAM against net heat flux through the edge corresponds very well to empirical scalings found experimentally for accessing the I mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Curvature of Radial Electric Field Aggravates Edge Magnetohydrodynamics Mode in Toroidally Confined Plasmas

We show that the radial electric field (E r ) plays a dual role in edge magnetohydrodynamics (MHD) activity. While E r shear (first spatial derivative of E r ) dephases radial velocity and displacement, and so is stabilizing, a new finding here is that E r curvature (second spatial derivative of E r ) tends to synchronize the radial velocity and displacement, and so destabilizes MHD. As a highlighted result, we analytically demonstrate that E r curvature can destabilize an otherwise stable kink mode, and so form a joint vortexkink mode. The synergetic effects of E r shear and E r curvature in edge MHD extend the familiar E × B shearing paradigm. This theory thus explains the experimental findings that a deeper E × B well may aggravate edge MHD, and so trigger the formation of the edge harmonic oscillation. Here, a simple criterion linking E r structure and the edge MHD activity is derived.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Project Title: Demonstration High Temperature Superconducting NonPlanar Stellarator Magnet with Advanced Manufactured Assemblies

This is the final report for the project “Demonstration High Temperature Superconducting Non- Planar Stellarator Magnet with Advanced Manufactured Assemblies”, funded by DOE, and performed by Type One Energy from September, 2020, to March 2024 involving the Fusion Technology Institute at the University of Wisconsin–Madison, the Plasma Science and Fusion Center (PSFC) at the Massachusetts Institute of Technology (MIT) and Commonwealth Fusion Systems (CFS) to design and fabricate the first non-planar HTS (REBCO) coil for a high-field stellarator based on the SPARC tokamak’s VIPER cable concept. Stellarators at high fields make high-temperature superconducting magnets necessary for a compact fusion device. But the asymmetric and non-planar nature of its components, especially the magnets make it difficult for scalable producibility. To address these challenges, two promising technologies have emerged: advanced manufacturing (AM) for the supporting plates for forming the magnets, and high-temperature superconducting (HTS) cables inside the plates. AM has advanced enough to produce stellarator components with the necessary geometric complexity, size, and the precision, leading to potentially significant reductions in production time, cost, and waste. The cost of HTS tape has decreased dramatically, and progress in HTS planar magnet development has reached a point where it can be proposed for application to complex 3D non-planar magnets. The main objective of this project is to develop, demonstrate and pre-commercialize a novel, non-planar HTS coil shape that remains superconducting to achieve production scalable reductions in time and cost and performance. The proposed technology is based on the novel concept of a precision sub-scale HTS nonplanar coil assembly. This project focuses on the design, fabrication, material optimization of cable design, and validation and demonstration of the high current carrying capability of superconducting magnets and their support in a complex 3D shape needed for application to stellarator magnetic plasma confinement. The specific objectives of this research program include: (1) The successful application of metal AM to build a precision sub-scale HTS nonplanar coil, (2) An HTS cable and cross-section design that can conform to the required nonplanar coil shape (bend radii as tight as 10-cm) and remains superconducting at an engineering current density of 1.35 kA/cm 2 at 77 K and 1 tesla at the conductor (5 kA in the cable). To achieve the above challenging goals, we have formed a multidisciplinary research team consisting of members from Type One Energy and UW-Madison, MIT PSFC and CFS with complementary skills and strong facilities. The team worked collaboratively on fundamental and applied research on the following three major technical areas: (1) Design, fabrication, and optimization of non-planar HTS Cable The ultimate goal of the project is to determine if commercial REBCO tapes and additive manufacturing can be used to fabricate high field (≥ 10T) non-planar coils with tight bending radii (≃ 100mm) and with a degradation of the critical current (Ic) smaller than 20% with respect to the expected performance. We started with shorter length cable to evaluate the scalability of the production process and eventually reached multiple turns for higher magnetic fields. Our findings suggest that a stellarator coil system of a relevant size, characterized by its asymmetric and non-planar components, can be fabricated using a formed cable in plate method. This system can be simulated using a large-scale modeling approach. The use of hybrid modeling 3 techniques will be pivotal in reducing the complexity of the model and in assessing expected performance in designs. (2) Modeling and simulation of the non-planar HTS Cable Multiphysics simulations are performed using the commercial software and are carried out in self-field conditions, involving 2D and 3D models and twisted around one slot of twist-pitched VIPER cable. Multiphysics simulations are mainly focused on the coil for the critical current evaluation, the magnetic field map, self-Lorentz forces and mechanical, and magnetothermal behavior and the quench dynamics. The detailed model and prediction of the superconducting performance of a stellarator-relevant demonstration cable from numerical simulations supports the results from the actual testing backing the results. A detailed description and results are provided in the later sections. (3) Design, fabrication, and optimization of support for the non-planar HTS Cable The team developed an additive manufactured (AM) coil positioning plate that formed into the required non-planar geometry (with bend radii as tight as 10-cm) and to acceptable tolerances required for a stellarator magnet: (+0.25-mm from ideal on dimensions of coil positioning plates and up to +1-mm from ideal for position of wound coil). The plate materials is also included in this selection process from fabrication and 3D printing perspective and commensurate with eventual application to a fusion reactor. From the cost effectiveness point of view, the HTS coil and plate has the potential to cost less than that made in conventional methods with less waste (<75% waste) reducing time (<50%) and cost (<50%), especially as the AM field matures. The application of advanced manufacturing in the construction of the support plates will also lead to cost reduction, as the cables can be easily replaced, thereby making the assembly modular. With such high primary cost and time savings, high current densities and magnetic field, the funded R&D work has validated the designs, proven the feasibility, and characterized the performance of the HTS coil and plate assembly, paving the way for a relevant-size stellarator coil system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tokamak ITG-KBM transition benchmarking with the mixed variables/pullback transformation electromagnetic gyrokinetic scheme

Electromagnetic gyrokinetic simulation of high temperature plasma is required to predict confinement in magnetic fusion devices and has posed challenges for existing codes. In this paper, we demonstrate successful global gyrokinetic simulation of the ion temperature gradient-driven mode-kinetic ballooning mode transition in a toroidal fusion plasma test case using the mixed variables/pullback transformation (MV/PT) scheme with the particle-in-cell codes XGC and ORB5, and compare to results from a conventional continuum code from the literature. Furthermore, the MV/PT scheme combines explicit time integration with mitigation of the well-known electromagnetic gyrokinetic “cancelation problem.” We calculate eigenmodes in the electrostatic and parallel vector potentials, and find good agreement in growth rate, real frequency, and the normalized plasma pressure of mode transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fast Ion Transport in the Three-Dimensional Reversed-Field Pinch

Here, we report on the first comprehensive experimental and numerical study of fast ion transport in the helical reversed-field pinch (RFP). Classical orbit effects dominate the macroscopic confinement properties. The strongest effect arises from growth in the dominant fast ion guiding-center island, but substantial influence from remnant subdominant tearing modes also plays a critical role. At the formation of the helical RFP, neutron flux measurements indicate a drastic loss of fast ions at sufficient subdominant mode amplitudes. Simulations corroborate these measurements and suggest that subdominant tearing modes strongly limit fast ion behavior. Previous work details a sharp thermal transport barrier and suggests the helical RFP as an Ohmically heated fusion reactor candidate; the enhanced transport of fast ions reported here identifies a key challenge for this scheme, but a workable scenario is conceivable with low subdominant tearing mode amplitudes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Flexible, integrated modeling of tokamak stability, transport, equilibrium, and pedestal physics

The STEP (Stability, Transport, Equilibrium, and Pedestal) integrated-modeling tool has been developed in OMFIT to predict stable, tokamak equilibria self-consistently with core-transport and pedestal calculations. STEP couples theory-based codes to integrate a variety of physics, including magnetohydrodynamic stability, transport, equilibrium, pedestal formation, and current-drive, heating, and fueling. The input/output of each code is interfaced with a centralized ITER-Integrated Modelling & Analysis Suite data structure, allowing codes to be run in any order and enabling open-loop, feedback, and optimization workflows. This paradigm simplifies the integration of new codes, making STEP highly extensible. STEP has been verified against a published benchmark of six different integrated models. Core-pedestal calculations with STEP have been successfully validated against individual DIII-D H-mode discharges and across more than 500 discharges of the H98,y2 database, with a mean error in confinement time from experiment less than 19%. STEP has also reproduced results in less conventional DIII-D scenarios, including negative-central-shear and negative-triangularity plasmas. Predictive STEP modeling has been used to assess performance in several tokamak reactors. Simulations of a high-field, large-aspect-ratio reactor show significantly lower fusion power than predicted by a zero-dimensional study, demonstrating the limitations of scaling-law extrapolations. STEP predictions have found promising scenarios for an EXhaust and Confinement Integration Tokamak Experiment, including a high-pressure, 80%-bootstrap-fraction plasma. ITER modeling with STEP has shown that pellet fueling enhances fusion gain in both the baseline and advanced-inductive scenarios. Finally, STEP predictions for the SPARC baseline scenario are in good agreement with published results from the physics basis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic neoclassical transport in the H-mode pedestal

Multi-species kinetic neoclassical transport through the QH-mode pedestal and scrape-off layer on DIII-D is calculated using XGC0, a 5D full-f particle-in-cell drift-kinetic solver with self-consistent neutral recycling and sheath potentials. Quantitative agreement between the flux-driven simulation and the experimental electron density, impurity density, and orthogonal measurements of impurity temperature and flow profiles is achieved by adding random-walk particle diffusion to the guiding-center drift motion. The radial electric field (E{sub r}) that maintains ambipolar transport across flux surfaces and to the wall is computed self-consistently on closed and open magnetic field lines and is in excellent agreement with experiment. The E{sub r} inside the separatrix is the unique solution that balances the outward flux of thermal tail deuterium ions against the outward neoclassical electron flux and inward pinch of impurity and colder deuterium ions. Particle transport in the pedestal is primarily due to anomalous transport, while the ion heat and momentum transport are primarily due to the neoclassical transport. The full-f treatment quantifies the non-Maxwellian energy distributions that describe a number of experimental observations in low-collisionallity pedestals on DIII-D, including intrinsic co-I{sub p} parallel flows in the pedestal, ion temperature anisotropy, and large impurity temperatures in the scrape-off layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗