Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic response”

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 343 records · Page 19

Nondestructive Evaluation of Carbon Fiber Reinforced Polymers

The American Society of Mechanical Engineers Boiler and Pressure Vessel Code requires repair and replacement of safety-related piping materials to meet the original Construction Code; however, no construction criteria currently exist for carbon fiber reinforced polymer (CFRP) materials in the nuclear industry. Although nondestructive examination (NDE) techniques for cast and wrought ferritic and austenitic steels are well established, licensees are increasingly deploying novel materials such as CFRP for which assessment of the use of NDE is required, as the inspectability of such materials and the influence of manufacturing processes on inspectability remain insufficiently understood. To address these gaps, Pacific Northwest National Laboratory (PNNL) evaluated the fabrication of several CFRP repair mockups and the effectiveness of NDE methods to support regulatory review of CFRP repairs in nuclear power plant applications. Representative flat-plate mockups containing realistic fabrication defects were manufactured and inspected using manual and automated tap testing, dynamic response spectroscopy (DRS), and ultrasonic testing (UT). The study identified significant fabrication variability, particularly in controlling defect size and epoxy thickness, which strongly influenced defect detectability by NDE methods. Tap testing was effective for shallow defects in thin epoxy configurations but unreliable for thicker repairs and deeper flaws. DRS demonstrated higher sensitivity to dry spots but produced unconfirmed indications for thicker plates, requiring further validation. Conventional UT, particularly at 1.0 MHz, showed the strongest overall capability for detecting a range of defects, although performance degraded with increasing repair thickness and complexity. The results highlight the need to better define critical defect sizes, develop reliable mockup fabrication methods, and validate NDE methods needed to support CFRP repairs.

36 MATERIALS SCIENCE↗

Fed-DeepONet: Stochastic Gradient-Based Federated Training of Deep Operator Networks

The Deep Operator Network (DeepONet) framework is a different class of neural network architecture that one trains to learn nonlinear operators, i.e., mappings between infinite-dimensional spaces. Traditionally, DeepONets are trained using a centralized strategy that requires transferring the training data to a centralized location. Such a strategy, however, limits our ability to secure data privacy or use high-performance distributed/parallel computing platforms. To alleviate such limitations, in this paper, we study the federated training of DeepONets for the first time. That is, we develop a framework, which we refer to as Fed-DeepONet, that allows multiple clients to train DeepONets collaboratively under the coordination of a centralized server. To achieve Fed-DeepONets, we propose an efficient stochastic gradient-based algorithm that enables the distributed optimization of the DeepONet parameters by averaging first-order estimates of the DeepONet loss gradient. Then, to accelerate the training convergence of Fed-DeepONets, we propose a moment-enhanced (i.e., adaptive) stochastic gradient-based strategy. Finally, we verify the performance of Fed-DeepONet by learning, for different configurations of the number of clients and fractions of available clients, (i) the solution operator of a gravity pendulum and (ii) the dynamic response of a parametric library of pendulums.

Moya, Christian↗

Parallel Broadband Femtosecond Reflection Spectroscopy at a Soft X-Ray Free-Electron Laser

X-ray absorption spectroscopy (XAS) and the directly linked X-ray reflectivity near absorption edges yield a wealth of specific information on the electronic structure around the resonantly addressed element. Observing the dynamic response of complex materials to optical excitations in pump–probe experiments requires high sensitivity to small changes in the spectra which in turn necessitates the brilliance of free electron laser (FEL) pulses. However, due to the fluctuating spectral content of pulses generated by self-amplified spontaneous emission (SASE), FEL experiments often struggle to reach the full sensitivity and time-resolution that FELs can in principle enable. Here, we implement a setup which solves two common challenges in this type of spectroscopy using FELs: First, we achieve a high spectral resolution by using a spectrometer downstream of the sample instead of a monochromator upstream of the sample. Thus, the full FEL bandwidth contributes to the measurement at the same time, and the FEL pulse duration is not elongated by a monochromator. Second, the FEL beam is divided into identical copies by a transmission grating beam splitter so that two spectra from separate spots on the sample (or from the sample and known reference) can be recorded in-parallel with the same spectrometer, enabling a spectrally resolved intensity normalization of pulse fluctuations in pump–probe scenarios. We analyze the capabilities of this setup around the oxygen K- and nickel L-edges recorded with third harmonic radiation of the free electron laser in Hamburg (FLASH), demonstrating the capability for pump–probe measurements with sensitivity to reflectivity changes on the per mill level.

