Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transient analysis”

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,045 records · Page 58

Discrete modeling of a transformer with ALEGRA

We report progress on a task to model transformers in ALEGRA using the “Transient Magnetics” option. We specifically evaluate limits of the approach resolving individual coil wires. There are practical limits to the number of turns in a coil that can be numerically modeled, but calculated inductance can be scaled to the correct number of turns in a simple way. Our testing essentially confirmed this “turns scaling” hypothesis. We developed a conceptual transformer design, representative of practical designs of interest, and that focused our analysis. That design includes three coils wrapped around a rectangular ferromagnetic core. The secondary and tertiary coils have multiple layers. The tertiary has three layers of 13 turns each; the secondary has five layers of 44 turns; the primary has one layer of 20 turns. We validated the turns scaling of inductance for simple (one-layer) coils in air (no core) by comparison to available independent calculations for simple rectangular coils. These comparisons quantified the errors versus reduced number of turns modeled. For more than 3 turns, the errors are <5%. The magnetic field solver failed to converge (within 5000 iterations) for >10 turns. Including the core introduced some complications. It was necessary to capture the core surfaces in thin grid sheaths to minimize errors in computed magnetic energy. We do not yet have quantitative benchmarks with which to compare, but calculated results are qualitatively reasonable.

42 ENGINEERING↗

Coherent Phonon Dynamics in Plasmonic Gold Tetrahedral Nanoparticle Ensembles

Coherent phonon modes supported by plasmonic nanoparticles offer prospective applications in chemical and biological sensing. Whereas characterization of these phonon modes often requires single particle measurements, synthetic routes to narrow size distributions of nanoparticles permit ensemble investigations. Recently, the synthesis of highly monodisperse gold tetrahedral nanoparticles with tunable edge length and corner sharpness has been developed. Herein, we characterize a size series of these nanoparticles in colloidal dispersion via transient absorption spectroscopy to examine their mechanical and plasmonic responses upon photoexcitation. Oscillations of transient absorption signals are observed in the plasmon resonance and correspond to the lowest order radial breathing modes of the nanoparticles, the frequencies of which are affected by edge length and truncation of the corners. Furthermore, quality factor values ranging from 10 to 15 are observed for the oscillations that convey potential utility in mass-sensing and plasmon-exciton coupling photonics schemes. Finite-difference time domain and finite element analysis calculations establish the specific optically-relevant phonon modes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tracking carrier and exciton dynamics in mixed-cation lead mixed-halide perovskite thin films

Mixed-cation lead mixed-halide perovskites simultaneously possess structural stability and high power conversion efficiency. A thorough study of both carrier and exciton dynamics is needed to understand the photophysical properties that underpin its superior photovoltaic performance. By utilizing a broadband transient absorption spectroscopy, we observe the carrier and exciton dynamics in a FA 0.85 Cs 0.15 Pb(I 0.97 Br 0.03 ) 3 (FCPIB) perovskite by simultaneously resolving the carrier and exciton contribution to the transient change of the absorption spectra, from which the carrier density and exciton oscillator strength can be determined. Our data reveal a quick and significant conversion of the photogenerated carriers to excitons, on top of the usual carrier recombination process. Moreover, the decay of carrier density shows a change of kinetics from a second-order recombination at high pump fluence to a third-order recombination at low pump fluence. Our analysis utilizes band anharmonicity, presents an independent determination of electronic temperature and quasi-Fermi energy, and reveals an interesting interplay among the processes of carrier cooling, exciton formation/decay and carrier recombination, all as a function of time after photoexcitation. Our work demonstrates the use of pump fluence as a knob to tune the relative populations of carriers and excitons in halide perovskite materials.

36 MATERIALS SCIENCE↗

Growing neoclassical tearing modes seeded via transient-induced-multimode interactions

