Search NASA⌕ Search

SEARCH · Search NASA

Results for “emission modeling”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The Emissions Model Intercomparison Project (Emissions-MIP): quantifying model sensitivity to emission characteristics

Abstract. Anthropogenic emissions of aerosols and precursor compounds are known to significantly affect the energy balance of the Earth–atmosphere system, alter the formation of clouds and precipitation, and have a substantial impact on human health and the environment. Global models are an essential tool for examining the impacts of these emissions. In this study, we examine the sensitivity of model results to the assumed height of SO2 injection, seasonality of SO2 and black carbon (BC) particulate emissions, and the assumed fraction of SO2 emissions that is injected into the atmosphere as particulate phase sulfate (SO4) in 11 climate and chemistry models, including both chemical transport models and the atmospheric component of Earth system models. We find large variation in atmospheric lifetime across models for SO2, SO4, and BC, with a particularly large relative variation for SO2, which indicates that fundamental aspects of atmospheric sulfur chemistry remain uncertain. Of the perturbations examined in this study, the assumed height of SO2 injection had the largest overall impacts, particularly on global mean net radiative flux (maximum difference of −0.35 W m−2), SO2 lifetime over Northern Hemisphere land (maximum difference of 0.8 d), surface SO2 concentration (up to 59 % decrease), and surface sulfate concentration (up to 23 % increase). Emitting SO2 at height consistently increased SO2 and SO4 column burdens and shortwave cooling, with varying magnitudes, but had inconsistent effects across models on the sign of the change in implied cloud forcing. The assumed SO4 emission fraction also had a significant impact on net radiative flux and surface sulfate concentration. Because these properties are not standardized across models this is a source of inter-model diversity typically neglected in model intercomparisons. These results imply a need to ensure that anthropogenic emission injection height and SO4 emission fraction are accurately and consistently represented in global models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Corrigendum to “The Emissions Model Intercomparison Project (Emissions-MIP): quantifying model sensitivity to emission characteristics” published in Atmos. Chem. Phys., 23, 14779–14799, 2023

Anthropogenic emissions of aerosols and precursor compounds are known to significantly affect the energy balance of the Earth–atmosphere system, alter the formation of clouds and precipitation, and have a substantial impact on human health and the environment. Global models are an essential tool for examining the impacts of these emissions. In this study, we examine the sensitivity of model results to the assumed height of SO 2 injection, seasonality of SO 2 and black carbon (BC) particulate emissions, and the assumed fraction of SO 2 emissions that is injected into the atmosphere as particulate phase sulfate (SO 4 ) in 11 climate and chemistry models, including both chemical transport models and the atmospheric component of Earth system models. We find large variation in atmospheric lifetime across models for SO 2 , SO 4 , and BC, with a particularly large relative variation for SO 2 , which indicates that fundamental aspects of atmospheric sulfur chemistry remain uncertain. Of the perturbations examined in this study, the assumed height of SO 2 injection had the largest overall impacts, particularly on global mean net radiative flux (maximum difference of -0.35?W?m-2), SO2 lifetime over Northern Hemisphere land (maximum difference of 0.8?d), surface SO 2 concentration (up to 59?% decrease), and surface sulfate concentration (up to 23?% increase). Emitting SO 2 at height consistently increased SO 2 and SO 4 column burdens and shortwave cooling, with varying magnitudes, but had inconsistent effects across models on the sign of the change in implied cloud forcing. The assumed SO 4 emission fraction also had a significant impact on net radiative flux and surface sulfate concentration. Because these properties are not standardized across models this is a source of inter-model diversity typically neglected in model intercomparisons. These results imply a need to ensure that anthropogenic emission injection height and SO 4 emission fraction are accurately and consistently represented in global models.

