Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first wall”

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 55 records · Page 3

Integrated modelling: Coupling of surface evolution and plasma-impurity transport

During the interaction of the scrape off layer (SOL) plasma with the first wall the evolution of both wall and plasma are tightly coupled: The erosion of the first wall leads to an impurity concentration in the plasma which affects the particle and power balance in the plasma. In turn the impurities, when leaving the plasma via transport, can form deposits and mixed materials, far away from their initial source location. These deposits can be eroded, allowing the impurities to stepwise migrate through the fusion device until they end up at a location where the plasma at the wall is cold enough and no further erosion occurs. To describe these processes an integrated model of surface evolution and plasma transport of impurities is needed. The WallDYN code achieves this required coupling of processes by parameterising the output of surface evolution- and plasma impurity-migration-codes by analytical models. For a given fixed background plasma it evolves the surface composition, derives impurity flux into and from the plasma and can from this derive the impurity densities in the plasma. This paper will show the importance of including this recycling of impurities at the wall in impurity migration modelling: The 13 CH 4 seeding experiment performed in the Wendelstein 7-X Stellerator is modelled using the recent extension of the WallDYN code to 3D plasma and wall geometries. A comparison with post mortem analysis of the 13 C deposition shows both qualitative and quantitative agreement with the WallDYN calculations.

36 MATERIALS SCIENCE↗

Thermal energy transformer

For use in combination with a heat engine, a thermal energy transformer is presented. It is comprised of a flux receiver having a first wall defining therein a radiation absorption cavity for converting solar flux to thermal energy, and a second wall defining an energy transfer wall for the heat engine. There is a heat pipe chamber interposed between the first and second walls having a working fluid disposed within the chamber and a wick lining the chamber for conducting the working fluid from the second wall to the first wall. Thermal energy is transferred from the radiation absorption cavity to the heat engine.

Berdahl, C. M.↗

Reassessing energy deposition for the ITER 5 MA vertical displacement event with an improved DINA model

The beryllium (Be) main chamber wall interaction during a 5 MA/1.8 T upward, unmitigated VDE scenario, first analysed in [J. Coburn et al., Phys. Scr. T171 (2020) 014076] for ITER, has been re-evaluated using the latest energy deposition analysis software. Updates to the DINA disruption model are summarized, including an improved numerical convergence for the 0D power balance, limitations on the safety factor within the plasma core, and the choice to maintain a constant plasma +halo poloidal cross-section. Such updates result in a broad halo region and higher radiated power fractions compared to previous models. The new scenario lasts for ~75 ms and deposits ~29 MJ of energy, with the radial distribution of parallel heat flux q‖(r)resembling an exponential falloff with an effective λ_q=75-198 mm. A maximum halo width w_h of 0.52 m at the outboard midplane is observed. SMITER field line tracing and energy deposition simulations calculate a q_{⊥,max} of ~83 MW/m^2 on the upper first wall panels (FWP). Heat transfer calculations with the MEMOS-U code show that the FWP surface temperature reaches ~1000 K, well below the Be melt threshold. Variations of this 5 MA scenario with Be im-purity densities from 0 to 3∙10^19 m^-3 also remain below the melt threshold despite differences in energy deposition and duration. These results are in contrast to the early study which predicted melt damage to the first wall [J. Coburn et al., Phys. Scr. T171 (2020) 014076], and emphasize the importance of accurate models for the halo width w_h and the heat flux distribution q‖(r)within that halo width. The 2020 halo model in DINA has been compared with halo current experiments on COMPASS, JET, and Alcator C-Mod, and the preliminary results build confidence in the broad halo width predictions. Results for the 5 MA VDE are compared with those for a 15 MA equivalent, generated using the new DINA model. At the higher current, significant melting of the upper FWP is to be expected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Airfoil having internally cooled wall with liner and shell

An airfoil includes an airfoil section that has an airfoil wall that surrounds a cavity. The airfoil wall includes a first wall section comprised of a monolithic ceramic shell, a second wall section comprised of a CMC liner that is bonded along an interface to an interior side of the monolithic ceramic shell, and a ceramic barrier coating that is bonded to an exterior side of the monolithic ceramic shell. The airfoil wall includes at least one of: inter-section cooling channels that extend along the interface and are bounded by the monolithic ceramic shell and the CMC liner, or intra-section cooling channels that extend within the CMC liner.

