Search NASA⌕ Search

SEARCH · Search NASA

Results for “kinetic distributions”

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 199 records · Page 11

A multi-dimensional kinetic-based upwind solver for the Euler equations

A multidimensional kinetic fluctuation-splitting scheme has been developed for the Euler equations. The scheme is based on an N-scheme discretization of the Boltzmann equation at the kinetic level for triangulated Cartesian meshes with a diagonal-adaptive strategy. The resulting Euler scheme is a cell-vertex fluctuation-splitting scheme where fluctuations in the conserved-variable vector Q are obtained as moments of the fluctuation in the Maxwellian velocity distribution function at the kinetic level. Encouraging preliminary results have been obtained for perfect gases on Cartesian meshes with first-order spatial accuracy. The present approach represents an improvement to the well-established dimensionally-split upwind schemes.

Eppard, W. M.↗

Explosive Soot Challenge (Final Report)

This project assembled a broad ensemble of modeling and experimentation tools to study the morphological and optical properties of detonation soots in explosive fireballs. A gram-scale hemispherical high explosive was studied in a low-pressure controlled environment using in-situ experimentation with diffusely illuminated visible absorption spectroscopy, particle sizing through light scattering techniques, and post-test collections with subsequent morphological analysis. Hydrocode modeling was performed to replicate the detonation flow observations, and subsequent aerosol kinetics models provided particle size distributions and extinction coefficients from the hydrocode results. Experimentally observed soot morphologies agreed with expectation from the literature - a bimodal distribution was found, brought upon by the particles growing to a size where their inertia and fluid wakes are non-negligible. The aerosol kinetics model did not replicate the observed bimodal size distribution for lack of a coagulation kernel to represent the behavior. To recover particulate optical properties, a spectrally resolved absorption spectroscopy method termed Spectral diffuse back-illuminated extinction imaging (SBI-EI) was developed and implemented on two explosive types. Inverting the absorption spectra using a Kramers-Kronig consistent method yielded the complex index of refraction for the soots produced by the explosives. This method resulted in an unrealistic index of refraction for one of the two explosives, and this is suspected to be due to the model neglecting scattering brought upon by the large particle sizes observed. In addition to the core work, three additional studies were performed in parallel. These investigated the impact of scattering on diffuse absorption spectroscopy, studied how soots oxidate and sublimate in a well-controlled shock tube, and laid the theoretical groundwork for a new collision kernel to replicate the bimodal size distribution from the observations. Summaries of these efforts are included at the end of this report.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Fission with Exotic Nuclei (Abbreviated Report)

Nuclear fission is a key mechanism involved in the synthesis of heavy elements in the Cosmos and is the primary explanation for the stability of superheavy elements. Nevertheless, our knowledge of fission remains extremely fragmented. Most experiments have been conducted only on a tiny number of stable actinide nuclei and are often incomplete, leading to gaps in our basic understanding of the process. For many radioactive isotopes, basic fission data such as the charge or mass distribution of the fragments is unknown. These gaps cannot always be filled by simulation alone. Common fission models contain too many free parameters and lack predictive power. In contrast, the fundamental theory of fission under development at LLNL is much more predictive, but its current computational cost is too high to be used extensively for data evaluations. A unique window of opportunity to resolve these limitations has recently opened: the U.S. nuclear science community is ramping up major experimental programs at the Facility for Rare Isotope Beams (FRIB, the DOE flagship facility in low-energy nuclear science), and techniques from machine learning have shown great potential to simplify the use of a fundamental, quantum-mechanical theory of fission. This project has two components. On the experimental side, we acquired and deployed at the HIGS facility a new dual Frisch-Grid ionization chamber to measure correlated fragment-mass, kinetic energy, and angular distributions of fission fragments from induced fission. This new device was used to perform measurements of charge, mass and total kinetic energy of fission fragments in the photofission of 238 U and eight gamma-ray beam energies between 6.2 and 13 MeV, which allowed extracting high-precision independent yields for this reaction. The device was also used to perform measurements of the same quantities in the neutron-induced fission of 234 U with monoenergetic beams of energy between 5 and 8 MeV. In parallel, we collaborated with a team at Commissariat à l’énergie atomique et aux énergies alternatives (CEA) to perform a series of measurements of fission yields in inverse kinematics for the two isotopes of 236 U and 240 Pu. The experiment took place at the Grand Accélérateur National d’Ions Lourds in France in June 2023. The deployment of the VAMOS spectrometer with a new array called PISTA allowed determining the excitation energy of the fissioning system within 1 Mega-electronvolts. The second component of the project involved using deep neural networks to build fast and reliable emulators of our current fission models. In an invited paper published in Frontier in Physics, we showed that autoencoders could successfully compress nuclear wavefunctions in nuclear density functional theory. We achieved a dimensionality reduction of the order of two orders of magnitude while keeping the error in the total energy to less than 0.01%. In a second paper submitted to Physical Review Letters in June 2023 with our collaborators at CEA, we showed that variational autoencoders can learn the collective degrees of freedom driving the fission process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multiscale Modeling of Nanoparticle Precipitation in Oxide Dispersion-Strengthened Steels Produced by Laser Powder Bed Fusion

