Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forward modeling”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Electron emission from fluorene (C 13 H 10 ) upon 3.5 MeV/u Si 8+ ion impact: double differential distributions

In this work, we have studied the electron emission from one of the polycyclic aromatic hydrocarbon (PAH) molecules namely, fluorene (C13H10), upon 3.5 MeV/u Si8+ ion impact. The experimentally measured absolute double differential cross sections (DDCS) are compared with the continuum distorted wave-eikonal initial state (CDW-EIS) model and the first Born approximation including correct boundary conditions (CB1). The measurements are carried out in the ejected e⁻-energy range of 1 eV–400 eV and in the angular range of 20°–160°. We have obtained the single differential and the total cross sections (TCSs) of e⁻-emission as well. The CB1 calculation largely underestimates the data. The CDW-EIS model, which is applied for the PAH molecule for the first time, provides an overall better agreement with the double differential, single differential and TCS data. The DDCS data for fluorene has also been compared with that for CH4 molecule, at a few angles. The forward–backward angular asymmetry shows a relatively flatter distribution compared to the theoretical predictions. The contribution due to the giant plasmon resonance could not be clearly observed except a mild indication in the asymmetry parameter. The angular distribution of the carbon KLL Auger electron cross section shows certain variations. The study of the KLL hyper-satellite component indicates the double K-ionization cross section is about 8.6% of the single K-ionization one.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Insights into the relation between vertical cloud structure and dynamics of three polar lows: Observations from COMBLE

The vertical structure of three polar lows is described using profiling radar, lidar, and passive remote sensors deployed at a coastal site in northern Norway. The data were collected as part of the Cold-air Outbreaks in the Marine Boundary Layer Experiment (COMBLE). These data are examined in the context of satellite imagery, operational weather radar reflectivity imagery, and output from the AROME-Arctic model. A comparison to satellite images shows that AROME-Arctic captured the polar lows well. All three polar lows form in marine cold-air outbreaks and are surrounded by open-cellular convection typical of such outbreaks. Two of the lows form in a forward shear environment, and one under reverse shear. Initially, the three polar lows are comma-shaped; two of them transition to be more spirali-form at their mature stage and have mostly cloud-free warm cores. The warm cores are the result of a warm-seclusion process. All three lows have stratiform precipitation bands, marked by little cloud liquid water, and rather high surface precipitation rate. The vertical drafts and turbulence in these stratiform clouds are generally weak. In conclusion, all three lows also have convective clouds, which have stronger vertical drafts, stronger turbulence, and pockets of high liquid water content.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Optimal control to handle variations in moisture content and reactor in-feed rate

The variations in feedstock characteristics, such as moisture and particle size distribution, lead to an inconsistent flow of feedstock from the biomass pre-processing system to the reactor in-feed system. These inconsistencies result in low on-stream times at the reactor in-feed equipment. This research develops an optimal process control method for a biomass pre-processing system comprised of milling and densification operations to provide the consistent flow of feedstock to a reactor's throat. This method uses a mixed-integer optimization model to identify optimal bale sequencing, equipment in-feed rate, and buffer location and size in the biomass pre-processing system. This method, referred to as the hybrid process control (HPC), aims to maximize throughput over time. We compare HPC with a baseline feed forward process control. Our case study based on switchgrass finds that HPC reduces the variation of a reactor's feeding rate by up to 100% without increasing the operating cost of the biomass pre-processing system for biomass with moisture ranging from 10 to 25%. Additionally, HPC reduces the cost of processing biomass by 0.36%–2.22%, and reduces processing time by 0.35%-2.24%. Furthermore, a biorefinery can adapt HPC to achieve its design capacity.

09 BIOMASS FUELS↗
DOE OSTI2021

Role of deep levels and barrier height lowering in current-flow mechanism in 150 μ m thick epitaxial n-type 4H–SiC Schottky barrier radiation detectors

