Search NASA⌕ Search

SEARCH · Search NASA

Results for “proxy”

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 541 records · Page 30

Paleoclimatic implications of glacial fluctuations in the Sierra Nevada del Cocuy, northern Andes, Colombia, during the Lateglacial and Holocene

The reconstruction of former mountain glaciers from geomorphic mapping and cosmogenic-nuclide surface-exposure dating provides a unique opportunity to infer patterns of past terrestrial climate variability. Tropical mountain glaciers are particularly valuable as there are comparatively few terrestrial climate proxies at equatorial latitudes relative to higher latitudes. As the single largest climate zone on Earth, the tropics play an outsized role in mediating global climate via the ocean-atmosphere transfer of latent heat and water vapor. Nonetheless, there remains a persistent gap in our understanding of how the tropics influenced – or were influenced by – the high-magnitude climate shifts of the Late Pleistocene, and whether this high-energy region simply responded to extratropical forcing or was itself a driver of global climatic change. To help address this knowledge gap, we analyzed geologic evidence for past glacial fluctuations in three adjacent valleys in the Sierra Nevada del Cocuy, the highest subrange of the Eastern Cordillera in the Colombian Andes, to provide a terrestrial record of atmospheric temperature during the latter part of Termination 1. Coupled with geomorphic mapping and paleo-snowline reconstructions, our beryllium-10 glacial chronology indicates that glaciers in the humid inner tropics underwent pronounced growth and gradual decay during the Antarctic Cold Reversal (14.5–12.8 ka) and Younger Dryas (12.8–11.7 ka) periods, respectively, following a trend that, according to directly dated moraine records from throughout both polar hemispheres, appears to have been global. While the specific mechanism(s) behind this large-scale behavior remains to be corroborated, we revisit the hypothesis that ocean-atmosphere heat transfer and water vapor flux are key drivers of abrupt Lateglacial temperature fluctuations. Subsequent to the Lateglacial, deglaciation of the Sierra Nevada del Cocuy accelerated during the Early Holocene, a pattern also observed in other tropical glacier records. More recently, the magnitude of snowline rise and glacier retreat over the last two centuries supports the view that modern tropospheric warming is anomalously strong at least relative to the last ∼16,000 years.

Andes↗

Downwind coning concept rotor for a 25 MW offshore wind turbine

The size of offshore wind turbines over the next decade is expected to continually increase due to reduced balance of station costs per MW and also the higher wind energy at increased altitudes that can lead to higher capacity factors. However, there are challenges that may limit the degree of upscaling which is possible. In this paper, a two-bladed downwind turbine system is upscaled from 13.2 MW to 25 MW, by redesigning aerodynamics, structures, and controls. In particular, three 25 MW rotors have been developed: V1 is the upscaled model, V2 is a partial redesigned model, and V3 is a fully redesigned model. Despite their radically large sizes, it is found that these 25 MW turbine rotors satisfy this limited set of design drivers at the rated condition and that larger blade lengths are possible with cone-wise load-alignment. In addition, flapwise morphing (varying the cone angle with a wind-speed schedule) is investigated in terms of minimizing mean and fluctuating root bending loads using steady inflow proxies for the maximum and damage equivalent load moments. The resulting series of 25 MW rotors, which are the largest ever designed, can be a useful baseline for additional development and assessment.

17 WIND ENERGY↗

Machine learning and shallow groundwater chemistry to identify geothermal prospects in the Great Basin, USA

This study discovers various geothermal prospects in the Great Basin, USA based on shallow groundwater chemical (geochemical) data. The geochemical data are expected to include hidden (latent) information that is a proxy for geothermal prospectivity. We processed the sparse geochemical data in the Great Basin at 14,341 locations including 18 attributes. Next, a non-negative matrix factorization with customized k-means clustering is applied to the geochemical data matrix that automatically finds three hidden geothermal signatures representing modestly, moderately, and highly confident geothermal prospects. The algorithm also evaluated the probability of occurrence of these types of resources through the studied region. There is a consistency between regional geothermal prospectivity as estimated by our ML methodology and the traditional play fairway analysis conducted over a portion of the study area. We also identify the dominant data attributes associated with each signature. Finally, our ML analyses allow us to reconstruct attributes from sparse into continuous over the study domain. The predicted continuous attributes can be used for future detailed geothermal explorations in the Great Basin.

