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 919 records · Page 51

Combining Machine Learning and Numerical Simulation for High-Resolution PM2.5 Concentration Forecast

Forecasting ambient PM2.5 concentrations with spatiotemporal coverage is key to alerting decision-makers of pollution episodes and preventing detrimental public exposure, especially in regions with limited ground air monitoring stations. The existing methods either rely on chemical transport models (CTMs) to forecast spatial distribution of PM2.5 with nontrivial uncertainty or statistical algorithms to forecast PM2.5 concentration time-series at air monitoring locations without continuous spatial coverage. In this study, we developed a PM2.5 forecast framework by combining the robust Random Forest algorithm with a publicly accessible global CTM forecast product – NASA’s Goddard Earth Observing System “Composition Forecasting” (GEOS-CF), providing spatiotemporally continuous PM2.5 concentration forecasts for the next five days at a 1-km spatial resolution. Our forecast experiment was conducted for a region in Central China including the populous and polluted Fenwei Plain. The forecast for the next two days had overall validation R2 of 0.76 and 0.64, respectively; the R2 was around 0.5 for the following three forecast days. Spatial cross-validation showed similar validation metrics. Our forecast model, with validation normalized mean bias close to zero, substantially reduced the large biases in GEOS-CF. The proposed framework requires minimal computational resources compared to running CTMs at urban scales, enabling near-real-time PM2.5 forecast in resource-restricted environments.

PM2.5↗

Transport of Mars atmospheric water into high northern latitudes during a polar warming

Several numerical experiments were conducted with a simplified tracer transport model in order to attempt to examine the poleward transport of Mars atmospheric water during a polar warming like that which occurred during the winter solstice dust storm of 1977. The flow for the transport experiments was taken from numerical simulations with a nonlinear beta-plane dynamical model. Previous studies with this model have demonstrated that a polar warming having essential characteristics like those observed during the 1977 dust storm can be produced by a planetary wave mechanism analogous to that responsible for terrestrial sudden stratospheric warmings. Several numerical experiments intended to simulate water transport in the absence of any condensation were carried out. These experiments indicate that the flow during a polar warming can transport very substantial amounts of water to high northern latitudes, given that the water does not condense and fall out before reaching the polar region.

Barnes, J. R.↗

A Liquid–Vapor Two-Phase Model of Direct Methanol Fuel Cells With Platinum Group Metal-Free Cathode Catalyst

Abstract This study has developed a two-dimensional, two-phase transport model to investigate the transport characteristics in direct methanol fuel cells (DMFCs) using platinum group metal (PGM)-free cathode catalysts. The model considered anisotropic properties of the gas diffusion layer (GDL) caused by current collector’s mechanical compression, the interfacial mass transfer of water and methanol between liquid and vapor, and unique properties of the cathode PGM-free catalyst layer. Results showed that the liquid methanol solution from the anode could provide sufficient water to hydrate the proton exchange membrane (PEM), and the relative humidity of the cathode air did not impact the membrane hydration. Fully hydrating the cathode air may deteriorate the fuel cell performance, especially when the operating temperature is close to 100 °C because the exponential increase of the saturated water pressure with temperature decreased the partial pressure of oxygen. The optimized operating temperature increased with the increase of air pressure and was about 80 °C at 1.5 atm cathode pressure. To achieve the US Department of Energy’s performance target of 300 mW/cm2 peak power density, catalytic activities of both the anode and cathode catalysts need to be improved by one order of magnitude compared with the state-of-the-art commercial catalysts.

Electrochemistry↗

Global Simulations of Phase State and Equilibration Time Scales of Secondary Organic Aerosols with GEOS-Chem

The phase state of secondary organic aerosols (SOA) can range from liquid through amorphous semisolid to glassy solid, which is important to consider as it influences various multiphase processes including SOA formation and partitioning, multiphase chemistry, and cloud activation. In this study, we simulate the glass transition temperature and viscosity of SOA over the globe using the global chemical transport model, GEOS-Chem. The simulated spatial distributions show that SOA at the surface exist as liquid over equatorial regions and oceans, semisolid in the midlatitude continental regions, and glassy solid over lands with low relative humidity. The predicted SOA viscosities are mostly consistent with the available measurements. In the free troposphere, SOA particles are mostly predicted to be semisolid at 850 hPa and glassy solid at 500 hPa, except over tropical regions including Amazonia, where SOA are predicted to be low viscous. Phase state also exhibits seasonal variation with a higher frequency of semisolid and solid particles in winter compared to warmer seasons. We calculate equilibration time scales of SOA partitioning (τ eq ) and effective mass accommodation coefficient (α eff ), indicating that τ eq is shorter than the chemical time step of GEOS-Chem of 20 min and αeff is close to unity for most locations at the surface level, supporting the application of equilibrium SOA partitioning. However, τ eq is prolonged and α eff is lowered over drylands and most regions in the upper troposphere, suggesting that kinetically limited growth would need to be considered for these regions in future large-scale model studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Characterizing the combined impact of nucleation-driven precipitation and secondary passivation on carbon mineralization

