Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model compounds”

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 1,351 records · Page 75

Next Generation Co-Molded One-Piece Automotive Parts

The purpose of this project was to determine the feasibility of co-molding Class A Sheet Molding Compound (SMC), structural SMC, and continuous fiber prepreg materials to produce a single piece co-molded automotive part, such as a hood. The combination of the three molding materials and the incorporation of selective design features (ribs, flanges, corrugations) was expected to eliminate the need for inner reinforcement panels, significantly reducing tooling costs and simplifying the manufacturing process. A multi-material solution would also result in significant weight savings. The scope of this project included the material characterization of the three materials, development of cure and flow simulation models, and validation of the models against parts molded with different combinations of the materials on an 11”x11” plaque tool with rib features. The intent of this project was to apply the learnings obtained from co-molding a part with simplified geometry to a Phase 2 project that would produce a single piece hood, co-molded with the same materials. Resins from INEOS Composites were provided to IDI to produce SMC and continuous fiber prepregs. Purdue University and INEOS Composites characterized the rheological and curing behavior of the different materials as well as the mechanical properties of the co-molded parts. This data was used by Purdue University to create simulation models. Models were created to predict both flow patterns and predict mechanical properties of different laminate constructions in and around the ribs. Michigan State University - Corktown validated Purdue’s models by co-molding the three materials in different combinations on a tool containing rib features provided by Century Tool. The co-molded parts were evaluated against the simulation models by Purdue, Corktown and INEOS. The project team demonstrated the following: • The ability to obtain a cohesive co-molded structure made up of Class A SMC, Structural SMC, and continuous fiber prepreg. • The ability to obtain a co-molded part with a Class A surface • Modeling of flow and fiber orientation • Modeling to predict mechanical properties of multi-material co-molded structures The predicted flow behavior and fiber orientation from simulations were then compared with the molded samples. Exact local orientation state was difficult to compare between microscopy and flow simulation, but captured general trends. Multi-material flow behavior was generally modeled well with SMC materials, but the introduction of woven material sheets requires further model development. Purdue compared the predicted versus actual mechanical properties, and INEOS Composites evaluated the surface appearance of unpainted and painted co-molded parts. The models developed by Purdue demonstrated that the mechanical properties can be predicted for parts with varying material configurations. The resulting models can be applied to future co-molded part design, tooling design, and molding conditions.

36 MATERIALS SCIENCE↗

Influence of the Molecular Structure on the Electrocatalytic Hydrogenation of Carbonyl Groups and H 2 Evolution on Pd

We investigated the electrocatalytic hydrogenation (ECH) of model aldehydes and ketones over carbon supported Pd in aqueous phase. We propose reaction mechanisms based on kinetic measurements and on spectroscopic and electrochemical characterization of the working catalyst. The reaction rates of ECH and of the H 2 evolution reaction (HER) vary with the applied electric potential following trends that strongly depend on the organic substrate. The intrinsic rates of hydrogenation and H 2 evolution are influenced, in opposing ways, by the sorption of the reacting organic substrate. Strong interactions, i.e., higher standard free energies of adsorption of the organic compound, induce high hydrogenation rates but decrease the overall currents. The fast hydrogenation kinetics produces a hydrogen-depleted environment that kinetically hinders the HER and the bulk phase transition of Pd to a H-rich bulk Pd hydride, which is triggered by the applied potential in the absence of reacting organic compounds. As consequence of strong organic-metal interactions, hydrogenation dominates at low overpotential. However, the coverages of organic substrates on the metal surface decrease and the rates of H 2 evolution surpass those of hydrogenation with increasingly negative electric potential. We determined the range of electric potential favoring hydrogenation on Pd and quantitatively deconvoluted the effects of the sorption of the organic compound, and of the rates of proton coupled electron transfers, on the kinetics of both ECH and HER. The results indicate that electrocatalysis offers hydrogenation pathways for polar molecules that are different and, in some cases, faster than those dominating in the absence of an external electric potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic and magnetic properties of the RuX3 (X = Cl, Br, I) family: two siblings—and a cousin?

