Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spectral shape”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Experimental investigation of turbulent wake flows in a helically wrapped rod bundle in presence of localized blockages

In nuclear sodium fast reactors, bundles of rods are tightly packed into a triangular lattice, enclosed in a hexagonal duct, and each pin is spirally wrapped with a thin wire. Flow blockages can potentially impact the local flow characteristics and heat transfer mechanisms in the bundle due to its small subchannel size. The effects of the blockage on the flow structures and heat transfer mechanisms are important aspects that require an accurate investigation. In this study, the flow-field characteristics in the vicinity of a blockage located in the exterior subchannel of rod bundles with helically wrapped wires were experimentally investigated. The velocity fields in the exterior subchannel were acquired by applying matched-index-of-refraction and time-resolved particle image velocimetry (TR-PIV) techniques for Reynolds numbers of Re 1 = 4000 and Re 2 = 17 000, i.e., equivalent to Rew 1 = 19 600 and Rew 2 = 83 200, respectively, based on the blockage width. The results from the TR-PIV measurements revealed an arch-shaped vortex with a large flow recirculation and a pair of counter-rotating vortices in the wake region downstream of the blockage, which is commonly observed in the wake flow of bluff bodies. The relative lateral distance and angle between the two vortices decreased when the Reynolds numbers increased. Profiles of maximum turbulence intensity along the shear layers illustrated the transition process including the growth, peak, and decay along the flow direction. From the spectral analysis of the turbulent velocities extracted at points along the shear layer, the Strouhal numbers ( St ) representing the vortex shedding frequency were found to be St = 0.25 and St = 0.56 for the left and right shear layers, respectively. Characteristics of shear layers generated by the blockage in the exterior subchannel were investigated via the two-point cross correlation of fluctuating velocities. The patiotemporal cross correlations of turbulent velocities, computed at points in the region where the left shear layer exhibited rolling effects and vortex breakdowns, were considerably wider and longer. The convection velocity U c was estimated to be ~0.82 U m to 0.93 U m . Proper orthogonal decomposition (POD) analysis was applied to the instantaneous velocity fields to extract the statistically dominant flow structures. It was found that POD modes 2–3 and 4–5 formed the pair modes when the corresponding POD temporal coefficients depicted sinusoidal shapes and exhibited nearly circular orbits in the phase space. Spectral analysis of the POD temporal coefficients confirmed the vortex shedding frequencies detected in the analysis of turbulent velocities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Assessment of four strain energy decomposition methods for phase field fracture models using quasi-static and dynamic benchmark cases

Abstract Strain energy decomposition methods in phase field fracture models separate strain energy that contributes to fracture from that which does not. However, various decomposition methods have been proposed in the literature, and it can be difficult to determine an appropriate method for a given problem. The goal of this work is to facilitate the choice of strain decomposition method by assessing the performance of three existing methods (spectral decomposition of the stress or the strain and deviatoric decomposition of the strain) and one new method (deviatoric decomposition of the stress) with several benchmark problems. In each benchmark problem, we compare the performance of the four methods using both qualitative and quantitative metrics. In the first benchmark, we compare the predicted mechanical behavior of cracked material. We then use four quasi-static benchmark cases: a single edge notched tension test, a single edge notched shear test, a three-point bending test, and a L-shaped panel test. Finally, we use two dynamic benchmark cases: a dynamic tensile fracture test and a dynamic shear fracture test. All four methods perform well in tension, the two spectral methods perform better in compression and with mixed mode (though the stress spectral method performs the best), and all the methods show minor issues in at least one of the shear cases. In general, whether the strain or the stress is decomposed does not have a significant impact on the predicted behavior.

Zhang, Shuaifang↗
DOE OSTI2022

Shallow Cumulus Properties as Captured by Adiabatic Fraction in High-Resolution LES Simulations

