SEARCH · Engineering Papers

Results for “Field Validation”

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 1,009 records · Page 56

All sources NASA NTRS DOE OSTI

DOE OSTI2022

First application of a digital mirror Langmuir probe for real-time plasma diagnosis

For the first time, a digital Mirror Langmuir probe (MLP) has successfully sampled plasma temperature, ion saturation current, and floating potential together on a single probe tip in real time in a radio-frequency driven helicon linear plasma device. This is accomplished by feedback control of the bias sweep to ensure a good fit to I-V characteristics with a high frequency, high power digital amplifier and field-programmable gate array (FPGA) controller. Measurements taken by the MLP were validated by a low speed I-V characteristic manually collected during static plasma conditions. Plasma fluctuations, induced by varying the axial magnetic field (f̃ = 10 Hz), were also successfully monitored with the MLP. Further refinement of the digital MLP pushes it towards a turn-key system that minimizes the time to deployment and lessens the learning curve, positioning the digital MLP as a capable diagnostic for the study of low radio-frequency plasma physics. These demonstrations bolster confidence in fielding such digital MLP diagnostics in magnetic confinement experiments with high spatial and adequate temporal resolution such as edge plasma, scrape-off layer, and divertor probes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2014

Algorithm Development and Validation of CDOM Properties for Estuarine and Continental Shelf Waters Along the Northeastern U.S. Coast

An extensive set of field measurements have been collected throughout the continental margin of the northeastern U.S. from 2004 to 2011 to develop and validate ocean color satellite algorithms for the retrieval of the absorption coefficient of chromophoric dissolved organic matter (aCDOM) and CDOM spectral slopes for the 275:295 nm and 300:600 nm spectral range (S275:295 and S300:600). Remote sensing reflectance (Rrs) measurements computed from in-water radiometry profiles along with aCDOM() data are applied to develop several types of algorithms for the SeaWiFS and MODIS-Aqua ocean color satellite sensors, which involve least squares linear regression of aCDOM() with (1) Rrs band ratios, (2) quasi-analytical algorithm-based (QAA based) products of total absorption coefficients, (3) multiple Rrs bands within a multiple linear regression (MLR) analysis, and (4) diffuse attenuation coefficient (Kd). The relative error (mean absolute percent difference; MAPD) for the MLR retrievals of aCDOM(275), aCDOM(355), aCDOM(380), aCDOM(412) and aCDOM(443) for our study region range from 20.4-23.9 for MODIS-Aqua and 27.3-30 for SeaWiFS. Because of the narrower range of CDOM spectral slope values, the MAPD for the MLR S275:295 and QAA-based S300:600 algorithms are much lower ranging from 9.9 and 8.3 for SeaWiFS, respectively, and 8.7 and 6.3 for MODIS, respectively. Seasonal and spatial MODIS-Aqua and SeaWiFS distributions of aCDOM, S275:295 and S300:600 processed with these algorithms are consistent with field measurements and the processes that impact CDOM levels along the continental shelf of the northeastern U.S. Several satellite data processing factors correlate with higher uncertainty in satellite retrievals of aCDOM, S275:295 and S300:600 within the coastal ocean, including solar zenith angle, sensor viewing angle, and atmospheric products applied for atmospheric corrections. Algorithms that include ultraviolet Rrs bands provide a better fit to field measurements than algorithms without the ultraviolet Rrs bands. This suggests that satellite sensors with ultraviolet capability could provide better retrievals of CDOM. Because of the strong correlations between CDOM parameters and DOM constituents in the coastal ocean, satellite observations of CDOM parameters can be applied to study the distributions, sources and sinks of DOM, which are relevant for understanding the carbon cycle, modeling the Earth system, and to discern how the Earth is changing.

CDOM Properties↗

NASA NTRS2000

LIS Validation at The KSC-ER

