Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Condition reports”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Techno-Economic, Feasibility, and Life Cycle Analysis of Renewable Propane: 2025 Update

To clarify the current and future landscape for renewable propane (RP) production, this work evaluates the value proposition of recovering RP from existing and planned hydroprocessed esters and fatty acids (HEFA) biorefineries and surveys emerging technologies under development or deployment. HEFA biorefineries co-produce a propane-rich fuel gas stream, normally used to meet HEFA process heat requirements, from which propane can be recovered and sold to create an additional revenue stream alongside liquid transportation fuels such as renewable diesel (RD) and sustainable aviation fuel (SAF). This report updates and extends a 2022 analysis of RP recovery from HEFA facilities by escalating capital and operating costs to 2024 prices, incorporating recent policy developments (including the Section 45Z Clean Fuel Production Credit), evaluating RP recovery for both RD- and SAF-focused HEFA facilities at two scales (3,000 and 75,000 barrels per day of feedstock), and quantifying the impact of RP recovery on HEFA liquid-fuel carbon intensity (CI) and associated tax credits using the 45ZCF-GREET model. For a 3,000 BPD RD-focused HEFA facility, approximately 3.5 million gallons per year (MGPY) of RP can be recovered; in this base case, the estimated payback period is 18 months based on the total installed cost of the RP recovery equipment and 36 months based on the total capital investment for the entire RP recovery project. The payback period is slightly shorter for the analogous SAF-focused configuration (approximately 4.3 MGPY RP). Sensitivity analysis shows that CAPEX magnitude, RP recovery plant scale, and CI-driven tax credit valuations are the dominant determinants of project viability. RP recovery may increase the CI of HEFA liquid fuels, which can reduce liquid-fuel tax credits (a key revenue stream for the HEFA biorefinery) and lengthen payback periods. However, RP recovery generally remains economically favorable across a wide range of plausible scenarios and market conditions. The report also summarizes emerging pathways that could expand future RP supply.

09 BIOMASS FUELS↗
DOE OSTI2022

Highly Active and Durable Air Electrodes for Reversible Protonic Ceramic Electrochemical Cells Enabled by an Efficient Bifunctional Catalyst

The commercialization of reversible protonic ceramic electrochemical cells is hindered by the lack of highly active and durable air electrodes exposed to high concentration of steam under operating conditions. We report findings that dramatically enhance the electrocatalytic activity and stability of a conventional (La 0.6 Sr 0.4 ) 0.95 Co 0.2 Fe 0.8 O 3-δ (LSCF) air electrode by a multiphase catalyst coating composed of a conformal Pr 1-x Ba x CoO 3-δ thin film and exsolved BaCoO 3-δ nanoparticles, are reported. At 600 °C, the catalyst coating decreases the polarization resistance of the LSCF air electrode by a factor of 25 (from 1.09 to 0.043 Ω cm 2 ) in air and the degradation rate by two orders of magnitude (from 1.0 × 10 -2 to 1.8 × 10 -4 Ω cm2 h -1 in humidified air with 30 vol% H2O). Further, a single cell with the catalyst-coated LSCF air electrode at 600 °C demonstrates a high peak power density of 1.04 W cm -2 in the fuel cell mode and a high current density of 1.82 A cm -2 at 1.3 V in the electrolysis mode. The significantly enhanced performance of the LSCF air electrode is attributed mainly to the high rate of surface oxygen exchange, fast surface proton diffusion, and the rapid H 2 O and O 2 dissociation on the catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Facile Access to Organostibines via Selective Organic Superbase Catalyzed Antimony‐Carbon Protonolysis