Abstract Motivated by reports of metallic behavior in the recently synthesized RuI 3 , in contrast to the Mott-insulating nature of the actively discussed α -RuCl 3 , as well as RuBr 3 , we present a detailed comparative analysis of the electronic and magnetic properties of this family of trihalides. Using a combination of first-principles calculations and effective-model considerations, we conclude that RuI 3 , similarly to the other two members, is most probably on the verge of a Mott insulator, but with much smaller magnetic moments and strong magnetic frustration. We predict the ideal pristine crystal of RuI 3 to have a nearly vanishing conventional nearest-neighbor Heisenberg interaction and to be a quantum spin liquid candidate of a possibly different kind than the Kitaev spin liquid. In order to understand the apparent contradiction to the reported resistivity ρ , we analyze the experimental evidence for all three compounds and propose a scenario for the observed metallicity in existing samples of RuI 3 . Furthermore, for the Mott insulator RuBr 3, we obtain a magnetic Hamiltonian of a similar form to that in the much-discussed α -RuCl 3 and show that this Hamiltonian is in agreement with experimental evidence in RuBr 3 .

36 MATERIALS SCIENCE↗

SDE_quark

This program provides a fixed-grid quadrature algorithm to compute integrals in the self-energy of the quark propagator within the Maris-Tandy model. The quark propagator both on t he spacelike real axis and at complex-valued momenta is determined from its Schwinger-Dyson equation (SDE). We first apply an iterative solver to find the quark propagator on the spacelike real axis. The propagator at complex-valued momenta is then computed from its self-energy based on this solution, where demanding integrals are encountered. In order to compute of these integrals, we apply customized variable transformations for the radial integral after subtracting the asymptotics. We subsequently apply an optional compound of quadrature rules for the angular integral.

Jia, Shaoyang↗

Nitrogen-doped char as a catalyst for wet oxidation of phenol-contaminated water

Catalytic wet oxidation (CWO) of aqueous effluents rich in organic compounds is a very promising technology for the treatment of liquid wastes from biomass conversion processes. CWO reactions occur through the formation of free radical species, produced in the presence of an oxidant, which act on organic contaminates in the effluent. Although the reaction is well known, there exists a lack of affordable catalysts to conduct this process at the lower temperatures and pressures in novel bioenergy processes. Here, this study assessed the catalytic effect of nitrogen-doped chars as such an option. Phenol in aqueous solution was used as a model waste effluent. Treatment was conducted at moderate temperatures (190 to 260 °C), oxygen partial pressure of 1 MPa, and reaction times of 15, 30, and 45 min in stainless steel and glass-lined tube reactors. High pressure liquid chromatography (HPLC) analyses of the products quantified phenol and by-product concentrations used in the calculation of reaction activation energy. The char catalyst was studied by X-ray photoelectron spectroscopy (XPS), transmission electron microscopy (TEM), and scanning electron microscopy (SEM) in order to gain insight into its structure and surface composition. The results indicate that nitrogen-doped char catalysts accelerate the oxidation of phenol by decreasing its reaction activation energy from 82.2 kJ/mol (non-catalyzed) to 40.4 kJ/mol (catalyzed). An analysis from first principles using density functional theory (DFT) was conducted to ascertain which N functional group has the most significant impact on free radical formation in the presence of oxygen. Among all the N functional groups studied, the dipyridinic functional groups showed the most promising characteristics to facilitate the formation of hydroxyl free radicals.

Advanced Oxidation Processes for Organic Pollutant↗

Numerical simulation for the design of induction heating based radio frequency reactor for ethylene production