Each year, thousands of lightning electric field disturbances are recorded and archived by the ground-based field mill (FM) network at the NASA Kennedy Space Center (KSC) and USAF Eastern Range (ER). The FM network has a range of several tens of kilometers, and a digital accuracy of 4 V/m. It has provided years of continuous lightning warning surveillance to KSC-ER space vehicle launch operations, and has undergone one major hardware upgrade since its inception in the early 1970s. Additional KSC lightning warning data is derived from a multistation radio time-of-arrival system called Lightning Detection and Ranging (LDAR). This system provides the location and space-time mapping of individual lightning channels (for both cloud and ground flashes). Additional lightning information for the KSC region is available from the National Lightning Detection Network (NLDN) and a 5-station local magnetic direction finder network. In this study, all of the above mentioned data are used to ground-validate data derived from the Lightning Imaging Sensor (LIS) onboard the Tropical Rainfall Measuring Mission (TRMM). The FM network can be used to retrieve the charges deposited in a lightning flash, provided the flash is within a few kilometers of the FM Network. Although it is rare to obtain a TRMM overpass of thunderstorms hat occur this close to the FM network, seven such storms have been found and examined in this study. We compare the times and locations of LIS optical pulses with the spatial-temporal character of the FM, LDAR, and magnetic direction finder data. We also inter-compare LIS optical pulse amplitude data with FM-derived charge magnitudes, number of LDAR radio sources, and peak current values from magnetic direction finder data. Generally speaking, LIS lightning locations and times agree favorably with the KSC ground-based systems for most cases, but little correlation appears to exist between optical pulse amplitude and any of charge, # LDAR sources, peak current), owing possibly to the effects of source complexity and/or cloud multiple scattering.

Koshak, William J.↗

DOE OSTI2022

Feasibility study of Mg storage in a bilayer silicene anode via application of an external electric field

With the goal of developing a Si-based anode for Mg-ion batteries (MIBs) that is both efficient and compatible with the current semiconductor industry, the current research utilized classical Molecular Dynamics (MD) simulation in investigating the intercalation of a Mg 2+ ion under an external electric field (E-field) in a 2D bilayer silicene anode (BSA). First principles density functional theory calculations were used to validate the implemented EDIP potentials. Our simulation shows that there exists an optimum E-field value in the range of 0.2–0.4 V Å –1 for Mg 2+ intercalation in BSA. To study the effect of the E-field on Mg 2+ ions, an exhaustive spread of investigations was carried out under different boundary conditions, including calculations of mean square displacement (MSD), interaction energy, radial distribution function (RDF), and trajectory of ions. Our results show that the Mg 2+ ions form a stable bond with Si in BSA. The effects of E-field direction and operating temperature were also investigated. In the X–Y plane in the 0°–45° range, 15° from the X-direction was found to be the optimum direction for intercalation. The results of this work also suggest that BSA does not undergo drastic structural changes during the charging cycles with the highest operating temperature being ~300 K

25 ENERGY STORAGE↗

NASA NTRS2024

Filament Wound Composite Analysis Using the NASA Multiscale Analysis Tool (NASMAT) and Finite Element Analysis

Fiber reinforced composite materials, owing to their tailorable thermomechanical and functional properties, allow one to produce a structure that is stronger, stiffer, and lighter than its metal counterpart while performing the same function, yielding a more efficient structure. This not only allows for the improvement of current technologies like aircraft structures, but also enables new technologies like gaseous hydrogen storage for mobility applications, which are otherwise impractical when manufactured using traditional metals due to weight and space restrictions or material embrittlement. However, the use of composites imposes greater design and manufacturing challenges on an engineer, since they are heterogenous, having a distinct structure across multiple length scale, behave generally anisotropically at the structural level and require complex manufacturing and processing methods. Capturing this complex behavior requires detailed numerical simulations, including the modeling of microstructural features like undulations, voids, and fiber alignment. In this paper, multiple repeating unit cells (RUCs), representing filament wound composites, are developed (via a script provided in the Appendix) and analyzed. The refinement of these RUCs is varied, and the analyses are performed using both the Abaqus finite element software and the NASA Multiscale Analysis Tool (NASMAT). A study is undertaken to compare the predicted effective elastic properties of the wound RUC to a laminate representation of the wound RUC, which neglects the undulations. Additionally, two different sets of periodic boundary conditions (PBCs) have been examined. One approximates the real boundary conditions using a standard approach and the other represents the PBCs exactly through the use of an offset. Lastly, a comparison of the local elastic stress fields is made among the models and approaches. Since wound structures are often approximated as laminated structures, it is important to understand the degree to which this assumption is valid, namely by first comparing the elastic constants and local elastic fields. This will provide, on the one hand, information concerning the bulk mechanical behavior and, on the other hand, insights concerning local load distributions and likely damage initiation sites.