Croteau, Paul F.↗

Low- n global ideal MHD instabilities in the CFETR baseline scenario

This article reports an evaluation on the linear ideal magnetohydrodynamic (MHD) stability of the China Fusion Engineering Test Reactor (CFETR) baseline scenario for various first-wall locations. The initial-value code NIMROD and eigen-value code AEGIS are employed in this analysis. Despite the distinctly different approaches in modeling the scrape off layer region, the dominant growth in each of the low-n (n=1–10) modes are consistent between the two codes. We report the higher-n modes are dominated by ballooning modes and localized in the pedestal region, while the lower-n modes have more prominent external kink components and broader mode profiles. The influence of the plasma-vacuum profile and wall shape are also examined using NIMROD. In the presence of a resistive wall, the low-n ideal MHD instabilities are further studied using AEGIS. For the designed first-wall location, the n=1 resistive wall mode is found to be unstable, which could be fully stabilized by uniform toroidal rotation above the 2.9% core Alfven speed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Lithium Divertor Targets and Walls for the ASTER Liquid Stellarator Reactor, Distributed Divertor

Stellarators may have advantages for certain liquid metal options as Plasma Facing Components (PFC) for divertor targets and first walls due to the wide range of possible magnetic configurations, which additionally are free of disruptions and fast field variations. In a previous work (V. Queral et al., IEEE Trans. Plasma Sci. 52, 2024), a concept of stellarator reactor (ASTER-CP) based on swirling Li-molten salts and liquid lithium floating on the molten salt as PFC was presented. The divertor matters were not studied then and, thus, they are being studied and experimentally tested now. The ASTER-CP reactor concept, the initial liquid metal experiments and potential concepts for the ASTER-CP divertor and first wall are reported. Concerning the experiments, several small scale experiments of galinstan in a small rotating cylinder under magnetic field have been produced, including one experiment with high viscosity galinstan-mixture for increased thickness of layer. An experiment of floating lithium on the molten salt LiCl-PbCl2 gave fast volatilization/decomposition of the molten salt. Particularly for divertors, the traditional free-flow, Capillary Porous Systems and ‘divertorlets’ have been studied for application to ASTER-CP. Surface waves (hot spots), lack of enough surface fluid turbulence and excessive fluid speed are the main issues found in fast free-flow. The perhaps original concept of Distributed Divertor and Equi-power Surface is tentatively proposed and studied, taking advantage of stellarator fields and low recycling regime.

distributed divertor↗

Q3 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q3: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis. • Benchmark between two first wall heat flux mapping methods, identify importance of various heat flux sources and physics impact of using fully coupled CESOL vs post-analysis evaluation. 2. Generate medium fidelity parametrized CAD. • Generate parametrized CAD components for the CAT example case via either user-defined modules called within the geometry generation or by defeatured/parametrized CAD, including DCLL blanket matched to divertor boundary and magnets. Define materials, labels, and boundary conditions for passing the mesh to CFD tools. 3. Demonstrate multiphysics magnet analysis. • Demonstrate magnet analysis workflow called from the FREDA workflow, and 4. Demonstrate nuclear analysis. • Add model to OpenFOAM and/or other codes possibly including Diablo to account for tritium diffusion in solids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Empirical scaling of the L–H threshold power for metal wall tokamaks using a multi-device database

