Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “enhanced kinetics”

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

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

At least 1,513 records · Page 84

Reactive Transport Modeling of Hydrogen Production from Serpentinization of Olivine

Hydrogen production from serpentinization of ultramafic rocks represents a promising natural pathway for generating carbon-free energy, yet its kinetics and controlling factors remain incompletely understood. A key challenge in advancing serpentinization research lies in the heterogeneity of porosity and permeability in rocks, which leads to nonuniform fluid velocity fields, as well as uncertainties in estimating reactive surface area and identifying appropriate mineral reaction equilibria. Additional complexities arise from the role of dissolved SiO 2 , Fe 2+ /Fe 3+ partitioning, and the limited effect of pH variations within the strongly alkaline regime on hydrogen yields. These challenges hinder straightforward extrapolation from laboratory tests to practical applications of hydrogen production from natural rocks. Here, in this work, we address these questions using a simulation-based reactive transport modeling framework calibrated against controlled laboratory experiments reported elsewhere. The model couples geochemical kinetics, multiphase flow, and mineralogical feedbacks, enabling systematic evaluation of how surface area, dissolved silica concentration, Fe redox state, temperature, and pressure govern serpentinization and H2 generation. We find that surface area exerts the strongest control on reaction rates and hydrogen yields, while Fe 2+ /Fe 3+ ratios act as secondary modulators. Elevated dissolved silica concentrations suppress hydrogen production but accelerate serpentine precipitation, whereas increasing pH beyond 12 within the strongly alkaline regime produces only marginal gains. Finally, we demonstrate that integrating targeted experiments with calibrated simulations offers a powerful and efficient approach for predicting hydrogen yields and assessing parameter trade-offs in industrial-scale applications. This integration can substantially reduce the experimental burden while improving predictive capability, thereby enhancing both the mechanistic understanding and the practical feasibility of hydrogen production from serpentinization.

08 HYDROGEN↗

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗

Finite-element-based simulations of electrodes for CO 2 cascade reduction reactions

The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the enhanced sulfur and coking tolerance of Ni-Co-rare earth oxide catalysts for the dry reforming of methane

In this study, sulfur and coking tolerance of Ni-based dry reforming catalysts were examined. Catalysts utilizing both Ce/Zr and Ce/La oxide supports, some with additional Co, were tested. Long-term reaction runs were conducted with and without sulfur in the feed. Catalysts were also characterized by STEM, XPS, XAFS and XANES and CO chemisorption. Only catalysts where Co was also present, and supported on the Ce-Zr oxide, were capable of extended sulfur tolerance at >20 ppm sulfur. This tolerance, along with a greatly reduced coking rate, is linked to Co in intimate contact with Ni, the mixture existing as clusters anchored and influenced electronically by the oxide support. The activation of methane takes place on these sites. Larger metal aggregates formed by ripening during reaction appear to be spectators. The measured activation energies for dry reforming suggest that CO 2 activation takes place at the oxide interface, and is a kinetically significant step.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational study of inlet injection for a Pre-Mixed, Shock-Induced Combustion (PM/SIC) engine

A computational simulation of reacting 2-D and 3-D flowfields in a model inlet section of a Pre-Mixed, Shock-Induced Combustion (PM/SIC) engine concept was performed. LARCK, a multi-dimensional Navier-Stokes code with finite-rate kinetics chemistry developed at NASA LaRC by J.A. White, was adapted for this simulation. The flow conditions in the simulation match those envisioned for the PM/SIC engine experiments currently planned at LaRC. The reacting flowfields were Mach 6.3 freestream air and Mach 2 hydrogen at various pressure and temperature conditions injected through a slot injector at the base of the inlet section. In the PM/SIC engine, fuel is injected at the inlet section upstream of the combustor, and reaction is initiated by the shock wave at the inlet which increases the gas temperature and pressure beyond the kinetic limits for reaction. Many challenges exist prior to establishing shock-controlled combustion as a practical engine concept. These challenges include fuel injection schemes that can provide proper fuel-air mixing without creating large losses in the inlet section, and control of the combustion process so that early ignition or combustion propagation through the inlet boundary layer does not occur. For this project, a parametrics study was carried out to model the fuel injection of hydrogen at different flow conditions. It was found that, as the fuel temperature and pressure were increased, the potential for pre-ignition was high at a short distance downstream of the slot injector. The next stage of this work will investigate injection techniques for enhancing mixing of fuel and air in a manner that prevents or reduces the potential for premature ignition observed numerically.

