Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “particle improvement”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Implementation of A New Microwave Scattering Database and A Forward Model for Active Microwave Sensors in CRTM

Radiative transfer models are extensively used for the assimilation of satellite observations into NWP models as well as retrieving geophysical products from satellite measurements. CRTM is a community model developed by NOAA JCSDA and widely used for different purposes requiring RT calculations. CRTM requires bulk scattering lookup tables in order to perform all-sky RT calculations. However, the current CRTM lookup tables for microwave frequencies were generated based on the Mie theory by assuming spherical frozen particles. The scattering lookup tables generated using the DDA technique has shown to largely improve the RT scattering calculations in the MW region. This presentation targets (i) the implementation and validation of a DDA database that was originally developed for the ARTS RT model into CRTM, and (ii) developing the CRTM active sensor module that takes advantage of the backscattering coefficients computed using the DDA method. The DDA database only provides single scattering properties of different habits, while CRTM requires bulk scattering properties. The CRTM cloud coefficients were previously generated based on the effective radius for representing the size of the particles. However, effective radius is neither measurable nor provided by the NWP models, thus need to be estimated from other geophysical variables such as water content. Therefore, in addition to calculating the CRTM bulk scattering properties from the DDA single scattering database, the CRTM was also largely modified to use cloud water content (kg.m-3), instead of effective radius, for performing the interpolation over size/mass of the particles. CRTM already requires water content as input, thus no extra variables are required for performing scattering calculations using the new ARTS DDA database. The CRTM scattering modules search for effective radius in cloud coefficient files and will use the cloud water content if the effective radius dimension is not found in the cloud coefficient files. Figure 1 shows the CRTM simulated brightness temperatures computed using different cloud coefficients versus ATMS observed values over Hurricane Irma on September 7, 2017 at 18:00 UTC. We used all the cloud water content values included in ERA5 with default CRTM/DDA habits for water, rain, snow, ice, hail, and graupel. ERA5 does not provide separate water content values for ice, hail, and graupel, thus the ice water content values were divided between ice, hail, and graupel clouds similar to what was explained in the previous section. In channels with a frequency lower than 90 GHz, emission from water and rain clouds can compensate for cloud scattering so that cloud contaminated Tbs are larger than corresponding clear sky Tbs. The DDA simulations for channels 1-7 largely perform better than the Mie simulations. The DDA simulations show a mix of small negative and positive simulated minus observed values, while the Mie results show large negative biases. The weighting functions for some of the ATMS temperature sounding channels (channels 9-15) peak mostly above the clouds, therefore the measured Tbs become less sensitive to clouds so that the results of both Mie and DDA become very similar. The Mie lookup tables generate excessive scattering for channel 16, but not enough scattering for the water vapor channels. In the specific case of Hurricane Maria, the DDA lookup tables do not generate enough scattering for channel 16, but the DDA results are much more consistent with observations for water vapor channels than for channel 16. It should be noted that the results may vary if we use other habits to represent snow, hail, and graupel in the DDA simulations. Although these results clearly show the advantage of the DDA database over the Mie dataset, different error sources such as error in the observations, displacement of clouds in the ERA5 reanalysis, and also lack of convective clouds or in general errors in the input atmospheric and cloud profiles contribute to the differences between the simulated and observed values. Aside from the improvements in the simulations, a major advantage of the new dataset is a large number of habits that can be used to tune the data assimilation systems to perform well in different weather conditions.

Isaac Moradi↗
NASA NTRS2014

Hot Hydrogen Testing of Tungsten-Uranium Dioxide (W-UO2) CERMET Fuel Materials for Nuclear Thermal Propulsion

