Search NASA⌕ Search

SEARCH · Search NASA

Results for “collisional processes”

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 37 records · Page 2

Boundary of the Distribution of Solar Wind Proton Beta versus Temperature Anisotropy

Abstract The frequency distribution of solar wind protons, measured in the vicinity of Earth’s orbit, is customarily plotted in ( β ∥ , T ⊥ / T ∥ ) phase space. Here, T ⊥ / T ∥ is the ratio of perpendicular and parallel temperatures, and β ∥ = 8 π nT ∥ / B 2 is the ratio of parallel thermal energy to background magnetic field energy, the so-called “parallel beta,” with ⊥ and ∥ denoting directions with respect to the ambient magnetic field. Such a frequency distribution, plotted as a two-dimensional histogram, forms a peculiar rhombic shape defined with an outer boundary in the said phase space. Past studies reveal that the threshold conditions for temperature anisotropy–driven plasma instability partially account for the boundary on the high- β ∥ side. The low- β ∥ side remains largely unexplained despite some efforts. Work by Vafin et al. recently showed that certain contours of collisional relaxation frequency, ν pp , when parameterized by T ⊥ / T ∥ and β ∥ , could match the overall shape of the left-hand boundary, thus suggesting that the collisional relaxation process might be closely related to the formation of the left-hand boundary. The present paper extends the analysis by Vafin et al. and carries out the dynamical computation of the collisional relaxation process for an ensemble of initial proton states with varying degrees of anisotropic temperatures. The final states of the relaxed protons are shown to closely match the observed boundary to the left of the ( β ∥ , T ⊥ / T ∥ ) phase space. When coupled with a similar set of calculations for the ensemble in the collective instability regime, it is found that the combined collisional/collective effects provide the baseline explanation for the observation.

Astronomy & Astrophysics↗

Breakups are complicated: an efficient representation of collisional breakup in the superdroplet method

Abstract. A key constraint of particle-based methods for modeling cloud microphysics is the conservation of total particle number, which is required for computational tractability. The process of collisional breakup poses a particular challenge to this framework, as breakup events often produce many droplet fragments of varying sizes, which would require creating new particles in the system. This work introduces a representation of collisional breakup in the so-called “superdroplet” method which conserves the total number of superdroplets in the system. This representation extends an existing stochastic collisional-coalescence scheme and samples from a fragment size distribution in an additional Monte Carlo step. This method is demonstrated in a set of idealized box model and single-column warm-rain simulations. We further discuss the effects of the breakup dynamic and fragment size distribution on the particle size distribution, hydrometeor population, and microphysical process rates. Box model experiments serve to characterize the impacts of properties such as coalescence efficiency and fragmentation function on the relative roles of collisional breakup and coalescence. The results demonstrate that this representation of collisional breakup can produce a stationary particle size distribution, in which breakup and coalescence rates are approximately equal, and that it recovers expected behavior such as a reduction in precipitate-sized particles in the column model. The breakup algorithm presented here contributes to an open-source pythonic implementation of the superdroplet method, PySDM, which will facilitate future research using particle-based microphysics.

de Jong, Emily (ORCID:0000000253104554)↗

Collisionless Larmor coupling and blob formation in a laser-plasma expanding into a magnetized ambient plasma

Collisionless Larmor coupling is a fundamental process in space and astrophysical plasmas that enables momentum transfer between an expanding plasma and a magnetized ambient medium. Here, in this paper, we report on the laboratory experimental study of Larmor coupling leading to the formation of a plasma blob associated with a laser-driven, super-Alfvénic plasma flow on the Large Plasma Device at the University of California, Los Angeles. The high-repetition rate enables systematic spatial and temporal scans of the plasma evolution using Doppler spectroscopy, as well as measurements of the magnetic field, electrostatic field, and self-emission of both debris and ambient ions using filtered imaging. We observe the self-focusing of the laser-produced plasma and the formation of a secondary diamagnetic cavity associated with a blob composed of background ions. Doppler spectroscopy reveals the transverse velocity distribution of the background ions, providing direct evidence of ion energization via Larmor coupling. The systematic spatial and temporal scans enabled by the high-repetition-rate experiment allow for a detailed characterization of the ion dynamics. These experimental observations are supported by numerical simulations that provide more insight into the kinetic-scale physics associated with blob formation as well as the role of the ambient plasma density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Femtosecond electronic structure response to high intensity XFEL pulses probed by iron X-ray emission spectroscopy

