SEARCH · Search NASA

Results for “MAGNETOHYDRODYNAMIC GENERATOR”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS2021

Magnetic Reconnection in a Sheared Magnetic Flux Tube: Slippage Versus Tearing

The process of magnetic reconnection in a flux tube can occur in a time-stationary fashion as slippage reconnection or in a time-dependent manner based on, for example, the tearing instability. However, it is not well known under which conditions a system can sustain slippage reconnection. Likewise, it is unclear whether systems can exhibit slippage reconnection and time-dependent reconnection simultaneously. In order to investigate these questions, we employ a set of 3D magnetohydrodynamic simulation. Using these simulations, we model a twisted flux tube with a spatially localized resistive region in the center of the simulation box, and an applied driver on one end plane of the simulation box. As a result, the Poynting flux injected at the boundary propagates to the resistive region and dissipates there. If the driver is nearly circular, we find a predominance of slippage reconnection, which decouples the applied driver entirely from one half of the flux tube. Second, a tearing-like instability occurs with slippage reconnection in the resistive region when the resistivity is strong enough and the cross section of the flux tube, as defined by the shape of the applied velocity perturbation, is elliptical enough. In this case we find the surprising feature that magnetic field perturbation generated by the tearing-like instability radiates away from the resistive region to influence the entire flux tube shape.

Hidetaka Kuniyoshi↗

DOE OSTI2024

Design and modeling of indirectly driven magnetized implosions on the NIF

The use of magnetic fields to improve the performance of hohlraum-driven implosions on the National Ignition Facility (NIF) is discussed. The focus is on magnetically insulated inertial confinement fusion, where the primary field effect is to reduce electron-thermal and alpha-particle loss from the compressed hotspot (magnetic pressure is of secondary importance). We summarize the requirements to achieve this state. The design of recent NIF magnetized hohlraum experiments is presented. These are close to earlier shots in the three-shock, high-adiabat (BigFoot) campaign, subject to the constraints that magnetized NIF targets must be fielded at room-temperature, and use ≲1 MJ of laser energy to avoid the risk of optics damage from stimulated Brillouin scattering. We present results from the original magnetized hohlraum platform, as well as a later variant that gives a higher hotspot temperature. In both platforms, imposed fields (at the capsule center) of up to 28 T increase the fusion yield and hotspot temperature. Integrated radiation-magneto-hydrodynamic modeling with the Lasnex code of these shots is shown, where laser power multipliers and a saturation clamp on cross-beam energy transfer are developed to match the time of peak capsule emission and the P2 Legendre moment of the hotspot x-ray image. The resulting fusion yield and ion temperature agree decently with the measured relative effects of the field, although the absolute simulated yields are higher than the data by 2.0−2.7×. The tuned parameters and yield discrepancy are comparable for experiments with and without an imposed field, indicating the model adequately captures the field effects. Self-generated and imposed fields are added sequentially to simulations of one BigFoot NIF shot to understand how they alter target dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1976

High-beta turbulence in two-dimensional magnetohydrodynamics

Equations of ideal magnetohydrodynamics are used to study incompressible turbulent flows in a specified geometry where all the field quantities vary with only two spatial dimensions. The procedures adopted are basically those of Kraichnan (1967), in which classical equilibrium ensembles are built around constants of the motion identified from the Fourier-transformed equations of motion. Once the constants of the motion are identified, the statistical formulation of the problem is presented in a phase space whose coordinates are the real and imaginary parts of the Fourier coefficients. Canonical ensembles are constructed in this phase space by classical arguments. The theory developed permits isolation of some qualitatively new gross physical effects which have so far not been calculated. One of the more novel physical effects is the appearance of macroscopic structures involving long-wavelength, self-generated, magnetic fields for a wide range of initial parameters.

Fyfe, D.↗

NASA NTRS1997

The Ambient and Perturbed Solar Wind: From the Sun to 1 AU

