Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “field test”

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 919 records · Page 51

Rear-Side Irradiance Simulation of Field PV Modules: Preprint

Assessing the durability of photovoltaic (PV) module backsheet is critical to increasing module lifetime. Historically, laboratory-based accelerated testing has insufficient in predicting large-scale failures of commercial polymeric materials. Additionally, there is growing concern that standard condition tests do not reflect non-uniformities in field exposure, and that certain modules experience more severe degradation due to their location. Anisotropy in field exposures is installation dependent and reflects different levels of exposure to irradiance due to mounting geometry, ground surface albedo, and climatic zone. "Bificial Radiance" [1], simulates the amount of reflected irradiance on the Polyethylene Naphthalate (PEN) backsheet of a PV array located in Maryland. In our work, site specific weather data are gathered for the entire length of exposure for the modules. The total full-spectrum dose for the incident irradiance on the backsheet was then determined. The simulation results are integrated with historical field survey data to better understand realworld outdoor degradation.

backsheet↗

Characterization of Losses in Superconducting Radio-Frequency Cavities by Combined Temperature and Magnetic Field Mapping

Superconducting radio-frequency (SRF) cavities are one of the fundamental building blocks of modern particle accelerators. To achieve the highest quality factors (10 10 -10 11 ), SRF cavities are operated at liquid helium temperatures. Magnetic flux trapped on the surface of SRF cavities during cool-down below the critical temperature is one of the leading sources of residual RF losses. Instruments capable of detecting the distribution of trapped flux on the cavity surface are in high demand in order to better understand its relation to the cavity material, surface treatments and environmental conditions. We have designed, developed, and commissioned two novel diagnostic tools to measure the distribution of trapped flux at the surface of SRF cavities. One is a magnetic field scanning system (MFSS) which uses cryogenic Hall probes and anisotropic magnetoresistance sensors that fit the contour of a 1.3 GHz cavity. The second setup is a stationary, combined magnetic and temperature mapping system which uses AMR sensors and carbon resistor temperature sensors, covering the surface of a 3 GHz SRF cavity. The MFSS system revealed a non-uniform distribution of trapped flux on the cavities' surface, dependent on the magnitude of the applied magnetic field during field-cooling below the critical temperature. The MFSS shows that magnetic field scanning as a function of the RF field indicates redistribution of trapped flux at some locations. About ~ 33% of hot-spots observed by temperature mapping during high power RF tests overlapped with the high B-field spots. Almost all high B-field spots observed after field-cool were found to be overlapped on grain boundaries. A clear correlation between the hot- spot formed after quench and local trapped flux was found, providing insight on the RF dissipation of trap vortices. The combined B&T map system showed that the T-map system is capable of detecting hot-spots and quench location on the surface of the 3 GHz SRF cavity. Different distribution of trapped flux was measured after different cool-down and residual B-field, but, no variation in magnetic field distribution was observed during quench, possibly due to magnetic sensors being far from the quench location.

Parajuli, Ishwari↗

Automated Quantification of Wind Turbine Blade Leading Edge Erosion from Field Images

Wind turbine blade leading edge erosion is a major source of power production loss and early detection benefits optimization of repair strategies. Two machine learning (ML) models are developed and evaluated for automated quantification of the areal extent, morphology and nature (deep, shallow) of damage from field images. The supervised ML model employs convolutional neural networks (CNN) and learns features (specific types of damage) present in an annotated set of training images. The unsupervised approach aggregates pixel intensity thresholding with calculation of pixel-by-pixel shadow ratio (PTS) to independently identify features within images. The models are developed and tested using a dataset of 140 field images. The images sample across a range of blade orientation, aspect ratio, lighting and resolution. Each model (CNN v PTS) is applied to quantify the percent area of the visible blade that is damaged and classifies the damage into deep or shallow using only the images as input. Both models successfully identify approximately 65% of total damage area in the independent images, and both perform better at quantifying deep damage. The CNN is more successful at identifying shallow damage and exhibits better performance when applied to the images after they are preprocessed to a common blade orientation.

Aird, Jeanie A.↗

Characterizing falling particle curtain receivers at commercially relevant scales: Research Performance Progress Report (RPPR-1)

