Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “IMPINGEMENT”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Cosmic rays in molecular clouds probed by H 2 rovibrational lines

Context. Low-energy cosmic rays (<1 TeV) play a fundamental role in the chemical and dynamical evolution of molecular clouds, as they control the ionisation, dissociation, and excitation of H 2 . Their characterisation is therefore important both for the interpretation of observations and for the development of theoretical models. However, the methods used so far for estimating the cosmic-ray ionisation rate in molecular clouds have several limitations due to uncertainties in the adopted chemical networks. Aims. Here, we refine and extend a previously proposed method to estimate the cosmic-ray ionisation rate in molecular clouds by observing rovibrational transitions of H 2 at near-infrared wavelengths, which are mainly excited by secondary cosmic-ray electrons. Methods. Combining models of interstellar cosmic-ray propagation and attenuation in molecular clouds with the rigorous calculation of the expected secondary electron spectrum and updated electron-H 2 excitation cross sections, we derive the intensity of the four H 2 rovibrational transitions observable in cold dense gas: (1—0)O(2), (1—0)Q(2), (1—0)S(0), and (1—0)O(4). Results. The proposed method allows the estimation of the cosmic-ray ionisation rate for a given observed line intensity and H 2 column density. We are also able to deduce the shape of the low-energy cosmic-ray proton spectrum impinging upon the molecular cloud. In addition, we present a look-up plot and a web-based application that can be used to constrain the low-energy spectral slope of the interstellar cosmic-ray proton spectrum. We finally comment on the capability of the James Webb Space Telescope to detect these near-infrared H 2 lines, which will make it possible to derive, for the first time, spatial variation in the cosmic-ray ionisation rate in dense gas. Besides the implications for the interpretation of the chemical-dynamic evolution of a molecular cloud, it will finally be possible to test competing models of cosmic-ray propagation and attenuation in the interstellar medium, as well as compare cosmic-ray spectra in different Galactic regions

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Investigation of flow regimes in arc plasma–gas interactions using a two-temperature arc in crossflow model

The perpendicular impingement of a gas stream on an electric arc, a configuration known as the arc in crossflow, is of primary relevance in the study of plasma-gas interactions, as well as in industrial applications such as circuit breakers and wire-arc spraying. The flow dynamics in the arc in crossflow are the result of coupled fluid-thermal-electromagnetic phenomena accompanied by large property gradients, which can produce significant deviations from Local Thermodynamic Equilibrium (LTE) among electrons and gas species. These characteristics can lead to the establishment of distinct flow regimes depending on the relative values of the controlling parameters of the system, such as inflow velocity, arc current, and inter-electrode spacing. A two-temperature non-LTE (NLTE) model is used to investigate the arc dynamics and the establishment of flow regimes in the arc in crossflow. The plasma flow model is implemented within a nonlinear Variational Multiscale (VMS) numerical discretization approach that is less dissipative, and hence better suited to capture unstable behavior, than traditional VMS methods commonly used in computational fluid dynamics simulations. The Reynolds and the Enthalpy dimensionless numbers, characterizing the relative flow strength and arc strength, respectively, are chosen as the controlling parameters of the system. Simulation results reveal the onset of dynamic behavior and the establishment of steady, periodic, quasi-periodic, and chaotic or potentially turbulent, regimes, as identified by distinct spatiotemporal fluctuations. The computational findings reveal the role of increasing the relative arc strength on enhancing flow stability by delaying the growth of fluctuating and unstable flow behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A graphene-based magnetoplasmonic metasurface for actively tunable transmission and polarization rotation at terahertz frequencies

The rise of terahertz (THz)-based technologies has prompted intense interest in methods for manipulating the properties of THz radiation. In this Letter, we experimentally demonstrate a graphene micro-ribbon-based metasurface that allows us to tailor and dynamically tune the transmission and polarization rotation spectra of impinging THz pulses. This is achieved by controlling the geometry of the ribbon pattern and the applied DC magnetic field, both of which strongly influence the magnetoplasmonic response of the metasurface. The properties of the metasurface can be accurately modeled using an effective medium approach that considers the capacitive interaction of the ribbons resulting from the THz electric field. Furthermore, our results highlight the promise of graphene-based THz devices that leverage the strongly enhanced light–matter interactions resulting from patterning.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Spectral solar irradiance on inclined surfaces: A fast Monte Carlo approach