The evolution of mineral reactive surface area is one of the primary phenomena controlling the progression and extent of mineral carbonation. The CO 2 mineralization begins with nucleation of crystals that provide initial surface area for subsequent growth of the mineral. However, many reactive transport models (RTMs) for CO 2 mineralization do not include the nucleation process. The few RTMs that do include it are yet to be validated against experimental data. Similarly, many RTMs ignore passivating effects of the secondary mineral, which coats the surface of the dissolving mineral, slow down the reaction process, and reduce the total extent of carbonation. Furthermore, the combined impact of nucleation and passivation on carbon mineralization is yet to be properly characterized. In this study, we consider the coupled effects of passivation and nucleation on the mineralization extent. The nucleation-driven precipitation model relies on the formation of nuclei to provide a surface area for crystal growth, while a new model is proposed to account for passivation effects. Our analysis shows that (i) omission of nucleation leads to overestimation of extent of mineralization, and (ii) omission of passivation leads to overestimation of host rock reactivity. The model was evaluated via comparison with CO 2 mineralization data from the literature and models that ignore these processes. We observed that including nucleation and passivation lead to closer predictions of the CO 2 mineralization extent. Therefore, this study highlights the importance of including the coupled nucleation-driven precipitation and secondary passivation in RTMs. Here, the findings from the study can be applied in various scientific and engineering applications such as petroleum production, cement carbonation, CO 2 sequestration, chemical weathering, and concrete degradation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Increasing aggregate size reduces single-cell organic carbon incorporation by hydrogel-embedded wetland microbes

Abstract Microbial degradation of organic carbon in sediments is impacted by the availability of oxygen and substrates for growth. To better understand how particle size and redox zonation impact microbial organic carbon incorporation, techniques that maintain spatial information are necessary to quantify elemental cycling at the microscale. In this study, we produced hydrogel microspheres of various diameters (100, 250, and 500 μm) and inoculated them with an aerobic heterotrophic bacterium isolated from a freshwater wetland (Flavobacterium sp.), and in a second experiment with a microbial community from an urban lacustrine wetland. The hydrogel-embedded microbial populations were incubated with 13C-labeled substrates to quantify organic carbon incorporation into biomass via nanoSIMS. Additionally, luminescent nanosensors enabled spatially explicit measurements of oxygen concentrations inside the microspheres. The experimental data were then incorporated into a reactive-transport model to project long-term steady-state conditions. Smaller (100 μm) particles exhibited the highest microbial cell-specific growth per volume, but also showed higher absolute activity near the surface compared to the larger particles (250 and 500 μm). The experimental results and computational models demonstrate that organic carbon availability was not high enough to allow steep oxygen gradients and as a result, all particle sizes remained well-oxygenated. Our study provides a foundational framework for future studies investigating spatially dependent microbial activity in aggregates using isotopically labeled substrates to quantify growth.

59 BASIC BIOLOGICAL SCIENCES↗

Impacts of Intercontinental Transport of Anthropogenic Fine Particulate Matter on Human Mortality

Fine particulate matter with diameter of 2.5 microns or less (PM2.5) is associated with premature mortality and can travel long distances, impacting air quality and health on intercontinental scales. We estimate the mortality impacts of 20 % anthropogenic primary PM2.5 and PM2.5 precursor emission reductions in each of four major industrial regions (North America, Europe, East Asia, and South Asia) using an ensemble of global chemical transport model simulations coordinated by the Task Force on Hemispheric Transport of Air Pollution and epidemiologically-derived concentration-response functions. We estimate that while 93-97 % of avoided deaths from reducing emissions in all four regions occur within the source region, 3-7 % (11,500; 95 % confidence interval, 8,800-14,200) occur outside the source region from concentrations transported between continents. Approximately 17 and 13 % of global deaths avoided by reducing North America and Europe emissions occur extraregionally, owing to large downwind populations, compared with 4 and 2 % for South and East Asia. The coarse resolution global models used here may underestimate intraregional health benefits occurring on local scales, affecting these relative contributions of extraregional versus intraregional health benefits. Compared with a previous study of 20 % ozone precursor emission reductions, we find that despite greater transport efficiency for ozone, absolute mortality impacts of intercontinental PM2.5 transport are comparable or greater for neighboring source-receptor pairs, due to the stronger effect of PM2.5 on mortality. However, uncertainties in modeling and concentration-response relationships are large for both estimates.