Sandia will construct a cold flow receiver test platform in order to characterize falling particle curtain receivers at commercially relevant scales. While Sandia has extensive experience in R&D of falling particle curtain receivers, most have been at pilot scale and smaller - on the order of 1 MWth with characteristic dimensions of nominally 1-2 m to adequately collect solar energy from the heliostat field at the National Solar Thermal Test Facility (NSTTF). However, scaling up receivers to commercially relevant scales (25 MWth and above) will require a thorough understanding of particle curtain dynamics at larger scales, especially longer drop heights, for design certainty. The goal of this project will be to construct a cold falling particle curtain test rig capable of simulating particle characteristics that are expected in a commercial scale CSP plant, namely the drop height, curtain thickness, and particle mass flow rate (normalized by length of curtain). This will enable data collection on curtain opacity and spread, both of which are correlated to receiver efficiency and reliable construction, for commercially relevant scales. It will also permit validation of numerical models that will enable detailed receiver characterization and design past currently validated scales.

14 SOLAR ENERGY↗

Analyzing the Validity of Brazilian Testing Using Digital Image Correlation and Numerical Simulation Techniques

Characterizing the mechanical behavior of rocks plays a crucial role to optimize the fracturing process in unconventional reservoirs. However, due to the intrinsic anisotropy and heterogeneity in unconventional resources, fracture process prediction remains the most significant challenge for sustainable and economic hydrocarbon production. During the deformation tracking under compression, deploying conventional methods (strain gauge, extensometer, etc.) is insufficient to measure the deformation since the physical attachment of the device is restricted to the size of the sample, monitoring limited point-wise deformation, producing difficulties in data retrieval, and a tendency to lose track in failure points, etc. Where conventional methods are limited, the application of digital image correlation (DIC) provides detailed and additional information of strain evolution and fracture patterns under loading. DIC is an image-based optical method that records an object with a camera and monitors the random contrast speckle pattern painted on the facing surface of the specimen. To overcome the existing limitations, this paper presents numerical modeling of Brazilian disc tests under quasi-static conditions to understand the full-field deformation behaviors and finally, it is validated by DIC. As the direct tensile test has limitations in sample preparation and test execution, the Brazilian testing principle is commonly used to evaluate indirectly the tensile strength of rocks. The two-dimensional numerical model was built to predict the stress distribution and full-field deformation on Brazilian disc under compression based on the assumptions of a homogenous, isotropic and linear elastic material. The uniaxial compression test was conducted using the DIC technique to determine the elastic properties of Spider Berea sandstone, which were used as inputs for the simulation model. The model was verified by the analytical solution and compared with the digital image correlation. The numerical simulation results showed that the solutions matched reasonably with the analytical solutions where the maximum deviation of stress distribution was obtained as 14.59%. The strain evolution (normal and shear strains) and displacements along the central horizontal and vertical planes were investigated in three distinguishable percentages of peak loads (20%, 40%, and 90%) to understand the deformation behaviors in rock. The simulation results demonstrated that the strain evolution contours consistently matched with DIC generated contours with a reasonable agreement. The changes in displacement along the central horizontal and vertical planes showed that numerical simulation and DIC generated experimental results were repeatable and matched closely. In terms of validation, Brazilian testing to measure the indirect tensile strength of rocks is still an issue of debate. The numerical model of fracture propagation supported by digital image correlation from this study can be used to explain the fracturing process in the homogeneous material and can be extended to non-homogeneous cases by incorporating heterogeneity, which is essential for rock mechanics field applications.

58 GEOSCIENCES↗

Characterizing falling particle curtain receivers at commercially relevant scales: Research Performance Progress Report (RPPR-1)

Sandia will construct a cold flow receiver test platform in order to characterize falling particle curtain receivers at commercially relevant scales. While Sandia has extensive experience in R&D of falling particle curtain receivers, most have been at pilot scale and smaller - on the order of 1 MW th with characteristic dimensions of nominally 1-2 m to adequately collect solar energy from the heliostat field at the National Solar Thermal Test Facility (NSTTF). However, scaling up receivers to commercially relevant scales (25 MW th and above) will require a thorough understanding of particle curtain dynamics at larger scales, especially longer drop heights, for design certainty. The goal of this project will be to construct a cold falling particle curtain test rig capable of simulating particle characteristics that are expected in a commercial scale CSP plant, namely the drop height, curtain thickness, and particle mass flow rate (normalized by length of curtain). This will enable data collection on curtain opacity and spread, both of which are correlated to receiver efficiency and reliable construction, for commercially relevant scales. It will also permit validation of numerical models that will enable detailed receiver characterization and design past currently validated scales.

