SEARCH · Search NASA

Results for “gas path dynamics”

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.

SourceResearch connectionFrom yearThrough year

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

18 records

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Path integral molecular dynamics: A high-fidelity approach to quantum dynamics of electrons

We investigate electron transport in the uniform electron gas using ring-polymer molecular dynamics (RPMD). Working in the weakly coupled, non-degenerate regime, we use RPMD to probe how the onset of quantum diffraction effects at high temperature reshapes electron–electron collisions and leads to a classical-to-quantum crossover in macroscopic transport properties. Static thermodynamics obtained with RPMD are consistent with the weak-coupling equation of state, confirming correct quantum Boltzmann sampling. Real-time transport extracted from mean square displacements exhibits the expected ballistic-to-diffusive transition and a systematic reduction of the electronic self-diffusivity as quantum effects strengthen, due to quantum diffraction modifying electron–electron collisions. Direct ring-polymer scattering simulations reveal diffractive “softening” of binary deflections, providing a micro-to-macro link between collision physics and diffusion. The present study establishes RPMD as a quantitative, trajectory-based tool for electron transport across the classical–quantum crossover and furnishes benchmarks for improving Coulomb-log interpolation models. We outline extensions to multi-component plasmas and a path to incorporate Fermi–Dirac statistics within path-integral dynamics.

Electronic transport↗

DOE OSTI2025

Fiber-Coupled Multipass NIR Sensor for In Situ, Real-Time Water Vapor Outgassing Monitoring

This work presents the recent development of a fiber-coupled multipass near-infrared (NIR) gas sensor used to monitor water vapor desorption of small material coupons. The gas sensor design employs a White cell topology to maximize the optical path length over a compact, hand-size footprint. Water vapor concentrations are quantified over a large dynamic range by simultaneously applying wavelength modulation and tunable diode laser absorption spectroscopy techniques. A custom headspace optimized for material desorption experiments is assembled using commercially available vacuum chamber components. We provide in situ measurements of water vapor desorption from two geometries of the industrially important silicone elastomer Sylgard-184 as a case study for sensor viability. To corroborate the results, the gas sensor data are compared to numerical simulations based on a triple-mode diffusion–sorption model, consisting of Henry, Langmuir, and Pooling modes.

gas sensor↗

DOE OSTI2025

Molecular axis distribution moments in ultrafast transient absorption spectroscopy: A path toward ultrafast quantum state tomography

In ultrafast time-resolved experiments with gas phase molecules, the alignment of the molecular axis relative to the polarization of the interacting laser pulses plays a crucial role in determining the dynamics following this light–matter interaction. The molecular axis distribution is influenced by the interacting pulses and is intrinsically linked to the electronic coherences of the excited molecules. However, in typical theoretical calculations of such interactions, the signal is either calculated for a single molecule in the molecular frame or averaged over all possible molecular orientations to compare with the experiment. Such averaging removes information about anisotropy in the molecular-axis distribution, even though anisotropic contributions can play a significant role in the measured experimental signal. Here, we calculate the laboratory frame transient electronic first-order polarization [P (1) ] spectra in terms of separated molecular frame and laboratory frame quantities. The laboratory frame polarizations are compared with orientation-averaged quantum master equation calculations, demonstrating that orientation-averaging captures only the isotropic contributions. We show that our formalism also allows us to evaluate the anisotropic contributions to the spectrum. Lastly, we discuss the application of this approach to achieve ultrafast quantum state tomography using transient absorption spectroscopy and field observables in nonlinear spectroscopy.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Epoxy-based vitrimeric semi-interpenetrating network/MXene nanocomposites for hydrogen gas barrier applications

