Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transport models”

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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

NAS Demand Predictions, Transportation Systems Analysis Model (TSAM) Compared with Other Forecasts

The current work incorporates the Transportation Systems Analysis Model (TSAM) to predict the future demand for airline travel. TSAM is a multi-mode, national model that predicts the demand for all long distance travel at a county level based upon population and demographics. The model conducts a mode choice analysis to compute the demand for commercial airline travel based upon the traveler s purpose of the trip, value of time, cost and time of the trip,. The county demand for airline travel is then aggregated (or distributed) to the airport level, and the enplanement demand at commercial airports is modeled. With the growth in flight demand, and utilizing current airline flight schedules, the Fratar algorithm is used to develop future flight schedules in the NAS. The projected flights can then be flown through air transportation simulators to quantify the ability of the NAS to meet future demand. A major strength of the TSAM analysis is that scenario planning can be conducted to quantify capacity requirements at individual airports, based upon different future scenarios. Different demographic scenarios can be analyzed to model the demand sensitivity to them. Also, it is fairly well know, but not well modeled at the airport level, that the demand for travel is highly dependent on the cost of travel, or the fare yield of the airline industry. The FAA projects the fare yield (in constant year dollars) to keep decreasing into the future. The magnitude and/or direction of these projections can be suspect in light of the general lack of airline profits and the large rises in airline fuel cost. Also, changes in travel time and convenience have an influence on the demand for air travel, especially for business travel. Future planners cannot easily conduct sensitivity studies of future demand with the FAA TAF data, nor with the Boeing or Airbus projections. In TSAM many factors can be parameterized and various demand sensitivities can be predicted for future travel. These resulting demand scenarios can be incorporated into future flight schedules, therefore providing a quantifiable demand for flights in the NAS for a range of futures. In addition, new future airline business scenarios are investigated that illustrate when direct flights can replace connecting flights and larger aircraft can be substituted, only when justified by demand.

Viken, Jeff↗
DOE OSTI2024

Comparative study of decay heat calculations with FLUKA and MCNP/CINDER2008

In designing a high-power tungsten target, decay heat driven temperature rise in the spallation volume is a safety concern during maintenance and in loss of coolant accidents. As tungsten hydrates and becomes volatile in steam at above 800 ◦C, it is important to keep the target temperature below this threshold when active cooling is unavailable. Decay heat in a target is calculated with particle transport simulation codes combined with transmutation codes. The calculated decay heat usually differs depending on the nuclear cross sections and the decay particle transport models built in the code architecture. In this paper, we calculated decay heat of a water-cooled tungsten target using popular particle transport codes, FLUKA and MCNP6® paired with CINDER2008 and compared the results. The target-moderator-reflector (TMR) system is modeled with a water-cooled solid tungsten target, water premoderators, liquid hydrogen cold moderators and beryllium reflectors. Water-cooled stainless-steel shielding is modeled around the TMR system. The tungsten volume is clad with a thin layer of erosion/corrosion resistant material. This study provides information about the uncertainty range in decay heat prediction of high-power spallation targets for hazard analysis.

Lee, Yong Joong↗
DOE OSTI2021

Bed form-induced hyporheic exchange and geochemical hotspots

Small-scale bed form topographies control hyporheic exchange and biogeochemical processes within aquatic sediments, which ultimately affect water quality and nutrient cycling at the watershed scale. The impact of three-dimensional and small-scale bed form topographies on hyporheic exchange and solute mixing is investigated in the present work. The effect of bed form morphologies on the development of zones of enhanced reaction rates (i.e., hotspots) is also studied. A computational fluid dynamics model to simulate river flow over bed forms is combined with a subsurface flow and multicomponent reactive solute transport model. A wide variety of bed form topographies are generated using geometric models by varying parameters controlling curvature as well as bed form wavelength and amplitude. The results in this research suggest that out-of-phase bed forms generate more complex hyporheic flow patterns which reduce the efficiency of solute transformations. Higher phase shifts in bed form shapes result in overall higher average velocity, larger zones of enhanced pressure and reaction rates, and higher amounts of solute exchange. Moreover, the bed form shapes control the reaction process for a wide range of sediment conductivities. Here this study advances the understanding of the effects of complex and small scale morphological features on hyporheic exchange processes including the rate and spatio-temporal distribution of reaction hotspots.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

MOOSE ProbML: Parallelized probabilistic machine learning and uncertainty quantification for computational energy applications

