Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase velocity”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

ZTF SN Ia DR2: Improved SN Ia colors through expanded dimensionality with SALT3+

Context. Type Ia supernovae (SNe Ia) are a key probe in modern cosmology, as they can be used to measure luminosity distances at gigaparsec scales. Models of their light curves are used to project heterogeneous observed data onto a common basis for analysis. Aims. The SALT model currently used for SN Ia cosmology describes SNe as having two sources of variability, accounted for by a color parameter c , and a “stretch” parameter x 1 . We extend the model to include an additional parameter we label x 2 , to investigate the cosmological impact of currently unaddressed light-curve variability. Methods. We constructed a new SALT model, that we dub “SALT3+”. This model was trained by an improved version of the SALTshaker code, using training data combining a selection of the second data release of cosmological SNe Ia from the Zwicky Transient Facility and the existing SALT3 training compilation. Results. We find additional, coherent variability in supernova light curves beyond SALT3. Most of this variation can be described as phase-dependent variation in g − r and r − i color curves, correlated with a boost in the height of the secondary maximum in i -band. These behaviors correlate with spectral differences, particularly in line velocity. We find that fits with the existing SALT3 model tend to address this excess variation with the color parameter, leading to less informative measurements of supernova color. We find that neglecting the new parameter in light-curve fits leads to a trend in Hubble residuals with x 2 of 0.039 ± 0.005 mag, representing a potential systematic uncertainty. However, we find no evidence of a bias in current cosmological measurements. Conclusions. We conclude that extended SN Ia light-curve models promise mild improvement in the accuracy of color measurements, and corresponding cosmological precision. However, models with more parameters are unlikely to substantially affect current cosmological results.

Kenworthy, W. D. (ORCID:0000000251535983)↗
DOE OSTI2023

Exploring the stellar populations of backsplash galaxies

Backsplash galaxies are those that traverse and overshoot cluster cores as they fall into these structures. They are affected by environment, and should stand out in contrast to the infalling population. We target galaxies in the vicinity of clusters (R≳R 200 ) and select a sample in projected phase space (PPS), from the compilation of Sampaio et al. based on SDSS data. We present a statistical analysis, comparing two regions in PPS, with the same projected distance to the cluster but different velocity. The analysis relies on the presence of variations in the stellar population content of backsplash galaxies. We find a lower limit in the fractional contribution of ∼5 per cent with respect to the general sample of infalling galaxies at similar group-centric distance when using single line strength analysis, or ∼15-30 per cent when adopting bivariate distributions. The stellar populations show a subtle but significant difference towards older ages, and a higher fraction of quiescent galaxies. We also compare this set with a general field sample, where a substantially larger difference in galaxy properties is found, with the field sample being consistently younger, metal poorer and with a lower fraction of quiescent galaxies. Noting that our ‘cluster’ sample is located outside of the virial radius, we expect this difference to be caused by pre-processing of the infalling galaxies in the overall higher density regions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

SILCC – VII. Gas kinematics and multiphase outflows of the simulated ISM at high gas surface densities

ABSTRACT We present magnetohydrodynamic (MHD) simulations of the star-forming multiphase interstellar medium (ISM) in stratified galactic patches with gas surface densities Σgas = 10, 30, 50, and 100 $\mathrm{M_\odot \, pc^{-2}}$. The silcc project simulation framework accounts for non-equilibrium thermal and chemical processes in the warm and cold ISM. The sink-based star formation and feedback model includes stellar winds, hydrogen-ionizing UV radiation, core-collapse supernovae, and cosmic ray (CR) injection and diffusion. The simulations follow the observed relation between Σgas and the star formation rate surface density ΣSFR. CRs qualitatively change the outflow phase structure. Without CRs, the outflows transition from a two-phase (warm and hot at 1 kpc) to a single-phase (hot at 2 kpc) structure. With CRs, the outflow always has three phases (cold, warm, and hot), dominated in mass by the warm phase. The impact of CRs on mass loading decreases for higher Σgas and the mass loading factors of the CR-supported outflows are of order unity independent of ΣSFR. Similar to observations, vertical velocity dispersions of the warm ionized medium (WIM) and the cold neutral medium (CNM) correlate with the star formation rate as $\sigma _\mathrm{z} \propto \Sigma _\mathrm{SFR}^a$, with a ∼ 0.20. In the absence of stellar feedback, we find no correlation. The velocity dispersion of the WIM is a factor ∼2.2 higher than that of the CNM, in agreement with local observations. For $\Sigma _\mathrm{SFR} \gtrsim 1.5 \times 10^{-2}\, \mathrm{M}_\odot \, \mathrm{yr}^{-1}\, \mathrm{kpc}^{-2}$ the WIM motions become supersonic.