CERMET fuel materials are being developed at the NASA Marshall Space Flight Center for a Nuclear Cryogenic Propulsion Stage. Recent work has resulted in the development and demonstration of a Compact Fuel Element Environmental Test (CFEET) System that is capable of subjecting depleted uranium fuel material samples to hot hydrogen. A critical obstacle to the development of an NCPS engine is the high-cost and safety concerns associated with developmental testing in nuclear environments. The purpose of this testing capability is to enable low-cost screening of candidate materials, fabrication processes, and further validation of concepts. The CERMET samples consist of depleted uranium dioxide (UO2) fuel particles in a tungsten metal matrix, which has been demonstrated on previous programs to provide improved performance and retention of fission products1. Numerous past programs have utilized hot hydrogen furnace testing to develop and evaluate fuel materials. The testing provides a reasonable simulation of temperature and thermal stress effects in a flowing hydrogen environment. Though no information is gained about radiation damage, the furnace testing is extremely valuable for development and verification of fuel element materials and processes. The current work includes testing of subscale W-UO2 slugs to evaluate fuel loss and stability. The materials are then fabricated into samples with seven cooling channels to test a more representative section of a fuel element. Several iterations of testing are being performed to evaluate fuel mass loss impacts from density, microstructure, fuel particle size and shape, chemistry, claddings, particle coatings, and stabilizers. The fuel materials and forms being evaluated on this effort have all been demonstrated to control fuel migration and loss. The objective is to verify performance improvements of the various materials and process options prior to expensive full scale fabrication and testing. Post test analysis will include weight percent fuel loss, microscopy, dimensional tolerance, and fuel stability.

Hickman, Robert↗
DOE OSTI2026

Performance of the LHCb muon detector in Run 3

In Run 3 of the LHC, the instantaneous luminosity at the LHCb interaction point has been increased by a factor of five, from to $\mathscr{L}$ = 4 x 10 32 cm -2 s -2 to $\mathscr{L}$ = 2 x 10 33 cm -2 s -1 . Several hardware interventions, including a complete overhaul of the readout electronics, have been carried out on the muon detector. The muon identification algorithms in the software trigger were improved with the aim of ensuring Run 2 performance under a higher particle rate. The operation and calibration of the upgraded muon detector are presented. The muon detection efficiency and muon identification performance are evaluated on data calibration samples collected during the year 2024. A muon identification efficiency above 90% with sub-percent hadron misidentification probability is achieved by exploiting the pattern of hits in the muon detector.

Gaseous detectors↗
DOE OSTI2025

Towards improved neutral exhaust in the HSX stellarator

To improve our knowledge of neutral exhaust in non-resonant divertors, a particle exhaust study is performed for the Helically Symmetric eXperiment (HSX). The chaotic magnetic field topology in the plasma edge is examined with connection lengths and Lyapunov exponents. Plasma flux tubes are identified and these flux tubes deposit particles onto plasma facing components. This helps determine the areas where neutrals are generated and provides guidance on where to put baffles. As shown using EMC3-EIRENE simulations, the introduction of baffles can help increase the neutral pressure at and away from the strike line locations. Furthermore, it is demonstrated that baffle design cannot only be based on placing structures at locations with small connection lengths.

Baffles↗
DOE OSTI2021

Collective Nanoparticle Dynamics Associated with Bridging Network Formation in Model Polymer Nanocomposites

The addition of nanoparticles (NPs) to polymers is a powerful method to improve the mechanical and other properties of macromolecular materials. Such hybrid polymer–particle systems are also rich in fundamental soft matter physics. Among several factors contributing to mechanical reinforcement, a polymer-mediated NP network is considered to be the most important in polymer nanocomposites (PNCs). Here, we present an integrated experimental–theoretical study of the collective NP dynamics in model PNCs using X-ray photon correlation spectroscopy and microscopic statistical mechanics theory. Silica NPs dispersed in unentangled or entangled poly(2-vinylpyridine) matrices over a range of NP loadings are used. Static collective structure factors of the NP subsystems at temperatures above the bulk glass transition temperature reveal the formation of a network-like microstructure via polymer-mediated bridges at high NP loadings above the percolation threshold. The NP collective relaxation times are up to 3 orders of magnitude longer than the self-diffusion limit of isolated NPs and display a rich dependence with observation wavevector and NP loading. A mode-coupling theory dynamical analysis that incorporates the static polymer-mediated bridging structure and collective motions of NPs is performed. It captures well both the observed scattering wavevector and NP loading dependences of the collective NP dynamics in the unentangled polymer matrix, with modest quantitative deviations emerging for the entangled PNC samples. Here, we identify an unusual and weak temperature dependence of collective NP dynamics, in qualitative contrast with the mechanical response. Hence, the present study has revealed key aspects of the collective motions of NPs connected by polymer bridges in contact with a viscous adsorbing polymer medium and identifies some outstanding remaining challenges for the theoretical understanding of these complex soft materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Rational design of mechanically robust Ni-rich cathode materials via concentration gradient strategy