Schottky barrier detectors (SBDs) require larger surface areas than conventional electronics to increase the detection efficiency although such SBDs manifest large diode ideality factors due to inhomogeneous areal distribution of surface barrier height (SBH). Inhomogeneous SBH distributions lead to various current flow mechanisms in SBDs, which need to be identified to optimize detector performance. In this Letter, we identify the current flow mechanism in large area Schottky barrier diodes for radiation detection fabricated on 150 μm thick n-4H–SiC epitaxial layers. The analysis of temperature-dependent forward current–voltage (I–V–T) characteristics of SBDs revealed two linear regions in current–voltage curves up to 450 K, one corresponding to the current flow through a low barrier patch, while the other corresponds to that of average barrier distribution. Applying a SBH distribution model to the reverse I–V–T characteristics, an activation energy of 0.76 eV for the current flow over the Schottky barrier was calculated. Additionally, the activation energy did not directly correspond to any of the defect levels observed from the deep level transient spectroscopy (DLTS). Above 450 K, a Schottky type barrier lowering suggested a current flow through a low barrier patch of ≈ 0.8 eV. The absence of any SBH lowering below 450 K indicated that the current corresponded to a neutrally charged trap level at ≈ 0.6 eV below the conduction band edge, which was consistent with DLTS measurements revealing the presence of an electron trap level Z1/2 at 0.59 eV below the conduction band edge.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Toward model-free temperature diagnostics of warm dense matter from multiple scattering angles

Warm dense matter plays an important role in astrophysical objects and technological applications, but the rigorous diagnostics of corresponding experiments is notoriously difficult. In this work, we present a model-free analysis of x-ray Thomson scattering (XRTS) measurements on isochorically heated graphite obtained at the Linac Coherent Light Source at multiple scattering angles. We demonstrate that the recent imaginary-time thermometry technique works for scattering data that have been measured in both forward and backward scattering geometry. This opens up the way toward a rigorous quantification of nonequilibrium effects in future experiments, where XRTS measurements are being obtained from multiple scattering angles from the same sample.

Equations of state↗
DOE OSTI2020

A New Capability of E4D For 3D Parallel Joint Inversion of DC Resistivity And Traveltime Data on Unstructured Mesh

A major challenge in interpreting geophysical data is how to derive consistent three-dimensional (3D) earth models of different physical properties from spatially and temporally limited measurements. Joint inversion with cross-gradient constraints is an approach to find such models by imposing structural similarities between different physical parameters. We have developed a parallel distributed-memory joint inversion code for direct-current (DC) resistivity and traveltime data using the cross-gradient constraint on unstructured mesh. The code utilizes existing E4D framework for parallel forward simulation, distributed storage and computation of the Jacobian matrix of forward operator, and parallel execution of matrix-vector multiplication during inversion. Besides, the joint inversion is solved by nonlinear conjugate gradient algorithm parallelized for DC resistivity and traveltime data. The joint inversion capability of E4D was tested using synthetic data from cross-borehole DC resistivity and traveltime data. The results indicate that the shape and size of the anomalies from the joint inversion are more reliable than those from separate inversions.

58 GEOSCIENCES↗
DOE OSTI2020

Validation of SAS4A/SASSYS-1 for Steady-State Single-Phase Natural Circulation

The development of numerical software for engineering applications requires the validation of the code against experimental benchmark datasets. To support the validation of SAS4A/SASSYS-1 in simulating the single-phase natural circulation within nuclear reactor systems, high-precision steady-state experiments are performed to capture single-phase natural circulation phenomena on an existing scaled facility with comprehensive instrumentation. Forced convection tests are designed precisely capturing the facility’s critical thermal-hydraulic parameters required for one-dimensional modeling, and a comprehensive single-phase natural circulation dataset is obtained with well-documented experimental uncertainty and facility description. Analysis and discussion based on the finalized dataset confirm the dataset’s ability in capturing dominant phenomena and address important physical interpretation of parameters. The dataset reported in this project provides a valuable benchmark resource for the validation of system analysis codes under single-phase natural circulation inside nuclear reactors. Validation of SAS4A/SASSYS-1 is then performed against the obtained benchmark dataset to confirm the capability of its physics model in simulating the steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators to quantify and restrict numerical errors from discretization. Numerical uncertainty originating from finite maximum pseudo-transient time is also quantified and restricted. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, directly quantifying the simulation output in a Monte Carlo approach. The composition of the output uncertainty is further determined by the estimators of Sobol’ indices through a variance-based sensitivity analysis. With uncertainty quantified for each individual condition, detailed comparison between the simulation results and experimental data is performed covering the whole dataset, which shows consistent agreement for all important quantities. The validation activity reported in this project demonstrates that SAS4A/SASSYS-1 can predict the primary parameters with satisfactorily accuracy under steady-state single-phase natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Uncertainty-quantified phenomenological optical potentials for single-nucleon scattering

