Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface mass variations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The Milky Way tomography with APOGEE: intrinsic density distribution and structure of mono-abundance populations

ABSTRACT The spatial distribution of mono-abundance populations (MAPs, selected in [Fe/H] and [Mg/Fe]) reflect the chemical and structural evolution in a galaxy and impose strong constraints on galaxy formation models. In this paper, we use APOGEE data to derive the intrinsic density distribution of MAPs in the Milky Way, after carefully considering the survey selection function. We find that a single exponential profile is not a sufficient description of the Milky Way’s disc. Both the individual MAPs and the integrated disc exhibit a broken radial density distribution; densities are relatively constant with radius in the inner Galaxy and rapidly decrease beyond the break radius. We fit the intrinsic density distribution as a function of radius and vertical height with a 2D density model that considers both a broken radial profile and radial variation of scale height (i.e. flaring). There is a large variety of structural parameters between different MAPs, indicative of strong structure evolution of the Milky Way. One surprising result is that high-α MAPs show the strongest flaring. The young, solar-abundance MAPs present the shortest scale height and least flaring, suggesting recent and ongoing star formation confined to the disc plane. Finally we derive the intrinsic density distribution and corresponding structural parameters of the chemically defined thin and thick discs. The chemical thick and thin discs have local surface mass densities of 5.62 ± 0.08 and 15.69 ± 0.32 M⊙pc−2, respectively, suggesting a massive thick disc with a local surface mass density ratio between thick to thin disc of 36 per cent.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Low-order aquifer thermal energy storage model for geothermal system simulation

This paper presents a low-order aquifer thermal energy storage (ATES) model for simulation of combined subsurface and above-surface energy systems. The model is included in the Modelica IBPSA Library, which is a free open-source library with basic models for building and district energy and control systems. The model uses a lumped-component method, in which the transient conductive-convective heat and mass transfer equation is radially discretized. To verify the accuracy of the model, we present an inter-model comparison from a simulation test suite. Results show that the Modelica ATES model is in good agreement, with a normalized mean bias error for yearly variation of aquifer temperatures of 1.6×10−2 and 9×10−5 at 1 m and 10 m distance from the well.

Maccarini, Alessandro↗
DOE OSTI2022

Using Sentinel-1 and GRACE satellite data to monitor the hydrological variations within the Tulare Basin, California

Abstract Subsidence induced by groundwater depletion is a grave problem in many regions around the world, leading to a permanent loss of groundwater storage within an aquifer and even producing structural damage at the Earth’s surface. California’s Tulare Basin is no exception, experiencing about a meter of subsidence between 2015 and 2020. However, understanding the relationship between changes in groundwater volumes and ground deformation has proven difficult. We employ surface displacement measurements from Interferometric Synthetic Aperture Radar (InSAR) and gravimetric estimates of terrestrial water storage from the Gravity Recovery and Climate Experiment (GRACE) satellite pair to characterize the hydrological dynamics within the Tulare basin. The removal of the long-term aquifer compaction from the InSAR time series reveals coherent short-term variations that correlate with hydrological features. For example, in the winter of 2018–2019 uplift is observed at the confluence of several rivers and streams that drain into the southeastern edge of the basin. These observations, combined with estimates of mass changes obtained from the orbiting GRACE satellites, form the basis for imaging the monthly spatial variations in water volumes. This approach facilitates the quick and effective synthesis of InSAR and gravimetric datasets and will aid efforts to improve our understanding and management of groundwater resources around the world.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Density streams in the disc winds of Classical T Tauri stars

