Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Local mixing”

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 109 records · Page 6

Ducted fuel injection

Various technologies presented herein relate to heating a fuel and enhancing fuel and charge-gas mixing inside a combustion chamber to enable minimal, or no, generation of soot and/or other undesired emissions during ignition and subsequent combustion of the locally premixed mixtures. To enable sufficient mixing of the fuel and charge-gas, a jet of fuel can be directed to pass through a bore of a duct causing charge-gas to be drawn into the bore creating turbulence to mix the fuel and the drawn charge-gas. The duct can be heated to provision heating of the fuel and charge-gas mixture. The duct can be located proximate to an opening in a tip of a fuel injector. An ignition assist component can be located downstream of the duct to facilitate ignition of the fuel/charge-gas mixture.

42 ENGINEERING↗

Pressure-Induced Metal-like Transport and Magnetoresistance in a Au 2+ –Au 3+ Halide Perovskite

The Cs 4 Au II Au III 2 Cl 12 perovskite (1), featuring AuCl 4 trimers separated by vacancies, enables the first high-pressure study of Au 2+/3+ mixed-valence. Our computational analysis of the gold frontier orbitals suggests that the Au 2+ →Au 3+ intervalence charge transfer (IVCT) occurs across the vacancies. Computational structures indicate that these vacancies rapidly shrink with pressure and the Au 2+ and Au 3+ coordination spheres become very similar at the phase transition to nearly cubic symmetry at ca. 15 GPa─enabling facile IVCT. Although the activation energy of conductivity of 0.73(4) meV and far-infrared absorption indicate a small but nonzero bandgap, ambient thermal energy drives the IVCT, affording metallic properties: prominent infrared reflectivity and transport values of 10 2 S·cm –1 . This prompted us to perform the first high-pressure studies of magnetoresistance (MR) and Hall effect in halide perovskites. At 16 GPa, the MR increases by 9.3% at 2 K and 9 T; this value is maintained up to 27 GPa, when a local distortion drives electronic localization. By globally fitting the MR and Hall resistance to a two-carrier model we quantify how the carrier densities and mobilities evolve with pressure. Thus, metal-like transport and MR in 1 is driven by a pressure-induced transition from localized to partially delocalized mixed-valence.

