Search NASA⌕ Search

SEARCH · Search NASA

Results for “scaling reactions”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Effects of gravity on sheared and nonsheared turbulent nonpremixed flames

The present numerical study is concerned with the fundamental physics of the multiway interaction between turbulence, chemical reaction, and buoyancy in a nonpremixed flame. The method of direct numerical simulation (DNS) is used to solve the instantaneous, three-dimensional governing equations. Because of the present supercomputer limitations, we consider two simple flow geometries, namely an initially uniform flow without shear (equivalent to grid-generated turbulence) and an initially uniform shear flow. In each flow, the fuel and oxidant initially exist as two separate streams. As the reactants mix, chemical reaction takes place and exothermic energy is released causing variations in density. In the presence of a gravity field, the spatial and temporal distributions of the induced buoyancy forces depend on the local density gradients and the direction of the gravitational acceleration. The effects of buoyancy include the generation of local shear, baroclinic production or destruction of vorticity, and countergradient heat and mass transport. Increased vorticity and small-scale turbulence promote further mixing and reaction. However, if the strain-rates become too high, local flame extinction can occur. Our objective is to gain an understanding of the complex interactions between the physical phenomena involved, with particular attention to the effects of buoyancy on the turbulence structure, flame behavior, and factors influencing flame extinction.

Elghobashi, Said↗
DOE OSTI2021

Influence of microstructure on replacement and porosity generation during experimental dolomitization of limestones

Replacement reactions commonly alter the multiscale pore structures of rocks during fluid-rock interactions. Analysis of these processes in various model fluid-rock systems during controlled laboratory experiments provides insights into the origins of microstructures found in natural materials. This study focused on understanding the effects of initial starting material permeability and resultant differences in transport pathways on porosity and mineralogical changes during limestone dolomitization.A series of replacement experiments (32–317 days in duration) have been conducted in which 1.59 cm (5/8 in.) diameter cores of two different limestones were reacted with saturated MgCl2 solutions at 200 °C. The Texas Cream (Austin Chalk) is a high-porosity, high-permeability limestone, whereas both the porosity and permeability of the Carthage Marble (Burlington Limestone) are relatively low. Altered limestones were imaged using scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM-EDX), Time-of-Flight Secondary Ion Mass Spectrometry (ToF–SIMS) and electron microprobe analysis (EMPA). A representative grain boundary of the low-porosity limestone was targeted for a focused ion beam (FIB) lift-out and characterized using transmission electron microscopy (TEM). These results were coupled with analyses of radial changes in the porosity distribution of the core derived from X-ray and neutron small- and ultra-small angle scattering ((U)SANS/(U)SAXS).The high-porosity/permeability limestone showed a four times faster bulk replacement rate than the lower-porosity/permeability material, and a different mechanism of porosity development. For the low-porosity limestone, a two-stage replacement occurred, with the reacted region of the core consisting of an inner rim in which the limestone was replaced by two calcite-dolomite solid solutions, and an outer rim in which the dolomite was replaced by magnesite. Elongated pores formed along grain boundaries at the initial limestone/dolomite reaction interface, and additional nanometer-scale porosity was formed at the secondary magnesite replacement rim. Grain boundaries were identified as preferential pathways for transport leading to dolomitization and a grain boundary diffusion rate was calculated based on microstructural characterization. In contrast, replacement in the high-porosity limestone was accompanied by porosity generation through replacement of individual grains by dolomite throughout the sample and, in longer runs, magnesite in outer parts of the core. These observations emphasize that both the mechanisms of the replacement reaction and the microstructure and chemistry of the replaced product are contingent on the initial structure of the starting material.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Kinetic model development for single step ethanol to butene rich olefin process over Cu-Y/Beta catalysts

Here, this study presents the first intrinsic kinetic model for the single-step conversion of ethanol to butene-rich olefins over bifunctional Cu-Y/Beta catalysts, addressing a critical gap in the design and scale-up of Sustainable Aviation Fuel (SAF) processes. The reaction network comprises ten global steps involving dehydrogenation, aldol condensation, hydrogenation, and dehydration reactions, distributed across Cu and rare-earth (Y) active sites. The model incorporates dual-site functionality (Cu and Y site) and explicitly accounts for key intermediates such as crotonaldehyde and butanal. Reaction rates are formulated using Langmuir–Hinshelwood–Hougen–Watson (LHHW) kinetics. Kinetic parameters are extracted by fitting the model to lab-scale packed-bed reactor data across a wide range of temperatures and space velocities, demonstrating strong agreement in ethanol conversion and product selectivity. The reaction kinetics developed in this work provide a foundational basis for constructing reactor models that enable process optimization and scale-up of ethanol-to-jet fuel technologies.