Abstract Shallow convective clouds are important players in Earth’s energy budget and hydrological cycle, and are abundant in the tropical and subtropical belts. They greatly contribute to the uncertainty in climate predictions due to their unresolved, complex processes that include coupling between the dynamics and microphysics. Analysis of cloud structure can be simplified by considering cloud motions as a combination of moist adiabatic motions like adiabatic updrafts and turbulent motions leading to deviation from adiabaticity. In this work, we study the sizes and occurrence of adiabatic regions in shallow cumulus clouds during their growth and mature stages, and use the adiabatic fraction (AF) as a continuous metric to describe cloud processes and properties from the core to the edge. To do so, we simulate isolated trade wind cumulus clouds of different sizes using the System of Atmospheric Modeling (SAM) model in high resolution (10 m) with the Hebrew University spectral bin microphysics (SBM). The fine features in the clouds’ dynamics and microphysics, including small near-adiabatic volumes and a thin transition zone at the edge of the cloud (∼20–40 m in width), are captured. The AF is shown to be an efficient measure for analyzing cloud properties and key processes determining the droplet-size distribution formation and shape during the cloud evolution. Physical processes governing the properties of droplet size distributions at different cloud regions (e.g., core, edge) are analyzed in relation to AF. Significance Statement 1) This study investigates the evolution of cumulus clouds (Cu) using a 10-m-resolution LES model with spectral bin microphysics. 2) The study improves the understanding of the mutual effects of adiabatic updrafts and lateral entrainment and mixing. 3) The study demonstrates the existence of an adiabatic core in nonprecipitating Cu. 4) Shapes of the droplet size distributions are closely related to the adiabatic fraction values. 5) Utilization of high resolution reveals the existence of physically significant small features in the cloud structure, such as a narrow cloud interface zone and small adiabatic volumes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Machine learning of superconducting critical temperature from Eliashberg theory

Abstract The Eliashberg theory of superconductivity accounts for the fundamental physics of conventional superconductors, including the retardation of the interaction and the Coulomb pseudopotential, to predict the critical temperature T c . McMillan, Allen, and Dynes derived approximate closed-form expressions for the critical temperature within this theory, which depends on the electron–phonon spectral function α 2 F ( ω ). Here we show that modern machine-learning techniques can substantially improve these formulae, accounting for more general shapes of the α 2 F function. Using symbolic regression and the SISSO framework, together with a database of artificially generated α 2 F functions and numerical solutions of the Eliashberg equations, we derive a formula for T c that performs as well as Allen–Dynes for low- T c superconductors and substantially better for higher- T c ones. This corrects the systematic underestimation of T c while reproducing the physical constraints originally outlined by Allen and Dynes. This equation should replace the Allen–Dynes formula for the prediction of higher-temperature superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Nanomechanical vibrational response from electrical mixing measurements

Driven nanomechanical resonators based on low-dimensional materials are routinely and efficiently detected with electrical mixing measurements. However, the measured signal is a non-trivial combination of the mechanical eigenmode displacement and an electrical contribution, which makes the extraction of the driven mechanical response challenging. Here, we report a simple yet reliable method to extract solely the driven mechanical vibrations by eliminating the contribution of pure electrical origin. This enables us to measure the spectral mechanical response as well as the driven quadratures of motion. This method is crucial for nanomechanical vibrations in the nonlinear regime, since the shape of the mechanical response depends on the physics at work. We further show how to calibrate the measured signal into units of displacement. Our method marks a key step forward in the study of nanoelectromechanical resonators based on low-dimensional materials in the nonlinear regime.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dynamical Correlations and Order in Magic-Angle Twisted Bilayer Graphene