Ethylene is a vital petrochemical compound produced in vast amounts yearly by manufacturers that have enough scale to overcome the inherent thermodynamic inefficiencies of the process. In order to address the inefficiencies that prevent smaller scale or intermittent production ethylene, investigation of new production methods are required. In this work, we investigate the use of a radio frequency (RF) based reactor system that generates heat internally as opposed to applying heat externally via steam or direct combustion of fossil fuels. In order to guide the design of an electromagnetic based reactor system, we have created a macroscale model capable of capturing heat generation at the susceptors from a produced electromagnetic field and subsequent chemical reactions. Several susceptor and induction coil designs were investigated to better understand the performance of the system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spin-coupling topology in the copper hexamer compounds A 2 Cu 3 O(SO 4 ) 3 ( A =Na, K)

The compounds A 2 Cu 3 O(SO 4 ) 3 (A=Na,K) are characterized by copper hexamers that are weakly coupled to realize antiferromagnetic order below TN≈3K. They constitute quantum spin systems with S=1 triplet ground states. In this work, we investigated the energy-level splittings of the copper hexamers by inelastic neutron scattering experiments covering the entire range of the magnetic excitation spectra. Additionally, the observed transitions are governed by very unusual selection rules that we ascribe to the underlying spin-coupling topology. This is rationalized by model calculations that allow an unambiguous interpretation of the magnetic excitations concerning both the peak assignments and the nature of the spin-coupling parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Lignin conversion to β-ketoadipic acid by Pseudomonas putida via metabolic engineering and bioprocess development

Bioconversion of a heterogeneous mixture of lignin-related aromatic compounds (LRCs) to a single product via microbial biocatalysts is a promising approach to valorize lignin. Here, Pseudomonas putida KT2440 was engineered to convert mixed p-coumaroyl– and coniferyl-type LRCs to β-ketoadipic acid, a precursor for performance-advantaged polymers. Expression of enzymes mediating aromatic O-demethylation, hydroxylation, and ring-opening steps was tuned, and a global regulator was deleted. β-ketoadipate titers of 44.5 and 25 grams per liter and productivities of 1.15 and 0.66 grams per liter per hour were achieved from model LRCs and corn stover-derived LRCs, respectively, the latter representing an overall yield of 0.10 grams per gram corn stover-derived lignin. Technoeconomic analysis of the bioprocess and downstream processing predicted a β-ketoadipate minimum selling price of $\$2.01$ per kilogram, which is cost competitive with fossil carbon-derived adipic acid ($\$1.10$ to 1.80 per kilogram). Overall, this work achieved bioproduction metrics with economic relevance for conversion of lignin-derived streams into a performance-advantaged bioproduct.

09 BIOMASS FUELS↗

Multi-Probe 24A Pre-Shot report

The Multi-Probe Radiography campaign will field a full day of experiments at the Omega EP laser facility. The day consists of shots alternating between Omega EP’s two short-pulse laser beams. The Backlighter beam will generate proton and deuteron beams from 500-800nm CH/CD film targets, and the Sidelighter beam will accelerate electrons to generate x-rays from 0.5mm x 25-50µm Ø Ta wire targets attached to Compound Parabolic Concentrator (CPC) cones. The ion beams shots will optimize CD foil thickness (maximizing ion energy/yield) for the transparency regime for use as the pitcher in a pitcher catcher neutron-generation platform. Several shots will include a LiF catcher within the NTA and this feeds into future neutron radiography work. X-ray shots will be used to characterize the platform established in FY2020B via measurement of electron spectra with and without the Ta wire, and radiography will be conducted at varying laser energies. This will validate PIC electron acceleration modeling while also demonstrating potential for using the x-ray platform for imaging of thicker objects.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Predicting receptor-ligand pairing preferences in plant-microbe interfaces via molecular dynamics and machine learning