We report the time-resolved femtosecond evolution of the K-shell X-ray emission spectra of iron during high intensity illumination of X-rays in a micron-sized focused hard X-ray free electron laser (XFEL) beam. Detailed pulse length dependent measurements revealed that rapid spectral energy shift and broadening started within the first 10 fs of the X-ray illumination at intensity levels between 10 17 and 10 18 W cm -2 . We attribute these spectral changes to the rapid evolution of high-density photoelectron mediated secondary collisional ionization processes upon the absorption of the incident XFEL radiation. These fast electronic processes, occurring at timescales well within the typical XFEL pulse durations (i.e., tens of fs), set the boundary conditions of the pulse intensity and sample parameters where the widely-accepted ‘probe-before-destroy’ measurement strategy can be adopted for electronic-structure related XFEL experiments.

99 GENERAL AND MISCELLANEOUS↗

Driven electrostatic phase space vortices in a 1D weakly dissipative Vlasov–Poisson system

In this paper, the effect of collisions on driven electrostatic phase space vortices is analyzed by means of Eulerian simulation for two different collision models. It was demonstrated recently that in the absence of collisions, at late times, steady state phase space vortices manifest to form a plateau in the resonant region of the particle velocity distribution function, due to trapping of particles supporting multiextrema giant phase space vortices (PSVs). In the presence of collisions, over long time, this multiextrema plateau are found to smooth out, since collisions drive the velocity distribution toward Maxwellian, irrespective of how weak the collisions are as long as they are non-zero. In these conditions, kinetic processes and collisionality are found to be in competition, and the evolution of the plasma is found, therefore, to be a result of nontrivial combination of these two effects. An attempt has been made by means of numerical simulations to study the effect of weak collisionality on the electrostatic driven phase space vortices with two types of collision operators: (1) Bhatnagar–Gross–Krook (Krook) collision operator, where the colliding particles can be treated as isolated pairs and, (2) Fokker–Planck (FP) type collision operator (Zakharov–Karpman) in one dimension, where many weak collisions lead to particle diffusion in velocity space. It is shown that depending on the collision model used, the nature of smoothing in velocity space of giant PSVs results in qualitatively very different phase space structures. However, irrespective of the collision model used, excess density fractions over 10% are retained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multi-fluid and kinetic models of partially ionized magnetic reconnection

Magnetic reconnection in partially ionized plasmas is a ubiquitous and important phenomenon in both laboratory and astrophysical systems. Here, simulations of partially ionized magnetic reconnection with well-matched initial conditions are performed using both multi-fluid and fully-kinetic approaches. Despite similar initial conditions, the time-dependent evolution differs between the two models. In multi-fluid models, the reconnection rate locally obeys either a decoupled Sweet–Parker scaling, where neutrals are unimportant, or a fully coupled Sweet–Parker scaling, where neutrals and ions are strongly coupled, depending on the resistivity. In contrast, kinetic models show a faster reconnection rate that is proportional to the fully-coupled, bulk Alfvén speed, $v^*_A$. In this work, these differences are interpreted as the result of operating in different collisional regimes. Multi-fluid simulations are found to maintain $ν_{ni}L/v^*_A$ ≳1, where $ν_{ni}$ is the neutral–ion collision frequency and L is the time-dependent current sheet half-length. This strongly couples neutrals to the reconnection outflow, while kinetic simulations evolve to allow $ν_{ni}L/v^*_A$ <1, decoupling neutrals from the reconnection outflow. Differences in the way reconnection is triggered may explain these discrepancies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computing the Paschen curve for argon with speed-limited particle-in-cell simulation

Upon inclusion of collisions, the speed-limited particle-in-cell (SLPIC) simulation method successfully computed the Paschen curve for argon. The 1D3V simulations modeled an electron cascade across an argon-filled capacitor, including electron-neutral ionization, electron-neutral elastic collisions, electron-neutral excitation, and ion-induced secondary electron emission. In electrical breakdown, the timescale difference between ion and electron motion makes traditional PIC methods computationally slow. To decrease this timescale difference and speed up computation, here we used SLPIC, a time-domain algorithm that limits the speed of the fastest electrons in the simulation. The SLPIC algorithm facilitates a straightforward, fully kinetic treatment of dynamics and collisions. SLPIC was as accurate as PIC, but ran up to 200 times faster. SLPIC accurately computed the Paschen curve for argon over three orders of magnitude in pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Benchmarking NIMROD continuum kinetic formulations through the steady-state poloidal flow