Laser Powder Bed Fusion (LPBF) enables the efficient production of near-net-shape oxide dispersion-strengthened (ODS) alloys, which possess superior mechanical properties due to oxide nanoparticles (e.g., yttrium oxide, Y-O, and yttrium-titanium oxide, Y-Ti-O) embedded in the alloy matrix. To better understand the precipitation mechanisms of the oxide nanoparticles and predict their size distribution under LPBF conditions, we developed an innovative physics-based multiscale modeling strategy that incorporates multiple computational approaches. These include a finite volume method model (Flow3D) to analyze the temperature field and cooling rate of the melt pool during the LPBF process, a density functional theory model to calculate the binding energy of Y-O particles and the temperature-dependent diffusivities of Y and O in molten 316L stainless steel (SS), and a cluster dynamics model to evaluate the kinetic evolution and size distribution of Y-O nanoparticles in as-fabricated 316L SS ODS alloys. The model-predicted particle sizes exhibit good agreement with experimental measurements across various LPBF process parameters, i.e., laser power (110–220 W) and scanning speed (150–900 mm/s), demonstrating the reliability and predictive power of the modeling approach. The multiscale approach can be used to guide the future design of experimental process parameters to control oxide nanoparticle characteristics in LPBF-manufactured ODS alloys. Additionally, our approach introduces a novel strategy for understanding and modeling the thermodynamics and kinetics of precipitation in high-temperature systems, particularly molten alloys.

Wang, Zhengming (ORCID:0000000241627112)↗

Stratosphere chemistry in a 2-D model with residual circulation

The objective of this research was to examine the effects of chemical perturbations on the stratosphere using models which can incorporate fully interactive radiative, chemical, and dynamical responses, in the context of a zonally averaged model. Model runs for the unperturbed, chlorine-perturbed and simultaneously chlorine-and CO2-perturbed cases were completed using the JPL-87 chemical kinetics data. The base case was analyzed and submitted for publication. The perturbed cases show substantial sensitivity of the predicted column ozone depletion to the perturbations affecting lower stratosphere temperature, but less to far dynamical perturbations. The column ozone distribution changed substantially when the kinetics data was changed. This implies a greater-than-expected uncertainty in predicted latitude distributions of ozone depletion, due to uncertainty about the accuracy and completeness of the chemical kinetics data set.

Guthrie, Paul D.↗

Shock propagation and the generation of magnetohydrodynamic wave fields in inhomogeneous molecular clouds