Composites↗

NASA NTRS1997

Snow Parameters Derived from Microwave Measurements During the BOREAS Winter Field Campaign

Passive microwave data have been used to infer the snow-covered area and snow water equivalent (SWE) over forested areas, but the accuracy of these retrieved snow parameters cannot be easily validated for heterogeneous vegetated regions. The Boreal Ecosystem-Atmosphere Study Winter Field Campaign provided the opportunity to study the effect of boreal forests on snow parameter retrieval in detail. Microwave radiometers (18, 37, and 92 GHz) were flown on board the Canadian National Aeronautical Establishment's Twin Otter. Flight lines covered both the southern study area near Prince Albert and the northern study area near Thompson, Canada. During the 1994 winter campaign, extensive ground-based snow cover information, including depth, density, and grain size, was collected along most of the flight lines, jointly by U.S. and Canadian investigators. Satellite data collected by the special sensor microwave imager are also used for comparison. Preliminary results reconfirmed the relationship between microwave brightness temperature and SWE. However, the effect of forest cover observed by the aircraft sensors is different from that of the satellite observations. This is probably due to the difference in footprint averaging. There were also several flight lines flown over Candle Lake and Waskesiu Lake to assess lake ice signatures. Preliminary results show the thickness of the lake ice may be inferred from the airborne microwave observations. The microwave signature relationship between lake ice and snow matches the results from radiative transfer calculations.

Chang, A. T. C.↗

DOE OSTI2023

A sharp immersed method for 2D flow-body interactions using the vorticity-velocity Navier-Stokes equations

Immersed methods discretize boundary conditions for complex geometries on background Cartesian grids. Here, this makes such methods especially suitable for two-way coupled flow-body problems, where the body mechanics are partially driven by hydrodynamic forces. However, for the vorticity-velocity form of the Navier-Stokes equations, existing immersed geometry discretizations for two-way coupled problems only achieve first order spatial accuracy near solid boundaries. Here we introduce a sharp-interface approach based on the immersed interface method to handle the one- and two-way coupling between an incompressible flow and one or more rigid bodies using the 2D vorticity-velocity Navier-Stokes equations. Our main contributions are three-fold. First, we develop and analyze a moving boundary treatment for sharp immersed methods that can be applied to PDEs with implicitly defined boundary conditions, such as those commonly imposed on the vorticity field. Second, we develop a two-way coupling methodology for the vorticity-velocity Navier-Stokes equations based on control-volume momentum balance that does not require the pressure field. Third, we show through extensive testing and validation that our resulting flow-body solver reaches second-order accuracy for most practical scenarios, and provides significant efficiency benefits compared to a representative first-order approach.

42 ENGINEERING↗

DOE OSTI2021

ATHENA: A unique radiation environment platform at the National Ignition Facility

This paper describes the ATHENA platform, an energy tuning assembly, which was developed to spectrally shape the National Ignition Facility (NIF) deuterium–tritium fusion neutron source to a thermonuclear (fusion) plus prompt fission neutron spectrum with a capability to act as a short-pulse neutron source. This unique, otherwise inaccessible radiation environment complements existing experimental facilities and capabilities. Here, the flexible ATHENA irradiation positions were modeled using an ensemble of Monte Carlo simulations with stochastic sampling of the nuclear cross-sections to characterize the radiation environments and uncertainty for the platform. Validation of the internal neutron spectrum produced from fielding ATHENA at NIF occurred through neutron flux unfolding with 20 measured activation products.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2025

Identifying High Ionic Conductivity Compositions of Ionic Liquid Electrolytes Using Features of the Solvation Environment