Microbiome assembly, structure, and dynamics significantly influence plant health. Secreted microbial signaling molecules initiate and mediate symbiosis by binding to structurally compatible plant receptors. For example, lipo-chitooligosaccharides (LCOs), produced by nitrogen-fixing rhizobial bacteria and various fungi, are recognized by plant lysin motif receptor-like kinases (LysM-RLKs), which activate the common symbiotic pathway. Accurately predicting these molecular interactions could reveal complementary signatures underlying the initial stages of endosymbiosis. Despite the breakthrough in protein-ligand structure prediction with deep learning-based tools, such as AlphaFold3, the large size and highly flexible nature of signaling compounds like LCOs present major challenges for detailed structural characterization and binding-affinity prediction. Typical structure-/physics-based methods of ligand virtual screening are designed for small, drug-like molecules, often rely on high-resolution, experimentally determined structures of the protein receptors, and rarely achieve sufficient sampling to obtain converged thermodynamic quantities with large ligands. In this study, we developed a hybrid molecular dynamics/machine learning (MD/ML) approach capable of predicting binding affinity rankings with high accuracy in systems involving large, flexible ligands, despite limited experimental structural information. Using coarse initial structural models, the predictions using the MD/ML workflow achieved strong alignment with experimental trends, particularly in the top-affinity tier for four legume LysM-RLKs (LYR3) binding to LCOs and a chitooligosaccharide. Furthermore, the MD-based conformation selection protocol provided critical structural insights into substrate specificity and binding mechanisms. This study demonstrates a powerful method to screen for challenging cognate ligand-receptors and advance our understanding of the molecular basis of microbial colonization in plants.

Lipo-chitooligosaccharides↗

3D Printed Parahydrophobic Surfaces as Multireaction Platforms

We report Parahydrophobic surfaces (PHSs) composed of arrays of cubic μ-pillars with a double scale of roughness and variable wettability were systematically obtained in one step and a widely accessible stereolithographic Formlabs 3D printer. The wettability control was achieved by combining the geometrical parameters (H = height and P = pitch) and the surface modification with fluoroalkyl silane compounds. Homogeneous distribution of F and Si atoms onto the pillars was observed by XPS and SEM-EDAX. A nano-roughness on the heads of the pillars was achieved without any post-treatment. The smallest P values lead to surfaces with static contact angles (CAs) >150° regardless of the H utilized. Interestingly, the relationship 0.6 ≤ H/P ≤ 2.6 obtained here was in good agreement with the H/P values reported for nano- and submicron pillars. Furthermore, experimental CAs, advancing and receding CAs, were consistent with the theoretical prediction from the Cassie–Baxter model. Structures covered with perfluorodecyltriethoxysilane with high H and short P lead to PHSs. Conversely, structures covered with perfluorodecyltrimethoxysilane exhibited a superhydrophobic behavior. Finally, several aqueous reactions, such as precipitation, coordination complex, and nanoparticle synthesis, were carried out by placing the reactive agents as microdroplets on the parahydrophobic pillars, demonstrating the potential application as chemical multi-reaction array platforms for a large variety of relevant fields in microdroplet manipulation, microfluidics systems, and health monitoring, among others.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quasilinear dispersion in electronic band structure and high Seebeck coefficient in CuFeS 2 -based thermoelectric materials

The earth-abundant natural mineral chalcopyrite CuFeS 2 is a potential n-type thermoelectric material because of its large Seebeck coefficient at high carrier concentrations. For a long time, the large Seebeck coefficient of CuFeS 2 has been attributed to a large electron effective mass, but the reasons for this and the unusual carrier concentration dependent behavior have rarely been discussed. Here, we systematically investigated the special transport behavior of CuFeS 2 and found the classical parabolic band model to be inadequate in explaining it. Our experimental and theoretical studies indicate that there are two flat electronic pockets at the and Z points of the Brillouin zone near the conduction band edge of CuFeS 2 that dominate the charge transport. These electronic pockets result from nonparabolic quasilinearly dispersing bands that give rise to a linear wave vector dependent energy (E ~ k) and a carrier density dependent effective mass (m* ~ m 0 n 1/3 ). Such a strong carrier concentration dependent carrier effective mass results in the high Seebeck coefficient of CuFeS 2 compound under a large carrier density. Here, the work demonstrates that quasilinearly dispersing bands can give strongly enhanced Seebeck coefficient, and could be useful in optimizing the properties of thermoelectric materials.