The overall objective of the proposed research was to use numerical solutions of the magnetohydrodynamic (MHD) equations along with comparisons of the computed results with observations to study the following topics: (1) ambient solar wind solutions that extend from the solar surface to 1 astronomical unit (AU), contain closed magnetic structures near the Sun, and are consistent with observed values; (2) magnetic and plasma structures in coronal mass ejections (CMES) as they propagate to the interplanetary medium; (3) relation of MHD shocks to CMEs in the interplanetary medium; (4) interaction of MHD shocks with structures (such as other shocks, corotating interaction regions, current sheets) in the interplanetary plasma; and (5) simulations of observed interplanetary structures. A planned close collaboration with data analysts served to make the model more relevant to the data. The outcome of this research program is an improved understanding of the physical processes occurring in solar-generated disturbances in the interplanetary medium between the Sun and 1 AU.

Steinolfson, R. S.↗

NASA NTRS2022

Dayside Transient Phenomena and Their Impact on the Magnetosphere and Ionosphere

Dayside transients, such as hot flow anomalies, foreshock bubbles, magnetosheath jets, flux transfer events, and surface waves, are frequently observed upstream from the bow shock, in the magnetosheath, and at the magnetopause. They play a significant role in the solar wind-magnetosphere-ionosphere coupling. Foreshock transient phenomena, associated with variations in the solar wind dynamic pressure, deform the magnetopause, and in turn generates field-aligned currents (FACs) connected to the auroral ionosphere. Solar wind dynamic pressure variations and transient phenomena at the dayside magnetopause drive magnetospheric ultra low frequency (ULF) waves, which can play an important role in the dynamics of Earth’s radiation belts. These transient phenomena and their geoeffects have been investigated using coordinated in-situ spacecraft observations, spacecraft-borne imagers, ground-based observations, and numerical simulations. Cluster, THEMIS, Geotail, and MMS multi-mission observations allow us to track the motion and time evolution of transient phenomena at different spatial and temporal scales in detail, whereas ground-based experiments can observe the ionospheric projections of transient magnetopause phenomena such as waves on the magnetopause driven by hot flow anomalies or flux transfer events produced by bursty reconnection across their full longitudinal and latitudinal extent. Magnetohydrodynamics (MHD), hybrid, and particle-in-cell (PIC) simulations are powerful tools to simulate the dayside transient phenomena. This paper provides a comprehensive review of the present understanding of dayside transient phenomena at Earth and other planets, their geoeffects, and outstanding questions.

Bow shock↗

NASA NTRS2022

The Imprint of Intermittent Interchange Reconnection on the Solar Wind

The solar wind is known to be highly structured in space and time. Observations from Parker Solar Probe have revealed an abundance of so-called magnetic switchbacks within the near-Sun solar wind. In this Letter, we use a high-resolution, adaptive-mesh, magnetohydrodynamics simulation to explore the disturbances launched into the solar wind by intermittent/bursty interchange reconnection and how they may be related to magnetic switchbacks. We find that repeated ejection of plasmoid flux ropes into the solar wind produces a curtain of propagating and interacting torsional Alfvénic waves. We demonstrate that this curtain forms when plasmoid flux ropes dynamically realign with the radial field as they are ejected from the current layer and that this is a robust effect of the 3D geometry of the interchange reconnection region. Simulated flythroughs of this curtain in the low corona reveal an Alfvénic patch that closely resembles observations of switchback patches, but with relatively small magnetic field deflections. Therefore, we suggest that switchbacks could be the solar wind imprint of intermittent interchange reconnection in the corona, provided an in situ process subsequently amplifies the disturbances to generate the large deflections or reversals of radial field that are typically observed. That is to say, our results indicate that a combination of low-coronal and inner-heliospheric mechanisms may be required to explain switchback observations.

corona↗

DOE OSTI2024

High-Fidelity Arc-Discharge Model for Hydrogen-Plasma-Smelting-Reduction of Iron Ore