Binary mixtures of ionic liquids with molecular solvents are gaining interest in electrochemical applications due to the improvement in their performance over neat ionic liquids. Dilution with suitable molecular solvents can reduce the viscosity and facilitate faster diffusion of ions, thereby yielding substantially higher ionic conductivity than that for a pure ionic liquid. Although viscosity and diffusion coefficients typically behave as monotonic functions of concentration, ionic conductivity often passes through a peak value at an optimum molar ratio of the molecular solvent to the ionic liquid. The ionic conductivity maximum is generally explained in terms of a balance between the ease of charge transport and the concentration of the charge carriers. In this work, fluctuation in the local environment surrounding an ion is invoked as a plausible explanation for the ionic conductivity mechanism with a binary mixture of 1-ethyl-3-methylimidazolium tetrafluoroborate and ethylene glycol as an example. The magnitude of the dynamism in the local environment is captured by measuring the spatial and temporal features of the solvation environment. Standard deviation in the number of ions in the solvation environment serves as a spatial feature, while the cage correlation lifetimes for oppositely charged ions within the first solvation shell serve as a temporal feature. Large standard deviations in the cluster ion population and short cage correlation lifetimes are indicators of highly dynamic ionic environment at the molecular level and consequently yield high ionic conductivity. Such compositions were found to be in good agreement with the optimum ionic liquid mole fractions obtained through experimental measurement. Short cage correlation lifetimes enable the identification of optimum mixture compositions using simulation trajectories significantly shorter than those required to implement the Nernst–Einstein or Einstein formalisms for calculating ionic conductivity. We validated the applicability of this approach across force fields and in six ionic liquid-molecular solvent electrolytes formed with combination of cations, anions, and solvents. We offer a computationally efficient approach of screening ionic liquid-molecular solvent binary mixture electrolytes to identify molar ratios that yield high ionic conductivity.

25 ENERGY STORAGE↗

DOE OSTI2019

Coherent elastic neutrino-nucleus scattering on Ar 40 from first principles

Coherent elastic neutrino scattering on the 40 Ar nucleus is computed with coupled-cluster theory based on nuclear Hamiltonians inspired by effective field theories of quantum chromodynamics. Our approach is validated by calculating the charge form factor and comparing it to data from electron scattering. We make predictions for the weak form factor, the neutron radius, and the neutron skin and estimate systematic uncertainties. The neutron-skin thickness of 40 Ar is consistent with results from density-functional theory. Precision measurements from coherent elastic neutrino-nucleus scattering could potentially be used to extract these observables and help to constrain nuclear models.

39 ≤ A ≤ 58↗

DOE OSTI2022

Effective bias expansion for 21-cm cosmology in redshift space

A near-future detection of the 21-cm signal from the epoch of reionization will provide unique opportunities to probe the underlying cosmology, provided that such cosmological information can be extracted with precision. To this end, we further develop effective field theory (EFT) inspired techniques for the 21-cm brightness temperature field during the epoch of reionization, incorporating renormalized bias and a treatment of redshift space distortions. Notably, we confirm that in redshift space, measures of the 21-cm brightness, e.g., the power spectrum, should have irreducible contributions that lack a bias coefficient and therefore contain direct, astrophysics-free information about the cosmological density field; in this work, we study this effect beyond linear order. To validate our theoretical treatment, we fit the predicted EFT Fourier-space shapes to the thesan suite of hydrodynamical simulations of reionization at the field level, where the considerable number of modes prevents overfitting. We find agreement at the level of a few percent between the 21-cm power spectrum from the EFT fits and simulations over the wave number range k ≲ 0.8 h/Mpc and neutral fraction x HI ≳ 0.4, which is imminently measurable by the Hydrogen Epoch of Reionization Array and future experiments. Furthermore, the ability of the EFT to describe the 21-cm signal extends to simulations that have different astrophysical prescriptions for reionization as well as simulations with interacting dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Impact of strong magnetization in cylindrical plasma implosions with applied B-field measured via x-ray emission spectroscopy

Magnetization is a key strategy for enhancing inertial fusion performance, though accurate characterization of magnetized dense plasmas is needed for a better comprehension of the underlying physics. Measured spectra from imploding Ar-doped D 2 -filled cylinders at the OMEGA laser show distinctive features with and without an imposed magnetic field. A multizone spectroscopic diagnosis leads to quantitative estimates of the plasma conditions, namely revealing a 50% core temperature rise at half mass density when a 30-T seed field is applied. Concurrently, experimental spectra align well with predictions from extended-magnetohydrodynamics simulations, providing strong evidence that the attained core conditions at peak compression are consistent with the impact of a 10-kT compressed field. These results pave the way for the validation of magnetized transport models in dense plasmas and for future magnetized laser implosion experiments at a larger scale. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Machine Learning Enabled Position Detection for 6.78 MHz UAV Wireless Power Transfer System