Abstract The selective formation of antimony‐carbon bonds via organic superbase catalysis under metal‐ and salt‐free conditions is reported. This novel approach utilizes electron‐deficient stibine, Sb(C 6 F 5 ) 3 , to give upon base‐catalyzed reactions with weakly acidic aromatic and heteroaromatic hydrocarbons access to a range of new aromatic and heteroaromatic stibines, respectively, with loss of C 6 HF 5 . Also, the significantly less electron‐deficient stibines, Ph 2 SbC 6 F 5 and PhSb(C 6 F 5 ) 2 smoothly underwent base‐catalyzed exchange reactions with a range of terminal alkynes to generate the stibines of formulae PhSb(C≡CPh) 2 , and Ph 2 SbC≡CR [R=C 6 H 5 , C 6 H 4 ‐NO 2 , COOEt, CH 2 Cl, CH 2 NEt 2 , CH 2 OSiMe 3 , Sb(C 6 H 5 ) 2 ], respectively. These formal substitution reactions proceed with high selectivity as only the C 6 F 5 groups serve as a leaving group to be liberated as C 6 HF 5 upon formal proton transfer from the alkyne. Kinetic studies of the base‐catalyzed reaction of Ph 2 SbC 6 F 5 with phenyl acetylene to form Ph 2 SbC≡CPh and C 6 HF 5 suggested the empirical rate law to exhibit a first‐order dependence with respect to the base catalyst, alkyne and stibine. DFT calculations support a pathway proceeding via a concerted σ‐bond metathesis transition state, where the base catalyst activates the Sb‐C 6 F 5 bond sequence through secondary bond interactions.

Culvyhouse, Jacob↗
DOE OSTI2024

Facile Access to Organostibines via Selective Organic Superbase Catalyzed Antimony‐Carbon Protonolysis

The selective formation of antimony-carbon bonds via organic superbase catalysis under metal- and salt-free conditions is reported. Here, this novel approach utilizes electron-deficient stibine, Sb(C 6 F 5 ) 3 , to give upon base-catalyzed reactions with weakly acidic aromatic and heteroaromatic hydrocarbons access to a range of new aromatic and heteroaromatic stibines, respectively, with loss of C 6 HF 5 . Also, the significantly less electron-deficient stibines, Ph 2 SbC 6 F 5 and PhSb(C 6 F 5 ) 2 smoothly underwent base-catalyzed exchange reactions with a range of terminal alkynes to generate the stibines of formulae PhSb(C≡CPh) 2 , and Ph 2 SbC≡CR [R=C 6 H 5 , C 6 H 4 -NO 2 , COOEt, CH 2 Cl, CH 2 NEt 2 , CH 2 OSiMe 3 , Sb(C 6 H 5 ) 2 ], respectively. These formal substitution reactions proceed with high selectivity as only the C 6 F 5 groups serve as a leaving group to be liberated as C 6 HF 5 upon formal proton transfer from the alkyne. Kinetic studies of the base-catalyzed reaction of Ph 2 SbC 6 F 5 with phenyl acetylene to form Ph 2 SbC≡CPh and C 6 HF 5 suggested the empirical rate law to exhibit a first-order dependence with respect to the base catalyst, alkyne and stibine. DFT calculations support a pathway proceeding via a concerted σ-bond metathesis transition state, where the base catalyst activates the Sb-C 6 F 5 bond sequence through secondary bond interactions.

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

Autonomous monitoring of algal biomass: Success stories and lessons learned from long-term field deployment

Autonomous, high-frequency monitoring of outdoor algal ponds is needed to quantify biomass productivity and detect culture decline in environments prone to contamination, grazers, and variable operating conditions. We report successes and lessons learned in translating a laboratory spectroradiometric monitoring approach to a multi-year autonomous field deployment at the Arizona Center for Algae Technology and Innovation (AzCATI). The system measures spectrally resolved pond reflectance by ratioing upwelling radiance from each raceway to simultaneous downwelling sky irradiance using fiber-coupled spectrometers. A physics-based reflectance model (ASHARP) is fit to each spectrum pair to estimate optical parameters, including a biomass-proxy coefficient (C a ) which enables near-real-time tracking of biomass accumulation and culture state at 2–5 min intervals. From May 2022 through September 2025 the platform operated continuously while scaling from two to six raceway ponds. Several strains of algae were monitored successfully, including the high productivity Tetraselmis striata and Picochlorum celeri. Transitioning data acquisition from a Windows laptop to a Raspberry Pi improved uptime from 57% (2022) to ~89% (2024–2025) and enabled routine real-time analysis. Further, we converted relative biomass estimates to absolute ash-free dry weight (AFDW) using experimentally-derived calibrations, providing field-relevant biomass predictions with conservative confidence bounds. These results demonstrate the feasibility of long-term, autonomous optical monitoring for well-mixed open-raceway algal cultivation and provide practical guidance for reliable field operation and scaling.

