Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SURROGATE MODELS”

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 1,621 records · Page 90

Coalescence of DNA Double Strand Breaks Induced by Galactic Cosmic Radiation is Modulated by Genetics in 15 Inbred Strains of Mice

In this manuscript we address the challenges associated with the ability to predict radiation sensitivity associated with exposure to either cosmic radiation or X-rays in a population study, by monitoring DNA damage sensing protein 53BP1 forming small nuclear radiation-induced foci (RIF) as a surrogate biomarker of DNA double strand breaks (DSB). 76 primary skin fibroblasts were isolated from 10 collaborative cross strains and five reference inbred mice (C57Bl/6, BALB/CByJ, B6C3, C3H and CBA/CaJ) and exposed to three different charged nuclei of increasing LET (350 MeV/n Si, 350 MeV/n Ar and 600 MeV/n Fe) and X-ray. Our data brings strong evidence against the classic "contact-first" model where DSBs are assumed to be immobile and repaired at the lesion site. In contrast, our model suggests nearby DSBs move into single repair unit characterized by large RIF before the repair machinery kicks in. Such model has the advantage of being much more efficient molecularly but is poorly suited to deal with cosmic radiation, where energy is concentrated along the particle trajectory, inducing a large density of DSBs along each particle track. In accordance with this model, RIF quantification after X-ray exposition showed a saturated dose response for early time points post-irradiation for all strains. Similarly, the high-LET response showed that RIF number matched the number of track per cell, not the number of expected DSB per cell (1). At the temporal level, we noted that the percentage of unrepaired high-LET tracks over a 48 hour time-course increased with LET, confirming that the DNA repair process becomes more difficult as more DSB coalesce into single RIF. There was also good agreement between persistent RIF levels measured in-vitro in the primary skin cultures and survival levels of T-cells and B-cells collected in blood samples from 10 CC strains 24 hours after 0.1 Gy whole-body dose of X-ray. This suggests that persistent RIF 24 hour post-IR is a good surrogate in-vitro biomarker for in-vivo radiation toxicity. Finally, at the genomic level, large differences in repair rates between strains for high-LET allowed us to identify suggestive genetic loci associated with radiation sensitivity. Interestingly, the two highest LETs provided the most strain variation with a common locus on Chromosome 10 highly enriched for DNA repair associated genes we discussed in detail.

Radiation-Induced Foci↗

Compression Analysis of Materials via in situ X-ray Computed Tomography

X-ray computed tomography (X-ray CT) is an analytical technique used in materials science to non– destructively characterize features in a variety of materials like polymers, metals, composites, and explosives. The non-destructive imaging allows for the analysis of features (voids and cracks), which give a fundamental understanding of material characteristics.This is more effective when combinedwith a load cell to expose the material to realworld stimuli to understand material characteristics and morphological behavior in situ. This study focused on using in situ X-ray CT to understand the material characteristics and morphological changes of sugar prills and a lattice under a compressive load. Sugar prills are used as a surrogate to analyze explosives in a safer and less expensiveway. The sugar prills were found to be brittle and shattered when compressed via the CT renderings and the stress/strain information. The lattices are used as both the basis and ground truth comparison for the development of models to predict the morphological and material characteristic changes that occur as a result of loading. The lattice was found to have weak regions due to low polymer interlayer fusion in these regions during the printing process. It was also found that under the initial compression the lattice demonstrated elastic behavior as its cellular structure was intact until the cellular structure started to collapse causing yielding in the lattice. Additionally, the analysis of the sugar prills and lattice demonstrate the information that can be gleaned from in situ X-ray CT experiments that would be lost in traditional pre- and post-mortem X-rayCT analysis.

36 MATERIALS SCIENCE↗

Active Learning for Directed Exploration of Complex Systems