Mechanical integrity issues such as particle cracking are considered one of the leading causes of structural deterioration and limited long-term cycle stability for Ni-rich cathode materials of Li-ion batteries. Indeed, the detrimental effects generated from the crack formation are not yet entirely addressed. Here, applying physicochemical and electrochemical ex situ and in situ characterizations, the effect of Co and Mn on the mechanical properties of the Ni-rich material are thoroughly investigated. As a result, we successfully mitigate the particle cracking issue in Ni-rich cathodes via rational concentration gradient design without sacrificing the electrode capacity. Our result reveals that the Co-enriched surface design in Ni-rich particles benefits from its low stiffness, which can effectively suppress the formation of particle cracking. Meanwhile, the Mn-enriched core limits internal expansion and improve structural integrity. The concentration gradient design also promotes morphological stability and cycling performances in Li metal coin cell configuration.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Upgrade of the Lyman-alpha diagnostic system on DIII-D for main chamber edge neutral studies

The LLAMA (Lyman Alpha Measurement Apparatus) pinhole camera diagnostic had previously been deployed on DIII-D to measure radial profiles of the Lyman-α (Ly-α) deuterium neutral line brightness across the plasma boundary in the lower chamber to infer neutral deuterium density and ionization rate profiles. This system has recently been upgraded with a new diagnostic head, named ALPACA, that also encloses two pinhole cameras and duplicates the LLAMA views in the upper chamber. Similar to LLAMA, ALPACA provides two times 20 lines of sight, viewing the plasma edge on the inboard and outboard sides with a radial resolution of ~2.5 cm (FWHM) and an effective time resolution of ~1 ms that allows for the investigation of inter-ELM dynamics. The extended Ly-α system provides better coverage to study neutrals in experiments with various plasma shapes utilizing both the upper and lower divertors. Furthermore, post-campaign calibration of the LLAMA diagnostic has successfully been demonstrated for the first time. This was facilitated by various upgrades to the calibration set-up and detailed measurements of the emissivity distribution of the Ly-α calibration source using a pinhole collimator. It was found that the sensitivity of the inboard LLAMA pinhole camera was reduced by a factor of 2.0 ± 0.2 over the course of six months of plasma operation in 2021. In conclusion, the upgraded Ly-α system, equipped with improved absolute calibration, will provide key input for neutral fueling and pedestal particle transport studies and for 2D edge transport code validation on the DIII-D tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Recovery of Scandium from Bauxite Residue Using Sulfation Roasting-Assisted Water Leaching Process

Here, this study presents a method for recovering scandium from Indonesian bauxite residue (red mud) through a sulfation roasting-assisted water leaching process. Ammonium sulfate, (NH 4 ) 2 SO 4 , was chosen as the sulfate source for its lower melting point, less corrosivity, decreased toxicity, and the potential for ammonia recovery. A modified mixing system was introduced to improve the distribution of (NH 4 ) 2 SO 4 on red mud particles by combining them in water. A series of experiments were conducted to optimize the sulfation roasting process by examining the effects of the (NH 4 ) 2 SO 4 /red mud ratio, temperature, and roasting time. Subsequently, water leaching was optimized while considering the effect of leaching time and liquid-to-solid (L/S) ratio. The leaching of scandium from the roasted red mud showed rapid kinetics, even when using only water. Under the roasting conditions of 1/1 (g/g) (NH 4 ) 2 SO 4 /red mud ratio, 700°C, 60 min roasting time, and the leaching conditions of 65°C, 60 min, and 10 mL/g (L/S), over 90% of scandium was recovered, along with 64.4% Al leaching and below 8% Fe leaching. These results were further corroborated through comprehensive analyses such as x-ray diffraction (XRD) and scanning electron microscopy-energy dispersive X-ray spectroscopy (SEM-EDS). This study offers insights into an alternative, more environmentally friendly method for Sc recovery from red mud, promoting the sustainable use of hazardous waste like red mud.