Herein, we report MXene-filled epoxy-based vitrimeric nanocomposites featuring a semi-interpenetrating network (S-IPN) to develop a hydrogen gas (H 2 ) barrier coating with self-healing characteristics for compressed H 2 storage applications. The reversible epoxy network was formed by synthesizing linear epoxy chains with pendent bis-hydroxyl groups using amino diol, which were then crosslinked with 1,4-benzenediboronic acid to generate dynamic boronic ester linkages. To achieve the S-IPN-type molecular arrangement, the epoxy chains were in situ crosslinked in the presence of poly(ethylene-co-vinyl alcohol) (EVOH), giving rise to a self-healing network (EEP) with a healing efficiency of 87%. Into the S-IPN vitrimer (EEP), a 2D platelet-type nanofiller MXene was incorporated to introduce a tortuous path for H 2 gas diffusion along with improved mechanical properties. The nanocomposite coating was applied to nylon 6 liner material, which is conventionally used in all-composite H 2 storage vessels. Further, the application of a 2 wt% MXene/EEP nanocomposite coating showed a permeability coefficient of 0.062 cm 3 mm m -2 d -1 atm -1 exhibiting ~96% reduction in gas permeability compared to uncoated nylon 6. The same nanocomposite exhibited a healing efficiency of 79%. Increasing the MXene loading to 10 wt% further reduced the permeability coefficient to 0.002 cm 3 mm m -2 d -1 atm -1 ; however, the healing efficiency decreased due to restricted chain mobility. In essence, the current work highlights the potential of vitrimeric S-IPN nanocomposite coatings for H 2 gas-barrier applications, enhancing safety and performance.

42 ENGINEERING↗

DOE OSTI2025

Success Path Method: Introduction to the Success Path Method Software Tool©

As part of its commitment to advancing safety and reliability assessment methodologies, Argonne National Laboratory pioneered the use of an evaluation method called the Success Path Method (SPM) to improve risk management for offshore oil and gas operations. The development of the SPM at Argonne has been driven by the need to improve existing risk assessment methodologies by focusing on the steps necessary for success rather than failure modes alone. This is particularly important for industrial environments like offshore facilities that perform multiple functions under a continuously evolving set of operational conditions – such as water depth and temperature, currents, and weather conditions. In these dynamic environments, the traditional Probabilistic Risk Assessment (PRA) approach is far too complex as it focuses on what can go wrong – which comprises an infinite failure space that must be fully explored and understood. By shifting the focus to a finite space of success paths, the SPM enables operators and decision makers to prioritize a manageable number of steps that must go right to ensure success. Building on its five decades of experience in safety assessments for the nuclear industry, Argonne made major adaptations to existing risk assessment methods utilizing features similar to fault trees that are traditionally used in PRA to map all pathways in which the system can malfunction. In contrast, SPM identifies the components and processes that must function correctly to achieve specific outcomes – such as preventing the uncontrolled release of hydrocarbons during drilling operations. The SPM framework integrates equipment, procedures, software, processes, and human actions to ensure that physical barriers meet critical safety functions in dynamic operational conditions. This approach helps identify failure modes and improve operational risk management by narrowing the focus to key success elements, which in turn reduces uncertainty and helps users understand, manage, and respond to failures.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Heat Transfer Experiments of a 1st Stage Blade Cascade for Supercritical CO2 Oxy-Combustion Turbine Application