The empirical scaling for the H-mode power threshold in tokamaks has been revisited using a database with threshold data from machines with a metallic first wall as part of International Tokamak Physics Activity (ITPA) task TC-26. The database contains discharges from ASDEX Upgrade (AUG) (W), JET (Be/W) and Alcator C-Mod (Mo). This was motivated by reports that in like-for-like discharges the power threshold was reduced by approximately 30% after the change from carbon based to metallic first wall materials on AUG (Ryter et al 2013 Nucl. Fusion 53 113003) and JET (Maggi et al 2014 Nucl. Fusion 54 023007). The database contains L–H transition data for all hydrogen isotopes and mixtures, including T and DT from the recent JET campaigns. Compared to the ITPA 2008 scaling (Martin et al 2008 J. Phys.: Conf. Ser. 123 012033), the metal wall scaling has a smaller magnetic field exponent but a larger density exponent. We present an additional parameter to capture the strong dependence of the L–H power threshold (approx. factor 2) on the magnetic configuration in the divertor on JET. The scaling recovers the approximate inverse isotope mass scaling of the threshold power. Alternative scalings involving the plasma current and poloidal magnetic field are explored. Despite the reduction in threshold observed earlier, the scalings based on the metal wall database do not necessarily extrapolate to a lower threshold for ITER compared to the ITPA 2008 scaling, especially at high density. The divertor configuration effect induces the largest uncertainty in the extrapolation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of neoclassical tearing modes and toroidal field ripple on lost alpha power in the SPARC tokamak

Here, using the SPIRAL Monte Carlo, full particle-orbit simulation code, we investigate the effects of neoclassical tearing modes (NTMs) and toroidal field (TF) ripple on alpha power losses during steady-state operation of the SPARC primary reference discharge. Model perturbations for TF ripple and the $m/n = 2/1$ and $3/2$ NTMs with exaggerated widths selected based on an H-mode plasma approaching thermal quench are added to a simulated SPARC magnetic equilibrium through which marker particles are tracked. The 3/2 and 2/1 NTMs are located at $\rho_{\mathrm{pol}} \sim 0.76$ and $\rho_\mathrm{pol} \sim 0.86$ respectively, well positioned to increase alpha particle transport into and within an outer lossy region of the plasma beyond $\rho_\mathrm{pol} \sim 0.8$ where over 95% of lost alpha particles are born. Total alpha power losses are shown to increase modestly from 1.73% lost at a minimum to 2.34% lost at a maximum, and alpha particle surface power densities form localized hotspots on the first-wall near the low-field side midplane due to NTMs and TF ripple. We establish a conservative upper limit for first-wall alpha surface power densities on a toroidally symmetric wall for typical, flattop operation and motivate the consideration of NTMs in the design of three dimensional limiter surfaces for SPARC.

alpha power losses↗

Fluid mixing apparatus using liquid fuel and high- and low-pressure fluid streams

A fluid mixing apparatus includes mixing conduits that extend through a fluid plenum. The fluid plenum, which surrounds a first wall defining a main passage fluidly coupled to a low-pressure fluid source, is surrounded itself by a second wall defining a high-pressure plenum fluidly coupled to a high-pressure fluid source. An insulated tube disposed at the inlet of the first wall delivers a third fluid. The mixing conduits fluidly couple the high-pressure plenum to the main passage, where the high-pressure fluid is mixed with low-pressure fluid and the third fluid. Optionally, the fluid plenum may house a fourth fluid that is injected through injection holes in the mixing conduits. The fluid mixing apparatus may be used to mix one or more fuels with high- and low-pressure air in a gas turbine combustor. Alternately, the fluid mixing apparatus may mix a fluid with high- and low-pressure water streams.

Berry, Joseph Jonathan↗

The new ITER baseline, research plan and open R&D issues

A new baseline (NB) has been proposed by the ITER Project to ensure a robust achievement of the Projects’ goals, in view of past challenges including delays incurred due to the Covid-19 pandemic, technical challenges in completing first-of-a-kind components and in nuclear licensing. The NB includes modifications to the configuration of the ITER device and its ancillaries (e.g. change from beryllium to tungsten as first wall material, modification of the heating and current drive mix, etc.) as well as additional testing of components (e.g. toroidal field coils) or phased installation (start with inertially cooled first wall before later installation of the final actively water-cooled components) to minimise operational risks. In the NB, the ITER research plan (IRP) will be divided into three main phases: (a) start of research operation, with 40 MW of ECH and 10 MW of ICH, which will focus on the demonstration of 15 MA operation in L-mode, commissioning of all required systems, including disruption mitigation, and the demonstration of H-mode plasma operation in deuterium; (b) DT-1, with 60–67 MW of ECH, 33 MW of neutral beam injection (NBI) and 10–20 MW of ICH, which will demonstrate robust operation in high confinement H-mode plasmas in DT up to Q ⩾ 10 and for burn durations of 300–500 s within an accumulated neutron fluence of ∼1% of the ITER machine’s lifetime total, and; (c) DT-2, with up to 67 MW of ECH, up to 49.5 MW of NBI and up to 20 MW of ICH, with the ITER tokamak and ancillaries in their final configuration to demonstrate routine operation in DT plasmas at high Q and the Q ⩾ 5 long-pulse and steady-state scenarios to the final neutron fluence and to perform R&D on nuclear fusion reactor issues. The logic, physics basis, modelling and experimental evaluations carried out to support the NB and the associated IRP are described. These include the impact of the tungsten wall on plasma scenarios and associated risk mitigation measures, as well as the optimisation of the tokamak components and ancillaries to minimise Project risks. Open R&D issues related to these evaluations and mitigation measures are also described together with experimental, modelling and validation activities required to address them.