Ahsan, Hamza [Pacific Northwest National Laborator↗
DOE OSTI2024

A Lake Biogeochemistry Model for Global Methane Emissions: Model Development, Site‐Level Validation, and Global Applicability

Abstract Lakes are important sentinels of climate change and may contribute over 30% of natural methane (CH 4 ) emissions; however, no earth system model (ESM) has represented lake CH 4 dynamics. To fill this gap, we refined a process‐based lake biogeochemical model to simulate global lake CH 4 emissions, including representation of lake bathymetry, oxic methane production (OMP), the effect of water level on ebullition, new non‐linear CH 4 oxidation kinetics, and the coupling of sediment carbon pools with in‐lake primary production and terrigenous carbon loadings. We compiled a lake CH 4 data set for model validation. The model shows promising performance in capturing the seasonal and inter‐annual variabilities of CH 4 emissions at 10 representative lakes for different lake types and the variations in mean annual CH 4 emissions among 106 lakes across the globe. The model reproduces the variations of the observed surface CH 4 diffusion and ebullition along the gradients of lake latitude, depth, and surface area. The results suggest that OMP could play an important role in surface CH 4 diffusion, and its relative importance is higher in less productive and/or deeper lakes. The model performance is improved for capturing CH 4 outgassing events in non‐floodplain lakes and the seasonal variability of CH 4 ebullition in floodplain lakes by representing the effect of water level on ebullition. The model can be integrated into ESMs to constrain global lake CH 4 emissions and climate‐CH 4 feedback.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A new process-based and scale-aware desert dust emission scheme for global climate models – Part I: Description and evaluation against inverse modeling emissions

Desert dust accounts for most of the atmosphere's aerosol burden by mass and produces numerous important impacts on the Earth system. However, current global climate models (GCMs) and land-surface models (LSMs) struggle to accurately represent key dust emission processes, in part because of inadequate representations of soil particle sizes that affect the dust emission threshold, surface roughness elements that absorb wind momentum, and boundary-layer characteristics that control wind fluctuations. Furthermore, because dust emission is driven by small-scale (~ 1 km or smaller) processes, simulating the global cycle of desert dust in GCMs with coarse horizontal resolutions (~ 100 km) presents a fundamental challenge. This representation problem is exacerbated by dust emission fluxes scaling nonlinearly with wind speed above a threshold wind speed that is sensitive to land-surface characteristics. Here, we address these fundamental problems underlying the simulation of dust emissions in GCMs and LSMs by developing improved descriptions of (1) the effect of soil texture on the dust emission threshold, (2) the effects of nonerodible roughness elements (both rocks and green vegetation) on the surface wind stress, and (3) the effects of boundary-layer turbulence on driving intermittent dust emissions. We then use the resulting revised dust emission parameterization to simulate global dust emissions in a standalone model forced by reanalysis meteorology and land-surface fields. We further propose (4) a simple methodology to rescale lower-resolution dust emission simulations to match the spatial variability of higher-resolution emission simulations in GCMs. The resulting dust emission simulation shows substantially improved agreement against regional dust emissions observationally constrained by inverse modeling. We thus find that our revised dust emission parameterization can substantially improve dust emission simulations in GCMs and LSMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Improving commercial truck fleet composition in emission modeling using 2021 US VIUS data

Commercial trucks are essential elements of the nation's supply chain system. Meanwhile, intensive truck movements contribute significantly to system externalities, such as energy use and air pollution. However, collecting detailed fleet composition and distribution of operational patterns remains a barrier to accurately accounting for these impacts. The recently released 2021 US Vehicle Inventory and Use Survey (US VIUS) fills a critical gap in understanding commercial truck fleet distributions, their operations, and business constraints at the national scale. This study aims to understand the latest US commercial vehicle fleet composition and operational characteristics using 2021 US VIUS data and calibrate the fleet inputs in regulatory emission models to assess the potential emission implications of the VIUS-derived fleet composition. The emission rates for commercial trucks and default fleet composition are collected from the U.S. EPA's MOtor Vehicle Emission Simulator (MOVES4). The 2021 US VIUS data is applied to improve fleet characteristics such as the long-haul fraction and the vehicle mileage accumulation rate. The study also investigates potential emission reduction benefits under various forecasted fleet electrification scenarios. The energy consumption and critical air pollutant rates by vehicle types are compared between MOVES4 and US VIUS fleets for both current and future scenarios to provide insights into the latest U.S. commercial vehicle fleet characteristics and their implications on energy and emissions. This study helps policymakers and practitioners advance the commercial fleet generation for emission models. It also deepens the understanding of the emission reduction potential of the commercial fleet under various fleet projections.

2021 US VIUS↗
DOE OSTI2025

Analyzing Tradeoffs Associated with the Manufacture of Refined Fuels by Comparing Volatile Organic Compound (VOC) Emissions Model Results with Carbon Impacts

A first-of-its kind, publicly-available model for estimating emissions of volatile organic compounds from operations along the supply chain of liquid fuels and refinery products is presented and demonstrated. This standardized model permits comprehensive, reproducible, and comparable evaluation of the supply chain/life cycle VOC emissions of liquid fuels and refinery products.

09 BIOMASS FUELS↗
DOE OSTI2026

Electromagnetic Analysis of ITER Electron Cyclotron Emission Model

The ITER electron cyclotron emission (ECE) diagnostic system is located at diagnostics shielding module 2 (DSM2), equatorial port 9 (EP9), to measure electron temperature profile and electron temperature fluctuations, and also assess nonthermal electron distributions via the oblique view. Therefore, ECE has both radial and oblique views with two hot sources for calibration and two couples of mirrors for two different optical views. It shall receive that the large electromagnetic (EM) loads up to 100-MN/m3 force density due to the eddy currents on the copper mirrors generated by the short transient plasma disruption. The simplified EM analysis model using magnetic field ( B ) data and flux variations (dB/dt) method based on the worst case of plasma disruption, MD_DW_EXP16MS_CATIII, has been developed to calculate the force and moments for each ECE component and support structure for EP9 DSM2 Bay2 and Bay3. The B and dB/dt methods use the constant B and dB/dt assumed for each local bay space, but the real B and dB/dt are basically varied along the radial direction of each bay. The complex global model with the local ECE components and support structure using Maxwell transient has been developed to do EM analysis to compare the results of the B and dB/dt models. The volumetric EM force density of the whole ECE components and support structure can be used for the subsequent structural analysis. In conclusion, combined with the thermal and nuclear loads and seismic and initial loads, the ECE integration analysis can be finalized for the operation case.

Fang, J. [Princeton Plasma Physics Laboratory (PPP↗
DOE OSTI2021

Comparison of simulated neutrino emission models with data on Supernova 1987A

Here we compare models of supernova (SN) neutrino emission with the Kamiokande II data on SN 1987A using the Bayesian approach. These models are taken from simulations and are representative of current one-dimensional SN models. We find that models with a brief accretion phase of neutrino emission are the most favored. This result is not affected by varying the overall flux normalization or considering neutrino oscillations. We also check the compatibility of the best-fit models with the data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

General kinetic ion-induced electron emission model for metallic walls applied to biased Z-pinch electrodes

A kinetic ion-induced electron emission (IIEE) model for general applications is developed to obtain the emitted electron energy spectrum for a distribution of ion impacts on a metallic surface. We assume an ionization cascade mechanism and use empirical models for the ion and electron stopping powers. The emission spectrum and the secondary electron yield (SEY) are validated for a variety of materials. The IIEE model is used to study the effect of IIEE on the plasma-material interactions of Z-pinch electrodes. Un-magnetized Boltzmann-Poisson simulations are performed for a Z-pinch plasma doubly bounded by two biased copper electrodes with and without IIEE at bias potentials from 0 to 9 kV. At the anode, the SEY decreases from 0 to 1 kV, but then increases at higher bias potentials. At the cathode, the SEY is much larger due to higher energy ion bombardment and grows with bias potential. As the bias potential increases, the emitted cathode electrons are accelerated to higher energies into the domain, collisionally heating the plasma. Above 1 kV, the heating is strong enough to increase the plasma potential. Despite SEY greater than 1, only a classical sheath forms as opposed to a space-charge limited or inverse sheath due to the emitted electron flux not reaching the space charge current saturation limits. Furthermore, the current in the emissionless cases saturates to a value lower than experiment. With IIEE, the current does not saturate and continues to increase with the 4 kV case, matching most closely with the experiment.

Carbon based materials↗
DOE OSTI2021

Modeled Emission of Hydroxyl and Ozone Reactivity from Evaporation of Fragrance Mixtures

Volatile chemical products (VCPs) account for increasing fractions of organic carbon emitted to the atmosphere, particularly in urban areas. Fragrances are potentially reactive components that are added to many VCPs. To better constrain these emissions, 11 commercially available liquid fragrance mixtures were characterized for their composition and their evaporation modeled. Here, emissions of mass, hydroxyl reactivity, and ozone reactivity were estimated by modeling under four different scenarios. Fragrance compounds were generally less than one-half the mass of fragrance mixtures, with the balance comprised of solvents and plasticizers and unresolved mass thought to be dominated by plasticizers. The results showed that terpenes and terpenoids account for nearly all of the emitted mass and reactivity while only comprising ~10% w/w on average of the liquid fragrance mixtures. Most of the reactivity is emitted within hours, with ozone reactivity evolving more rapidly than OH reactivity and comprised almost entirely of terpenes. Limonene, a common fragrance constituent, dominates the reactivity of emitted carbon. Generally, 20–40% of the potential hydroxyl reactivity contained in the fragrance mixture does not evaporate on time scales sufficient to have an impact on local or regional air quality.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Simulation of radon-222 with the GEOS-Chem global model: Emissions, seasonality, and convective transport

Radon-222 (222Rn) is a short-lived radioactive gas naturally emitted from land surface, and has long been used to assess convective transport in atmospheric models. In this study, we simulate 222Rn using the GEOS-Chem chemical transport model with aims to improve our understanding of 222Rn emissions and surface concentration seasonality, and to characterize convective transport associated with two Goddard Earth Observing System (GEOS) meteorological products, MERRA and GEOS-FP. We evaluate four available 222Rn emission scenarios by comparing model results with surface observations at 51 global surface sites. The default emission scenario in GEOS-Chem yields a moderate agreement with global surface observations (< 70% data within a factor of 2) and a large underestimate of wintertime surface 222Rn concentrations at Northern Hemisphere mid- and high-latitudes due to an oversimplified formulation of 222Rn emission fluxes (1 atom cm-2 s-1 over land with a reduction by a factor of 3 under freezing conditions). We compose a new global 222Rn emission scenario based on Zhang et al. (2011) and show its potential to improve simulated surface 222Rn concentrations and seasonality. The regional components of this emission scenario include spatially and temporally varying emission fluxes derived from previous measurements of soil radium content and soil exhalation models, which are key factors to 222Rn emission flux rates. However, large model underestimates of surface 222Rn concentrations still exist in Asia, suggesting unusually high regional 222Rn emissions. We propose a conservative up-scaling factor of 1.2 for 222Rn emission fluxes in China, as also constrained by the observed deposition fluxes of 210Pb (decay daughter of 222Rn). With this modification, the model shows better agreement with the observations in Europe and North America (>80% data within a factor of 2), and reasonable agreement in Asia (close to 70%). Further constraints on 222Rn emissions would require additional observations of surface 222Rn concentrations and emission fluxes in central U.S., Canada, Africa, and Asia. We also compare and assess convective transport in model simulations driven by MERRA and GEOS-FP using observed 222Rn vertical profiles during northern mid-latitude summertime and from two short-term airborne campaigns. While the simulations with both GEOS products are able to capture the observed vertical gradient of 222Rn concentrations in the lower troposphere (0-4 km), neither correctly represents the level of convective detrainment, resulting in biases in the middle and upper troposphere. Compared to GEOS-FP, MERRA leads to stronger convective transport of 222Rn, which is partially compensated by its weaker large-scale vertical advection, resulting in similar global vertical distributions of 222Rn concentrations between the two simulations.

Zhang, Bo↗
DOE OSTI2025

Impact of a dynamic grid mix and climate on operational carbon emissions modeling for different building typologies and climate zones

Calculating operational carbon emissions through a building’s lifecycle is complex due to the dynamic nature of influencing factors such as climate and energy grid mix. This paper introduces a novel methodology for modeling 30-year operational carbon impacts of buildings and applies this method to mid-rise office and residential typologies across various US climate zones. The method accounts for these temporal variabilities using new and scarcely cited data sources. Key findings indicate that future changes in the climate, while impactful, play a relatively modest role in operational carbon emissions compared to significant reductions with modeling scenarios using the projected decarbonization of the electricity grid. Here, the study also finds that using annual, month-hourly, or hourly grid emission factors have a minimal impact on carbon accounting, except in certain climates and program types where emission patterns do not align with a building’s energy consumption. Warmer climates like Miami, Florida and Tucson, Arizona, which rely heavily on cooling, demonstrate larger variations in carbon emissions when using higher temporal resolution emission factors. Ultimately, this study underscores the critical role of grid decarbonization in reducing long-term emissions and the importance of incorporating this variable in life cycle assessment (LCA) modeling.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Development of Methane Emissions Model to Assess Fuel Recovery Potential at Gas Well Sites Using On Site Compression

The U.S. natural gas production and consumption has increased 85.5% since 2005 primarily due to the unconventional production methods of horizontal drilling and hydraulic fracturing. Natural gas used as a fuel has a lower greenhouse gas (GHG) footprint than coal and petroleum due to lower Carbon Dioxide (CO2) emissions when combusted. However, the “greener” benefit to natural gas may be negated by leaks in production and transmission systems. Methane (CH4), the primary hydrocarbon in natural gas, has an estimated Global Warming Potential (GWP) of 28-36 over 100 years, meaning it can absorb 28-36 more energy than CO2 which has a GWP of 1.0. Natural gas well sites are prone to methane emissions, or leaks and irregular gas releases, vented to atmosphere throughout production and transmission. The U.S. Department of Energy (DOE) and the National Energy Technology Laboratory (NETL) has recently granted West Virginia University (WVU) funding under agreement DE-FOA-0002005, to “Advance technologies to mitigate methane emissions and increase the efficiency of the natural gas transportation infrastructure”. As part of this funding WVU was tasked with identifying and quantifying sources of methane emissions at unconventional well sites, processing this data, and developing a system to recapture these emissions. A 0-D Simulink model was developed, utilizing standardized methodologies, data from previously conducted studies, as well as collected data from well sites in the Marcellus shale play region. The model was developed to predict emission rates from various components at natural gas well sites as well as the potential to utilize these emissions as fuel for the natural gas powered compressor engines on-site. This model was utilized to run high, medium, and low cases for four identified emission sources, engine size, pneumatic controller count, liquid level production which dictates tank emissions, and compressor packing vent emissions. Due to discrepancies in transient tank emission data, a high and low emission factor for tanks was used, resulting in two sets of 81 executed cases, and 162 unique cases of total site emissions and potential for fuel consumption. Each of the cases were run over 86,400 seconds at a 1 Hz, representative of a full 24 hour day of operation. The fuel consumption offset an average of 557% of fuel consumption on an energy density basis across all 81 cases with the high tank emission factor with a maximum offset of 2334%. The fuel consumption offset was an average of 82.9% for all 81 cases with the low tank emission factor with a maximum offset of 337%. This study highlights flaws in the use of publicly available methane number calculations to determine natural gas’s suitability as an engine fuel as well as the lack of public data for transient liquid storage tank emissions.

03 NATURAL GAS↗
DOE OSTI2024

Energy-dependent implementation of secondary electron emission models in continuum kinetic sheath simulations

The plasma-material interactions present in multiple fusion and propulsion concepts between the flow of plasma through a channel and a material wall drive the emission of secondary electrons. This emission is capable of altering the fundamental structure of the sheath region, significantly changing the expected particle fluxes to the wall. The emission spectrum is separated into two major energy regimes, a peak of elastically backscattered primary electrons at the incoming energy, and cold secondary electrons inelastically emitted directly from the material. The ability of continuum kinetic simulations to accurately represent the secondary electron emission is limited by relevant models being formulated in terms of monoenergetic particle interactions which cannot be applied directly to the discrete distribution function. As a result, rigorous implementation of energy-dependent physics is often neglected in favor of simplified, constant models. We present here a novel implementation of semi-empirical models in the boundary of continuum kinetic simulations which allows the full range of this emission to be accurately captured in physically-relevant regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Parameterized Neutrino Emission Model to Study Mass Ejection in Failed Core-collapse Supernovae

Some massive stars end their lives as failed core-collapse supernovae (CCSNe) and become black holes (BHs). Although in this class of phenomena the stalled supernova (SN) shock is not revived, the outer stellar envelope can still be partially ejected. This occurs because the hydrodynamic equilibrium of the star is disrupted by the gravitational mass loss of the protoneutron star (PNS) due to neutrino emission. We develop a simple parameterized model that emulates PNS evolution and its neutrino emission and use it to simulate failed CCSNe in spherical symmetry for a wide range of progenitor stars. Our model allows us to study mass ejection of failed CCSNe where the PNS collapses into a BH within ~100 ms and up to ~10 6 s. We perform failed CCSNe simulations for 262 different pre-SN progenitors and determine how the energy and mass of the ejecta depend on progenitor properties and the equation of state (EOS) of dense matter. In the case of a future failed CCSN observation, the trends obtained in our simulations can be used to place constraints on the pre-SN progenitor characteristics, the EOS, and on PNS properties at BH formation time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A Broad Grid of 2D Kilonova Emission Models

Depending upon the properties of their compact remnants and the physics included in the models, simulations of neutron star mergers can produce a broad range of ejecta properties. The characteristics of this ejecta, in turn, define the kilonova emission. To explore the effect of ejecta properties, we present a grid of 2-component 2D axisymmetric kilonova simulations that vary mass, velocity, morphology, and composition. The masses and velocities of each component vary, respectively, from 0.001 to 0.1 M ⊙ and 0.05 to 0.3c, covering much of the range of results from the neutron star merger literature. The set of 900 models is constrained to have a toroidal low electron fraction (Y e ) ejecta with a robust r-process composition and either a spherical or lobed high-Y e ejecta with two possible compositions. We simulate these models with the Monte Carlo radiative transfer code SuperNu using a full suite of lanthanide and 4th row element opacities. We examine the trends of these models with parameter variation, show how it can be used with statistical tools, and compare the model light curves and spectra to those of AT2017gfo, the electromagnetic counterpart of GW170817.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Opacities of singly and doubly ionized neodymium and uranium for kilonova emission modeling

ABSTRACT Even though the electromagnetic counterpart AT2017gfo to the binary neutron star merger GW170817 is powered by the radioactive decay of r-process nuclei, only few tentative identifications of light r-process elements have been made so far. One of the major limitations for the identification of heavy nuclei is incomplete or missing atomic data. While substantial progress has been made on lanthanide atomic data over the last few years, for actinides there has been less emphasis, with the first complete set of opacity data only recently published. We perform atomic structure calculations of neodymium (Z = 60) as well as the corresponding actinide uranium (Z = 92). Using two different codes [flexible atomic code (fac) and hartree–fock-relativistic (hfr)] for the calculation of the atomic data, we investigate the accuracy of the calculated data (energy levels and electric dipole transitions) and their effect on kilonova opacities. For the fac calculations, we optimize the local central potential and the number of included configurations and use a dedicated calibration technique to improve the agreement between theoretical and available experimental atomic energy levels (AELs). For ions with vast amounts of experimental data available, the presented opacities agree quite well with previous estimations. On the other hand, the optimization and calibration method cannot be used for ions with only few available AELs. For these cases, where no experimental nor benchmarked calculations are available, a large spread in the opacities estimated from the atomic data obtained with the various atomic structure codes is observed. We find that the opacity of uranium is almost double the neodymium opacity.

Flörs, A. (ORCID:0000000320242819)↗