Cu-Y/Beta catalyst↗
NASA NTRS1973

Effects of reaction control system jet simulation on the stability and control characteristics of a 0.015-scale space shuttle orbiter model in the Ames Research Center 3.5-foot hypersonic wind tunnel

An experimental investigation was performed in the Ames Research Center 3.5-Foot Hypersonic Wind Tunnel to obtain detailed effects which interactions between the RCS jet flow field and the local orbiter flow field have on orbiter hypersonic stability and control characteristics. Six-component force data were obtained through an angle-of-attack range of 15 to 35 deg with 0 deg angle of sideslip. The test was conducted with yaw, pitch and roll jet simulation at a free-stream Mach number of 10.3. These data simulate two SSV reentry flight conditions at Mach numbers of 28.3 and 10.3. Fuselage base pressures and pressures on the nonmetric RCS pods were obtained in addition to the basic force measurements. Model 42-0 was used for these tests.

Dziubala, T. J.↗
NASA NTRS1974

Effects of reaction control system jet simulation on the stability and control characteristics of a 0.015 scale space shuttle orbiter model tested in the Langley Research Center unitary plan wind tunnel

An experimental investigation was performed in the Langley Research Center Unitary Plan Wind Tunnel (Test 0A70) to obtain the detailed effects that RCS jet flow interactions with local orbiter flow field have on supersonic stability and control characteristics of the space shuttle orbiter. Six-component force data were obtained through an angle-of-attack range from 15 to 35 degrees at angles of sideslip of 0, +5, and -5 degrees. The test was conducted with yaw jet simulation at free-stream Mach numbers of 2.5 and 4.6, simulating SSV re-entry flight conditions at these Mach numbers. In addition to the basic force measurements, fuselage base pressures and pressures on the non-metric RCS pods were obtained.

Daileda, J. J.↗
DOE OSTI2021

Quantification of heterogeneous, irreversible lithium plating in extreme fast charging of lithium-ion batteries

Realization of extreme fast charging (XFC, ≤15 minutes) of lithium-ion batteries is imperative for the widespread adoption of electric vehicles. However, dramatic capacity fading is associated with XFC, limiting its implementation. To quantitatively elucidate the effects of irreversible lithium plating and other degradation mechanisms on the cell capacity, it is important to understand the links between lithium plating and cell degradation at both the local and global (over the full cell) scales. Here, we study the nature of local lithium plating after hundreds of XFC cycles (charging C-rates ranging from 4C to 9C) in industrially-relevant pouch cells using spatially resolved X-ray diffraction. Our results reveal a spatial correlation at the mm scale between irreversible lithium plating on the anode, inactive lithiated graphite phases, and local state-of-charge of the cathode. In regions of plated lithium, additional lithium is locally and irreversibly trapped as lithiated graphite, contributing to the loss of lithium inventory (LLI) and to a local loss of active anode material. The total LLI in the cell from irreversibly plated lithium is linearly correlated to the capacity loss in the batteries after XFC cycling, with a non-zero offset originating from other parasitic side reactions. Finally, at the global (cell) scale, LLI drives the capacity fade, rather than electrode degradation. We anticipate that the understanding of lithium plating and other degradation mechanisms during XFC gained in this work will help lead to new approaches towards designing high-rate batteries in which irreversible lithium plating is minimized.

25 ENERGY STORAGE↗
DOE OSTI2023

Effects of sub‐mm cylindrical voids on detonation performance in PBX 9501

Abstract Internal features of varying scale and geometry are always present in explosives systems. Below a critical length, dependent on the explosive, these features can operate as a driving force for energy concentration and reaction, known as hotspots. At larger length scales, internal features result in jetting and detonation wave shaping, allowing for bulk work to be done by the explosive as is seen in shape charges. To date, a large volume of work has been performed to simulate hotspot ignition and large‐scale wave shaping. However, little experimental data exists on the effects of features in intermediate length scales, 0.1 to 1 millimeter. It has been observed in many tests that these small‐scale features can influence the high explosive (HE) performance and in some cases cause substantial damage to adjacent systems. This work provides quantitative data on the effects of machined voids moderately above more typical hotspot lengths, 0.3–0.8 mm in diameter, in PBX 9501 pressed to 1.785 g cm −3 ±2.5 mg cm −3 . Streak imaging was used to visualize void collapse, jet velocity, re‐initiation time, and wave shape evolution. Delay in detonation front propagation time was found to be linearly dependent on the void diameter and jet velocity was found to be independent within the tested range and resolution. Cut‐back experiments were used to investigate wave shape distortion and evolution downstream of the void, showing consistent growth and decay shapes across all void sizes. Simulations using CTH, a hydrocode by Sandia National Laboratory, were used to investigate void collapse showing agreement with the trends of experimental results but yielded inaccurate wave shape development delay values. Both simulation and experimental results identified several re‐initiation mechanisms with jetting and subsequent double shocking of localized HE being the dominant mechanisms for the void sizes that were studied.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