We develop a simple one-dimensional model for the interaction of a steady, thin, planar shock wave with a nonrigid cloud which may be in motion relative to the surrounding medium, and we apply the model to shocks impinging on, and propagating through, molecular clouds. Both 'adiabatic' (gamma = 5/3) and radiative (gamma = 1) shocks are considered and we allow for the presence of a uniform magnetic field directed either parallel or perpendicular to the shock normal. The former field orientation is equivalent to the hydrodynamic case, and the latter involves only fast MHD shocks. We focus on the manner in which such shocks can generate internal kinetic motions in the cloud on a range of size and density scales through the direct acceleration of cores and clumps by shocks transmitted into them and through the generation of an MHD wavefield via the reflection of the incident shock at clump boundaries. We find that stronger incident Mach numbers and smaller density contrasts lead to more efficient cloud acceleration, as do isothermal intercloud shocks and small intercloud magnetic field strengths. The acceleration efficiency is insensitive to the adiabatic index and the magnetic field strength in the cloud itself. For typical parameter choices, the direct acceleration of clouds and clumps by strong shocks is found to be substantial and could at least in part account for their observed velocity dispersions. If the shocks are moderately weak, the final velocity of the cloud is linearly related to its initial velocity, with higher acceleration giving shallower slopes (i.e., final velocity distributions which are less sensitive to the initial distribution). Compared to the kinetic energy of the postshock cloud, the energy given to the wavefield at each encounter is small, and the heating of the interclump medium by the dissipation of this wavefield is found to be insufficient to balance the cooling rate in the cloud as a whole (although it may be important in particular regions), even if this medium is warm, unless it is also extremely tenuous (n approximately less than 0.1/cu cm). Nevertheless, the correction for the velocity imparted to the cloud leads to a substantial increase in the critical incident Mach numner for wave emission over that reported by Spitzer for the rigid case. The implications of our model for shock-induced star formation are discussed briefly.

Miesch, Mark S.↗

Stability of superthermal strahl electrons in the solar wind

We present a kinetic stability analysis of the solar wind electron distribution function consisting of the Maxwellian core and the magnetic-field aligned strahl, a superthermal electron beam propagating away from the sun. We use an electron strahl distribution function obtained as a solution of a weakly collisional drift-kinetic equation, representative of a strahl affected by Coulomb collisions but unadulterated by possible broadening from turbulence. This distribution function is essentially non-Maxwellian and varies with the heliospheric distance. The stability analysis is performed with the Vlasov–Maxwell linear solver leopard. We find that depending on the heliospheric distance, the core-strahl electron distribution becomes unstable with respect to sunward-propagating kinetic-Alfvén, magnetosonic, and whistler modes, in a broad range of propagation angles. The wavenumbers of the unstable modes are close to the ion inertial scales, and the radial distances at which the instabilities first appear are on the order of 1 au. However, we have not detected any instabilities driven by resonant wave interactions with the superthermal strahl electrons. Instead, the observed instabilities are triggered by a relative drift between the electron and ion cores necessary to maintain zero electric current in the solar wind frame (ion frame). Contrary to strahl distributions modelled by shifted Maxwellians, the electron strahl obtained as a solution of the kinetic equation is stable. Our results are consistent with the previous studies based on a more restricted solution for the electron strahl.

79 ASTRONOMY AND ASTROPHYSICS↗

Elucidating the Impact of Cis – Trans Organic Structure Directing Agent Isomer Ratios on the Aluminum Distribution Within SSZ-39

Despite their widespread use, the mechanisms governing the synthesis of zeolite catalysts are still poorly understood. A notable example of this problem is the uncertainty surrounding the influence of synthesis conditions on the placement of Al atoms in the zeolite framework which determines the active sites available for catalytic species. In this work, the role of the cis to trans isomer ratio of the OSDA N,N-dimethyl-3-5-dimethylpiperidinium on the energetics of 26 distinct Al pair distributions in SSZ-39 is examined both in the presence and absence of Na using density functional theory calculations. The initial orientation of the OSDA was found to have a significant impact on the final energies present, necessitating the screening of a large number of initial orientations with force field calculations and single point DFT calculations. Ground state energies were found to vary significantly with the ratio of cis to trans OSDAs with a Boltzmann distribution revealing the most likely Al pair distributions shift from sharing the same 8 membered rings to sharing the same double six membered rings to having no shared subunits as one increases the amount of cis OSDA present within the framework. The presence of Na was found to favor Al pair distributions where both Als occupied the same 6-membered ring. When an implicit solvent model was used to evaluate ground state energies the ideal Na sites shifted from 6-membered rings to empty SSZ-39 cages while OSDA positions and orientations remained largely the same. To provide insight on how kinetic factors may influence Al distributions, formation energies we calculated for connected double six membered rings. Further, these formation energies revealed a preference for Al pairs to occupy the same 4-membered ring which indicates kinetic and thermodynamic control may lead to different Al distributions in SSZ-39.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A kinetic model for the thermal nitridation of SiO2/Si