36 MATERIALS SCIENCE↗

Effects of Selenium Doping in Zinc Telluride from First Principles

We report point defects and compensation in zinc telluride have been subjects of study for several decades due to the compound’s scientific/technological applications. Despite this, there remains some disagreement in the literature as to the fine details of its defect chemistry. Recently, it has also come under renewed interest due to the emerging Cd( 1–x )Zn x Te( 1–y )Se y material family, used in nuclear radiation detector applications. In this article, we use state-of-the-art computational tools to re-examine the native defect chemistry of ZnTe and explore how it changes upon addition of selenium using density functional theory with screened hybrid exchange–correlation functionals, thermochemical models of the equilibration environment, and statistical mechanics models of the point defect ensemble.

II-VI semiconductors↗

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kondo effect in ferromagnetic quantum critical CeRh 6 ⁢Ge 4

The mechanism of a pressure-induced quantum critical point in the heavy fermion ferromagnet CeRh 6 ⁢Ge 4 has attracted interest, as ferromagnetic quantum criticality in a clean itinerant Ce compound is typically avoided. The localized versus itinerant character of the 4⁢𝑓 electrons is a key aspect for understanding this behavior. We investigated the electronic structure of the 4⁢𝑓 shell in CeRh 6 ⁢Ge 4 using core-level photoelectron and x-ray absorption spectroscopy, demonstrating the hybridization of Ce 4⁢𝑓 with the conduction electrons. Linearly polarized x-ray absorption reveals a temperature-dependent linear dichroism consistent with the crystal-electric-field sequence as inferred from the static susceptibility. This dichroism cannot be described by an ionic full-multiplet model alone, but is reproduced by including the Kondo effect within a single-impurity Anderson model in the noncrossing approximation. The Kondo effect mixes higher-lying crystal-field states into a resulting multiorbital ground state with 4⁢𝑓 occupancy, 𝑛 𝑓 ∼ 0.9. Deviations at low temperatures between the measured linear dichroism and calculated dichroism suggest an orbital-dependent Kondo effect. A scenario in which there is a multiorbital ground state and orbital-dependent Kondo hybridization should be a starting point for a model of pressure-induced criticality in CeRh 6 ⁢Ge 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Enhancing superconductivity of Y 5 Rh 6 Sn 18 by atomic disorder

Here, we investigate the effect of enhancement of superconducting transition temperature $T_c$ by nonmagnetic atom disorder in the filled skutterudite-related tetragonal Y 5 Rh 6 Sn 18 compound doped with Ca. We documented experimentally that Y 5 Rh 6 Sn 18 and its Ca-doped alloys are electronically inhomogeneous at the nanoscale, when Ca content is smaller than ~1.5 in the Y 5– x Ca x Rh 6 Sn 18 system, while for $\textit{x}$ ≥ 1.5 much stronger chemical phase inhomogeneity is observed with an about 20% volume fraction of the second 3:4:13 cubic phase, which we interpret as a fluctuation of stoichiometry within the bulk sample. Then the enhancement of $T_c$ vs $\textit{x}$ could be modeled by a mechanism proposed in recent theoretical reports for increasing the mean-field transition temperature $T_c$ in the presence of nonmagnetic disorder. The increase in disorder with doping increases the sample inhomogeneity and causes a systematic increase in $T_c$, while for two-phase $\textit{x}$ ≥ 1.5 samples the critical temperature $T^{\star}_c ≈ 2 × T_c$ rapidly increases. Based on band structure calculations performed under pressure, we demonstrate how the change in density of states would affect $T_c$ of Y 5 Rh 6 Sn 18 . We obtained the Grüneisen parameter $γ_G$ larger for the inhomogeneous phase with respect to $γ_G$ of the bulk $T_c$ phase and attribute the enhancement of $T_c$ to larger stiffening of the high-temperature $T^{\star}_c$ phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Thin film combinatorial sputtering of Al-Ce alloys: Investigating the phase separation of as-deposited solid solutions and determining the coefficient of thermal expansion