Gonzalez, D. E.↗

Measurement of the double-differential cross section of muon-neutrino charged-current interactions with low hadronic energy in the NOvA Near Detector

The NOvA collaboration reports cross-section measurements for $\nu_{\mu}$ charged-current interactions with low hadronic energy (maximum kinetic energy of 250 MeV for protons and 175 MeV for pions) in the NOvA Near Detector. The results are presented as a double-differential cross section as a function of the direct observables of the final-state muon kinematics. Results are also presented as a single-differential cross section as a function of the derived square of the four-momentum transfer, $Q^{2}$, and as a function of the derived neutrino energy. The data correspond to an accumulated 8.09$\times10^{20}$ protons-on-target (POT) in the neutrino mode of the NuMI beam, with a narrow band of neutrino energies peaked at 1.8 GeV. The analysis provides a sample of neutrino-nucleus interactions with an enhanced fraction of quasi-elastic and two-particle-two-hole (2p2h) interactions. This enhancement allows quantitative comparisons with various nuclear models. We find strong disagreement between data and theory-based models in various regions of the muon kinematic phase space, especially in the forward muon direction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Stabilizing Ru-Formyl CO 2 Reduction Intermediates by Formation of Lewis Acid–Base Adducts: A Combined 2DIR and NMR Study

This work demonstrates how the stability of a Ru(II) formyl complex (cis-[Ru(bpy) 2 (CO)(CHO)][PF 6 ]), which is known to decarbonylate to form the metal hydride complex, is greatly enhanced by the addition of a Lewis acidic cation (Li + or Mg 2+ ). Multinuclear NMR is used to measure equilibrium constants of adduct formation, and then ultrafast 2DIR spectroscopy is applied to determine how the electronic structure and local molecular dynamics change upon complexation of the formyl ligand with the Lewis acid. We attribute the stabilization of the complex to the formation of a carbene-like structure where the bond strength between the Ru and the formyl group increases. The change in electronic structure is evidenced by steady-state NMR spectroscopy and DFT calculations. 2DIR measurements further confirm both the complexation of the formyl ligand with the Lewis acid and the change in electronic structure, where changes in vibrational frequency, relaxation kinetics, and a slowdown in the time scale of rotation about the metal-formyl bond are observed in the presence of Lewis acids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ideal conductor/dielectric model (ICDM): A generalized technique to correct for finite-size effects in molecular simulations of hindered ion transport

Molecular simulations serve as indispensable tools for investigating the kinetics and elucidating the mechanism of hindered ion transport across nanoporous membranes. In particular, recent advancements in advanced sampling techniques have made it possible to access translocation timescales spanning several orders of magnitude. In our prior study [Shoemaker et al., J. Chem. Theory Comput. 18, 7142 (2022)], we identified significant finite size artifacts in simulations of pressure-driven hindered ion transport through nanoporous graphitic membranes. We introduced the ideal conductor model, which effectively corrects for such artifacts by assuming the feed to be an ideal conductor. In the present work, we introduce the ideal conductor dielectric model (Icdm), a generalization of our earlier model, which accounts for the dielectric properties of both the membrane and the filtrate. Using the Icdm model substantially enhances the agreement among corrected free energy profiles obtained from systems of varying sizes, with notable improvements observed in regions proximate to the pore exit. Moreover, the model has the capability to consider secondary ion passage events, including the transport of a co-ion subsequent to the traversal of a counter-ion, a feature that is absent in our original model. We also investigate the sensitivity of the new model to various implementation details. The Icdm model offers a universally applicable framework for addressing finite size artifacts in molecular simulations of ion transport. It stands as a significant advancement in our quest to use molecular simulations to comprehensively understand and manipulate ion transport processes through nanoporous membranes.

Chemistry↗

Kinetic analysis of the template effect in ribooligoguanylate elongation

The paper presents kinetic studies on the reaction of elongation of the 3-prime-5-prime-linked ribooligoguanylates with guanosine 5-prime-phospho-2-methylimidazolide (2-MelmpG) in the presence or absence of a complementary template, the polycytidylic acid. In the absence of poly(C), the reaction leads to three isomeric oligomers that are elongated by one monomer unit: the 3-prime-5-prime linked, the 2-prime-5-prime linked, and the pyrophosphate, formed in a ratio of 1:2:5. In the presence of the template, the reaction is 20-fold faster and yields products (n + 1), (n + 2), (n + 3), etc., as long as 2-MelmpG is available. The formation of the natural, 3-prime-5-prime-linked isomer, is enhanced selectively by 140-fold at 37 C, and its relative yield increases with decreasing temperature.