ITER↗

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetic mesh generation and field line reconstruction for scrape-off layer and divertor modeling in stellarators

The design of divertor targets and baffles for optimal heat and particle exhaust from magnetically confined fusion plasmas requires a combination of fast, low-fidelity models (such as EMC3-Lite (Feng 2022 Plasma Phys. Control. Fusion 64 125012)) for scoping studies and high-fidelity ones (such as EMC3-EIRENE (Feng et al 2014 Contrib. Plasma Phys. 54 426–31)) for verification. Both of those approaches benefit from a magnetic flux tube mesh for fast interpolation and mapping of field line segments (Feng et al 2005 Phys. Plasmas 12 052505). A new automated mesh generator for unstructured quadrilateral flux tubes with adaptive refinement is presented and integrated into FLARE (Frerichs 2024 Nucl. Fusion 64 106034). For HSX with an extended first wall, it is found that several layers of flux tubes can span the entire half field period before splitting is required. This is an advantage over the traditional setup of the EMC3-EIRENE mesh where careful construction of several sub- domains is required already for the much tighter present first wall. In particular, there is no longer the need to manually construct a suitable outer boundary for the mesh. The divide and conquer paradigm with unstructured mesh layout offers a powerful alternative for fast head load approximation that is suitable for integration into optimization workflows. Further examples for W7-X and CTH demonstrate the versatile application range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Integral analysis of the effect of material dimension and composition on tokamak neutronics *

Abstract The neutronics performance of a tokamak has been identified as an important factor in designing a fusion power plant. The design of the tokamak should not only meet operational parameters such as sufficient tritium breeding, but also safety parameters such as low structural material activation. This paper investigates the impacts of the neutronics metrics for the ARC-class tokamak, a compact tokamak with an immersion blanket, by perturbing the first five layers of structural material—first wall, inner vacuum vessel, coolant salt channel, neutron multiplier, and outer vacuum vessel. The goal of this work is to provide insight into shaping and scaling the flux on each layer to obtain optimized operational and safety metrics through quantification of the responses from each perturbation. Results show that increased first wall thickness can increase the tritium breeding ratio (TBR) in specific configurations with high 6 Li enrichments and that vacuum vessels decrease TBR for low- 6 Li enrichment configurations. It was also found that the neutron multiplier can either increase or decrease TBR depending on the configuration. The response of metrics to the change in layer thickness and enrichment also varies depending on the vacuum vessel material. The integral impacts of 6 Li enrichment, layer thicknesses, and vacuum vessel material choice are investigated and presented in this paper.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Wall heating by subcritical energetic electrons generated by the runaway electron avalanche source *