Here, this paper presents the development and demonstration of massively parallel probabilistic machine learning (ML) and uncertainty quantification (UQ) capabilities within the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open-source computational platform for parallel finite element and finite volume analyses. In addressing the computational expense and uncertainties inherent in complex multiphysics simulations, this paper integrates Gaussian process (GP) variants, active learning, Bayesian inverse UQ, adaptive forward UQ, Bayesian optimization, evolutionary optimization, and Markov chain Monte Carlo (MCMC) within MOOSE. It also elaborates on the interaction among key MOOSE systems — Sampler, MultiApp, Reporter, and Surrogate — in enabling these capabilities. The modularity offered by these systems enables development of a multitude of probabilistic ML and UQ algorithms in MOOSE. Example code demonstrations include parallel active learning and parallel Bayesian inference via active learning. The impact of these developments is illustrated through five applications relevant to computational energy applications: UQ of nuclear fuel fission product release, using parallel active learning Bayesian inference; very rare events analysis in nuclear microreactors using active learning; advanced manufacturing process modeling using multi-output GPs (MOGPs) and dimensionality reduction; fluid flow using deep GPs (DGPs); and tritium transport model parameter optimization for fusion energy, using batch Bayesian optimization. These capabilities are part of the MOOSE framework.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

How are Heterogeneous Nucleation Rate Observations Influenced by Instrument Resolution?

Experimental measurements of the heterogeneous nucleation rate rely on counting the number of nuclei with time. However, the size of a thermodynamically stable nucleus is often a few nanometers in diameter and is below the resolution of most (in situ) measurement techniques that provide a statistically valid sample. Due to the finite resolution of the instruments and analysis methods, it is challenging to capture the incipient nuclei and the subsequent evolution of nuclei density over time. In this work, we demonstrate the impact of instrument resolution on observed nuclei densities by comparing numerical modeling with experimental results. Further, to achieve this, we implemented heterogeneous nucleation within the pore-scale reactive transport modeling framework using classical nucleation theory (CNT). We compared the modeling results with nucleation rates measured using X-ray nanotomography (XnT) and evaluated how these impact the apparent values of the prefactor and interfacial energy based on CNT and the crystal growth rate. Specifically, we applied a resolution threshold (artificial resolution limit) in the model during nuclei counting to resemble an experimental resolution, ranging from 15 to 500 nm. The findings reveal that the instrument resolution significantly impacts the apparent prefactor and interfacial energy. Both apparent prefactor and interfacial energy decrease with a decrease in the instrument resolution. While deviation in the prefactor due to resolution is anticipated, those in the interfacial energy are unexpected. The approach described here allows one to correct apparent nucleation rates that depend on the instrument’s resolution to derive “intrinsic” CNT parameters for the prefactor and interfacial energy.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Genomes to ecosystem function: Targeting critical knowledge gaps in methanogenesis and translation to updated global biogeochemical models

Natural freshwater temperate wetland systems currently represent the largest natural source of atmospheric methane, but are relatively understudied using systems biology tools (e.g. meta-omics) compared to other high producing methane systems (e.g. peat, tropical, or reconstructed wetlands). Using field investigations at the NOAA operated sentinel site on Lake Erie, methane producing activities and responses to geochemical conditions will be determined along seasonal and spatial gradients (Objective 1). Here, a combination of high-throughput activity and gas measurements, combined with high-resolution systems biology and analytical methods, will provide in depth knowledge of the microbiological, chemical, and physical constraints on methane production in wetlands. Using laboratory microcosms, the formation of anoxic microsites and their capacity to facilitate methane production in wetland soils will be simulated (Objective 2). This objective will validate the findings from the field investigations, offering a more controlled environment for teasing out the role of different yet interrelated variables. Lastly, these field and laboratory data will be used for multi-scale, process-level evaluation of an ecosystem biogeochemical model that accommodates these newly identified processes and parameterizes representation of these processes along relevant environmental gradients (Objective 3).This research will identify multiple interacting geochemical, ecological, and metabolic constraints that are poorly understood, oversimplified, or missing in global biogeochemical methane models. This proposal targets the role of oxygen limitation on methane processes in soil domains, to improve reactive transport models of microbial carbon cycling across terrestrial-aquatic soils and generate data on nutrient cycling activities in Great Lake wetlands. This information could provide new insights into the microbial controllers of Lake Erie eutrophication.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2020

Aircraft System Identification from Multisine Inputs and Frequency Responses