Kanavarioti, Anastassia↗

Turbulence-Resolving Simulations of Hurricane Laura (2020): Insights Into Extreme Winds and Eyewall Turbulence

Extreme weather events pose significant risks to coastal and offshore energy infrastructure. In this work, we investigate the structure of mean winds and turbulence near the surface ( z < 300 m) that occurred during Hurricane Laura in 2020 on the US Gulf Coast. To this end, we perform turbulence-resolving simulations (..delta..x = 33.33 m) that encompass the entire eyewall of Hurricane Laura by introducing a novel method that we will refer to as Moving-Mesoscale to Static-LES Integrated Coupling (MOSAIC). The simulation results are validated against surface, boundary-layer, and mid-tropospheric observations. Our analysis of the extreme wind conditions near the surface indicates that the mean wind and turbulence profiles vary greatly at and near the eyewall, resulting in extreme values of mean wind speed (U > 50 m * s-1) and turbulence kinetic energy (k ~ 50 m2 * s-2) at altitudes above 50 m . To investigate in detail the nature of turbulence within the eyewall of the storm, we conduct higher-resolution simulations (..delta..x = 11.11 m) of a portion of the hurricane. We provide a comprehensive statistical description of turbulence in the eyewall region, highlighting the need to employ ..delta..x ~ 10 m grid spacing to capture the integral length-scales near the surface, second- and higher-order central moments (i.e., skewness and kurtosis), and spectral coherence in the flow accurately. These numerical simulations provide the most detailed description of mean winds and turbulent conditions within the eyewall of a historical tropical cyclone to date, illustrating how high-resolution simulations can enhance limited turbulence measurements in hurricanes and giving insight into the extreme conditions threatening offshore and coastal infrastructure.

17 WIND ENERGY↗

In-situ ETEM study of plasma-facing tungsten nanofuzz oxidation at atmospheric pressure: Microstructure evolution and substrate-free oxidation kinetics

To enable sustainable carbon-free fusion energy, managing reactor structural material degradation during normal operation as well as accident scenarios is vital. Tungsten (W) plasma-facing materials (PFMs) are susceptible to aggressive high-temperature oxidation during air-ingress fusion reactor accidents, yet there's a lack of oxidation kinetic data for irradiated tungsten. Here, in this study, we utilize atmospheric environmental transmission electron microscopy (ETEM) to present the first kinetic data for substrate-free W nanofuzz oxidation at 400 ºC and 500 ºC in 1 bar dry air. Comparison with pristine bulk W during the early parabolic stage suggests an irradiationdecelerated oxidation for W nanofuzz. Our time-resolved in-situ characterization reveals a durable amorphous surface oxide, likely promoted by high-flux He + irradiation-induced surface defects, serving as an effective passivating layer that impedes nanofuzz oxidation onset. This surface oxide layer also interfaces well with newly formed orthorhombic WO 3 , facilitated by stress relief through He bubble shrinkage, providing lasting passivating protection throughout the nanofuzz parabolic oxidation. This new finding challenges conventional notions of irradiation's negative impact on metal oxidation, and calls for advanced characterization to enhance our understanding of fusion energy materials degradation, informed by further accident modeling.

36 MATERIALS SCIENCE↗

Calculations of High-Temperature Jet Flow Using Hybrid Reynolds-Average Navier-Stokes Formulations

Two multiscale-type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds-averaged Navier Stokes equations. The first scheme is a hybrid Reynolds-averaged- Navier Stokes/large-eddy-simulation model using the two-equation k(epsilon) model with a Reynolds-averaged-Navier Stokes/large-eddy-simulation transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier Stokes model in which the unresolved kinetic energy parameter f(sub k) is allowed to vary as a function of grid spacing and the turbulence length scale. This parameter is estimated based on a novel two-stage procedure to efficiently estimate the level of scale resolution possible for a given flow on a given grid for partially averaged Navier Stokes. It has been found that the prescribed scale resolution can play a major role in obtaining accurate flow solutions. The parameter f(sub k) varies between zero and one and is equal to one in the viscous sublayer and when the Reynolds-averaged Navier Stokes turbulent viscosity becomes smaller than the large-eddy-simulation viscosity. The formulation, usage methodology, and validation examples are presented to demonstrate the enhancement of PAB3D's time-accurate turbulence modeling capabilities. The accurate simulations of flow and turbulent quantities will provide a valuable tool for accurate jet noise predictions. Solutions from these models are compared with Reynolds-averaged Navier Stokes results and experimental data for high-temperature jet flows. The current results show promise for the capability of hybrid Reynolds-averaged Navier Stokes and large eddy simulation and partially averaged Navier Stokes in simulating such flow phenomena.