Al x Ce 100–x thin films with a composition range of ~75.0 < x < 99.5 at% (36.5 < x < 97.5 wt%) were synthesized via combinatorial co-sputtering from an Al and an Al 50 Ce 50 target. The crystal structure, phase fraction, film morphology, electrical resistivity, and temperature-dependent coefficients of thermal expansion (CTE) are all correlated to the Al x Ce 100–x composition. The as-deposited films form a metastable solid-solution, and annealing leads to the formation of the thermodynamically stable two-phase system of Al and the α-Al 11 Ce 3 intermetallic. Temperature dependent x-ray diffraction (XRD) reveals that the two phases expand independently of one another, and the thin film Al temperature-dependent CTE is similar to bulk Al. The thin film Al 11 Ce 3 intermetallic phase has a nearly constant CTE of ~1.5 × 10 –5 /°C within the temperature range studied (25–550 °C). To confirm the thin film Al 11 Ce 3 results, bulk stoichiometric Al 11 Ce 3 and +/- 1 wt% Ce samples were prepared and the CTE of each was measured with the same conditions. A Rietveld analysis of the bulk data enabled an estimation of the CTE in each of the 3 orthorhombic lattice parameters, which displayed anisotropic behavior. The thin film and bulk CTE measurements were in very good agreement. Estimations of the temperature dependent CTE of the two-phase alloys are made via the Reuss and Voigt models. Finally, by demonstrating the efficacy of the approach, more complex multi-component rapid materials discovery of low CTE Al-alloys can be pursued via the combinatorial thin film synthesis and XRD measurement.

36 MATERIALS SCIENCE↗

Chamber Flux and Porewater Concentration of CH4, CO2 and N2O, 2018, Columbia River bank at the Hanford site, WA, USA

This is the corresponding observations data for the paper: Jorge A Villa; Garrett J Smith; Yang Ju; Lupita Renteria; Jordan C Angle; Evan Arntzen; Samuel F Harding; Huiying Ren; Xingyuan Chen; Audrey H Sawyer; Emily B Graham; James C Stegen; Kelly C Wrighton; Gil Bohrer (2020) Methane and nitrous oxide porewater concentrations and surface fluxes of a regulated river. Science of the Total Environment.Greenhouse gas (GHG) emissions from rivers are a critical missing component of current global GHG models. Their exclusion is mainly due to a lack of in-situ measurements and a poor understanding of the spatiotemporal dynamics of GHG production and emissions, which prevents optimal model parametrization. We combined simultaneous observations of porewater concentrations along different beach positions and depths, and surface fluxes of methane and nitrous oxide at a plot scale in a large regulated river during three water stages: rising, falling, and low. Our goal was to gain insights into the interactions between hydrological exchanges and GHG emissions and elucidate possible hypotheses that could guide future research on the mechanisms of GHG production, consumption, and transport in the hyporheic zone (HZ). Results indicate that the site functioned as a net source of methane. Surface fluxes of methane during river water stages at three beach positions (shallow, intermediate and deep) correlated with porewater concentrations of methane. However, fluxes were significantly higher in the intermediate position during the low water stage, suggesting that low residence time increased methane emissions. Vertical profiles of methane peaked at different depths, indicating an influence of the magnitude and direction of the hyporheic mixing during the different river water stages on methane production and consumption. The site acted as either a sink or a source of nitrous oxide depending on the elevation of the water column. Nitrous oxide porewater concentrations peaked at the upper layers of the sediment throughout the different water stages. River hydrological stages significantly influenced porewater concentrations and fluxes of GHG, probably by influencing heterotrophic respiration (production and consumption processes) and transport to and from the HZ. Our results highlight the importance of including dynamic hydrological exchanges when studying and modeling GHG production and consumption in the HZ of large rivers.

54 ENVIRONMENTAL SCIENCES↗