The interplay of dynamical correlations and electronic ordering is pivotal in shaping phase diagrams of correlated quantum materials. In magic-angle twisted bilayer graphene, transport, thermodynamic, and spectroscopic experiments pinpoint a competition between distinct low-energy states with and without electronic order, as well as between localized and delocalized charge carriers. In this study, we utilize dynamical mean-field theory on the topological heavy fermion model of twisted bilayer graphene to investigate the emergence of electronic correlations and long-range order in the absence of strain. We contrast moment formation, Kondo screening, and ordering on a temperature basis and explain the nature of emergent correlated states based on three central phenomena: (i) the formation of local spin and valley isospin moments around 100 K, (ii) the ordering of the local isospin moments around 10 K preempting Kondo screening, and (iii) a cascadic redistribution of charge between localized and delocalized electronic states upon doping. At integer fillings, we find that low-energy spectral weight is depleted in the symmetric phase, while we find insulating states with gaps enhanced by exchange coupling in the zero-strain ordered phases. Doping away from integer filling results in distinct metallic states: a “bad metal” above the ordering temperature, where scattering off the disordered local moments suppresses electronic coherence, and a “good metal” in the ordered states with coherence of quasiparticles facilitated by isospin order. This finding reveals coherence from order as the microscopic mechanism behind the Pomeranchuk effect observed experimentally by Rozen [] and by Saito []. Upon doping, there is a periodic charge reshuffling between localized and delocalized electronic orbitals leading to cascades of doping-induced Lifshitz transitions, local spectral weight redistributions, and periodic variations of the electronic compressibility ranging from nearly incompressible to negative. Our findings highlight the essential role of charge transfer, hybridization, and ordering in shaping the electronic excitations and thermodynamic properties in twisted bilayer graphene and provide a unified understanding of the most puzzling aspects of scanning tunneling spectroscopy, transport, and compressibility experiments. Published by the American Physical Society 2024

Physics↗
DOE OSTI2021

Grain size dependencies of intergranular solute segregation in nanocrystalline materials

Grain boundary segregation is of increasing importance for stabilizing nanocrystalline alloys, and while segregation is typically modeled as a scale-independent phenomenon, at the finest nanocrystalline grain sizes there are several sources of grain size dependence. Here we provide quantitative detail on two grain size dependencies not previously elaborated, both of which arise from the spectrality of intergranular segregation sites. The first of these arises because a change in grain size affects which portions of the segregation spectrum are occupied by the solutes, with no change in the shape or nature of the spectrum itself. The second of these arises from the presence of higher-order grain junctions, which have unique segregation behaviors and an increasing relative presence at finer grain sizes. This paper uses molecular statics to quantify these size dependencies in detail for the case of an Al(Mg) binary alloy. Together, these two effects combined have a significant effect on the effective McLean segregation energy of more than 6 kJ/mol (from approximately -18.5 kJ/mol at 5 nm to -12.5 kJ/mol at more than 40 nm), with the majority of grain size dependence arising from the spectrality of segregation sites.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantum control of exciton wave functions in 2D semiconductors

Excitons—bound electron-hole pairs—play a central role in light-matter interaction phenomena and are crucial for wide-ranging applications from light harvesting and generation to quantum information processing. A long-standing challenge in solid-state optics has been to achieve precise and scalable control over excitonic motion. We present a technique using nanostructured gate electrodes to create tailored potential landscapes for excitons in 2D semiconductors, enabling in situ wave function shaping at the nanoscale. Our approach forms electrostatic traps for excitons in various geometries, such as quantum dots, rings, and arrays thereof. We show independent spectral tuning of spatially separated quantum dots, achieving degeneracy despite material disorder. Owing to the strong light-matter coupling of excitons in 2D semiconductors, we observe unambiguous signatures of confined exciton wave functions in optical reflection and photoluminescence measurements. This work unlocks possibilities for engineering exciton dynamics and interactions at the nanometer scale, with implications for optoelectronic devices, topological photonics, and quantum nonlinear optics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Efficient Low-Order Refined Preconditioners for High-Order Matrix-Free Continuous and Discontinuous Galerkin Methods