To explain the observed nitrogen distributions in thermally nitridated SiO2 films, a kinetic model is proposed in which the nitridation process is simulated, using the first-order chemical kinetics and Arrhenius dependence of the diffusion and reaction rates on temperature. The calculations show that initially, as the substrate reacts with diffusing nitrogen, a nitrogen-rich oxynitride forms at the SiO2-Si interface, while at nitridation temperatures above 1000 C, an oxygen-rich oxynitride subsequently forms at the interface, due to reaction of the substrate with an increasingly concentrated oxygen displaced by the slower nitridation of the SiO2. This sequence of events results in a nitrogen distribution in which the peak of the interfacial nitrogen concentration occurs away from the interface. The results are compared with the observed nitrogen distribution. The calculated results have correctly predicted the positions of the interfacial nitrogen peaks at the temperatures of 800, 1000, and 1150 C. To account for the observed width of the interfacial nitrogen distribution, it was found necessary to include in the simulations the effect of interfacial strain.

Vasquez, R. P.↗

Barkas effect in strongly magnetized plasmas

Strongly magnetized plasmas, which are characterized by the particle gyrofrequency exceeding the plasma frequency, exhibit novel transport properties. For example, recent work showed that the friction force on a test charge moving through a strongly magnetized plasma not only consists of the typical stopping power component but also includes components perpendicular to the test charge's velocity. However, these studies only considered test charges that have the same sign as the charge of the plasma particles. Here, we extend these calculations to the case of charges with opposite signs (such as an ion interacting with strongly magnetized electrons). This is done with both a novel generalized Boltzmann kinetic theory and molecular dynamics simulations. It is found that the friction force changes dramatically depending on the sign of the interacting charges. Likewise, the stopping power component for oppositely charged particles decreases in magnitude compared with like-charged particles, and the perpendicular components increase in magnitude. Moreover, the difference between the two cases increases as the gyrofrequency becomes larger compared with the plasma frequency. Finally, the electrical resistivity is calculated from the friction force, where it is found that strong magnetization in conjunction with oppositely charged interactions significantly decreases the parallel resistivity and increases the perpendicular resistivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Overview of recent SCALE activities for Non-LWR inventory and decay heat analysis

In 2019, the US Nuclear Regulatory Commission initiated a project for the development and assessments of non-light-water reactor (non-LWR) accident progression using the SCALE and MELCOR simulation tools. SCALE simulations are used to generate nuclide inventories, full-core power distributions, decay heat, and kinetics parameters to initialize MELCOR simulations of severe accident scenarios. Five non-LWR concepts were studied: high-temperature gas-cooled reactor (HTGR), heat pipe reactor (HPR), high-temperature fluoride salt-cooled reactor (FHR), molten salt-fueled reactor (MSR), and sodium-cooled fast reactor (SFR). This paper summarizes the SCALE results obtained in 2021 for the first three non-LWR concepts, compares characteristics and results to common LWRs, and provides the strategy for the analysis of the remaining two non-LWRs. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Operando Neutron Imaging of Lithium Flux and Gradient Cathode Design for Enhanced Kinetics in High‐Loading All‐Solid‐State Li─S Batteries

All-solid-state Li–sulfur batteries (ASSLSBs) are considered promising candidates for next-generation energy storage owing to their inherent safety, high energy density, and abundant sulfur resources. However, slow redox kinetics greatly limit sulfur utilization during solid-solid sulfur reactions, leading to significant challenges to achieve efficient performance in high-mass-loading ASSLSBs. Here, operando neutron image is employed to directly visualize, for the first time, that sluggish Li + transport kinetics and the uneven distribution of Li + during cathodic reactions are critical factors limiting sulfur conversion. To address this issue, gradient cathode architectures comprising three and five layers are designed, in which catholyte concentrations are strategically varied to optimize Li-ion flux and enhance ionic conductivity of the whole composite cathode electrode. Operando neutron imaging distinctly visualizes and confirms that three-layer gradient approach significantly enhances Li-ion mobility, resulting in more uniform redox reactions and greatly improved sulfur utilization compared to traditional non-gradient structures. Consequently, the three-layer gradient cathode achieves superior rate performance and reduced electrode polarization at high sulfur mass loadings of 4.5 and 6.0 mg cm −2 . Furthermore, the applicability and scalability of this design are demonstrated in a five-layer gradient cathode architecture, achieving an impressive discharge specific capacity increase from 656 mAh g −1 (three-layer gradient) to 1232 mAh g −1 at 1/20 C for ultra-high sulfur loading of 7.5 mg cm −2 . In conclusion, this innovative gradient cathode design offers substantial advancements in understanding and overcoming Li-ion transport limitations, paving the way toward practical, high-energy-density ASSLSBs.