Nonlinear extended magnetohydrodynamics simulations demonstrating seeding of neoclassical tearing modes (NTMs) via magnetohydrodynamics-transient-induced multimode interactions are presented. Simulations of NTMs are enabled by two recent NIMROD code developments: the implementation of heuristic neoclassical stresses and the application of transient magnetic perturbations (MPs) at the boundary. NTMs are driven unstable by the inherently pressure driven kinetic bootstrap current, which arises due to collisional viscosity between passing and trapped electrons. These simulations use heuristic closures that model the neoclassical electron and ion stresses. NTM growth requires a seed island, which is generated by a transiently applied MP in simulations. The capability is demonstrated using kinetic-based reconstructions with flow of a DIII-D ITER Baseline Scenario discharge (La Haye et al., in Proceedings IAEA FEC, 2020). The applied MP seeds a 2/1 NTM that grows in two phases: a slow growth phase followed by a faster robust growth phase. Additionally, an evolving sequence of higher order core modes are excited at first. Power transfer analysis shows that nonlinear interactions between the core modes and the 2/1 helps drive the initial slow growth. Once the induced 2/1 magnetic island reaches a critical width, the NTM transitions to faster robust growth, which is well described by the nonlinear modified Rutherford equation. Furthermore, this work highlights the role of nonlinear mode coupling in seeding NTMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Toward the origins of binding energy shifts and “satellites” formation during plasma-XPS measurements

In-plasma x-ray photoelectron spectroscopy (plasma-XPS) emerges as a powerful platform for real-time, in situ chemical analysis under conditions relevant to semiconductor processing and other plasma-enabled technologies. This study investigates the origins of binding energy (BE) shifts and the formation of “satellite” peaks observed during plasma-XPS measurements across conductive, dielectric, and gas-phase systems. Using a standard laboratory-based ambient pressure XPS apparatus coupled with an alternating current (AC)-driven capacitively coupled plasma source, we demonstrate that metastable surface species, such as transient Au oxides, can be detected during plasma exposure, revealing chemical states that are hardly accessible using conventional ultrahigh vacuum (UHV) XPS. In dielectric samples (e.g., undoped diamond, sapphire), we observe pressure- and plasma-type-dependent BE shifts of more than 50 eV, attributed to x-ray-induced and plasma-mediated surface charging. These shifts are mitigated at higher pressures/plasmas or in electronegative plasmas (e.g., O 2 ), the latter due to enhanced charge compensation mechanisms involving slow negative ions. For gas-phase species, AC-plasma excitation leads to spectral broadening and the emergence of “satellite” peaks with energy separations of a few electron volts, linked to oscillating local plasma potentials in the probing volume. Furthermore, these findings highlight the complex and important interplay between plasma parameters, surface charging, and local electric fields in shaping XPS spectra. Overall, plasma-XPS emerges as a critical metrological tool for probing transient surface chemistry, with implications for semiconductor processing, material synthesis, and plasma diagnostics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Physics-Infused AI/ML Based Digital-Twin Framework for Flow-Induced-Vibration Damage Prediction in a Nuclear Reactor Heat Exchanger

This report summarizes some of the ongoing work related to the development of an expert-elicitation-digital-twin framework for real time damage state prediction in heat exchanger components of a nuclear reactor. The framework is targeted towards predicting damage associated with coupled low cycle fatigue (associated with regular heat-up, cool-down and power operation transients) and high cycle fatigue (associated with flow induced vibration transients). The overall framework will be based on a NoSQL based database, physics-infused-geometry-dependent virtual-sensor data, different AI/ML techniques-based data-driven-predictive-model applications (Apps) and real-time plant sensor measurements available through few existing sensors. Towards this overall goal, this report updates some of the ongoing work, such as on implementation of a NoSQL Database (such as MongoDB), FE based heat transfer analysis of a heat exchanger (e.g. of a PWR steam generator) for generating geometry-dependent virtual sensor data and evaluation of various AI/ML models such as based on multivariate linear regression, ensembled decision-tree based Random-Forest and Gradient-Boosting regression and high-dimensional-kernel-function-transformation based Support-Vector-Machine regression models. The AI/ML models were evaluated for predicting multi-time-series thermal states at thousands of 3D point-clouds

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Analysis of Molten Salt Thermal Energy Storage Tank Sizing and Preheating

