SEARCH · Engineering Papers

Results for “evaporation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54

All sources NASA NTRS DOE OSTI

NASA NTRS2014

A Thermal Analysis of a Hot-Wire Probe for Icing Applications

This paper presents a steady-state thermal model of a hot-wire instrument applicable to atmospheric measurement of water content in clouds. In this application, the power required to maintain the wire at a given temperature is used to deduce the water content of the cloud. The model considers electrical resistive heating, axial conduction, convection to the flow, radiation to the surroundings, as well as energy loss due to the heating, melting, and evaporation of impinging liquid and or ice. All of these parameters can be varied axially along the wire. The model further introduces a parameter called the evaporation potential which locally gauges the maximum fraction of incoming water that evaporates. The primary outputs of the model are the steady-state power required to maintain a spatially-average constant temperature as well as the variation of that temperature and other parameters along the wire. The model is used to understand the sensitivity of the hot-wire performance to various flow and boundary conditions including a detailed comparison of dry air and wet (i.e. cloud-on) conditions. The steady-state power values are compared to experimental results from a Science Engineering Associates (SEA) Multi-Element probe, a commonly used water-content measurement instrument. The model results show good agreement with experiment for both dry and cloud-on conditions with liquid water content. For ice, the experimental measurements under read the actual water content due to incomplete evaporation and splashing. Model results, which account for incomplete evaporation, are still higher than experimental results where the discrepancy is attributed to splashing mass-loss which is not accounted in the model.

Thermal Analysis↗

NASA NTRS2017

Molecular Composition Analysis of Distant Targets

This document is the Final Report for NASA Innovative Advanced Concepts (NIAC) Phase I Grant 15-NIAC16A-0145, titled Molecular Composition Analysis of Distant Targets. The research was focused on developing a system concept for probing the molecular composition of cold solar system targets, such as Asteroids, Comets, Planets and Moons from a distant vantage, for example from a spacecraft that is orbiting the target (Hughes et al., 2015). The orbiting spacecraft is equipped with a high-power laser, which is run by electricity from photovoltaic panels. The laser is directed at a spot on the target. Materials on the surface of the target are heated by the laser beam, and begin to melt and then evaporate, forming a plume of asteroid molecules in front of the heated spot. The heated spot glows, producing blackbody illumination that is visible from the spacecraft, via a path through the evaporated plume. As the blackbody radiation from the heated spot passes through the plume of evaporated material, molecules in the plume absorb radiation in a manner that is specific to the rotational and vibrational characteristics of the specific molecules. A spectrometer aboard the spacecraft is used to observe absorption lines in the blackbody signal. The pattern of absorption can be used to estimate the molecular composition of materials in the plume, which originated on the target. Focusing on a single spot produces a borehole, and shallow subsurface profiling of the targets bulk composition is possible. At the beginning of the Phase I research, the estimated Technology Readiness Level (TRL) of the system was TRL-1. During the Phase I research, an end-to-end theoretical model of the sensor system was developed from first principles. The model includes laser energy and optical propagation, target heating, melting and evaporation of target material, plume density, thermal radiation from the heated spot, molecular cross section of likely asteroid materials, and estimation of the absorption profile at a distant spectrometer. Results obtained by executing simulations based on the model provide compelling evidence that the concept of remote laser evaporative molecular absorption spectroscopy is feasible. In this document, technical details of the model are presented, and results of simulations are described that indicate the utility of the proposed sensor system. Additionally, an asteroid rendezvous mission is analyzed, with a survey of system requirements to accomplish molecular composition analysis of the asteroid. Based on positive theoretical results obtained during Phase I, the estimated TRL of the system is now TRL-2. This document also describes potential future research and experimentation that could push the system to TRL-4 within 2 years. Steps required for construction of a laboratory prototype are described. An experiment to test predictions of the theory is described, based on the laboratory prototype setup.

asteroid↗

NASA NTRS2023

Land Surfaces at the Tipping-Point for Water and Energy Balance Coupling