A frequency-domain approach is described for estimating parameters, such as stability and control derivatives, in aircraft flight dynamic models from measured input and output data. The approach uses orthogonal phase-optimized multisines for moving the aircraft control effectors, Fourier analysis for computing multiple-input multiple-output frequency responses, and a maximum likelihood estimator called frequency response error (FRE) for determining values and uncertainties for the model parameters. The approach is demonstrated using flight test data for two subscale airplanes: the T-2 generic transport model and the X-56A aeroelastic demonstrator. Results and comparisons with the output-error method indicated that the approach produced accurate estimates of stability and control derivatives and their uncertainties from flight test data.

Grauer, Jared A.↗
DOE OSTI2022

Impact of Concurrent Solubilization and Fines Migration on Fracture Aperture Growth in Shales during Acidized Brine Injection

Acidic hydraulic fracturing fluid chemically and physically alters shale rock fabric during injection and shut-in, creating a “reaction-altered zone” along the fracture faces. To better characterize the variable thickness and composition of this reaction-altered zone under advective flow, we take a coupled experimental and modeling approach. A fluidic cell, with six fiducial markers, is first fabricated to keep the rock sample in place during the core floods and to allow image alignment of acquired images. Then, we conduct a series of reactive core floods in a clay-rich siliceous Wolfcamp shale sample with 10 wt % carbonate, using a synthetic fracturing fluid under no confining stress and at room temperature. High-resolution computed tomography (CT) scans are periodically conducted to observe the spatial alteration of the fracture network. We then perform scanning electron microscopy (SEM-EDS) on the two orthogonal surfaces (fracture surface and freshly cut profile face) to generate high-resolution elemental maps that show the change in mineralogy, both with distance along a given flow path along a fracture surface and with depth from the fracture surface into the shale matrix. These results are contrasted against a two-dimensional (2D) advection-diffusion reaction model developed previously for batch reactions between shale and synthetic fracturing fluids. The model simulates the geochemical interaction occurring at the fracture/matrix interface and penetrating into the shale matrix during the reactive core flood. Both model and experimental results show that the acidic brine is neutralized during the core flood, corresponding to an increase in fracture aperture as a function of fluid volume injected with the greatest change near the inlet. SEM-EDS scans reveal significant dissolution of carbonates on the fracture surface without pyrite oxidation. The reactive transport model indicates that carbonate depletion into the shale interior should be observable, yet SEM-EDS shows no discernible loss of carbonate in the orthogonal profile face. The combination of these observations suggests an additional fracture evolution mechanism in the reactive system, i.e., fines migration. We show that fines migration enhances the access of fracturing fluid to the matrix resulting in a more pronounced fracture widening. Finally, we conclude that coupled mineral dissolution and fines migration govern fracture aperture growth during acidized brine injection. In this work, we effectively show the underlying risk of relying solely on models that do not include an important (transport) process that can alter the system significantly and propose a combined chemomechanical mechanism for fracture evolution appropriate for this shale mineralogy.

58 GEOSCIENCES↗
DOE OSTI2020

Mapping the Concentration Distribution of Contaminant Plumes to the Computational Grid of the Plateau to River Model (Ver. 8.3)

This environmental calculation file (ECF) outlines the process conducted for defining the initial condition groundwater plumes for select contaminants of interest (COIs) to the Plateau River model version 8.3 (P2Rv8.3) domain (CP-57037, Model Package Report: Plateau to River Groundwater Model Version 8.3). Various data sources were combined into a single set of files defining plume conditions at time 1/1/2018. The resulting files are intended for use as inputs for transport modeling to support forecasts of the Central Plateau contaminant plume fate and transport.

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

Time-dependent water vapor desorption isotherm model of hardened cement paste

Highlights: • Short-term and long-term desorption isotherms of hardened cement paste are reviewed. • Water content at RH depends on both thermodynamic and microstructural equilibria. • An extended GAB equation is proposed to fit the cavitation drop during desorption. • A dynamic time-dependent sorption isotherm model is shown and used for simulations. • Numerical results capture anomalous drying in the same RH as seen experimentally. This study aims to enhance moisture transport modeling by elucidating the so-called anomalous water diffusion in cementitious materials. Water desorption isotherms are measured from samples at various drying stages using techniques with different durations to quantify the impact of the drying duration on the water sorption ability. A single water sorption isotherm does not solely give a relationship between the water content and relative humidity but also represents the state of microstructure. The continuous evolution of a desorption isotherm due to drying-induced microstructural rearrangement is demonstrated. For numerical modeling, the microstructural alteration can be explicitly considered through a dynamic desorption isotherm model, which governs the local thermodynamic equilibrium at the capillary meniscus. This approach is implemented into a multiphase transport code, whose ability to predict drying is validated using literature data. Finally, the effect of prolonged drying on the colloidal nature of the calcium-silicate-hydrate gel is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Global Impacts of Marine Methanethiol Emissions and Chemistry in the Atmosphere