15 GEOTHERMAL ENERGY↗

Towards multi-fidelity deep learning of wind turbine wakes

We report engineering wake models that accurately predict wake in a computationally efficient manner are very important for tasks such as layout optimization and control of wind farms. In this paper, we explore an application of deep learning (DL) to learn the wake model from hierarchies of physics-based approaches ranging from analytical models to an approximate form of the Reynolds-averaged Navier-Stokes equations. We first illustrate the application of principal component analysis to obtain a lower-dimensional representation that allows a computationally tractable training and deployment of DL models. Then, the DL model is trained to learn the mapping from input parameter space to the principal components, which are then used to reconstruct the three-dimensional flow field. Additionally, we investigate a composite framework consisting of two neural networks to learn the correlation between low- and high-fidelity data with Gauss and curl models treated as proxies for low- and high-fidelity models, respectively. The prediction from both DL models matches well with the high-fidelity data with a maximum relative percentage error for the kinetic energy flux of <1%. This work opens up possibilities for data-efficient construction of surrogate models for wake prediction that can be used to study the influence of wind speed and yaw angles on wind farm power production.

17 WIND ENERGY↗

TROPOMI observations allow for robust exploration of the relationship between solar-induced chlorophyll fluorescence and terrestrial gross primary production

Solar-induced chlorophyll fluorescence (SIF) observed by satellites has advanced the monitoring of terrestrial photosynthesis regionally and globally. The relationship between SIF and gross primary production (GPP) at leaf, canopy, and ecosystem scales has received tremendous attention in recent years. It remains controversial whether the SIF-GPP relationship at the ecosystem scale is universal or dependent upon vegetation type. New SIF observations from the TROPOspheric Monitoring Instrument (TROPOMI) with unprecedented high spatial and temporal resolution provide a new opportunity to elucidate the SIF-GPP relationship. Here, we examine the SIF-GPP relationship for seven major vegetation types across the U.S. with TROPOMI SIF and in-situ GPP data for 83 eddy covariance flux sites. We find that TROPOMI SIF shows a strong and consistent relationship with tower-based GPP at both satellite footprint and grid-cell levels. The slope of the SIF-GPP relationship is similar among all the vegetation types except croplands, demonstrating a nearly universal (converging to ~13.5 g C m -2 d -1 /W m -2 μm -1 sr -1) rather than vegetation type-specific SIF-GPP relationship. This confirms that TROPOMI SIF can be used as a proxy for GPP across a wide variety of vegetation types and can also be used to quantify GPP by avoiding uncertainty associated with land cover maps. The C 4 crops have a much higher slope than the C 3 crops, and therefore croplands tend to have a higher slope than C 3 -dominated vegetation types (e.g., forests, shrublands, savannas). We also find that the TROPOMI SIF is well correlated with GPP under normal or wetter conditions, while their relationship becomes weaker under water stress. In conclusion, our TROPOMI-based study could improve our understanding of the SIF-GPP relationship at the ecosystem scale and advance the mapping of GPP globally with SIF observations from space.

54 ENVIRONMENTAL SCIENCES↗

Decoupling of greenness and gross primary productivity as aridity decreases

Ecosystem primary productivity is a key ecological process influencing many ecosystem services, including carbon storage. Thus, clarifying how primary productivity in terrestrial ecosystems responds to climatic variability can reveal key mechanisms that will drive future changes in the global carbon budget. Satellite products of canopy greenness are widely used as proxies for vegetation productivity to evaluate how ecosystems respond to climate variability. However, to what degree inter-annual variations in productivity are consistent with greenness and how this relationship varies spatially remains unclear. Here we investigated the strength of the coupling between inter-annual variations in leaf area index (LAI, a measure of greenness) and ecosystem gross primary productivity (GPP) derived from eddy covariance towers, i.e., the r2 of the LAI-GPP relationship. Overall, inter-annual GPP and LAI were highly coupled (i.e., high r2) in arid grasslands, but were fully decoupled in mesic evergreen broadleaf forests, indicating that this relationship varies strongly along aridity gradients. A possible mechanism of the spatial variation in the LAI-GPP relationship is that the tradeoff between ecosystem structure (LAI) and physiology (photosynthesis per unit leaf area) becomes stronger in more humid climates. Land models overestimated the r2 of LAI-GPP correlation for most ecosystem types and failed to capture the spatial pattern along aridity gradients. We conclude that relying on greenness products for evaluating inter-annual changes in vegetation productivity may bias assessments, especially in tropical rainforest ecosystems. Our findings may also reconcile observed disparities between responses in greenness and GPP during drought in Amazon forests.

