Search NASA⌕ Search

SEARCH · Search NASA

Results for “divertors”

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 109 records · Page 6

Resilient stellarator divertor characteristics in the Helically Symmetric eXperiment

Resilient divertor features connected to open chaotic edge structures in the Helically Symmetric eXperiment are investigated. For the first time, an expanded vessel wall was considered that would give space for implementation of a physical divertor target structure. The analysis was done for four different magnetic configurations with very different chaotic plasma edges. A resilient plasma wall interaction pattern was identified across all configurations. This manifests as qualitatively very similar footprint behavior across the different plasma equilibria. Overall, the resilient field lines of interest with high connection length L C lie within a helical band along the wall for all configurations. This resiliency can be used to identify the best location of a divertor. The details of the magnetic footprint’s resilient helical band is subject to specific field line structures which are linked to the penetration depth of field lines into the plasma and directly influence the heat and particle flux patterns. The differences arising from these details are characterized by introducing a new metric, the minimum radial connection min $(\delta_N)$ of a field line from the last closed flux surface. The relationship, namely the deviation from a scaling law, between min $(\delta_N)$ and L C of the field lines in the plasma edge field line behavior suggests that the field lines are associated with structures such as resonant islands, cantori, and turnstiles. This helps determine the relevant magnetic flux channels based on the radial location of these chaotic edge structures and the divertor target footprint. These details will need to be taken into account for resilient divertor design.

divertor↗

Optimizing detachment control using the magnetic configuration of divertors

Abstract As tokamak research moves to reactor conditions, the control of a stable, optimally-detached divertor plasma has become increasingly relevant. Simple predictions of such detachment control have been performed previously using the detachment location sensitivity (DLS) model (Lipschultz et al 2016 Nucl. Fusion 56 056007). In this study the DLS model is extended and combined with SOLPS-ITER simulations of isolated divertor grids to study the effects of alternate divertor magnetic field properties on detachment control. The DLS model predicts that divertors can achieve easier access to detachment through a long connection length, a high total flux expansion, and a high average magnetic field in the divertor compared with that at the x -point. SOLPS-ITER simulations suggest an even stronger impact of total flux expansion and connection length on detachment access than the DLS model predicts. In terms of detachment evolution, both simulation and modelling show a high gradient in the total magnetic field and low B p o l B are able to more easily keep a detachment front at a desired poloidal location. Regions with no stable detachment front locations can arise for high magnetic field gradients directed towards the target. Significant differences have been found between impurity scan simulations and the DLS model. These differences may be attributed to sources and sinks of power and electron pressure, and should be explored further in future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of edge harmonic oscillations on the divertor heat flux in NSTX

Previously, ELM-free and inter-ELM divertor peak heat flux reduction induced by an edge harmonic oscillation (EHO) was observed in NSTX (K.F. Gan et al., 2017 Nucl. Fusion 57 126053). This paper introduces new analysis of the EHO impact on the divertor heat flux. It was found that enhanced edge turbulence significantly altered the divertor heat flux footprint in NSTX H-modes, relative to ELM-free discharges. When the background edge fluctuation level was low, the EHO significantly increased the heat flux width and decreased the divertor peak heat flux. When the background edge fluctuation level was high, the EHO actually increased the divertor peak heat flux. As a result, it was also found that the heat flux width increased with increasing frequency of the EHO.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Assessment of filtered cameras for quantitative 2D analysis of divertor conditions during detachment in JET L-mode plasmas