Oceanic emissions of dimethyl sulfide (DMS) have long been known to influence aerosol particle composition, cloud condensation nuclei (CCN) concentration, and Earth’s radiative budget. However, the impact of oceanic emissions of methanethiol (MeSH), a sulfur compound produced by the same oceanic precursor as DMS, has been relatively less explored. The gas-phase oxidation of MeSH has a higher effective yield of SO 2 and a shorter oxidative lifetime compared to DMS, highlighting the relevance of this pathway for the modeled representation of particle formation, growth, and CCN abundance in the marine atmosphere. Here, we use the global chemical transport model GEOS-Chem to explore possible scenarios representative of specific environmental conditions and MeSH emission schemes based on previous experimental studies. We further implement and test previously reported chemical mechanisms for MeSH oxidation, along with additional improvements, highlighting key uncertainties and sensitivities for regional and global sulfur budgets. We place our results in the context of recent modeling updates to DMS chemistry and cloud processing, which further impact SO 2 production in the marine atmosphere in parallel with MeSH oxidation. Within the overall marine sulfur budget, our findings highlight that MeSH plays a significant role in SO 2 production in the marine atmosphere, contributing to regional surface layer concentration increases of up to 40–60%. These results point to the importance of MeSH for efforts aimed at improving the modeled representation of sulfur spatiotemporal patterns relevant to air quality predictions and climate impact assessments.

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

Deconvolving the roles of E × B shear and pedestal structure in the energy confinement quality of super H-mode experiements

Analysis of “super H-mode” experiments on DIII-D has put forward that high plasma toroidal rotation, not high pedestal, plays the essential role in achieving energy confinement quality H 98y2 >>1. Recently, super H-mode experiments with variable input torque have confirmed that high rotation shear discharges have very high levels of H 98y2 (>1.5), independent of the pedestal height, and that high pedestal discharges with low rotation shear have levels of H 98y2 only slightly above 1 (≤1.2). Although some increase in stored energy with higher pedestal occurs, the energy confinement quality mainly depends on the toroidal rotation shear, which varies according to different levels of injected neutral beam torque per particle. Here, quasi-linear gyrofluid modeling achieves a good match of the experiment when including the E×B shear; without including plasma rotation, the modeling predicts a confinement quality consistent with the empirical observation of H 98y2 ~1.2 at low rotation. Nonlinear gyrokinetic transport modeling shows that the effect of E×B turbulence stabilization is far larger than other mechanisms, such as the so-called hot-ion stabilization (T i /T e ) effect. Consistent with these experimental and modeling results are previous simulations of the ITER Baseline Scenario using a super H-mode pedestal solution, which showed the potential to exceed the Q=10 target if the pedestal density could be increased above the Greenwald limit. A close look at these simulations reveals that the predicted energy confinement quality is below 1 even at the highest pedestal pressure. The improvement in Q at higher pedestal density is due to the improved fusion power generation at the higher core density associated with higher pedestal density, not to an improved energy confinement quality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1988

Modeling magnetospheric plasma; Proceedings of the First Huntsville Workshop on Magnetosphere/Ionosphere Plasma Models, Guntersville, AL, Oct. 14-16, 1987

The conference presents papers on the global modeling of magnetospheric plasma processes, the modeling of the midlatitude ionosphere and plasmasphere, the modeling of the auroral zone and boundary layer, the modeling of the polar magnetosphere and ionosphere, and the modeling of the plasma sheet and ring current. Particular attention is given to the kinetic approach in magnetospheric plasma transport modeling, self-consistent neutral point current and fields from single particle dynamics, preliminary statistical survey of plasmaspheric ion properties from observations by DE 1/RIMS, and a model of auroral potential structures based on dynamics explorer plasma data. Other topics include internal shear layers in auroral dynamics, quantitative parameterization of energetic ionospheric ion outflow, and open flux merging in an expanding polarcap model.

Moore, T. E.↗
NASA NTRS1994

Kinetic modeling of O(+) upflows resulting from E x B convection heating in the high-latitude F region ionosphere