54 ENVIRONMENTAL SCIENCES↗

Improving estimates of sub-daily gross primary production from solar-induced chlorophyll fluorescence by accounting for light distribution within canopy

Solar-induced chlorophyll fluorescence (SIF) has long been regarded as a proxy for photosynthesis and has shown superiority in estimating gross primary production (GPP) compared to traditional vegetation indices, especially in evergreen ecosystems. However, current SIF-based GPP estimations regard the canopy as a large leaf and seldom consider the impact of interactions among light, canopy structure, and leaf physiology. In this study, we proposed GPP estimation models with different descriptions of light–structure–physiology interactions (including the layered model, the two-leaf model, and the layered two-leaf model) and compared their performances with the big-leaf model using half-hourly (or hourly) observations at evergreen needleleaf forest sites. First, we found that the big-leaf model underestimated GPP, especially at noon. All models showed higher accuracy than that of the big-leaf model. Second, we investigated the diurnal dynamics of GPP estimations in each canopy layer and found that models with a two-leaf assumption captured the diurnal variations in GPP better than that with the layered assumption. Here we also deduced that the poor performance of the big-leaf model was related to its overestimation of the overall light stress on the redox state of PSII reaction centers (qL). Finally, we noticed that the qL at the canopy scale had lower sensitivity to light change than the single-leaf qL and that the light response of canopy-scale qL was influenced by the leaf area index during seasonal cycles. Overall, this study describes methods to accurately estimate sub-daily GPP from SIF in evergreen needleleaf forests and demonstrates that the interactions among light, canopy structure, and leaf physiology regulate the SIF-GPP relationship at the canopy scale. Further, it indicates the need to consider the description of light distribution within the canopy in next-generation terrestrial biosphere models, even if they incorporate SIF to constrain their parameterization. Thus, upscaling the established leaf-scale mechanistic SIF-GPP relationship or findings to canopy-scale applications still requires much work, especially when there are significant changes in environmental conditions and their within-canopy distributions.

59 BASIC BIOLOGICAL SCIENCES↗

Enhancing analytical merits of laser-induced breakdown spectroscopy of hydrogen isotopes using an orthogonal double-pulsing scheme

Accurate detection and quantification of hydrogen isotopes in solid materials are vital for diverse applications, including fusion energy, hydrogen storage, and tritium production. Laser-induced breakdown spectroscopy (LIBS) is a well established, rapid, standoff method for this purpose, but it faces challenges related to the analytical merits required for isotopic analyses. In this study, we enhance the analytical and detection capabilities of traditional single-pulse LIBS by implementing an orthogonal double pulsing approach, focusing on the analysis of a range of 2 H concentrations in Zircaloy-4 substrates (acting as a proxy for 3 H). The double-pulse experiments employed an orthogonal re-heating configuration with two nanosecond Nd:YAG lasers. We systematically evaluated critical parameters affecting the signal intensity in double-pulse LIBS, including interpulse delay, ambient gas pressure, and heating laser energy. Finally, our results demonstrate that employing an orthogonal double-pulse scheme significantly enhances 2 H α emission while minimizing line broadening and self-absorption, ultimately improving the technique’s analytical capabilities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Combining organic amendments with enhanced rock weathering shifts soil carbon storage in croplands

