Search NASA⌕ Search

SEARCH · Search NASA

Results for “full”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Comprehensive evaluation of commercially scalable atomic-layer-deposited alumina coating impact on full cell battery performance across varied test conditions

Atomic Layer Deposition (ALD) has emerged as a strategic enhancement method for lithium-ion battery (LIB) materials offering potential benefits and durability benefits for industrial battery production. However, the translation from laboratory achievements to commercial-scale applications has been limited. Here, this study aims to bridge this gap by comprehensively evaluating the effects of commercially scalable Al 2 O 3 ALD coatings using full pouch cell performance as a means to assess the ALD impact. We utilized large-scale slot-die coating techniques to ensure consistent electrode quality and tested four configurations of pouch cells to analyze the individual effects of ALD coating on anode and cathode electroactive materials. Our extensive testing matrix included long-term cycling, fast discharge, fast charge, leakage current, and high voltage tests. While at lower C-rates (<~1C), the influence of Al 2 O 3 coatings on cell performance is not significant. Fast charging conditions reveal that the anode ALD coating significantly enhances performance via a passivating effect, while on the cathode, it is detrimental, potentially due to increased resistance of the thin interfacial layer formed during the ALD processing. Leakage current and high-voltage tests show that the application of ALD coatings on either anode or cathode effectively minimizes side reactions at the electrode-electrolyte interface. Additionally, ALD coatings significantly mitigate concentrated and localized lithium plating on the anodes. These insights provide a valuable understanding of the potential of ALD technologies in LIB manufacturing to tailor cell performance, paving the way for safer, more efficient, and cost-effective battery solutions.

25 ENERGY STORAGE↗
DOE OSTI2025

Numerical simulation of coupled THM behaviour of full-scale EBS in backfilled experimental gallery in the Horonobe URL

Bentonite-based engineered barrier system (EBS) is a key component of many repository designs for the geological disposal of high-level radioactive waste. Given the complexity and interaction of the phenomena affecting the barrier system, coupled thermo-hydro-mechanical (THM) numerical analyses are a potentially useful tool for a better understanding of their behaviour. In this context, a Task (the Horonobe EBS experiment) was undertaken to study, using numerical analyses, the thermo-hydro-mechanical (and thermo-hydro) interactions in bentonite based engineered barriers within the international cooperative project DECOVALEX 2023. One full-scale in-situ experiment and four laboratory experiments, largely complementary, were selected for modelling. The Horonobe EBS experiment is a temperature-controlled non-isothermal experiment combined with artificial groundwater injection. The Horonobe EBS experiment consists of the heating and cooling phases. Six research teams performed the THM or TH (depended on research team approach) numerical analyses using a variety of computer codes, formulations and constitutive laws. Finally, for each experiment, the basic features of the analyses are described and the comparison between calculations and laboratory experiments and field observations are presented and discussed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Hydrogen production from full-strength corn stover fermentation effluent in single-chamber replaceable-cathode microbial electrolysis cells

Lignocellulosic residual biomass generated by the agricultural sector is an abundant feedstock for biohydrogen production via dark fermentation. However, this process is intrinsically inefficient, converting only ~30% of the reductant energy into H2 and leaving substantial amounts of reduced byproducts. These byproducts, mostly found in the fermentation effluents, can be further valorized in microbial electrolysis cells (MECs) to enhance the overall H2 recovery. However, current MEC configurations are typically dual- or single-chamber systems, yet both suffer from key inefficiencies. Dual-chamber systems rely on proton exchange membranes that are costly and prone to rapid biofouling, whereas single-chamber, membraneless systems are limited by reduced productivity due to H2 recycling and methanogenic consumption In this study, three single-chamber, 50-mL replaceable-cathode microbial electrolysis cells (RC-MECs) were 3D-printed and equipped with a physical separator to isolate anode and cathode compartments and limit H2 migration. Full-strength milled corn stover (MCS) fermentation effluent (COD of ~23.8 g-COD/L) was treated in fed-batch mode over two operational periods spanning 21 (Run 1) and 80 (Run 2) days. The RC-MECs exhibited comparable performance in both runs: after biofilm maturation, current densities exceeded 100 A/m²_cathode, COD removal reached up to 43%. Notably, extended RC-MECs operation led to a substantial methanogenic activity with the CH4 fraction in the cathode gas increasing to as high as 80% of the total biogas. Additions of a methanogenesis inhibitor 2-bromoethanesulfonate (2-BES) produced transient increases in hydrogen yields (11.51 and 5.12 L-H2/L_reactor/day in Runs 1 and 2, respectively); however, sustained 2-BES addition in subsequent cycles reduced total biogas production, decreased COD removal, and led to volatile fatty acid accumulation. Overall, single-chamber MECs can treat high-strength dark fermentation effluents while improving H2 recovery, but methanogenesis remains a key bottleneck, and complete long-term inhibition may be operationally unsustainable.