14 SOLAR ENERGY↗

HERMES 2019 Upgrade Characterization - Quick Look Report (UUR)

Characterizing the electromagnetic environment created by the recently upgraded HERMES III pulser is the primary objective addressed by this report. The pulser upgrade design was intended not only to provide higher electric field and magnetic field illumination amplitudes but also to provide a more collimated beam, focused along the center-line of the courtyard. For high-voltage pulsers to be useful as qualification systems, they must have the requisite peak values and uniform, predictable peak field spatial patterns. To assess the peak electric field values and uniformity of the spatial pattern of the electric fields for the upgraded HERMES system, a 1 6 full-power shot illumination of 9 electric and 9 magnetic field sensors in the courtyard of the HERMES III facility was conducted. The data acquired from this test showed an increase an increase in electric field amplitude, but still an asymmetry in the spatial distributions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Using isotope shift for testing nuclear theory: The case of nobelium isotopes

We calculate field isotope shifts for nobelium atoms using nuclear charge distributions which come from different nuclear models.We demonstrate that comparing calculated isotope shifts with experiment can serve as a testing ground for nuclear theories. It also provides a way of extracting parameters of nuclear charge distribution beyond nuclear root mean square (rms) radius, e.g., parameter of quadrupole deformation β. The measurements of at least two atomic transitions is needed to disentangle the contributions of the changes in deformation and nuclear rms radius into field isotopic shift. We argue that a previous interpretation of the isotope measurements in terms of δ < r 2 > between 252,254 No isotopes should be amended when nuclear deformation is taken into account. As a result, we calculate isotope shifts for other known isotopes and for hypothetically metastable isotope 286 No for which the predictions of nuclear models differ substantially.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Magnetic Field Mapping of a 2.5 T Fixed-Field HTS Gantry Magnet for Proton Therapy

We present results from testing a high-temperature superconducting (HTS) magnet prototype for proton therapy. This magnet is specifically designed for a novel rotating gantry capable of delivering the entire proton beam energy range (70225 MeV) while maintaining a fixed magnetic field in the superconducting magnets. The gantry layout simplifies the magnet design, enabling the use of straight, flat racetrack Bi-2223 (DI-BSCCO) coil technology and operation at higher temperatures (1015 K). The magnet has a non-linear field distribution for bending and focusing the proton beams. To validate this feature, we developed a system for measuring the magnetic field distribution in the magnet aperture. We present the design of this hall probe array and experimental results from two different magnet tests at 4.2 K in a liquid helium bath. These results are compared with the simulated field distribution and discussed in the context of the required field quality for the application.

Mosat, M↗

A windowed mean trajectory approximation for condensed phase dynamics

We propose a trajectory-based quasi-classical method for approximating dynamics in condensed phase systems. Building upon the previously developed optimized mean trajectory approximation that has been used to compute linear and nonlinear spectra, we borrow some ideas from filtering trajectory methods to obtain a novel semiclassical method for the dynamical propagation of density matrices. This new approximation is tested rigorously against standard multistate electronic models, spin-boson models, and models of the Fenna–Matthews–Olson complex. For dissipative systems, the current method is significantly better or as good as many other semiclassical methods available, especially at low temperatures and for off-diagonal density matrix elements, whereas for scattering models, the current method bears similar limitations as mean-field propagation schemes. All results are tested against the numerically exact hierarchical equations of motion method. In conclusion, the new method shows excellent agreement across various parameter regimes with numerically exact results, highlighting the robustness and accuracy of our approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental evaluation of localized creep deformation in grade 91 steel weldments