Electrification and use of renewable hydrogen is currently a necessity for decarbonizing the iron-and-steel industry. In this regard, hydrogen plasma smelting reduction (HPSR) is a novel pathway that is being explored for reduction of iron ore. HPSR provides several decarbonization merits compared to conventional blast furnaces. Firstly, the use of renewable hydrogen drastically reduces the CO2 emissions compared to the use of coke. Secondly, renewable electricity in the form of a thermal plasma for making reactive hydrogen species (radicals, ions) are more efficient at reducing iron ore compared to neutral H2. Thirdly, a molten product compatible with downstream processes is obtained from the intense heat transfer from the plasma. However, the scale-up of this technology requires fundamental exploration of hydrogen plasma dynamics and its interaction with complex solid material that include phase changing iron-ore and slag. In this work, we present a first principles continuum scale model for thermal plasmas in Ar/H2 gas mixtures typically used for HPSR. The thermal plasma governing equations for mass, momentum and energy with Lorentz force and Joule heating source terms are solved along with electromagnetic equations for electrostatic and magnetic vector potential. Our solver will be based on Pele, a suite of reacting flow solvers designed for advanced scientific computing architectures (Henry De Frahan et al., Proceedings of SIAM Parallel Processing, 13-25, 2024), and will utilize adaptive mesh generation for enhanced resolutions at locations of intense physicochemical interactions. This study will present the impact of Ar to H2 ratios on excited/dissociated hydrogen species concentrations, plasma temperature and conductivity along with the impact of outgassed species (water, metal vapor, O, OH radicals) from ore surface on gas phase chemistry. Furthermore, the heat and species flux to the surface will be quantified as a function of applied voltages in a transferred arc configuration.

hydrogen plasma↗

NASA NTRS2011

Multi-MW Closed Cycle MHD Nuclear Space Power Via Nonequilibrium He/Xe Working Plasma

Prospects for a low specific mass multi-megawatt nuclear space power plant were examined assuming closed cycle coupling of a high-temperature fission reactor with magnetohydrodynamic (MHD) energy conversion and utilization of a nonequilibrium helium/xenon frozen inert plasma (FIP). Critical evaluation of performance attributes and specific mass characteristics was based on a comprehensive systems analysis assuming a reactor operating temperature of 1800 K for a range of subsystem mass properties. Total plant efficiency was expected to be 55.2% including plasma pre-ionization power, and the effects of compressor stage number, regenerator efficiency and radiation cooler temperature on plant efficiency were assessed. Optimal specific mass characteristics were found to be dependent on overall power plant scale with 3 kg/kWe being potentially achievable at a net electrical power output of 1-MWe. This figure drops to less than 2 kg/kWe when power output exceeds 3 MWe. Key technical issues include identification of effective methods for non-equilibrium pre-ionization and achievement of frozen inert plasma conditions within the MHD generator channel. A three-phase research and development strategy is proposed encompassing Phase-I Proof of Principle Experiments, a Phase-II Subscale Power Generation Experiment, and a Phase-III Closed-Loop Prototypical Laboratory Demonstration Test.

Litchford, Ron J.↗

DOE OSTI2024

Density and Magnetic Field Asymmetric Kelvin‐Helmholtz Instability

Abstract The Kelvin‐Helmholtz (KH) instability can transport mass, momentum, magnetic flux, and energy between the magnetosheath and magnetosphere, which plays an important role in the solar‐wind‐magnetosphere coupling process for different planets. Meanwhile, strong density and magnetic field asymmetry are often present between the magnetosheath (MSH) and magnetosphere (MSP), which could affect the transport processes driven by the KH instability. Our magnetohydrodynamics simulation shows that the KH growth rate is insensitive to the density ratio between the MSP and the MSH in the compressible regime, which is different than the prediction from linear incompressible theory. When the interplanetary magnetic field (IMF) is parallel to the planet's magnetic field, the nonlinear KH instability can drive a double mid‐latitude reconnection (DMLR) process. The total double reconnected flux depends on the KH wavelength and the strength of the lower magnetic field. When the IMF is anti‐parallel to the planet's magnetic field, the nonlinear interaction between magnetic reconnection and the KH instability leads to fast reconnection (i.e., close to Petschek reconnection even without including kinetic physics). However, the peak value of the reconnection rate still follows the asymmetric reconnection scaling laws. We also demonstrate that the DMLR process driven by the KH instability mixes the plasma from different regions and consequently generates different types of velocity distribution functions. We show that the counter‐streaming beams can be simply generated via the change of the flux tube connection and do not require parallel electric fields.

Astronomy & Astrophysics↗

NASA NTRS1985

A Leaky Waveguide Model for MHD Wave Driven Winds from Coronal Holes