The results of internally cooled 1st stage blade (S1B) cascade testing in a supercritical CO2 environment is presented. The turbine blade design has been previously established for the end application of an oxy-combustion turbine operating in the Allam-Fetvedt cycle with turbine inlet conditions of 305 bar and 1150°C. The internally cooled blade features leading edge (LE) region impingement cooling, mid-section ribbed serpentine passages, and a pin-finned trailing edge (TE) region before cooling ejection holes. The geometry for the tested blade cascade has a cooled central blade with un-cooled blades on either side to match flowpath areas of the actual turbine. The flowpath reuses internal components previously employed for mid-section region ribbed serpentine passage experiments that established Nusselt number enhancement ratios over a range of Reynolds numbers from 100,000-400,000. New components include flow conditioning plates upstream and downstream of the blade cascade to adequately represent the flow field and blade external heat transfer coefficient profiles for the actual turbine. The cooled central blade utilizes uniform crystal temperature sensors (UCTS) with six sensors each on the blade pressure and suction surfaces distributed radially and from LE to TE. The post-processed UCTS quantified the maximum wall temperature seen at each installed sensor location. The test procedure consisted of establishing supercritical CO2 cooling flow temperature and flow rate and maintaining it throughout the test. The flow rate aims to match that for the actual in-service turbine blade design and is maintained through an orifice restriction to keep the pressure differential between internal cooling flow and external hot flow nearly constant. For the sCO2 flow path external to the blade, temperatures were ramped throughout the test via control of the test loop’s natural gas burner heater. The maximum temperature seen was 468°C and held constant for a duration of 10 minutes at which the blade metal temperature was predicted to be at its maximum before ramping down. For the turbine blade design for service inlet conditions, external flow path computational fluid dynamics (CFD) results and an internal cooling 1-D thermal and hydraulic flow network model using experimentally validated correlations served as thermal finite element (FE) boundary conditions to predict blade metal temperatures. These predicted temperatures were subsequently utilized in a structural FE model to predict blade life ratings dictated by Haynes 282 creep strength data, having a strong dependence on temperature. The boundary conditions experienced during testing are used in the same workflow and compared to the experimental results, with the goal of validating the analysis methodology and providing insight on the uncertainty in local metal temperature predictions.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2024

Design of a Three-Layer Impactor for Material Strength Measurements at Multi-Megabar Pressures

Off-Hugoniot loading methods for dynamic compression have been a topic of interest for many years as a means of accessing high pressures without melting a material. Simulations and an initial test experiment of a three-layer gas-gun impactor suitable for use in measuring dynamic material strength at multi-megabar pressures are described. As shown through ARES hydrocode simulations, an impactor consisting of a low-, medium-, and high-impedance material, such as TPX, Al, and Ta or Pt, can produce a quasi-isentropic loading path up to 4.6 Mbar pressure with a low homologous temperature (T/T melt = 0.1–0.25). These impactors can also be tailored to create a pressure drive suitable for a Lagrangian sound speed experiment. Simulated in-situ velocities were analyzed to determine the change in calculated strength due to non-uniformity in the impactor, which provides a measure of error sensitivity. An initial test of a three-layer projectile onto a flat target at 4.4 km/s impact velocity is reported. Finally, while the impactor incurred some delamination during firing, the measured sample/window interface velocity generally agrees with the simulated peak velocity and subsequent release, indicating the feasibility of this approach for measuring material strength.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Shock-induced bubble jets: a dual perspective of bubble collapse and interfacial instability theory

Interactions between shock waves and gas bubbles in a liquid can lead to bubble collapse and high-speed liquid jet formation, relevant to biomedical applications such as shock wave lithotripsy and targeted drug delivery. This study reveals a complex interplay between acceleration-induced instabilities that drive jet formation and radial accelerations causing overall bubble collapse under shock wave pressure. Using high-speed synchrotron X-ray phase contrast imaging, the dynamics of micrometre-sized air bubbles interacting with laser-induced underwater shock waves are visualised. These images offer full optical access to phase discontinuities along the X-ray path, including jet formation, its propagation inside the bubble, and penetration through the distal side. Jet formation from laser-induced shock waves is suggested to be an acceleration-driven process. A model predicting jet speed based on the perturbation growth rate of a single-mode Richtmyer–Meshkov instability shows good agreement with experimental data, despite uncertainties in the jet-driving mechanisms. The jet initially follows a linear growth phase, transitioning into a nonlinear regime as it evolves. To capture this transition, a heuristic model bridging the linear and nonlinear growth phases is introduced, also approximating jet shape as a single-mode instability, again matching experimental observations. Upon piercing the distal bubble surface, jets can entrain gas and form a toroidal secondary bubble. Linear scaling laws are identified for the pinch-off time and volume of the ejected bubble relative to the jet’s Weber number, characterising the balance of inertia and surface tension. At low speeds, jets destabilise due to capillary effects, resulting in ligament pinch-off.

Drops and Bubbles: Bubble dynamics↗