Enhanced rock weathering (ERW) involves applying crushed silicate minerals to cropland soils to remove carbon dioxide and stabilize the global climate. If practiced widely, ERW has the potential to mitigate climate change and improve soil health and crop productivity. However, most ERW studies emphasize inorganic carbon (IC) chemistry, using model-based estimates and short-term mesocosms. Limited field data exist on how ERW interacts with organic amendments to affect organic carbon (C) cycling in soils. In a three-year field study in conventionally managed, irrigated maize fields, we monitored how key soil variables responded to crushed rock-alone, and in combination with compost and/or biochar. We measured weathering indicators (pH, major cations, and IC contents) and organic fractions, including particulate organic matter (POM), mineral-associated organic matter (MAOM), microbial biomass C, and water-extractable organic C. Rock-alone treatments increased weathering proxies (pH and IC) and showed an increasing trend in POM and MAOM, relative to control. In contrast, combining crushed rock with organic amendments resulted in lower soil organic C and nitrogen (N) concentrations (in both POM and MAOM) compared to organic amendments alone, though IC increased in the rock+compost treatment. Combining rock with both compost and biochar (compost/biochar) significantly lowered MAOM-N compared to compost/biochar alone. Overall, co-applying rock with organic inputs may promote weathering and C accrual but slow the accrual rate of organic C and N relative to organic amendments alone. Quantifying these trade-offs over multiple years and scales is critical to integrating ERW with existing soil health practices and climate mitigation strategies.

Biological and medical sciences↗

Effects of warming on bacterial growth rates in a peat soil under ambient and elevated CO 2

Boreal peatlands are important global carbon reservoirs that are particularly vulnerable to predicted climate changes such as increasing CO 2 and temperature. Since microbial activities regulate the balance of carbon sequestered into soil organic matter or remineralized to CO 2 , characterizing their response to these environmental factors is critical to predicting how peatland ecosystems will affect climate-carbon cycle feedbacks. Here we examined in-situ taxon-specific variation in microbial growth under long-term elevated CO 2 and across a gradient of warming treatments in a northern Minnesota peat bog using quantitative stable isotope probing with 18 O-water. Across temperatures, bacterial taxa were grouped according to the excess atom fraction 18 O (EAF) of their genomes, a proxy for DNA replication and hence, growth. Taxon-specific growth across CO 2 and temperature treatments clustered into relatively few response patterns. While a large portion of taxa showed little to no growth under ambient CO 2 , many of the same taxa grew rapidly under elevated CO 2 . We found support for phylogenetic conservation of response patterns among Acidobacteria and Proteobacteria, the two most abundant phyla in our data. Our results suggest certain taxa may be primed for new climate conditions and have a greater influence on carbon cycling with implications for future climate mitigation strategies.

16S amplicon sequencing, Carbon Dioxide (CO2), pea↗

Intensity measure adequacy assessment for nonlinear site response using Information Theory

Nonlinear site response modeling is a crucial aspect of Probabilistic Seismic Hazard Analysis. Site amplification models routinely rely on a rock intensity measure to characterize the strength of the bedrock motion. However, the adequacy of such intensity measures towards predicting amplifications across the oscillator period range has not been investigated in the literature. This paper analyzes the adequacy of rock intensity measures using state of the art criterion established in Performance-Based Earthquake Engineering and techniques from Information Theory. The efficiency and the sufficiency of several rock intensity measure are assessed. It was found that spectral accelerations at low periods usually are adequate for predicting amplifications across the period range. This supports the current practice of using Peak Ground Acceleration in Ground Motion Models. However, for extremely soft sites, which demonstrate nonlinear effects well into the period range, it is better practice to ensure that amplification and spectral acceleration share the same oscillator period. Lastly, for predicting the peak shear strain (an important parameter that controls nonlinearity of site response), Peak Ground Velocity is generally adequate, and this conclusion is in line with the definition of proxy shear strain.

42 ENGINEERING↗

Correlation between infrared spectra features and coverage of different adsorption sites for the NO/Pd(111) system

This study investigates the adsorption of nitric oxide (NO) on Pd(111) surfaces at 200 K using infrared reflection absorption spectroscopy (IRRAS). Peak positions and areas are used as proxies for tracking coverage of different species, serving as groundwork for in-situ experiments that need real-time tracking of surface species. Here, we derived a mathematical model correlating NO coverage with dose, consistent with the Molecular Langmuir Model. Distinct correlations are observed between coverage and spectral features (peak position and peak area for the NO vibration): at lower coverages (θ<0.6), the total integrated peak areas are linearly corelated to coverage; at higher coverages (θ>0.6), the peak center for a compressed-hollow site becomes the primary coverage indicator, showing a linear relationship between wavenumber and coverage. Note that this study addresses coverages below saturation. These findings refine our understanding of NO adsorption and establish a foundation for real-time in situ transient kinetics studies of chemical reactions involving NO on Pd(111).