Magnetohydrodynamic (MHD) waves, driven by the large scale convective motions of the photosphere are suggested as a possible source of additional acceleration for the stellar wind. Most of the turbulent power in a coronal hole is carried by MHD waves with periods of a few hundred seconds or longer. This is evident from direct observations of turbulence in the solar photosphere, as well as in situ observations of turbulence in the solar wind. But waves with periods this long have wavelengths which are typically as large as the transverse scale of the coronal hole flux tube itself. For these waves boundary effects are important and the coronal hole must be treated as a waveguide. The propagation of MHD waves using this waveguide approach is discussed. The simple model presented demonstrates that coronal holes can act as waveguides for MHD waves. For typical solar parameters the waves are compressible and can generate a wave tensile force which tends to cancel at least part of the wave pressure force. This effect tends to decrease the efficiency of MHD wave acceleration.

Davila, J. M.↗

NASA NTRS2006

GRMHD Simulations of Jet Formation with RAISHIN

We have developed a new three dimensional general relativistic magnetohydrodynamic (GRMHD) code, RAISHIN, using a conservative, high-resolution shock capturing scheme. Numerical fluxes are calculated using the Harten, Lax, & van Leer (HLL) approximate Riemann solver scheme. The flux-interpolated, constrained transport scheme is used to maintain a divergence-free magnetic field. We describe code performance on some test problems in both special and general relativity. Our new GRMHD code has proven to be accurate to second order and has successfully passed several numerical test problems including highly relativistic and magnetized tests in both special and general relativity. We have performed several simulations of non-rotating and rotating black hole systems with a geometrically thin accretion disk. The simulations show the formation of jets driven by the Lorentz force and the gas pressure. It appears that the rotating black hole creates an additional faster, and more collimated outflow inside a broader, slower outflow that is also generated by the rotating accretion disk around a non-rotating black hole. The kinematic jet structure could thus be a sensitive function of black hole rotation.

Mizuno, Y.↗

NASA NTRS2005

Warm Molecular Gas Traced with CO J = 7 --> 6 in the Galaxy's Central 2 Parsecs: Dynamical Heating of the Circumnuclear Disks

We present an 11" resolution map of the central 2 pc of the Galaxy in the CO J = 7 --> 6 rotational transition. The CO emission shows rotation about Sgr A* but also evidence for noncircular turbulent motion and a clumpy morphology. We combine our data set with available CO measurements to model the physical conditions in the disk. We find that the molecular gas in the region is both warm and dense, with T approx. 200-300 K and n(sub H2) approx. (5-7) x 10(exp 4) cm(exp -3). The mass of warm molecular gas we measure in the central 2 pc is at least 2000 M(solar), about 20 times the UV-excited atomic gas mass, ruling out a UV heating scenario for the molecular material. We compare the available spectral tracers with theoretical models and conclude that molecular gas is heated with magnetohydrodynamic shocks with v approx. 10-20 km s(exp -1) and B approx. 0.3- 0.5 mG. Using the conditions derived with the CO analysis, we include the other important coolants, neutral oxygen and molecular hydrogen, to estimate the total cooling budget of the molecular material. We derive a mass-to-luminosity ratio of approx. 2-3 M(solar)(L(solar)exp -1), which is consistent with the total power dissipated via turbulent decay in 0.1 pc cells with v(sub rms) approx. 15 kilometers per second. These size and velocity scales are comparable to the observed clumping scale and the velocity dispersion. At this rate, the material near Sgr A* is dissipating its orbital energy on an orbital timescale and cannot last for more than a few orbits. Our conclusions support a scenario in which the features near Sgr A* such as the circumnuclear disk and northern arm are generated by infalling clouds with low specific angular momentum.

galatic nuclei↗

DOE OSTI2025

Numerical investigation of two-dimensional buoyancy-driven eddies in liquid metal magnetohydrodynamic flows in breeding blankets *