Physics-based simulation codes are widely used in science and engineering to model complex systems that would be infeasible to study otherwise. Such codes provide the highest-fidelity representation of system behavior, but are often so slow to run that insight into the system is limited. For example, conducting an exhaustive sweep over a d-dimensional input parameter space with k-steps along each dimension requires k(sup d) simulation trials (translating into k(sup d) CPU-days for one of our current simulations). An alternative is directed exploration in which the next simulation trials are cleverly chosen at each step. Given the results of previous trials, supervised learning techniques (SVM, KDE, GP) are applied to build up simplified predictive models of system behavior. These models are then used within an active learning framework to identify the most valuable trials to run next. Several active learning strategies are examined including a recently-proposed information-theoretic approach. Performance is evaluated on a set of thirteen synthetic oracles, which serve as surrogates for the more expensive simulations and enable the experiments to be replicated by other researchers.

oracles↗

Effects of fuel composition and octane sensitivity on polycyclic aromatic hydrocarbon and soot emissions of gasoline–ethanol blend surrogates

The sooting propensity of a fuel is closely coupled with the fuel composition and chemistry. A detailed understanding of their effects is, therefore, needed to develop next-generation fuels which can minimize particulate emissions. With this overarching goal, the present work numerically investigates the effects of fuel composition and octane sensitivity (S) on polycyclic aromatic hydrocarbons (PAH) and soot emissions, for four-component gasoline-ethanol blend surrogates comprising isooctane, n-heptane, toluene, and ethanol. A partially-premixed counterflow flame is chosen as the canonical configuration for this study and simulations are performed using CHEMKIN-Pro-employing a kinetic mechanism developed by Park et al. (2017). In addition, a, detailed soot model based on the sectional method is used to capture the spatial characteristics of soot emissions. The kinetic mechanism and soot model are validated using available experimental data for various targets. A total of 86 TPRF-E mixtures, spanning a wide range of concentration of each component, and a wide range of S are analyzed. The effect of each non-paraffinic fuel component on the resultant PAH and soot emissions is investigated. PAH and soot emissions are found to vary significantly depending upon the blend composition. Additionally, based on the parametric sweeps, a regression analysis is carried out to identify global parameters that govern the formation of PAHs and soot. The analysis shows that both toluene content and S have a prominent effect on the formation of PAHs and soot, with toluene content having a stronger impact. Moreover, larger PAHs have higher dependency on toluene content and S. Furthermore, a detailed analysis is carried out to understand the physical and chemical phenomena associated with the observed trends of PAH and soot emissions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Benchmarking universal machine learning interatomic potentials for rapid analysis of inelastic neutron scattering data

The accurate calculation of phonons and vibrational spectra remains a significant challenge, requiring highly precise evaluations of interatomic forces. Traditional methods based on the quantum description of the electronic structure, while widely used, are computationally expensive and demand substantial expertise. Emerging universal machine learning interatomic potentials (uMLIPs) offer a transformative alternative by employing pre-trained neural network surrogates to predict interatomic forces directly from atomic coordinates. This approach dramatically reduces computation time and minimizes the need for technical knowledge. In this paper, we produce a phonon database comprising nearly 5000 inorganic crystals to benchmark the performance of several leading uMLIPs. We further assess these models in real-world applications by using them to analyze experimental inelastic neutron scattering data collected on a variety of materials. Through detailed comparisons, we identify the strengths and limitations of these uMLIPs, providing insights into their accuracy and suitability for fast calculations of phonons and related properties, as well as the potential for real-time interpretation of neutron scattering spectra. Our findings highlight how the rapid advancement of AI in science is revolutionizing experimental research and data analysis.

inelastic neutron scattering↗

Dakota, A Multilevel Parallel Object-Oriented Framework for Design Optimization Parameter Estimation Uncertainty Quantification and Sensitivity Analysis: Version 6.12 User's Manual

The Dakota toolkit provides a flexible and extensible interface between simulation codes and iterative analysis methods. Dakota contains algorithms for optimization with gradient and nongradient-based methods; uncertainty quantification with sampling, reliability, and stochastic expansion methods; parameter estimation with nonlinear least squares methods; and sensitivity/variance analysis with design of experiments and parameter study methods. These capabilities may be used on their own or as components within advanced strategies such as surrogate-based optimization, mixed integer nonlinear programming, or optimization under uncertainty. By employing object-oriented design to implement abstractions of the key components required for iterative systems analyses, the Dakota toolkit provides a flexible and extensible problem-solving environment for design and performance analysis of computational models on high performance computers. This report serves as a user's manual for the Dakota software and provides capability overviews and procedures for software execution, as well as a variety of example studies.