The surface water and energy balances can be coupled or uncoupled depending on whether the evaporation regime is water-limited or energy-limited. As the landscape loses soil moisture during drydowns, a transition between the regimes may occur, which signifies a nonlinear change in water-energy-carbon coupling. Regions that switch often between these two regimes, that is, are dominated by neither regime, are particularly vulnerable to climate variability and change. To robustly identify these tipping points, we identify drydown events based on global soil moisture data sets from remote sensing. The event identification does not rely on precipitation information and is robust with respect to measurement noise. Then, the soil moisture thresholds delineating the evaporation regime transitions are determined by Sequential Monte Carlo Sampling and a two-stage parametrization strategy. Based on the estimated soil moisture thresholds across the globe, we estimate observation-based water availability indices which quantify the nonlinear controls of soil moisture on evaporation. This framework is tested and applied globally using Soil Moisture Active Passive soil moisture retrievals. Combined with a new tippling-point metric that describes the frequency of evaporation regime transitions, we identify regions that switch often between different evaporation regimes at the global scale. Given unit shifts in soil moisture, these regions will experience the most change in how their surface water and energy are coupled.

Jianzhi Dong↗

DOE OSTI2020

Hot Gravitons and Gravitational Waves From Kerr Black Holes in the Early Universe

Any abundance of black holes that was present in the early universe will evolve as matter, making up an increasingly large fraction of the total energy density as space expands. This motivates us to consider scenarios in which the early universe included an era that was dominated by low-mass ($M < 5\times 10^8$ g) black holes which evaporate prior to primordial nucleosynthesis. In significant regions of parameter space, these black holes will become gravitationally bound within binary systems, and undergo mergers before evaporating. Such mergers result in three potentially observable signatures. First, any black holes that have undergone one or more mergers will possess substantial angular momentum, causing their Hawking evaporation to produce significant quantities of high-energy gravitons. These products of Hawking evaporation are predicted to constitute a background of hot ($\sim$eV-keV) gravitons today, with an energy density corresponding to $\Delta N_{\rm eff} \sim 0.01-0.03$. Second, these mergers will produce a stochastic background of high-frequency gravitational waves. And third, the energy density of these gravitational waves can be as large as $\Delta N_{\rm eff} \sim 0.3$, depending on the length of time between the mergers and evaporation. These signals are each potentially within the reach of future measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

A conservative, interface-resolved, compressible framework for the modeling and simulation of liquid/gas phase change

This paper presents a method for simulating evaporation in a compressible, interface-resolved framework appropriate for modeling problems of engineering interest. In order to achieve robustness and broad applicability, the method has been designed to discretely enforce consistent mass and thermal energy transport at the phase interface, to globally conserve mass, momentum, and energy, and to be capable of modeling compressible and incompressible systems. Verification is performed via the Sod-shock test, one-dimensional heat conduction, evaporation from a planar interface, and evaporation of three-dimensional droplets. Convergence with increasing mesh resolution is demonstrated in all tested configurations, and conservation is maintained near machine precision for a translating droplet. Conservation and accurate phase change rates are preserved at the low numerical resolutions commonly encountered in engineering calculations. Following verification, the method is validated by comparison to an empirical correlation for evaporating droplets in high temperature crossflow, and the presentation concludes with the simulation of an iso-octane spray at conditions representative of gasoline direct injection. In conclusion, successful verification, validation, and demonstrated practical utility suggest the method to be an accurate, efficient, and robust approach for the study of phase change in engineering systems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Estimating Future Surface Water Availability Through an Integrated Climate‐Hydrology‐Management Modeling Framework at a Basin Scale Under CMIP6 Scenarios