Hydrogen Production↗
DOE OSTI2025

The influence of kinematics of deformation on polycrystalline halite dynamic recrystallization: Full-field simulation of simple shear versus pure shear

Rock salt, composed mainly of halite, flows viscoplastically over a wide range of geological conditions, strongly impacting the dynamic evolution of sedimentary basins and orogens. Understanding how dislocation creep, which involves dislocation glide, intracrystalline recovery and dynamic recrystallization, influences the microstructure and rheology of halite under various deformation kinematics and temperatures is crucial for enhancing knowledge of salt flow dynamics. Here, this study employs a full-field numerical simulation method to compare the viscoplastic deformation of polycrystalline halite under simple shear and pure shear conditions up to a natural strain of ε = 1.5 at temperatures ranging from 100 °C to 300 °C. The results are presented in terms of crystallographic preferred orientation (CPO), grain shape preferred orientation (SPO), subgrain boundary direction, grain size and strain rate distribution. The results indicate that the crystallographic anisotropy of individual halite crystals is transferred to the polycrystalline scale, resulting in strain localization, particularly in simple shear simulations. The kinematics of deformation affect the evolution and distribution of high strain-rate bands, determining the direction of intragranular substructures and the morphology of strain-induced grain boundaries, with minimal impact on grain size. The intensity of grain boundary migration increases with temperature, significantly influencing grain morphology and size, thereby obscuring strain localization, while having little effect on CPOs. At low strain (ε < 1.0), CPOs relative to both the maximum shortening direction and the grain SPO are similar regardless of the deformation kinematics. At high strain (ε > 1.0), simple shear CPOs exhibit three stronger {100} maxima with a monoclinic symmetry relative to the grain SPO compared to the six {100} maxima with an orthotropic symmetry relative to the grain SPO generated under pure shear. Therefore, microstructures and CPOs can serve as indicators of the strain path in polycrystalline halite under various conditions, aiding in determining the shear sense and elucidating the deformation kinematics of salt structures.

58 GEOSCIENCES↗
DOE OSTI2025

Microstructure and CPO evolution of dynamically recrystallized olivine during complex deformation conditions: a full-field numerical modeling approach

The rheological properties of mantle rocks are strongly dependent on their crystallographic preferred orientation (CPO). Olivine CPO, defined by the orientation of seismically fast [100] axes parallel to flow direction, is also thought to be a dominant contributor to seismic anisotropy in the Earth's upper mantle. However, the amount of deformation needed to overprint a new CPO on a pre-existing fabric and the impact of the inherited CPOs on the transient microstructure evolution, remain unknown. This study employs a full-field numerical approach (VPFFT-ELLE) to explore the dynamic recrystallization and microstructural evolution of olivine polycrystalline aggregates under complex deformation conditions. We test four combinations of successive pure shear and simple shear boundary conditions. Findings indicate that inherited CPOs influence subsequent deformation in a manner dependent on the kinematic relationship between successive stages. In all cases, a minor strain increment (ε ∼0.3–0.6) is sufficient to erase the previous microstructure and CPO. However, when deformation conditions change dramatically (e.g., stretching direction changes orthogonally), the intensity of the new CPO developed is significantly lower and strain distribution are specially altered. During a transient strain stage, pre-existing microstructures undergo extensive reworking, especially when deformation conditions are changed dramatically, such as switching from simple shear to a pure shear condition with a parallel shortening direction relative to the stretching direction. We estimate the significance of these results in interpreting observations of seismic velocity anisotropy, concluding that P-wave seismic anisotropy is significantly and positively correlated with the evolution of olivine CPO with deformation history. This research underscores the transient nature of microstructural rearrangement in olivine aggregates and the necessity for caution in interpreting seismic anisotropy in regions with complex deformation histories, as inherited CPOs can influence current fabric development and induce deviation to the present deformation conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A full-scope, high-fidelity simulator-based hardware-in-the-loop testbed for comprehensive nuclear power plant cybersecurity research