DOE OSTI2025

Numerical assessment of triply periodic minimal surfaces for direct air capture of carbon dioxide

Direct air capture (DAC) systems often consist of packing material wetted by a capture fluid that reacts with CO 2 in the airstream. The efficiency of the contactor is determined by a complex relationship of fluid dynamics, heat and mass transfer, contactor geometry, and chemical properties. The efficiency of the contactor must be balanced with other factors, primarily pressure drop through the system. Triply periodic minimal surfaces (TPMS) are a class of differential surfaces that have been explored in multiple engineering applications and have been shown to exhibit excellent performance when used in heat exchangers. Their tortuous path provides a high surface-to-volume ratio and favorable trade-off between contact area and pressure drop. In this work, a gyroid-type TPMS contactor was evaluated using computational fluid dynamics for a variety of geometric parameters to explore the potential benefit of TPMS shapes for DAC applications. A thin-film model was employed to model the flow and distribution of the capture solvent, allowing efficient simulations of TPMS structures at scale by eliminating the need for a computationally intensive interface capturing method. A liquid-gas mass transfer model was implemented in the commercial software STAR-CCM+ and used to predict the CO 2 capture efficiency and study the trade-off between capture performance and pressure drop through analysis of capture rates, mass transfer coefficients, and other relevant variables. TPMS contactors with a variety of geometric parameters and two capture solvent options were investigated to determine the effect of design choices on the operational performance of DAC systems. In conclusion, results showed that while contactor geometry is the dominant factor in efficiency and pressure drop, the physiochemical properties of the solvent are an important secondary influence on the contactor performance.

CFD↗

DOE OSTI2025

Reaction-induced departures from continuum Navier–Stokes turbulence

Reactive hydrodynamic turbulence is an inherently multiscale phenomenon, characterized by the separation between energy-containing, viscous, and molecular length and time scales. The separation between the viscous scale (the Kolmogorov scale) and the molecular mean free path ostensibly justifies a macroscopic description of reactive turbulence via the Navier–Stokes (NS) equations. However, here we use molecular-level simulations to demonstrate that exothermic bimolecular reactions can cause the NS description of turbulence to break down in the near-continuum regime. Sufficiently energetic heat-releasing reactive collisions strongly distort the Maxwell–Boltzmann velocity distribution function, modifying not only the macroscopic chemical rate law but the kinetic-energy-transfer processes as well. This translational nonequilibrium ultimately introduces significant departures from the NS description in the kinetic energy spectra at scales orders of magnitude larger than both the molecular mean free path and the Kolmogorov length scale. These departures prove substantial enough to meaningfully alter integrated quantities, including the overall turbulence kinetic energy itself.

molecular gas dynamics↗

DOE OSTI2024

Cosmic-Ray Feedback on Bistable Interstellar Medium Turbulence

Abstract While cosmic rays ( E ≳ 1 GeV) are well coupled to a galaxy’s interstellar medium (ISM) at scales of L > 100 pc, adjusting stratification and driving outflows, their impact on small scales is less clear. Based on calculations of the cosmic-ray diffusion coefficient from observations of the grammage in the Milky Way, cosmic rays have little time to dynamically impact the ISM on those small scales. Using numerical simulations, we explore how more complex cosmic-ray transport could allow cosmic rays to couple to the ISM on small scales. We create a two-zone model of cosmic-ray transport, with the cosmic-ray diffusion coefficient set at the estimated Milky Way value in cold gas but smaller in warm gas. We compare this model to simulations with a constant diffusion coefficient. Quicker diffusion through cold gas allows more cold gas to form compared to a simulation with a constant, small diffusion coefficient. However, slower diffusion in warm gas allows cosmic rays to take energy from the turbulent cascade anisotropically. This cosmic-ray energization comes at the expense of turbulent energy which would otherwise be lost during radiative cooling. Finally, we show our two-zone model is capable of matching observational estimates of the grammage for some transport paths through the simulation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Wakefield generation in hydrogen and lithium plasmas at FACET-II: Diagnostics and first beam-plasma interaction results