We report that estimates for 2D distributions of electron temperature, T e , electron density, n e , and atomic deuterium density, n o , in the JET divertor volume have been inferred from deuterium Balmer line intensity ratios obtained from tomographic reconstructions of divertor camera measurements. This enables also investigation of ionization, S ion , and recombination, S rec , rates. The analysis shows a decrease of T e to 0.5–1.0 eV throughout the outer divertor during detachment in low-confinement (L-mode) plasmas. Simultaneously, the high-n e region and the n 0 distribution in the outer divertor are observed to elongate and shift from the outer strike point towards the X-point. The observations are in qualitative agreement and follow the same sequence with modelling predictions of EDGE2D-EIRENE simulations of a density scan. While the method was found to provide good representation of the evolution of volumetric recombination during detachment, in agreement with the simulations, the movement of the ionization front upstream could not be followed due to lack of spatial overlap between the ionization region and the necessary emission distributions. Consequently, the representation of the ionization conditions and the particle balance in the detached outer divertor are compromised.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measurements of plasma parameters in the divertor island of Wendelstein 7-X through line-ratio spectroscopy on helium

We report electron temperature T e and density n e have been measured in the divertor island of Wendelstein 7-X by means of a new thermal He-beam system utilizing line-ratio spectroscopy. This system is equipped with a poloidal arrangement of five gas valves that allows to infer 2D plasma profiles T e ( R,Z ) and n e ( R,Z ) across the magnetic island forming the island divertor. A local maximum in T e ( R,Z ) was measured along the flux surfaces located between the island center and the outer separatrix. In the island center, characterized by closed field lines, a clear local minimum in T e was measured, indicating that no direct heat deposition occurs in this region and that the perpendicular heat flux from the core plasma into this domain is small. Increasing the island size leads to a reduction of the connection length in the scrape-off layer (SOL) and to a shift of the island center towards the divertor target. In this scenario, the T e ( R,Z ) are reduced in the entire island domain. Density profiles are found flat across the divertor island for both island size scenarios. These findings are supported by similar observations made with a reciprocating probe plunging through the same SOL island at the mid-plane of W7-X representing the upstream position in this 3D divertor geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Understanding detachment of the W7-X island divertor

The fundamental behavior of the W7-X island divertor under detached conditions, which has been theoretically predicted with the EMC3-Eirene code, is re-examined here under the experimental conditions achieved so far and compared with the first experimental results. Both simulations and experiments cover a range of divertor configurations and plasma parameters, and show the following common trends: (1) with rising impurity radiation, the target heat load decreases 'uniformly' over the entire target surface in the sense that both the peak and average heat loads can drop by an order of magnitude. Impurity radiation (mainly from intrinsic carbon) occurs primarily at the plasma edge and the resulting negative impact on the stored energy is less than 10%. (2) When the total radiation exceeds a critical level, the target particle flux (the recycling flux Γrecy) begins to fall and can drop by a factor of 3–5 at high radiation levels without an obvious indication of significant volume recombination. (3) While Γ recy decreases, the divertor neutral pressure continues to build up and reaches a maximum, at which point Γ recy has declined significantly. (4) During detachment, the electron temperature at the last closed flux surface falls in a way that is not quantitatively understandable from parallel classical heat conduction processes. This paper presents a physical explanation of the numerical/experimental results described above. Furthermore, using the EMC3-Eirene code as a diagnostic tool, we are able, apparently for the first time, to provide a full quantitative analysis of each transport channel in the island divertor, aiming to clarify how the island divertor plasma self-regulates to maintain particle, energy, and momentum balance under detached conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Toroidal modeling of 3D perturbations generated by current filaments in scrape-off layer in tokamak with biased divertor targets