97 MATHEMATICS AND COMPUTING↗

Dakota A Multilevel Parallel Object-Oriented Framework for Design Optimization Parameter Estimation Uncertainty Quantification and Sensitivity Analysis (V.6.14) (User's Manual)

The Dakota toolkit provides a flexible and extensible interface between simulation codes and iterative analysis methods. Dakota contains algorithms for optimization with gradient and nongradient-based methods; uncertainty quantification with sampling, reliability, and stochastic expansion methods; parameter estimation with nonlinear least squares methods; and sensitivity/variance analysis with design of experiments and parameter study methods. These capabilities may be used on their own or as components within advanced strategies such as surrogate-based optimization, mixed integer nonlinear programming, or optimization under uncertainty. By employing object-oriented design to implement abstractions of the key components required for iterative systems analyses, the Dakota toolkit provides a flexible and extensible problem-solving environment for design and performance analysis of computational models on high performance computers. This report serves as a users manual for the Dakota software and provides capability overviews and procedures for software execution, as well as a variety of example studies.

97 MATHEMATICS AND COMPUTING↗

Dakota, A Multilevel Parallel Object-Oriented Framework for Design Optimization, Parameter Estimation, Uncertainty Quantification, and Sensitivity Analysis: Version 6.13 User's Manual

The Dakota toolkit provides a flexible and extensible interface between simulation codes and iterative analysis methods. Dakota contains algorithms for optimization with gradient and nongradient-based methods; uncertainty quantification with sampling, reliability, and stochastic expansion methods; parameter estimation with nonlinear least squares methods; and sensitivity/variance analysis with design of experiments and parameter study methods. These capabilities may be used on their own or as components within advanced strategies such as surrogate-based optimization, mixed integer nonlinear programming, or optimization under uncertainty. By employing object-oriented design to implement abstractions of the key components required for iterative systems analyses, the Dakota toolkit provides a flexible and extensible problem-solving environment for design and performance analysis of computational models on high performance computers. This report serves as a user’s manual for the Dakota software and provides capability overviews and procedures for software execution, as well as a variety of example studies.

97 MATHEMATICS AND COMPUTING↗

Dakota, A Multilevel Parallel Object-Oriented Framework for Design Optimization, Parameter Estimation, Uncertainty Quantification, and Sensitivity Analysis: Version 6.15 User's Manual

The Dakota toolkit provides a flexible and extensible interface between simulation codes and iterative analysis methods. Dakota contains algorithms for optimization with gradient and nongradient-based methods; uncertainty quantification with sampling, reliability, and stochastic expansion methods; parameter estimation with nonlinear least squares methods; and sensitivity/variance analysis with design of experiments and parameter study methods. These capabilities may be used on their own or as components within advanced strategies such as surrogate-based optimization, mixed integer nonlinear programming, or optimization under uncertainty. By employing object-oriented design to implement abstractions of the key components required for iterative systems analyses, the Dakota toolkit provides a flexible and extensible problem-solving environment for design and performance analysis of computational models on high performance computers.

97 MATHEMATICS AND COMPUTING↗

Dakota, A Multilevel Parallel Object-Oriented Framework for Design Optimization, Parameter Estimation, Uncertainty Quantification, and Sensitivity Analysis (V.6.16 User's Manual)

The Dakota toolkit provides a flexible and extensible interface between simulation codes and iterative analysis methods. Dakota contains algorithms for optimization with gradient and nongradient-based methods; uncertainty quantification with sampling, reliability, and stochastic expansion methods; parameter estimation with nonlinear least squares methods; and sensitivity/variance analysis with design of experiments and parameter study methods. These capabilities may be used on their own or as components within advanced strategies such as surrogate-based optimization, mixed integer nonlinear programming, or optimization under uncertainty. By employing object-oriented design to implement abstractions of the key components required for iterative systems analyses, the Dakota toolkit provides a flexible and extensible problem-solving environment for design and performance analysis of computational models on high performance computers. This report serves as a user's manual for the Dakota software and provides capability overviews and procedures for software execution, as well as a variety of example studies.

97 MATHEMATICS AND COMPUTING↗

Heart Pump Design for Cleveland Clinic Foundation

Through a Lewis CommTech Program project with the Cleveland Clinic Foundation, the NASA Lewis Research Center is playing a key role in the design and development of a permanently implantable, artificial heart pump assist device. Known as the Innovative Ventricular Assist System (IVAS), this device will take on the pumping role of the damaged left ventricle of the heart. The key part of the IVAS is a nonpulsatile (continuous flow) artificial heart pump with centrifugal impeller blades, driven by an electric motor. Lewis is part of an industry and academia team, led by the Ohio Aerospace Institute (OAI), that is working with the Cleveland Clinic Foundation to make IVAS a reality. This device has the potential to save tens of thousands of lives each year, since 80 percent of heart attack victims suffer irreversible damage to the left ventricle, the part of the heart that does most of the pumping. Impeller blade design codes and flow-modeling analytical codes will be used in the project. These codes were developed at Lewis for the aerospace industry but will be applicable to the IVAS design project. The analytical codes, which currently simulate the flow through the compressor and pump systems, will be used to simulate the flow within the blood pump in the artificial heart assist device. The Interdisciplinary Technology Office heads up Lewis' efforts in the IVAS project. With the aid of numerical modeling, the blood pump will address many design issues, including some fluid-dynamic design considerations that are unique to the properties of blood. Some of the issues that will be addressed in the design process include hemolysis, deposition, recirculation, pump efficiency, rotor thrust balance, and bearing lubrication. Optimum pumping system performance will be achieved by modeling all the interactions between the pump components. The interactions can be multidisciplinary and, therefore, are influenced not only by the fluid dynamics of adjacent components but also by thermal and structural effects. Lewis-developed flow-modeling codes to be used in the pump simulations will include a one-dimensional code and an incompressible three-dimensional Navier-Stokes flow code. These codes will analyze the prototype pump designed by the Cleveland Clinic Foundation. With an improved understanding of the flow phenomena within the prototype pump, design changes to improve the performance of the pump system can be verified by computer prior to fabrication in order to reduce risks. The use of Lewis flow modeling codes during the design and development process will improve pump system performance and reduce the number of prototypes built in the development phase. The first phase of the IVAS project is to fully develop the prototype in a laboratory environment that uses a water/glycerin mixture as the surrogate fluid to simulate blood. A later phase of the project will include testing in animals for final validation. Lewis will be involved in the IVAS project for 3 to 5 years.

Source record↗

Development of risk mitigation guidance for sensor placement inside mechanically ventilated enclosures – Phase 1

Guidance on Sensor Placement was identified as the top research priority for hydrogen sensors at the 2018 HySafe Research Priority Workshop on hydrogen safety in the category Mitigation, Sensors, Hazard Prevention, and Risk Reduction. This paper discusses the initial steps (Phase 1) to develop such guidance for mechanically ventilated enclosures. This work was initiated as an international collaborative effort to respond to emerging market needs related to the design and deployment equipment for hydrogen infrastructure that is often installed in individual equipment cabinets or ventilated enclosures. The ultimate objective of this effort is to develop guidance for an optimal sensor placement such that, when integrated into a facility design and operation, will allow earlier detection at lower levels of incipient leaks, leading to significant hazard reduction. Reliable and consistent early warning of hydrogen leaks will allow for the risk mitigation by reducing or even eliminating the probability of escalation of small leaks into large and uncontrolled events. To address this issue, a study of a real-world mechanically ventilated enclosure containing GH2 equipment was conducted, where CFD modeling of the hydrogen dispersion (performed by AVT and UQTR, and independently by the JRC) was validated by the NREL Sensor laboratory using a Hydrogen Wide Area Monitor (HyWAM) consisting of a 10-point gas and temperature measurement analyzer. In the release test, helium was used as a hydrogen surrogate. Expansion of indoor releases to other larger facilities (including parking structures, vehicle maintenance facilities and potentially tunnels) and incorporation into QRA tools, such as HyRAM is planned for Phase 2. It is anticipated that results of this work will be used to inform national and international standards such as NFPA 2 Hydrogen Technologies Code, Canadian Hydrogen Installation Code (CHIC) and relevant ISO/TC 197 and CEN documents.

08 HYDROGEN↗

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↗

ASK: Adversarial Soft k-Nearest Neighbor Attack and Defense

K-Nearest Neighbor (kNN)-based deep learning methods have been applied to many applications due to their simplicity and geometric interpretability. However, the robustness of kNN-based deep classification models has not been thoroughly explored and kNN attack strategies are underdeveloped. In this paper, we first propose an Adversarial Soft kNN (ASK) loss for developing more effective kNN-based deep neural network attack strategies and designing better defense methods against them. Our ASK loss provides a differentiable surrogate of the expected kNN classification error. It is also interpretable as it preserves the mutual information between the perturbed input and the in-class-reference data. We use the ASK loss to design a novel attack method called the ASK-Attack (ASK-Atk), which shows superior attack efficiency and accuracy degradation relative to previous kNN attacks on hidden layers. We then derive an ASK-Defense (ASK-Def) method that optimizes the worst-case ASK training loss. Experiments on CIFAR-10 (ImageNet) show that (i) ASK-Atk achieves ≥13% (≥ 13% ) improvement in attack success rate over previous kNN attacks, and (ii) ASK-Def outperforms the conventional adversarial training method by ≥ 6.9% (≥ 3.5% ) in terms of robustness improvement. Relevant codes are available at https://github.com/wangren09/ASK .

97 MATHEMATICS AND COMPUTING↗

Microcomputed X-Ray Tomographic Imaging and Image Processing for Microstructural Characterization of Explosives

Microstructural characterization of composite high explosives (HEs) has become increasingly important over the last several decades in association with the development of high fidelity mesoscale modeling and an improved understanding of ignition and detonation processes. HE microstructure influences not only typical material properties (e.g., thermal, mechanical) but also reactive behavior (e.g., shock sensitivity, detonation wave shape). A detailed nondestructive 3D examination of the microstructure has generally been limited to custom-engineered samples or surrogates due to poor contrast between the composite constituents. Highly loaded (>90 wt%) HE composites such as plastic-bonded explosives (PBX) are especially difficult. Here, we present efforts to improve measurement quality by using single and dual-energy microcomputed X-ray tomography and state-of-the-art image processing techniques to study a broad set of HE materials. Some materials, such as PBX 9502, exhibit suitable contrast and resolution for an automatic segmentation of the HE from the polymer binder and the voids. Other composite HEs had varying levels of success in segmentation. Post-processing techniques that used commercially available algorithms to improve the segmentation quality of PBX 9501 as well as zero-density defects such as cracks and voids could be easily segmented for all samples. Aspects of the materials that lend themselves well to this type of measurement are discussed.

36 MATERIALS SCIENCE↗

Simulation-based inference for parameter estimation of complex watershed simulators

High-resolution, spatially distributed process-based (PB) simulators are widely employed in the study of complex catchment processes and their responses to a changing climate. However, calibrating these PB simulators using observed data remains a significant challenge due to several persistent issues, including the following: (1) intractability stemming from the computational demands and complex responses of simulators, which renders infeasible calculation of the conditional probability of parameters and data, and (2) uncertainty stemming from the choice of simplified representations of complex natural hydrologic processes. Here, we demonstrate how simulation-based inference (SBI) can help address both of these challenges with respect to parameter estimation. SBI uses a learned mapping between the parameter space and observed data to estimate parameters for the generation of calibrated simulations. To demonstrate the potential of SBI in hydrologic modeling, we conduct a set of synthetic experiments to infer two common physical parameters – Manning's coefficient and hydraulic conductivity – using a representation of a snowmelt-dominated catchment in Colorado, USA. We introduce novel deep-learning (DL) components to the SBI approach, including an “emulator” as a surrogate for the PB simulator to rapidly explore parameter responses. We also employ a density-based neural network to represent the joint probability of parameters and data without strong assumptions about its functional form. While addressing intractability, we also show that, if the simulator does not represent the system under study well enough, SBI can yield unreliable parameter estimates. Approaches to adopting the SBI framework for cases in which multiple simulator(s) may be adequate are introduced using a performance-weighting approach. The synthetic experiments presented here test the performance of SBI, using the relationship between the surrogate and PB simulators as a proxy for the real case.

54 ENVIRONMENTAL SCIENCES↗

Multi-stage heat release of multi-component fuels: Insights and implications for advanced engine operation

Multi-stage heat release (MSHR) is a unique phenomenon typically seen in lean/diluted and low- to intermediate-temperature combustion. Despite its relevance to advanced engine operation, the MSHR of multi-component fuels has barely been quantified. This study aims to characterize the MSHR of multi-component fuels in a rapid compression machine (RCM) at conditions representative of advanced combustion engines. New experimental data are first reported in the RCM at an equivalence ratio of 0.4, pressure of 60 bar and temperatures from 702 to 795 K for a research grade, multi-component gasoline surrogate, termed PACE-20. Here, experiments confirm the existence of MSHR for PACE-20 at all temperatures, and reveal the strong inhibiting effect of temperature on MSHR, where increasing temperature inhibits MSHR by suppressing first stage heat release and promoting second and third stage heat release. A detailed chemical kinetic model is also adopted to model the experiments, with good agreement observed. The response of MSHR characteristics to changes in different engine operating parameters (i.e., temperature, pressure, equivalence ratio and CO 2 dilution level) and fuel compositions (i.e., mole fraction of n-pentane, n-heptane, isooctane, cyclopentane, 1-hexene, ethanol and aromatics in PACE-20) is further evaluated via statistical analysis coupling quasi-random sampling, extensive computer experiments and high-dimensional model representation. The change in MSHR is characterized through 8 quantities of interest (QOIs), i.e., duration and extent of each heat release stage, and their standard deviations. The analysis highlights the dependence of the QOIs on the individual parameters as well as their interplays, with temperature exhibiting the strongest impact among all parameters. Furthermore, it is demonstrated that if MSHR is to be facilitated and combustion phasing is to be extended to enable engine operation at higher compression ratios, it is recommended to operate the engine at low intake temperatures, boosted intake pressures and high CO 2 dilutions (i.e., high EGR levels) with gasoline fuels containing more n-pentane and less aromatics.

42 ENGINEERING↗

Shock Tube and Flame Speed Measurements of 2,4,4-Trimethyl-1-Pentene: A Co-Optima Biofuel

Abstract The combustion of 2,4,4-trimethyl-1-pentene (diisobutylene, C8H16), which is a biofuel and a component of surrogate fuels, is examined in this work. Carbon monoxide time–histories and ignition delay times are collected behind reflected shock waves utilizing a shock tube and mid-infrared laser absorption spectroscopy. Measurements were obtained near 10 atm pressure during stoichiometric oxidation of 0.15%C8H16/O2/Ar. Simulated results from chemical kinetic models are provided, and sensitivity analyses are used to discuss differences between models for both ignition delay times and carbon monoxide formation. In addition, laminar burning speeds are obtained at 1 atm, 428 K, and equivalence ratios, phi, between 0.91 and 1.52 inside a spherical chamber facility. Measured burning speeds are found to be less than that of ethanol over the equivalence ratio span. Burning speed measurements are compared to predictions of chemical kinetic mechanisms and are in agreement for the richest conditions; however, at lean conditions, the model predicts a far slower-burning speed. The maximum burning speed occurs at an equivalence ratio of 1.08 with a magnitude of 0.70 m/s. The current work provides the crucial experimental data needed for assessing the feasibility of this biofuel and for the development of future combustion chemical kinetics models.

Energy & Fuels↗