Optical-model potentials (OMPs) continue to play a key role in nuclear reaction calculations. However, the uncertainty of phenomenological OMPs in widespread use—inherent to any parametric model trained on data—has not been fully characterized, and its impact on downstream users of OMPs remains unclear. Here we assign well-calibrated uncertainties for two representative global OMPs, those of Koning-Delaroche and Chapel Hill '89, using Markov-chain Monte Carlo for parameter inference. By comparing the canonical versions of these OMPs against the experimental data originally used to constrain them, we show how a lack of outlier rejection and a systematic underestimation of experimental uncertainties contributes to bias of, and overconfidence in, best-fit parameter values. Our updated, uncertainty-quantified versions of these OMPs address these issues and yield complete covariance information for potential parameters. Scattering predictions generated from our ensembles show improved performance both against the original training corpora of experimental data and against a new “test” corpus comprising many of the experimental single-nucleon scattering data collected over the last twenty years. Finally, we apply our uncertainty-quantified OMPs to two case studies of application-relevant cross sections. We conclude that, for many common applications of OMPs, including OMP uncertainty should become standard practice. Furthermore, to facilitate their immediate use, digital versions of our updated OMPs and related tools for forward uncertainty propagation are included as Supplemental Material.

150 ≤ A ≤ 189↗
DOE OSTI2023

Electron scattering and neutrino physics

Abstract A thorough understanding of neutrino–nucleus scattering physics is crucial for the successful execution of the entire US neutrino physics program. Neutrino–nucleus interaction constitutes one of the biggest systematic uncertainties in neutrino experiments—both at intermediate energies affecting long-baseline deep underground neutrino experiment, as well as at low energies affecting coherent scattering neutrino program—and could well be the difference between achieving or missing discovery level precision. To this end, electron–nucleus scattering experiments provide vital information to test, assess and validate different nuclear models and event generators intended to test, assess and validate different nuclear models and event generators intended to be used in neutrino experiments. Similarly, for the low-energy neutrino program revolving around the coherent elastic neutrino–nucleus scattering (CEvNS) physics at stopped pion sources, such as at ORNL, the main source of uncertainty in the evaluation of the CEvNS cross section is driven by the underlying nuclear structure, embedded in the weak form factor, of the target nucleus. To this end, parity-violating electron scattering (PVES) experiments, utilizing polarized electron beams, provide vital model-independent information in determining weak form factors. This information is vital in achieving a percent level precision needed to disentangle new physics signals from the standard model expected CEvNS rate. In this white paper, we highlight connections between electron- and neutrino–nucleus scattering physics at energies ranging from 10 s of MeV to a few GeV, review the status of ongoing and planned electron scattering experiments, identify gaps, and lay out a path forward that benefits the neutrino community. We also highlight the systemic challenges with respect to the divide between the nuclear and high-energy physics communities and funding that presents additional hurdles in mobilizing these connections to the benefit of neutrino programs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

1.2-kV Vertical GaN PIN Rectifier With Ion-Implanted Floating Guard Rings