Katinas, Christopher Michael [Sandia National Labo↗
DOE OSTI2025

Energy efficiency in industrial drying: A hybrid ultrasonic system with a novel dynamic optimization framework

Drying processes are among the most energy-consuming operations in industrial and manufacturing settings, demanding strategic selection, design, and control for enhanced efficiency. Advancing drying technologies is critical for improving sustainability, lowering energy use, reducing carbon emissions, and minimizing waste. This study explores two innovative strategies aimed at transforming drying processes into sustainable, low-carbon systems by reducing energy consumption, minimizing waste, and maintaining a strong emphasis on preserving product quality. The first strategy showcases a sub-pilot scale hybrid ultrasonic-convective dryer for agrifood products. This technology, powered by electricity (process electrification), integrates non-thermal ultrasonic dehydration with convective heating and is presented as a sustainable and energy-efficient solution that enhances eco-friendly practices. The second strategy involves introducing and implementing a novel, multiobjective, mixed integer dynamic optimization technique to determine the optimal time-dependent process parameter values for the drying operation. This optimization technique yields operating conditions that are piecewise constant in time aiming to maximize the energy efficiency of the hybrid ultrasonic-convective dryer while ensuring strict adherence to product quality constraints. By adopting the hybrid ultrasonic-convective dryer, a notable 35% improvement in energy efficiency was achieved compared to conventional hot-air drying systems for drying apple slices. The proposed optimization framework further enhanced energy efficiency by nearly 14% over the most efficient process on the identical testbed, under static operating conditions. The reported enhancements have been experimentally validated. Regarding drying time (thereby improving production yield), the developed hybrid ultrasonic-convective dryer demonstrates as much as a 41% reduction in total processing time, which is further optimized by an additional 10% using our proposed optimization framework. The research outcomes have profound implications for the design and operation of drying systems, encompassing crucial aspects such as process electrification, cost-effectiveness, energy savings, time efficiency, product yield, product quality, and process automation.

Dynamic optimization↗
DOE OSTI2021

Review of analytical studies on seismic fluid-structure interaction of base-supported cylindrical tanks

Analytical solutions for seismic fluid-structure interaction (FSI) of tanks are used by analysts and engineers for a number of purposes, including preliminary sizing, risk assessment, and verification of numerical models. Analytical solutions for seismic FSI of tanks were first developed in the 1930s and have been extended and/or modified in the following nine decades to accommodate different seismic inputs and boundary conditions. This report reviews published analytical solutions for rigid and flexible, base-supported, cylindrical tanks subjected to unidirectional horizontal motion of a small amplitude. The studies parsed the fluid-structure (FSI) response into an impulsive and a convective component, and addressed frequencies of lateral motions of the tank (i.e., impulsive frequencies), hydrodynamic pressures, frequencies of waves (i.e., convective frequencies), wave heights, and reactions at the base. The solutions were categorized into exact and approximate solutions based on the methodologies used. Herein, FSI responses are calculated and compared for a range of tank dimensions using the analytical solutions from the different studies. The responses presented here are normalized to be unitless, and can be used for tanks with different dimensions and mechanical properties, and subjected to different input motions. A sample steel water tank is analyzed using the normalized solutions for two earthquake ground motions with different frequency contents. Ignoring tank flexibility may significantly underestimate seismic responses of the tank and its contained fluid. The approximate solutions, which have been widely used for decades, may not be sufficiently accurate in the modern era. Calculation errors in the reviewed studies are identified and corrected, and differences are quantified for the sample tank.

42 ENGINEERING↗
DOE OSTI2020

Tracking nickel oxide reduction in solid oxide cells via ex-situ ptychographic nano-tomography

In solid oxide cell (SOC) electrodes, nickel oxide reduction is a fundamental step occurring during the initial stages of cell operation. The reduction procedure is considered to be crucial in producing the final fuel electrode microstructure. This work studies the reduction reaction occurring in a SOC fuel electrode by means of X-ray ptychographic nanotomography. The same microstructure is analyzed in its pristine state and after 1 h of reduction at 850 °C, in dry hydrogen. The microstructure evolution is discussed in terms of changes in volume fractions, surface areas, and connectivity of the different constituent materials. The pore network in the unreduced cermet is revealed to be fully isolated, implying that reduction proceeds inwards from the edges of the sample towards the center. For the given experimental conditions we report that reduction is complete after 1 h of treatment but the resulting nickel network is not yet percolated and detached from the yttria-stabilized zirconia.

25 ENERGY STORAGE↗
DOE OSTI2025

Mechanical properties and strengthening mechanisms of cold sprayed Al-Mg alloys

This paper investigates the mechanical properties and strengthening mechanisms of cold sprayed Al-Mg alloys as a function of Mg content. Cold spray (CS) is a solid-state deposition process that can be used as a repair technique for structural components. However, the limited ductility of cold sprayed materials is an unsettled issue in CS community. Herein, reproducible, and significantly increased elongation to failure (∼15%) in as-deposited conditions is reported in CS Al-Mg binary alloys with low Mg content (≤ 2.5 wt%). Furthermore, reducing Mg content down to 2.5 wt% in CS Al-Mg binary alloys demonstrated comparable mechanical strength to its wrought counterparts. Strengthening mechanisms were analyzed in wrought and CS Al-Mg alloys based upon microstructural characteristics, such as grain size, Mg solute concentration, and dislocation density. A quantitative evaluation of yield strength revealed that grain boundary strengthening is predominant in CS Al-Mg alloys as compared with their wrought counterparts. When comparing among CS Al-Mg binary alloys of different compositions, Mg solute concentration, and dislocation density governed mechanical strength. In conclusion, the underlying mechanism for the variation of mechanical properties and strengthening mechanisms is discussed in terms of these microstructural characteristics.

Additive manufacturing↗
DOE OSTI2023

Non-Centrosymmetric Sr 2 IrO 4 Obtained Under High Pressure

Sr 2 IrO 4 with strong spin-orbit coupling (SOC) and Hubbard repulsion (U) hosts Mott insulating states. The similar crystal structure, magnetic and electronic properties, particularly the d-wave gap observed in Sr 2 IrO 4 enhanced the analogies to cuprate high-$T_c$ superconductor, La 2 CuO 4 . The incomplete analogy was due to the lack of broken inversion symmetry phases observed in Sr 2 IrO 4 . Here, under high pressure and high temperature conditions, we report a non-centrosymmetric Sr 2 IrO 4 . The crystal structure and its noncentrosymmetric character were determined by single crystal X-ray diffraction and high-resolution scanning transmission electron microscopy (HR-STEM). The magnetic characterization confirms the Ir 4+ with $\textit{S}$ = 1/2 at low temperature in Sr 2 IrO 4 with magnetic ordering occurred at around 86 K, where a larger moment is observed than the ambient pressure Sr 2 IrO 4 . Moreover, the resistivity measurement shows three-dimensional Mott variable-range hopping existed in the system. Further, this non-centrosymmetric Sr 2 IrO 4 phase appears to be a unique material to offer further understanding of high-$T_c$ superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Estimation of the Kelvin–Helmholtz Unstable Boundary

The Kelvin–Helmholtz (KH) instability is one of the most important mechanisms of the viscous-like interaction between the solar wind and the magnetosphere (MSP), which transport the mass, energy, momentum, and magnetic flux. Thus, it is important to examine whether the magnetopause boundary is KH unstable or not. Based on the KH onset conditions, this report proposes to use a matrix to identify the most KH unstable direction based on the in situ measurements of the density, velocity, and magnetic field in the MSP and magnetosheath. The range of the KH unstable direction can be easily estimated based on the eigenvalues of the matrix. The eigenvectors of the matrix provide a new boundary normal coordinate system, which could be useful for 2-D KH instability simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Magnetic field screening in hydrogen-rich high-temperature superconductors

In the last few years, the superconducting transition temperature, T c , of hydrogen-rich compounds has increased dramatically, and is now approaching room temperature. However, the pressures at which these materials are stable exceed one million atmospheres and limit the number of available experimental studies. Superconductivity in hydrides has been primarily explored by electrical transport measurements, whereas magnetic properties, one of the most important characteristic of a superconductor, have not been satisfactory defined. Here, we develop SQUID magnetometry under extreme high-pressure conditions and report characteristic superconducting parameters for Im-3m-H 3 S and Fm-3m-LaH 10 —the representative members of two families of high-temperature superconducting hydrides. We determine a lower critical field H c1 of ~0.82 T and ~0.55 T, and a London penetration depth λ L of ~20 nm and ~30 nm in H 3 S and LaH 10 , respectively. The small values of λ L indicate a high superfluid density in both hydrides. These compounds have the values of the Ginzburg-Landau parameter κ ~12–20 and belong to the group of “moderate” type II superconductors, rather than being hard superconductors as would be intuitively expected from their high T c s.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Structural determinants for pH-dependent activation of a plant metacaspase

Arabidopsis thaliana metacaspase 9 (AtMC9) plays roles in clearing dead cells, forming xylem vessels, and regulating immunity and programmed cell death in plants. The protease's activation is controlled by pH levels, but the exact structural mechanism behind this has not been elucidated. In this work, high-resolution crystal structures for AtMC9 at active (pH 5.5 and pH 4.2) and inactive (pH 7.5) conditions are reported. The three structures are similar except for local conformations where their hydrogen bonding interactions with solvents are mediated through the protonation of specific titratable amino acid residues' side chains. By combining structural analysis, molecular dynamics simulations under constant pHs, and biochemical assays coupled with site-directed mutagenesis, we show that the regulation of AtMC9 activation involves multiple titratable glutamate and histidine residues across the three domains of p20, linker, and p10. Specifically, deprotonated Glu112, His193, and His208 can suppress AtMC9 proteolytic activity, while protonation of Glu255 and His307 at acidic pH may promote it. This study provides valuable insights into the pH-dependent activation of AtMC9 and could potentially lead to improving crops with enhanced immunity and controlled cell death, ultimately increasing agricultural productivity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

“Breathing” organic cation to stabilize multiple structures in low-dimensional Ge-, Sn-, and Pb-based hybrid iodide perovskites

Low-dimensional hybrid inorganic–organic perovskites are excellent candidates for stable optoelectronic devices. The dimensionality of these perovskites depends largely on the organic and inorganic compositions, as well as the synthetic conditions. We report five new hybrid iodides, (ETU) 4 Ge 5 I 18 , (ETU)GeI 4 , (ETU)SnI 4 , (ETU)PbI 4 , and (ETU) 3 Pb 2 I 10 using only one type of organic cation, namely, S-(2-aminoethyl)isothiouronium (ETU). (ETU)GeI 4 and (ETU)SnI 4 belong to the (110)-oriented structure-type with “3 × 3” sawtooth corrugated layers and crystallize in a structure with the orthorhombic space group Pbca. (ETU) 4 Ge 5 I 18 crystallizes in a structure with the triclinic space group P$\bar{1}$with combining macron], featuring a 2D layered structure with combinations of corner, edge, and face-sharing [GeI 6 ] octahedra. For the Pb-based series, (ETU)PbI 4 has the conventional (100) – oriented 2D type whereas (ETU) 3 Pb 2 I 10 has a unique 0D structure. Remarkably, the unstable 2D orange-phase (ETU)PbI 4 transforms to a stable 0D yellow phase (ETU) 3 Pb 2 I 10 , accompanied by the reduction of the C–S–C angle of the organic cation ETU. The optical band gaps are largely regulated by the diverse types of structure and are in the range of 1.8 eV to 2.8 eV. (ETU)SnI 4 is the only material showing notable photoluminescence at room-temperature. Our work showcases the flexibility of the organic cation in determining the structural dimensionality and provides a new strategy in generating new hybrid materials.

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