Astronomy & Astrophysics↗
DOE OSTI2025

Blending as an In-Tank Waste Processing Strategy: Year 1 Progress Report

The Hanford Site stores approximately 56 million gallons of radioactive legacy defense waste in underground tanks. Direct feed of southeast quadrant tank waste to a high-level waste treatment facility is currently being explored to accelerate tank treatment activities. This study seeks to de-risk waste transfers by favorably altering the resuspension and critical velocities of wastes through strategic waste blending. Wastes rich in aluminum- and zirconium-bearing solids pose significant risks due to their high settled strength and density. A series of tests explores how solid phase particle size and suspending phase chemistry alter waste transport properties, quantified through the just-suspended mixing speed (NJS), for pure and blended aluminum- and zirconium-bearing waste solid simulants. The study utilizes a range of gibbsites (three sizes spanning 10 to 90 µm) and zirconia solids (three sizes spanning 10 to 100 µm). As expected, the NJS of pure solids largely follows predictions based on literature correlations; blends of gibbsite and zirconium-bearing solids of differing particle sizes exhibit NJS maxima attributed to dense packing of the mixed-size settled solids. The findings of this study will assist in developing strategies to optimize the transport of waste solids for the safe and efficient delivery of tank waste.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Ultrasonic characterization of Ti-5Al-5Mo-5V-3Cr

Metastable β titanium alloys such as Ti-5Al-5Mo-5V-3Cr are of increasing interest due to their excellent corrosion resistance, high strength, and in particular, the ability to manipulate microstructure to a high degree to control the mechanical properties of the material. It is therefore important to consider how this manipulation of mechanical properties will affect the inspection of this material with ultrasonic testing both for the purposes of flaw detection and material characterization. The connection between material properties and ultrasonic characterization techniques will be discussed and then dissected in the context of Ti-5553. This metastable alloy was shown to exhibit remarkable variation in microstructure as well as ultrasonic velocity and attenuation demonstrating these techniques as potential candidates for characterizing Ti-5553 and Ti-5553 as a good candidate for aiding in future studies aimed at learning more about the role of secondary phases in ultrasonic wave propagation.

Sunderman, Ruth↗
DOE OSTI2020

Topological Dirac Semimetal Phase in Bismuth Based Anode Materials for Sodium-Ion Batteries

Bismuth has recently attracted interest in connection with Na-ion battery anodes due to its high volumetric capacity. It reacts with Na to form Na 3 Bi which is a prototypical Dirac semimetal with a nontrivial electronic structure. Density-functional-theory based first-principles calculations are playing a key role in understanding the fascinating electronic structure of Na 3 Bi and other topological materials. In particular, the strongly-constrained-and-appropriately-normed (SCAN) meta-generalized-gradient-approximation (meta-GGA) has shown significant improvement over the widely used generalized-gradient-approximation (GGA) scheme in capturing energetic, structural, and electronic properties of many classes of materials. Here, we discuss the electronic structure of Na 3 Bi within the SCAN framework and show that the resulting Fermi velocities and s-band shift around the Γ point are in better agreement with experiments than the corresponding GGA predictions. SCAN yields a purely spin-orbit-coupling (SOC) driven Dirac semimetal state in Na 3 Bi in contrast with the earlier GGA results. Our analysis reveals the presence of a topological phase transition from the Dirac semimetal to a trivial band insulator phase in Na 3 Bi x Sb 1-x alloys as the strength of the SOC varies with Sb content, and gives insight into the role of the SOC in modulating conduction properties of Na 3 Bi.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Sound Velocities in Vanadium Reveal Complex Elastic Behavior at High Pressures