Spatially resolved measurement of localized creep deformation in heterogeneous creep resistant steel weldments is crucial but challenging for lifetime assessments of critical steam components in power plants. Here, experimental approaches were established to quantitatively evaluate commonly observed localized creep deformation in multi-pass Grade 91 steel weldments. An in-situ digital image correlation (DIC) system was utilized with a creep testing frame to monitor and measure both full-field strain and localized strain accumulation across the weldments during long-term creep testing at elevated temperatures. The in situ DIC method measured not only the creep deformation behavior of the weld metal, heat affected zone (HAZ), and base metal, but also creep strain evolution for each sub-region within the HAZ itself, including coarse-grained HAZ, fine-grained HAZ (FGHAZ), and intercritical HAZ (ICHAZ). The DIC results revealed that local creep strain in the ICHAZ reached up to 90% strain before final rupture, whereas nominal creep strain measured by the standard extensometer of the tested cross-weld specimen was below 10%, indicative of Type IV cracking of the Grade 91 weld. Microstructural analyses revealed that the faster creep degradation/deformation in the HAZ was caused mainly by accelerated matrix grain recrystallization/growth and a reduced pinning effect from the segregated and coarsened precipitates in the FGHAZ and ICHAZ. The ultimate creep rupture occurred in the ICHAZ owing to its lowest creep resistance induced by the largest recrystallized grain size, the lowest fraction of coincidence site lattice, and the lowest local strain energies/dislocation densities.

36 MATERIALS SCIENCE↗

Vertically Resolved Convective–Stratiform Echo-Type Identification and Convectivity Retrieval for Vertically Pointing Radars

Using data from the airborne HIAPER Cloud Radar (HCR), a partitioning algorithm (ECCO-V) that provides vertically resolved convectivity and convective versus stratiform radar-echo classification is developed for vertically pointing radars. The algorithm is based on the calculation of reflectivity and radial velocity texture fields that measure the horizontal homogeneity of cloud and precipitation features. The texture fields are translated into convectivity, a numerical measure of the convective or stratiform nature of each data point. The convective–stratiform classification is obtained by thresholding the convectivity field. Subcategories of low, mid-, and high stratiform, shallow, mid-, deep, and elevated convective, and mixed echoes are introduced, which are based on the melting-layer and divergence-level altitudes. As the algorithm provides vertically resolved classifications, it is capable of identifying different types of vertically layered echoes, and convective features that are embedded in stratiform cloud layers. Its robustness was tested on data from four HCR field campaigns that took place in different meteorological and climatological regimes. The algorithm was adapted for use in spaceborne and ground-based radars, proving its versatility, as it is adaptable not only to different radar types and wavelengths, but also different research applications.

54 ENVIRONMENTAL SCIENCES↗

Transport Properties of Magnetized High-Energy Density Plasma (Final Project Report)

This report summarizes results of project DE-SC0016159 "Transport Properties of Magnetized High-Energy-Density Plasma,'' which ran from 7/15/2016 - 7/14/2021. The primary objectives of the work, as stated in the original proposal, were to develop a theory that describes transport coefficients in magnetized high energy density plasmas and to test the theory with molecular dynamics (MD) simulations. The science challenge is that strong ion coupling, strong magnetization of electrons, and partial degeneracy of electrons, cause the system to be in a regime that is not well described by current theory. The particular processes that were to be investigated include: ion stopping power, diffusion, electron-ion temperature relaxation, thermal conduction and viscosity. All of the primary research objectives were accomplished during the course of this work. In addition, some unexpected results were discovered along the way that led to new and productive research directions. This work resulted in 17 publications in well-respected peer-reviewed journals (primarily Physics of Plasmas and Physical Review E), and 2 more are being prepared for publication. A few of these were selected as Editor's Picks and one was published as a Rapid Communication. Highlights of the research results include: Transport phase space: The first identification of the parameter space that defines fundamental transport regimes in terms of the Coulomb coupling and magnetization parameters. Mean force kinetic theory: Systematic derivation of a kinetic theory for strongly coupled plasmas based on a new expansion parameter of the BBGKY hierarchy. This provided the derivation of a theory that we had previously posed phenomenologically, and also revealed a new term that captures the equation of state properties at all coupling strengths. Transverse friction force: Discovery that the friction force on a test charge in a strongly magnetized plasma includes a component that is perpendicular to the motion of the test charge in the plane defined by the velocity and magnetic field vectors. Gyrofriction force: Discovery that the friction force on a test charge in a plasma that is both strongly magnetized and strongly coupled includes a component of the force in the direction of the Lorentz force. Kinetic theory for warm dense matter: Extension of the mean force kinetic theory to include partial electron degeneracy, so that it applies to warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantifying and Optimizing the Energy Benefits of Mass Timber Construction