ammonium sulfate↗
DOE OSTI2024

Turbulence-reduced high-performance scenarios in Wendelstein 7-X

In the Wendelstein 7-X (W7-X) stellarator, turbulence is the dominant transport mechanism in most discharges. This leads to a 'clamping' of ion temperature over a wide range of heating power, predominantly flat density profiles where hollow profiles driven by neoclassical thermo-diffusion would be expected and by rapid impurity transport in injection experiments. Significantly reduced turbulent transport is observed in the presence of strong core density gradients found transiently after core pellet injection and irregularly after boronisation or boron pellet injection. Density peaking is also achieved in a controlled manner in purely neutral beam heated discharges where particle transport analysis reveals an abrupt reduction in the main-ion particle flux leading to significant density profile peaking not explained by the NBI particle source alone. The plasmas exhibit a heat diffusivity of around $\chi = 0.25 \pm 0.1\,\mathrm{m}^2\ \mathrm{s}^{-1}$ at mid radius, a factor of around 4 lower than ECRH dominated discharges. Despite the improved confinement, the achieved ion temperature is limited by broader heat deposition and the lower power-per-particle given the higher density. This is overcome with limited reintroduction of ECRH power, where the low heat diffusivity diffusivity is maintained, the density rise supressed and ion temperatures above the clamping limit are achieved. The applicability of these plasmas for a high performance scenario on transport relevant time scales is assessed, including initial predictions for planned heating upgrades of W7-X, based on a range of assumptions about particle transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Graph Neural Network for Object Reconstruction in Liquid Argon Time Projection Chambers

The Exa.TrkX project presents a graph neural network (GNN) technique for low-level reconstruction of neutrino interactions in a Liquid Argon Time Projection Chamber (LArTPC). GNNs are still a relatively novel technique, and have shown great promise for similar reconstruction tasks in the Large Hadron Collider (LHC). Graphs describing particle interactions are formed by treating each detector hit as a node, with edges describing the relationships between hits. We utilise a multi-head attention message passing network which performs graph convolutions in order to label each node with a particle type.We present an updated variant of our GNN architecture, with several improvements. After testing the model on more realistic simulation with regions of unresponsive wires, the target was modified from edge classification to node classification in order to increase robustness. Removing edges as a classification target opens up a broader possibility space for edge-forming techniques; we explore the model’s performance across a variety of approaches, such as Delaunay triangulation, kNN, and radius-based methods. We also extend this model to the 3D context, sharing information between detector views. By using reconstructed 3D spacepoints to map detector hits from each wire plane, the model naively constructs 2D representations that are independent yet fully consistent.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Mesoscale Effects of Composition and Calendering in Lithium-Ion Battery Composite Electrodes

Macrohomogeneous battery models are widely used to predict battery performance, necessarily relying on effective electrode properties, such as specific surface area, tortuosity, and electrical conductivity. While these properties are typically estimated using ideal effective medium theories, in practice they exhibit highly non-ideal behaviors arising from their complex mesostructures. In this paper, we computationally reconstruct electrodes from X-ray computed tomography of 16 nickel–manganese–cobalt-oxide electrodes, manufactured using various material recipes and calendering pressures. Due to imaging limitations, a synthetic conductive binder domain (CBD) consisting of binder and conductive carbon is added to the reconstructions using a binder bridge algorithm. Reconstructed particle surface areas are significantly smaller than standard approximations predicted, as the majority of the particle surface area is covered by CBD, affecting electrochemical reaction availability. Finite element effective property simulations are performed on 320 large electrode subdomains to analyze trends and heterogeneity across the electrodes. Significant anisotropy of up to 27% in tortuosity and 47% in effective conductivity is observed. Electrical conductivity increases up to 7.5× with particle lithiation. We compare the results to traditional Bruggeman approximations and offer improved alternatives for use in cell-scale modeling, with Bruggeman exponents ranging from 1.62 to 1.72 rather than the theoretical value of 1.5. We also conclude that the CBD phase alone, rather than the entire solid phase, should be used to estimate effective electronic conductivity. This study provides insight into mesoscale transport phenomena and results in improved effective property approximations founded on realistic, image-based morphologies.