Compressional (VP) and shear (VS) wave velocities of polycrystalline vanadium were measured simultaneously up to 11.5 GPa at room temperature using ultrasonic interferometry in a multi-anvil press. Complex softening behavior in VS and resulting shear moduli are discovered, possibly revealing a precursor to the reported phase transition within 30–60 GPa. The current data enables a comprehensive assessment of the elastic and mechanical properties of vanadium at high pressures, including bulk and shear moduli, Young’s modulus, Poisson’s ratio, and Pugh’s ratio. Through fitting to the 3rd-order finite strain equations, the elastic moduli and their pressure derivatives were determined to be K S0 = 151 (2) GPa, G 0 = 46.9 (8) GPa, K$^{’}_{S0}$ = 3.47 (5), and G$^{’}_{0}$ = 0.62 (1). These experimental results allow us to compare with and benchmark the existing Steinberg–Guinan models for extrapolations to extreme pressure and temperature conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Signatures of Coherent Phonon Transport in Ultralow Thermal Conductivity Two-Dimensional Ruddlesden–Popper Phase Perovskites

An emerging class of methylammonium lead iodide (MAPbI 3 )-based Ruddlesden-Popper (RP) phase perovskites, BA 2 MA n-1 Pb n I 3n+1 (n = 1-7), exhibit enhanced stability to environmental conditions relative to MAPbI 3 , yet still degrade at elevated temperatures. Here, we experimentally determine the thermal conductivities of these layered RP phases for n = 1-6, where n defines the number of repeated perovskite octahedra per layer. We measure thermal conductivities of 0.37 ± 0.13/0.12, 0.17 ± 0.08/0.07, 0.21 ± 0.05/0.04, and 0.19 ± 0.04/0.03 W/m∙K in thin films of n = 1-4 and 0.08 ± 0.06/0.04, 0.06 ± 0.04/0.03, 0.06 ± 0.03/0.03, and 0.08 ± 0.07/0.04 W/m∙K in single crystals of n = 3-6. With the exception of n = 1, these thermal conductivities are lower than the range of 0.34-0.50 W/m∙K reported for single-crystal MAPbI 3 . Reduced-order lattice dynamics modeling suggests that the initially decreasing trend of thermal conductivity in similarly oriented perovskites with increasing n may result from the transport properties of coherent phonons, emergent from the superstructure, that do not scatter at the interfaces of organic butylammonium chains and perovskite octahedra. Reduced group velocity of coherent phonons in n = 3-6, a consequence of band flattening in the phonon dispersion, is primarily responsible for their ultralow thermal conductivities. Similar effects on thermal conductivity have been experimentally demonstrated in deposited superlattices, but never in naturally defined materials such as RP phases. GIWAXS measurements reveal that higher n RP phase thin films are less orientationally controlled and therefore possess apparently elevated thermal conductivities relative to single crystals of the same n.

2D perovskites↗
DOE OSTI2022

On the integration of molecular dynamics, data science, and experiments for studying solvent effects on catalysis

Computational workflows that combine molecular dynamics (MD) simulations and emerging data-centric (DC) methods can accelerate the screening and analysis of solvent systems experimentally and computationally. Here, MD simulations provide atomic positions and velocities of reactant, solvent, and catalyst materials that can be manipulated into data representations that in turn can be used by DC techniques to conduct predictive modeling, feature extraction, and experimental design. For liquid-phase catalytic applications, emerging DC techniques such as Convolutional and Graph Neural Networks (CNN/GNN), Topological Data Analysis (TDA), and Active Learning (AL) can leverage MD and experimental data to quickly predict solvent effects on reaction outcomes. For instance, in recent studies, 3D solvent environments obtained with MD have been exploited by CNNs to predict experimental reaction rates for homogeneous acid-catalyzed lignocellulosic processes. In this perspective, we discuss basic principles of DC methods and how these can be combined with MD to enable high-throughput screening of solvent selection for diverse catalysis applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Dynamics of Downdrafts Around a Growing Convective Cloud: A Numerical Study