Deschene, Christina R. [Stanford University, CA (U↗

Investigation of high-pressure turbulent jets using direct numerical simulation

Direct numerical simulations of free round jets at a Reynolds number (Re D ) of 5000, based on jet diameter (D) and jet-exit bulk velocity (U e ), are performed to study jet turbulence characteristics at supercritical pressures. The jet consists of nitrogen (N 2 ) that is injected into N 2 at the same temperature. To understand turbulent mixing, a passive scalar is transported with the flow at unity Schmidt number. Two sets of inflow conditions that model jets issuing from either a smooth contraction nozzle (laminar inflow) or a long pipe nozzle (turbulent inflow) are considered. By changing one parameter at a time, the simulations examine the jet-flow sensitivity to the thermodynamic condition (characterized in terms of the compressibility factor ($Z$) and the normalized isothermal compressibility), inflow condition and ambient pressure (p ∞ ) spanning perfect- to real-gas conditions. The inflow affects flow statistics in the near field (containing the potential core closure and the transition region) as well as further downstream (containing fully developed flow with self-similar statistics) at both atmospheric and supercritical p ∞ . The sensitivity to inflow is larger in the transition region, where the laminar-inflow jets exhibit dominant coherent structures that produce higher mean strain rates and higher turbulent kinetic energy than in turbulent-inflow jets. Decreasing $Z$ at a fixed supercritical p ∞ enhances pressure and density fluctuations (non-dimensionalized by local mean pressure and density, respectively), but the effect on velocity fluctuations depends also on the local flow dynamics. When $Z$ is reduced, large mean strain rates in the transition region of laminar-inflow jets significantly enhance velocity fluctuations (non-dimensionalized by local mean velocity) and scalar mixing, whereas the effects are minimal in jets from turbulent inflow.

42 ENGINEERING↗

Direct Observations of Solute Dispersion in Rocks With Distinct Degree of Sub‐Micron Porosity

Abstract The transport of chemical species in rocks is affected by their structural heterogeneity to yield a wide spectrum of local solute concentrations. To quantify such imperfect mixing, advanced methodologies are needed that augment the traditional breakthrough curve analysis by probing solute concentration within the fluids locally. Here, we demonstrate the application of asynchronous, multimodality imaging by X‐ray computed tomography (XCT) and positron emission tomography (PET) to the study of passive tracer experiments in laboratory rock cores. The four‐dimensional concentration maps measured by PET reveal specific signatures of the transport process, which we have quantified using fundamental measures of mixing and spreading. We observe that the extent of solute spreading correlate strongly with the strength of subcore‐scale porosity heterogeneity measured by XCT, while dilution is enhanced in rocks containing substantial sub‐micron porosity. We observe that the analysis of different metrics is necessary, as they can differ in their sensitivity to the strength and forms of heterogeneity. The multimodality imaging approach is uniquely suited to probe the fundamental difference between spreading and mixing in heterogeneous media. We propose that when multi‐dimensional data is available, mixing and spreading can be independently quantified using the same metric. We also demonstrate that one‐dimensional transport models have limited predictive ability toward the internal evolution of the solute concentration, when the model is solely calibrated against the effluent breakthrough curves. The data set generated in this study can be used to build realistic digital rock models and to benchmark transport simulations that account deterministically for rock property heterogeneity.

Kurotori, Takeshi [Department of Chemical Engineer↗

Cosolvent-tuned interactions in ionic liquids: A vibrational and quantum-chemical study of ethylene glycol ratio effects

Ionic liquids (ILs) are attractive media for CO 2 capture but remain limited by viscosity and cost. Blending ILs with ethylene glycol (EG) is a practical route to mitigate these constraints, yet the molecular origins of cosolvent effects and their dependence on composition are not well resolved. We combine Fourier-transform infrared (FT-IR) spectroscopy with quantum-chemical (DFT) analysis to elucidate how the IL:EG molar ratio modulates intermolecular interactions and electronic structure. Computed vibrational frequencies enable mode assignment and deconvolution of overlapping bands, revealing systematic, ratio-dependent shifts and broadenings in (i) EG O–H stretching, (ii) cation and EG C–H stretchings (imidazolium C2–H, C4–H, C5–H, methyl and ethyl groups, -CH2 of EG), (iii) anion signature modes (e.g., CN motifs), and (iv) EG C–O and C–C stretchings, consistent with the redistribution of hydrogen-bonding networks. Molecular electrostatic potential (MESP) maps quantify attenuation of extreme potential regions with increasing EG, indicating progressive screening of cation–anion electrostatic interactions. Quantum Theory of Atoms in Molecules (QTAIM) identifies emergent bond critical points between EG and the IL ions, while Reduced Density Gradient–Noncovalent Interaction (RDG–NCI) analysis differentiates strong directional hydrogen bonds from dispersive contacts across compositions. Together, these results show that EG fraction controls a switch from predominantly ion–ion to mixed ion–EG coordination, altering local polarity and polarizability that underlie the observed FT-IR trends. The framework provides composition–structure–spectrum relationships that can guide rational selection of IL:EG ratios to balance favorable molecular interactions with practical performance targets in scalable CO 2 capture systems.

DAC↗

Using molecular dynamics to predict the solidus and liquidus of mixed oxides ( Th , U ) O 2 , ( Th , Pu ) O 2 and ( Pu , U ) O 2

Molecular dynamics (MD) was used to establish a mechanistic basis for the experimentally observed reduction in liquidus and solidus temperatures below the melting point of the end-members for the mixed oxides (Th,U)O 2 , (Th, Pu) O 2 and (Pu, U) O 2 . This dip is found at additions of the oxide with higher melting point to the oxide with the lower melting point. There are many causes suggested for the dip; here the distribution of the cation Frenkel energy for the mixed oxides caused by the local environment is proposed as a contributor. Furthermore, a variant of the moving interface method which yields information on the position of the solidus and liquidus boundaries, is used to predict the phase diagrams of these systems.

36 MATERIALS SCIENCE↗

Experimental study of vorticity-strain interactions in turbulent premixed counterflow flames

Here, the effects of heat release on interactions between vorticity (ω) and strain rate (s) in turbulent premixed CH 4 /O 2 /N 2 counterflow flames are investigated using simultaneous OH laser-induced fluorescence (LIF) and tomographic particle image velocimetry (TPIV) measurements. A comparison between the flames and a corresponding turbulent non-reacting variable density N 2 -vs-products counterflow reveals the impact of heat release on vorticity-strain rate alignment statistics. Vorticity and strain rate statistics in the flames and non-reacting flow are conditioned on distance from the local flame front and gas mixing layer interface (GMLI) contours, respectively. The magnitude, alignment, and spatial distribution of the vorticity and principal strain rates (s 1 , s 2 , s 3 ) are rather different when heat release is present. Density variations without heat release enhance the ω-s 2 alignment while significantly reducing the ω-s 3 alignment and modestly reducing the ω-s 1 alignment. In contrast, heat release at the flame front further reduces the ω-s 1 alignment but increases the ω-s 3 alignment and suppresses the preferential ω-s 2 alignment. Furthermore, increasing turbulence diminishes the effect of heat release on this preferential alignment. In regions with the largest vorticities, both the reacting and non-reacting counterflows show an increase in the probability of ω-s 2 alignment. All counterflow cases have a net positive vortex-stretching contribution to the enstrophy production with a peak production rate at the flame front or GMLI, but the peak values depend on the density variation, heat release, and turbulence level. Elucidation of the complex interplay between these factors contributes to the understanding of the dynamics of turbulence-flame interactions.

42 ENGINEERING↗

Perturbative analysis of low-frequency instabilities in high-field ST40 experiments

Abstract Linear and nonlinear perturbative analyses of low-frequency instabilities are presented for ST40 spherical torus high field plasmas, with simultaneous neutral beam heating due to two tangential sources at 25 keV and 55 keV. Initial plasma profiles on which the TRANSP and NUBEAM codes were based were obtained using an integrated analysis of several diagnostics including line-of-sight integrated and volume average measurements, as well as limited profile information from a charge-exchange-recombination spectrometer. Magnetohydrodynamics (MHD) analysis using the NOVA/NOVA-K codes suggests that the modes in the experiments are core-localized n = 1 eigenmodes with mixed acoustic and electromagnetic Alfvénic polarization. The global transport combined with the perturbative mode analysis indicates that there is greater tendency for bursty chirping response as the level of background micro-turbulence decreases in time and those predictions are shown to be consistent with the experimental data. This finding is interpreted in terms of the suppression of coherent phase-space structures that support chirping when the resonant fast ion dynamics is dominated by frequent stochastic, orbit-decorrelating events.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Techno-economic optimization of renewable urea production for sustainable agriculture and CO 2 utilization

Abstract Urea is the most used nitrogen fertilizer due to its ease of storage, transportation, and application. It is made by combining ammonia and carbon dioxide (CO 2 ), both of which are produced predominantly from fossil fuels at present. The recent momentum behind ammonia production using renewable-powered electrolysis offers an opportunity to both make urea in a more sustainable way and utilize CO 2 from external sources. In this work, we present a techno-economic optimization model to minimize the cost of making urea in this way. The model allows for time-varying chemical production in response to renewable variability by simultaneously optimizing production facility design and hourly operation. We performed a case study for Minnesota considering the use of byproduct CO 2 from bioethanol production. We found that the present-day levelized cost of renewable urea is between $268 mt −1 and $413 mt −1 at likely implementable production scales up to 250 000 mt yr −1 . This is within the range of historical conventional urea prices while offering at least 78% carbon intensity reduction. Projecting to 2030, there is a clear economic case for renewable urea production with levelized cost as low as $135 mt −1 due to technology improvement and electrolysis manufacturing expansion, facilitating a urea production scale increase to 525 000 mt yr −1 . Optimal facilities use wind energy, with hydrogen and ammonia production operating in a flexible, time-varying way to minimize battery and hydrogen storage capacities. Urea production operates near steady state due to the relatively low cost of intermediate ammonia buffer storage. A mix of imported methane and locally produced hydrogen are used to provide heat for steam consumed in the urea synthesis.

Chemistry↗

Interpolating between small- and large-$\mathscr{g}$ expansions using Bayesian model mixing

Bayesian model mixing (BMM) is a statistical technique that can be used to combine models that are predictive in different input domains into a composite distribution that has improved predictive power over the entire input space. We explore the application of BMM to the mixing of two expansions of a function of a coupling constant $\mathscr{g}$ that are valid at small and large values of $\mathscr{g}$ respectively. This type of problem is quite common in nuclear physics, where physical properties are straightforwardly calculable in strong and weak interaction limits or at low and high densities or momentum transfers, but difficult to calculate in between. Interpolation between these limits is often accomplished by a suitable interpolating function, e.g., Padé approximants, but it is then unclear how to quantify the uncertainty of the interpolant. We address this problem in the simple context of the partition function of zero-dimensional Φ 4 theory, for which the (asymptotic) expansion at small $\mathscr{g}$ and the (convergent) expansion at large $\mathscr{g}$ are both known. We consider three mixing methods: linear mixture BMM, localized bivariate BMM, and localized multivariate BMM with Gaussian processes. Here, we find that employing a Gaussian process in the intermediate region between the two predictive models leads to the best results of the three methods. The methods and validation strategies we present here should be generalizable to other nuclear physics settings.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

RF mixing modules for superconducting qubit room temperature control systems v1.0

We have developed a series of compact RF mixing boards by integrating I/Q quadrature mixing, IF(intermediate frequency)/LO(local oscillator)/RF power level adjustments, and DC (direct current) bias fine tuning on a 40 mm × 80 mm 4-layer PCB (printed circuit board) board with EMI (electromagnetic interference) shielding.

Huang, Gang↗

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗

Photoinduced Anion Segregation in Mixed Halide Perovskites

Alloyed lead halide perovskites have taken a dominant role in the quest for third-generation solar cells. This is due to optimal light-harvesting properties, which can be tuned across the visible spectrum by mixing halide (X = Cl – , Br – , and I – ) anions and A + cations (A + = FA + , MA + , and Cs + ). Durability issues related to ion movement within the perovskite lattice, however, impede large-scale commercialization. Uniformly mixed halide perovskites [e.g., APb(I 1– x Br x ) 3 ] reversibly segregate into narrow bandgap I-rich and wide bandgap Br-rich domains during continuous visible illumination. Subsequent I-rich domains reduce local open circuit voltages and decrease mixed halide perovskite solar cell power conversion efficiencies. In this review, we assess the known effects of halide segregation on the structural and optical properties of mixed halide materials, discuss ongoing research to suppress the phenomenon, and finally provide a mechanistic overview of its underlying origins.

36 MATERIALS SCIENCE↗

Standoff Infrared Measurements of Chemical Plume Dynamics in Complex Terrain Using a Combination of Active Swept-ECQCL Laser Spectroscopy with Passive Hyperspectral Imaging

Chemical plume detection and modeling in complex terrain present numerous challenges. We present experimental results from outdoor releases of two chemical tracers (sulfur hexafluoride and Freon-152a) from different locations in mountainous terrain. Chemical plumes were detected using two standoff instruments collocated at a distance of 1.5 km from the plume releases. A passive long-wave infrared hyperspectral imaging system was used to show time- and space-resolved plume transport in regions near the source. An active infrared swept-wavelength external cavity quantum cascade laser system was used in a standoff configuration to measure quantitative chemical column densities with high time resolution and high sensitivity along a single measurement path. Both instruments provided chemical-specific detection of the plumes and provided complementary information over different temporal and spatial scales. The results show highly variable plume propagation dynamics near the release points, strongly dependent on the local topography and winds. Effects of plume stagnation, plume splitting, and plume mixing were all observed and are explained based on local topographic and wind conditions. Measured plume column densities at distances ~100 m from the release point show temporal fluctuations over ~1 s time scales and spatial variations over ~1 m length scales. The results highlight the need for high-speed and spatially resolved measurement techniques to provide validation data at the relevant spatial and temporal scales required for high-fidelity terrain-aware microscale plume propagation models.

47 OTHER INSTRUMENTATION↗

Influence of mixing on the generation of nanoparticles in cement systems

In this work, we study the effect of mixing intensity on the size and concentration of nanoparticles in quartz and cement systems. We vary the local shear rate in the studied suspensions induced by mixing by either changing solid volume fraction or mixing speed. Our results obtained on quartz, cement and quartz-cement suspensions provide a master curve relating the concentration of nanoparticles in the pore solution as a function of the local shearing conditions computed at a micromechanical level. Our results suggest, in accordance with literature, that high mixing intensity is at the origin of a detachment of nano-hydrates from the surface of cement particles. These nanoparticles can be stabilized and therefore detected in the pore-solution only in the presence of superplasticizers, which inhibit their agglomeration to other particles.

36 MATERIALS SCIENCE↗

Local Power: Comparing County-Level Renewable Energy Potential to Consumption Using the SLOPE Platform

Wide-scale deployment of renewable energy technologies has the potential to significantly reduce greenhouse gas emissions and mitigate the effects of climate change. Many communities have ambitious clean energy goals with targets for locally generated renewable energy. To inform state and local clean energy planning, analysts from the Joint Institute for Strategic Energy Analysis (JISEA) and National Renewable Energy Laboratory (NREL) used data from NREL's State and Local Planning for Energy (SLOPE) platform to compare annual technical generation potential of renewable energy technologies to modeled electricity consumption in every county of the contiguous United States. Annual costs were calculated to produce a 20% share of electricity consumed annually from each technology to examine localized cost effectiveness of a diversified mix of generation sources. For example, combining distributed and utility-scale wind and solar generation can offset the need for storage and nonintermittent fossil resources to achieve high deployment of renewables. This county-level analysis provides insight into where localized renewable energy generation could cost-effectively match annual electricity consumption.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

First observations of edge instabilities in strongly shaped negative triangularity plasmas on DIII-D

Throughout the 2023 DIII-D negative triangularity (NT) campaign, unexpected magnetic fluctuations indicative of an edge localized instability were observed in the majority of strongly shaped (–0.55 ≤ δ ≤ –0.45) discharges. Moreover, these same fluctuations were found to be absent in weaker NT plasma shapes. Based on a database study, three separate mode behaviors are observed based primarily on the edge safety factor and normalized plasma beta β N . The main regime is the bursty state, in which short (~5 ms) magnetic bursts occur every 5–10 ms. Each burst drives a small reduction in the edge electron temperature (T e ) as measured by electron cyclotron emissions (ECE), which then rapidly recovers between the bursts. The quiescent regime appears to have governing mechanics similar to the bursty regime, but with multiple unique mode signatures observed in the magnetics. These additional modes greatly reduce the bursty nature of the plasma. Finally, a saturated state is observed in a number of discharges, particularly at early time slices. The saturated regime is often short-lived (maximum 500 ms), but can have a large impact on the dynamics of the plasma edge region, causing significant drops in electron temperature as observed in the ECE diagnostic. 2D ECE imaging (ECE-I) is utilized, showing large kink-like perturbations localized near the last closed flux surface that extend into the SOL indicative of an external kink/peeling behavior. Additionally, evidence of phase mixing and observed tearing-parity perturbations localized near the edge rational surfaces is observed in the quiescent regime plasmas. As a result, the various modes do not appear to have a large impact on global confinement but may play a role in regulating the transport of the NT edge and the observed edge gradients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Relative importance of high-latitude local and long-range-transported dust for Arctic ice-nucleating particles and impacts on Arctic mixed-phase clouds

Dust particles, serving as ice-nucleating particles (INPs), may impact the Arctic surface energy budget and regional climate by modulating the mixed-phase cloud properties and lifetime. In addition to long-range transport from low-latitude deserts, dust particles in the Arctic can originate from local sources. However, the importance of high-latitude dust (HLD) as a source of Arctic INPs (compared to low-latitude dust, LLD) and its effects on Arctic mixed-phase clouds are overlooked. In this study, we evaluate the contribution to Arctic dust loading and INP population from HLD and six LLD source regions by implementing a source-tagging technique for dust aerosols in version 1 of the US Department of Energy's Energy Exascale Earth System Model (E3SMv1). Our results show that HLD is responsible for 30.7 % of the total dust burden in the Arctic, whereas LLD from Asia and North Africa contributes 44.2 % and 24.2 %, respectively. Due to its limited vertical transport as a result of stable boundary layers, HLD contributes more in the lower troposphere, especially in boreal summer and autumn when the HLD emissions are stronger. LLD from North Africa and East Asia dominates the dust loading in the upper troposphere with peak contributions in boreal spring and winter. The modeled INP concentrations show better agreement with both ground and aircraft INP measurements in the Arctic when including HLD INPs. The HLD INPs are found to induce a net cooling effect (-0.24 W m -2 above 60 °N) on the Arctic surface downwelling radiative flux by changing the cloud phase of the Arctic mixed-phase clouds. The magnitude of this cooling is larger than that induced by North African and East Asian dust (0.08 and -0.06 W m -2 , respectively), mainly due to different seasonalities of HLD and LLD. Uncertainties of this study are discussed, which highlights the importance of further constraining the HLD emissions.

54 ENVIRONMENTAL SCIENCES↗