Long-Term Evaluation of Remedial Technology Performance in the Laboratory: Multi-year Experimental Test Plan

An understanding of the long-term effectiveness of remediation technologies is central to sustainable environmental cleanup. Long-term experiments are valuable for reducing uncertainty and predicting remediation outcomes at scales required for regulatory compliance. However, these types of tests can be costly and challenging to interpret. Therefore, contaminated sites often rely on short-term laboratory experiments that may not account for the potentially significant effects of gradual, time- dependent processes governing contaminant retention, release, and species transformation. For example, short-term lab experiments (from months to a year) conducted with sediments from the unsaturated and saturated zones at the Hanford Site play an important role in initial evaluations of remediation technologies but cannot capture the full extent of time-dependent reactions, leaving uncertainties in field-scale deployment. Here, long-term testing will be conducted on select technologies based on their performance in short-term testing. The overall objective is to directly address the challenges described above by generating and analyzing data on long-term (2–10 years) efficacy of selected remediation technologies, integrating this understanding into models, and providing critical input for remediation planning, monitoring, and 5-year review cycles for field-implemented remedies. Specific objectives include: 1. Evaluating long-term efficiency of promising remedies under site-specific conditions. 2. Generating robust parameters for modeling, reducing uncertainty in predictive simulations. 3. Advancing integrated monitoring by combining geochemical and geophysical observations. 4. Informing field-scale implementation by incrementally advancing technologies, identifying failure mechanisms early, and prioritizing robust, cost-effective remedies. Through systematic evaluation of technologies in laboratory-scale column experiments, integrated monitoring, and modeling support, this project is designed to bolster confidence in the long-term robustness of selected technologies. The ultimate outcome is the identification and deployment of more reliable, cost-effective remedies that safeguard human health and the environment while reducing the uncertainties that have historically hindered cleanup progress at Hanford Site. An experimental approach was developed and initiated for long-term testing potentially up to 10 years. The table below summarizes the experimental approach developed for testing select technologies and presented in this multi-year experimental test plan.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1971

Pyrotechnic hazards classification and evaluation program. Run-up reaction testing in pyrotechnic dust suspensions

A preliminary investigation of the parameters included in run-up dust reactions is presented. Two types of tests were conducted: (1) ignition criteria of large bulk pyrotechnic dusts, and (2) optimal run-up conditions of large bulk pyrotechnic dusts. These tests were used to evaluate the order of magnitude and gross scale requirements needed to induce run-up reactions in pyrotechnic dusts and to simulate at reduced scale an accident that occurred in a manufacturing installation. Test results showed that propagation of pyrotechnic dust clouds resulted in a fireball of relatively long duration and large size. In addition, a plane wave front was observed to travel down the length of the gallery.

Source record↗
DOE OSTI2021

Electrochemical ion insertion from the atomic to the device scale

Electrochemical ion insertion involves coupled ion–electron transfer reactions, transport of guest species and redox of the host. The hosts are typically anisotropic solids with 2D conduction planes but can also be materials with 1D or isotropic transport pathways. These insertion compounds have traditionally been studied in the context of energy storage but also find extensive applications in electrocatalysis, optoelectronics and computing. Recent developments in operando, ultrafast and high-resolution characterization methods, as well as accurate theoretical simulation methods, have led to a renaissance in the understanding of ion-insertion compounds. In this Review, we present a unified framework for understanding insertion compounds across timescales and length scales ranging from atomic to device levels. Using graphite, transition metal dichalcogenides, layered oxides, oxyhydroxides and olivines as examples, we explore commonalities in these materials in terms of point defects, interfacial reactions and phase transformations. We illustrate similarities in the operating principles of various ion-insertion devices, ranging from batteries and electrocatalysts to electrochromics and thermal transistors, with the goal of unifying research across disciplinary boundaries.