Here, this work reports on the fabrication and properties of a homojunction gallium nitride (GaN) p-i-n (PIN) rectifier fabricated on a free-standing GaN substrate. Uniform device performance is achieved with breakdown voltage (BV) >1.2 kV and low ON-resistance × area ( R ON A ). The statistics of the BV measurements show 58.5% of devices achieve BV >1.3 kV, and 71.1% of devices achieve BV >1.2 kV, as attributed to high quality and control in both epitaxial growth and device process. At room temperature, R ON A is 0.23 mΩ∙ cm 2 at a current density ( J ) of 6.9 kA/cm2. The corresponding Baliga’s figure of merit is >5.97 GW/cm 2 . Temperature-dependent reverse I – V measurements were performed and show a positive temperature coefficient of 0.42 V/K, indicating the avalanche capability of reverse breakdown. Further analysis with the Poole-Frenkel model on the temperature-dependent measurement suggested that a trap-assisted tunneling process contributed to the reverse leakage current. Floating guard rings (FGRs) formed by nitrogen implantation serve as an effective edge termination technique in these GaN PIN rectifiers, resulting in uniform performance in both forward and reverse bias.

42 ENGINEERING↗
DOE OSTI2020

Computational analysis of deployable wind turbine systems in defense operational energy applications.

The U.S. military has been exploring pathways to reduce the logistical burden of fuel on virtually all their missions globally. Energy harvesting of local resources such as wind and solar can help increase the resilience and operational effectiveness of military units, especially at the most forward operating bases where the fuel logistics are most challenging. This report considers the potential benefits of wind energy provided by deployable wind turbines as measured by a reduction in fuel consumption and supply convoys to a hypothetical network of Army Infantry Brigade Combat Team bases. Two modeling and simulation tools are used to represent the bases and their operations and quantify the impacts of system design variables that include wind turbine technologies, battery storage, number of turbines, and wind resource quality. The System of Systems Analysis Toolkit Joint Operational Energy Model serves as a baseline scenario for comparison. The Hybrid Optimization of Multiple Energy Resources simulation tool is used to optimize a single base within the larger Joint Operational Energy Model. The results of both tools show that wind turbines can provide significant benefits to contingency bases in terms of reduced fuel use and number of convoy trips to resupply the base. The match between the turbine design and wind resource, which is statistically low across most of the global land area, is a critical design consideration. The addition of battery storage can enhance the benefits of wind turbines, especially in systems with more wind turbines and higher wind resources. Wind turbines may also provide additional benefits to other metrics such as resilience that may be important but not fully considered in the current analysis. ACKNOWLEDGEMENTS The authors would like to thank the following individuals for their helpful support, feedback and review to improve this report: U.S. Department of Energy Wind Energy Technologies Office, Patrick Gilman and Bret Barker; Idaho National Laboratory, Jake Gentle and Bradley Whipple; The National Renewable Energy Laboratory, Robert Preus and Tony Jimenez; Sandia National Laboratories, Alan Nanco, Dennis Anderson, and Hai Le. In addition, numerous discussions with military and industry stakeholders over the year were invaluable in focusing the efforts represented in this report.

17 WIND ENERGY↗
DOE OSTI2022

Modified Calcium Manganites for Thermochemical Energy Storage Applications

CaAl 0.2 Mn 0.8 O 3-δ (CAM28) and CaTi 0.2 Mn 0.8 O 3-δ (CTM28) are perovskite metal oxides developed for high-temperature thermochemical energy storage (TCES) applications, e.g., in support of air Brayton power generation. Previous reports for these compounds focus on the equilibrium non-stoichiometry (δ) as a function of temperature and oxygen partial pressure (pO 2 ) and the endotherm (or exotherm) accompanying changes in δ resulting from thermal reduction (or re-oxidation). Herein, we report results for elemental substitution and doping (Al, Co, Fe, La, Sr, Ti, Y, Zn, and Zr) of calcium manganites (CM) that establish the preference for CAM28 and CTM28. Techniques employed include conventional (screening and equilibrium) and ballistically heated multi-cycle thermogravimetric analysis (TGA), conventional and high temperature ( in-situ ) X-ray diffraction (XRD), and differential scanning calorimetry (DSC). Forward-looking results for A-site Y-doped materials, e.g., Ca 0.9 Y 0.1 MnO 3-δ (CYM910), establish a route to increasing the reduction enthalpy relative to CAM28 and CTM28, albeit at the expense of increased reduction temperatures and raw materials costs. A thermodynamic model presented for CAM28, but extendable to related materials, provides values for the reaction enthalpy and extent of reduction as a function of temperature and oxygen partial pressure for use in design efforts. Taken as a whole, the results support the choice of Al-doped CaMnO 3-δ as a low-cost material for TCES in a high temperature air Brayton application, but point the way to achieving higher stored energy densities that could lead to overall cost savings.