human health↗

Development of Scalable Reactive Transport Framework for PFAS

Research on PFAS chemicals is extensive and covers all aspects, including analytical quantification, determination of properties, toxicology, ex situ treatment, and in situ remediation. PFAS chemicals are very stable and thus persistent/recalcitrant in the environment. Although there are many unknowns about PFAS chemicals, degradation pathways, reaction rates, etc., many different sorption, oxidation, reduction, biological, and innovative treatment approaches are being developed. Sorption with activated carbon is currently the only fully available in situ treatment technology for PFAS-impacted groundwater. Given the wide array of PFAS chemicals and transformation products, remediation may need multi-step treatment trains to fully address the PFAS contamination. The work here provides kinetic reaction modules that represent an initial cut at functionality representing PFAS migration and reaction in groundwater aquifer flow and transport models. One reaction kinetics module provides a method to model kinetically limited adsorption using a mass transfer model. The second reaction module represents biological transformation of 8:2 FTOH and daughter species, illustrating how a complex reaction pathway network can be represented. Both reaction modules allow for spatially variable parameter values so that a variety of remediation approaches (e.g., a PRB or volumetric treatment or variations in geochemical conditions) can be simulated. The intent with these PFAS reaction modules is to provide tools for practitioners to aid in the selection, design, and assessment of potential in situ PFAS remediation strategies. It is anticipated that, as PFAS remediation technologies and scientific understanding advances, these modules would be refined or replaced to match new knowledge.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Geochemical Assessment for Carbon Sequestration in the Conasauga Group, Northwest Georgia, USA

Sedimentary geological formations are known to be great candidates for geological carbon sequestration. Published studies suggest the southeast of the United States contains many formations suitable for carbon storage. The Cassville 1 Stratigraphic Borehole well could act as a potential carbon reservoir for nearby energy resource facilities in Georgia, United States. Although studies have shown that porous formations are adequate for geological carbon sequestration, it is important to understand possible geochemical reactions between CO 2 and the targeted geological formation before injecting any fluids. In this study, a sandstone sample from the Cassville 1 well is being considered for geological carbon sequestration in the Conasauga Group in Northwest Georgia. Here, the collected sandstone sample, consisting of quartz, K-feldspar, micas, kaolinite, and carbonate minerals such as calcite and dolomite, has a 6% porosity. Leveraging the formation composition and porosity, a one-dimensional continuum reactive transport model was built using CrunchFlow to assess possible geochemical reactions between injected CO 2 and the geological formation. Simulation results show that the carbonate minerals, calcite and dolomite, dissolve during the injection period of 10,000 days, increasing formation porosity from 6% to as much as 30%. The rate and extent of carbonate mineral dissolution and resulting porosity increase are highly sensitive to mineral reactive surface area values. No evidence of mineral precipitation was observed, suggesting that dissolution reactions will control porosity evolution during the CO 2 injection period.

42 ENGINEERING↗

Bayesian inference of in-medium baryon-baryon scattering cross sections from HADES proton flow data

Within a Bayesian statistical framework using a Gaussian Process emulator for an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model simulator of heavy-ion reactions at intermediate energies, we infer from the HADES proton flow data the posterior probability distribution functions of in-medium baryon-baryon scattering cross section modification factor X with respect to free-space and the corresponding incompressibility K of nuclear matter as well as their correlation function. In conclusion, the mean value of X is found to be $X$ = $1.32_ {+0.28} ^ {-0.40} $ at 68% confidence level assuming the nuclear incompressibility K will not exceed 400 MeV, providing circumstantial evidence for enhanced baryon-braryon scattering cross sections in hot and dense nuclear matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Principal Component Analysis of azimuthal flow in intermediate-energy heavy-ion reactions