The divertor biasing technique offers a promising alternative to control the edge localized mode (ELM) as well as the divertor heat load in tokamaks, as compared to the resonant magnetic perturbation (RMP) generated by magnetic coils. The linear resistive MHD code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) is employed to study the plasma response to the n = 2 (n is the toroidal mode number) field perturbation in HL-2A, produced by the divertor biasing current filaments flowing in the scrape-off layer region. The response field due to biasing currents is compared with the RMP field utilized for ELM control in HL-2A. The strength of the radial resonant field perturbation, produced by the biasing currents of 100 A level, is found to be comparable with RMP generated by several kA of ELM control coil currents for the reference plasma. The plasma normal displacement near the X-point and the associated neoclassical toroidal viscosity torque are also computed to be similar between these two techniques. The modeling results thus strongly suggest that the biasing technique can be applied to control ELMs. Moreover, the biasing currents produce field perturbations, including the plasma response, that are localized more near the plasma edge (compared to the RMP counterpart), thus reducing the chance of mode locking associated with core perturbations. Particle orbit tracing also reveals that the biasing current produced magnetic perturbation tends to widen the heat deposition region and induce the strike point splitting of the ion saturation flow on the outer divertor surface, consistent with experimental observations in HL-2A. These toroidal modeling results confirm the possibility of ELM control and plasma exhaust solution by the divertor biasing technique.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Divertor heat load estimates on NSTX and DIII-D using new and open-source 2D inversion analysis code

A thermography inversion algorithm has been developed in the open-source Python-based computer code, HYPERION, to calculate the heat flux incident on plasma-facing components (PFCs) in axisymmetric tokamaks. The chosen mesh size at the surface significantly affects the calculated transient heat flux results. The calculated transient heat flux will exceed the real value when the mesh size tends to zero but will underestimate the real value when the mesh size is large. A criterion for determining the appropriate mesh size for the transient heat flux calculation will be discussed. The numerical scheme for HYPERION uses a 2D fully implicit finite-difference approach, allowing temperature-dependent thermal properties of PFC materials. The inversion algorithm is benchmarked against established heat flux calculation codes, TACO and THEODOR, based on thermography data from NSTX and DIII-D respectively. The primary benefits of HYPERION compared to TACO and THEODOR are that it is open-source and it allows for the optimization of mesh thickness along the substrate. The algorithm also accounts for the thermal properties of thin surface layers that characteristically form on PFCs due to plasma-material interactions. The agreement between HYPERION and THEODOR is excellent, as the percent difference between the codes is ~5% on average in the case of the DIII-D data for moderate to high heat flux. Verification tests with TACO show slightly higher average percent differences of 8% and 12%. In using HYPERION to study filaments in heat flux, the initial results indicate that small ELMs filaments significantly broaden the divertor heat flux, and decrease divertor peak flux. Compared to the inter-ELM, the small ELM filaments decrease the divertor peak surface temperature. With intermittent divertor filaments, the divertor heat flux width is comparable with that found in L-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

FIREFLY: heat load and particle exhaust approximations for rapid evaluation of divertor designs

The divertor in a magnetic confinement fusion reactor is an essential component for power dissipation and particle removal. The FIREFLY package for rapid evaluation of divertor designs is presented as an extension of the FLARE code for field line reconstruction from a flux tube mesh. First, divertor loads are approximated with a simplified heat transport model. Neutralized particles are then sampled from the resulting load distribution, and the EIRENE code is used to track molecules and atoms in a plasma background while accounting for dissociation, charge exchange and ionization. Particles are removed on pumping surfaces in order to estimate the exhaust efficiency for a given divertor geometry. Optimization of the divertor geometry for more efficient particle exhaust is explored by using W7-X as an example, and the sensitivity to model parameters for the plasma background in the proxy calculations is evaluated.

mesh generation, magnetic field lines, scrape-off ↗

Thermal Hydraulics Analysis of a Divertor Monoblock Using SALAMANDER

Divertors are critical components in magnetic confinement fusion devices. One of the Divertor's crictial roles is absorbing the highest heat flux from the plasma. The Divertor accomplishes this through the use of divertor monoblocks and cooling channels. This work presents a thermal hydraulic analysis of a divertor monoblock using the Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering and Research (SALAMANDER). SALAMANDER is an open source tool developed at Idaho National Laboratory for conducting multiphysics and multscale analysis of Plasma Facing Components (PFCs). For the monoblock, solid heat conduction is modeled in the block with a constant heat flux of 1E7 W/m^2 applied to the top of the block to represent peak heat flux from the plasma. The cooling channel (which is pressurized water at a Reynolds number of ~1,000,000) is modeled using a k-epsilon turbulence model with standard wall functions. The heat transfer between the solid and fluid domain is modeled using the Dittus-Boelter correlation for the convective heat transfer coefficient. This work aims to show the importance of a multiphysics modeling approach when performing simulations on PFCs.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of Drift Effects in UEDGE Simulations of NSTX-U Edge Plasma With Lithium Divertors