Miller, James E.↗
DOE OSTI2024

Promoting the regulatory acceptance of combined ion and neutron irradiation for material degradation in nuclear reactors

The Advanced Materials and Manufacturing Technologies (AMMT) program within the Department of Energy (DOE) Office of Nuclear Energy has developed its current recommendation for promoting the use of combined ion irradiation and neutron irradiation for the accelerated qualification of materials to be deployed in nuclear reactors. This plan is intended to provide a collaborative path forward that can be adopted by academia, national laboratories, and industry, and has been developed with input from the regulatory research arm of the U.S. Nuclear Regulatory Commission (NRC). To deploy new materials or materials manufactured with new technologies, such as additive manufacturing, materials must be evaluated for reactor-induced degradation from the combination of harsh temperatures, corrosive environments, and radiation fields. However, rapid deployment of materials necessitates accelerated testing methods rather than relying on years of neutron irradiation in a material test reactor. Ion irradiation has demonstrated success in reproducing material microstructure and select property evolution resulting from neutron irradiation with three to four orders of magnitude reduction in time and cost, making it an ideal candidate for accelerated irradiation testing. This presentation provides context governing both the scientific and regulatory aspects of the proposed goal. The discussion is aimed at a broad audience including researchers from industry, national laboratories, and academia. The recommended path forward is presented as a conceptual framework of specific steps. In brief, the strategy entails developing an integrated ion and neutron irradiation test plan for the material property of interest based on the fundamental tenet of the linkage of microstructure and properties in materials. Physics-based modeling interprets ion irradiation data and predicts neutron irradiation microstructure and properties with uncertainty bounds. The first round of testing is sufficient for an initial licensing application using a risk-informed approach, while a minimum required neutron irradiation test plan reduces cost and time requirements. A surveillance program with witness specimens in-reactor provides additional data over time to improve model predictions to higher damage levels and further reduce uncertainty bounds, which can be used for license extensions or longer lifetimes in new license applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Complexity growth in integrable and chaotic models

We use the SYK family of models with N Majorana fermions to study the complexity of time evolution, formulated as the shortest geodesic length on the unitary group manifold between the identity and the time evolution operator, in free, integrable, and chaotic systems. Initially, the shortest geodesic follows the time evolution trajectory, and hence complexity grows linearly in time. We study how this linear growth is eventually truncated by the appearance and accumulation of conjugate points, which signal the presence of shorter geodesics intersecting the time evolution trajectory. By explicitly locating such “shortcuts” through analytical and numerical methods, we demonstrate that: (a) in the free theory, time evolution encounters conjugate points at a polynomial time; consequently complexity growth truncates at O($\sqrt{N}$), and we find an explicit operator which “fast-forwards” the free N-fermion time evolution with this complexity, (b) in a class of interacting integrable theories, the complexity is upper bounded by O(poly(N)), and (c) in chaotic theories, we argue that conjugate points do not occur until exponential times O(e N ), after which it becomes possible to find infinitesimally nearby geodesics which approximate the time evolution operator. Finally, we explore the notion of eigenstate complexity in free, integrable, and chaotic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Coupling High-Throughput Experiments and Regression Algorithms to Optimize PGM-Free ORR Electrocatalyst Synthesis