36 MATERIALS SCIENCE↗

Sensitivity Analysis of Modal Parameters of a Jacket Offshore Wind Turbine to Operational Conditions

Accurate estimation of offshore wind turbine (OWT) modal parameters has a prominent effect on the design loads, lifetime prediction, and dynamic response of the system. Modal parameters can vary during the operation of OWTs. This paper studies the variation and sensitivity analysis of an OWT’s modal parameters with respect to operational and environmental conditions. Three finite element models of a jacket-supported OWT at the Block Island Wind Farm are created within the OpenSees, SAP2000, and OpenFAST platforms and validated using experimental measurements. The OpenFAST model is used to simulate the modal parameters of the turbine under various wind speed, rotor speed, power, yaw angle, mean sea level, blade pitch angle, and soil spring values. The model-predicted modal parameters of the first fore–aft (FA) and side–side (SS) modes are compared to those identified from experimental measurements. Results from the simulations show that the first FA natural frequency and damping ratio mostly depend on the rotor speed and wind speed, respectively, while yaw angle and mean sea level do not have a visible effect. It is observed that there is about 8% stiffening in the first FA frequency and an aerodynamic damping of 7.5% during the operation of the OWT.

17 WIND ENERGY↗

ESCARGOT: Mapping Vertical Phase Spiral Characteristics Throughout the Real and Simulated Milky Way

Abstract The recent discovery of a spiral pattern in the vertical kinematic structure in the solar neighborhood provides a prime opportunity to study nonequilibrium dynamics in the Milky Way from local stellar kinematics. Furthermore, results from simulations indicate that even in a limited volume, differences in stellar orbital histories allow us to trace variations in the initial perturbation across large regions of the disk. We present ESCARGOT , a novel algorithm for studying these variations in both simulated and observed data sets. ESCARGOT automatically extracts key quantities from the structure of a given phase spiral, including the time since perturbation and the perturbation mode. We test ESCARGOT on simulated data and show that it is capable of accurately recovering information about the time since the perturbation occurred as well as subtle differences in phase spiral morphology due to stellar locations in the disk at the time of perturbation. We apply ESCARGOT to kinematic data from data release 3 of the Gaia mission in bins of guiding radius. We show that similar structural differences in morphology occur in the Gaia phase spirals as a function of stellar orbital history. These results indicate that the phase spirals are the product of a complex dynamical response in the disk with large-scale coupling between different regions of phase space.

Darragh-Ford, Elise (ORCID:0000000288005652)↗

Impact of volcanic eruptions on CMIP6 decadal predictions: a multi-model analysis

Abstract. In recent decades, three major volcanic eruptions of different intensity have occurred (Mount Agung in 1963, El Chichón in 1982 and Mount Pinatubo in 1991), with reported climate impacts on seasonal to decadal timescales that could have been potentially predicted with accurate and timely estimates of the associated stratospheric aerosol loads. The Decadal Climate Prediction Project component C (DCPP-C) includes a protocol to investigate the impact of volcanic aerosols on the climate experienced during the years that followed those eruptions through the use of decadal predictions. The interest of conducting this exercise with climate predictions is that, thanks to the initialisation, they start from the observed climate conditions at the time of the eruptions, which helps to disentangle the climatic changes due to the initial conditions and internal variability from the volcanic forcing. The protocol consists of repeating the retrospective predictions that are initialised just before the last three major volcanic eruptions but without the inclusion of their volcanic forcing, which are then compared with the baseline predictions to disentangle the simulated volcanic effects upon climate. We present the results from six Coupled Model Intercomparison Project Phase 6 (CMIP6) decadal prediction systems. These systems show strong agreement in predicting the well-known post-volcanic radiative effects following the three eruptions, which induce a long-lasting cooling in the ocean. Furthermore, the multi-model multi-eruption composite is consistent with previous work reporting an acceleration of the Northern Hemisphere polar vortex and the development of El Niño conditions the first year after the eruption, followed by a strengthening of the Atlantic Meridional Overturning Circulation the subsequent years. Our analysis reveals that all these dynamical responses are both model- and eruption-dependent. A novel aspect of this study is that we also assess whether the volcanic forcing improves the realism of the predictions. Comparing the predicted surface temperature anomalies in the two sets of hindcasts (with and without volcanic forcing) with observations we show that, overall, including the volcanic forcing results in better predictions. The volcanic forcing is found to be particularly relevant for reproducing the observed sea surface temperature (SST) variability in the North Atlantic Ocean following the 1991 eruption of Pinatubo.