Analytic stability boundaries for compressional and global Alfvén eigenmodes driven by fast ions. II. Interaction via Landau resonance

We report that conditions for net fast ion drive are derived for beam-driven, co-propagating, sub-cyclotron compressional (CAE) and global (GAE) Alfvén eigenmodes driven by the Landau resonance with super-Alfvénic fast ions. Approximations applicable to realistic neutral beam distributions and mode characteristics observed in spherical tokamaks enable the derivation of marginal stability conditions for these modes. Such conditions successfully reproduce the stability boundaries found from numerical integration of the exact expression for local fast ion drive/damping. Coupling between the CAE and GAE branches of the dispersion due to finite ω/ω ci and |k ∥ /k ⊥ | is retained and found to be responsible for the existence of the GAE instability via this resonance. Finally, encouraging agreement is demonstrated between the approximate stability criterion, simulation results, and a database of NSTX observations of co-CAEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Validation of heat transport modeling using directly driven beryllium spheres

Recent experiments involving directly driven beryllium spheres are reported. Plasma conditions are measured using Thomson scattering with the probe beam pointed 200, 300, and 400 μ m from the surface of the sphere. Laser coupling is assessed using calorimeters that collect scattered light placed at various locations within the target chamber. Laser intensities of 1 0 14 W / c m 2 and 2.5 × 1 0 14 W / c m 2 are chosen to minimize unmodeled laser-plasma interactions (LPIs) that lead to laser-target decoupling. Two-dimensional simulations are compared to the interpreted data using the radiation-hydrodynamics code Lasnex. Heat transport is simulated using flux-limited Spitzer–Harm with both high ( f = 0.15) and low ( f = 0.03) flux limiters and the Schurtz–Nicolai-Busquet (SNB) model. At 1 0 14 W / c m 2 , all three heat transport models agree well with the measurement, demonstrating that the heat flux is local at low intensities near the measurement locations. At 2.5 × 1 0 14 W / c m 2 , the SNB and high flux model roughly match the plasma conditions but predict 2% uncoupled light compared to 10% measured. Additionally, the use of drive multipliers to match the measured coupled light does not alter the agreement between measured and simulated plasma conditions, suggesting that decoupling due to LPI is unlikely to alter this agreement. The low flux model cannot match the plasma conditions and results in 19% scattered light. The use of a resonant absorption model can be used to bring the simulated scattered light into agreement, but the simulated plasma conditions are still in disagreement with the measurement. For this reason, the low flux model is rejected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Solid organic-coated ammonium sulfate particles at high relative humidity in the summertime Arctic atmosphere