ABSTRACT Spectral and photometric variability of the Classical T Tauri stars RY Tau and SU Aur from 2013 to 2022 is analysed. We find that in SU Aur the H α line’s flux at radial velocity RV = −50 ± 7 km s−1 varies with a period P = 255 ± 5 d. A similar effect previously discovered in RY Tau is confirmed with these new data: P = 21.6 d at RV = −95 ± 5 km s. In both stars, the radial velocity of these variations, the period, and the mass of the star turn out to be related by Kepler’s law, suggesting structural features on the disc plane orbiting at radii of 0.2 au in RY Tau and 0.9 au in SU Aur, respectively. Both stars have a large inclination of the accretion disc to the line of sight – so that the line of sight passes through the region of the disc wind. We propose there is an azimuthal asymmetry in the disc wind, presumably in the form of ‘density streams,’ caused by substructures of the accretion disc surface. These streams cannot dissipate until they go beyond the Alfven surface in the disc’s magnetic field. These findings open up the possibility to learn about the structure of the inner accretion disc of CTTS on scales less than 1 au and to reveal the orbital distances related to the planet’s formation.

Petrov, P. P. (ORCID:0000000155514963)↗
DOE OSTI2020

Surrogate Suspension-Resuspension Study

The objective of this study was to conduct a series of tests looking at the deposition and resuspension of aerosol particles deposited onto multiple representative substrate surfaces for a range of particle sizes under varying environmental conditions. The benefit of this study is to provide additional insight into the understanding of early time resuspension from different mechanisms and compare to existing literature. The resuspension methods utilized in this study were full-scale and the substrates were representative of real- world ground level conditions. Multiple experiments were conducted to assess the impact on resuspension from the varying substrates and mechanisms. The results of this study show variations in the size distribution of aerosol as a function of height from the source resuspension factors. Additionally, the aerosolized mass concentration and resuspension factor were evaluated. The maximum resuspension factor was found to be on the order of 10 -4 m -1 which is higher than most resuspension factors found in literature but represents idealized conditions due to the well constrained experimental setup.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Composition, Activity, and Stability of IrO x Oxygen Evolution Reaction Electrocatalysts

The oxygen evolution reaction (OER) is integral to several electrochemical energy conversion and storage technologies, including carbon dioxide reduction to value added fuels, nitrogen reduction to ammonia, reversible fuel cells, rechargeable metal−air batteries, and water electrolysis to produce hydrogen. Iridium oxide (IrO x ) is widely recognized as the benchmark OER catalyst for acidic environments. Despite widespread use of IrO x catalysts, most notably in proton-exchange membrane water electrolyzers (PEMWEs), a comprehensive understanding of the physicochemical properties of commercial catalysts and the impact of these properties on both the activity and stability of these catalysts is lacking. Here, we study commercial IrO x catalysts with different physicochemical properties, three nominally considered amorphous and three rutile, to elucidate how structural and compositional variations affect OER activity and stability. Utilizing standardized aqueous electrochemical protocols, time-resolved dissolution quantification using inductively-coupled plasma mass spectrometry, and physicochemical characterization, including multiple synchrotron X-ray techniques, we systematically correlate catalyst properties with OER performance and degradation behavior aided by principal component analysis (PCA). Our results demonstrate the general trend of amorphous IrO x having higher intrinsic activity but limited stability and crystalline rutile IrO 2 having lower activity but enhanced stability against dissolution. The trends within the amorphous and rutile catalyst groups correlate with inherent material properties, including phase composition and structure, crystallinity, particle size, surface area, and surface structure/chemistry. Notably, we identify a rutile catalyst with the largest crystallite/ domain sizes, moderate surface area, a small fraction of hydrous phase, and a favorable pore structure (trimodal distributions of pore sizes ranging from 2−5 nm) that exhibits the best balance between activity and stability among the six catalysts studied here. These findings illustrate a fundamental structure-governed trade-off between activity and stability and highlight the critical role of surface chemistry modification and structure engineering in IrO x catalyst optimization.

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

Satellite-based evaluation of AeroCom model bias in biomass burning regions