Bilbao, Roberto (ORCID:0000000307294980)↗

Wind turbine drivetrains: state-of-the-art technologies and future development trends

Abstract. This paper presents the state-of-the-art technologies and development trends of wind turbine drivetrains – the system that converts kinetic energy of the wind to electrical energy – in different stages of their life cycle: design, manufacturing, installation, operation, lifetime extension, decommissioning and recycling. Offshore development and digitalization are also a focal point in this study. Drivetrain in this context includes the whole power conversion system: main bearing, shafts, gearbox, generator and power converter. The main aim of this article is to review the drivetrain technology development as well as to identify future challenges and research gaps. The main challenges in drivetrain research identified in this paper include drivetrain dynamic responses in large or floating turbines, aerodynamic and farm control effects, use of rare-earth material in generators, improving reliability through prognostics, and use of advances in digitalization. These challenges illustrate the multidisciplinary aspect of wind turbine drivetrains, which emphasizes the need for more interdisciplinary research and collaboration.

17 WIND ENERGY↗

Roper Resonance Structure and Exploration of Emergent Hadron Mass from CLAS Electroproduction Data

The N(1440)1/2 + nucleon resonance, first identified in 1964 by Roper and collaborators in analyses of πN hadroproduction data, has continued to provide pivotal insights that serve to advance our understanding of nucleon excited states. In this contribution, we present results from studies of the structure of the Roper resonance based on exclusive πN and π + π − p electro production data measured with the CLAS detector at Jefferson Lab. These analyses have revealed the Roper resonance as a complex interplay between an inner core of three dressed quarks and an external meson–baryon cloud. Analyses of the CLAS results on the evolution of the Roper resonance electroexcitation amplitudes with photon virtuality Q 2 , within the framework of the Continuum Schwinger Method, have conclusively demonstrated the capability to gain insight into the strong interaction dynamics responsible for generating more than 98% of hadron mass. Further extension of such studies to higher Q 2 , through experiments currently underway with the CLAS12 detector and in the future with a potential CEBAF energy upgrade to 22 GeV, offers the only foreseeable opportunity to explore the full range of distances where the dominant portion of hadron mass and resonance structure emerge.