In this paper, we design preconditioners for the matrix-free solution of high-order continuous and discontinuous Galerkin discretizations of elliptic problems based on finite element method--spectral element method (FEM-SEM) equivalence and additive Schwarz methods. The high-order operators are applied without forming the system matrix, making use of sum factorization for efficient evaluation. The system is preconditioned using a spectrally equivalent low-order (p = 1) finite element operator discretization on a refined mesh. The low-order refined mesh is anisotropic and not shape regular in the polynomial degree of the high-order operator, requiring specialized solvers to treat the anisotropy. We make use of an element-structured, geometric multigrid V-cycle with ordered ILU(0) smoothing. The preconditioner is parallelized through an overlapping additive Schwarz method that is robust in h and p. The method is extended to interior penalty and Bassi and Rebay (BR2) discontinuous Galerkin discretizations, for which it is also robust in the size of the penalty parameter. Finally, numerical results are presented on a variety of examples, verifying the uniformity of the preconditioner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

GRB 110213A: A Study of Afterglow Electromagnetic Cascade Radiation

We obtained well-sampled optical photometry of GRB 110213A, including Swift/UVOT and XRT. Combining our data from those of other ground-based telescopes, we present 15 optical multicolor light curves showing similar shapes with two peaks. In contrast, in the X-ray band, only a single peak is observed between the two optical peaks. Temporal and spectral analysis of GRB 110213A shows that the X-rays differ from the optical for Phases I–III (before the second peak of the optical band at ~5.6 ks). Moreover, they have the same spectral behavior at late times (Phases IV–VI). These data indicate that the optical and X-ray emission are dominated by different components. The synchrotron-supported pair cascade emission is included in the standard external forward-shock model, which is dominated by synchrotron radiation and synchrotron self-Compton (SSC). We find that the optical bands of GRB 110213A are dominated by the cascade emission from synchrotron radiation + SSC at the early stage, while the primary synchrotron + SSC radiation dominates the X-ray band. At late stages, both the X-ray and optical bands are dominated by emission from primary synchrotron + SSC radiation. The cascade component can reasonably explain the first optical peak. In contrast, the primary synchrotron + SSC emission mainly contributes to the second peak.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Third Fermi Large Area Telescope Catalog of Gamma-Ray Pulsars

We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems colocated with LAT sources also likely harbor gamma-ray MSPs. This catalog thus reports roughly 340 gamma-ray pulsars and candidates, 10% of all known pulsars, compared to ≤11 known before Fermi. Half of the gamma-ray pulsars are young. Of these, the half that are undetected in radio have a broader Galactic latitude distribution than the young radio-loud pulsars. The others are MSPs, with six undetected in radio. Overall, ≥236 are bright enough above 50 MeV to fit the pulse profile, the energy spectrum, or both. For the common two-peaked profiles, the gamma-ray peak closest to the magnetic pole crossing generally has a softer spectrum. The spectral energy distributions tend to narrow as the spindown power $\dot{E}$ decreases to its observed minimum near 10 33 erg s -1 , approaching the shape for synchrotron radiation from monoenergetic electrons. We calculate gamma-ray luminosities when distances are available. Our all-sky gamma-ray sensitivity map is useful for population syntheses. The electronic catalog version provides gamma-ray pulsar ephemerides, properties, and fit results to guide and be compared with modeling results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Gravitational effective field theory islands, low-spin dominance, and the four-graviton amplitude

Abstract We analyze constraints from perturbative unitarity and crossing on the leading contributions of higher-dimension operators to the four-graviton amplitude in four spacetime dimensions, including constraints that follow from distinct helicity configurations. We focus on the leading-order effect due to exchange by massive degrees of freedom which makes the amplitudes of interest infrared finite. In particular, we place a bound on the coefficient of the R 3 operator that corrects the graviton three-point amplitude in terms of the R 4 coefficient. To test the constraints we obtain nontrivial effective field-theory data by computing and taking the large-mass expansion of the one-loop minimally-coupled four-graviton amplitude with massive particles up to spin 2 circulating in the loop. Remarkably, we observe that the leading EFT coefficients obtained from both string and one-loop field-theory amplitudes lie in small islands. The shape and location of the islands can be derived from the dispersive representation for the Wilson coefficients using crossing and assuming that the lowest-spin spectral densities are the largest. Our analysis suggests that the Wilson coefficients of weakly-coupled gravitational physical theories are much more constrained than indicated by bounds arising from dispersive considerations of 2 → 2 scattering. The one-loop four-graviton amplitudes used to obtain the EFT data are computed using modern amplitude methods, including generalized unitarity, supersymmetric decompositions and the double copy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Resonant propagation of extreme-ultraviolet pulses through strongly driven high-density media