Lithium is a low-Z material, and lithium-based plasma-facing components (PFCs) are planned for the National Spherical Torus Experiment Upgrade (NSTX-U) to explore potential benefits for divertor power exhaust and core plasma management. NSTX-U is a medium-sized spherical tokamak with up to 12 MW of auxiliary heating, capable of generating reactor-relevant plasma conditions. This work presents boundary plasma simulations for NSTX-U with lithium PFCs using the UEDGE code, incorporating full magnetic and 𝐄 ×𝐁 drift physics. The simulations show that drifts strongly influence heat and particle transport: they enhance convective transport, broaden the scrape-off layer heat-flux width 𝜆 𝑞 , and reduce the anomalous heat diffusivity 𝜒 required to reproduce predicted SOL heat-flux width. 𝐄 ×𝐁 drifts provide poloidal transport, while ∇𝐵 (which includes both gradB and curvature) drifts provide radial heat and particle transport. Lithium transport is also affected by drifts, with lithium ions migrating from the outer divertor to the inner divertor through the private flux region (PFR) following the 𝐄 ×𝐁 drifts flow, lowering upstream impurity lithium densities. UEDGE is self-consistently coupled with the Wall-Li model to study plasma lithium PFC interactions depending on the local lithium sourcing based on local plasma conditions and lithium surface temperature. In these simulations, lithium evaporation shows a vapor-shielding effect that reduces divertor heat flux and increases radiative losses once surface temperatures exceed 450°C. This research work provides a first step toward self-consistent modeling of lithium PFCs in NSTX-U, demonstrating the impact of drift-driven plasma transport in SOL and divertor regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Liquid lithium divertor analysis using coupled plasma material interaction model

A liquid lithium divertor can improve performance of future fusion devices by creating efficient power exhaust and improving the energy confinement via pumping of the hydrogen isotopes. In addition, significantly higher heat fluxes can be handled if controlled vapor shielding is used to redistribute the divertor heat flux over a wider area. Design and optimization of such a system calls for an analysis model which includes a strong two-way coupling between the plasma and divertor material. The incoming plasma heat and particle flux will affect the divertor surface temperature, which is a defining factor of the lithium evaporative and sputtered flux going into the plasma. Results of the coupled model based on the plasma edge code SOLPS-ITER and the computational fluid dynamics (CFD) code ANSYS-CFX will be presented for different configurations. An analytical slab flow model is used as a heat transfer boundary condition for SOLPS, defining particle flux from the wall via calculation of the surface temperature. At the final step, results of the SOLPS analysis are verified using a 3D CFD magnetohydrodynamics (MHD) analysis which uses heat and particle flux from SOLPS as a boundary condition. In addition to plasma heat flux, both analytical and CFD temperature models include several plasma material interaction effects, such as lithium evaporation, condensation and sputtering based on deuterium target flux. New adatom sputtering model based on the available experimental data is presented. Analytical model is expanded to include free surface axisymmetric configurations. Results of parametric studies of the divertor configurations with different lithium inlet temperature and velocity will be presented leading to the optimal design resulting in the lowest possible lithium contamination in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Quantifying heat and particle flux to primary and secondary divertors for various types of edge-localized-modes