We report here the results of modeling work aimed at understanding the development of ionospheric O(+) field-aligned upflows that develop in response to high-latitude E x B drift induced frictional heating. The model used is a collisional semikinetic model which includes ion-neutral resonant charge exchange and polarization collisions as well as Coulomb self-collisions. It also includes the process of chemical removal of O(+) as well as all of the macroscopic forces: ambipolar electric, gravity, magnetic mirror, and centripetal. Model results show the development of several types of non-Maxwellian velocity distributions including toroids at low altitude, distributions with large heat flow in the perpendicular component at intermediate altitudes, and distributions with a separate upflowing population or upward superthermal tail at high altitudes. Whenever the convection electric field increases from a small value (less than 25 mV/m) to a large value (100-200 mV/m) in 6 min or less large upflows develop with parallel drift speeds which peak (below 1000 km) at values between 500 m/s and 2 km/s, parallel fluxes which peak between 6.0 x 10(exp 8) and 3.2 x 10(exp 9)/sq cm/s, and parallel per particle heat flows which peak between 8.0 x 10(exp -9) and 8.0 x 10(exp -8) ergs cm/s. The higher values in these ranges occur for a cooler neutral atmosphere, with a larger convection electric field that is turned on quickly. The model produces field-aligned O(+) flow speeds that are larger than those produced by a 20-moment generalized transport model but smaller then those produced by an isotropic hydrodynamic model for comparable values of the convection turn on times. The model results compare favorably with some topside satellite and radar data.

Wilson, G. R.↗
NASA NTRS2000

The Interaction Between Dynamics and Chemistry of Ozone in the Set-Up Phase of the Northern Hemisphere Polar Vortex

Understanding stratospheric ozone loss to the point of accurately predicting ozone in the future requires correctly distinguishing chemical from transport-induced changes in ozone. For example, evaluating the impact of chlorine reduction in controlling stratospheric ozone loss requires estimating the amount of ozone lost in chemical reactions in the polar winter and spring. The Northern Hemisphere winter polar region is a particularly crucial and interesting area because it appears that the Northern vortex may currently be poised near the threshold of extreme ozone destruction such as that which now occurs annually in the Antarctic ozone "hole." In this presentation we explore the interaction of ozone transport and chemistry through the Northern late summer and fall seasons as the vortex circulation becomes established. This phase of the seasonal cycle determines the starting point for heterogeneous processes and chlorine-driven loss that take control in the winter vortex. Using a combination of profile data from POAM, HALOE, and in situ measurements, we show that relatively low ozone at high latitudes in the middle stratosphere is associated with vortex airmasses and that these ozone abundances evolve photochemically from characteristically higher values at the end of the summer. The zonal variance of ozone mixing ratio also increases greatly at this time consistent with increasing wave-driven transport. Comparison with a three-dimensional chemistry-transport model is used to generalize the findings from the limited set of observations and quantify the relative roles of transport and chemistry in determining the ozone mixing ratio distributions.

Kawa, S.R.↗
NASA NTRS2007

Initial Verification of GEOS-4 Aerosols Using CALIPSO and MODIS: Scene Classification

A-train sensors such as MODIS and MISR provide column aerosol properties, and in the process a means of estimating aerosol type (e.g. smoke vs. dust). Correct classification of aerosol type is important because retrievals are often dependent upon selection of the right aerosol model. In addition, aerosol scene classification helps place the retrieved products in context for comparisons and analysis with aerosol transport models. The recent addition of CALIPSO to the A-train now provides a means of classifying aerosol distribution with altitude. CALIPSO level 1 products include profiles of attenuated backscatter at 532 and 1064 nm, and depolarization at 532 nm. Backscatter intensity, wavelength ratio, and depolarization provide information on the vertical profile of aerosol concentration, size, and shape. Thus similar estimates of aerosol type using MODIS or MISR are possible with CALIPSO, and the combination of data from all sensors provides a means of 3D aerosol scene classification. The NASA Goddard Earth Observing System general circulation model and data assimilation system (GEOS-4) provides global 3D aerosol mass for sulfate, sea salt, dust, and black and organic carbon. A GEOS-4 aerosol scene classification algorithm has been developed to provide estimates of aerosol mixtures along the flight track for NASA's Geoscience Laser Altimeter System (GLAS) satellite lidar. GLAS launched in 2003 and did not have the benefit of depolarization measurements or other sensors from the A-train. Aerosol typing from GLAS data alone was not possible, and the GEOS-4 aerosol classifier has been used to identify aerosol type and improve the retrieval of GLAS products. Here we compare 3D aerosol scene classification using CALIPSO and MODIS with the GEOS-4 aerosol classifier. Dust, smoke, and pollution examples will be discussed in the context of providing an initial verification of the 3D GEOS-4 aerosol products. Prior model verification has only been attempted with surface mass comparisons and column optical depth from AERONET and MODIS.