25 ENERGY STORAGE↗

Gaining Perspective on Unconventional Well Design Choices through Play-level Application of Machine Learning Modeling

The recent development of unconventional oil and gas (O&G) reservoirs has led to an abundant hydrocarbon supply, both domestically and globally. However, there is a continued push to develop new and innovative approaches to improve exploration and extraction efficiencies and overall well productivity moving forward. Substantial improvements in unconventional O&G development are expected through optimized well completion and stimulation strategies aimed at maximizing well productivity. Optimizing well designs will require tailoring to the distinctive geologic conditions present for any newly placed well. To better evaluate the impact of well design attributes and their associated interactions on productivity in a major unconventional play, multivariate machine learning-based models that use empirical datasets were developed. A gradient boosted regression tree (GBRT) algorithm was applied. GBRT has been narrowly investigated for O&G applications but enables straightforward parametric importance and influence evaluation, as well as assessment of parameter interaction effects. Models were trained on well design and locational parameters that serve as a proxy for variable geologic conditions to estimate two types of productivity indicator response variables strongly correlated to estimated ultimate recovery (EUR). The dataset utilized consists of over 7,000 well observations that cover the majority of the productive region of the Marcellus Shale. Model performance was evaluated and algorithm parameters tuned by analyzing the goodness-of-fit for simulated results against observed data in a cross-validation approach. Models were found capable of 73–79 percent prediction accuracy on held out testing data of gas equivalent production and can be used to inform future well design and placement decisions for increasing EUR per well and improving overall field-level recovery. Study results indicate that Marcellus well performance improves most with upscaling perforated interval lengths and water and proppant volumes per foot; but relative productivity improvements are spatially dependent across the play. Finally, optimal combinations of water and proppant on well performance were found to vary depending on well location, emphasizing the utility of data-driven models capable of broad application across a play of interest for informing tailored well design approaches prior to their field deployment.

04 OIL SHALES AND TAR SANDS↗

Morphotype-resolved characterization of microalgal communities in a nutrient recovery process with ARTiMiS flow imaging microscopy

Microalgae-driven nutrient recovery represents a promising technology for phosphorus removal from wastewater while simultaneously generating biomass that can be valorized to offset treatment costs. As full-scale processes come online, system parameters including biomass composition must be carefully monitored to optimize performance and prevent culture crashes. In this study, flow imaging microscopy (FIM) was leveraged to characterize microalgal community composition in near real-time at a full-scale municipal wastewater treatment plant (WWTP) in Wisconsin, USA, and population and morphotype dynamics were examined to identify relationships between water chemistry, biomass composition, and system performance. Two FIM technologies, FlowCam and ARTiMiS, were evaluated as monitoring tools. ARTiMiS provided a more accurate estimate of total system biomass, and estimates derived from particle area as a proxy for biovolume yielded better approximations than particle counts. Deep learning classification models trained on annotated image libraries demonstrated equivalent performance between FlowCam and ARTiMiS, and convolutional neural network (CNN) classifiers proved significantly more accurate when compared to feature table-based dense neural network (DNN) models. Across a two-year study period, Scenedesmus spp. appeared most important for phosphorus removal, and were negatively impacted by elevated temperatures and increase in nitrite/nitrate concentrations. Chlorella and Monoraphidium also played an important role in phosphorus removal. For both Scenedesmus and Chlorella, smaller morphological types were more often associated with better system performance, whereas larger morphotypes likely associated with stress response(s) correlated with poor phosphorus recovery rates. Furthermore, these results demonstrate the potential of FIM as a critical technology for high-resolution characterization of industrial microalgal processes.

59 BASIC BIOLOGICAL SCIENCES↗

Development of Corynebacterium glutamicum as a monoterpene production platform