The ability of atmospheric aerosols to impact climate through water uptake and cloud formation is fundamentally determined by the size, composition, and phase (liquid, semisolid, or solid) of individual particles. Particle phase is dependent on atmospheric conditions (relative humidity and temperature) and chemical composition and, importantly, solid particles can inhibit the uptake of water and other trace gases, even under humid conditions. Particles composed primarily of ammonium sulfate are presumed to be liquid at the relative humidities (67 to 98%) and temperatures (–2 to 4 °C) of the summertime Arctic. Under these atmospheric conditions, we report the observation of solid organic-coated ammonium sulfate particles representing 30% of particles, by number, in a key size range (<0.2 µm) for cloud activation within marine air masses from the Arctic Ocean at Utqiagvik, AK. The composition and size of the observed particles are consistent with recent Arctic modeling and observational results showing new particle formation and growth from dimethylsulfide oxidation to form sulfuric acid, reaction with ammonia, and condensation of marine biogenic sulfate and highly oxygenated organic molecules. Aqueous sulfate particles typically undergo efflorescence and solidify at relative humidities of less than 34%. Therefore, the observed solid phase is hypothesized to occur from contact efflorescence during collision of a newly formed Aitken mode sulfate particle with an organic-coated ammonium sulfate particle. With declining sea ice in the warming Arctic, this particle source is expected to increase with increasing open water and marine biogenic emissions.

54 ENVIRONMENTAL SCIENCES↗