25 ENERGY STORAGE↗
DOE OSTI2020

Discovering the nature of neutrinos: How nuclear theory will advance this grand endeavor. Final report

Why is there more matter than antimatter in the visible universe? Are neutrinos their own antiparticles? These are some of the most profound questions in science. The detection of a (0vββ) neutrinoless doublebeta decay and the measurement of its half-life is one of the most important experimental endeavors in science because it will answer these questions. Paramount to these measurements are accurate predictions with quantified uncertainties of the 0vββ-decay nuclear transitions, especially for the germanium-76 and xeon-136 nuclei, two candidate materials for building the next generation of DOE funded 0vββ-decay experiments. These predictions will guide the Office of Science’s down selection of experiments as well as aid in determining how much material is needed by these experiments. We studied a specific difficult to assess source of uncertainty that will impact coming microscopic predictions, uncertainty which arises from an incomplete yet systematically improvable model of interactions between nucleons and of nucleons with the other particles involved in decay process, referred to as the theory uncertainty. We employed a surrogate model of the germanium-76 nucleus, and its decay products, just detailed enough to adequately mimic the decaying nucleons need for a reasonable measure of this uncertainty and nothing more, avoiding the immense computational expense required in quantifying uncertainties of the full calculation. We found that microscopic calculations of this process can be expected to achieve relative theory uncertainties around 30%, including an unconstrained contribution that had the potential to spoil microscopic calculations. This results in a nearly order of magnitude on previously estimated uncertainties for this rare decay process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Solid State Detectors Dedicated to Fast Timing

Development and implementation of 4D-trackers will play a key role at future machines. These detectors will reduce backgrounds, improve track reconstruction and Level-1triggering, and enable new capabilities for particle identification and reconstruction of long-lived particles. All of these pose unique challenges and opportunities to detector design. In this talk I will overview the currently available technologies and areas of active R&D, and proposed disruptive technologies to achieve the goals for future experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Progress Toward the First Search for Bound Neutron Oscillation into Antineutron in a Liquid Argon TPC

This note presents current progress for a neutron-antineutron oscillation ($n–\overline{n}$) search in MicroBooNE paving the way for the first search analysis of such process in a Liquid Argon Time Projection Chamber (LArTPC). Convolutional Neural Network (CNN) and Boosted Decision Tree (BDT) algorithms were used to select signal $n–\overline{n}$ events over cosmogenic backgrounds. The CNN-only, BDT-only, and the combined (CNN+BDT) methods were demonstrated on the Monte-Carlo signal and background events. Validation of the CNNonly and the BDT-only methods was carried out on a small dataset of MicroBooNE Run1 off-beam data, setting the starting point toward further improvement of the analysis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Interface Box for Dual Power Amplifier Modulators in Main Injector

The Proton Improvement Plan II (PIP-II) aims to enhance the current Fermilab particle accelerator complex to intensify and accelerate more beam. To achieve this, a second power amplifier will be added to the Main Injector (MI) cavities. However, the control unit could only read signals coming from the previous power supplies, so the interface box was designed to facilitate the issue.

Ayala, Diana↗
DOE OSTI2024

Physics Design of Proton/Ion SRF Accelerators

Superconducting Radio Frequency (SRF) technology has become foundation of the modern high energy particle accelerators. The technology enables a higher accelerating field with minimal power dissipation and a stronger compact magnetic field. These advancements have broadened the operational horizon of the particle accelerators making a higher duty factor and Continuous Wave (CW) beam operations both feasible and economically viable. The intensity frontier research founded on generation of intense beam of exotic particles including but not limited to neutrinos, muons, kaons and neutrons as well as practical applications involving transmutation of nuclear reactor waste, rare isotopes generation, sub-critical nuclear power generation etc. could be accomplished efficiently only by using superconducting proton/ions particle accelerators. As a result, many newly constructed or under construction high energy accelerator facilities around the world such as Proton Improvement Plan-II (PIP-II) at Fermilab, FRIB at MSU, ESS at Sweden are utilizing SRF accelerators. This talk presents the design philosophy of SRF proton/ion particle accelerator, with a focus on optimizing its performance while addressing physics and operational challenges including the fault scenarios, machine availability and preservation of crucial beam parameters.