Abstract Lead–lithium flows are key features in the design of tokamak breeding blanket concepts such as the dual-coolant lead–lithium (DCLL). Since they flow under magnetic fields, they are affected by magnetohydrodynamic (MHD) effects. The neutron flux originating in the tokamak plasma heats the breeding blanket channels in a non-uniform manner, inducing buoyancy forces in the liquid metal. Buoyancy may become a source of quasi-two-dimensional (Q2D) turbulence, and the appearance of eddies may affect the transport of heat and tritium across the blanket. Blankets characterized by high-speed liquid metal flows (such as DCLL) will need ceramic insulating walls to reduce the MHD-related pressure drop in the channels. In our simulations, we have used the Q2D model proposed by Sommeria and Moreau (SM82) which is especially suitable for modeling electrically insulating channel flows. On top of that, we have modelled buoyancy forces in the momentum equation using the Oberbeck–Boussinesq approximation. In this work, we include a validation of the implemented Q2D model in buoyancy-driven cases and identify a thermohydraulic configuration that promotes the generation of eddies and the accumulation of tritium. We also show the results provided by our post-processing tool based on the bi-dimensional fast Fourier transform for eddy detection and characterization. We complete our investigation by performing an initial assessment of how eddies can accumulate tritium in the breeding blanket. We conclude by discussing the relevance of preventing accumulation of tritium within the blanket and suggest a possible solution.

Suarez, Daniel (ORCID:0000000251740100)↗

NASA NTRS2022

Quantifying Errors in 3D CME Parameters Derived from Synthetic Data Using White-Light Reconstruction Techniques

Current efforts in space weather forecasting of CMEs have been focused on predicting their arrival time and magnetic structure. To make these predictions, methods have been developed to derive the true CME speed, size, position, and mass, among others. Difficulties in determining the input parameters for CME forecasting models arise from the lack of direct measurements of the coronal magnetic fields and uncertainties in estimating the CME 3D geometric and kinematic parameters after eruption. White-light coronagraph images are usually employed by a variety of CME reconstruction techniques that assume more or less complex geometries. This is the first study from our International Space Science Institute (ISSI) team “Understanding Our Capabilities in Observing and Modeling Coronal Mass Ejections”, in which we explore how subjectivity affects the 3D CME parameters that are obtained from the Graduated Cylindrical Shell (GCS) reconstruction technique, which is widely used in CME research. To be able to quantify such uncertainties, the “true” values that are being fitted should be known, which are impossible to derive from observational data. We have designed two different synthetic scenarios where the “true” geometric parameters are known in order to quantify such uncertainties for the first time. We explore this by using two sets of synthetic data: 1) Using the ray-tracing option from the GCS model software itself, and 2) Using 3D magnetohydrodynamic (MHD) simulation data from the Magnetohydrodynamic Algorithm outside a Sphere code. Our experiment includes different viewing configurations using single and multiple viewpoints. CME reconstructions using a single viewpoint had the largest errors and error ranges overall for both synthetic GCS and simulated MHD white-light data. As the number of viewpoints increased from one to two, the errors decreased by approximately 4° in latitude, 22° in longitude, 14° in tilt, and 10° in half-angle. Our results quantitatively show the critical need for at least two viewpoints to be able to reduce the uncertainty in deriving CME parameters. We did not find a significant decrease in errors when going from two to three viewpoints for our specific hypothetical three spacecraft scenario using synthetic GCS white-light data. As we expected, considering all configurations and numbers of viewpoints, the mean absolute errors in the measured CME parameters are generally significantly higher in the case of the simulated MHD white-light data compared to those from the synthetic white-light images generated by the GCS model. We found the following CME parameter error bars as a starting point for quantifying the minimum error in CME parameters from white-light reconstructions: Δθ (latitude)=6° +2° -3° , Δϕ (longitude)=11° +18° -6° , Δγ (tilt)=25° +8° -7° , Δx (half-angle)=10° +12° -6° , Δh (height)=0.6 +1.2 -0.4 R ⨀ , and Δκ (ratio)=0.1 +0.03 -0.02 .

Coronal mass ejections↗

NASA NTRS2022

Quantifying Errors in 3D CME Parameters Derived From Synthetic Data Using White-Light Reconstruction Techniques

Coronal mass ejections↗

NASA NTRS1994

Dynamics of the magnetic shearing instability and magnetohydrodynamic turbulence in accretion disks. 1: Vertical magnetic field