Nuclear power plant (NPP) cybersecurity research often relies on hardware-in-the-loop (HIL) testbeds that integrate real hardware components into simulated environments. These testbeds allow researchers to identify vulnerabilities, evaluate attack impacts, and test security measures in a controlled setting. Furthermore, previous HIL testbeds lacked fidelity to accurately represent real nuclear systems, limiting the scope of cybersecurity analysis. This study presents the creation of a HIL testbed, devised upon a full-scope, high-fidelity NPP simulator, to facilitate realistic and comprehensive cybersecurity research. To demonstrate its capabilities, the control logic for the steam generator water level was migrated from the simulator to an external programmable logic controller. As a practical application of the developed testbed, supply chain attack scenarios were simulated by injecting malicious code into the controller logic, and the effects of manipulating sensor inputs and control commands were observed. While this HIL testbed provides more detailed simulations, enhanced realism, and wider applicability compared to other options utilizing a less complex simulator, it is also more intricate and costly. For this reason, we include a detailed comparison with some alternative architectures to aid fellow researchers and practitioners in the selection of a suitable HIL architecture based on specific research objectives.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Quantitative DNA Stable Isotope Probing Identifies Active Microorganisms Assimilating Volatile Fatty Acids in Full-Scale Enhanced Biological Phosphorus Removal Processes

Enhanced biological phosphorus removal (EBPR) systems often rely on exogenous carbon sources, such as volatile fatty acids (VFAs), to achieve higher P removal. Here, we employed DNA quantitative stable isotope probing (qSIP) using two VFAs, acetate and propionate, in cyclic anaerobic/aerobic incubations to assess their effects on P cycling and microbial activity with biomass from two full-scale EBPR water resource-recovery facilities that utilize VFA addition. We found that anaerobic VFA uptake preferences differed within known groups of PAOs, such as Candidatus Accumulibacter and Tetrasphaera-affiliated members (e.g., Ca. Phosphoribacter), between the two biomasses. The combination of qSIP with metagenomics identified isotopically labeled phages that were linked to active PAOs, highlighting their potential roles in modulating EBPR community composition and activity. The highest levels of anaerobic labeling from acetate were in genomes belonging to Saccharimonadales and Rickettsiales, which are generally host-associated with bacteria and eukaryotes, respectively. Furthermore, this finding highlights the possibility of cross-feeding between PAO hosts and their parasites or predators, as well as the role of so-far uncharacterized organisms participating in carbon cycling under EBPR conditions. Collectively, these results expand our understanding of the ecological interactions involved in communities anaerobically uptaking VFAs and cycling P that are central to EBPR.

Polyphosphate accumulating organisms↗
DOE OSTI2025

Linear-Scaling Local Natural Orbital-Based Full Triples Treatment in Coupled-Cluster Theory

We present an efficient, asymptotically linear-scaling implementation of the canonically O(N 8 ) coupled-cluster method with singles, doubles, and full triples excitations (CCSDT) method. We apply the domain-based local pair natural orbital (DLPNO) approach for computing CCSDT amplitudes. Our method, called DLPNO–CCSDT, uses the converged coupled-cluster amplitudes from a preceding DLPNO–CCSD(T) computation as a starting point for the solution of the CCSDT equations in the local natural orbital basis. To simplify the working equations, we t1-dress our two-electron integrals and Fock matrices, allowing our equations to take on the form of CCDT. With appropriate parameters, our method can recover more than 99.99% of the total canonical CCSDT correlation energy. In addition, we demonstrate that our method consistently yields sub-kJ mol –1 errors in relative energies when compared to canonical CCSDT, and, likewise, when computing the difference between CCSDT and CCSD(T). Finally, to highlight the low scaling of our algorithm, we present timings on linear alkanes (up to 30 carbons and 730 basis functions) and water clusters (up to 131 water molecules and 3144 basis functions).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