Welton, Ellsworth J.↗
NASA NTRS2017

Daily Ambient Air Pollution Metrics for Five Cities: Evaluation of Data Fusion-Based Estimates and Uncertainties

Spatiotemporal characterization of ambient air pollutant concentrations is increasingly relying on the combination of observations and air quality models to provide well-constrained, spatially and temporally complete pollutant concentration fields. Air quality models, in particular, are attractive, as they characterize the emissions, meteorological, and physiochemical process linkages explicitly while providing continuous spatial structure. However, such modeling is computationally intensive and has biases. The limitations of spatially sparse and temporally incomplete observations can be overcome by blending the data with estimates from a physically and chemically coherent model, driven by emissions and meteorological inputs. We recently developed a data fusion method that blends ambient ground observations and chemical transport-modeled (CTM) data to estimate daily, spatially resolved pollutant concentrations and associatedcorrelations. In this study, we assess the ability of the data fusion method to produce daily metrics (i.e., 1-hr max, 8-hr max, and 24-hr average) of ambient air pollution that capture spatiotemporal air pollution trends for 12 pollutants (CO, NO2, NOx, O3, SO2, PM (sub10), PM (sub 2.5), and five PM (sub 2.5) components) across five metropolitan areas (Atlanta, Birmingham, Dallas, Pittsburgh, and St. Louis), from 2002 to 2008. Three sets of comparisons are performed: (1) the CTM concentrations are evaluated for each pollutant and metropolitan domain, (2) the data fusion concentrations are compared with the monitor data, (3) a comprehensive cross-validation analysis against observed data evaluates the quality of the data fusion model simulations across multiple metropolitan domains. The resulting daily spatial field estimates of air pollutant concentrations and uncertainties are not only consistent with observations, emissions, andmeteorology, but substantially improve CTM-derived results for nearly all pollutants and all cities, with the exception of NO2 for Birmingham. The greatest improvements occur for O3 and PM (sub 2.5). Squared spatiotemporal correlation coefficients range between simulations and observations determined using cross-validation across all cities for air pollutants of secondary and mixed origins are R-squared equal to 0.88-0.93 (O3), 0.81-0.89 (SO4), 0.67-0.83 (PM (sub 2.5)), 0.52-0.72 (NO3), 0.43-0.80 (NH4), 0.32-0.51 (OC), and 0.14-0.71 (PM (sub 10)). Results for relatively homogeneous pollutants of secondary origin, tend to be better than those for more spatially heterogeneous (larger spatial gradients) pollutants of primary origin (NOx, CO, SO2 and EC). Generally, background concentrations and spatial concentration gradients reflect interurban airshed complexity and the effects of regional transport, whereas daily spatial pattern variability shows intra urban consistency in the fused data. With sufficiently high CTM spatial resolution, traffic-related pollutants exhibit gradual concentration gradients that peak toward the urban centers. Ambient pollutant concentration uncertainty estimates for the fused data are both more accurate and smaller than those for either the observations or the model simulations alone.

spatiotemporal fusion↗
DOE OSTI2020

Bridging Scales to Model Reactive Diffusive Transport in Porous Media

Published by ECS. Two novel scale-bridging algorithms to model reaction-diffusion transport in porous media are presented. The algorithms are based on direct numerical simulations and couple the information of a micro-scale model, which accounts for the large field of view provided by micro X-ray computed tomography (X-ray CT), and a nano-scale model, which locally resolves transport in the fine structure extracted from nano X-ray CT. The micro-scale model is discretized in the through-plane direction into a 1D grid, where effective properties and internal boundaries are determined based on the results from the nano-scale model. The validated algorithms are used to examine transport of oxygen in precious group metal-free electrodes considering both zero- and first-order kinetics. Unlike conventional methods, the results show that the effective diffusivity is not a passive property but increases in regions where the reaction-rate coefficient is large. The proposed algorithms account for the multiscale coupling of reaction-diffusion transport and material microstructure, thus improving the predictions compared to conventional methods.

42 ENGINEERING↗