Recently, the magnetic shearing instability (MSI) has been proposed as a dynamical mechanism for angular momentum transport in accretion disks (Balbus & Hawley 1991; Hawley & Balbus 1991). In this paper, the nonlinear dynamics of MSI modes in the presence of a vertical magnetic field B(sub 0) is discussed. In particular, the saturation levels of the fluctuating fields, the angular momentum flux, and the energy dissipation mechanism, are examined in detail. It is shown that MSI induces strong magnetohydrodynamic (MHD) turbulence in a range of wavenumbers 1/H is less than K is less than or equal to Omega/V(sub A)(sub 0)), where H is the thickness, Omega is the rotation frequency of the disk, and V(sub A)(sub 0) is the Alfven velocity. Despite the fact that the linear growth rate of MSI is maximal at small-scale (i.e., k is approximately Omega/V(sub A(sub 0)), angular momentum transport due to MSI turbulence is dominated by the magnetic Reynolds stress driven by large-scale modes (k is approximately 1/H). It is shown that the amplitude of low k(sub r) MSI eddies is limited primarily by subscale shear flow instability. Thus, dominant MSI cells are quasi-isotropic. In a stationary state, the effective Shakura-Sunyaev 'alpha' value is predicted to be of order V(sub A)(sub 0)0/C(sub s). In addition, the veritcal magnetic-field-induced MSI cells convert vertical magnetic field B(sub 0) into azimuthal magnetic field B(sub theta) in the disk. The generation of azimuthal magnetic field in turn introduces new physical processes, such as dynamo activity and azimuthal MSI turbulence. We conclude that it is not possible to decouple vertical MSI saturation from azimuthal MSI evolution. Low-frequency MSI cells are shown to co-exist with high-frequency radial buoyancy or internal waves. We show that modulational interaction between waves on these two frequency ranges is usually weak in the case when mean magnetic field is vertical. Thus, MSI and internal wave dynamics must be treated on an equal footing.

Zhang, W.↗

NASA NTRS2020

A Double Disturbed Lunar Plasma Wake

Under nominal solar wind conditions, a tenuous wake forms downstream of the lunar nightside. However, the lunar plasma environment undergoes a transformation as the Moon passes through the Earth's magnetotail, with hot subsonic plasma causing the wake structure to disappear. We investigate the lunar wake response during a passing coronal mass ejection (CME) on March 8, 2012 while crossing the Earth's magnetotail using both a magnetohydrodynamic (MHD) model of the terrestrial magnetosphere and a three-dimensional hybrid plasma model of the lunar wake. The CME arrives at 1 AU around 10:30 UT and its impact is first detected inside the geomagnetic tail after 11:10 UT by the Acceleration, Reconnection, Turbulence and Electrodynamics of the Moon's Interaction with the Sun (THEMIS-ARTEMIS) satellites in lunar orbit. A global magnetospheric MHD simulation using Wind data for upstream conditions with the OpenGGCM model reveals the magnetosheath compression to the lunar position from 11:20–12:00 UT, accompanied by multiple flux rope or plasmoid-like features developing and propagating tailward. MHD results support plasma changes observed by the THEMIS-ARTEMIS satellites. Lunar-scale simulations using the Amitis hybrid code show a short and misaligned plasma wake during the Moon's brief entry into the magnetosheath at 11:20 UT, with plasma expansion into the void being aided by the higher plasma temperatures. Sharply accelerated flow speed and a compressed magnetic field lead to an enhanced electric field in the lunar wake capable of generating sudden changes to the nightside near-surface electric potential.

A. P. Rasca↗

DOE OSTI2024

Time-resolved biphase signatures of quadratic nonlinearity observed in coupled Alfvén eigenmodes on the DIII-D tokamak

We report the detection of nonstationary quadratic coupling between toroidicity-induced Alfvén eigenmodes (TAEs) on sub-millisecond time scales. Identification of phase coherency between multiple TAEs and nonlinearly generated modes is facilitated by wavelet-based bicoherence analysis of time-series from inductive coils, taken from a DIII-D discharge heated by neutral beam injection (NBI). Characterization of nonlinear three-wave interaction is inferred by stationary local bispectrum phase (biphase) and confirmed via bandpass filtering. Biphase dynamics associated with prominent bispectral features are well-resolved in time and consistent with transient quadratic coupling. Onset and duration of nonlinearity are correlated with enhanced amplitude of participating TAEs; coincident changes in amplitude are observed for modes at difference frequency |f TAE,1 - f TAE,2 |.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