Abstract Climate change and increasing water demand due to population growth pose serious threats to surface water availability. The biggest challenge in addressing these threats is the gap between climate science and water management practices. Local water planning often lacks the integration of climate change information, especially with regard to its impacts on surface water storage and evaporation as well as the associated uncertainties. Using Texas as an example, state and regional water planning relies on the use of reservoir “Firm Yield” (FY)—an important metric that quantifies surface water availability. However, this existing planning methodology does not account for the impacts of climate change on future inflows and on reservoir evaporation. To bridge this knowledge gap, an integrated climate‐hydrology‐management (CHM) modeling framework was developed, which is generally applicable to river basins with geographical, hydrological, and water right settings similar to those in Texas. The framework leverages the advantages of two modeling approaches—the Distributed Hydrology Soil Vegetation Model (DHSVM) and Water Availability Modeling (WAM). Additionally, the Double Bias Correction Constructed Analogues method is utilized to downscale and incorporate Coupled Model Intercomparison Project Phase 6 GCMs. Finally, the DHSVM simulated naturalized streamflow and reservoir evaporation rate are input to WAM to simulate reservoir FY. A new term—“Ratio of Firm Yield” (RFY)—is created to compare how much FY changes under different climate scenarios. The results indicate that climate change has a significant impact on surface water availability by increasing reservoir evaporation, altering the seasonal pattern of naturalized streamflow, and reducing FY.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Two-Phase Cooling for High Packaging Density Planar Inverter: Cooperative Research and Development Final Report, CRADA Number CRD-15-00592

The project objective was to utilize two-phase heat transfer to enable compact and power-dense inverters. The high heat transfer capacity and near isothermal characteristics of two-phase cooling can enable higher power densities, reducing the size and cost of the power electronics system. This project focused on developing efficient, passive (no pump or compressor) two-phase cooling strategies that consist of an air-cooled condenser and evaporator. This report describes the modeling work conducted to design the evaporator and air-cooled condenser for a prototype inverter development work underway in the John Deere. Prototypes of the final evaporator and condenser designs were then fabricated via 3D printing for testing. Experiments on the evaporator confirmed that this component can meet the heat load and junction temperatures requirements for Si IGBT (insulated-gate bipolar transistor) inverter. Experiments with the air-cooled condenser were not possible due to leaks through the condenser walls.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2021

Analysis of Tank 38H (HTF-38-21-24, -25) and Tank 43H (HTF-43-21-26, -27) Samples for Support of the Enrichment Control and Corrosion Control Programs

Feed limits have been established for the 2H-Evaporator system to ensure nuclear criticality is not possible and corrosion is minimized. These limits are protected by the Enrichment Control Program (ECP) and the Corrosion Control Program (CCP) that require periodic sampling and analysis to confirm that the waste supernate composition stays within the limits. Savannah River Remediation (SRR) obtained samples from two different heights within each of the two waste tanks supporting the 2H-Evaporator operations on March 26, 2021. The Tank 38H (evaporator drop tank) and Tank 43H (evaporator feed tank) samples were received by the Savannah River National Laboratory (SRNL) Shielded Cells on March 26, 2021. Analysis of these samples provides information necessary for determining compliance with the ECP and CCP. The sample characterization was requested via a Technical Task Request (TTR) and conducted based on a Task Technical and Quality Assurance Plan (TTQAP).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Journal paper

Laser technology is being widely studied for controlled energy deposition for a range of applications, including flow control, ignition, combustion, and diagnostics. The absorption and scattering of laser radiation by liquid droplets in aerosols affects propagation of the laser beam in the atmosphere, while the ignition and combustion characteristics in combustion chambers are influenced by the evaporation rate of the sprayed fuel. In this work, we present a mathematical model built on OpenFOAM for laser heating and evaporation of a single droplet in the diffusion-dominated regime taking into account absorption of the laser radiation, evaporation process, and vapor flow dynamics. The developed solver is validated against available experimental and numerical data for heating and evaporation of ethanol and water droplets. The two main regimes—continuous and pulsed laser heating—are explored. For continuous laser heating, the peak temperature is higher for larger droplets. For pulsed laser heating, when the peak irradiance is close to transition to the boiling regime, the temporal dynamics of the droplet temperature does not depend on the droplet size. With the empirical normalization of time, the dynamics of the droplet shrinkage and cooling are found to be independent of droplet sizes and peak laser intensities.

42 ENGINEERING↗

DOE OSTI2024

High performance heat exchangers for next generation hybrid cooler (CRADA NFE-19-07891 Final Report)