Estimating spectral plane-of-array (POA) solar irradiance on inclined surfaces is an important step in the design and performance evaluation of both photovoltaic and concentrated solar plants. This work introduces a fast, line-by-line spectral, Monte Carlo (MC) radiative transfer model approach to simulate anisotropic distributions of shortwave radiation through the atmosphere as photon bundles impinge on inclined surfaces. This fast Monte Carlo approach reproduces the angular distribution of solar irradiance without the undesirable effects of spatial discretization and thus computes detailed POA irradiance values on surfaces at any orientation and also when surfaces are subjected to the anisotropic ground and atmospheric scattering. Polarization effects are also easily incorporated into this approach that can be considered as direct numerical simulation of the physics involved. Here, we compare our Monte Carlo radiative transfer model with the most widely used empirical transposition model, Perez4, under various conditions. The results show that the Perez4 model reproduces the more detailed Monte Carlo simulations with less than 10% deviation under clear skies for all relevant surface tilt and azimuth angles. When optically thin clouds are present, observed deviations are larger, especially when the receiving surface is strongly tilted. Deviations are also observed for large azimuth angle differences between the receiving surface and the solar position. When optically thick clouds are present, the two models agree within 15% deviation for nearly all surface orientation and tilt angles. The overall deviations are smaller when compared with cases for optically thin clouds. The Perez4 model performs very well (6.0% deviation) in comparison with the detailed MC simulations for all cases, thus validating its widespread use for practical solar applications. When detailed atmospheric profiles and cloud optical properties are available, the proposed fast Monte Carlo radiative model reproduces accurate spectral and angular POA irradiance levels for various atmospheric and cloud cover conditions, surface orientations, and different surface and ground properties.

42 ENGINEERING↗
DOE OSTI2022

Characterization of self-magnetic pinch (SMP) radiographic diode performance on RITS-6 at Sandia National Laboratories. I. Diode dynamics, DC heating to extend radiation pulse

Radiographic diodes focus on an intense electron beam to a small spot size to minimize the source area of energetic photons for radiographic interrogation. The self-magnetic pinch (SMP) diode has been developed as such a source and operated as a load for the six-cavity radiographic integrated test stand (RITS-6) inductive voltage adder driver. While experiments support the generally accepted conclusion that a 1:1 aspect diode (cathode diameter equals anode–cathode gap) delivers optimum SMP performance, such experiments also show that reducing the cathode diameter, while reducing spot size, also results in reduced radiation dose, by as much as 50%, and degraded shot reproducibility. Analysis of the effective electron impingement angle on the anode converter with time made possible by a newly developed dose-rate array diagnostic indicates that fast-developing oscillations of the angle are correlated with early termination of the radiation pulse on many of the smaller-diameter SMP shots. This behavior as a function of relative cathode size persists through experiments with output voltages and currents up to 11.5 MV and 225 kA, respectively, and with spot sizes below approximately few millimeters. Since simulations to date have not predicted such oscillatory behavior, considerable discussion of the angle behavior of SMP shots is made to lend credence to the inference. We report there is clear anecdotal evidence that DC heating of the SMP diode region leads to stabilization of this oscillatory behavior. This is the first of two papers on the performance of the SMP diode on the RITS-6 accelerator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A reaction mechanism for plasma electrolysis of AgNO 3 forming silver nanoclusters and nanoparticles