NOCI-F Electronic Couplings in Assemblies of Indolonaphthyridine Molecules: From Dimers to the Full Stack

Key electronic processes related to molecular excitonic states of finite stacks of indolonaphthyridine molecules are analyzed via the non-orthogonal configuration interaction with fragments (NOCI-F) method. Indolonaphthyridine is an organic chromophore that can undergo several electronic photoexcitation-related intermolecular processes, such as exciton and electron transfer. The structures studied here are noncrystalline arrangements built as either ordered stacks of indolonaphthyridine or stacks extracted from molecular dynamics simulations including thermal disorder. Taking dimers or trimers from either model, we performed CASSCF and NOCI-F calculations to quantify the intermolecular electronic couplings governing singlet fission, excited singlet and triplet diffusion, and hole and electron diffusion processes. Also, comparing the results for the different models, we studied the effect of structural disorder and distortion on these couplings. Finally, we present a newly developed, advanced postanalysis tool. It takes the NOCI-F data as input to carry out a multifragment full Hamiltonian procedure that involves the complete stack, providing physical information not available from the dimer/trimer models, hence giving access to additional insight into the material’s properties.

coupling reactions↗
DOE OSTI2026

Local Pair Natural Orbital-Based Coupled-Cluster Theory through Full Quadruples (DLPNO–CCSDTQ)

In this work, we implement a local pair natural orbitalbased coupled-cluster method through the full treatment of quadruple excitations (CCSDTQ). The domain-based local pair natural orbital (DLPNO) approach, which has successfully been applied to lower levels of coupled-cluster theory, is utilized in our algorithm, and thus our algorithm is called DLPNO-CCSDTQ. For simplicity in the working equations and in the implementation, we t 1 -dress the twoelectron integrals as well as Fock matrix elements. Our method can recover CCSDTQ-CCSDT and CCSDTQ-CCSDT(Q) energy differences on the order of 0.01−0.05 kcal mol −1 , even at a loose quadruples natural orbital (QNO) occupation number cutoff of 3.33 × 10 −6 . To highlight the capabilities of our code and its potential future applications, we showcase computations that would be intractable with canonical CCSDTQ, such as the benzene dimer, (H 2 O) 17 , and adamantane. With sufficient computing resources, computations up to 15 heavy atoms (40 atoms overall) may be feasible for fully bonded 3D systems.

Cluster chemistry↗
DOE OSTI2025

Many-Body Basis Set Amelioration Method for Incremental Full Configuration Interaction

Incremental full configuration interaction (iFCI) is a polynomial-cost electronic structure method that systematically approaches the FCI limit by employing the method of increments to solve the Schrödinger equation through a many-body expansion. This article introduces the many-body basis set amelioration (MBBSA) method, which is designed to allow iFCI to be applicable to larger atomic orbital basis sets. MBBSA uses a series of inexpensive iFCI calculations to approximate the correlation energy that would be found using a more expensive, highly accurate iFCI calculation. Here, when compared to standard iFCI computations on smaller molecules in triple-zeta and larger basis sets, MBBSA provides approximations to the total and relative energies within chemical accuracy. MBBSA exhibits a reduced cost of between 60-92% when compared to standard iFCI calculations, with larger systems experiencing the largest benefit. Tests of MBBSA on two reactions that involve highly correlated systems, the automerization of cyclobutadiene and a Criegee intermediate reaction, show that MBBSA has practical utility for studying realistic chemistries.

Basis sets↗
DOE OSTI2025

A Solvatochromic Near Infrared Fluorophore Sensitive to the Full Amyloid Beta Aggregation Pathway