Between 15% and 30% of edge-localized-mode (ELM) heat flux can be deposited to regions outside the main divertor including the far scrape-off-layer (SOL), private flux region (PFR), and secondary divertor inner target. Analysis shown here demonstrates that type-I and small ELM plasmas are transported to the PFR and to the secondary inner divertor, which is magnetically isolated from the outer divertor, leading to well-defined heat flux profiles and with peak values that can surpass those at the secondary outer target. Such features are consistently observed for pedestal collisionalities from $ν$$^{*}_{e}$ ~1.5 to $ν$$^{*}_{e}$ ~3.9. Heat flux profiles of the examined ELM types feature rippled structures due to bursts in the outer far-SOL region but not on the secondary inner target, causing long decay lengths in the time-averaged ELM profiles. The contribution of each ELM type to the total time-averaged power deposited to the secondary divertor has been evaluated, showing that grassy ELMs contribute ~8%, small ELMs ~67%, and type-I ELMs ~85%. These findings imply that small ELMs may yet pose a concern for future machines if some regions of the main wall are not designed to withstand significant heat and particle fluxes. Due to the low intra-ELM heat flux contribution, however, the grassy ELM regime is an attractive option for an ELMing scenario in future machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulations of the churning mode: Toroidally symmetric plasma convection and turbulence around the X-points in a snowflake divertor

Using a reduced MHD model, extended to include field-aligned thermal conduction, we present numerical simulations of the churning mode (CM): a toroidally symmetric, non-linear plasma vortex in the vicinity of the null points in a snowflake (SF) divertor [D. D. Ryutov et al., Phys. Scr. 89, 088002 (2014)]. Simulations are carried out across a range of inter-null separations, d xx ⁠, and inter-null orientations, θ⁠, primarily in conditions relevant to the MAST-U tokamak. We find that, when d xx is small, the CM induces additional transport across the X-points when β pm ≳8%, where β pm is the ratio of the plasma pressure in the null region to poloidal magnetic pressure at the midplane. This transport also increases approximately linearly as d xx is reduced. A diffusive model of this transport is shown to predict the total transport across the null points, where diffusion coefficients of up to ∼ 10 2 m 2 s −1 centered on a small region around the X-points are used. However, the CM also results in significant changes to the flux surfaces in the null region, which is not captured by this diffusive model. The changes in magnetic geometry mean the fractional exhaust power delivered to each divertor leg is highly sensitive to β pm ⁠, d xx ⁠, and θ⁠. For small values of θ⁠, the CM can induce a change in topology, redirecting exhaust power from a secondary divertor leg on the high field side to one on the low field side. Similar behavior is found in the fraction of exhaust power going to the inner and outer divertor. Such changes in the flux surfaces may not be captured by Grad–Shafranov solvers and so may be a source of error in the magnetic reconstruction of SF experiments. We consistently find that the fractional exhaust power going to a secondary divertor leg on the high field side is small, consistent with SF experiments.

Computer simulation↗

Divertor detachment characterization in negative triangularity discharges in DIII-D via 2D edge-plasma transport modeling