Saini, Arun↗
DOE OSTI2024

Two-moment Neutrino Flavor Transformation with Applications to the Fast Flavor Instability in Neutron Star Mergers

Abstract Multi-messenger astrophysics has produced a wealth of data with much more to come in the future. This enormous data set will reveal new insights into the physics of core-collapse supernovae, neutron star mergers, and many other objects where it is actually possible, if not probable, that new physics is in operation. To tease out different possibilities, we will need to analyze signals from photons, neutrinos, gravitational waves, and chemical elements. This task is made all the more difficult when it is necessary to evolve the neutrino component of the radiation field and associated quantum-mechanical property of flavor in order to model the astrophysical system of interest—a numerical challenge that has not been addressed to this day. In this work, we take a step in this direction by adopting the technique of angular-integrated moments with a truncated tower of dynamical equations and a closure, convolving the flavor-transformation with spatial transport to evolve the neutrino radiation quantum field. We show that moments capture the dynamical features of fast flavor instabilities in a variety of systems, although our technique is by no means a universal blueprint for solving fast flavor transformation. To evaluate the effectiveness of our moment results, we compare to a more precise particle-in-cell method. Based on our results, we propose areas for improvement and application to complementary techniques in the future.

Astronomy & Astrophysics↗
DOE OSTI2022

Evaluating seasonal and regional distribution of snowfall in regional climate model simulations in the Arctic

In this study, we investigate how the regional climate model HIRHAM5 reproduces the spatial and temporal distribution of Arctic snowfall when compared to CloudSat satellite observations during the examined period of 2007–2010. For this purpose, both approaches, i.e., the assessments of the surface snowfall rate (observation-to-model) and the radar reflectivity factor profiles (model-to-observation), are carried out considering spatial and temporal sampling differences. The HIRHAM5 model, which is constrained in its synoptic representation by nudging to ERA-Interim, represents the snowfall in the Arctic region well in comparison to CloudSat products. The spatial distribution of the snowfall patterns is similar in both identifying the southeastern coast of Greenland and the North Atlantic corridor as regions gaining more than twice as much snowfall as the Arctic average, defined here for latitudes between 66 and 81°N. Excellent agreement (difference less than 1%) in the Arctic-averaged annual snowfall rate between HIRHAM5 and CloudSat is found, whereas ERA-Interim reanalysis shows an underestimation of 45% and significant deficits in the representation of the snowfall rate distribution. From the spatial analysis, it can be seen that the largest differences in the mean annual snowfall rates are an overestimation near the coastlines of Greenland and other regions with large orographic variations as well as an underestimation in the northern North Atlantic Ocean. To a large extent, the differences can be explained by clutter contamination, blind zone or higher resolution of CloudSat measurements, but clearly HIRHAM5 overestimates the orographic-driven precipitation. The underestimation of HIRHAM5 within the North Atlantic corridor south of Svalbard is likely connected to a poor description of the marine cold air outbreaks which could be identified by separating snowfall into different circulation weather type regimes. By simulating the radar reflectivity factor profiles from HIRHAM5 utilizing the Passive and Active Microwave TRAnsfer (PAMTRA) forward-modeling operator, the contribution of individual hydrometeor types can be assessed. Looking at a latitude band at 72–73°N, snow can be identified as the hydrometeor type dominating radar reflectivity factor values across all seasons. The largest differences between the observed and simulated reflectivity factor values are related to the contribution of cloud ice particles, which is underestimated in the model, most likely due to the small sizes of the particles. The model-to-observation approach offers a promising diagnostic when improving cloud schemes, as illustrated by comparison of different schemes available for HIRHAM5.

54 ENVIRONMENTAL SCIENCES↗