Mokeev, Victor I. [Thomas Jefferson National Accel↗

Numerical Model Development of a Variable-Geometry Attenuator Wave Energy Converter: Preprint

Because the wave energy industry is still in its infancy, an optimal design for wave energy converters (WECs) has yet to be established; more work is needed to explore various cost reduction pathways. The primary cost-reduction pathway considered for this work is the optimization of the geometric profile on an attenuator WEC to maximize power production while, at the same time, minimizing capital expenditures through the use of variable-geometry modules. In this investigation, the variable geometry modules consist of inflatable bags placed on either side of a base central steel cylinder that would be inflated in low moderate sea states to maximize power capture and then deflated in moderate-extreme sea states to minimize wave loading. The numerical model and simulation of the attenuator WEC were developed and completed using WEC-Sim, which is an open-source code that is appropriate for use in evaluating the dynamic response of the different WEC models in operational seas. The power production estimates were obtained from the Wave Energy Prize (WEP) sea states, which are representative of U. S. deployment sites, to calculate the average climate capture width that is used in the WEP ACE calculation. Preliminary capital expenditure costs were obtained assuming the base central steel cylinder mass was equal to the fluid displaced mass, minus the mass of the variable-geometry bags. The additional weight required to offset the additional buoyancy from the variable-geometry bags was assumed to come from the addition of seawater ballast. The variable-geometry attenuator model was found to have a similar power capture efficiency as a fixed-body model, but is expected to have a lower characteristic capital expenditure given its more streamlined profile, which demonstrates that variable-geometry modules may provide a realistic cost-reduction pathway to help design a more cost-competitive WEC.

50 EE - Wind and Water Power Program - Water (EE-4↗

Probing and manipulating magnetic and 2D quantum materials using ultrafast laser and high harmonic sources

The ability to probe the full dynamic response of quantum materials on the length- and time-scales (Å to attoseconds on up) fundamental to charge, spin and phonon interactions is leading to a host of new discoveries. The coupled interactions between charges, spins, orbital and lattice degrees of freedom are key to determining the state of a material, whether metallic, insulating, superconducting or magnetic. Under thermal equilibrium conditions these states can be tuned by varying the temperature, pressure, chemical doping or dimensionality. However, their inherent complexity and rich phase landscape make them challenging to understand or manipulate deterministically. Fortunately, ultrafast light sources have undergone remarkable advances in recent years, achieving what was merely a dream three decades ago, i.e., full coherent control of light fields spanning the THz to the X-ray regions. Here, these new capabilities are providing powerful new tools for coherently manipulating and probing quantum materials using light. Recent advances in high harmonic generation (HHG) have provided new techniques for uncovering new phases and couplings in 2D and magnetic materials, using femtosecond lasers to coherently manipulate and tune the properties and interactions in the material, and HHG as a probe to map out these new phases and couplings. There are several beautiful recent observations of exciting new states in materials using light, but in some cases, the limited experimental data points made it challenging to validate an interpretation. The new ability of ultrafast HHG-based photoemission and magneto-optic spectroscopies to precisely probe multiple phases, electronic and magnetic order parameters, when combined with scanning the laser excitation fluence, is allowing us to map how material phases, magnetic states, and electron-phonon-spin couplings in these delicate materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Reliable Power Rating of Perovskite PV Modules: Preprint

As the perovskite technology is ramping up into commercialization, reliable and accurate power rating of large-size perovskite modules becomes a prominent aspect for its future deployment in the PV market. It is known that the performance calibration of perovskite PV devices is very challenging due to its complex dynamic response during a conventional current-voltage (IV) measurement. PV researchers have previously proposed several steady-state performance calibration methods to reliably extract PV efficiencies, but mostly focus on small area research-type cells. In this paper, we emphasize the importance of reliable performance calibration on large-size perovskite modules. Extending the NREL Cell and Module Performance (CMP) group’s steady-state performance calibration protocol (i.e., Asymptotic PMAX Scan) for perovskite cells to modules, we justify the necessity of reporting steady-state efficiencies for perovskite cells and discuss the challenges of applying this protocol to modules. We also present our protocol for Maximum Power Point Tracking (MPPT), which is a technique often used for performance calibration of perovskite cells and modules, and show a comparison between MPPT and Asymptotic PMAX. Using MPPT we demonstrate the interplay between metastability and degradation in perovskite modules, and emphasize the necessity to develop preconditioning protocols for stabilizing these devices. Our aim is to promote development of consensus protocols for performance calibration of perovskite modules, and to advance their credible power ratings, which will be beneficial to the growth of perovskite technology in the PV market.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Networked HIL Simulation System for Modeling Large-scale Power Systems

This paper presents a network hardware-in-the-loop (HIL) simulation system for modeling large-scale power systems. Researchers have developed many HIL test systems for power systems in recent years. Those test systems can model both microsecond-level dynamic responses of power electronic systems and millisecond-level transients of transmission and distribution grids. By integrating individual HIL test systems into a network of HIL test systems, we can create large-scale power grid digital twins with flexible structures at required modeling resolution that fits for a wide range of system operating conditions. This will not only significantly reduce the need for field tests when developing new technologies but also greatly shorten the model development cycle. In this paper, we present a networked OPAL-RT based HIL test system for developing transmission-distribution coordinative Volt-VAR regulation technologies as an example to illustrate system setups, communication requirements among different HIL test systems, and connection mechanisms. Impacts of communication delays, information exchange cycles, and computing delays are illustrated. Simulation results show that the performance of a networked HIL test system is satisfactory.

digital twin, distribution system, co-simulation, ↗

In Situ Wind and Turbulence Measurements in a Field of Full-Size Parabolic Trough Solar Collectors

Concentrated Solar Power (CSP) is a promising method for using Solar power for electricity generation with thermal energy storage. One of the primary drivers of structural design costs of CSP collector structures is wind loading. To date, the design of these structures has relied on data from wind tunnels that do not adequately capture the dynamic effects observed at scale. NREL initiated a field measurement campaign at the operational Nevada Solar One (NSO) powerplant that uses parabolic troughs as solar collectors. The aim of the project is a detailed characterization of prevailing wind and turbulence conditions and resulting operational loads on parabolic troughs, providing insights on structural dynamic response, and generating a first-of-a-kind, comprehensive, high-resolution wind-loading dataset available for validating simulations of wind loading on collector structures. The measurements at NSO consist of Sonic anemometers on masts at different heights to characterize the incoming flow and conditions within and above the trough field, in combination with a Doppler Lidar scanning the horizontal plane above the troughs. The wind measurements at NSO have been continuously operating since October 2021 and provide a year-long dataset characterizing wind and turbulence conditions. The structural load measurements start in November 2022 and will complement the wind measurements. In this poster, we present first results of the wind measurement campaign by highlighting days with different atmospheric flow conditions. We identified three main factors altering the flow over the parabolic troughs: Wind speed, wind direction, and the angle of the parabolic troughs. The highest loads are expected when the wind blows perpendicular to the trough rows. In this case, the first rows experience the highest loads and block the subsequent rows, creating conditions with decreased wind speed and enhanced turbulence within the trough field. Also, turbulent length scales change. These affected conditions produce unique load cases on the structures that will be captured by the load measurements. If the wind blows along the trough rows, wind and turbulence conditions are less influenced by the troughs.

atmospheric turbulence↗

In Situ Wind and Turbulence Measurements in a Field of Full-Size Parabolic Trough Collectors

Concentrated Solar Power (CSP) is a promising method for using Solar power for electricity generation with thermal energy storage. One of the primary drivers of structural design costs of CSP collector structures is wind loading. To date, the design of these structures has relied on data from wind tunnels that do not adequately capture the dynamic effects observed at scale. NREL initiated a field measurement campaign at the operational Nevada Solar One (NSO) powerplant that uses parabolic troughs as solar collectors. The aim of the project is a detailed characterization of prevailing wind and turbulence conditions and resulting operational loads on parabolic troughs, providing insights on structural dynamic response, and generating a first-of-a-kind, comprehensive, high-resolution wind-loading dataset available for validating simulations of wind loading on collector structures. The measurements at NSO consist of Sonic anemometers on masts at different heights to characterize the incoming flow and conditions within and above the trough field, in combination with a Doppler Lidar scanning the horizontal plane above the troughs. The wind measurements at NSO have been continuously operating since October 2021 and provide a year-long dataset characterizing wind and turbulence conditions. The structural load measurements start in November 2022 and will complement the wind measurements. In this poster, we present first results of the wind measurement campaign by highlighting days with different atmospheric flow conditions. We identified three main factors altering the flow over the parabolic troughs: Wind speed, wind direction, and the angle of the parabolic troughs. The highest loads are expected when the wind blows perpendicular to the trough rows. In this case, the first rows experience the highest loads and block the subsequent rows, creating conditions with decreased wind speed and enhanced turbulence within the trough field. Also, turbulent length scales change. These affected conditions produce unique load cases on the structures that will be captured by the load measurements. If the wind blows along the trough rows, wind and turbulence conditions are less influenced by the troughs.

concentrated solar power↗

Unstructured mesh-based multi-physics reactor transient simulation with iMC code

In this study, we present a time-dependent Monte Carlo simulation result about a complex fuel-loaded light water reactor system with thermal-hydraulic feedback using the KAIST iMC code. The intra-pin temperature distribution is calculated based on Monte Carlo-tallied detailed power distribution and the finite element heat transfer method. The temperature effect of fuel, absorber, and coolant are adequately reflected into the dynamic response by the on-the-fly Doppler broadened cross-section correction. For the time-dependent Monte Carlo transport simulation, the predictor-corrector quasi-static Monte Carlo method is used for this study. We tested the coupled analysis framework with a 17-by-17 fuel assembly loaded with the centrally-shielded burnable absorber (CSBA) fuel element. The temperature feedback effect on the material cross-section and the moderator density incurred negative reactivity in response to the external positive reactivity insertion, showing the inherent stability of the reactor system. The pin-wise power and temperature change during the system transient are presented and discussed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Performance Evaluation of Multi-Vendor Grid-Forming Inverters for Grid-Connected Operation Through Hardware Experimentation

This paper presents the functional performance evaluation tests of multiple (three) commercial grid-forming (GFM) inverters when they operate in parallel with the grid through hardware experiments. The goal of these tests is to explore and benchmark the GFM inverters' functionalities and dynamic response when they are operated in parallel with power grids. Both steady-state (changing the inverter's frequency and voltage droop) and transient (adding step changes to grid's frequency/voltage) tests are performed for each GFM inverter with the same testing circuit and testing protocol. The key findings are summarized as follows: 1) The GFM inverters can be dispatched through frequency and voltage droop intercepts to output the target power; and 1) the GFM inverters automatically respond to system frequency and voltage events to output the needed power; 3) all the GFM inviters shows stability issues when absorbing reactive power from the grid.

frequency droop↗

Evaporator Temperature Transient Testing of a High-Performance Sodium Filled Heat Pipe

Heat pipes are two-phase heat transfer devices that enable passive removal of heat from the reactor core to the power conversion system in heat pipe-cooled microreactor designs. Experimental investigations of heat pipe transients are needed for technology demonstration, verification and validation of numerical codes, and the establishment of regulatory requirements. The Single Primary Heat Extraction and Removal Emulator (SPHERE) facility at Idaho National Laboratory (INL) serves as a platform for evaluating the dynamic response of high-temperature heat pipes under a variety of operating conditions. The present work details the experimental investigation of a high-performance, defined as over 2 kW sodium heat pipe subjected to rapid input power fluctuations induced by sudden changes in the evaporator temperature setpoint. In addition, the heat pipe was subjected to an asymmetrical heat load where a subset of heaters operated at 30% and 70% below their nominal power. These experimental conditions were chosen to simulate thermal and operational stresses expected to be encountered in microreactors to provide data on heat pipe behavior during such important transient events. Key data and performance metrics, including time series of temperatures and strains, axial temperature profiles, thermal response times, and heat transfer capabilities, the thermal output over thermal input, were reported and discussed. The results highlight the resilience of heat pipes, revealing their potential to maintain thermal stability and efficiency under varying power loads. Lastly, the paper concludes with a discussion on the significance of the results and their implications for future research.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Rapid eco‐phenotypic feedback and the temperature response of biomass dynamics

Biomass dynamics capture information on population dynamics and ecosystem-level processes (e.g., changes in production over time). Understanding how rising temperatures associated with global climate change influence biomass dynamics is thus a pressing issue in ecology. The total biomass of a species depends on its density and its average mass. Consequently, disentangling how biomass dynamics responds to increasingly warm and variable temperatures ultimately depends on understanding how temperature influences both density and mass dynamics. Here, we address this issue by keeping track of experimental microbial populations growing to carrying capacity for 15 days at two different temperatures, and in the presence and absence of temperature variability. We develop a simple mathematical expression to partition the contribution of changes in density and mass to changes in biomass and assess how temperature responses in either one influence biomass shifts. Moreover, we use time-series analysis (Convergent Cross Mapping) to address how temperature and temperature variability influence reciprocal effects of density on mass and vice versa. We show that temperature influences biomass through its effects on density and mass dynamics, which have opposite effects on biomass and can offset each other. We also show that temperature variability influences biomass, but that effect is independent of any effects on density or mass dynamics. Last, we show that reciprocal effects of density and mass shift significantly across temperature regimes, suggesting that rapid and environment-dependent eco-phenotypic dynamics underlie biomass responses. Overall, our results connect temperature effects on population and phenotypic dynamics to explain how biomass responds to temperature regimes, thus shedding light on processes at play in cosmopolitan and abundant microbes as the world experiences increasingly warm and variable temperatures.

54 ENVIRONMENTAL SCIENCES↗