Alzheimer's disease has long been associated with the aggregation of amyloid beta peptides (Aβ42) into macroscale plaques, although specific neurodegenerative agents have not been definitively identified. Much evidence has pointed to the soluble nanoscale oligomers that form early in the Aβ42 aggregation pathway, but there is little understanding of these structures, their mechanisms of formation, or how they grow into plaques. Here, we show that a solvatochromic fluorophore with near-infrared (NIR) emission can track synthetic Aβ42 aggregation through environment-sensitive spectral shifts from the earliest time points through plaque formation. This azide-functionalized phosphine oxide azetidine rhodol (Phazr-N3) shows large polarity-dependent changes in fluorescence emission, with maxima shifting from 630 nm in toluene to 703 nm in aqueous buffer, and a maximum quantum yield of 62%. Upon induction of Aβ42 aggregation, we observe immediate solvatochromic changes in Phazr-N3 fluorescence, with multiple apparent phases over 12 h, and which culminate before the onset of any major fluorescence changes of conformation-specific aggregation fluorophore thioflavin T. Solution anisotropy measurements show a low micromolar affinity of Phazr-N3 for disordered, free Aβ42 in solution, and real-time measurements are consistent with rapid liquid-liquid phase separation and slow dehydration of the growing aggregate. Spectral imaging of synthetic plaques stained in the presence of live cells and lipid-binding protein albumin shows over 4000-fold Phazr-N3 fluorescence intensity above background under no-wash conditions, and over 100-fold intensity above coplated microglial cells or a large excess of albumin. This use of a solvatochromic probe with structure-independent binding to free Aβ42 offers real-time, minimally invasive insight into the full Aβ42 aggregation pathway.

Wang, Zeming↗
DOE OSTI2024

WUS324: Multiscale Full Waveform Inversion Approaching Convergence Improves Waveform Fits While Imaging Seismic Structure of the Western United States