Molten Salt Thermal energy storage (TES) is essential for concentrating solar power (CSP) plants, enhancing their capacity factor and dispatchability to ultimately reduce electricity costs. This technology is becoming more important for other technologies including nuclear power, thermal-electrical storage, and concentrating solar thermal (CST). While molten salt TES tanks have been successfully deployed in lower-temperature parabolic-trough plants, the current state-of-the-art central receiver systems operating at temperatures up to 565 degrees Celsius have experienced costly and frequent failures. These failures are associated with the immaturity of the technology and problems in tank design, fabrication, commissioning and, sometimes, aggressive operation. A primary contributing factor is the absence of a dedicated technical standard for high-temperature molten nitrate TES tanks; designers must currently rely on existing codes like API 650 and ASME Section II, which are insufficient because they fail to account for the unique demands of high temperatures, thermal cycling, and transient conditions.

25 ENERGY STORAGE↗

Ultrafast dynamics of CN radical reactions with chloroform solvent under vibrational strong coupling

Polariton chemistry may provide a new means to control molecular reactivity, permitting remote, reversible modification of reaction energetics, kinetics, and product yields. A considerable body of experimental and theoretical work has already demonstrated that strong coupling between a molecular vibrational mode and the confined electromagnetic field of an optical cavity can alter chemical reactivity without external illumination. However, the mechanisms underlying cavity-altered chemistry remain unclear in large part because the experimental systems examined previously are too complex for detailed analysis of their reaction dynamics. Here, we experimentally investigate photolysis-induced reactions of cyanide radicals with strongly-coupled chloroform (CHCl 3 ) solvent molecules and examine the intracavity rates of photofragment recombination, solvent complexation, and hydrogen abstraction. Here we use a microfluidic optical cavity fitted with dichroic mirrors to facilitate vibrational strong coupling (VSC) of the C–H stretching mode of CHCl 3 while simultaneously permitting optical access at visible wavelengths. Ultrafast transient absorption experiments performed with cavities tuned on- and off-resonance reveal that VSC of the CHCl 3 C–H stretching transition does not significantly modify any measured rate constants, including those associated with the hydrogen abstraction reaction. This work represents, to the best of our knowledge, the first experimental study of an elementary bimolecular reaction under VSC. We discuss how the conspicuous absence of cavity-altered effects in this system may provide insights into the mechanisms of modified ground state reactivity under VSC and help bridge the divide between experimental results and theoretical predictions in vibrational polariton chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Load Following Analysis of the Holos-Quad 10MWe Micro-Reactor with Plant Dynamics Code