Principal Component Analysis (PCA) via Singular Value Decomposition (SVD) of large datasets is an adaptive exploratory method to uncover natural patterns underlying the data. Several recent applications of the PCA-SVD to event-by-event single-particle azimuthal angle distribution matrices in ultra-relativistic heavy-ion collisions at RHIC-LHC energies indicate that the sine and cosine functions chosen a priori in the traditional Fourier analysis are naturally the most optimal basis for azimuthal flow studies according to the data itself. We perform PCA-SVD analyses of mid-central Au+Au collisions at $E$ $beam$ / $A$ =1.23 GeV simulated using an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model to address the following two questions: (1) if the principal components of the covariance matrix of nucleon azimuthal angle distributions in heavy-ion reactions around 1 GeV/nucleon are naturally sine and/or cosine functions and (2) what if any advantages the PCA-SVD may have over the traditional flow analysis using the Fourier expansion for studying the EOS of dense nuclear matter. In conclusion, we find that (1) in none of our analyses the principal components come out naturally as sine and/or cosine functions, (2) while both the eigenvectors and eigenvalues of the covariance matrix are appreciably EOS dependent, the PCA-SVD has no apparent advantage over the traditional Fourier analysis for studying the EOS of dense nuclear matter using the azimuthal collective flow in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Along the Channel Gradients Impact on the Spatioactivity of Gas Diffusion Electrodes at High Conversions during CO 2 Electroreduction

Results of a 2-D transport model for a gas diffusion electrode performing CO 2 reduction to CO with a flowing catholyte are presented, including the concentration gradients along the flow cell, spatial distribution of the current density and local pH in the catalyst layer. The model predicts that both the concentration of CO 2 and the buffer electrolyte gradually diminish along the channels for a parallel flow of gas and electrolyte as a result of electrochemical conversion and nonelectrochemical consumption. At high single-pass conversions, significant concentration gradients exist along the flow channels leading to large local variations in the current density (>150 mA/cm 2 ), which becomes prominent when compared to ohmic losses. In addition, concentration overpotentials change dramatically with CO 2 flow rate, which results in significant differences in outlet concentrations at high conversions. The outlet concentration of CO attains a maximum of 80% along with 5% CO 2 and 15% H 2 , although the maximum single-pass conversion is limited to below 60% due to homogeneous consumption by the electrolyte. Fundamental and practical implications of our findings on electrochemical CO 2 reduction are discussed with a focus on the trade-off between high current density operation and high single-pass conversion efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Photon detector response function methodology using MCNP and shift hybrid radiation transport code for wide-area contamination assay applications

Here, radiation transport modeling using the Monte Carlo N-Particle (MCNP) radiation transport code and Monte Carlo code, Shift, were employed to model detector responses for a variety of wide-area photon contamination scenarios. In this study, 2" × 2" and 3" × 3" cylindrical NaI(Tl) scintillation detector configurations at source detector-distances of 0.5 cm, 1 cm, 2.54 cm, 10 cm, and 30 cm were modeled. Media of soil, concrete, and steel were evaluated for contamination depths ranging from surface to a depth of an infinite thickness in each medium for photon energies ranging from 20 keV to 3 MeV, which correspond to the energies that current detectors can discern. Monoenergetic photon surface contamination detector responses for each of the media, source–detector distances, and detectors were estimated using MCNP v6.2. Shift was harnessed for improved variance reduction of particle transport in highly attenuating media to obtain average cell fluxes in the two MCNP NaI(Tl) scintillation detector configurations. Average cell flux values in Shift were coupled with detector responses from MCNP to convert average cell flux in a void to energy distribution of pulses in the NaI(Tl) scintillation detector crystal of interest. An optimized detector response function methodology was developed by coupling the MCNP radiation transport method with the Consistent Adjoint Driven Importance Sampling (CADIS) hybrid radiation transport method built into Shift to significantly decrease the runtime of thousands of MCNP pulse height simulations. The methodology may be utilized to quickly and accurately facilitate the assessment of a broad range of wide-area environmental contamination assay and decommissioning cleanup applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Deployment of portable, modular gas samplers as part of an atmospheric tracer experiment

Underground nuclear explosions release noble gases into the atmosphere that can be detected to support international monitoring efforts. Atmospheric transport models help predict the movement of these gases over long distances, but struggle to predict the movement in the atmosphere local to the release. A field experiment was designed to monitor the movement of 127 Xe within a 5-km radius. Four gas samplers were deployed as part of this experiment to collect atmospheric samples at various distances from the release point. In conclusion, these samples were then analyzed in a near-field lab using a NaI detector and in an off-site lab using gamma-gamma coincidence and beta-gamma coincidence counting.

Radioxenon↗

The impact of mineral reactive surface area variation on simulated mineral reactions and reaction rates