Plasma wakefield acceleration provides ultrahigh acceleration gradients of tens of GeV/m, providing a novel path toward efficient, compact, TeV-scale linear colliders, and high brightness free electron lasers. Critical to the success of these applications is demonstrating simultaneously high gradient acceleration, high energy transfer efficiency, and preservation of emittance, charge, and energy spread. Experiments at the FACET-II National User Facility at SLAC National Accelerator Laboratory aim to achieve all of these milestones in a single-stage plasma wakefield accelerator, providing a 10 GeV energy gain in a < 1 m plasma with high energy transfer efficiency. Such a demonstration depends critically on diagnostics able to measure emittance with mm mrad accuracy, energy spectra to determine both percent level energy spread, and broadband energy gain and loss, incoming longitudinal phase space, and matching dynamics. This paper discusses the experimental setup at FACET-II, including the incoming beam parameters from the FACET-II linac, plasma sources, and diagnostics developed to meet this challenge. Initial progress on the generation of beam ionized wakes in meter-scale hydrogen gas is discussed as well as commissioning of the plasma sources and diagnostics. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Chemistry of Sugar Formation in the Gas Phase: Following the Activated Aldehyde

Sugars are produced by living organisms, and are required building blocks for life as we know it, which raises the foundational question of how sugars formed in a prebiotic environment. The abiotic formose reaction produces sugars from formaldehyde, but our understanding of its initiation step remains murky, with chemists invoking the concept of an “activated aldehyde” to seed this reaction. Singlet hydroxycarbenes, high-energy isomers of aldehydes, were recently reported to facilitate sugar formation under cold, nonaqueous conditions relevant to interstellar environments. Here, we generate singlet methylhydroxycarbene ( 1 CH 3 –C̈–OH) from the photodissociation of pyruvic acid and experimentally measure its gas-phase reaction with d 4 -acetaldehyde using multiplexed photoionization mass spectrometry. The C 4 H 4 D 4 O 2 isomer d 4 -acetoin is the sole product, which we kinetically link to the reactant CH 3 –C̈–OH, and attribute to a carbonyl-ene formation mechanism. We see no evidence of 3-hydroxybutanal, the C–H insertion product expected in carbene chemistry. Using automated exploration we calculate stationary points on the potential energy surface and report master equation rate coefficients from T = 20–600 K, providing quantitative kinetics of this fast reaction for use in chemical models. The prereactive complex in this reaction is stabilized by both hydrogen bonding and electrophilic carbene C═O interactions. These effects create a short-range dynamical bottleneck for the reaction besides the long- and midrange barrierless bottlenecks. Combined with recent reports of 1 HC̈OH production from methanol photodissociation and pyruvic acid production in cold irradiated ices, this work provides evidence that singlet hydroxycarbene + aldehyde chemistry is a feasible path to prebiotic sugar formation.

aldehydes↗

DOE OSTI2025

Evolution of intermittent filaments in the scrape-off layer of NSTX