Abdol-Hamid, Khaled S.↗

Enhanced oxidation resistance of (Mo 95 W 5 ) 85 Ta 10 (TiZr) 5 refractory multi-principal element alloy up to 1300°C

Refractory-metal-based alloys are a potential replacement of current nickel-based superalloys due to their excellent mechanical strength at extremely high temperatures. However, severe oxidation in a high-temperature working environment limits their application. To address this challenge, a two-step coating process (including a Mo precoat and a Si-B pack cementation) was applied to an innovative refractory multi-principal element alloy (RMPEA) (Mo 95 W 5 ) 85 Ta 10 (TiZr) 5 . The coating is composed of an aluminoborosilica glass layer on top of a RMPEA-Si-B multilayered structure. Here, the coating effectively protects the RMPEA from oxidation in high-temperature environments, as demonstrated by phase-stable operation at 10–20% higher temperatures over state-of-the-art systems without any forced-cooling system. Following an isothermal exposure at 1300 °C, the weight change of the coated sample follows a paralinear kinetics with a minor weight loss of 4.2 mg/cm 2 after 50 h. Thermal cycling tests between 1300 °C and room temperature in air resulted in the total weight gain of only 2.6 mg/cm 2 after 450 cycles. The coating shows an excellent adherence to the substrate with a boride layer acting as a barrier that maintains the coating integrity. This two-step Mo-Si-B coating method can be adapted to provide environmental resistance to a wide range of RMPEA.

36 MATERIALS SCIENCE↗

Enhancing the Sensitivity of Redox Detection in Search-for-Life Missions Using Chemical Vapor Deposition (CVD) Graphene Electrodes

The primary objective for this effort is to fabricate and characterize redox-active molecules using CVD graphene as working electrodes using standard "on-chip" lithographic techniques. We then plan to move away from “on-chip” graphene electrodes and study methods for CVD graphene transfer that can also accommodate in-space manufacturing (ISM), enabling the use of this technology for crewed mission operations of the future.

graphene↗

First Simultaneous Measurement of Differential Muon-Neutrino Charged-Current Cross Sections on Argon for Final States with and without Protons Using MicroBooNE Data

We report the first double-differential neutrino-argon cross section measurement made simultaneously for final states with and without protons for the inclusive muon neutrino charged-current interaction channel. The proton kinematics of this channel are further explored with a differential cross section measurement as a function of the leading proton’s kinetic energy that extends across the detection threshold. These measurements use data collected with the MicroBooNE detector from 6.4 × 10 20 protons on target from the Fermilab booster neutrino beam with a mean neutrino energy of ∼ 0.8 GeV . Extensive data-driven model validation utilizing the conditional constraint formalism is employed. This motivates enlarging the uncertainties with an empirical reweighting approach to minimize the possibility of extracting biased cross section results. The extracted nominal flux-averaged cross sections are compared to widely used event generator predictions revealing severe mismodeling of final states without protons for muon neutrino charged-current interactions, possibly from insufficient treatment of final state interactions. These measurements provide a wealth of new information useful for improving event generators which will enhance the sensitivity of precision measurements in neutrino experiments. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Deposition and Characterization of Thin Films on Metallic Substrates