The evaporative cooling process has been successfully deployed in multiple energy conversion processes such as power generation, process cooling, heating, ventilation & cooling (HVAC), and commercial and industrial refrigeration. The heat exchanger is a very common yet extremely critical component of such systems. The device is important since the performance of the overall cooling system heavily relies on the thermal-hydraulic efficiency of the heat exchanger. The possible wet operation as evaporative cooling heat exchanger complicates the issue. Indeed, based on the hybrid nature of operation (fully dry, partially wet, fully wet) the efficiency can vary significantly and can cause performance variation for the whole system. Understanding the importance of the issue, Oak Ridge National Laboratory (Contractor) and Baltimore Aircoil Company (Participant) propose an activity focused on the analysis of the existing heat exchanger technology used in air and evaporative cooling processes and design, development and demonstration of a novel heat exchanger design based on porous materials. The potential candidates include metallic wire meshes and metal foams carefully designed for the proposed application. Prior research has indicated that these materials have great potential for deployment in thermal systems, particularly heat sinks and heat exchangers. It is expected that the successful completion of the project will lead to novel heat exchanger technology for hybrid systems (capable of dry and wet operation) by exploring newly emerging materials and their potential for heat transfer application. The overall goal is to develop a next generation heat exchanger technology which can be deployed in evaporative cooling systems. The impact of the such development is significant since, it will not only make the overall system highly efficient but will also, reduce the total refrigerant charge in the heat exchanger, a critical aspect for the deployment of refrigerants with lower Global Warming Potential (GWP).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1994

Dynamics and Thermodynamics of Dense and Dilute Clusters of Drops

Collective drop behavior is encountered in sprays produced for a variety of applications. The existence of clusters of drops indicates that the interaction among the drops is important in determining the dynamics of the drops because the drop proximity changes the flow around the individual drops in a way that affects the drag on each drop. Additionally, if there is a phase change between the liquid drops and the gaseous surroundings (either evaporation or condensation), this will also influence the flow around the individual drops; and phase change is also affected by the drop proximity. What will be discussed are studies of collections of drops where the dynamics and thermodynamics of the drops are coupled. This means the modification of the drag force due to the drop number density, the evaporation rate and the slip velocity are all taken into account. In turn, the enhancement in the evaporation rate due to the slip velocity , and the hindering of evaporation due to the drop proximity are also part of the models.

thermodynamics↗

DOE OSTI2024

Molecular beam scattering from flat jets of liquid dodecane and water

Abstract Molecular beam experiments in which gas molecules are scattered from liquids provide detailed, microscopic perspectives on the gas–liquid interface. Extending these methods to volatile liquids while maintaining the ability to measure product energy and angular distributions presents a significant challenge. The incorporation of flat liquid jets into molecular beam scattering experiments in our laboratory has allowed us to demonstrate their utility in uncovering dynamics in this complex chemical environment. Here, we summarize recent work on the evaporation and scattering of Ne, CD 4 , ND 3 , and D 2 O from a dodecane flat liquid jet and present first results on the evaporation and scattering of Ar from a cold salty water jet. In the evaporation experiments, Maxwell–Boltzmann flux distributions with a cos θ angular distribution are observed. Scattering experiments reveal both impulsive scattering (IS) and trapping followed by thermal desorption (TD). Super‐specular scattering is observed for all four species scattered from dodecane and is attributed to anisotropic momentum transfer to the liquid surface. In the IS channel, rotational excitation of the polyatomic scatterers is a significant energy sink, and these species accommodate more readily on the dodecane surface compared to Ne. Our preliminary results on cold salty water jets suggest that Ar atoms undergo some vapor‐phase collisions when evaporating from the liquid surface. Initial scattering experiments characterize the mechanisms of Ar interacting with an aqueous jet, allowing for comparison to dodecane systems. Key Points Molecular beam scattering from flat liquid jets is a powerful technique to elucidate mechanistic detail at the gas–liquid interface. Previous dodecane scattering experiments have uncovered angularly‐resolved thermal desorption fractions and energy transfer at the interface for several small molecule scatterers. Preliminary results on scattering from cold salty water reveal mechanisms of interaction between argon and an aqueous jet.

Yang, Walt↗

DOE OSTI2020

Numerical and experimental analysis of precooled desiccant system