We examine the dynamics of cloud-edge downdrafts over the growth phase of isolated cumuli, combining Eulerian and Lagrangian analyses. As in previous studies, our results show that growing cumuli are surrounded by downdrafts linked to cloud-scale quasi-toroidal circulations at all times at middle and upper cloud levels consistent with the thermal chain description of convective clouds. These toroidal circulations are responsible for the most intense cloud-edge downdrafts in our simulations. In the upper cloud half, roughly 30%–50% of the upward mass flux is typically compensated within a radius of about twice the updraft radius in quasi-laminar simulations forced by a warm bubble in an initially quiescent flow. In a turbulent cloud forced by surface fluxes, this compensation fraction is around 10%–30% over the same region. In contrast to the buoyancy-centered view of subsiding shells, Eulerian and Lagrangian vertical momentum budget analyses show that the most intense cloud-edge downdrafts in the turbulent setup, and after spin-up of the toroidal circulation in the quasi-laminar experiments, are predominantly mechanically forced (i.e., driven by dynamic pressure accelerations). This is consistent throughout the entire growth phase of the cumulus clouds and across tests with varying assumptions, including drier and moister environments. Despite dynamic pressure perturbations being the main driver of toroidal downdrafts, the downdraft speed (relative to the corresponding updraft velocity) exceeds the prediction of the non-buoyant Hill's spherical vortex—a simple model frequently used for cloud-scale circulations—by more than 30%.