This report presents the load following analysis of the “Holos-Quad” 10MWe micro-reactor design developed by HolosGen LLC. This analysis is focused on the Holos-Quad micro-reactor ability to match the changing grid demand at 10%/min rate. The control mechanisms for the plant are identified, simulated, and compared. Based on the comparison results, control strategy for load following of the Holos-Quad micro-reactor is developed. The control strategy and load following capabilities are demonstrated in a full-range down-and-up reactor power transient from 100% to 0% and back to 100% load at 10%/min rate. All calculations are carried out with Argonne’s Plant Dynamics Code.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Data, model inputs, and analysis scripts associated with a manuscript on stream intermittency controls across spatial scales in Pacific Northwest watersheds

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript "Hydroclimatic Memory and Watershed Template Shape Stream Intermittency: Multi-scale Attribution Using Process-based Simulation and Explainable ML" by Niroula et al. (2026), submitted to Water Resources Research (WRR). The study investigates the dominant controls on stream intermittency across local, reach, and watershed scales using a coupled process-based simulation and explainable machine-learning framework. Long-term daily simulations from the Advanced Terrestrial Simulator (ATS) were used to generate wetness states and ponded-depth responses over river-corridor cells. These ATS outputs were then aggregated across scales and used to train XGBoost (eXtreme Gradient Boosting) models. SHAP (SHapley Additive exPlanations) was applied to quantify the relative importance of hydroclimatic forcings, watershed template attributes, and antecedent-memory effects in shaping intermittency behavior. The analysis is carried out for three contrasting Pacific Northwest watersheds: Oak Creek (OCW), American River Watershed (ARW), and H.J. Andrews (HJA). Across these testbeds, the package contains ATS-ready watershed inputs, ATS run configuration and selected output files, model-evaluation data products, intermittency-analysis datasets, machine-learning target-feature tables, SHAP outputs, and notebooks used to organize, analyze, and visualize results. At a high level, the package documents a workflow in which ATS provides the physically based simulation backbone and explainable machine learning is used as a post-processing attribution tool. The contents are intended to support interpretation of the manuscript figures and results, provide context for how intermittency metrics were generated at multiple scales, and preserve the key artifacts needed to understand and reuse the analysis workflow. The package contains a high-level directory summary file (`summary.txt`) and four main content folders (1) `evaluation_plots` contains evaluation figures and supporting evaluation datasets; (2) `intermittency_plots` contains intermittency-focused analysis notebook and prepared datasets; (3) `ml-training-and-shap_values_plots` contains ML training inputs, SHAP outputs, and figure-generation notebooks; and (4) `watershed_mesh_and_ats_input` contains ATS model setup materials, forcing inputs, geometry, and selected run files. More specifically, the `evaluation_plots` folder contains the notebook used for ATS evaluation plotting and site-specific evaluation datasets. These include evapotranspiration and water-balance products for three watersheds, as well as an Oak Creek field-measurement discharge file. The `intermittency_plots` folder contains the notebook used for intermittency analysis and the prepared datasets used to analyze intermittent and non-intermittent wetness behavior across the study watersheds. The `ml-training-and-shap_values_plots` folder contains notebooks and outputs for the machine-learning and explainability workflow. This includes the main XGBoost and SHAP notebook(s), a beeswarm plotting notebook, target-feature tables for machine-learning training, SHAP summary tables, and per-sample SHAP value archives. The `watershed_mesh_and_ats_input` folder contains ATS-related watershed inputs and supporting materials. This includes mesh and shape products, ATS-readable LAI and meteorological forcing inputs, selected ATS spinup and transient-run files, and a watershed workflow example notebook. Subdirectories are organized by watershed where applicable.All files are .cpg (codepage files), .csv (comma-separated values), .dbf (database files), .exo (Exodus mesh format), .h5 (HDF5 format), .ipynb (Jupyter notebooks), .pkl (Python pickle), .prj (projection files), .sh (shell scripts), .shp (shapefile geometry), .shx (shapefile index), .txt (text files), or .xml (markup data).

Advanced Terrestrial Simulator↗

Comparison of pyrometry and thermography for thermal analysis of thermite reactions

This article examines the thermal behavior of a laser ignited thermite composed of aluminum and bismuth trioxide. Temperature data were collected during the reaction using a four-color pyrometer and a high-speed color camera modified for thermography. The two diagnostics were arranged to collect data simultaneously, with similar fields of view and with similar data acquisition rates, so that the two techniques could be directly compared. Results show that at initial and final stages of the reaction, a lower signal-to-noise ratio affects the accuracy of the measured temperatures. Both diagnostics captured the same trends in transient thermal behavior, but the average temperatures measured with thermography were about 750 K higher than those from the pyrometer. This difference was attributed to the lower dynamic range of the thermography camera’s image sensor, which was unable to resolve cooler temperatures in the field of view as well as the photomultiplier tube sensors in the pyrometer. Overall, while the camera could not accurately capture the average temperature of a scene, its ability to capture peak temperatures and spatial data make it the preferred method for tracking thermal behavior in thermite reactions.

36 MATERIALS SCIENCE↗

Viral‐mediated delivery of morphogenic regulators enables leaf transformation in Sorghum bicolor (L.)