Monoterpenes are commonly known for their role in the flavors and fragrances industry and are also gaining attention for other uses like insect repellant and as potential renewable fuels for aviation. Corynebacterium glutamicum, a Generally Recognized as Safe microbe, has been a choice organism in industry for the annual million ton-scale bioproduction of amino acids for more than 50 years; however, efforts to produce monoterpenes in C. glutamicum have remained relatively limited. In this study, we report a further expansion of the C. glutamicum biosynthetic repertoire through the development and optimization of a mevalonate-based monoterpene platform. In the course of our plasmid design iterations, we increased flux through the mevalonate-based bypass pathway, measuring isoprenol production as a proxy for monoterpene precursor abundance and demonstrating the highest reported titers in C. glutamicum to date at 1504.6 mg/L. Our designs also evaluated the effects of backbone, promoter, and GPP synthase homolog origin on monoterpene product titers. Monoterpene production was further improved by disrupting competing pathways for isoprenoid precursor supply and by implementing a biphasic production system to prevent volatilization. With this platform, we achieved 321.1 mg/L of geranoids, 723.6 mg/L of 1,8-cineole, and 227.8 mg/L of linalool. Furthermore, we determined that C. glutamicum first oxidizes geraniol through an aldehyde intermediate before it is asymmetrically reduced to citronellol. Additionally, we demonstrate that the aldehyde reductase, AdhC, possesses additional substrate promiscuity for acyclic monoterpene aldehydes.

59 BASIC BIOLOGICAL SCIENCES↗

Nanoscopic lipid domains determined by microscopy and neutron scattering

Biological membranes are highly complex supramolecular assemblies, which play central roles in biology. However, their complexity makes them challenging to study their nanoscale structures. To overcome this challenge, model membranes assembled using reduced sets of membrane-associated biomolecules have been found to be both excellent and tractable proxies for biological membranes. Due to their relative simplicity, they have been studied using a range of biophysical characterization techniques. Here, in this review article, we will briefly detail the use of fluorescence and electron microscopies, and X-ray and neutron scattering techniques used over the past few decades to study the nanostructure of biological membranes.

59 BASIC BIOLOGICAL SCIENCES↗

Optimization of nonlinear turbulence in stellarators

We present new stellarator equilibria that have been optimized for reduced turbulent transport using nonlinear gyrokinetic simulations within the optimization loop. The optimization routine involves coupling the pseudo-spectral GPU-native gyrokinetic code GX with the stellarator equilibrium and optimization code DESC . Since using GX allows for fast nonlinear simulations, we directly optimize for reduced nonlinear heat fluxes. To handle the noisy heat flux traces returned by these simulations, we employ the simultaneous perturbation stochastic approximation (SPSA) method that only uses two objective function evaluations for a simple estimate of the gradient. We show several examples that optimize for both reduced heat fluxes and good quasi-symmetry as a proxy for low neoclassical transport. Finally, we run full transport simulations using the T3D stellarator transport code to evaluate the changes in the macroscopic profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

How does ion temperature gradient turbulence depend on magnetic geometry? Insights from data and machine learning

Magnetic geometry has a significant effect on the level of turbulent transport in fusion plasmas. Here, we model and analyse this dependence using multiple machine learning methods and a dataset of >200 000 nonlinear gyrokinetic simulations of ion-temperature-gradient turbulence in diverse non-axisymmetric geometries. The dataset is generated using a large collection of both optimised and randomly generated stellarator equilibria. At fixed gradients and other input parameters, the turbulent heat flux varies between geometries by several orders of magnitude. Trends are apparent among the configurations with particularly high or particularly low heat flux. Regression and classification techniques from machine learning are then applied to extract patterns in the dataset. Due to a symmetry of the gyrokinetic equation, the heat flux and regressions thereof should be invariant to translations of the raw features in the parallel coordinate, similar to translation invariance in computer vision applications. Multiple regression models including convolutional neural networks (CNNs) and decision trees can achieve reasonable predictive power for the heat flux in held-out test configurations, with highest accuracy for the CNNs. Using Spearman correlation, sequential feature selection and Shapley values to measure feature importance, it is consistently found that the most important geometric lever on the heat flux is the flux surface compression in regions of bad curvature. The second most important geometric feature relates to the magnitude of geodesic curvature. These two features align remarkably with surrogates that have been proposed based on theory, while the methods here allow a natural extension to more features for increased accuracy. The dataset, released with this publication, may also be used to test other proposed surrogates, and we find that many previously published proxies do correlate well with both the heat flux and stability boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