Global models are widely used to simulate biomass burning aerosol (BBA). Exhaustive evaluations on model representation of aerosol distributions and properties are fundamental to assess health and climate impacts of BBA. Here we conducted a comprehensive comparison of Aerosol Comparisons between Observations and Models (AeroCom) project model simulations with satellite observations. A total of 59 runs by 18 models from three AeroCom Phase-III experiments (i.e., biomass burning emissions, CTRL16, and CTRL19) and 14 satellite products of aerosols were used in the study. Aerosol optical depth (AOD) at 550 nm was investigated during the fire season over three key fire regions reflecting different fire dynamics (i.e., deforestation-dominated Amazon, Southern Hemisphere Africa where savannas are the key source of emissions, and boreal forest burning in boreal North America). The 14 satellite products were first evaluated against AErosol RObotic NETwork (AERONET) observations, with large uncertainties found. But these uncertainties had small impacts on the model evaluation that was dominated by modeling bias. Through a comparison with Polarization and Directionality of the Earth’s Reflectances measurements with the Generalized Retrieval of Aerosol and Surface Properties algorithm (POLDER-GRASP), we found that the modeled AOD values were biased by –93 % to 152 %, with most models showing significant underestimations even for the state-of-the-art aerosol modeling techniques (i.e., CTRL19). By scaling up BBA emissions, the negative biases in modeled AOD were significantly mitigated, although it yielded only negligible improvements in the correlation between models and observations, and the spatial and temporal variations in AOD biases did not change much. For models in CTRL16 and CTRL19, the large diversity in modeled AOD was in almost equal measures caused by diversity in emissions, lifetime, and the mass extinction coefficient (MEC). We found that in the AeroCom ensemble, BBA lifetime correlated significantly with particle deposition (as expected) and in turn correlated strongly with precipitation. Additional analysis based on Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) aerosol profiles suggested that the altitude of the aerosol layer in the current models was generally too low, which also contributed to the bias in modeled lifetime. Modeled MECs exhibited significant correlations with the Ångström exponent (AE, an indicator of particle size). Comparisons with the POLDER-GRASP-observed AE suggested that the models tended to overestimate the AE (underestimated particle size), indicating a possible underestimation of MECs in models. The hygroscopic growth in most models generally agreed with observations and might not explain the overall underestimation of modeled AOD. Our results imply that current global models contain biases in important aerosol processes for BBA (e.g., emissions, removal, and optical properties) that remain to be addressed in future research.

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

The impact of environmental effects on active galactic nuclei: A decline in the incidence of ionized outflows

Active galactic nuclei (AGNs) have generally been considered to be less frequent in denser environments due to a lower number of galaxy-galaxy interactions and/or the removal of their gas-rich reservoirs by the dense intergalactic medium. However, recent observational and theoretical works suggest that the effect of ram-pressure stripping acting on galaxies in dense environments might reduce the angular momentum of their gas, causing it to infall towards the super massive black hole at their centre, activating the AGN phase. In this work we explore the connection between environment and nuclear activity by evaluating the variation in the incidence of ionized outflows, a common phenomenon associated with nuclear activity, in AGNs across different environments. We select a sample of approximately 3300 optical AGNs from the Sloan Digital Sky Survey Data Release 13, which we match with a group catalogue to identify galaxies in isolation or residing in groups. We further probe their environments through the projected distance to the central galaxy of the group or cluster and the projected surface density to the fifth neighbour (δ 5 ). The presence of ionized outflows is determined through the modelling of the [O III] λ5007 emission line. We find that at lower masses (< 10 10.3 M ⊙ ), the fraction of ionized outflows is significantly lower in satellite AGNs (~7%) than in isolated (~22%) AGNs, probably due to their different AGN luminosity, L [O III] , in this stellar mass range. The fraction of outflows decreases at distances closer to the central galaxy of the group or cluster for all satellite AGNs; however, only the lower-mass ones (10 9 - 10 10.3 M ⊙ ) display a significant decline with δ 5 . Although this study does not include AGNs in the densest regions of galaxy clusters, our findings suggest that AGNs in dense environments accrete less gas than those in the field, potentially due to the removal of the gas reservoirs via stripping or starvation, leading to a negative connection between environment and AGN activity. Based on our results, we propose that the observed change in the incidence of outflows when moving towards denser regions of groups and clusters could contribute to the higher gas metallicities of cluster galaxies compared to field galaxies, especially at lower masses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Chemical Mixing Induced by Internal Gravity Waves in Intermediate-mass Stars