25 ENERGY STORAGE↗

Turbulence modeling to aid tidal energy resource characterization in the Western Passage, Maine, USA

Numerical models combined with field measurements are regularly used to characterize tidal energy resources at potential energetic sites. However, most existing works only focus on the tidal hydrodynamic characteristics, and turbulence parameters are often not reported because of the lack of high-quality turbulence measurements and the limitations of numerical models in resolving turbulent eddies. In this study, we used FVCOM - a hydrostatic primitive equation (HPE) model - to characterize the tidal energy resource in the Western Passage, Maine, USA, by taking care of the essential macro-scale turbulence properties. We observed an excellent model performance using the Mellor-Yamada Level 2.5 Turbulence Model; estimating the spatial and vertical distribution of the turbulent kinetic energy and intensity added a new perspective to the site ranking for tidal energy converter (TEC) deployments. In addition, we also examined the role of channel geometry and bathymetry, such as headlands and underwater sills, in enhancing turbulent eddies around potential TEC siting locations. Ultimately, the detailed analysis of the turbulent flow characteristics has changed the site-ranking results and demonstrated that the regional-scale HPE models could be used for the relative understanding of more or less turbulent sites for a refined resource assessment.

16 TIDAL AND WAVE POWER↗

UV photodissociation dynamics of the acetone oxide Criegee intermediate: experiment and theory

Here, the photodissociation dynamics of the dimethyl-substituted acetone oxide Criegee intermediate [(CH 3 ) 2 COO] is characterized following electronic excitation to the bright 1 ππ* state, which leads to O ( 1 D) + acetone [(CH 3 ) 2 CO, S 0 ] products. The UV action spectrum of (CH 3 ) 2 COO recorded with O ( 1 D) detection under jet-cooled conditions is broad, unstructured, and essentially unchanged from the corresponding electronic absorption spectrum obtained using a UV-induced depletion method. This indicates that UV excitation of (CH 3 ) 2 COO leads predominantly to the O ( 1 D) product channel. A higher energy O ( 3 P) + (CH 3 ) 2 CO (T 1 ) product channel is not observed, although it is energetically accessible. In addition, complementary MS-CASPT2 trajectory surface-hopping (TSH) simulations indicate minimal population leading to the O ( 3 P) channel and non-unity overall probability for dissociation (within 100 fs). Velocity map imaging of the O ( 1 D) products is utilized to reveal the total kinetic energy release (TKER) distribution upon photodissociation of (CH 3 ) 2 COO at various UV excitation energies. Simulation of the TKER distributions is performed using a hybrid model that combines an impulsive model with a statistical component, the latter reflecting the longer-lived (>100 fs) trajectories identified in the TSH calculations. The impulsive model accounts for vibrational activation of (CH 3 ) 2 CO arising from geometrical changes between the Criegee intermediate and the carbonyl product, indicating the importance of CO stretch, CCO bend, and CC stretch along with activation of hindered rotation and rock of the methyl groups in the (CH 3 ) 2 CO product. Detailed comparison is also made with the TKER distribution arising from photodissociation dynamics of CH 2 OO upon UV excitation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonlinear dynamics of multiple Alfvén modes driven by trapped energetic ions in tokamaks: A triplet paradigm