This paper presents a novel supervised machine learning (SML) approach for accurate position detection of the receiver coil in wireless power transfer (WPT) systems using only secondary-side electrical measurements, with applications in autonomous unmanned aerial vehicle (UAV) charging. The proposed method trains a supervised learning model to map measured secondary-side voltage and current features to the receiver’s spatial position with high precision. This enables an autonomous UAV to determine its location relative to the primary coil center, the optimal position for maximizing wireless charging efficiency. The sensing method is fully integrated into a standard WPT system, utilizing the same primary and secondary coils for both power transfer and position detection, thereby eliminating additional sensing hardware. The use of a 6.78 MHz operating frequency enhances positional sensitivity, as high-frequency near-field electromagnetic fields respond strongly to small spatial variations. Experimental validation is performed on a 30 W scaled prototype featuring a 210 mm × 140 mm primary coil, a 50 mm × 80 mm receiver coil, and a 15 mm air gap. Results demonstrate reliable position estimation and a strong correlation between predicted position and optimal coil alignment. This integrated framework unifying position detection and wireless charging offers a promising foundation for future autonomous electric vertical takeoff and landing (eVTOL) systems, enabling compact, hardware-efficient, and high-accuracy charging solutions.

Colak, Kerim [New York University]↗

DOE OSTI2023

Status of the MQXFB Nb 3 Sn Quadrupoles for the HL-LHC

Here, the cold powering test of the first two prototypes of the MQXFB quadrupoles (MQXFBP1, now disassembled, and MQXFBP2), the Nb 3 Sn inner triplet magnets to be installed in the HL-LHC, has validated many features of the design, such as field quality and quench protection, but has found performance limitations. In fact, both magnets showed a similar phenomenology, characterized by reproducible quenches in the straight part inner layer pole turn, with absence of training and limiting the performance at 93% (MQXFBP1) and 98% (MQXFBP2) of the nominal current at 1.9 K, required for HL-LHC operation at 7 TeV. Microstructural inspections of the quenching section of the limiting coil in MQXFBP1 have identified fractured Nb 3 Sn sub-elements in strands located at one specific position of the inner layer pole turn, allowing to determine the precise origin of the performance limitation. In this paper we outline the strategy that has been defined to address the possible sources of performance limitation, namely coil manufacturing, magnet assembly and integration in the cold mass.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

On the Development of Compact HTS Coil Modules for Large Bore High Field Superconducting Magnets

Lower cost, high current density superconducting coil modules producing higher magnetic fields and cooled affordably are crucial for obtaining cost-effective, compact commercial fusion reactors. Accessibility to low cost, higher field magnets (>30 T) is also critical for the discovery of new quantum phenomena in materials, cosmic frontier and other topics in basic science research. The Princeton Plasma Physics Laboratory (PPPL) is working with Princeton University to develop unique large bore, compact superconducting magnets to support science experiments including the development of new instrumentation for condensed matter physics and Axion dark matter search in the cosmic frontier. Core elements of these experiments are unique for access to lower cost, simple fabrication of compact superconducting magnets that can be cooled affordably, while integrated with dedicated science instruments. Conductor qualification and coil design concepts are discussed in support of needs for these experiments. PPPL has the unique expertise and experimental facilities to design, construct and test subscale coil modules for these projects. Compact high field coil modules were fabricated and tested to validate coil design concepts and coil performance. Finally, the design and model coil integration challenges are discussed to identify performance risks and demonstrate feasibility for deploying full scale large bore compact superconducting magnets for cost effective operations of multiple laboratory experiments.

axion dark matter↗

DOE OSTI2021

Hardware Design of Medium Voltage SiC-Based Modular Multilevel Converters for Grid-Tied Applications

