Search NASA⌕ Search

SEARCH · Search NASA

Results for “tokamak control”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Corroborating VNA and thermal measurements of transmission loss on the DIII-D ECH waveguide system

Electron cyclotron heating (ECH) and current drive (ECCD) will play a large role in tokamak-based fusion reactors. At the DIII-D tokamak, 110 GHz microwaves injected into the plasma can provide core heating and current drive as well as impurity control, neoclassical tearing mode mitigation, and breakdown assistance. Understanding the physics of these processes relies on accurate estimates of injected ECH power. DIII-D’s ECH system consists of six MW-class Microwave Power Products (MPP) gyrotron microwave sources. Operating the gyrotrons far from the tokamak removes them from magnetic field interference, so 31.75 mm inner-diameter corrugated waveguides transmit the microwave power the 80 m from the gyrotrons to steerable launchers in the tokamak chamber. Estimates of injected power rely on knowing the generated power at the source and then subtracting transmission loss. Conventional transmission loss measurements based on calorimetric dummy loads are onerous and only possible during extended maintenance periods. This work examines two tools that provide more flexibility for the transmission loss measurements. Furthermore, a resistive temperature detector (RTD) array installed along a waveguide measures heat lost to the transmission line, and low power time domain reflectometry (TDR) measurements with a vector network analyzer (VNA) allows loss measurements without burdensome hardware modifications.

ECH↗
DOE OSTI2022

Non-disruptive error field identification based on magnetic island healing

Here a technique to identify intrinsic error fields (EFs) in tokamaks with minimized risk of disruption is demonstrated on the DIII-D tokamak. The method extends the conventional driven magnetic island ‘compass scan’ approach by modifying asynchronous control waveforms to enable prompt healing of the island instability. Healing of the island is achieved by reducing the imposed non-axisymmetric coil current and raising the density (here via gas fueling). The method is also shown to support multiple island threshold measurements per pulse, thus reducing the number of dedicated pulses necessary to conduct an EF identification. Non-linear modeling with the TM1 code reproduces the experimental results and approximately recovers the critical density required for island healing. Island healing is explained in the non-linear modeling by an increase in the viscous coupling between the static island and the nearby flowing plasma, thus healing the island as it accelerates into the plasma frame. Due to both simplicity and risk minimization, this technique is suitable for plasma-based EF identification in the early commissioning stages of future disruption-averse tokamaks such as ITER and SPARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

DIII-D research advancing the physics basis for optimizing the tokamak approach to fusion energy

DIII-D physics research addresses critical challenges for the operation of ITER and the next generation of fusion energy devices. This is done through a focus on innovations to provide solutions for high performance long pulse operation, coupled with fundamental plasma physics understanding and model validation, to drive scenario development by integrating high performance core and boundary plasmas. Substantial increases in off-axis current drive efficiency from an innovative top launch system for EC power, and in pressure broadening for Alfven eigenmode control from a co-/counter-I p steerable off-axis neutral beam, all improve the prospects for optimization of future long pulse/steady state high performance tokamak operation. Fundamental studies into the modes that drive the evolution of the pedestal pressure profile and electron vs ion heat flux validate predictive models of pedestal recovery after ELMs. Understanding the physics mechanisms of ELM control and density pumpout by 3D magnetic perturbation fields leads to confident predictions for ITER and future devices. Validated modeling of high-Z shattered pellet injection for disruption mitigation, runaway electron dissipation, and techniques for disruption prediction and avoidance including machine learning, give confidence in handling disruptivity for future devices. For the non-nuclear phase of ITER, two actuators are identified to lower the L–H threshold power in hydrogen plasmas. With this physics understanding and suite of capabilities, a high poloidal beta optimized-core scenario with an internal transport barrier that projects nearly to Q = 10 in ITER at ~8 MA was coupled to a detached divertor, and a near super H-mode optimized-pedestal scenario with co-I p beam injection was coupled to a radiative divertor. The hybrid core scenario was achieved directly, without the need for anomalous current diffusion, using off-axis current drive actuators. Also, a controller to assess proximity to stability limits and regulate β N in the ITER baseline scenario, based on plasma response to probing 3D fields, was demonstrated. Finally, innovative tokamak operation using a negative triangularity shape showed many attractive features for future pilot plant operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Active control of resistive wall mode via modification of external tearing index