We show that by combining strong-field dressing and resonant propagation of XUV pulses, the transition of absorption lines from their natural Lorentzian profiles through Fano and complex multipeak shapes all the way back to broadened near-Lorentzian profiles can be achieved, in the limit of optically thick samples. The final stage of this spectral modification can be understood in terms of a significant temporal stretching and delay of the resonant XUV pulse as it propagates through the dense gas, which alters the ultrafast absorption that is modified by the time-synchronized few-femtosecond laser pulse. We first demonstrate this concept in numerical calculations, both using a model system and through a fully coupled solution of the time-dependent Schrödinger equation and the Maxwell wave equation. The applicability and generality of the underlying mechanism is then illustrated in proof-of-principle attosecond transient absorption measurements in a helium gas with a widely varying atomic density. These results provide insights into the interaction of ultrashort laser pulses with dense media and its coherent control. Published by the American Physical Society 2025

He, Yu (ORCID:0000000152120176)↗
DOE OSTI2024

Assessment of Machine Learning Wall Modeling Approaches for Large Eddy Simulation of Gas Turbine Film Cooling Flows: An a Priori Study

Here, in this work, a priori analysis of machine learning (ML) strategies is carried out with the goal of data-driven wall modeling for large eddy simulation (LES) of gas turbine film cooling flows. High-fidelity flow datasets are extracted from wall-resolved LES (WRLES) of flow over a flat plate interacting with the coolant flow supplied by a single row of 7-7-7 shaped cooling holes inclined at 30 degrees with the flat plate at different blowing ratios (BR). The WRLES are performed using the high-order Nek5000 spectral element computational fluid dynamics (CFD) solver. Light gradient boosting machine (LightGBM) is employed as the ML algorithm for the data-driven wall model. Parametric tests are conducted to systematically assess the influence of a wide range of input flow features (velocity components, velocity gradients, pressure gradients, and fluid properties) on the accuracy of ML wall model with respect to prediction of wall shear stress. In addition, the use of spatial stencil and time delay is also explored within the ML wall modeling framework. It is shown that features associated with gradients of the streamwise and spanwise velocity components have a major impact on the prediction fidelity of wall model, while the effect of gradients of wall-normal velocity component is found to be negligible. Moreover, adding flow feature information from an x-y-z spatial stencil significantly improves the ML model accuracy and generalizability compared to just using local flow features from the matching location. Overall, highest prediction accuracy is achieved when both spatial stencil and time delay features are incorporated within the data-driven wall modeling paradigm.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Large-eddy simulation of non-vaporizing sprays using the spectral-element method

Predictive simulations of high-pressure sprays require accurate representation of the turbulent gaseous flow field generated by liquid jet. Typically, the accuracy that can be obtained with low-order numerical methods (e.g. finite volume, finite element) is limited by stability issues in fine grids and the order of convergence of the method. In this work, we resolve the turbulent flow field in an Eulerian manner using the high-order spectral element method, coupled with a Lagrangian parcels approach to model the atomizing liquid jet. Large eddy simulations of single-hole sprays under non-evaporative conditions were conducted and compared against experimental data from Margot et al. (2008) and Spray A data from the Engine Combustion Network. The sensitivity of liquid penetration and droplet sizes to different breakup model parameters was studied. The effect of different numerical parameters, such as polynomial order of the solution (grid resolution), on liquid penetration was also analyzed. The method achieved grid-independent results using p-refinement, achieving finer resolution (by a factor of x 1.7 - x 3.5) in the gas-phase solution than in state-of-the-art simulations using the finite-volume method. Here the results showed good agreement with experimental data, demonstrating the ability of the current method to accurately capture liquid penetration and the shape of the spray.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Affinity propagation clustering of full-field, high-spatial-dimensional measurements for robust output-only modal identification: A proof-of-concept study