Internal gravity waves can cause mixing in the radiative interiors of stars. We study this mixing by introducing tracer particles into 2D hydrodynamic simulations. Following the work of Rogers & McElwaine, we extend our study to different masses (3, 7, and 20 M ⊙ ) and ages (ZAMS, midMS, and TAMS). The diffusion profiles of these models are influenced by various parameters such as the Brunt–Väisälä frequency, density, thermal damping, the geometric effect, and the frequencies of waves contributing to these mixing profiles. We find that the mixing profile changes dramatically across age. In younger stars, we noted that the diffusion coefficient increases toward the surface, whereas in older stars the initial increase in the diffusion profile is followed by a decreasing trend. We also find that mixing is stronger in more massive stars. Hence, future stellar evolution models should include this variation. In order to aid the inclusion of this mixing in 1D stellar evolution models, we determine the dominant waves contributing to these mixing profiles and present a prescription that can be included in 1D models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Pore-Scale Transport Effects in Electrochemical CO 2 Reduction on Gold via Coupled Microkinetic-Transport Modeling

A pore-resolved modeling framework is developed to quantify how pore-scale transport affects the intrinsic microkinetics of CO 2 -to-CO on Au. A DFT-informed microkinetic model is coupled self-consistently to a Generalized-Modified Poisson–Nernst–Planck (GMPNP) transport description in a single, electrolyte-filled cylindrical pore, allowing local concentrations and potential to feed back into site-specific reaction rates. FIB-SEM is used to determine pore sizes within realistic electrode materials. Across pore diameters, d p = 10–6000 nm, the surface-averaged CO 2 reduction rate is systematically reduced relative to the ideal microkinetic baseline where mass transport is not accounted for; the effectiveness factor 𝜂 𝑠,CO 2 , which quantifies this ratio, decreases rapidly at more negative potentials and is about 1% near −1.0 V vs SHE due to reactant depletion. Spatial maps reveal pore-bulk alkalization that emerges at higher cathodic bias, with a small, near-wall pH dip due to electrostatic repulsion of hydroxide at the cathode interface. For a fixed aspect ratio L p /d p , narrower pores exhibit larger 𝜂 𝑠,CO 2 by shortening diffusion paths, whereas variations in the aspect ratio L p /d p play a secondary role. A dimensionless analysis (surface/bulk Damköhler numbers) delineates operating regimes. In conclusion, this work offers a concept for incorporating microkinetic models into homogenized porous-electrode models through effectiveness factors and pore-size distribution.

Au-catalyst↗
DOE OSTI2022

Visualizing the interplay of Dirac mass gap and magnetism at nanoscale in intrinsic magnetic topological insulators

In intrinsic magnetic topological insulators, Dirac surface-state gaps are prerequisites for quantum anomalous Hall and axion insulating states. Unambiguous experimental identification of these gaps has proved to be a challenge, however. Here, we use molecular beam epitaxy to grow intrinsic MnBi 2 Te 4 thin films. Using scanning tunneling microscopy/spectroscopy, we directly visualize the Dirac mass gap and its disappearance below and above the magnetic order temperature. We further reveal the interplay of Dirac mass gaps and local magnetic defects. We find that, in high defect regions, the Dirac mass gap collapses. Ab initio and coupled Dirac cone model calculations provide insight into the microscopic origin of the correlation between defect density and spatial gap variations. Finally, this work provides unambiguous identification of the Dirac mass gap in MnBi 2 Te 4 and, by revealing the microscopic origin of its gap variation, establishes a material design principle for realizing exotic states in intrinsic magnetic topological insulators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Turbulence and Cloudiness in Cumulus-Topped Marine Boundary Layers