This paper presents a design of a 10 kV silicon carbide (SiC) MOSFET based modular multilevel converter (MMC). The MMC functions as the interconnection between a 24 kV dc bus and a 13. SkV distribution grid. This paper focuses on the hardware design of the MMC and its submodules. The gate driver, auxiliary power supply and bus bur in submodules are major components that need to meet high dv/dt immunity and high insulation voltage requirements. The gate driver achieves over 200 kV/us common mode transient immunity (CMTI). Its associated auxiliary power supply is based on an on-board transformer with over 30 kV insulation capability and 1.4 pF barrier capacitance. The busbar achieves less than 11 kV/cm electric field in the air to avoid partial discharges. Validation tests of the MMC submodule were performed by double pulse tests at 5 kV and buck converter tests at3 kV, 52A (peak). Reliable operations were achieved in validation tests.

auxiliary power supply↗

DOE OSTI2026

Asi Nuclear Energy Sensors Data Portal Chatbot And Data Structuring Tool

The Idaho National Laboratory (INL) is advancing the development of an AI-powered chatbot and data structuring tool specifically designed to accelerate data mining processes for sensor-related information and seamlessly integrate the results into the ASI Sensors Data Portal (https://nes.energy.gov/). By doing so, the software aims to enhance the accessibility, usability, and organization of sensor data for nuclear energy applications. The software initial phase focuses on retrieving comprehensive datasets, prioritizing the past five years of publicly available information from the Office of Scientific and Technical Information (OSTI). These datasets will be meticulously processed to ensure compatibility, employing cleaning and preprocessing steps to eliminate irrelevant, incomplete, or corrupted information, thus establishing a robust foundation for subsequent AI use. The data will serve as the backbone for training an AI model and chatbot, which will act as an interactive tool enabling users to ask complex, context-specific questions and receive accurate, validated answers derived from constrained literature. In parallel, the project incorporates a data structuring process supported by AI to organize sensor information from multiple sources into a standardized format. This structured data will include detailed sensor specifications, such as measurement range, applications, accuracy, and operating conditions, generated and documented with AI. These specifications will be systematically integrated into the sensor portal. To maintain the highest levels of accuracy and relevance, all AI-generated outputs will be reviewed and validated by subject matter experts (SMEs), with additional fields or parameters added as needed. Future stages of the project aim to expand the dataset beyond OSTI to include other sources and potentially incorporate unclassified controlled information (UCI) with restricted access protocols to address security and confidentiality requirements.

Mapes, NormanJ. [Idaho National Laboratory (INL), ↗

DOE OSTI2020

Experimental Investigation on Separation Behavior of Heavy-Oil Emulsion for Polymer Flooding on Alaska North Slope

The first-ever field pilot on Alaska North Slope (ANS) to validate using polymer floods for heavy-oil enhanced oil recovery is currently ongoing. One of the major concerns of the operator is the effect of polymer on oil/water-separation efficiency after polymer breakthrough. This work investigates the influence of polymer on the separation behavior of heavy-oil emulsions and evaluates the performance of emulsion breakers (EBs). In this study, two types of heavy-oil emulsions were prepared and tested at 20 and 50% water cut (WC), respectively. The bottle test method was used in the experiments, in which the separated water volume with time, the separated water quality, and the volume fraction of phases were recorded. Results showed that polymer accelerated the oil/water separation acting as an emulsion inhibitor at 20% WC but tended to impede the water separation at 50% WC. Regardless of WC, polymer resulted in poor water quality and the formation of a stable intermediate oil in water (o/w) emulsion, because of the increased viscosity of the water phase. The performance of EBs showed a complex dependency on the WC, the type of demulsifier and dosage, and the polymer concentration. Despite the varied conditions encountered in the heavy-oil/water/polymer/demulsifier system, a compound EB achieved satisfactory demulsification performance, showing the highest potential for deployment in the current ANS polymer flooding pilot. Here, we systematically studied the potential influence of polymer breakthrough on the separation behavior of heavy-oil emulsion on ANS for the first time. The findings of this study will provide practical guidance in advance for produced fluid treatment of the ongoing first-ever polymer flooding pilot on ANS.

54 ENVIRONMENTAL SCIENCES↗