Modification of the external tearing index, Δ ext ′ , by magnetic feedback is analytically investigated for the purpose of controlling the resistive plasma resistive wall mode (RP-RWM). The matching method is pursued by deriving expressions for the close-loop Δ ext ′ and by linking it to the counterpart from the inner layer. Various feedback coil configurations are found to generally reduce Δ ext ′ and stabilize the RWM, with either proportional or derivative control. Feedback modification of Δ ext ′ is found to be generally independent of the inner layer resistive interchange index D R , confirming that feedback action primarily modifies the solution in the outer ideal region for the RP-RWM. Exception occurs when either the inner layer favorable curvature effect becomes sufficiently large or the feedback action is sufficiently strong to introduce a rotating RP-RWM in the static plasma, leading to complex-valued close-loop Δ ext ′ . The perturbed magnetic energy dissipation in the outer region, associated with the eddy current in the resistive wall, is identified as the key physics reason for feedback induced complex Δ ext ′ . Similar results are also obtained for active control of the external kink instability, whose open-loop growth rate is significantly reduced by inclusion of the plasma resistivity. Within the single poloidal harmonic approximation, which is most suitable for the matching approach, external active coils combined with poloidal sensors are often found to be more efficient for feedback stabilization of the mode at large proportional gain values. This counter-intuitive result is explained as the lack of (non-resonant) poloidal harmonics for proper description of the feedback coil geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Editorial for the 28th Workshop on Magnetohydrodynamic (MHD) Stability Control

This special issue contains four papers which were presented at the 28th Workshop on Magnetohydrodymanic (MHD) Stability Control. The workshop was held at Auburn University and on Zoom from 3–5 October 2024. Furthermore, this annual workshop covers the active and passive control of MHD instabilities such as kink, tearing and edge-localized modes in tokamaks, reversed field pinches, spherical tori or other magnetic confinement devices.