Turbulence and cloudiness in cumulus-topped marine boundary layers are studied using data from the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) site. The analysis includes eight periods of nonprecipitating shallow cumulus clouds, spanning 141 h and encompassing 603 individual clouds. On average, the cumulus had bases at 558 m, were 99 m thick, had a chord length of 500 m, and exhibited an hourly base-layer cloudiness of 12%. Changes in cloud fraction were primarily driven by variations in cloud number not chord length. High-resolution Doppler lidar and cloud radar observations were used to estimate updraft and downdraft mass fluxes in the cumulus base layer under both clear and cloudy regions. Cloudy updrafts contributed only ∼23% of the total updraft mass flux at cloud base, indicating that most upward mass flux originated from clear-air eddies. Updraft strength, rather than updraft fraction, was found to predominantly control both clear and cloudy updraft mass fluxes. Cloud-base cloudiness strongly correlated with cloudy updraft mass flux but showed negligible correlation with clear-air updraft mass flux. Clear-air downdraft mass flux exhibited a strong relationship with the ratio of surface buoyancy to inversion strength. Mesoscale analysis revealed that moist patches had greater cloudiness, updraft mass flux, and vertical velocity variance compared to dry patches. Additionally, mesoscale base-layer cloud fraction was highly and significantly correlated with cloudy updraft mass flux in both dry and moist environments. Results presented herein have implications for cumulus parameterization development along with climatological and model evaluation studies conducted at the ENA site.

cumulus clouds↗
DOE OSTI2025

Implications of reduced-complexity aerosol thermodynamics on organic aerosol mass concentration and composition over North America

Atmospheric organic aerosol (OA) mass concentrations can be affected by water uptake through its impact on the gas–particle partitioning of semi-volatile compounds. Current chemical transport models (CTMs) neglect this process. We have implemented the Binary Activity Thermodynamics model coupled to a volatility basis set partitioning scheme in the GEOS-Chem CTM, providing an efficient reduced-complexity OA model that predicts relative-humidity-dependent mixing and partitioning thermodynamics, while limiting the impact on computational efficiency. We provide a quantitative assessment of this water-sensitive OA treatment, focusing on a subdomain over North America. The updated OA scheme predicts a spatiotemporal mean enhancement in surface-level OA mass concentration of 145 % for January 2019 and 76 % for July 2019 compared to GEOS-Chem's most advanced OA scheme. The temporal mean surface-level OA organic mass concentration can increase by up to ∼590 % for January 2019 and ∼280 % for July 2019, with the greatest enhancements occurring over the ocean. The updated OA scheme also quantifies the OA-associated water content. The simulations show how different OA precursors and related OA surrogates contribute and respond to water uptake, including that due to changes in temperature and relative humidity over the diurnal cycle in selected winter and summer months. These results are independent of future CTM improvements involving updates to chemical reaction schemes and emission inventories. Our water-sensitive OA scheme allows for a better representation of the seasonal and regional variations in OA mass concentration in CTMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Electrochemistry of PUREX: R is for Reduction and Ion Transfer