Recent advancements in monocot transformation, using leaf tissue as explant material, have expanded the number of grass species capable of transgenesis. However, the complexity of vectors and reliance on inducible excision of essential morphogenic regulators have so far limited widespread application. Plant RNA viruses, such as Foxtail Mosaic Virus (FoMV), present a unique opportunity to express morphogenic regulator genes, such as Babyboom (Bbm), Wuschel2 (Wus2), Wuschel-like homeobox protein 2a (Wox2a) and the GROWTH-REGULATING FACTOR 4 (GRF4) GRF-INTERACTING FACTOR 1 (GIF1) fusion protein transiently in leaf explant tissues. Furthermore, altruistic delivery of conventional and viral vectors could provide opportunities to simplify vectors used for leaf transformation—facilitating vector optimization and reducing reliance on morphogenic regulator gene integration. In this study, both viral and conventional T-DNA vectors were tested for their ability to promote the formation of embryonic calli, a critical step in leaf transformation protocols, using Sorghum bicolor leaf explants. Although conventional leaf transformation vectors yielded viable embryonic calli (43.2 ± 2.9%: GRF4-GIF1, 50.2 ± 3%: Bbm/Wus2), altruistic conventional vectors employing the GRF4-GIF1 morphogenic regulator resulted in improved efficiencies (61.3 ± 4.7%). Altruistic delivery was further enhanced with the use of viral vectors employing both GRF4-GIF1 and Bbm/Wus2 regulators, resulting in 75.1 ± 2.3% and 79.2 ± 2.5% embryonic calli formation, respectively. Embryonic calli generated from both conventional and viral vectors produced shoots expressing fluorescent reporters, which were confirmed using molecular analysis. This work provides an important proof-of-concept for the use of both altruistic vectors and viral-expressed morphogenic regulators for improving plant transformation.

59 BASIC BIOLOGICAL SCIENCES↗

Non-dimensional performance and safety parameters for heat pipes

The use of heat pipes in safety-critical systems such as nuclear microreactors dictates the development of generalized, practical, scalable performance and safety parameters. Traditional dimensional metrics, while informative, lack the universality required for comparative analysis across varying designs and operating regimes. Here, this work introduces a comprehensive set of non-dimensional parameters to characterize heat pipe performance and safety, including capillary performance, effective thermal conductivity, response time, exergetic efficiency, allowable temperature gradients, allowable rate of temperature change, priming coefficients, and factor of safety. A reference heat pipe design representative of microreactor applications was analyzed via the developed parameters using both traditional analytical models and Sockeye simulations under transient and steady-state conditions. Sodium, potassium, and water were evaluated as working fluids to demonstrate the applicability of the framework across a broad temperature range. The proposed non-dimensional parameters effectively captured key thermal-hydraulic behaviors and safety concerns, as was demonstrated via Sockeye simulations. This framework supports the development of design optimization strategies, operational protocols, and safety assurance practices for advanced reactor systems and other high-reliability applications.

42 - ENGINEERING↗

Study of the film boiling heat transfer and two-phase flow interface behavior using image processing

Here, in this study, based on a small-scale quench test facility, the two-phase flow interface behavior during quench transients is visualized and analyzed utilizing an image processing framework. The high-fidelity experimental results obtained for two-phase flow in the current framework can support various studies both in the time domain and in the frequency domain. In particular, visualization of the data obtained from different heating surfaces under different test conditions are used to perform a full-scale transient 2-D vapor film reconstruction. The liquid-vapor interface variations in various heat transfer regimes as well as at the initial film breakup point can be directly obtained through the processed data. Moreover, the temporal variation of the interfacial wave frequency approaching quench is investigated in detail. Based on the high-resolution data obtained for the liquid-vapor interface, the detailed phase velocity and temperature profiles are obtained through theoretic analysis, based on which the film boiling heat transfer coefficient (HTC) can be determined. In addition, an improved film boiling HTC model is developed considering the effects of wall superheat, liquid subcooling temperature, vapor film thickness as well as fluid properties. The model is found to predict film boiling HTC well within 15% error.

42 ENGINEERING↗

Highlighting the Impact of Yaw Control by Parsing Atmospheric Conditions Based on Total Variation: Preprint

Identification of atmospheric conditions within a multivariate atmospheric dataset is a necessary step in the validation of wind plant control strategies. Most often, operating conditions are characterized in terms of aggregated observations and assume that the atmosphere is 'quasi-steady'. Aggregation of observations without regard to covariance between time series discounts the dynamical nature of the atmosphere and is not sufficiently representative of wind plant operating conditions. Identification and characterization of continuous time periods with atmospheric conditions that have a high value for analysis or simulation sets the stage for more advanced model validation and the development of real-time control and operation strategies. Controlling observational data for statistical stationarity highlights significant enhancements to the power production of waked turbines under wake steering wind plant control. Considering the combined energy ratio of the turbine to which yaw offsets are prescribed as well as the waked turbine show only moderate improvements when filtering for total validation. The same quality control that highlights wake steering also serves to underpin the off-nominal operation of the controlled turbine.