A CVD method was successfully developed to produce conversion coatings on aluminum alloys surfaces with reproducible results with a variety of precursors. A well defined protocol to prepare the precursor solutions formulated in a previous research was extended to other additives. It was demonstrated that solutions prepared following such a protocol could be used to systematically generate protective coatings onto aluminum surfaces. Experiments with a variety of formulations revealed that a refined deposition protocol yields reproducible conversion coatings of controlled composition. A preliminary correlation between solution formulations and successful precursors was derived. Coatings were tested for adhesion properties enhancement for commercial paints. A standard testing method was followed and clear trends were identified. Only one precursors was tested systematically. Anticipated work on other precursors should allow a better characterization of the effect of intermetallics on the production of conversion/protective coatings on metals and ceramics. The significance of this work was the practical demonstration that chemical vapor deposition (CVD) techniques can be used to systematically generate protective/conversion coating on non-ferrous surfaces. In order to become an effective approach to replace chromate-based pre- treatment processes, namely in the aerospace or automobile industry, the process parameters must be defined more precisely. Moreover, the feasibility of scale-up designs necessitates a more comprehensive characterization of the fluid flow, transport phenomena, and chemical kinetics interacting in the process. Kinetic characterization showed a significantly different effect of magnesium-based precursors when compared to iron-based precursors. Future work will concentrate on refining the process through computer simulations and further experimental studies on the effect of other transition metals to induce deposition of conversion/protective films on aluminum and other metallic substrates.

Gatica, Jorge E.↗

GRMHD and GRPIC Simulations

We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code by using a conservative, high-resolution shock-capturing scheme. The numerical fluxes are calculated using the HLL approximate Riemann solver scheme. The flux-interpolated constrained transport scheme is used to maintain a divergence-free magnetic field. We have performed various 1-dimensional test problems in both special and general relativity by using several reconstruction methods and found that the new 3D GRMHD code shows substantial improvements over our previous code. The simulation results show the jet formations from a geometrically thin accretion disk near a nonrotating and a rotating black hole. We will discuss the jet properties depended on the rotation of a black hole and the magnetic field configuration including issues for future research. A General Relativistic Particle-in-Cell Code (GRPIC) has been developed using the Kerr-Schild metric. The code includes kinetic effects, and is in accordance with GRMHD code. Since the gravitational force acting on particles is extreme near black holes, there are some difficulties in numerically describing these processes. The preliminary code consists of an accretion disk and free-falling corona. Results indicate that particles are ejected from the black hole. These results are consistent with other GRMHD simulations. The GRPIC simulation results will be presented, along with some remarks and future improvements. The emission is calculated from relativistic flows in black hole systems using a fully general relativistic radiative transfer formulation, with flow structures obtained by GRMHD simulations considering thermal free-free emission and thermal synchrotron emission. Bright filament-like features protrude (visually) from the accretion disk surface, which are enhancements of synchrotron emission where the magnetic field roughly aligns with the line-of-sight in the co-moving frame. The features move back and forth as the accretion flow evolves, but their visibility and morphology are robust. We would like to extend this research using GRPIC simulations and examine a possible new mechanism for certain X-ray quasi-periodic oscillations (QPOs) observed in blackhole X-ray binaries.

Nishikawa, K.-I.↗

ALMA Observations of the Galactic Center: SiO Outflows and High Mass Star Formation Near Sgr A

Using ALMA observations of the Galactic center with a spatial resolution of 2.61" x 0.97 ", we detected 11 SiO (5-4) clumps of molecular gas in the within 0.6pc (15") of Sgr A*, interior of the 2-pc circumnuclear molecular ring. Three SiO (5-4) clumps closest to Sgr A* show the largest central velocities of approximately 150 kilometers per second and broadest asymmetric linewidths with total linewidths FWZI approximately 110-147 kilometers per second. Other clumps are distributed mainly to the NE of the ionized minispiral with narrow linewidths of FWHM approximately 11-27 kilometers per second. Using CARMA data, LVG modeling of the broad velocity clumps, the SiO (5-4) and (2-1) line ratios constrain the column density N(SiO) approximately 10(exp 14) per square centimeter, and the H2 gas density n(sub H2) = (3-9) x 10(exp 5) per cubic centimeter for an assumed kinetic temperature 100-200K. The SiO (5-4) clumps with broad and narrow linewidths are interpreted as highly embedded protostellar outflows, signifying an early stage of massive star formation near Sgr A* in the last 104 years. Additional support for the presence of YSO outflows is that the luminosities and velocity widths lie in the range detected from protostellar outflows in star forming regions in the Galaxy. Furthermore, SED modeling of stellar sources along the N arm show two YSO candidates near SiO clumps supporting in-situ star formation near Sgr A*. We discuss the nature of star formation where the gravitational potential of the black hole dominates. In particular, we suggest that external radiative pressure exerted on self-shielded molecular clouds enhance the gas density, before the gas cloud become gravitationally unstable near Sgr A*.

Yusef-Zadeh, F.↗