The PUREX (Plutonium, Uranium, Reduction, EXtraction) process is the foundation of all industrial activities that involve the recycling of plutonium and uranium from used nuclear fuel. With over 70 years of research, engineering, and operations, it can be argued that little opportunity is left for further discovery and new fundamental understandings regarding PUREX, in general, and tri-n-butyl phosphate (TBP), in particular. Through use of an advanced electroanalytical approach, molecular-level insights regarding the back-extraction of Pu(IV) by reductive separation is reported. Studies of model PUREX systems, with lanthanide ions extracted into 20 % TBP in n-dodecane, are used to demonstrate the approach. These studies reveal knowledge pertinent to mass transfer aspects of the biphasic reductive stripping reaction central to the separation of U and Pu in PUREX. Two lanthanide ions—Ce(IV) and Yb(III)—were extracted from their concentrated solutions in aqueous 3 M HNO 3 to expressly provide electrically conducting third phases to facilitate electrochemical data acquisition. To simplify the diabolically-complicated reductive separation aspects of PUREX, the controlled polarization of electrode surfaces was used to reduce Ce(IV) and Yb(III) in place of chemical reductants. Electroanalyses of the biphasic systems with Ce(IV) and Yb(III) using three-phase electrode differential pulse voltammetry demonstrate the effects of mass transfer across water-oil interfaces on the electrode potentials for the reduction to Ce(III) and Yb(II). These reductions trigger nitrate transfer out of the oil phase for Ce and proton transfer into the organic phase for Yb. This is in correlation with the transfer of the charge-neutral ion pairs Ce·3NO 3 and Yb·2NO 3 out of the oil phases. Here, the differences in ion transfer steps reflect significant variations in the equilibrium speciation of multinuclear Ce(IV) and mononuclear Yb(III) solvates of TBP in the oil phases and the loading of the organic phase.

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

Dynamic signatures of electronically nonadiabatic coupling in sodium hydride: a rigorous test for the symmetric quasi-classical model applied to realistic, ab initio electronic states in the adiabatic representation

Sodium hydride (NaH) in the gas phase presents a seemingly simple electronic structure making it a potentially tractable system for the detailed investigation of nonadiabatic molecular dynamics from both computational and experimental standpoints. The single vibrational degree of freedom, as well as the strong nonadiabatic coupling that arises from the excited electronic states taking on considerable ionic character, provides a realistic chemical system to test the accuracy of quasi-classical methods to model population dynamics where the results are directly comparable against quantum mechanical benchmarks. Here, using a simulated pump–probe type experiment, this work presents computational predictions of population transfer through the avoided crossings of NaH via symmetric quasi-classical Meyer–Miller (SQC/MM), Ehrenfest, and exact quantum dynamics on realistic, ab initio potential energy surfaces. The main driving force for population transfer arises from the ground vibrational level of the D 1 Σ + adiabatic state that is embedded in the manifold of near-dissociation C 1 Σ + vibrational states. When coupled through a sharply localized first-order derivative coupling most of the population transfers between t = 15 and t = 30 fs depending on the initially excited vibronic wavepacket. While quantum mechanical effects are expected due to the reduced mass of NaH, predictions of the population dynamics from both the SQC/MM and Ehrenfest models perform remarkably well against the quantum dynamics benchmark. Additionally, an analysis of the vibronic structure in the nonadiabatically coupled regime is presented using a variational eigensolver methodology.

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

Electrochemical Corrosion and Catalysis Dynamics of Tin Oxide during Water Oxidation

Metal oxide corrosion severely limits anodic electrocatalysis, particularly at high potentials in acidic environments, where degradation pathways remain poorly defined. This study establishes explicit connections between corrosion and electrocatalysis on tin oxide during water oxidation by examining the roles of lattice defects, reactive oxygen species, interfacial pH variations, and speciation of corroded tin in acid. We first demonstrate the presence of structural defects such as oxygen vacancies and substoichiometric Sn(II) species by integrating electron paramagnetic resonance spectroscopy, ultraviolet photoelectron spectroscopy, and Mott–Schottky analysis. Kohn–Sham density functional theory calculations reveal that explicit water structures thermodynamically stabilize reaction intermediates and lower reaction overpotentials. Moreover, we propose that water dissociation leads to hydrogen-bonding networks formed by H* and OH* intermediates, which may span the entire catalyst surface and decrease the interfacial pH to drive corrosion. In contrast, the electrochemical generation of reactive oxygen species is shown to play a minor role in catalyst corrosion during water oxidation using inductively coupled plasma mass spectrometry coupled with selective chemical scavengers. Square-wave voltammetry combined with rotating ring-disk electrodes is used to reveal that under open-circuit conditions, only Sn(IV) cations chemically dissolve from tin oxide, while both Sn(IV) and Sn(II) species electrochemically corrode during water oxidation. Our results unveil a dynamic and complicated interplay between corrosive and catalytic pathways on metal oxide electrocatalysts: a decrease in interfacial pH due to water oxidation exacerbates Sn(II)/Sn(IV) corrosion. Subsequently, the electrochemical corrosion of Sn(II)/Sn(IV) facilitates product formation from lattice oxygen, while the redeposition of corroded Sn(II) as Sn(IV) can enable oxygen exchange with water. By elucidating the roles of defects and interfacial chemistry, this work provides a roadmap for engineering improved electrocatalysts that balance activity and stability, a critical step toward scalable and durable energy technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effect of soot concentration on radiative transport in a 2 m JP-8 pool fire