In this work, continuum kinetic formulations are employed as a mechanism to include closure physics in an extended magnetohydrodynamics model. Two continuum kinetic approaches have been implemented in the plasma fluid code NIMROD [Sovinec et al., “Nonlinear magnetohydrodynamics with high-order finite elements,” J. Comput. Phys. 195, 355 (2004)] including a Chapman–Enskog-like (CEL) formulation and a more conventional df approach. Ion kinetic closure schemes are employed to describe the neoclassical flow properties in axisymmetric toroidal geometry. In particular, predictions for steady-state values of poloidal flow profiles in tokamak geometry are provided using both the df formulation and two different solution techniques for the CEL approach. These results are benchmarked against analytic theory predictions as well as results from the drift kinetic code DK4D. The continuum kinetic formulations employed here show agreement with both the analytic theory and DK4D results, and offer a novel velocity space representation involving higher-order finite elements in pitch angle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measurements of the imploding plasma sheath in triple-nozzle gas-puff z pinches

Gas-puff z-pinch implosions are characterized by the formation of a dense annular plasma shell, the sheath, that is driven to the axis by magnetic forces and therefore subject to the magneto-Rayleigh–Taylor instability. Here, the conditions within these sheaths are measured on the 1-MA COBRA generator at Cornell University for various gas species and initial fill densities. The gas-puff loads are initialized by a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet. Thomson scattering and laser interferometry provide spatially resolved flow, temperature, and electron density profiles midway through the implosion, while extreme ultraviolet pinhole cameras record the evolution of the plasma column and photoconducting diodes measure x-ray emission. Analysis of the scattering spectra includes a means of discriminating between thermal and non-thermal broadening to test for the presence of hydrodynamic turbulence. Two types of sheath profiles are observed, those with sharp discontinuities at the leading edge and those with smooth gradients. In both cases, non-thermal broadening is generally peaked at the front of the sheath and exhibits a characteristic decay length that roughly scales with the sheath ion mean free path. Here, we demonstrate that this non-thermal broadening term is inconsistent with laminar velocity gradients and is more consistent with dissipative turbulence driven by unstable plasma waves in a collisionless shock. The resulting differences in sheath profile are then set by the sheath ion collisionality in a manner consistent with recent 1D kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of a central barrier potential on asymmetry-induced transport in a coaxial Malmberg–Penning trap

Experiments are presented measuring the radial particle flux produced when a central barrier potential is applied to the central ring of a coaxial Malmberg-Penning trap at the same time as the voltages producing the main asymmetry potential $\phi_1(r)\cos{(kz)}\cos{(l\theta-\omega t)}$. When a negative DC voltage is applied to the entire central ring, the peak flux produced by the main asymmetry is reduced by a factor $e^{(V_{eff}/1.2)}$ where $V_{eff}$ is the effective applied potential in volts. This barrier potential also increases the frequency of this peak flux. When asymmetric $\pm$ voltages are applied to the two azimuthally-divided halves of the central ring, DC voltages increase the radial flux while AC voltages decrease it. It is shown that these asymmetric barrier voltages produce their own transport and that the resulting modification of the plasma may be a factor in changing the observed flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bounce-averaged theory in arbitrary multi-well plasmas: solution domains and the graph structure of their connections

Bounce-averaged theories provide a framework for simulating relatively slow processes, such as collisional transport and quasilinear diffusion, by averaging these processes over the fast periodic motions of a particle on a closed orbit. This procedure dramatically increases the characteristic time scale and reduces the dimensionality of the modelled system. The natural coordinates for such calculations are the constants of motion (COM) of the fast particle motion, which by definition do not change during an orbit. However, for sufficiently complicated fields – particularly in the presence of local maxima of the electric potential and magnetic field – the COM are not sufficient to specify the particle trajectory. In such cases, multiple domains in COM space must be used to solve the problem, with boundary conditions enforced between the domains to ensure continuity and particle conservation. Previously, these domains have been imposed by hand, or by recognising local maxima in the fields, limiting the flexibility of bounce-averaged simulations. Here, we present a general set of conditions for identifying consistent domains and the boundary condition connections between the domains, allowing the application of bounce-averaged theories in arbitrarily complicated and dynamically evolving electromagnetic field geometries. We also show how the connections between the domains can be represented by a directed graph, which can help to succinctly represent the trajectory bifurcation structure.

fusion plasma↗

Synthesis of observed and simulated rain microphysics to inform a new Bayesian statistical framework for microphysical parameterization in climate models (Final Scientific Report)