Non-linear dynamics of multiple infernal Alfvén eigenmodes—a subset of global Alfvén eigenmodes in tokamak plasmas with extended low-shear central core [Marchenko et al., Phys. Plasmas 16, 092502 (2009)]—is studied. The analysis is carried out for a mode triplet with toroidal mode-numbers n=1, 2, 3. It was assumed that the n = 1 mode was linearly unstable due to precession resonance with trapped fast ions, whereas the other modes were linearly damped. The modes were coupled due to a non-linearity in a bounce-averaged drift kinetic equation for the distribution function of fast ions. Nonlinear equations for the mode amplitudes and phases are derived and solved numerically. It is found that the temporal evolution of the amplitudes and the phase (responsible for the frequency chirping) of the modes exhibit Hopf bifurcations to stable limit cycles. This can explain a synchronous cyclic destabilization of multiple modes in Alfvén avalanches (sudden growth of amplitudes of the mode cluster with different n and approximately equal frequency spacing) in NSTX and bursting modes in MAST—events, which resulted in enhanced loss of fast ions.

Marchenko, V. S. (ORCID:0000000261480533)↗

Tracing inter-Coulombic decay of molecular dimers

We have conducted an experimental study on the photo double ionization (PDI) of carbon dioxide dimers and oxygen dimers, while focusing on the dissociation dynamics upon single photon absorption. The results in terms of the kinetic energy and angular distributions of the charged particles show unambiguous experimental evidence of intermolecular Coulombic decay (ICD) in carbon dioxide dimers. In the oxygen dimer, the results show that ICD is accompanied by knock-off ionization mechanisms.

74 ATOMIC AND MOLECULAR PHYSICS↗

Accuracy limitations for composition analysis by XPS using relative peak intensities: LiF as an example

Although precision in XPS can be excellent, allowing small changes to be easily observed, obtaining accurate absolute elemental composition of a solid material from relative peak intensities is generally much more problematical, involving many factors: background removal; differing analysis depths at different photoelectron kinetic energies; possible angular distribution effects; calibration of the instrument transmission function, and variations of the distribution of the photoelectron intensity between “main” peaks (those usually used for analysis) and associated substructure following the main peak, as a function of the chemical bonding of the elements concerned. The last item, coupled with the use of photoionization cross-sections and/or relative sensitivity factors, is the major subject of this paper, though it is necessary to consider the other items also, using LiF as a test case. The results show that the above issues, which are relevant to differing degrees in most XPS analyses, present significant challenges to highly accurate XPS quantification.

Brundle, Christopher R.↗

Sulfuric Acid New Particle Formation Field Campaign Report

This campaign aimed to measure gaseous sulfuric acid concentrations at the surface level using a newly developed sulfuric acid dimethylamine-reactive condensation particle counter (SAD-RCPC). Sulfuric acid is a key nucleation precursor and is the likely driving compound for new particle formation at the Atmospheric Radiation Measurement (ARM) user facility’s Southern Great Plains (SGP) observatory. The SAD-RCPC will be deployed in the future at ARM SGP (AFC010010: "Vertically Resolved Nucleation Precursors at SGP") scheduled for 2024. The campaign was in conjunction with Eleanor Browne's campaign (AFC010075: Boundary-layer gradients in new particle formation) and took place from 05/7/23 to 5/21/23. The measured sulfuric acid concentrations by the SAD-RCPC will be compared to concentrations measured by Browne's chemical ionization mass spectrometer (CIMS). The SAD-RCPC consisted of a glass flow reactor coupled to a 1-nm diethylene glycol condensation particle counter (DEG CPC). Sampled air was mixed with a high concentration of dimethylamine vapor (~1 ppbv). Sulfuric acid within the sampled gas nucleated with dimethylamine to form detectable 1-nm particles. The measured particle concentration was converted to sulfuric acid concentration using nucleation reaction kinetics. The pulse height distribution of particles detected by the DEG CPC was observed using a multi-channel analyzer. Large pre-existing particles were separated from the freshly formed particles in the SAD-RCPC from the pulse height distribution (PHD). During the short campaign, seven new particle formation events were observed. The SAD-RCPC detected an increase in sulfuric acid concentration, but this also coincided with an increase in pre-existing <3 nm particle concentrations. In addition, Browne's CIMS detect sulfuric acid concentrations in the range of 10 5 -10 6 cm -3 .

54 ENVIRONMENTAL SCIENCES↗