Determination of the model order is a challenging problem in system identification, especially in output-only or operational modal identification where some modes are weakly excited. While existing methods such as the stabilization diagram method (spectral information) are effective, they do not scale to high-dimensional data, which is usually needed for high-fidelity characterization of structural dynamics and has been made available in the emerging full-field measurement techniques using optical methods such as photogrammetry and laser vibrometers. In this proof-of-concept study we present a new non-parametric, data-driven approach for robust output-only identification of high-spatial-dimensional modal parameters of basic structures by efficiently processing and interactively exploiting the full-field measurement (i.e., very dense spatial measurement points). Specifically, we first over-estimate the system model once, producing a pool of candidate modes associated with their modal frequencies and full-field, high-spatial-dimensional mode shapes. This is accomplished by a data-driven method termed affinity propagation clustering (APC), where the active clusters, which are the active modes in our formulations, emerge from the "message-passing" procedure and does not require a pre-determination of the cluster number (mode or model order). Next, rather than using the spectral information to distinguish the physical and spurious modes in the stabilization diagram method, we exploit and visualize the spatial, full-field mode shape associated with each candidate mode to do so. We conduct extensive experiments on basic structural models with comparisons to a few existing methods. The results indicate that the new method is computationally efficient for identifying high-spatial-dimensional modal parameters, and robust to identify weak modes by exploiting the full-field measurement. We also discuss its applicability and limitations for structures with complex geometry (shapes).

42 ENGINEERING↗
DOE OSTI2020

Development of a High-Rate Front-End ASIC for X-Ray Spectroscopy and Diffraction Applications

We developed a new front-end application-specific integrated circuit (ASIC) to upgrade the Maia X-ray microprobe. The ASIC instruments 32 configurable channels that perform either positive or negative charge amplification, pulse shaping, peak amplitude, and time extraction along with buffered analog storage. At a gain of 3.6 V/fC, 1-μs peaking time, and a temperature of 248 K, an electronic resolution of 13 and 10 e- rms was measured with and without a silicon drift detector (SDD) sensor, respectively. A spectral resolution of 170-eV full-width at half-maximum (FWHM) at 5.9 keV was obtained with an 55Fe source. The channel linearity was better than ± 1% with rate capabilities up to 40 kcps. The ASIC was fabricated in a commercial 250-nm process with a footprint of 6.3 mm × 3.9 mm and dissipates 167 mW of static power.

Vernon, Emerson↗
DOE OSTI2021

Uncertainty of SW Cloud Radiative Effect in Atmospheric Models Due to the Parameterization of Liquid Cloud Optical Properties

Clouds are largely responsible for the spread of climate models predictions. Here we focus on the uncertainties in cloud shortwave radiative effect due to the parameterization of liquid cloud single scattering properties (SSPs) from liquid water content (LWC) and droplet number concentration (N), named parameterization of cloud optical properties. Uncertainties arise from not accounting for the droplet size distribution (DSD)—which affects the estimation of the effective radius (r eff ) and modulates the r eff -dependency of the SSPs—and from averaging SSPs over wide spectral bands. To assess these uncertainties a series of r eff -dependent SSPs parameterizations corresponding to various DSDs and spectral averaging methods are derived and implemented in a radiative code. Combined with the DSD-dependent estimation of r eff they are used to compute the bulk radiative properties (reflectance, transmittance, absorptance) of various clouds (defined in terms of LWC and N), including a homogeneous cloud, more realistic case studies, and outputs of a climate model. The results show that the cloud radiative forcing can vary up to 20% depending on the assumed DSD. Likewise, differences up to 20% are obtained for heating rates. The estimation of r eff is the main source of uncertainty, while the SSPs parameterization contributes to around 20% of the total uncertainty. Spectral averaging is less an issue, except for atmospheric absorption. Overall, global shortwave cloud radiative effect can vary by 6 W m –2 depending on the assumed DSD shape, which is about 13% of the best observational estimate.

54 ENVIRONMENTAL SCIENCES↗