Reactive transport modeling is an essential tool to simulate complex geochemical reactions in porous media that can impact formation properties including porosity and permeability. However, simulating these reactions is challenging due to uncertainties in model parameters, particularly mineral surface areas. Imaging has emerged as a powerful means of estimating model parameters including porosity and mineral abundance, accessibility and accessible surface area. However, these parameters, particularly mineral accessible surface area, vary with image resolution. This work aims to enhance understanding of the impact of image resolution and other means of estimating mineral reactive surface area on simulated mineral reactions and reaction rates. Mineral surface areas calculated from images with resolutions of 0.34 μm and 5.71 μm were used to simulate mineral reactions in the context of geologic CO 2 sequestration in the Paluxy formation at the continuum scale. Additional simulations were carried out using BET surface areas collected from the literature and geometric surface areas. Simulations were run for 7300 days and mineral volume fractions and effluent ion concentrations tracked and compared. Variations in mineral surface areas measured from images are within 1 order of magnitude and yield similar simulation results, indicating the impact of image resolution on simulated reactions and reaction rates is minimum for the resolutions and sample considered. In comparison, surface areas obtained from BET and geometric approaches are 1–4 orders of magnitude higher than image-obtained surface areas and result in greater simulated reaction rates and extents. Minerals with high reaction rates (calcite and siderite) are most impacted by surface area values at short times where simulated mineral volume fractions at longer times agree relatively well, even for simulations with several orders of magnitude variation in surface area. Phases with lower reaction rates, such as K-feldspar and muscovite, are predominantly impacted over longer times where variations in surface areas impact reaction extents and porosity evolution. Overall, variations in surface areas due to image resolution are small and result in little variation in simulated reactions and reaction rates while there are significant variations in simulation results when other surface area estimates are used. These variations, however, depend on the reactivity of the mineral phase where surface areas of fast-reacting phases largely impact simulations at short (hours to days) timescales and surface areas of slower-reactive phases impact longer term simulations (years).

58 GEOSCIENCES↗

A New Formulation of the Filter-Error Method for Aerodynamic Parameter Estimation in Turbulence

A new formulation of the filter-error method for estimating aerodynamic parameters in nonlinear aircraft dynamic models during turbulence was developed and demonstrated. The approach uses an estimate of the measurement noise covariance to identify the model parameters, their uncertainties, and the process noise covariance, in a relaxation method analogous to the output-error method. Prior information on the model parameters and uncertainties can be supplied, and a post-estimation correction to the uncertainty was included to account for colored residuals not considered in the theory. No tuning parameters, needing adjustment by the analyst, are used in the estimation. The method was demonstrated in simulation using the NASA Generic Transport Model, then applied to the subscale T-2 jet-engine transport aircraft flight. Modeling results in different levels of turbulence were compared with results from time-domain output error and frequency- domain equation error methods to demonstrate the effectiveness of the approach.

Grauer, Jared A.↗

Ion and Water Dynamics in the Transition from Dry to Wet Conditions in Salt-Doped PEG

The influence of the water content on ion and water transport mechanisms in polymer membranes under low to moderate hydration conditions remains poorly understood. In this study, we combine ion and water diffusivity (PFG-NMR) measurements with atomistic molecular dynamics simulations to better understand transport processes in hydrated salt-doped poly(ethylene glycol). Above the water percolation threshold, the experimental and simulated diffusivities are in good agreement with the free volume transport models. At low hydration levels, unlike dry systems, ion diffusion cannot be described by polymer segmental dynamics alone. We rationalize such observations by the interplay between ion–water and ion–polymer solvation of cations and between ion–water and cation–anion interactions for anions. Further, we demonstrate that a two-state model combining ion–water solvation and free volume transport can describe water dynamics across the entire hydration range of interest. Here, our findings provide a more encompassing analysis of ion and water transport in hydrated polyelectrolytes, specifically in the low hydration regime.

36 MATERIALS SCIENCE↗

Historical total ozone radiative forcing derived from CMIP6 simulations

Radiative forcing (RF) time series for total ozone from 1850 up to the present day are calculated based on historical simulations of ozone from 10 climate models contributing to the Coupled Model Intercomparison Project Phase 6 (CMIP6). In addition, RF is calculated for ozone fields prepared as an input for CMIP6 models without chemistry schemes and from a chemical transport model simulation. A radiative kernel for ozone is constructed and used to derive the RF. The ozone RF in 2010 (2005–2014) relative to 1850 is 0.35 W m -2 [0.08–0.61] (5–95% uncertainty range) based on models with both tropospheric and stratospheric chemistry. One of these models has a negative present-day total ozone RF. Excluding this model, the present-day ozone RF increases to 0.39 W m -2 [0.27–0.51] (5–95% uncertainty range). The rest of the models have RF close to or stronger than the RF time series assessed by the Intergovernmental Panel on Climate Change in the fifth assessment report with the primary driver likely being the new precursor emissions used in CMIP6. The rapid adjustments beyond stratospheric temperature are estimated to be weak and thus the RF is a good measure of effective radiative forcing.

54 ENVIRONMENTAL SCIENCES↗