This paper presents a numerical and experimental study of four novel, desiccant air-conditioning systems based on pre-cooling of the ambient air. The purpose of the study is to: 1) establish which configuration is most profitable for the future studies through numerical and experimental methods; 2) verify the concept of precooled system and confirm that it is working in experimental conditions (for this purpose two prototypes of evaporative cooling units were developed and compared for different operation arrangements with desiccant wheel); 3) Analyze practical aspects of the systems, which includes water consumption analysis, economic analysis and establishing operational modes during a cooling season in temperate climate conditions (the proposed solution has a significant advantage during operation in a temperate climate, since it can utilize its first evaporative cooler the only source of cooling energy when the humidity of outdoor air allows to to meet the latent cooling load). Four arrangements of the system with different dew point evaporative air cooler prototypes were selected and developed for comparative analysis. The system was simulated as a whole, also, each component was individually tested on measurement station. develop the optimal arrangement of the system through numerical and experimental analysis. Four different performance indicators were considered: thermal COP (Coefficient Of Performance), the ERR (Energy Efficiency Ratio), coefficient a describing the potential in compensating indoor cooling loads, water consumption ratio. Also costs of the systems were compared through economic analysis. It was established that arrangement which is characterized by best-averaged performance in terms of all considered factors is a unit combined with two different types of evaporative coolers (marked as System B)- this is caused by the favorable arrangement of two units, where one is better suited to the pre-cooler role, while the other is better suited as the post-cooler.

42 ENGINEERING↗

DOE OSTI2024

Characterizing and improving the performance of molten-salt-steam heat exchangers in concentrating solar power plants

Shell-and-tube heat exchangers (HXs) for steam generation from molten salts in concentrating solar power (CSP) plants experience thermal fatigue due to significant temperature gradients and inherent transient operation. Molten salt-steam HX design lifespans exceed actual lifespans, and, as a consequence, designers overpredict plant profitability and operators neglect appropriate prescriptions to optimize these lifetimes. Here, this study refines HX lifespan estimates with data benchmarked against thermal-fluid mechanical modeling of stress and accumulated fatigue. Reduced-order thermal models of the molten salt-steam, shell-and-tube evaporator and superheater predict transient temperature profiles along the two HXs salt-steam flow paths. The modeled evaporator and superheater temperature profiles enable assessment of cyclic stresses within the HX tubesheets, where molten-salt HX failures are most common. Evaporator and superheater performance data from a current 110 MW elec commercial CSP plant provide a basis for validating the reduced-order HX models. HX life predictions derived from stochastic failure distributions serve as inputs for simulating and optimizing existing plant operations. The impact of the updated lifespans on overall plant revenue depends on operating scenarios. This study suggests that typical ramping rates for a CSP plant with a high-temperature Rankine cycle result in an evaporator and superheater life of approximately 10 and 25 years, respectively, compared to the design target of 30 years. Reduced HX lifespans decrease operational plant revenue on average by 4.6-5.1%. Furthermore, there may be as many as four HX replacements over the 30-year lifetime of the plant; and, purchase agreement loss due to failure to meet contractual production requirements can have ramifications that include the risk of bankruptcy.

14 SOLAR ENERGY↗

DOE OSTI2025

Changes in atmospheric circulation amplify extreme snowfall fueled by Arctic sea ice loss over high-latitude land

Arctic sea-ice retreat has been linked to increased winter precipitation and heavy snowfall over land, likely due to a combination of enhanced evaporation from ice-free Arctic marginal seas (AMS) and changes in atmospheric circulation. However, their relative roles and contributions remain uncertain. Here, we show that a greater proportion of AMS evaporative moisture reached high-latitude land during the cold seasons from 1980–1989 to 2012–2021. Atmospheric circulation changes added an additional 13 % increase in the AMS moisture contribution, accounting for 11 % of the total increase in AMS-sourced land precipitation. Notably, 46 % of the increase in AMS-sourced extreme snowfall is attributed to circulation-driven landward moisture transport, representing an 84 % increase beyond the effect of enhanced AMS evaporation alone. Further analysis indicates that both the rise in Arctic moisture and the atmospheric circulation shifts are primarily driven by anthropogenic forcing. These findings highlight how atmospheric circulation changes amplify extreme snowfall fueled by AMS evaporation, underscoring the synergistic effects of Arctic sea ice loss and circulation change on high-latitude winter precipitation.