Pardo, Lianet Hernández [Goethe Univ., Frankfurt (↗
DOE OSTI2025

Equation of state measurement of detonation carbon condensates using optical microscopy and interferometry

Thermochemical models of detonation that estimate performance (e.g., detonation velocity, energy delivery, etc.), are based on assumptions that carbon condensates (soot) formed during detonation is largely similar to bulk carbon. However, soot constituents can range from amorphous carbon to nanodiamond and include other material phases. Since thermodynamic properties of the soot such as compressibility are imperative for accurate thermochemical modeling of detonation reaction chemistry, experimental measurements of the equation of state (EOS) which determine the compressibility are vital. Due to the mixed-phase nature of detonation soot, typical methods to measure the EOS (e.g., x-ray diffraction) are untenable. In this study, the high-pressure EOS up to 20 GPa was determined for detonation soot collected from PBX 9502, Composition B (Comp B), Hexanitrostilbene (HNS), and LX-21 high explosives by employing a direct volume technique using optical microscopy and interferometry in a diamond anvil cell. Comp B soot was determined to be the least compressible [K 0 = 57.9(17) GPa] with HNS soot [K 0 = 53.7(15) GPa], LX-21 soot [K 0 = 45.8(51) GPa], and PBX 9502 soot [K 0 = 28.2(27) GPa] being more compressible, likely due to differences in nanodiamond content as compared to amorphous carbon and graphite content.

Amorphous materials↗
DOE OSTI2025

XRISM spectroscopy on orbital modulation of Fe Ly α lines in Cygnus X-3

To understand physical processes such as mass transfer and binary evolution in X-ray binaries, the orbital parameters of the system are fundamental and crucial information. Cygnus X-3 is a high-mass X-ray binary composed of a compact object of unknown nature and a Wolf–Rayet star, which is of great interest in the context of wind-fed mass accretion and binary evolution. Here we present XRISM/Resolve high-resolution spectroscopy, focusing on the Fe Ly$\alpha$ lines when the source was in its hypersoft state. We perform an orbital phase-resolved spectral analysis of the lines to study the orbital modulation of the emission and absorption lines. It is found that the emission lines reflect the orbital motion of the compact object whose estimated velocity amplitude is $430^{+150}_{-140}$ km s$^{-1}$, while the absorption lines show a variation that can be interpreted as originating from the stellar wind. We discuss possible mass ranges for the binary components using the mass function with the estimated value of the velocity amplitude in this work, combined with the relation between the mass-loss rate and the orbital period derivative and the empirical mass and mass-loss rate relation for Galactic Wolf–Rayet stars. They are constrained to be (1.3–5.1) and (9.3–12) M$_\odot$, respectively, for the assumed inclination angle of $i =25^\circ$, and become more relaxed, (1.3–24) and (9.3–16) M$_\odot$, for $i = 35^\circ$. Thus, it remains unclear whether the system harbors a black hole or a neutron star.

X-rays: binaries↗
DOE OSTI2022

Heterogeneous phase transformation pathways in additively manufactured Al-Ce-Mn alloys

Heat treatment of additively manufactured Al-Ce based multicomponent alloys leads to complex microstructure evolution. In this research, the ability to extend the phase transformation theories involving nucleation of a product phase from a heterogeneous multi-phase microstructure typical to that of additively manufactured samples is explored. The Al-10Ce-8Mn (wt%) was used as a model alloy system. Under additive manufacturing conditions different solidification microstructures were obtained due to spatial and temporal variations of thermal gradients (G) and liquid-solid interface velocities (R) within a given melt pool. Near the melt pool boundary (high G and low R, referred as MPB region), initially, Al 20 Mn 2 Ce forms from the liquid followed by a eutectic of FCC Al and Al 11 Ce 3 . In the melt pool interiors (low G and high R referred as ES region) a eutectic structure between FCC Al and Al 20 Mn 2 Ce is observed. During subsequent heat treatments, the MPB and ES regions transform into different sets of microstructures. In the MPB region, a fine globular microstructure containing FCC Al, Al 11 Ce 3 , Al 6 Mn, and Al 12 Mn results from the decomposition of Al 20 Mn 2 Ce. In the ES region a faceted Al 51 Mn 7 Ce 4 plate phase results from the decomposition of Al 20 Mn 2 Ce. The formation of the Al 51 Mn 7 Ce 4 phase within the eutectic microstructure at the boundaries of FCC Al and Al 20 Mn 2 Ce has not been reported in the literature. Further, these two distinct phase transformation pathways are rationalized based on the role of driving force on the nucleation of (Al 6 Mn) and/or metastable intermetallic (Al 51 Mn 7 Ce 4 ) phases at the interface of aluminum (FCC) and the non-equilibrium intermetallic (Al 20 Mn 2 Ce) phases.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Influence of Solidification Velocity on Partitioning and Precipitate Evolution in a Multicomponent, Soft Magnetic Alloy Processed by Laser Melting

Multicomponent alloys have recently been shown to overcome the structural-functional barrier associated with soft magnetic materials. Here, single laser tracks have been used to investigate the effects of solidification velocity on the Al 7 Fe 32.6 Co 27.7 Ni 27.7 Ta 5 (at pct) alloy. Ta-rich cell wall partitioning was shown to affect the final distribution and size of the precipitating L1 2 phase, after aging. These findings demonstrate that solidification rate control offers the potential to tailor localized performance in alloys with size-dependent precipitation properties.

Brookins, Jeffrey M. [Oak Ridge National Laborator↗
DOE OSTI2025

Benchmark microgravity experiments and computations for 3D dendritic-array stability in directional solidification

In this study, we present a comprehensive quantitative analysis of stability bands for dendritic arrays during directional solidification of a transparent succinonitrile-0.46 wt % camphor alloy, spanning a broad range of pulling velocities. Taking advantage of the microgravity environment aboard the International Space Station where most convection effects are suppressed, we obtain unique measurements that quantify the stable primary spacing range of spatially extended three-dimensional dendritic array structures under purely diffusive growth conditions. Through carefully designed velocity jump experiments and detailed examination of sub-grain boundary dynamics, we characterize key instabilities, including elimination and tertiary branching, shedding new light on the mechanisms governing dynamic dendritic spacing selection in extended 3D arrays. Phase field simulations are performed to characterize the stability limits of dendritic array structures for quantitative comparison with the flight experiments. Although the simulations capture general trends, significant deviations are noted at the upper stability boundary, indicating the influence of additional, unexplored factors. These findings contribute to a deeper understanding of dendritic growth dynamics and offer valuable benchmark data that could aid in refining predictive models and improving control of dendritic microstructures in metallurgical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Flow-induced oscillation patterns for two tandem cylinders with turbulence stimulation and variable stiffness and damping

Herein, the oscillation patterns of two 1-DOF cylinder oscillators, undergoing VIV and galloping, are investigated in a free surface water channel for 3.2 × 104 ≤ Re ≤ 1.2 × 10 5 . The cylinders are arranged in tandem and supported by springs for a range of different spring stiffness and damping parameters. The efficiency of a current energy converter (CEC), based on flow-induced oscillation (FIO) of multiple cylinders in tandem, is critically related to the cylinder oscillation patterns. Due to the limited number of studies on the FIO for multiple cylindrical or prismatic bodies, their oscillation patterns have neither been identified nor classified. The surfaces of the cylinders are modified by turbulence stimulation to enhance FIO. Three different center-to-center spacing, five stiffness, and six damping ratios, for a total of 90 sets of experiments were conducted. The current velocity range is from 0.34 m/s to 1.32 m/s. From more than 2000 tests, five major patterns, nine sub-patterns are identified and classified. The patterns are defined based on the amplitudes, frequencies and the phase angle differences between the two cylinders. The characteristics and mechanics of each oscillation pattern are explored and explained from the perspective of fluid-structure interaction (FSI). By systematically varying the parameters, the underlying hydrodynamic mechanisms, including the coupling level between vortices and cylinders, the shielding effect, the wake effect, and the stability states are revealed. Few preliminary observations on the connection between oscillation pattern and harnessed power by the tandem cylinders are reported.

42 ENGINEERING↗
DOE OSTI2020

Continuous and discontinuous transitions in the depinning of two-dimensional dusty plasmas on a one-dimensional periodic substrate

We report Langevin dynamical simulations are performed to study the depinning dynamics of two-dimensional dusty plasmas on a one-dimensional periodic substrate. From the diagnostics of the sixfold coordinated particles P 6 and the collective drift velocity V x , three different states appear, which are the pinning, disordered plastic flow, and moving ordered states. It is found that the depth of the substrate is able to modulate the properties of the depinning phase transition, based on the results of P 6 and V x , as well as the observation of hysteresis of V x while increasing and decreasing the driving force monotonically. When the depth of the substrate is shallow, there are two continuous phase transitions. When the potential well depth slightly increases, the phase transition from the pinned to the disordered plastic flow states is continuous; however, the phase transition from the disordered plastic flow to the moving ordered states is discontinuous. When the substrate is even deeper, the phase transition from the pinned to the disordered plastic flow states changes to discontinuous. When the depth of the substrate further increases, as the driving force increases, the pinned state changes to the moving ordered state directly, so that the disordered plastic flow state disappears completely.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Urbanization-induced land and aerosol impacts on sea-breeze circulation and convective precipitation

Abstract. Changes in land cover and aerosols resulting from urbanization may impact convective clouds and precipitation. Here we investigate how Houston urbanization can modify sea-breeze-induced convective cloud and precipitation through the urban land effect and anthropogenic aerosol effect. The simulations are carried out with the Chemistry version of the Weather Research and Forecasting model (WRF-Chem), which is coupled with spectral-bin microphysics (SBM) and the multilayer urban model with a building energy model (BEM-BEP). We find that Houston urbanization (the joint effect of both urban land and anthropogenic aerosols) notably enhances storm intensity (by ∼ 75 % in maximum vertical velocity) and precipitation intensity (up to 45 %), with the anthropogenic aerosol effect more significant than the urban land effect. Urban land effect modifies convective evolution: speed up the transition from the warm cloud to mixed-phase cloud, thus initiating surface rain earlier but slowing down the convective cell dissipation, all of which result from urban heating-induced stronger sea-breeze circulation. The anthropogenic aerosol effect becomes evident after the cloud evolves into the mixed-phase cloud, accelerating the development of storm from the mixed-phase cloud to deep cloud by ∼ 40 min. Through aerosol–cloud interaction (ACI), aerosols boost convective intensity and precipitation mainly by activating numerous ultrafine particles at the mixed-phase and deep cloud stages. This work shows the importance of considering both the urban land and anthropogenic aerosol effects for understanding urbanization effects on convective clouds and precipitation.

54 ENVIRONMENTAL SCIENCES↗