Over the past decades, significant improvement has been achieved in the performance of platinum group metal-free (PGM-free) materials as an alternative to Pt-based electrocatalysts for oxygen reduction reaction (ORR). However, further progress in ORR activity requires evaluation of precursors and synthesis approaches. In response to this challenge, we generated a first of its kind experimental data set of 36 samples using high-throughput synthesis and activity measurements. Several control parameters (e.g., Fe precursor identity, the precursor content, and pyrolysis temperature) were varied. We then developed several state-of-the-art machine learning (ML) based regression models to predict ORR activity, dependent on selected synthesis variables. Through an iterative algorithm, higher prediction accuracy (smaller root-mean-square error) was achieved. We identified that gradient boosting regression (GBR) and support vector regression (SVR), among several methods, work best for this data set. Aided by our ML-based surrogate models, we decided to alter catalyst synthesis conditions, which resulted in a 36% increase in measured ORR activity in comparison to the maximum ORR mass activity value of 21.9 A/g catalyst in the original data set. Overall, this combined experiment and machine learning approach represents a promising path forward toward developing highly efficient next-generation ORR electrocatalysts and, more generally, functional materials.

25 ENERGY STORAGE↗
DOE OSTI2022

Modeling and experimental validation of dynamical effects in Bragg coherent x-ray diffractive imaging of finite crystals

Bragg coherent diffractive imaging (BCDI) is a noninvasive microscopy technique that can visualize the shape and internal lattice deviations of crystals with nanoscale spatial resolution and picometer deformation sensitivity. Its strain imaging capability relies on Fourier transform–based iterative phase retrieval algorithms, which are mostly developed under the kinematical approximation. Such approximation prohibits the application of BCDI on larger crystals, which are commonly seen in most emerging functional materials. Understanding the dynamical effect in BCDI, as well as developing a validated method for modeling BCDI at the dynamical diffraction limit, is crucial for applying BCDI to hierarchical systems that contain micron-sized crystals and grains. Thus we report a comparative study on the impact of dynamical diffraction effects by comparing the reconstruction results from two measurements of the same crystal. Forward simulation is implemented to show subtle changes of interference fringes in the diffraction pattern due to the dynamical diffraction, and is compared directly with the experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2021

EXPLAINABLE AND TRUSTWORTHY DIAGNOSTICS ACHIEVABLE THROUGH PROCESS-BASED AUTOMATED REASONING

An approach has been developed that incorporates domain knowledge to obtain a more explainable and trustworthy equipment health monitoring diagnosis than might otherwise be obtained from a purely data-driven method. Physics-based models serve to constrain the realizable solution space and render a more trusted diagnosis. An automated reasoning algorithm performs backward chaining to infer a diagnosis that is consistent with logic statements that have been evaluated as true. This diagnosis is made explainable to an operator by providing the forward chaining path that elucidates for inspection and validity testing those truths implied by the diagnosis.

automated reasoning↗
DOE OSTI2021

Multisensor Agile Adaptive Sampling of Convective Storms Driven by Real-time Analytics

Convective storms vertically transport water vapor and condensate from Earth’s surface to the upper troposphere. Life on Earth is fundamentally linked to this transport which determines the hydrological cycle, and the intensity of severe weather responsible for the destruction of life and property. Despite advances in high-resolution modeling and better observational capabilities, the scientific community continues to be confronted with knowledge gaps about convective storms that limit our predictive capabilities. The ongoing developments in the high-resolution Energy Exascale Earth System Model (E3SM), large eddy simulations, and AI-based analytics to evaluate uncertainties are expected to provide a comprehensive framework for new scientific discovery. The model-experiment (MODEX) approach suggests that the aforementioned advancements in model development and AI-based inference techniques should be complemented by similar advancements in the experimental (observational) side so that the former does not outstrip the ability of the latter to provide meaningful constraints. What are the recent advancements in observations that will provide the necessary leap forward in improving our predictive capabilities? To address this question, we propose a new experimental paradigm called Multisensor Agile Adaptive Sampling (MAAS) that capitalizes on advancements in communications (5G), computational resources (edge/fog computing), sensor capabilities, and machine learning (ML) and AI techniques (Kollias et al., 2020). The MAAS framework allows for the collection of higher spatiotemporal resolution and quality observations of convective storms than is traditionally possible. The MAAS framework is scalable and applicable to atmospheric observatories such as those operated by the Department of Energy (DoE) Atmospheric Radiation Measurement (ARM) facility.

54 ENVIRONMENTAL SCIENCES↗