36 MATERIALS SCIENCE↗
NASA NTRS2022

Building and Breaking Carbon Composites with REACTER

Carbon-based composites have become indispensable materials in aerospace and other high-performance applications, yet obtaining a detailed, nanoscale understanding of their morphology and failure mechanisms using only experimental methods remains a difficult challenge. REACTER is a versatile computational modeling tool for atomistic molecular dynamics simulations designed to model chemical reactions at the speed and length scales of classical force fields. In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites: carbon nanotube (CNT) composites and carbon fiber reinforced polymers (CFRP). A network of CNTs was grown dynamically using the new ‘create atoms’ feature of REACTER. The CNT filler was embedded into a polyarylacetylene (PAA) matrix by simulated in situ polymerization to obtain the final composite model. To generate the second class of carbon composite, fully carbonized (graphitic) carbon fiber morphologies were created by the method of Desai et al. [1], but using the advanced reaction constraints framework of REACTER. Two fiber models were created, representing a circular carbon fiber core and a flat surface, and similarly infiltrated with resin to obtain the final CFRP structure. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third order DFT-based tight-binding simulations, via a reaction constraint on the total potential energy of the involved atoms.

Polymer↗
DOE OSTI2025

On-line Inductively Coupled Plasma Mass Spectrometry Reveals Material Degradation Dynamics of Au and Cu Catalysts during Electrochemical CO 2 Reduction

A significant challenge in commercializing electrochemical CO 2 reduction (CO 2 R) is achieving catalyst durability. Here, in this study, online inductively coupled mass spectrometry (ICP–MS) was used to investigate catalyst degradation via nanoparticle detachment and/or dissolution into metal ions under CO 2 R operating conditions in 0.1 M KHCO 3 . We developed an experimental framework with ex situ characterization to validate the online ICP–MS method for in situ evaluation of degradation from metal foils. By varying the applied potential and microenvironment (CO 2 vs N 2 -saturated electrolyte), we gained insights into the degradation of Au and Cu foils under CO 2 R and hydrogen evolution reaction (HER) conditions. While both Au and Cu foils were observed to be stable to dissolution in these regimes, degradation via nanoparticle detachment from the foil surface at the femtogram scale was observed as a function of reaction conditions, providing new insights into material degradation mechanisms. When applying potential steps at −0.1 and −1.0 V vs the reversible hydrogen electrode (RHE), Au was found to degrade via nanoparticle detachment under CO 2 R operating conditions more than under HER conditions, while Cu was found to degrade via nanoparticle detachment in similar amounts during both reactions. Au lost ∼1.8× more mass and ∼7.5× more nanoparticles than Cu under CO 2 R operating conditions. This study demonstrates the use of online ICP–MS to gain insight into the degradation of Au and Cu, the importance of studying unconventional degradation mechanisms such as nanoparticle detachment, and that online ICP–MS can be further utilized to gain fundamental understanding of catalyst durability for a variety of reaction systems.