In plasma-driven solution electrolysis (PDSE), gas-phase plasma-produced species interact with an electrolytic solution to produce, for example, nanoparticles. An atmospheric pressure plasma jet (APPJ) directed onto a liquid solution containing a metallic salt will promote reduction of metallic ions in solution, generating metallic clusters that nucleate to form nanoparticles. In this article, results from a computational investigation are discussed of a PDSE process in which a radio-frequency APPJ sustained in helium impinges on a silver nitrate solution, resulting in growth of silver nanoparticles. A reaction mechanism was developed and implemented in a global plasma chemistry model to predict nanoparticle growth. To develop the reaction mechanism, density functional theory was used to generate probable silver growth pathways up to Ag 9 . Neutral clusters larger than Ag 9 were classified as nanoparticles. Kinetic reaction rate coefficients for thermodynamically favorable growth pathways were estimated based on an existing, empirically determined base reaction mechanism for smaller Ag particle interactions. These rates were used in conjunction with diffusion-controlled reaction rate coefficients that were calculated for other Ag species. The role of anions in reduction of Agn ions in forming nanoparticles is also discussed. Oxygen containing impurities or admixtures to the helium, air entrainment into the APPJ, and dissociation of saturated water vapor above the solution can produce additional reactive oxygen species in solution, resulting in the production of anions and O 2 – in particular. For a given molarity, delivering a sufficient fluence of reducing species will produce similar nanoparticle densities and sizes for all applied power levels. In conclusion, comparisons are made to alternate models for nanoparticle formation, including charged nanoparticles and use of direct current plasmas.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Twin density and twin thickness evolution in sputtered Al–Mg alloys

Twinned Al–Mg alloys have been reported. However, the role of Mg solute in facilitating the formation of growth twins remains unclear. By using a precession-assisted crystal orientation mapping technique (PACOM) coupled with transmission electron microscopy (known as ASTAR), we examined the evolution of twin boundaries in Al, Al–1Mg, and Al–2.2Mg (at. %) films. The twinned grain fraction elevates with increasing film thickness until it reaches a peak when the film thickness is 120–160 nm. The Al–Mg alloys exhibited greater twinned grain fractions than pure Al. To investigate the fluctuation of twinned grain fraction, two types of twin boundaries were classified including intergranular and intragranular twins. The initial increase in twin density is attributed to the impingement of twinned grains during island coalescence and the twinned grains are more likely to survive during the grain growth process. Whereas the decrease in twinned grain fraction in thicker films is related to the removal of intragranular twins, and a lack of formation mechanisms of new twins.

Physics↗
DOE OSTI2024

Production of warm ions in electron beam generated E × B plasma

Several recent experiments have demonstrated low-damage processing of 2D materials, such as graphene and single crystal diamond, using electron beam (e-beam) generated plasmas with applied crossed electric and magnetic (E × B) fields. Here, the low damage of these sensitive materials is commonly attributed to the low energy of ions incident to the substrate surface and the ion confinement in E x B fields. In this work, measurements of atom and ion velocity distribution functions in an e-beam E × B plasma at sub-mTorr argon pressures using a laser-induced fluorescence diagnostic revealed the presence of a warm population of ions with temperatures of ~ 1 eV that are sufficient to destroy the ion confinement in E × B fields and drive the ion flux by cross field diffusion in the direction opposite to the applied electric field, toward the plasma-bounded walls or substrate. Thus, it is this nearly ambipolar diffusion process that is responsible for the flux of charged particles impinging on the wall/substrate surface.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Development of a compact magnetic spectrometer for use at the OMEGA Laser Facility and the National Ignition Facility

Measurement of proton spectra is an important diagnostic for a variety of high energy density physics experiments. Current diagnostics are either not designed to capture the spectrum of low-energy protons or are unsuitable for high debris experiments. To bridge the gap, a new CR-39 based compact magnetic spectrometer (MagSpec) has been developed to measure proton spectra in the 1–20 MeV energy range, with a particular focus on the low-energy (1–6 MeV) spectrum, for use in experiments at the OMEGA Laser Facility and the National Ignition Facility (NIF). In the MagSpec diagnostic, protons of different energies are dispersed as they pass through a magnetic field before impinging on a differentially filtered CR-39 surface, resulting in a spatial distribution of CR-39 tracks that corresponds to the energy spectrum. In this paper, we discuss details of the design and implementation of MagSpec on the NIF and OMEGA.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

