Size and oxidation state tracking of dynamic Rh catalysts on rutile TiO2 by ambient-pressure XPS
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
In the present work, a first-of-its-kind three-dimensional (3D) large-eddy simulation (LES) study is conducted to numerically investigate the combustion dynamics as well as aero-thermal phenomena in a full-scale nonpremixed hydrogen–air rotating detonation engine (RDE) (with a diverging-shaped lower-end wall), when integrated with nozzle guide vanes (NGV) acting as the turbine stator. The wall-modeled LES framework incorporates hydrogen–air detailed chemical kinetics and adaptive mesh refinement (AMR). A comparative analysis is carried out for two operating conditions with different fuel/air mass flow rates but global equivalence ratio of unity, and considering RDE configurations without and with stator. The LES model is validated against available experimental data for the low mass flux condition with respect to detonation wave speed/height, wave dynamics, and axial static pressure distribution. Numerical results indicate significant deflagrative combustion occurring in the fill region near the inner wall due to formation of recirculation zones in the injection near-field driven by the backward facing step. The leading detonation wave is found to be trailed by an azimuthal reflected-shock combustion (ARSC) wave, consistent with experimental observations, which consumes unburned vitiated reactants that leak through the main detonation wave. The main detonation wave characteristics, such as detonation wave speed/height and combustion efficiency, do not change appreciably with the presence of NGV. A novel combustion diagnostic technique based on chemical explosive mode analysis (CEMA) is employed to quantify the fraction of heat release occurring in the detonative mode versus deflagrative mode for the simulated conditions. The exit flow is found to be nearly fully subsonic and supersonic for the low and high mass flux conditions, respectively. Further analysis of the exit flow profiles shows that the presence of NGV renders the flow more axial and significantly impacts the exit Mach number and total pressure, while the total temperature shows negligible change. In addition, the low mass flux operating point, despite exhibiting more deflagrative losses within the combustor, yields overall lower pressure drop from plenum to exhaust, which is mainly attributed to lower pressure drop across the injectors. Lastly, the rotating detonation engine-nozzle guide vanes (RDE-NGV) configuration exhibits higher total pressure loss compared to rotating detonation engine (RDE) without stator across both the mass flux conditions. In conclusion, this study extends the state-of-the-art in numerical modeling of pressure gain combustion (PGC) systems by demonstrating high-fidelity 3D reacting LES of full-scale RDE-NGV systems relevant to RDE-turbine integration for stationary power generation.
Fission gas bubbles in UO 2 nuclear fuel have been observed to exhibit pressures in excess of the equilibrium bubble pressure; however, the cause of bubble over-pressurization has not yet been demonstrated. The mechanical interaction between a bubble and the surrounding matrix or grain boundary depends on the internal pressure of the bubble and local stress state, such that over-pressurized bubbles are thought to be responsible for fragmentation and pulverization, when exposed to a temperature ramp. Here, in this work, we investigate the role of U interstitials, produced through irradiation, in over-pressurizing bubbles by using a combined molecular dynamics (MD) and cluster dynamics approach. Firstly, the energies for the capture of interstitials and vacancies by bubbles have been determined from MD as a function of the ratio of gas atoms to vacancies that make up the bubble. Secondly, these reaction energies have been implemented in the cluster dynamics code Centipede to predict bubble over-pressurization as a function of temperature for typical fission rates. It was found that there is a transition from low pressure bubbles (at high temperatures) to high pressure bubbles (at lower temperatures). The cause of this behavior was shown to be the creation of irradiation-induced interstitials that are highly mobile relative to vacancies at low temperature; whereas, vacancies are sufficiently mobile at high temperatures to limit bubble pressures. This result supports the hypothesis that over-pressurized bubbles form during steady-state operation and that this behavior is highly sensitive to the local pellet temperature.
Passive ambient noise monitoring is an emerging tool in environmental seismology, leveraging the ambient seismic field to assess temporal variations in shallow subsurface properties. This review focuses on the potential and challenges of using scattered coda waves from noise correlation functions to monitor critical zone dynamics. The sensitivity of seismic velocities to various environmental factors, including precipitation, snowmelt, atmospheric pressure, and groundwater fluctuations, underscores the method’s versatility. While coda waves excel in detecting subtle changes due to their scattered nature, ballistic waves provide higher spatial resolution, albeit with challenges in source stability. Advances in seismic sensing, including distributed acoustic sensing and low-cost geophone networks, have enabled high-resolution monitoring of hydrological processes, subsurface deformation, and seismic hazards. Integrating seismic data with hydrological models provides insights into water storage, pore pressure changes, and soil moisture dynamics. However, limitations in spatial resolution, calibration with ground truth data, and coupled effects between environmental factors remain key challenges. This review emphasizes the importance of interdisciplinary approaches in refining methodologies, enhancing sensor deployments, and addressing data gaps. Passive seismic monitoring offers opportunities to understand critical zone processes and their broader impacts on seismic hazards and environmental sustainability.
Integrating monitoring data to efficiently update reservoir pressure and CO 2 plume distribution forecasts presents a significant challenge in geological carbon storage (GCS) applications. Inverse modeling techniques are commonly used to fuse observational data and refine reservoir model parameters, thereby improving state variable forecasts. However, these techniques often rely on linear or Gaussian assumptions, which can limit their effectiveness in accurately predicting state variables. Moreover, simulating large-scale three-dimensional (3D) GCS problems is computationally expensive, making iterative runs in inverse problems prohibitive. To address these challenges, we propose a conditional generative model utilizing the score-based diffusion method for real-time 3D pressure and saturation field distribution predictions. Our approach involves solving the score function with a mini-batch-based Monte Carlo estimator to generate labeled data. This data is subsequently employed to train a fully connected neural network, enabling it to learn the conditional sample generator within a supervised learning framework. This method enables the rapid generation of a large ensemble of predictions, facilitating comprehensive uncertainty quantification of state variables. Here we applied our method to forecast the dynamic 3D distributions of pressure and saturation fields over a 30-year injection period. The statistical assessment with low root mean square error (RMSE) values demonstrates that our method can accurately predict the spatiotemporal distributions of both pressure and saturation fields. Moreover, the developed conditional generative model shows high computational efficiency by generating 100 ensemble forecasts of 3D state variables in less than 10 min. The consistency between ensemble averages and ground truth values further illustrates the model’s capability to capture state variable dynamics during the CO 2 plume injection process. Notably, the ground truth values fall within the ensemble forecasts, indicating that our uncertainty quantification effectively captures variability and potential noise in the observations. Thus, the developed conditional generative model proves to be a more efficient, accurate, and practical tool for GCS applications, facilitating timely risk analysis and informed decision-making.
Developing advanced low-temperature combustion engines with ammonia-biofuel blends requires a comprehensive understanding of low-temperature flame dynamics and kinetic interactions between ammonia and oxygenated fuels at elevated pressures. This work aims to study the dynamics and kinetics of non-premixed Dimethyl Ether (DME)/Ammonia (NH 3 ) cool and warm flames, and their reignition to hot flames. A counterflow burner is employed to establish DME/NH 3 cool/warm flames at pressures up to 5 atm. The extinction limits of cool flame and the reignition limits of warm flame to hot flame are measured by varying NH 3 concentrations and compared to simulations to quantitatively examine the effects on DME/NH 3 flames. It is found that NH 3 inhibits low-temperature DME oxidation and results in lower cool flame extinction limits. Warm flames in the presence of NH 3 are observed for the first time, revealing a non-monotonic effect of NH 3 addition: a small amount of NH 3 presence enhances warm flame chemistry and promotes reignition to hot flames, while a high NH 3 concentration weakens the warm flame. This trend is further explained by 0-D PSR kinetic simulations and 1-D S-curve flame dynamic calculations. Three flame transition regimes between cool flames (CF), warm flames (WF), and hot flames (HF) by different levels of NH 3 additions at a specific strain rate are identified, namely WFHF reignition, WF-CF transition, and WF extinction. Reaction sensitivity analyses of OH at low temperatures show that NH3 inhibits DME oxidation through OH consumption via H-abstraction and the kinetic couplings of RO 2 /NH 2 , RO 2 /NO x , R/NO x , and O 2 QOOH/NO x further suppress the low-temperature branching. At intermediate-temperatures, NH 2 /NO x /HO 2 coupling promotes warm flames via the pathway NH 2 → H 2 NO → HNO → NO by converting O 2 → HO 2 → OH. At even higher NH₃ concentrations, radical termination reactions of NH 2 + NO/NO 2 and excessive OH consumption via H-abstraction inhibit the flame. The insights into the kinetic coupling between NH 3 and low-temperature chemistry at elevated pressure and its impact on the dynamics of cool-warm-hot flame transitions will contribute to advancing combustion technologies with reduced emissions and improved energy-efficiency.
The high-pressure, temperature phase diagram and transport properties of materials are of broad interest to planetary sciences and high-energy-density sciences and applications. Alumina (Al 2 O 3 ) or its various forms (e.g., solutions with other oxides or silicates) are important constituents in Earth’s and super-Earths’ mantles, common window materials in dynamic compression experiments, and standard pressure calibrators in diamond-anvil-cell experiments. Its structures and transport properties are of particular importance but have not been well studied at above 100 GPa pressures. Based on extensive first-principles molecular dynamics calculations, we obtain atomic level insights on structural differences and accurate results on electrical conductivity (σ dc ), thermal conductivity (κ), and reflectivity of Al 2 O 3 in several phases at pressures of ~130–1300 GPa and temperatures of 4000–20,000 K. We find the solid-to-liquid state changes of Al 2 O 3 are accompanied by an insulator-to-semimetal transition with enhanced σ dc and κ, which is similar to other silicates and oxides and can facilitate magnetic field generation in deep interiors of rocky planets. We also find the solid-to-solid transition (from Rh 2 O 3 (II) to CaIrO 3 ) is accompanied by negligible changes in σ dc and κ at 160 GPa, contrary to previous expectations, whereas another higher pressure phase transition (CaIrO 3 to U 2 S 3 at 450 GPa) is associated with increases in σ dc and κ. Furthermore, we show the transport properties can be enhanced by oxygen vacancy in Al 2 O 3 , but are not significantly affected by forming solutions with MgSiO 3 , particularly when in the solid states.
Pressure solution is an important process in the evolution of sedimentary rocks, which provide storage space for most of our petroleum resources. It directly influences the generation, migration, and storage of petroleum fluids in subsurface sedimentary rocks. Here, in this paper, we develop a pore-scale, mechanochemical model to demonstrate a possible positive feedback between the local porosity and pore surface stress, in which a higher local porosity causes a higher local pore surface stress, thus enhancing pressure solution and consequently further increasing the local porosity. Pore surface stress represents stress on a solid grain adjacent to a pore. Specifically, the pore-scale, mechanochemical model directly simulates the stress distribution over solid and pore surfaces using a finite element model. The dissolution of solids at the solid-pore interfaces under a far-from-equilibrium condition is simulated using a first-order kinetics model that accounts for the local stress distribution. The updated pore geometry, caused by pore surface dissolution, is then used in the stress simulation in the next numerical iteration. Two types of porous media, the Oriskany sandstone and an artificial porous medium with spherical pores, were tested in the mechanochemical simulation. The positive stress-porosity feedback during pressure solution was observed in both samples. In addition, the model quantitatively illustrated the distribution of local mineral dissolution rates on all pore surfaces, as well as its relation to the effective mineral dissolution rate of the entire sample. Based on the comparison between the two porous media, the local mineral dissolution was regulated by pore space distribution, geometry, and coalescence during pressure solution. This work is the first that uses direct, pore-scale numerical simulation to demonstrate the positive stress-porosity feedback during pressure solution, which has the potential to advance the understanding of the mechanical-chemical (MC) coupling in many geological processes that are relevant to subsurface energy systems, such as the recovery of petroleum hydrocarbons and geothermal energy.
Abstract Fractures and faults represent planes of weakness and compliance in rock masses that serve as focal points for both microearthquakes and fluid transport, with seismicity and permeability evolution closely linked. Contact stiffness is highly stress‐sensitive and directly influences permeability. We explore the co‐evolution of specific stiffness and permeability of rough fractures under normal stress and shear offset using numerical simulations. Individual rough fractures are represented by variable amplitude (Root mean square) and wavelength ( λ ) using a granular mechanics model. Contacting rough surfaces are mated, offset in shear, and then compacted in displacement mode. The compacting fractures generate stress‐dependent changes in contact porosity, which govern both permeability and stiffness evolution. We establish a universal dimensionless relationship linking specific stiffness and permeability that inherently incorporates the effects of surface roughness, shear offset, and microcracking. The observed cracking effect—where local stress redistribution and pressure‐driven microcrack propagation dynamically alter the aperture field—introduces a nonlinear permeability response at high stress. Increased roughness amplitude and larger shear offsets reduce stiffness while dampening permeability sensitivity to stress, demonstrating a strong interplay between surface texture and hydro‐mechanical behavior. While the model captures this behavior effectively, deviations emerge at very low porosities due to extreme aperture sensitivity in this limit.
Homogeneous nucleation, a textbook transition path for phase transitions, is typically understood on thermodynamic grounds through the prism of classical nucleation theory. However, recent studies have suggested the applicability of classical nucleation theory to systems far from equilibrium. Here, in this article, we formulate a purely mechanical perspective of homogeneous nucleation and growth, elucidating the criteria for the properties of a critical nucleus without appealing to equilibrium notions. Applying this theory to active fluids undergoing motility-induced phase separation, we find that nucleation proceeds in a qualitatively similar fashion to equilibrium systems, with concepts such as the Gibbs–Thomson effect and nucleation barriers remaining valid. We further demonstrate that the recovery of such concepts allows us to extend classical theories of nucleation rates and coarsening dynamics to active systems upon using the mechanically derived definitions of the nucleation barrier and surface tensions. Three distinct surface tensions—the mechanical, capillary, and Ostwald tensions—play a central role in our theory. While these three surface tensions are identical in equilibrium, our work highlights the distinctive role of each tension in the stability of active interfaces and the nucleation and growth of motility-induced phases.
Tungsten (W) is the leading candidate for divertor target plates because of its high melting point (>3000°C), thermal conductivity, and ultimate tensile stress. While W and its alloys are the only solid materials that can survive the high heat fluxes incident on the divertor, W’s low-ductility high ductile-to-brittle transition temperature of ~600°C and relatively low recrystallization temperature (RT) of ~1300°C pose structural (among other) challenges. The objective of this work is to estimate the thermal-fluid and thermal-structural performance of the helium (He)-cooled T-tube divertor, which was originally developed by the Advanced Reactor Innovation and Evaluation Study (ARIES) using numerical simulations. Here, predictions of temperature distributions across the plasma-facing structural component and surface pressures from computational fluid dynamics simulations are used to determine stress distributions using commercial structural finite element modeling software over a range of fusion-relevant conditions. The maximum allowable incident heat fluxes are determined based on the temperature limits imposed by the ITER elastic Structural Design Criteria for In-vessel Components (SDC-IC) and the maximum RT over a range of He mass flow rates and presented in the form of performance design charts. Our recent work found that thermal- structural criteria accounting for the low ductility of W in a finger-type modular divertor constrain the maximum incident heat fluxes to values well below the ITER specifications, and those based on considering only the RT demonstrate that integrated thermal-fluid and elastic structural performance evaluation are required for accurate assessment of divertor performance. This novel analysis of the T-tube considers how nonuniform and transient incident heat fluxes affect its thermal-fluid and thermal-structural performance, as well as the effect of volumetric heating, which can be as great as 27% of the power incident on the divertor surface. The W tile of the T-tube, with its relatively large plasma-facing area of ~15 cm 2 , will likely experience significant spatial variations in incident heat flux. This work therefore assesses whether steady-state incident heat flux profiles with a peak of 10 MW/m 2 and maximum heat flux gradients of 200 MW/m 2 per m exceed the structural limits imposed by the ITER elastic SDC-IC and the maximum RT over a range of fusion-relevant conditions. The effect of transient heat fluxes typical of plasma detachment and reattachment from the target plate due, for example, to gas injection are also evaluated
The electric-field induced second harmonic (E-FISH) signal sensing is commonly used as a nonlinear optical technique to probe electric fields inside a plasma region. Cross-beam E-FISH is being investigated to improve spatial resolution by defining the interaction volume via a controlled geometry of two overlapping noncolinear optical beams. This drastic reduction in interaction length of the electric field and laser region results in a significant signal reduction. Here, to overcome this signal reduction, we introduce coherent amplification of the cross-beam E-FISH signal by mixing the low E-FISH signal with a phase-locked bright local oscillator. We demonstrate enhancement of the signal. By introducing a local oscillator we can now derive the polarity of the measured electric field through the phase of the homodyne signal. To illustrate the technique, we, for the first time, to the best of our knowledge, measure the magnitude and the direction of the electric field in a cold atmospheric pressure plasma jet, which dynamically follows the profile of the applied bias current.
High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large eddy simulations due to their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reactive Euler equations encountered in high-speed combustion. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of DG approach. Thus, the framework is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. The numerical method is implemented within the spectral element solver Nek5000. Validation cases are conducted for both non-reactive and reactive discontinuous flows to demonstrate the solver capability. In particular, canonical one-dimensional and two-dimensional detonation simulations are performed and the high-order numerical results are validated against available literature data.
Idaho National Labs serves as the spearhead for many innovative energy solutions to the world's energy crisis. One such solution is the INL's Microreactor which is designed to deploy to extreme/remote environments where other sources of power are either unavailable or unreliable. In order to best design these energy solutions for their operational environments, it is crucial to understand how the design, components, and materials will respond to the environmental conditions. One key material in these innovative designs is a nuclear-grade graphite known as PCEA. This study examines the behavior of PCEA graphite under dynamic loading, similar to that which may occur in extreme environments. The objective is to characterize the dynamic behavior and produce an accurate, reliable constitutive model suitable for use in simulation tools such as INL's MOOSE. Graphite specimens were tested using a Split Hopkinson Pressure Bar (SHPB) to administer the dynamic compressive load. The SHPB was charged at various pressures to produce a range of strain rates on the material in compression. Data was acquired via strain gauges on the SHPB setup, from which stress, strain, and time data were collected. Analysis revealed the stress-strain behavior of the material as well as insights into the material behavior's relationship to strain rate. Further work must continue to characterize the various other dynamic behaviors of the material which will combine to create a substantially trustworthy constitutive model for this grade of nuclear-grade graphite. Ultimately, this will allow for realistic simulation of the material in reactor designs, allowing for prediction of design weaknesses and leading to improved designs for increased resilience, security, and reliability.
The pressure-dependent and initial-density-dependent dynamic strengths of three brittle powders, SiC, B4C, and CaF2, are determined via the Tamped Richtmeyer-Meshkov instability (RMI) method.
Density-functional theory based molecular-dynamics simulations were used to investigate high-pressure chemical reactions in liquid mixtures of CO 2 with several elements (Si, Mn, and Fe) at high temperatures of 2000-3000 K. Our ab initio simulations indicate that these reactant elements can reduce CO 2 to C at high pressures (20 GPa) leading to the formation of C-C chains, with Si by far the most effective carbon-reducing agent. A combined chemical analysis using Bader charge analysis and Crystal Orbital Hamilton Population (COHP) on simulation snapshots shows that significant charge transfer from the reducing element to the C atoms creates instability in the C-O covalent bonds. COHP analysis further shows that Mn/Fe-O and Mn/Fe-C bonding interactions are weaker compared to the Si counterparts. These results further our understanding of the redox chemistry of CO 2 at conditions relevant to planetary mantle interiors and demonstrate the effectiveness of high pressure in the reduction of CO 2 directly to solid carbon.
Phase transitions are omnipresent in modern condensed matter physics and its applications. In solids, first-order phase transformations typically occur by nucleation and growth under nonequilibrium conditions. Under constant external conditions, e.g., constant annealing temperature and pressure, the nucleation and growth dynamics are often thought of as spatially and temporally independent. Here, in situ Bragg X-ray photon correlation spectroscopy (XPCS) reveals nanoscale spatial and dynamical heterogeneity in the perovskite-to-brownmillerite topotactic phase transformation in La 0.7 Sr 0.3 CoO 3 thin films annealed under constant reducing conditions over a time span of multiple hours. Specifically, a time scale associated with domain growth remains stable, with a corresponding domain wall speed of v d = 6 ± 0.5 × 10 –4 nm/s (2 ± 0.2 nm/h), while a slower time scale, associated with temperature-driven depinning of domains, leads to accelerating dynamics with time scales following an aging power law with exponent −2.2 ± 0.5. This experiment demonstrates that Bragg XPCS is a powerful tool to study nanoscale dynamics in structural phase transformations, with the ability to extract quantitative average values related to nanodomain motion in situ. Furthermore, the results are relevant for phase engineering of phase-change devices, as they show that nanoscale dynamics, linked to domain and domain-wall motion, can continuously evolve and speed up with time, even hours after the initiation of the phase transformation, with potential repercussions on electrical performance.
The influence of fuel composition on self-excited combustion instabilities in a high-pressure combustor operated with ammonia, hydrogen, and natural gas fuels is characterized with high-frequency pressure measurements and imaging of the flame structure. A micromix multi-stage injector with 19 elements is used to introduce the fuel blend premixed with heated air into an optically accessible combustor operated at ~1.1 MPa. As hydrogen is substituted for natural gas, longitudinal thermoacoustic instabilities are observed in the combustor with pressure fluctuation amplitudes as large as 8% of the mean chamber value. In the absence of natural gas, limit-cycle instability magnitude is largely insensitive to ammonia addition. However, fuel compositions with > 20% natural gas result in positive correlation between instability amplitude and ammonia concentration. Spatial distribution of heat release in the combustor evaluated from OH* chemiluminescence imaging reveals axial growth in the regions of heat release fluctuation as hydrogen decreases, thus correlating flame length to instability amplitude. Distinct transitions from the first harmonic of the fundamental longitudinal acoustic mode (~1100 Hz) to the fundamental mode (~550 Hz) are observed with root mean square pressure fluctuations exceeding 7% of the mean chamber value. Here, these cases correspond to hydrogen mole fractions ≤ 50%. Analysis of the phase relationship between pressure fluctuations upstream and downstream of the flame zone indicates acoustic coupling of the injector as a key contributor to instability growth.