17 WIND ENERGY↗

Operating Characteristics of Metal Hydride-Based Solar Energy Storage Systems

Thermochemical energy storage systems, based on a high-temperature metal hydride coupled with a low-temperature metal hydride, represent a valid option to store thermal energy for concentrating solar power plant applications. The operating characteristics are investigated for a tandem hydride bed energy storage system, using a transient lumped parameter model developed to identify the technical performance of the proposed system. The results show that, without operational control, the system undergoes a thermal ratcheting process, causing the metal hydride concentrations to accumulate hydrogen in the high-temperature bed over time, and deplete hydrogen in the low temperature. This unbalanced system is compared with a ’thermally balanced’ system, where the thermal ratcheting is mitigated by thermally balancing the overall system. The analysis indicates that thermally balanced systems stabilize after the first few cycles and remain so for long-term operation, demonstrating their potential for practical thermal energy storage system applications.

25 ENERGY STORAGE↗

ANTARES: a new IRSN computational chain for safety analysis

Since some years, the IRSN has developed a new project to simulate accurately the interaction between neutronic and thermal hydraulic phenomena. This project called ANTARES (Advanced Neutronics and Thermal-hydraulic for the Analysis of Reactor Safety) couples the thermal hydraulic code CATHARE-3 with the 3D neutronic nodal code PARCS. The purpose of this project is to take the advantage of the complementary of these 2 codes, by the knowledge of the 3D core power distribution in CATHARE-3 (computed by PARCS) and improving the thermal hydraulic boundaries conditions in PARCS. In this way, the IRSN diversifies and enlarges its available tools to perform safety analysis with improved accuracy. The current technique usually adopted in France for the safety demonstrations for Reactivity Insertion Accident (RIA) is the so-called 'conservative' approach, which consists of reducing all the feedback (Doppler and moderator effects) and in modifying some physical quantities in such a way to increase a power peak. For this reason, these facilities ('penalties') have been implemented in ANTARES. In this paper we will give two examples of accidental transients that can be simulated with ANTARES: a Main Steam Line Break (MSLB) and an inadvertent boron dilution event. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

High fidelity simulations of contaminant dispersion in an urban environment with comparison to magnetic resonance imaging measurements

The dispersion of a contaminant in an urban environment has the potential to impact a large population of people. In this work, a complex urban canopy flow based on the Oklahoma City downtown business district circa 2003 is studied using Magnetic Resonance Imaging (MRI) and high-fidelity Large Eddy Simulations (LES). MRI is a novel experimental technique that can provide high-resolution measurements in four dimensions (three spatial and temporal) for lab scale models. The experiments and simulations use the same geometry and boundary conditions providing a one-to-one comparison of the two methods. Results are presented on the time-averaged velocity and concentration fields, the temporal dynamics of the concentration plumes for a transient release, and a novel Cloud Identification Algorithm that can separate plumes produced by periodic contaminant releases used for ensemble averaging over many releases. The MRI and LES datasets both include millions of measurement voxels and the comparisons highlight the complex 3D nature of the flow including strong vertical velocities in spanwise street canyons and flow acceleration in streamwise street canyons. The concentration fields are qualitatively similar albeit the LES shows larger dispersion. A quantitative analysis with performance measures compares the datasets pointwise and demonstrates that the two 3D datasets are similar with respect to many measures including a fractional bias of 0.02 (ideal=0.0), correlation coefficient of 0.87 (ideal = 1.0), and the fraction points within a factor of 2 is 0.98 (ideal = 1.0). Plume analysis compares the arrival and residence time of contaminant and is found to vary significantly with location within the urban environment with arrival times between 0 and 1.25 and differences within the contaminant cloud less than 10% at most locations.

54 ENVIRONMENTAL SCIENCES↗