Hubbard et al. (2024) investigated conjugate heat transfer in a 7.9 m JP-8 pool fire, utilizing both heavily sooting and non-sooting fuel models. Their findings indicated that all models produced pool radiative heat fluxes within the measurement uncertainty for large-scale outdoor experiments. This study aims to further investigate the influence of soot concentration on radiative transport. To this end, Large Eddy Simulations (LES) of a 2 m diameter JP-8 pool fire were conducted, employing four turbulent combustion chemistry models: the Eddy Dissipation Concept (EDC) and three variants of flamelet model. Additionally, three variations of the one-equation subgrid scale turbulence kinetic energy model, and multiple mesh resolutions were utilized. The primary validation metrics were pool radiative heat flux and plume radiative fraction, with simulation data compared to experimental results published by Blanchat et al. (2010). The EDC model produced high soot concentrations near the pool surface, resulting in a significant shielding effect. In contrast, the flamelet combustion model has different kinetic limiters for soot production, leading to soot formation higher in the flame, where it can be more readily oxidized. Flamelet models incorporating both acetylene and benzene contributions to soot growth demonstrated better agreement with experimental data. To assess the sensitivity of pool radiative heat flux to soot mass concentration, we used a one-dimensional approximation over a wide range of soot concentrations and compared it to three-dimensional solutions obtained using Discrete Ordinates. Furthermore, this analysis illustrates the impacts of soot shielding and Turbulence Radiation Interactions (TRI), with TRI emission enhancement factors predicted to be in the range of 3–5 near the edges of the plume.

Computational fluid dynamics↗
DOE OSTI2022

Investigation of Coupled Processes in Fractures and the Bordering Matrix via a Micro‐Continuum Reactive Transport Model

Abstract In multi‐mineral fractured rocks, the altered porous layer on the fracture surface resulting from preferential dissolution of the fast‐reacting minerals can have profound impacts on subsequent chemical‐physical alteration of the fractures. This study adopts the micro‐continuum approach to provide further understanding of reactive transport processes in the altered layer (AL), and mass exchanges with the bordering matrix and fracture. The modeling framework couples the Darcy‐Brinkman‐Stokes (DBS) solver in COMSOL Multiphysics and the geochemical modeling capability of CrunchFlow. Three‐dimensional steady state simulations with systematically varied chemical‐physical parameters of the AL were performed to examine the impacts of individual factors and processes. Our simulation results confirm previous observations that dissolution of the fast‐reacting mineral (i.e., calcite) is largely controlled by diffusion across the AL. We also show that dissolution of the slow‐reacting mineral (i.e., dolomite), which controls AL development and fracture enlargement, increases with surface area and has a complex dependence on different local rate‐limiting processes. In particular, advection can result in evident spatial variations in the local dissolution rates of dolomite, although it does not affect the bulk chemistry significantly. The difference in the spatial patterns between simulations with and without advection in the AL is more noticeable in the locations with smaller apertures, with up to 20% difference in local reaction rates. Therefore, it is important to include a full depiction of advection, diffusion, and reactions for accurately capturing local dynamics that control long‐term fracture evolution.

58 GEOSCIENCES↗