The International Mass Timber Alliance (IMTA) is a global organization of industry leaders, engineers, scientists, and associations dedicated to advancing mass timber construction. Its mission is to generate and disseminate scientific data supporting the development of standardized construction and energy efficient practices that promote the adoption of mass timber worldwide. IMTA collaborated with Oak Ridge National Laboratory (ORNL) to leverage ORNL’s expertise in building envelope modeling and testing to evaluate how mass timber construction can reduce peak heating and cooling demand, lower overall energy use, and improve resilience during power outages. A previous study of 80 mass timber buildings in Finland found measured energy use up to 50% lower than predicted by simulation. This project aimed to validate and extend those findings for U.S. buildings through analytical modeling, laboratory testing, and full-scale building evaluations. The research focused on the thermal performance of low-embodied-energy wall assemblies, such as cross-laminated timber (CLT) panels and log walls, with particular attention to the effects of thermal inertia on indoor comfort and energy performance. While mass timber’s structural and fire-resistance properties are well documented, its whole-building thermal behavior has received limited attention. Field data, simulation results, and resilience testing from this study will inform future modeling practices, design guidelines, and construction practices by quantifying the unique thermal and demand-flexibility benefits of mass timber construction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Proxy-Based Bayesian Inversion Of Poroelastic Simulations To Interpret Strain Tensor Data Measured During Well Testing

The long runtimes of 3D poroelastic numerical simulators makes it impractical to interpret deformation datasets using many inversion schemes. Recent advances in instrumentation have made it possible to measure the strain tensor during well testing, but the lack of robust inversion methods is limiting the ability to interpret these data. We have developed an inversion workflow that reduces the number of computations required to complete a Bayesian inversion using DREAMzs. The workflow trains a KNN model using output from the poroelastic simulator, and then uses the KNN model as a proxy for the simulator during inversion. The workflow also includes a strategy for ensuring the results from the proxy model converge to the results from the simulator, ensuring the accuracy of the final results. An idealized example configured to represent a well test in a deep aquifer is used to verify that the workflow correctly identifies parameters and characterizes noise. Field data measured using strainmeters during an injection test at an oil reservoir in Oklahoma are used to evaluate performance with a real dataset. The workflow identified 265 history matching solutions out of 1240 total simulation runs (21% acceptance ratio), and the results are used to characterize posterior parameter distribution and evaluate the prediction uncertainty. This approach makes it feasible to invert strain data measured during well testing and this has the potential to improve the characterization of aquifers and reservoirs.

Roudini, Soheil↗

Efficient near-field ptychography reconstruction using the Hessian operator

X-ray ptychography is a powerful and robust coherent imaging method providing access to the complex object and probe (illumination). Ptychography reconstruction is typically performed using first-order methods due to their computational efficiency. Higher-order methods, while potentially more accurate, are often prohibitively expensive in terms of computation. In this study, we present a mathematical framework for reconstruction using second-order information derived from an efficient computation of the bilinear Hessian and Hessian operator. The formulation is provided for Gaussian-based models, enabling the simultaneous reconstruction of the object, probe, and object positions. Synthetic data tests, along with experimental near-field ptychography data processing, demonstrate a ten-fold reduction in computation time compared to first-order methods. The derived formulas for computing the Hessians, along with the strategies for incorporating them into optimization schemes, are well-structured and easily adaptable to various ptychography problem formulations.