This project aimed to investigate a new approach to bulk microphysics schemes. We identified the need to move beyond the fixed structural assumptions and approximation of existing schemes. For example, most schemes assume some functional form for the rain drop size distribution (e.g. a gamma or exponential distribution). Most schemes then evolve some number of statistical moments of that distribution via various processes, such as evaporation, sedimentation, collision-coalescence, and collisional breakup. These microphysical processes, in turn, are typically some fixed functional form derived from either some other (more detailed) model, or using some single-particle rates that are then integrated over the size distribution. While these process rate functions may have some free parameters to adjust, in other cases doing so is impossible (e.g. when the process rates are analytical or piecewise solutions to some target function). Our goal was to build and test a scheme that assumed no size distribution form, and used a series of power laws as the basis for the process rates. Moments of the size distribution would be predicted, but no underlying size distribution would be specified. Any number or choice of size distribution moments could be used, and any number of power laws could be employed to model the process rates. Thus, our approach is flexible, and can seamlessly scale across levels of complexity, as demanded by the data. Bayesian inference then would provide the formalism to estimate the model parameters and structure, allowing for robust uncertainty quantification. We proposed testing this framework in idealized simulations using bin microphysical schemes as a data source. We also proposed using real data to inform our microphysics scheme, and also that we would integrate our scheme into WRF.

58 GEOSCIENCES↗

Prompt and non-prompt J/$ψ$ production at midrapidity in Pb–Pb collisions at $\sqrt{s_{\textrm{NN}}}$ = 5.02 TeV

The transverse momentum (p T ) and centrality dependence of the nuclear modif cation factor R AA of prompt and non-prompt J/ψ, the latter originating from the weak decays of beauty hadrons, have been measured by the ALICE collaboration in Pb–Pb collisions at $\sqrt{s_{\textrm{NN}}}$ = 5.02 TeV. The measurements are carried out through the e + e - decay channel at midrapidity (|y| < 0.9) in the transverse momentum region 1.5 < p T < 10 GeV/c. Both prompt and non-prompt J/ψ measurements indicate a significant suppression for p T > 5 GeV/c, which becomes stronger with increasing collision centrality. The results are consistent with similar LHC measurements in the overlapping p T intervals, and cover the kinematic region down to p T = 1.5 GeV/c at midrapidity, not accessible by other LHC experiments. The suppression of prompt J/ψ in central and semicentral collisions exhibits a decreasing trend towards lower transverse momentum, described within uncertainties by models implementing J/ψ production from recombination of $c$ and $\overline{c}$ quarks produced independently in different partonic scatterings. At high transverse momentum, transport models including quarkonium dissociation are able to describe the suppression for prompt J/ψ. For non-prompt J/ψ, the suppression predicted by models including both collisional and radiative processes for the computation of the beauty-quark energy loss inside the quark-gluon plasma is consistent with measurements within uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measurement of beauty production via non-prompt D 0 mesons in Pb-Pb collisions at $ \sqrt{{\textrm{s}}_{\textrm{NN}}} $ = 5.02 TeV

The production of non-prompt D 0 mesons from beauty-hadron decays was measured at midrapidity (|y| < 0.5) in Pb-Pb collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_{NN}}$= 5.02 TeV with the ALICE experiment at the LHC. Their nuclear modification factor (R AA ), measured for the first time down to p T = 1 GeV/c in the 0–10% and 30–50% centrality classes, indicates a significant suppression, up to a factor of about three, for p T > 5 GeV/c in the 0–10% central Pb-Pb collisions. The data are described by models that include both collisional and radiative processes in the calculation of beauty-quark energy loss in the quark-gluon plasma, and quark recombination in addition to fragmentation as a hadronisation mechanism. The ratio of the non-prompt to prompt D 0 -meson R AA is larger than unity for p T > 4 GeV/c in the 0–10% central Pb-Pb collisions, as predicted by models in which beauty quarks lose less energy than charm quarks in the quark-gluon plasma because of their larger mass.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dynamics of deuterium retention and desorption from plasma-facing materials in fusion reactor-relevant conditions

Hydrogen isotopes retention and desorption during and after discharges in fusion devices are still not well understood due to the complex device conditions and limitations of in-situ diagnostics and measurements. Here, we simulated well-diagnosed recent experiments at the DIII-D facility to benchmark our ITMC-DYN integrated package of modeling deuterium diffusion, retention, and desorption during and after D discharge irradiation. Modeling results were compared with detail experimental data of D desorption fluxes for various irradiation conditions. We predicted the temporal evolution of free and trapped D distribution in tungsten (W) plasma-facing material (PFM). Effects of key parameters namely diffusion coefficient, recombination rate, trapping energies against different defect types, were examined in these simulations. Existing experimental data of these parameters in literature varies significantly which makes it harder to identify key mechanisms and physics responsible for hydrogen isotope retention and desorption. The purpose of this work is to accurately simulate recent well-diagnosed reactor experiments given the uncertainties in such parameters and identify mechanisms responsible for the retention and desorption. We implemented the best identified diffusion, recombination, and trapping parameters in ITMC-DYN package that integrate both various collisional and thermal processes. We predicted, for example, that sample cooling between discharges in DIII-D operations can significantly affect the spatial distribution of trapped D in W under reactor irradiation conditions. Correct prediction of desorption spectra from samples irradiated during 10 DIII-D discharges showed that up to 35% of D can be retained in high binding energy defects such as vacancy clusters or voids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Antihydrogen chemistry