Abstract We report a new model of radially anisotropic crustal and upper mantle structure of the western United States (WUS324) obtained from full waveform inversion of earthquake data. We ran three multiscale inversion stages beyond model WUS256 (Rodgers et al., 2022, https://doi.org/10.1029/2022jb024549 ) allowing them to approach convergence to fit a larger data set to a shorter minimum period of 16 s. WUS324 is based on 324 total iterations from its starting model, significantly more (16 times) than previous studies. Waveform misfit reductions are 66%–70% for both the inversion data and an independent validation data set providing confidence in the predictive power of the model. WUS324 provides much better fits and reveals shear wavespeed, v S , structure of this large region with more detail than previous waveform tomography models. We show representative images demonstrating the resolution of diverse seismic structure across this highly heterogeneous region including oceanic lithosphere, subducting slabs and continental magmatism.

58 GEOSCIENCES↗
DOE OSTI2026

Regularization by denoising diffusion models for solving inverse PDE problems with application to full waveform inversion

Partial differential equation (PDE)-governed inverse problems are fundamental across various scientific and engineering applications; yet they face significant challenges due to nonlinearity, ill-posedness, and sensitivity to noise. Here, we introduce a computational framework, regularization by denoising using diffusion models for partial differential equations (RED-DiffEq), by integrating physics-driven inversion and data-driven learning. RED-DiffEq leverages pretrained diffusion models as a regularization mechanism for PDE-governed inverse problems. We apply RED-DiffEq to solve the full waveform inversion problem in geophysics, a challenging seismic imaging technique that seeks to reconstruct high-resolution subsurface velocity models from seismic measurement data. Our method shows enhanced accuracy and robustness compared to benchmark methods. Additionally, it exhibits strong generalization and domain decomposition capacity, enabling the inversion of more complex velocity models with larger domains than those used in training the diffusion model. Our framework can also be directly applied to diverse PDE-governed inverse problems.

Shan, Siming [Yale University, New Haven, CT (Unit↗
DOE OSTI2026

Multi-resolution enhancement for full-spectrum neural representations

Scientific data acquisition continues to outpace storage and analysis capabilities, making voxel-basedrepresentations increasingly intractable. Implicit neural representations (INRs) offer a promising solutionby encoding signals through coordinate-based neural networks, serving as surrogates of data, withcomputational and storage requirements scaling with network complexity rather than data dimensionality.However, smaller INRs struggle to faithfully represent multiscale structures, high-frequency informationand fine textures that constitute a large proportion of scientific measurements. We propose WIEN-INR, atheoretically guided hierarchical INR framework that distributes modelling across resolution scales andenables improved representation capacity through a novel enhancement network to recover subtle details.This multiscale architecture allows smaller networks to retain the full spatial-frequency content of thesignal as well as preserve training efficiency and lower storage cost. Evaluated on distinct raw experimentalmeasurements across scales and complexities, WIEN-INR represents a practical step towards a broaderadoption of neural representations in scientific workflows, delivering compact, robust and high-fidelityrepresentations.

Ni, Yuan [SLAC National Accelerator Laboratory (SL↗
DOE OSTI2024

Virtual node graph neural network for full phonon prediction

Understanding the structure-property relationship is crucial for designing materials with desired properties. The past few years have witnessed remarkable progress in machine-learning methods for this connection. However, substantial challenges remain, including the generalizability of models and prediction of properties with materials-dependent output dimensions. Here we present the virtual node graph neural network to address the challenges. By developing three virtual node approaches, we achieve Γ-phonon spectra and full phonon dispersion prediction from atomic coordinates. We show that, compared with the machine-learning interatomic potentials, our approach achieves orders-of-magnitude-higher efficiency with comparable to better accuracy. This allows us to generate databases for Γ-phonon containing over 146,000 materials and phonon band structures of zeolites. Additionally, our work provides an avenue for rapid and high-quality prediction of phonon band structures enabling materials design with desired phonon properties. The virtual node method also provides a generic method for machine-learning design with a high level of flexibility. In this study, the authors present a virtual node graph neural network to enable the prediction of material properties with variable output dimensions. This method offers fast and accurate predictions of phonon band structures in complex solids.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Full-polarization millimeter wavelength variability of Sagittarius A * during the 2018 EHT campaign

Context. Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, provides a unique laboratory to study accretion dynamics and plasma processes near the event horizon. Aims. We investigated the variability and polarization properties of Sgr A* using ALMA observations during the 2018 Event Horizon Telescope campaign. Methods. We analyzed high-cadence full-polarization light curves from ALMA at millimeter wavelengths, performed time-series analysis, and investigated the temporal behavior during an X-ray flare observed by Chandra on 2018 April 24. The variability characteristics are compared with expectations from standard accretion flow models. Results. We find low variability in total intensity (σ/μ < 10%), but significantly higher variability in linear and circular polarization (∼30% and ∼50%, respectively). A time-series analysis reveals red-noise variability, with power spectral densities between −2 and −3 across all Stokes parameters. Polarized intensity shows stable intra-day timescales, while total intensity exhibits more variable timescales, suggesting distinct emission regions, with polarization likely arising from a coherent structure. On April 24, a statistically significant inter-band delay in polarized intensity coincides with a near-simultaneous X-ray and millimeter peak that deviates from the typical delayed flare scenario. This event also features enhanced millimeter variability and coherent polarization loop evolution. The observed simultaneity challenges standard models of transient synchrotron emission with cooling delays, favoring instead a scenario of continuous energy injection in an optically thin region. Conclusions. Our results offer new constraints on the physical mechanisms driving variability in Sgr A*, and provide key observational input for refining theoretical models of accretion and plasma behavior in the vicinity of supermassive black holes.

Galaxy: center↗
DOE OSTI2026

Full-wave simulation of ion cyclotron range of frequency heating in a mirror device

We are adapting the all-orders spectral algorithm (AORSA) [Jaeger Phys. Plasmas 8 (2001)] from tokamak geometry to axisymmetric magnetic mirror geometry. We will use the modified AORSA to study ion cyclotron range of frequency (ICRF) heating in the Wisconsin HTS Axisymmetric Mirror (WHAM) magnetic mirror device [D. Endrizzi J. Plasma Phys. 5 (2023)]. ICRF power will be used in WHAM to accelerate high-energy neutral beam injected deuterium ions with ion cyclotron absorption at the second to fourth harmonics of deuterium. At these harmonics, for the spatial scales present in WHAM, full-wave simulations most accurately capture wave propagation, including transmission, reflection, and absorption at the cyclotron layer. The antenna is modeled as a single strap with m = 0 excitation, where m is the azimuthal mode number. We present some preliminary results for wave coupling and heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