Anthropogenic forcing↗

DOE OSTI2024

Component innovations for lower cost mechanical vapor compression

Despite significant capital and operating costs, mechanical vapor compression (MVC) remains the preferred technology for challenging brine concentration applications. Here, this work seeks to assess the dependence of MVC costs on feedwater salinity and desired water recovery and to quantify the value of improved component performance or reduced component costs for reducing the levelized cost of water (LCOW) of MVC. We built a cost optimization model coupling thermophysical, heat and mass transfer, and technoeconomic models to optimize and identify low cost MVC system designs as a function of feedwater salinity and water recovery. The LCOW ranges over 3.6 to 6.1 $\$$/m 3 for seawater feed salinities of 25–150 g/kg and water recoveries of 40–80 %. We then perform sensitivity analysis on parameter inputs to isolate irreducible costs and determine high value component innovation targets. The LCOW was most sensitive to evaporator material costs and performance, including the overall heat transfer coefficient in the evaporator. Process and material innovations such as polymer-composite evaporator tubes that reduce evaporator costs by 25 % without reducing heat transfer performance by more than 10 % would result in MVC cost reductions of 8 %.

42 ENGINEERING↗

DOE OSTI2019

Simulating a Mesoscale Convective System Using WRF with a New Spectral Bin Microphysics: Part 1: Hail vs Graupel

A modified Fast Spectral Bin Microphysics scheme (FSBM-2) embedded into the Weather Research and Forecasting Model (WRF) is used to simulate a mesoscale deep convective system observed during the Midlatitude Continental Convective Clouds Experiment (MC3E). FSBM-2 uses modified source codes as compared to the current FSBM (FSBM-1). In contrast to FSBM-1, FSBM-2 can simulate hail of several cm in diameter and includes additional processes such as spontaneous breakup of raindrops and aerosol regeneration by drop evaporation. It is shown that allowing large hail particles of diameters exceeding about 1 cm substantially increases the agreement between the simulated and observed squall-line structures both in the convective and stratiform regions. In contrast, if graupel particles are used to represent high-density hydrometeors in convective areas, the ratio of convective-to-stratiform areas diverges from the ratio seen in observations and maximum radar reflectivities are substantially underestimated. Analysis of snow size distributions in the stratiform area shows an important link between the ice crystals formed by homogeneous freezing in the convective area to ice particle number concentration in the stratiform region. Simulated raindrop size distributions in the stratiform area below the melting level from FSBM-2 show a good agreement with observations. The regeneration of CCN by droplet evaporation and detrainment of these CCN to the stratiform region increases the concentration of CCN there up to 8-9 km altitude; the additional CCN penetrate clouds and produce new droplets, leading to some convection intensification. Key words: Mesoscale convective systems (MCS), squall lines, spectral bin microphysics, large hail, in-cloud nucleation, size distributions of hydrometeors in MCS.

Shpund, Kobby↗

DOE OSTI2023

Effects of Wind Shear and Aerosol Conditions on the Organization of Precipitating Marine Stratocumulus Clouds

Here this study examines how wind shear affects precipitating marine stratocumulus clouds under different cloud droplet number concentrations (N d ). We performed a series of large eddy simulations (LES) of nocturnal marine stratocumulus clouds using Cloud Model 1 (CM1). The simulations show that N d is the dominant factor for cloud cellular organization in this cloud regime rather than wind shear. Low N d characterizes the open cellular structure with a high in-cloud liquid water path (LWP). When wind shear is increased, the cloud fraction tends to decrease along with LWP, suggesting the cloud top region is significantly influenced by the entrainment and mixing of dry air from the free troposphere. We also examine cold pools in open and closed cellular clouds. Open-cell clouds produce larger and deeper cold pools compared to closed-cell clouds. Interestingly, cold pools can exist without surface precipitation and are produced by evaporation of light precipitation (drizzle) below the cloud base with weak downdrafts. The evaporation of raindrops and drizzle play an important role in initiating new convection, particularly where colliding outflows occur downstream of the cloud. This secondary convection contributes to the development and maintenance of the cloud cellular organization and formation.

54 ENVIRONMENTAL SCIENCES↗