Carlsson, Marcus [Lund Univ. (Sweden)] (ORCID:0000↗

Novel Concepts for High Gradient Acceleration (Final Technical Report)

We have conducted an intensive, pioneering program to demonstrate novel concepts for achieving high gradient acceleration at frequencies from the conventional microwave bands up to the millimeter wave /THz bands. High gradient accelerators hold the promise of smaller and less costly accelerators for applications ranging from the largest scale accelerators used for discovery science down to the smallest accelerators used for industrial, homeland security and medical applications. The research consisted of two major research thrusts: 1.) Structure-based wakefield accelerator (SWFA) research in collaboration with the Argonne Wakefield Accelerator research group and 2.) Millimeter Wave / THz high gradient acceleration in collaboration with SLAC. The specific goals of the research program were: Design novel metallic metamaterial structures that increase the beam-wave coupling for the accelerator mode and reduce the effect of high order modes; test novel metamaterial structures to achieve higher output power, > 1 GW at X-Band (11.7 GHz), in test at the Argonne Wakefield Accelerator (AWA); determine the experimental breakdown threshold for nanosecond-scale pulses at X-Band in testing at the AWA; test a 110 GHz accelerator structure with a field emission electron gun, built at SLAC, using pulses from a 1 MW, 110 GHz gyrotron; design, build and test a 110 GHz quasi-optical, resonant-ring pulse compressor to compress microsecond pulses from the 1 MW gyrotron into > 20 MW, 5 ns output pulses for accelerator structure testing. The proposed research program built on our successes in our research program including: Generation of 510 MW, 2.1 ns (FWHM) pulses at 11.7 GHz from a metallic metamaterial structure in test at the Argonne Wakefield Accelerator using a train of eight 65 MeV electron bunches spaced at 1.3 GHz with a total charge of 280 nC. The metamaterial structure consisted of 100 copper unit cells each consisting of a “wagon-wheel” plate and a spacer plate with a total structure length of 0.2 m. The 510 MW pulse generated an on-axis wakefield of 130 MV/m that could be used to accelerate a trailing witness bunch. Demonstration of coupling of an unprecedented rf power level of 575kW into a 110 GHz accelerator structure using a quasi-optical setup. The standing structure consisted of a central copper cavity located between two matching cavities fed by a TM01 mode. The 6 ns input pulses were sliced from 3 microsecond pulses from the gyrotron using a laser-driven silicon switch. We obtained an unprecedented high gradient up to 230MV/m corresponding to a peak surface electric field of more than 520 MV/m.

43 PARTICLE ACCELERATORS↗

Developing and Testing a Novel Stochastic Ice Microphysics Parameterization for Cloud and Climate Models Using ARM Field Campaign Data (Final Progress Report)

The major goals of this project were: 1) to use field campaign measurements from DOE’s Atmospheric Radiation Measurement (ARM) program to characterize variability of important parameters describing properties of ice particles in the atmosphere; 2) based on this observational analysis, to develop a parameterization scheme for weather and climate models that stochastically varies these parameters, and implement the new scheme into a weather model called the Weather Research and Forecasting model (WRF); 3) to use WRF coupled with the new stochastic scheme to simulate ARM field campaign thunderstorm cases and analyze how accounting for this parameter variability affects the model simulations. This work was performed jointly between the National Center for Atmospheric Research, University of Oklahoma, and University of Utah. To accomplish these goals, we extended an approach previously developed to characterize the variability in the size distribution of ice particles to parameters that are explicitly represented in models (i.e., relationships between ice particle mass and size, and between particle fall velocity and size). Our project was, to our knowledge, the first to apply observationally-constrained estimates of this parameter variability describing mass-size and fall velocity-size in a modeling framework. Our results showed efficacy of the approach, evaluated using ARM observations. Similarly, to our knowledge, work in this project was the first to propose and evaluate in detail a stochastic approach for unresolved turbulent mixing in high-resolution model simulations against detailed, benchmark large eddy simulations and ARM observations. Results showed some promising behavior, particularly with increased mixing and dilution of air in thunderstorm cores with surrounding environmental air, bringing the stochastic simulations closer to the benchmark large eddy simulations; however, results were somewhat degraded using stochastic mixing compared to observations from the AMIE/DYNAMO field campaign. This project also further refined and applied a modeling methodology called “piggybacking” that can robustly separate dynamical and thermodynamic impacts of model changes, and comparison studies of different models based on cases developed from ARM observations. Finally, this project directly supported three graduate students who completed their PhDs as well as a postdoctoral research fellow.

54 ENVIRONMENTAL SCIENCES↗