Yang, J. [Princeton Plasma Physics Laboratory (PPP↗
DOE OSTI2026

Power and particle exhaust in the ST-E1 fusion power plant

Power exhaust challenges and potential solutions for a 5 m major radius, low-aspect ratio burning tokamak have been explored. 1D edge plasma models have been used to screen for access to detachment using short and long outer divertor legs in double and single null configurations, using Ar as the primary impurity and assuming tungsten plasma-facing components (PFCs). These show that detachment access can be accessed for all but the most conservative assumptions on scrape-off layer (SOL) width and power, but that trade-offs will be required between magnet engineering and the size of the acceptable window of as-yet uncertain plasma parameters. SOLPS-ITER was used to further model selected plasma scenarios, confirming that Ar seeding can be used to achieve dissipative divertor scenarios with peak deposited heat fluxes below 15 MWm -2 . Initial scoping of first wall loads and positioning of limiters has been carried out, showing the feasibility of protecting the breeding blanket wall during steady state without impeding tritium breeding. Initial PFC technology selection is also presented, identifying this as a critical area where further work is needed to find an attractive solution for helium-cooled PFCs that can handle high heat fluxes without excessive power requirements. Key questions and trade-offs for concept development have been identified, including: how to achieve high radiation for reduction of SOL power without core performance degradation; whether power exhaust can be well-controlled in a double null plasma; mechanical design and materials challenges of high-heat flux PFCs; and control of material erosion, redeposition and tritium retention.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Advances in physics and applications of 3D magnetic perturbations on the J-TEXT tokamak

In the last two years, three major technical improvements have been made on J-TEXT in supporting of the expanded operation regions and diagnostic capabilities. (1) The successful commission of the 105 GHz/500 kW/1 s electron cyclotron resonance heating (ECRH) system increasing the core electron temperature from 0.9 keV up to around 1.5 keV. (2) The poloidal divertor configuration with an X-point in the high-field side has been achieved. In particular, the 400 kW electron cyclotron wave has also been successfully injected into the diverted plasma. (3) A 256-channel electron cyclotron emission imaging diagnostic system and two sets of four-channel Doppler backscattering diagnostics have been successfully developed on J-TEXT, allowing detailed measurement of the electron temperature and density fluctuations for turbulence and MHD research. The locked mode (LM), especially the 2/1 LM, is one of the biggest threats to the plasma operation. Both the thresholds of 2/1 and 3/1 LM are observed to vary non-monotonically on electron density. The electrode biasing was applied successfully to unlock the LM from either a rotating or static resonant magnetic perturbation (RMP) field. In the presence of 2/1 LM, three kinds of standing wave (SW) structures have been observed to share a similar connection to the island structure, i.e. the nodes of the SWs locate around the O- or X-points of the 2/1 island. The control and mitigation of disruption is essential to the safe operation of ITER, and it has been systematically studied by applying a RMP field, massive gas injection (MGI) and shattered pellet injection on J-TEXT. When the RMP-induced 2/1 LM is larger than a critical width, the MGI shutdown process can be significantly influenced. If the phase difference between the O-point of LM and the MGI valve is +90° (or –90°), the penetration depth and the assimilation of impurities can be enhanced (or suppressed) during the pre-thermal quench (TQ) phase and result in a faster (or slower) TQ. A secondary MGI can also suppress the runaway electron (RE) generation, if the additional high-Z impurity gas arrives at the plasma edge before TQ. When the secondary MGI has been applied after the formation of the RE current plateau, the RE current can be dissipated, and the dissipation rate increases with the injected impurity quantity but saturates with a maximum of 28 MA s –1 . The non-local transport is experimentally observed in the ion transport channel. The electron thermal diffusivity significantly increases with the ECRH power. Theoretical work shows that significant intrinsic current can be driven by electromagnetic turbulence, and the robust formation mechanism of the E × B staircase is identified from the Hasegawa–Wakatani system.

divertor↗
DOE OSTI2022

Disruption Prediction and Avoidance in High Beta Long Pulse KSTAR Plasmas

This research is a concentrated and coordinated experimental and physics analysis program by a small U.S. team, focused on the topical area of “Long Pulse Control” as stated in the proposal solicitation. The overall goal of the research is to produce world-class analysis regarding forecasting and avoidance of disruptions in a long-pulse, superconducting tokamak at high stability parameters (e.g. normalized beta) with direct, quantifiable figures of merit (i.e. plasma disruption rate - “disruptivity”) and its reduction based on improved stability physics and control. The research targets the backbone of the problem – the basic, main issues – to produce this quantitative analysis by: (1) Creating and applying a new physics-based paradigm to understand how to eliminate disruptions: disruption event characterization, and event chain analysis – a direct manifestation of the approach suggested in the DOE Transients Report, (2) Ensuring that the components required to produce high-quality measurements and analysis needed for disruption event characterization for the tokamak in question (here, KSTAR) exist (e.g. pressure profile peaking, high-accuracy pitch angle data, compensated low frequency toroidal mode number n = 1 mode measurement) to create a primarily physics-based disruption forecasting capability, and (3) Addressing disruption avoidance for the first time in KSTAR by attacking primary disruption ‘event chains’, utilizing the most important associated plasma measurements (dynamic n = 1 field). The focus of this work will be avoidance of macroscopic instabilities, especially global modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Role of turbulent separatrix tangle in the improvement of the integrated pedestal and heat exhaust issue for stationary-operation tokamak fusion reactors

Abstract The magnetic separatrix surface is designed to provide the final and critical confinement to the hot stationary-operation core plasma in modern tokamak reactors in the absence of an external magnetic perturbation (MP) or transient magneto-hydrodynamic perturbation, while diverting the exhaust heat to divertor plates. All the stationary operational boundary plasma studies and reactor designs have been performed under this assumption. However, there has been a long-standing suspicion that a stationary-operation tokamak plasma even without external MPs or edge localized modes (ELMs) activities may not have a stable closed separatrix surface, especially near the magnetic X-point. Here, the first gyrokinetic numerical observation is reported that the divertor separatrix surface, due to homoclinic tangles caused by intrinsic electromagnetic turbulence, is not a stable closed surface in a stationary operation phase even without MPs or ELMs. Unlike the MP- or ELM-driven homoclinic tangles that could cause deleterious effects to core confinement or divertor plates, it is found that the micro-turbulence driven homoclinic tangles could connect the divertor plasma to the pedestal plasma in a constructive way by broadening the divertor heat-exhaust footprint and weakening the pedestal slope to the ELM-safe direction. Micro-turbulent homoclinic tangles can open a new research direction in understanding and controlling these two most troublesome and non-locally connected edge-plasma issues in a tokamak fusion reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Sensitivities of time-dependent temperature profile predictions for NSTX with the multi-mode model

The Multi-Mode Model (MMM) for turbulent transport was applied to a large set of well-analyzed discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate its sensitivities to a wide range of plasma conditions. MMM calculations were performed for hundreds of milliseconds in each discharge by performing time-dependent predictive simulations with the 1.5D tokamak integrated modeling code TRANSP. A closely related study (Lestz et al 2025 Plasma Phys. Control. Fusion 67 105029) concluded that MMM predicted electron and ion temperature profiles that were in reasonable agreement with NSTX observations, generally outperforming a different reduced transport model, TGLF. This finding motivates the more thorough investigation of the characteristics of the MMM predictions conducted in this work. The simulations with MMM have electron energy transport dominated by electron temperature gradient modes in the examined discharges with relatively low plasma β (ratio of kinetic plasma pressure to magnetic field pressure) and high collisionality, transitioning to a mixture of different modes for higher β and lower collisionality. The thermal ion diffusivity predicted by MMM is much smaller than the neoclassical contribution, in line with previous experimental analysis of NSTX. Nonetheless, the electron and ion temperature profiles are coupled via collisional energy exchange and thus sensitive to which transport channels are predicted. The time-dependent simulations with MMM are robust to the simulation start time, converging to remarkably similar temperature profiles later during the discharge. MMM typically overpredicts confinement relative to NSTX observations, leading to the prediction of overly steep temperature profiles. Plasmas with spatially broader temperature profiles, higher plasma β, and longer energy confinement times tend to be predicted by MMM with better agreement with the experiment. As a result, these findings provide useful context for understanding the regime-dependent tendencies of MMM in anticipation of self-consistent, time-dependent predictive simulations of NSTX-U discharges with these same modeling tools.

MMM↗
DOE OSTI2022

Modeling plasma response to non-axisymmetric magnetic field perturbations in tokamak boundaries

The primary goal of this project is to provide physics insight into models that elucidate the mechanisms responsible for controlling the performance of H-mode plasmas and edge localized mode (ELM) stability when subjected to intrinsic and applied non-axisymmetric magnetic perturbations. The overall goal is to facilitate the development of physics-based criteria for the design of magnetic coils to be used for ELM suppression and H-mode pedestal control in future burning plasma devices such as ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Impact of Variable Perpendicular Transport Coefficients in WEST Simulations Using SolEdge-HDG

Plasma–wall interaction is one of the key research topics on the way to controlled fusion. To study the best operational designs with reduced heat and particle fluxes onto tokamak plasma facing components (PFCs) comprehensive plasma simulations are required. A recent implementation of a hybridized discontinuous Galerkin scheme into a new version of SolEdge code has the advantage of using magnetic equilibrium-free mesh. This allows us to conduct pioneering 2-D transport simulations of a full discharge in the WEST tokamak. In this work, we implemented plasma transport coefficients as functions of coordinate in the poloidal plane and neutral diffusion as a function of neutral mean free path. Moreover, the perpendicular convection flux terms were added to the code. Using the new features, a few test cases were investigated. Finally, the influence of nonconstant transport coefficients on the simulated particle and heat fluxes onto the WEST tokamak PFCs are demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

US-UK fusion energy collaborations in the digital space

The US and UK share the vision for fusion as a vital part of the clean energy future. This vision is reflected in the respective national plans in the form of the Bold Decadal Vision for Commercial Fusion (BDV) in the US and the Spherical Tokamak for Energy Production (STEP) program in the UK. Digital tools such as simulation and control frameworks, design tools, AI/ML, high performance computing (HPC), and virtual reality (VR) will play an important role in developing, diagnosing, operating, and further improving burning-plasma-class fusion power plants. Therefore, a collaborative approach, involving both the public and private sectors, to developing these digital tools can accelerate the path to fusion energy commercialization. In this report we discuss previous and ongoing collaborations and opportunities to expand these collaborations into new areas. We conclude this report with near term actions and a vision of the collaboration to the Joint Coordinating Committee.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Impurity screening behavior of the high-field side scrape-off layer in near-double-null configurations: prospect for mitigating plasma-material interactions on RF actuators and first-wall components

The impurity screening response of the high-field side (HFS) scrape-off layer (SOL) to localized nitrogen injection is investigated on Alcator C-Mod for magnetic equilibria spanning lower-single null, double-null and upper single null configurations under otherwise identical plasma conditions. L-mode, EDA H-mode and I-mode discharges are investigated. HFS impurity screening is found to depend on magnetic flux balance and the direction of B×?B relative to the most active divertor. Impurity ‘plume’ emission patterns indicate that both parallel and perpendicular (E×B) flows in the SOL contribute to the ‘flushing’ of impurities towards the active divertor, thereby affecting the overall impurity screening behavior. Despite the fact that the HFS scrape-off layer is extremely narrow in near-double-null configurations, this SOL is able to screen locally injected nitrogen at least as effectively as the low-field side (LFS) SOL – up to a factor of 10 more effective, depending on specific plasma conditions and whether the magnetic geometry produces parallel flows that work with or against E×B flows. For situations in which the E×B drift of the impurity ions opposes parallel flow toward the primary divertor, HFS impurity screening is found to be least effective. When E×B drifts assist parallel flow toward the primary divertor, HFS impurity screening is found to be very effective. These data support the idea of placing RF actuators and close-fitting wall components on the high-field side of the tokamak. With this configuration, near-double-null magnetic topologies may be used for active control of plasma parameters at the antenna/plasma interface for optimal RF coupling, mitigate the generation of local impurities from plasma-material interactions and, taking advantage of favorable plasma flows and good screening properties of the HFS SOL, further minimize the impact of wall-born impurity sources on the plasma core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Toward active disruption avoidance via real-time estimation of the safe operating region and disruption proximity in tokamaks

This paper describes a real-time capable algorithm for identifying the safe operating region around a tokamak operating point. The region is defined by a convex set of linear constraints, from which the distance of a point from a disruptive boundary can be calculated. The disruptivity of points is calculated from an empirical machine learning predictor that generates the likelihood of disruption. While the likelihood generated by such empirical models can be compared to a threshold to trigger a disruption mitigation system, the safe operating region calculation enables active optimization of the operating point to maintain a safe margin from disruptive boundaries. The proposed algorithm is tested using a random forest disruption predictor fit on data from DIII-D. The safe operating region identification algorithm is applied to historical data from DIII-D showing the evolution of disruptive boundaries and the potential impact of optimization of the operating point. Real-time relevant execution times are made possible by parallelizing many of the calculation steps and implementing the algorithm on a graphics processing unit. Lastly, a real-time capable algorithm for optimizing the target operating point within the identified constraints is also proposed and simulated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Drift kinetic theory of the NTM magnetic islands in a finite beta general geometry tokamak plasma

In (Imada et al 2019 Nucl. Fusion 59 046016 and references therein) a new 4D drift kinetic nonlinear theory, valid in the limit of a low beta, small inverse aspect ratio, circular cross section, toroidal geometry, to describe the plasma response to the neoclassical tearing mode (NTM) magnetic perturbation is derived. In (Dudkovskaia et al 2021 Plasma Phys. Control. Fusion 63 054001) this theory is reduced in a low collisionality limit, which allows a dimensionality reduction to a 3D problem to efficiently resolve the collisional dissipation layer in the vicinity of the trapped-passing boundary. (Dudkovskaia et al 2021 Plasma Phys. Control. Fusion 63 054001) adopts an improved model for the magnetic drift frequency, which reduces the threshold magnetic island half-width from $8.73\rho_{b i}$, where $\rho_{b i}$ is the trapped ion banana orbit width, to $1.46 \rho_{b i}$, making it in closer agreement with experimental observations for the large aspect ratio tokamak equilibrium. In the present paper, the theory is extended to a high beta, arbitrary tokamak geometry to capture the plasma shaping effects on the NTM threshold physics with the focus on the non-zero triangularity discharges that are known to have a strong impact on the plasma MHD stability. First, it is found that the higher triangularity plasma is more prone to NTMs which is in agreement with the 2/1 tearing mode onset relative frequency measurements in DIII-D. Second, the NTM threshold dependence on the tokamak inverse aspect ratio obtained in (Dudkovskaia et al 2021 Plasma Phys. Control. Fusion 63 054001) is refined and extended to a finite aspect ratio limit. Third, the NTM threshold dependence on poloidal beta is obtained and benchmarked against the EAST threshold island width measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Development of advanced vacuum technologies for extending plasma pulse duration on EAST

Advanced vacuum technologies, including pumping, fueling and wall conditioning, have been successfully developed or upgraded to efficiently control the fuel and impurity particles to extend the plasma pulse duration in the experimental advanced superconducting tokamak (EAST). To improve the particle exhaust rate cryopumps with a 60% increase in pumping speed and ∼2 times increase in saturation capacity have been developed, and molecular pumps with a ∼30% increase in pumping speed have been upgraded. In order to monitor the molecular pump status while avoiding bearing faults and overload accidents, a fault detection system has been built which can offer an early warning to avoid more losses within the fusion device. A series of fueling technologies have been developed including gas injection system, supersonic molecular beam injector, pellet injector (PI), massive gas injector and shattered pellet injector, installed at the midplane and divertor positions at different ports to improve fueling uniformity and efficiency. Meanwhile, routine wall conditioning such as electric and hot N 2 baking, ion cyclotron wall conditioning and glow discharge cleaning have been successfully developed to remove the impurity particles from the inner component and materials. The low Z material wall coating and real-time powder injection during plasma discharge are also designed and applied to further improve particle control capability. Finally, by using these advanced vacuum related technologies, good vacuum (<2 × 10 −6 Pa) and wall conditions are realized, and the fuel and impurity particles can be effectively and stably controlled, which promotes the achievement of the record plasma of ∼1056 s pulse duration with the line-averaged electron density of 1.8 × 10 19 m −3 on EAST. They provide a very important reference for vacuum system design and operation for future fusion devices.

EAST↗
DOE OSTI2021

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↗