Abstract Subcritical energetic electrons (SEEs) produced by the runaway electron (RE) avalanche source at energies below the runaway threshold are found to be the primary contributor to surface heating of plasma-facing components (PFCs) during final loss events. This finding is supported by theoretical analysis, computational modeling with the Kinetic Orbit Runaway electrons Code (KORC), and qualitative agreement with DIII-D experimental observations. The avalanche source generates significantly more secondary electrons below the runaway threshold, which thermalize rapidly when well-confined. However, during a final loss event, the RE beam impacts the first wall, and SEEs are deconfined before they can thermalize. Additionally, because the energy deposition length decreases faster than energy, the deposited energy density, and thus the maximum PFC surface temperature change, is larger for SEEs than REs. KORC simulations employ an analytic first wall to model particle deconfinement onto a non-axisymmetric wall composed of individual tiles. PFC surface heating is calculated using a 1D model extended to include an energy-dependent deposition length scale. Simulations of DIII-D qualitatively agree with infrared (IR) imaging only when SEEs from the avalanche source are included. These results demonstrate that SEEs are the dominant contributor to PFC surface heating and indicate that the avalanche source plays a critical role in the PFC damage caused during final loss events. The prominence of SEEs also has important implications for interpreting IR imaging, one of the primary diagnostics for RE-wall interaction diagnosis, despite REs dominating the energy and current density. This result improves predictions of wall damage due to post-disruption REs to estimate material lifetime and design RE mitigation systems for ITER and future reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Q1 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q1: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis: $/circ$ Initially run CESOL with SOLPS and map heat flux using simple HEAT method: ▪ Generate SOLPS grid for CAT reference case, ▪ Implement and apply ‘lower’ fidelity method of using HEAT-like analytic method to map SOLPS heat flux from charged particles to first wall, and ▪ Provide initial heat flux + neutron loading for further engineering analysis. 2. Demonstrate multiphysics magnet analysis: $/circ$ Demonstrate validation of the Elmer workflow by comparing induced stresses due to EM forces generated by TF coils using multiple codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Corrigendum: Simulations of stand-off runaway electron beam termination by tungsten particulates for tokamak disruption mitigation (2024 Nucl. Fusion 64 056019)

In the original article, we inadvertently omitted a reference which previously introduced the concept of tungsten injection for 'Disruption Mitigation in Tokamak Reactor via Reducing the Seed Electrons of Avalanche.' In that reference, the authors have proposed injection of a tungsten cylinder of 80 mm by 8 mm dimensions, by rail guns, to absorb the runaway seed population. This can be contrasted with another interesting idea of using small tungsten pellets coated with a low Z material for depleting the runaway seeds as an option for upgrade of the ITER Disruption Mitigation System. Interested readers are referred to those two papers on the tungsten injection for runaway seed removal/reduction, and the related issue of not shortening the current quench excessively. The application focus of was not on reducing the runaway seeds in the initial Ohmic-to-runaway current conversion phase, but on safe termination of a fully formed and likely decaying runaway beam that is about to scrape off against the first wall, which can be accelerated by vertical displacement events. This scheme thus serves as a last line of defense against potential wall damage. The choice here is locally released tungsten particulates, similar to the previous idea of using tungsten particulates as a dust shield for the divertor. The most significant finding of is that strong pitch angle scattering, in addition to energy attenuation and absorption, of high-energy runaways by tungsten particulates, provides a transport mechanism by which runaways damage of the first wall can be mitigated by reducing and spreading the runaway wall load. We regret the omission in the original article, and thank the authors of the cited article for bringing this matter to our attention. We also find the possibility of using the injected tungsten rod from to terminate the runaways, as opposed to the cloud of tungsten particulates in intriguing. In that context, it is of interest to mention another possibility of a retractable tungsten metal arm that can be swung out for deployment in runaway termination, which would remove the need for post-mitigation recovery of tungsten debris in injection schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydrogen isotope analysis in W-tiles using fs-LIBS

Abstract Laser-Induced Breakdown Spectroscopy (LIBS) is a promising technology for in-situ analysis of Plasma-Facing Components in magnetic confinement fusion facilities. It is of major interest to monitor the hydrogen isotope retention i.e. tritium and deuterium over many operation hours to guarantee safety and availability of the future reactor. In our studies we use ultraviolet femtosecond laser pulses to analyze tungsten (W) tiles that were exposed to a deuterium plasma in the linear plasma device PSI-2, which mimics conditions at the first wall. A high-resolution spectrometer is used to detect the Balmer- $$\alpha$$ α transition of the surface from implanted hydrogen isotopes (H and D). We use Calibration Free CF-LIBS to quantify the amount of deuterium stored in W. This proof-of-principle study shows the applicability of femtosecond lasers for the detection of low deuterium concentration as present in first wall material of prevailing fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