A survey of antimatter reactions is presented, including the formation of the antihydrogen atom and anionic, cationic, and molecular species by collisional and radiative processes. Our approach is rooted in the detailed knowledge available for many matter counterpart (hydrogenic) reactions, due to their importance in controlling early Universe chemistry. We point out that the availability of trapped antihydrogen at densities similar to those pertaining to the epoch of hydrogen chemistry will soon be available. In addition, using modern atomic physics techniques, it should be feasible to control antimatter in the laboratory to facilitate antihydrogen chemistry. Our purpose is to summarize what is known from hydrogen chemistry that is of relevance for antimatter and to indicate, based on possible reaction rates, which processes may be fruitful to pursue to create new antimatter entities as probes of fundamental symmetries. We include antihydrogen, positrons, and antiprotons in our discussion and additionally the electron due to its propensity to form positronium and perhaps to participate in certain reactions. We attempt to indicate whether further theoretical/computational work is necessary to add to the assessment of reaction rates, and we discount processes where the projected rates are too low to be of interest, given foreseeable experimental capabilities.

74 ATOMIC AND MOLECULAR PHYSICS↗

Hot Electron Scaling and Energy Coupling In Nonlinear Laser Plasma Interactions

Energy conversion from high power lasers to energetic hot electrons is a fundamental phenomenon in laser-plasma interactions. At laser intensities above a few times 10 14 W/cm 2 , nonlinear laser-plasma interactions lead to the generation of hot electrons with a broad energy distribution at typical temperatures of 10’s – 100’s kilo-electron-volt (keV). The fraction of laser energy transferred to hot electrons, their energy distribution and energy coupling efficiency to the targets are strongly affected by the laser intensity, wavelength and the plasma condition. Most investigations to date that are related to the fusion of high-energy-density laboratory plasmas (HEDLP) are carried out with relatively moderate laser intensity and lower laser wavelength to minimize hot electron generation, as electron preheating can degrade the fuel assembly. An alternative fusion scheme known as shock ignition (SI) has been proposed, requiring a high intensity spike pulse at the end of a fuel assembly phase to generate a strong shock to ignite the fuel. An initial shock pressure of 0.3 Gbar is required for ignition-scale targets using the SI scheme, requiring a relatively high on-target laser intensity of 5-10X10 15 W/cm 2 for the spike pulse. At such intensities, copious hot electrons can be produced by laser-plasma-instabilities (LPI), e.g., Stimulated Raman Scattering (SRS) and Two-Plasmon Decay (TPD). While very energetic electrons (>200 keV) could preheat the fuel, electrons with moderate energies may actually be beneficial for the SI scheme. Electrons in the 50-100 keV range can be stopped in the compressed high-density outer ablator layer leading to an increase in the ignitor shock strength and therefore enhancing overall efficiency in the SI scheme. Considering this importance of total energy and temperature of hot electrons, it is necessary to understand the LPI in the SI high intensity regime. Characterizing the hot electron source and energy coupling is key to investigating the viability of the SI concept. The overall goal of the HEDLP award was to perform a detailed study of the scaling of hot electron generation and energy coupling on high power laser systems, with intensity, wavelength, and target ablator material relevant to the nonlinear laser-plasma interaction at SI relevant conditions. This was a collaborative project between the Center for Energy Research at the University of California, San Diego (UCSD) and General Atomics. The experiment was performed at the Omega Laser Facility at the Laboratory for Laser Energetics (LLE) in the University of Rochester. Specific objectives of the work were to: understand the basic physics of laser-plasma interactions and instabilities at high laser intensities (up to 10 17 W/cm 2 ) and the effects of low-Z ablator materials and laser wavelengths (IR pulse at 1.054 µm versus UV pulse at 0.351 µm) and intensity scaling, fully characterize the resultant hot electron energy distribution, validate hot electron energy deposition via collisional and kinetic processes in the target and the induced shocks in support of electron-assisted SI concept.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