An implicit solution for Asay foil trajectories generated by separable, sustained-production ejecta source models

We present a simple implicit solution for the time-dependent trajectory of a thin Asay foil ejecta diagnostic for the general case where the impinging ejecta cloud is generated by a source function characterized by an arbitrary (sustained) time dependence and a time-independent (stationary) particle velocity distribution. In the limit that the source function time dependence becomes a delta function, this solution—which is amenable to rapid numerical calculations of arbitrary accuracy—exactly recovers a previously published solution for the special case of instantaneous ejecta production. We also derive simple expressions for the free-surface arrival (catch-up) time as well as the true ejecta areal mass accumulation on the accelerating foil and place bounds on the level of error incurred when applying instant-production mass solutions to a sustained-production trajectory. We demonstrate these solutions with example calculations for hypothetical source functions spanning a wide range of ejecta production durations, velocity distributions, and temporal behaviors. These calculations demonstrate how the foil trajectory is often insensitive to the temporal dependence of the source function, instead being dominated by the velocity distribution. We quantify this insensitivity using a “compatibility score” metric. Under certain conditions, one may capitalize upon this insensitivity to obtain a good approximation of the second integral of the velocity distribution from the observed foil trajectory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Rarefied xenon flow in orificed hollow cathodes

A parametric study is conducted to quantify the effect of the keeper electrode geometry on the xenon neutral flow quantities within orificed hollow cathodes, prior to ignition. The keeper impinges directly on the flow out of the cathode orifice and its geometry influences the product between the pressure in the orifice–keeper region and the cathode-to-keeper distance. A representative cathode is simulated using the Direct Simulation Monte Carlo method. The numerical model is first validated with computational results from the literature. A parametric study is then conducted. Parameters include the cathode pressure–diameter in the range of 1–5 Torr cm and the following geometric ratios (and ranges): cathode orifice-to-inner radii (0.1–0.7), keeper orifice-to-cathode orifice radii (1–5), and keeper distance-to-cathode-orifice diameter (0.5–10). It is found that, if both keeper and cathode have identical orifice radii, the flow remains subsonic in the orifice-to-keeper region. In most cases, however, the flow becomes underexpanded and supersonic, and the static pressure within the orifice-to-keeper region is, on average, 4% that of the upstream pressure value. The orifice–keeper region pressure increases with either a decrease in the keeper orifice diameter or an increase in the distance between cathode and keeper, in agreement with literature data. Both trends are explained through conservation laws. A statistical study of numerical results reveals that the ratio of ignition-to-nominal mass flow rates has a most probable value of 50, which suggests that heaterless cathode ignition at a minimum DC voltage may be achieved by increasing the input mass flow rate by a factor of 50.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Design and Analysis of an Advanced Three-Point Bend Test Approach for Miniature Irradiated Disk Specimens

Candidate tungsten armor materials in a magnetic confinement fusion device must be able to withstand thermal variation that leads to internal stresses caused by the impinging heat load. In addition, the thermomechanical properties of these materials are degraded by irradiation-induced defect accumulation. Fission reactor–based irradiation data are used to predict the fusion neutron damage and property change. This study examines the motivation and design of a custom-designed three-point bend test for neutron-irradiated disk specimens that are 3 mm in diameter to be able to define the flexural strength of advanced tungsten materials, alloys, and composites—and to the extent that embrittlement occurs after neutron irradiation. The theory provided shows a calculation for the flexural deflection and shear deflection due to the small-geometry constraints. A finite element deformation analysis is performed to evaluate the mechanical stress field of disk bend specimens. The stress values above 80% of the maximum stress are concentrated in 2.4 mm of the 3.0-mm length of the centerline across the tungsten disk diameter. A bend test fixture has been designed and fabricated to enable testing of these specimens with precisely engineered tolerance and minimal machine compliance. This fixture will be able to be placed inside a universal testing frame at elevated temperatures for the mechanical property evaluation of future neutron-irradiated disk specimens.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Analysis of Dynamic Behavior of the Target and the Deuterium-Tritium Ice in Magnetic-Field Assisted Implosions