Yan, Katherine [Stanford Univ., CA (United States)↗
DOE OSTI2023

Methanediol from cloud-processed formaldehyde is only a minor source of atmospheric formic acid

Atmospheric formic acid is severely underpredicted by models. A recent study proposed that this discrepancy can be resolved by abundant formic acid production from the reaction ( 1 ) between hydroxyl radical and methanediol derived from in-cloud formaldehyde processing and provided a chamber-experiment-derived rate constant, k 1 = 7.5 × 10 −12 cm 3 s −1 . High-level accuracy coupled cluster calculations in combination with E,J -resolved two-dimensional master equation analyses yield k 1 = (2.4 ± 0.5) × 10 −12 cm 3 s −1 for relevant atmospheric conditions ( T = 260–310 K and P = 0–1 atm). We attribute this significant discrepancy to HCOOH formation from other molecules in the chamber experiments. More importantly, we show that reversible aqueous processes result indirectly in the equilibration on a 10 min. time scale of the gas-phase reaction HCHO + H 2 O ⇌ HOCH 2 OH (2) with a HOCH 2 OH to HCHO ratio of only ca . 2%. Although HOCH 2 OH outgassing upon cloud evaporation typically increases this ratio by a factor of 1.5–5, as determined by numerical simulations, its in-cloud reprocessing is shown using a global model to strongly limit the gas-phase sink and the resulting production of formic acid. Based on the combined findings in this work, we derive a range of 1.2–8.5 Tg/y for the global HCOOH production from cloud-derived HOCH 2 OH reacting with OH. The best estimate, 3.3 Tg/y, is about 30 times less than recently reported. The theoretical equilibrium constant K eq (2) determined in this work also allows us to estimate the Henry’s law constant of methanediol (8.1 × 10 5 M atm −1 at 280 K).

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

A Possible Observation of Sigma-nn Continuum Structures and A Bound Sigma-NN State Using the (e, e 'K+) Reaction

The E12-17-003 (e, e'K+) experiment was carried out in the experimental Hall A facility of Jefferson Lab in November 2018. The experiment aimed at providing the experimental data for the unknown Lambda-n interaction by measuring the bound state or the resonance state as indicated by the HypHI experiment. The so obtained Lambda-n interaction was assumed to solve the existing charge symmetry breaking problem in the case of Lambda-N interaction. The experiment used Hall A high resolution spectrometers (both left and right spectrometers) positioned symmetrically at constant angle 13.2 deg. each. A 25 cm long aluminum cylinder filled with H-3 gas was used as a production target. The major calibration data were taken with the H target contained in an identical aluminum cylinder. The electroproduction of Lambda and Sigma^0 via p(e, e'K+)Lambda, Sigma^0 reactions was used to calibrate the absolute energy scale with the known masses of the Lambda and Sigma^0. Then by using the (e, e'K+) Sigma-reaction, two possible Lambda-nn resonance states and one NN bound state were observed. The observed states have the energy resolution about 1.6 MeV (FWHM), however, great r statistics are required to solidly confirm the observed states. The A = 3 and 4 Sigma bound states have been predicted long ago but only the A = 4 Sigma hypernucleus, that is He-4/Sigma, was found in the (K-, pi-) reaction. A careful Monte Carlo study was conducted to study the A (nuclear mass number) dependence on the missing mass resolution by using the identical experimental conditions. The simulated and experimentally obtained Lambda and Sigma^0 channels agreed within 100 keV in sigma, where sigma is the standard deviation. The intrinsic missing mass resolution of A = 3 (3/Lambda-n) resonance was predicted about sigma = 0.67 MeV with a natural width of about 0.6 MeV. However, due to the low statistics the precision does not permit sufficient constrain on the determination of the Lambda-n interaction. The reason of having such low statistics is due to the use of unoptimized Hall A system which was definitely feasible but not ideal for this experiment. In addition, the cross section was found to be much smaller than expected.

Pandey, Bishnu↗
DOE OSTI2022

Solar Receiver with Integrated Thermal Storage for a Supercritical Carbon Dioxide Power Cycle

This project endeavors to design and demonstrate an integrated CSP system that meets the Department of Energy’s 6¢/kW-hr LCOE goal [1]. This design is based on the 10 MWe Supercritical Transformational Electric Power (STEP) engine specifications, but adaptable to any sCO 2 power module up to 10 MW e . The enabling technology to be tested is a new class of metal hydride based thermal energy storage materials which can store thermal energy in the hydrogen bonds between the metal hydride and hydrogen. This thermochemical energy storage solution provides a much greater energy density than sensible and latent energy storage enabling an integrated system design that closely linked the metal hydride to the power block. This allows the entire system to be mounted up tower to provide cost reductions. A novel metal hydride, capable of operating above 760 °C, was developed during the course of this project and thoroughly tested at a lab scale. The material’s hydrogen capacity, reaction enthalpy, thermal conductivity, and cyclic stability were measured. A demonstration test article was built that can store 2kWh of thermal energy and was coupled to Brayton Energy’s CO 2 loop where the performance of the metal hydride under commercial conditions was measured. An integrated model was created that uses the inputs from all the different sub-models in order to return an overall system LCOE using SAM. This project greatly advanced metal hydride based thermochemical energy storage technology and demonstrated its use. Previously there were no metal hydrides which were capable of operating at the proposed temperature limits (760°C), there is now a high temperature metal hydride with documented performance/stability and a path for low-cost production. Additionally, the metal hydride was tested at a scale >650x larger than what had previously been done in the lab. While not at commercial scale, it was an excellent demonstration of the technology and provided knowledge that can be leveraged for other applications.

08 HYDROGEN↗
DOE OSTI2024

Shallow Rate-Redox Potential Scaling in Aqueous Molecular Oxygen Reduction Electrocatalysis Across a Family of Iron Macrocycles

Rate-overpotential scaling relationships have been employed widely to understand trends in oxygen reduction reaction (ORR) electrocatalysis by dissolved metal macrocycles in organic electrolytes. Similar scaling relationships remain unknown for surface-adsorbed ORR electrocatalysts in the acidic aqueous environments germane to proton-exchange membrane (PEM) fuel cells. Herein, we examine ORR catalysis in aqueous perchloric acid media for a structurally diverse array of iron macrocycle complexes adsorbed on Vulcan carbon black. The macrocycles encompass Fe– N 4 , Fe–N 2 N' 2 and Fe–N x C 4-x motifs bearing pyrrolic, pyridinic, and N-heterocyclic carbene (NHC) moieties in the primary ligation sphere, giving rise to a 670 mV range in Fe(III/II) redox potentials, E Fe(III/II) . Experimental Tafel data in the micropolarization regime were extrapolated to the E Fe(III/II) to furnish estimated per-site-normalized current density (j per-site ) values that span ~4.6 orders of magnitude across the family of compounds. Despite the structural diversity of this family of compounds, extrapolated j per-site values correlate with the Fe(III/II) redox potentials in a roughly log-linear fashion with a shallow scaling factor of approximately 145 mV/decade. Further, these findings highlight that negative shifts in E Fe(III/II) lead to diminishing returns in catalytic rate promotion and suggest that changes to the primary ligating environment in a macrocycle are insufficient to break fundamental rate-potential scaling relationships in aqueous ORR catalysis. Together these studies motivate the further development of higher-potential iron complexes that employ motifs beyond the equatorial ligation plane to enhance ORR catalysis.

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

Toward an Understanding of Linear Scaling Relations through Energy Decomposition Analysis

The discovery of linear scaling relations has fundamentally changed the field of heterogeneous catalysis. The scaling relations have been rationalized based on the d-band theory, specifically a separation of sp and d electron contributions to adsorption energies. Within the framework of energy decomposition analysis, a full understanding of such a separation would require one to further break down the adsorption energy into distinct energy components such as electrostatics, polarization, charge transfer, and van der Waals interactions, and to examine the sp and d contributions to each of them. As a step in this direction, we analyzed the interaction energy between CH x (x = 1–4) adsorbates and fcc(100) transition metal surfaces (M = Cu, Ag, Au, Rh, and Pt), with the surfaces represented both as slabs in plane-wave density functional theory (pw-DFT) calculations and as atomic clusters in atomic-orbital basis density functional theory (ao-DFT) calculations. Through an absolutely localized molecular orbital (ALMO) based energy decomposition analysis of the ao-DFT adsorption energy, each of the interaction energy components (electrostatics, polarization, van der Waals, and charge transfer) was found to follow its own scaling relations, with an intricate interplay among these energy components yielding the overall scaling relations for the total adsorption energies. Using the recently introduced ALMO-based polarization and charge-transfer analysis schemes, we further dissected polarization into metal surface and adsorbate contributions, and charge transfer into metal → adsorbate and adsorbate → metal contributions. The contributions from the sp and d electrons of the metal to these terms were further quantified, and the dominant role of the metal d electrons was reaffirmed. These results shed light on how CHx adsorbates interact with metal surfaces and further reveal the physical origin of the scaling relations.

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

Hydrological Controls on Riverbed Methane Emissions: A Numerical Investigation of Hydrodynamic and Ebullitive Mechanisms from Site to Basin Scales

Methane production and emission from riverbed sediments constitute a significant yet underexplored component of the river methane budget. Despite their importance, a mechanistic and quantitative understanding of these processes under dynamic flow conditions remains elusive, particularly at the basin scale. In this study, we investigate the multiscale (site-to-basin) mechanisms governing methane emissions from riverbeds, focusing on the interplay between biogeochemical processes and hydrological dynamics. We developed a reactive transport model that integrates methane production with hydrodynamic and ebullitive pathways and coupled it with a basin-scale groundwater flow model to simulate methane emissions across scales. By incorporating key factors such as vertical hydrological exchange flow (VHEF), particulate organic carbon availability, sediment properties, and temperature sensitivity, our model captures the spatiotemporal variability of methane ebullition at both the site and basin scales. Our results reveal that methane ebullition is strongly modulated by VHEF, with lower flow rates promoting methane accumulation and subsequent ebullition. Sensitivity analyses underscore nonlinear responses of methane fluxes to environmental drivers. Furthermore, these findings emphasize the critical roles of sediment composition, hydrological exchange, and environmental conditions in shaping methane dynamics. Overall, this work advances our understanding of riverine methane emissions and offers insights for guiding climate models and mitigation strategies.

Fluxes↗