Edge fluid modeling of the first divertor-plasma detachment experiments in negative triangularity (NT) discharges in the DIII-D tokamak is presented using the 2D multifluid edge transport code UEDGE, including cross-field particle drifts. Experiments were performed where the lower single-null magnetic equilibrium had a strong NT (δ≈−0.5), that is, where the magnetic X-point is at a larger major radius than the core magnetic axis. Divertor-plasma detachment was induced by increasing the core plasma density in DIII-D via intrinsic gas puffing. Here density scans are performed with UEDGE to reach a detached plasma and to quantitatively recover the experimental rollover of the ion saturation current on the outer divertor target plate. The simulations cover experiments with both signs of the toroidal magnetic field, B T , where the ion magnetic Grad-B drifts are directed into (forward B T ) and out of (reverse B T ) the divertor region. Consistent with experiments with neutral beam power injection, the NT simulations reproduce: 40% higher density is needed to reach detachment onset with forward B T compared with reverse B T , and the absence of deep detachment is found with reverse B T . Similarly, comparison between Ohmic discharges in NT and positive triangularity (PT) shaping confirms that a substantially higher density is needed to achieve detachment in NT than in PT, with NT requiring an line-average density of at least the Greenwald density or higher. Simulation results suggest that higher densities are needed to reach detachment in negative compared to PT because these discharges have a shorter midplane-to-target distance along the total magnetic field B, a shorter outer divertor poloidal leg length (0.06 m vs 0.2 m), and reduced radial transport [near-scrape-off layer (SOL) D ⊥ /χ ⊥ =0.3/0.5 vs D ⊥ /χ ⊥ =1.0/1.0, all in (m 2 /s)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental investigation of a closed vapour box module for a divertor-like configuration in Magnum-PSI

Efficient management of extreme heat fluxes in the divertor region to extend the lifetime of the components remains a critical challenge for the realization of nuclear fusion-based power plants. Among the alternative concepts explored for the divertor region, the use of liquid metals, particularly lithium, is of interest due its ability to dissipate the incoming plasma heat flux through the vapour shielding effect (VS). In this work, we experimentally investigated a ‘closed’ configuration of a dedicated Vapour Box Module (VBM) in the linear plasma device Magnum-PSI. The goal of the experiments is to simulate the vapour box divertor environment conditions and assess its performance in terms of power mitigation and redistribution and lithium confinement. Initial testing without Li demonstrated the efficacy of a closed VBM structure in inducing detachment via neutral gas accumulation. Apertures which enabled non-condensing gas to be effectively pumped while ensuring lithium condensed on the inner surfaces were therefore added. With a lithium capillary porous structure target used, lithium is directly vaporized by the plasma, forming a dense lithium vapour cloud that interacts with the incoming plasma. This resulted in a significant reduction of the target temperature of at least 48%, together with a temperature locking effect, a phenomenon typically observed in the VS regime. Lithium vapour confinement within the VBM was strongly correlated with the wall temperature. Relatively cold walls promoted Li re-condensation and therefore improved Li confinement, although with the expected trade-off of increased hydrogenic retention on lithium-wetted surfaces. As the wall temperature increased, the confinement efficiency decreased, consistent with reduced Li re-condensation and thermally activated Li–H chemistry and remobilization at the walls. Diagnostic measurements through embedded thermocouples and calorimetry revealed that lithium vaporization and re-condensation processes also playedsignificant roles in plasma power dissipation. The results advance the case for a closed divertor chamber with direct lithium evaporation from the strike-points as a viable method to manage divertor heat fluxes in future fusion reactors.

Romano, Fabio [Dutch Institute for Fundamental Ene↗

ELM sputter erosion modeling of a tungsten coated small angle slot divertor in DIII-D

We modeled plasma edge localized mode (ELM) sputter erosion for a Small Angle Slot divertor with a tungsten coated region (SAS-VW), designed for experiments in the DIII-D tokamak, and proposed for use in future advanced tokamaks. The simulations use a free-streaming, 1000 eV, C +6 and D +1 ELM impingement model, with SOLPS-ITER, ITMC-DYN, and REDEP/WBC code packages for background plasma, material response, and erosion/redeposition respectively. The results show ELM'ing plasma gross and net tungsten erosion fluxes of the mixed-material C/W surface peaking at the slot entrance region, and an order of magnitude higher than for non-ELMs. The per-pulse erosion, however, remains low, of order 0.5 nm, due to expected moderate ELM frequencies and duration in DIII-D. Here, the ELMs result in a ~25x higher peak sputtered W current leaving the divertor slot region, towards the core plasma, compared to the ELM-free plasma case. The time-integrated escape current, however, may not significantly affect core plasma high-Z contamination concerns, for a 1% ELM duty factor, but may be an issue for higher frequency ELMs. In general, the modeling results appear favorable for effective testing of the SAS-VW divertor in DIII-D, and extrapolation to innovative divertor designs in future ITER-like and DEMO fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