The use of strong magnetic fields to augment the output energy of inertial confinement fusion experiments at the National Ignition Facility is of high interest. It offers the potential of reducing electron thermal conduction and increasing hot-spot alpha heating with little to no change in hohlraum behavior. In these magnetically assisted ignition experiments, the ultimate goal is to add a B-field in the form of a pulse ranging from 25 to 60 T to a high-performing hohlraum implosion several microseconds before impingement of the laser beams. This requires eliminating metallic components in the target and replacing them with electrically nonconducting materials. However, the strong eddy currents generated by the rapidly increasing high B-field, which were calculated to be as high as 2000 K, can heat the hohlraum. Here in this paper, we examine the transient effects of this rapid temperature change on the behavior of the target as well as the fuel layer composed typically of deuterium and tritium. Using simulations and calculations for limiting case scenarios, we find that the effect of the heating is not restrictive toward the performance of the target or the quality of the deuterium-tritium ice.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Integral Nuclear Data and Benchmarking Needs for Fusion Energy Systems

Fusion energy systems are currently being designed and optimized using radiation transport codes. To deal with the unique environment inside a fusion-based system, many of these designs incorporate novel materials able to withstand the high radiation fields, ensure adequate cooling and thermal protection, and produce tritium. Validation plays a vital role in building trust in the predictive power of these models and computational methods. Validation of a code consists of modeling documented real-world experiments and comparing the code-predicted response to the measured response. Adequate validation requires measured responses from real-world experiments, also known as integral data, that mimic the system being designed, including materials, impinging radiation, and temperature, among other variables. The most trusted integral data are experimental responses that have been through a rigorous benchmarking process that develops a recommended computational model and evaluates all experimental uncertainties. Finally, there are a few research groups around the world that have been producing integral data for fusion applications, but a substantial investment is needed to address the unique validation needs of the fusion community.

Fusion↗
DOE OSTI2021

Absolute ion density measurements in the afterglow of a radiofrequency atmospheric pressure plasma jet

Atmospheric pressure plasma jets (APPJs) can enable large fluxes of reactive species impinging on treated substrates in many applications. The relative importance of ionic versus neutral reactive species fluxes produced by APPJs is to date mainly deduced from modeling due to the challenge of measuring absolute ion densities/fluxes. In this study, a calibration was performed to determine the absolute densities of cold ions in the afterglow region of an RF-driven atmospheric pressure plasma jet (APPJ) in Ar + 1% O 2 by molecular beam mass spectrometry. Here, positive ion densities were found to be about 3-4 orders of magnitudes lower than the densities of the dominant reactive neutral species (O, O 2 (a 1 Δ g ), O 3 ) in the afterglow region of the APPJ at a distance larger than 2 mm from the visible plasma plume. The results illustrate the ability of molecular beam mass spectrometry to enable a quantitative evaluation of both neutral and ionic species in many applications using APPJs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Atomic oxygen density determination in the effluent of the COST reference source using in situ effective lifetime measurements in the presence of a liquid interface

Abstract Spatially resolved, absolute densities of atomic oxygen are measured for several helium-based admixtures in the effluent of the COST Reference Microplasma Jet using two-photon absorption laser induced fluorescence (TALIF). Admixtures investigated include a helium-only admixture, four helium/oxygen admixtures (0.1 % , 0.5 % , 0.6 % , and 1.0 % oxygen), and a helium/water admixture (2500 ppm water), chosen to coincide with previously published characterizations of plasma-treated liquid. The atomic oxygen TALIF signal is calibrated for density using the noble gas xenon, which possesses a very similar two-photon excitation and fluorescence scheme. Measurements are conducted for the jet operating in ambient air with both an open effluent and a liquid surface present, allowing for comparison with liquid phase measurements conducted under similar conditions. The presence of a water surface does not appear to alter the background chemistry in the effluent but reduces O densities close to the liquid interface when compared to a similar distance from the nozzle in an open effluent case. This may be the result of a reduction in flow velocity caused by the liquid obstructing the gas flow. Additionally, measurements near the liquid surface revealed a region of atomic oxygen well outside of where the core of the effluent impinges on the liquid. This is likely relevant for applications as it considerably expands the surface area subject to O absorption. Critically, in situ measurements of the effective lifetimes of the laser-excited 3p 3 P J state of atomic oxygen were recorded in the effluent by employing a picosecond (ps) laser and a nanosecond (ns) ICCD. By experimentally determining the contribution from collisional quenching via the in situ effective lifetime measurements, significant improvements in the accuracy of the atomic oxygen density calibration were made, with differences of approximately 30 % from existing methods of estimating quenching rates at atmospheric pressure. Finally, spatially resolved atomic oxygen densities allow for an investigation of O formation and extinction pathways in the effluent and a comparison between admixtures.