Filamentary structures naturally arise from background turbulence in the scrape-off layer (SOL) of plasmas, leading to significant particle and heat transport that can degrade overall plasma confinement. This enhanced transport can contribute to unacceptably high heat loads on plasma-facing components. As such, understanding the physics of SOL plasma filaments is critical for predicting and mitigating their effects in future fusion devices. On the National Spherical Torus Experiment (NSTX), plasma filaments—commonly referred to as “blobs”—were investigated using the gas-puff imaging (GPI) diagnostic in the edge and SOL regions. The analysis involved identifying, segmenting, and tracking the characteristic contours of the blobs in each frame of the GPI video sequences. Their evolution was characterized through shape descriptors, velocity, and angular velocity derived from their contour coordinates. The results indicate that as the blob area increases, their shapes become more concave and less circular, suggesting reduced structural stability in larger blobs. This result aligns with previous theoretical results where it was shown that larger blobs are more susceptible to instabilities [Krasheninnikov et al., J. Plasma Phys. 74, 679–717 (2008) and D'Ippolito et al., Phys. Plasmas 18, 060501 (2011)]. A positive correlation was observed between radial velocity and radial position, suggesting radially outward acceleration of the filaments, potentially driven by decreasing viscous drag toward the far SOL. Interestingly, blobs in background SOL turbulence exhibited minimal spinning in contrast to filaments originating from edge localized modes, which show substantial rotation during their paths [Lampert et al., Phys. Plasmas 29, 102502 (2022)]. Statistical analysis of the solidity and total curvature shape descriptors, along with their temporal evolution, revealed relatively broad, near-Gaussian distributions. This suggests that blob morphology is strongly influenced by stochastic turbulent processes in the surrounding plasma environment. Blob parameters were also compared with bulk plasma and radial profile measurements. Notable trends were found between blob rotation and poloidal velocity with collisionality and line-integrated density. These findings contribute to a deeper understanding of blob dynamics and provide valuable insights for refining SOL turbulence models.

Covariance and correlation↗

DOE OSTI2024

Thermodynamic stability of a spin microemulsion in Rashba spin-orbit-coupled bosons

Recent finite-temperature numerical simulations have unveiled a quantum “spin” microemulsion analog, found by raising the temperature of a stripe supersolid phase in a Rashba spin-orbit coupled Bose gas. This microemulsion state is a highly correlated, isotropic normal fluid where atoms self-arrange based on their internal pseudospin into patterns that resemble bicontinuous microemulsions. This finding leaves several open questions regarding the broader accessibility of this phase in experiments. Here, we use equilibrium finite-temperature numerical simulations based on a coherent-state path integral representation to perform a computational investigation into the thermodynamic stability of the spin microemulsion state. Numerical simulations emphasize the requirement of a nearly, but not perfectly, isotropic spin-orbit coupling in order to achieve the microemulsion phase in cold-atom experiments. Moreover, the microemulsion state exists independent of miscibility of the pseudospin components and for a wide range of pseudospin population imbalance, suggesting a high degree of flexibility in choosing the atom and hyperfine states in an experimental realization. Lastly, we demonstrate this feasibility by mimicking a Rashba spin-orbit-coupled 87 Rb experiment in an isotropic harmonic trap, where we confirm the microemulsion's existence via its density profile and equilibrium quasimomentum distribution.

Complex Langevin dynamics↗

DOE OSTI2026

Evolution of size-selected Pt cluster catalysts on prototypical oxide supports

The current quest for new pathways into sustainable, efficient and durable energy conversion technologies makes the used for a fundamental understanding of the atomic-scale phenomena underlying catalytic processes ever more pressing. In this context, characterizing catalyst particles in situ provides valuable information about the evolution of their composition, structure, oxidation state and charge state during an ongoing process. To disentangle the influence of individual parameters – temperature, pressure, gas composition, cluster size, as well as support acidity, redox state and defect density – it is crucial to control them precisely and separately in experiments. At the example of size-selected Pt n clusters – i.e. sub-nm particles defined to the exact number of atoms – on flat oxide supports, we follow their rich evolution phenomena via (synchrotron-based) X-ray photoelectron spectroscopy (XPS) and scanning tunnelling microscopy (STM) during temperature ramps and in various gas environments. Here, we present our experience with these highly defined, yet complicated-to-create samples in ultra-high vacuum (UHV) and at mbar pressures. We discuss their stability during transport to synchrotrons and under reaction conditions on three prototypical oxide supports and show the various phenomena that can be disentangled. Ptn clusters on the non-reducible silicon dioxide, SiO 2 , remain size-selected and show remarkable stability. Doping strongly influences the size-dependent binding energy shifts and changes the Pt response to oxidative and reaction conditions, which we attribute to different cluster geometries: p-type doping leads to wetting, enhanced sinter resistance and a diminished response to oxidative environments, compared to clusters on n-type doped samples. We compare the system with our previous findings of a similar change in dimensionality for Pt 20 clusters on the reducible ceria, CeO 2 (111), support, induced by modulation of its O vacancy density. The Pt n /CeO 2 system is particularly interesting for strategies to stabilize and redisperse Pt clusters dynamically. Finally, we study the evolution of Pt n clusters on another reducible oxide support, magnetite, Fe 3 O 4 (001), in 0.1 mbar alternating redox conditions at RT and elevated temperature. In analogy to findings previously reported for Pt/TiO 2 (110), the clusters either become encapsulated by a thin oxide film via strong metal–support interaction (SMSI) or deeply buried in the magnetite. Overall, our approach of following the evolution of size-selected clusters on oxide supports leads to fundamental atomic-scale insights on nano-scale catalyst materials, on our path to sustainable, dynamic and self-repairing catalysts.

Falling, Lorenz J. [Technical Univ. of Munich (Ger↗

DOE OSTI2025

The Multiple Paths to Merger of Unequal-mass Black Hole Binaries in the Disks of Active Galactic Nuclei

The accretion disks that power active galactic nuclei (AGN) are thought to house populations of stars and compact objects; after forming binaries these compact objects may merge, begetting gravitational waves (GWs) such as those detected by LIGO and VIRGO. We present a comprehensive study of the early evolution of binaries within AGN disks as their orbits are influenced by the surrounding gas, focusing on eccentric and unequal-mass binaries. Nearly equal-mass binaries behave similarly to their equal-mass counterparts: Prograde binaries inspiral, albeit somewhat slowly, and have their eccentricities damped; retrograde binaries inspiral ∼2–3 times faster than their prograde counterparts, and those with near-equal masses are driven quickly toward near-unity eccentricities. However, the primaries in retrograde binaries with mass ratios of m 2 /m 1 ≲ 0.4 experience significantly weaker headwinds and retain substantial accretion disks that help damp binary eccentricities, slowing binary inspirals. Additionally, we find that while accretion drives prograde binaries toward equal masses thanks to the exchange of material between the primary and secondary accretion disks, retrograde binaries are driven slowly toward more extreme mass ratios. Prograde binaries, and generally those with low mass ratios, likely accrete for multiple e-folding timescales before merger. On the other hand, high-mass-ratio retrograde binaries may merge before accreting substantially, potentially approaching merger with detectable eccentricity. Future ground-based GW observatories, with their broader frequency coverage, should be particularly useful for studying these populations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Development, Verification, and Validation of an OpenFOAM-Based Solver for Modeling Inertial Fusion Energy Chambers

Our work seeks to introduce a computational tool tailored to the physics of inertial fusion energy chambers, in particular, those concepts based on thick liquid walls. In this approach, the structural materials are protected by several neutron mean-free-paths of renewable liquid and thus will be able to survive much longer than un-shielded walls, with virtually all structures lasting for the life of the plant and enabling the use of commercially available and qualified materials. The OpenFOAM-based solver named rhoCentralFoam has been used as a starting point. rhoCentralFoam belongs to the standard OpenFOAM solver toolset. It is a high-speed, explicit compressible flow solver with shock-capturing capability. While the main features have been retained, the solver had to be restructured to make use of tabular data for equations of states, a necessary addition to model the complex thermo-physical properties of ionized gasses. This entailed the need to change the independent state variables used by the solver, resulting in a new thermodynamic library and slightly different solution algorithm. Moreover, a radiation heat transfer model based on the P-1 approximation was added to the solver. The solver is verified against an analytical solution from the Sedov-Taylor-Neumann test problem to showcase the ability of the hydrodynamic solvers to handle strong shocks, whereas the P-1 model was verified using a simple one-dimensional problem with an analytical solution. Additionally, a validation case involving shock-wave propagation through jet array is presented, and the results are compared with experimental data from the open literature. Lastly, in order to showcase the utility of the solver for practical cases, we applied the refined solver to two representative scenarios: gas venting within the HYLIFE-II chamber and the compression of the gas following the partial ablation of the liquid wall.

Chamber dynamics↗