laser diagnostics↗
DOE OSTI2025

Two-dimensional mapping of absolute OH densities in an atmospheric pressure plasma effluent via planar laser-induced fluorescence: effects of He/H 2 O and He/O 2 mixtures in N 2 and air, with and without solid targets

Planar laser-induced fluorescence (LIF) was employed to measure the absolute density of hydroxyl radicals (OH) in the effluent of the COST Reference Microplasma Jet for two feed gas mixtures: He/H 2 O and He/O 2 . Experiments were conducted with the effluent propagating into air and N 2 environments. For the He/H 2 O case, measurements were also performed with the effluent impinging on a solid target at varying distances from the jet nozzle. Calibration of the OH-LIF signal from the COST-Jet was achieved by comparing it to a reference signal generated by the photofragmentation of H 2 O 2 . Results demonstrated that OH densities were sustained longer when the effluent propagates in a nitrogen environment compared to air, particularly with water added to the feed gas. The broader OH distribution in N 2 suggests slower consumption due to the absence of oxygen, which accelerates OH depletion in air via reactions involving O 2 and HO 2 . Even when water was not added to the feed, as in the He/O 2 case, appreciable OH densities were observed, due to gas impurities and reactive species interactions with atmospheric humidity, forming reaction fronts that delineate the gas flow. Two-dimensional fluid dynamics simulations elucidated the influence of atmospheric gas entrainment and solid targets on the OH distribution. Experimental trends were further compared with a zero-dimensional chemistry model to explore OH production and consumption mechanisms in air and nitrogen environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effect of Helium Flux on Near-surface Helium Accumulation in Plasma-Exposed Tungsten

We report results of object kinetic Monte Carlo (OKMC) simulations aimed at understanding the effect of helium flux on the near-surface helium accumulation in plasma-facing tungsten, which is initially defect-free and has a W(100) surface orientation. These OKMC simulations are performed at 933 K for fluxes ranging from 10 22 to 4 × 10 25 He/m 2 s, with 100 eV helium atoms impinging on a W(100) surface up to a maximum fluence of 4×10 19 He/m 2 . In the near-surface region, helium clusters interact elastically with the free surface. The interaction is attractive and results in the drift of mobile helium clusters towards the surface as well as increased trap mutation rates. The associated kinetics and energetics of the above-mentioned processes obtained from molecular dynamics simulations are also considered. The OKMC simulations indicate that as the flux decreases, the retention of implanted helium decreases, and its depth distribution shifts to deeper below the surface in initially pristine tungsten. Furthermore, the fraction of retained helium diffusing into the bulk increases as well, so much so that for 10 22 He/m 2 s, almost all of the retained helium diffused into the bulk with minimal/negligible near-surface helium accumulation. At a given flux, with increasing fluence, the fraction of retained helium initially decreases and then starts to increase after reaching a minimum. The occurrence of the retention minimum shifts to higher fluences as the flux decreases. Although the near-surface helium accumulation spreads deeper into the material with decreasing flux and increasing fluence, the spread appears to saturate at depths between 80 and 100 nm. Finally, we present a detailed analysis of the influence of helium flux on the size and depth distribution of total helium and helium bubbles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