Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase velocity”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Seismic Anisotropy of Temperate Ice in Polar Ice Sheets

Here we present a series of simple shear numerical simulations of dynamic recrystallization of two-phase nonlinear viscous materials that represent temperate ice. First, we investigate the effect of the presence of water on the resulting microstructures and, second, how water influences on P wave (V p ) and fast S wave (V s ) velocities. Regardless the water percentage, all simulations evolve from a random fabric to a vertical single maximum. For a purely solid aggregate, the highest V p quickly aligns with the maximum c-axis orientation. At the same time, the maximum c-axis development reduces V s in this orientation. When water is present, the developed maximum c-axis orientation is less intense, which results in lower V p and V s . At high percentage of water, V p does not align with the maximum c-axis orientation. If the bulk modulus of ice is assumed for the water phase (i.e., implying that water is at high pressure), we find a remarkable decrease of V s while V p remains close to the value for purely solid ice. These results suggest that the decrease in V s observed at the base of the ice sheets could be explained by the presence of water at elevated pressure, which would reside in isolated pockets at grain triple junctions. Under these conditions water would not favor sliding between ice grains. However, if we consider that deformation dominates over recrystallization, water pockets get continuously stretched, allowing water films to be located at grain boundaries. This configuration would modify and even overprint the maximum c-axis-dependent orientation and the magnitude of seismic anisotropy.

58 GEOSCIENCES↗
DOE OSTI2026

Phase-field model of freeze casting

Directional solidification of water-based solutions has emerged as a versatile technique for templating hierarchical porous materials. However, the underlying mechanisms of pattern formation remain incompletely understood. In this work, we present a detailed derivation and analysis of a quantitative phase-field model for simulating this nonequilibrium process. The phase-field model extends the thin-interface formulation of dilute binary alloy solidification with antitrapping to incorporate the highly anisotropic energetic and kinetic properties of the partially faceted ice-water interface. This interface is faceted in the basal plane normal to the ⟨0001⟩ directions and atomically rough in other directions within the basal plane. On the basal plane, the model reproduces a linear or nonlinear relationship between the interface growth rate and the kinetic undercooling that can be linked to experimental measurements. In both cases, spontaneous parity breaking of the solidification front is observed when the preferred growth direction is aligned with the temperature gradient. This phenomenon leads to the formation of partially faceted ice lamellae that drift laterally in one of the ⟨0001⟩ directions. Here, we demonstrate that the drifting velocity of the ice lamellae is controlled by the kinetics on the basal plane and converges as the thickness of the diffuse solid-liquid interface decreases. Furthermore, we examine the effect of the form of the kinetic anisotropy, which is chosen here such that the inverse of the kinetic coefficient varies linearly from a finite value in the ⟨0001⟩ directions to zero in all other directions within the basal plane, consistent with the assumption that the interface grows in local thermodynamic equilibrium in this plane. Our results indicate that the drifting velocity of ice lamellae is not affected by the slope of this linear relation, and the radius and undercooling at the tip of an ice lamella converge at relatively small slope values. Consequently, the phase-field simulations remain quantitative with computationally tractable choices of both the interface thickness and the slope assumed in the form of the kinetic anisotropy.

Materials science↗
DOE OSTI2022

Structural Damage and Recrystallization Response of Garnet Crystals to Intense Electronic Excitation

Here, the susceptibility to irradiation-induced damage, the fundamental mechanism of the relaxation kinetics, and the corresponding recrystallization effect related to the cation radius ratio are comprehensively investigated for Y 3 Al 5 O 12 and Gd 3 Ga 5 O 12 garnet crystals under 645 MeV Xe 35+ irradiation with different fluences of 5 × 10 11 –3 × 10 12 ions cm –2 . Regarding different lattice distortion and swelling levels, the observed microstructure transformations to disordered and amorphous phases, and corresponding hillock dimensions, consistently confirm that Gd 3 Ga 5 O 12 has a higher susceptibility to radiation damage than Y 3 Al 5 O 12 . Combined with iTS model calculations, although Y 3 Al 5 O 12 has higher atomic temperature and energy deposition than Gd 3 Ga 5 O 12 under the same ion velocity and electronic energy loss, the relatively high thermal conductivity and specific heat coefficient of Y 3 Al 5 O 12 crystals enhance the conduction and dissipation of deposition energy, and Gd 3 Ga 5 O 12 , with a higher cation-radius-ratio (r A /r B ), is more easily damaged to amorphous phase due to the less favorable kinetics of ordering and recovery of a melted track region to the crystalline phase. Additionally, the significant bandgap modification in spectral ranges of 5.88–6.75 eV for Y 3 Al 5 O 12 and 4.83–5.41 eV for Gd 3 Ga 5 O 12 , and the enhancement of defect-assisted-related luminescence are achieved, providing a basis to design novel optoelectronic devices in microelectronics fields.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Kinetic Model for the Reduction of Cu II Sites by NO + NH 3 and Reoxidation of NH 3 -Solvated Cu I Sites by O 2 and NO in Cu-SSZ-13

In this work, a kinetic model is developed for the reduction of Cu II sites by NO + NH 3 and the reoxidation of NH 3 -solvated Cu I sites by O 2 and NO in Cu-SSZ-13. Fourier transform infrared (FTIR) spectroscopy and spatially resolved capillary inlet mass spectrometry (SpaciMS) measurements during transient reactor experiments are utilized to identify the rate parameters associated with NO + NH 3 RHC (reduction half-cycle), proposed to occur via two distinct pathways involving adsorbed NH 3 and gas-phase NH 3 . The resulting NO + NH 3 RHC model is validated using spatiotemporal N 2 measurements covering a wide range of temperatures (200–450 °C) and space velocities (53,000–640,000 h –1 ). N 2 O formation is observed and modeled during NO + NH 3 RHC, with quantitative validation under standard selective catalytic reduction (SCR) conditions. Experimentally measured enthalpic and entropic changes associated with O 2 adsorption on NH 3 -solvated Cu I (ZCu(NH 3 ) 2 ) complexes [ Kamasamudram, K. Catal. Today 2010 , 151 (3–4), 212-222], along with activation energies estimated computationally for the intercage diffusion of ZCu(NH 3 ) 2 complexes [ Paolucci, C. Science 2017 , 357 (6 354), 898-903], are incorporated into a mean field kinetic model for the low-temperature oxidation half-cycle (OHC). Significant NH 3 release is observed during the isothermal oxidation of Cu I sites, attributed to desorption of NH 3 ligands from NH 3 -solvated Cu II dimers (Z 2 Cu 2 (NH 3 ) 4 O 2 ). Reduction of these dimeric complexes leads to the consumption of one NO/Cu II , contradicting the expected reduction stoichiometry. Inclusion of a global Arrhenius rate for the NO titration of Z 2 Cu 2 (NH 3 ) 4 O 2 complexes provides accurate representations of standard SCR on reduced and oxidized catalysts, predicting transient NO and NH 3 consumption between 150 and 250 °C as a function of hydrothermal aging. Deactivation of low-temperature standard SCR by NH 3 is observed at high NH 3 pressures, modeled via the formation of superoxo amino (ZCu(NH 3 ) 3 OO*) complexes during NH 3 titration of Z 2 Cu 2 (NH 3 ) 4 O 2 complexes [ Negri, C. J. Am. Chem. Soc. 2020 , 142 (37), 15884-15896]. The redox kinetic model presented here provides a foundational description of active site redox during low-temperature standard SCR, combining the recent kinetic, spectroscopic, and computational findings on the mechanism of standard SCR over Cu-SSZ-13.

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

Trumpler 16-26: a new centrifugal magnetosphere star discovered via SDSS/APOGEE H -band spectroscopy

ABSTRACT We report the discovery of a new example of the rare class of highly magnetized, rapidly rotating, helium enhanced, early B stars that produce anomalously wide hydrogen emission due to a centrifugal magnetosphere (CM). The star is Trumpler 16-26, a B1.5 V member of the Trumpler 16 open cluster. A CM was initially suspected based on hydrogen Brackett series emission observed in SDSS/APOGEE H-band spectra. Similar to the other stars of this type, the emission was highly variable and at all times remarkable due to the extreme velocity separations of the double peaks (up to 1300 km s−1.) Another clue lay in the TESS light curve, which shows two irregular eclipses per cycle when phased with the likely 0.971 8115-d rotation period, similar to the behaviour of the well-known CM host star σ Ori E. To confirm a strong magnetic field and rotation-phase-locked variability, we initiated a follow-up campaign consisting of optical spectropolarimetry and spectroscopy. The associated data revealed a longitudinal magnetic field varying between −3.1 and +1.6 kG with the period found from photometry. The optical spectra confirmed rapid rotation (v sin i = 195 km s−1), surface helium enhancement, and wide, variable hydrogen emission. Tr16-26 is thus confirmed as the 20th known, the fourth most rapidly rotating, and the faintest CM host star yet discovered. With a projected dipole magnetic field strength of Bd > 11 kG, Tr16-26 is also among the most magnetic CM stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Parametric Study of Panel PCM–Air Heat Exchanger Designs

Heat exchangers, devices for the transfer of heat between two or more working fluids, are extensively used in cooling applications and heating applications. Heat exchangers in buildings are typically components of space-conditioning systems, as well as of water-heating applications. Heat exchangers are also sometimes used in applications that require storage and release of energy at specific times. Phase change materials (PCMs) enhance these heat-exchange processes, given their ability to melt and solidify at a fixed range of temperatures, absorbing or releasing significant amounts of latent heat. Five different configurations of PCM–air heat exchangers for thermal control in buildings are analyzed in this work. The heat exchangers were fitted with PCM encapsulated in plastic and composite pouches of various shapes, and packaged in stackable panel layers. Three-dimensional computational fluid dynamics (CFD) modeling of coupled incompressible fluid and conjugate heat transfer were performed on the designs. The phase change process was numerically modelled using the apparent heat capacity method. Steady-state CFD simulations provided quantification of pressure drop as a function of air flow velocity. Transient simulation results describe the thermal evolution of PCM in the pouches, helping to determine the best performing configuration with respect to total thermal charging time.

10 SYNTHETIC FUELS↗
DOE OSTI2021

Seismic anisotropy, dominant slip systems and phase transitions in the lowermost mantle

SUMMARY The presence of seismic anisotropy at the base of the Earth's mantle is well established, but there is no consensus on the deformation mechanisms in lower mantle minerals that could explain it. Strong anisotropy in magnesium post-perovskite (pPv) has been invoked, but different studies disagree on the dominant slip systems at play. Here, we aim to further constrain this by implementing the most recent results from atomistic models and high-pressure deformation experiments, coupled with a realistic composition and a 3-D geodynamic model, to compare the resulting deformation-induced anisotropy with seismic observations of the lowermost mantle. We account for forward and reverse phase transitions from bridgmanite (Pv) to pPv. We find that pPv with either dominant (001) or (010) slip can both explain the seismically observed anisotropy in colder regions where downwellings turn to horizontal flow, but only a model with dominant (001) slip matches seismic observations at the root of hotter large-scale upwellings. Allowing for partial melt does not change these conclusions, while it significantly increases the strength of anisotropy and reduces shear and compressional velocities at the base of upwellings.

58 GEOSCIENCES↗
DOE OSTI2024

Microchannel reactive distillation for the conversion of aqueous ethanol to ethylene

Here we demonstrate the proof-of-concept for microchannel reactive distillation for alcohol-to-jet application: combining ethanol/water separation and ethanol dehydration in one unit operation. Ethanol is first distilled into the vapor phase, converted to ethylene and water, and then the water co-product is condensed to shift the reaction equilibrium. Process intensification is achieved through rapid mass transfer—ethanol stripping from thin wicks using novel microchannel architectures—leading to lower residence time and improved separation efficiency. Energy savings are realized with integration of unit operations. For example, heat of condensing water can offset vaporizing ethanol. Furthermore, the dehydration reaction equilibrium shifts towards completion by immediate removal of the water byproduct upon formation while maintaining aqueous feedstock in the condensed phase. For aqueous ethanol feedstock (40% w ), 71% ethanol conversion with 91% selectivity to ethylene was demonstrated at 220 °C, 600 psig, and 0.28 h -1 wt hour space velocity. 2.7 stages of separation were also demonstrated, under these conditions, using a device length of 8.3 cm. This provides a height equivalent of a theoretical plate (HETP), a measure of separation efficiency, of ~3.3 cm. By comparison, conventional distillation packing provides an HETP of ~30 cm. Thus, 9.1× reduction in HETP was demonstrated over conventional technology, providing a means for significant energy savings and an example of process intensification. Finally, preliminary process economic analysis indicates that by using microchannel reactive distillation technology, the operating and capital costs for the ethanol separation and dehydration portion of an envisioned alcohol-to-jet process could be reduced by at least 35% and 55%, respectively, relative to the incumbent technology, provided future improvements to microchannel reactive distillation design and operability are made.

10 SYNTHETIC FUELS↗
DOE OSTI2023

The impact of the $\text{WHIM}$ on the $\text{IGM}$ thermal state determined from the low- z Lyman $\alpha$ forest

At z ≲ 1, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas (T ~ 10 4 K) in the intergalactic medium (IGM) into warm–hot phase (WHIM) with T > 10 5 K, resulting in a significant deviation from the previously tight power-law IGM temperature–density relationship, T = T 0 (P/$\bar{p}$) γ-1 ⁠. This study explores the impact of the WHIM on measurements of the low-z IGM thermal state, [T 0 , γ], based on the b–NH 1 distribution of the Ly α forest. Exploiting a machine learning-enabled simulation-based inference method trained on Nyx hydrodynamical simulations, we demonstrate that [T 0 , γ] can still be reliably measured from the b–NH 1 distribution at z = 0.1, notwithstanding the substantial WHIM in the IGM. To investigate the effects of different feedback, we apply this inference methodology to mock spectra derived from the IllustrisTNG and Illustris simulations at z = 0.1. The results suggest that the underlying [T 0 , γ] of both simulations can be recovered with biases as low as |Δlog(T 0 /K)| ≲ 0.05 dex, |Δγ| ≲ 0.1, smaller than the precision of a typical measurement. Given the large differences in the volume-weighted WHIM fractions between the three simulations (Illustris 38 percent, IllustrisTNG 10 percent, and Nyx 4 per cent), we conclude that the b–N H1 distribution is not sensitive to the WHIM under realistic conditions. Finally, we investigate the physical properties of the detectable Ly α absorbers, and discover that although their T and Δ distributions remain mostly unaffected by feedback, they are correlated with the photoionization rate used in the simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Seismic Tomographic Modeling of the Crust and Upper Mantle beneath Israel and the Middle East: Improved Resolution through Optimized Model Parameterization

Accurate regional seismic travel-time (RSTT) predictions rely on regional phases (e.g., Pg, Lg, Pn, Sn ) to account for 3D effects in the crust and upper mantle that are not captured by 1D models traditionally used for real-time location. The RSTT prediction model accounts for regional-scale crust and upper mantle structure globally by incorporating regional seismic phases into its travel-time calculations. Previous versions of the RSTT model have used a constant grid cell size of 1°. To improve the tomographic accuracy of recovering velocity structure at regional scales, we, in this study, perform data-driven grid refinement on the RSTT model down to a 0.125° grid (~14 km) in pursuit of two main goals: (1) to test the limits of RSTT capability and accuracy of determined velocity structure through variable grid refinement and (2) to image smaller structures in Israel and the Middle East and illuminate upper mantle dynamics operating in this complex tectonic area. We investigate the effects of model parameterization as grid cell size decreases and the trade-offs between recovered velocity structures. Our final dataset includes 4751 events and 499 stations that recorded 79,344 Pn and 7489 Pg . The variable grid refinement method allows recovery of finer-scale velocity structures and reduces travel-time residuals in areas with the highest data coverage. At smaller grid cell sizes, longer paths need to be upweighted to stabilize the inversion. Results illuminate tectonic features undefined in coarser grid-size models; in particular, we observe mantle perturbations related to the subduction zone around the Cyprian arc and crustal anomalies near the Dead Sea fault and throughout the Anatolian plate.

58 GEOSCIENCES↗
DOE OSTI2020

Chasing Accreted Structures within Gaia DR2 Using Deep Learning

In previous work, we developed a deep neural network classifier that only relies on phase-space information to obtain a catalog of accreted stars based on the second data release of Gaia (DR2). In this paper, we apply two clustering algorithms to identify velocity substructure within this catalog. We focus on the subset of stars with line-of-sight velocity measurements that fall in the range of Galactocentric radii $r\in [6.5,9.5]\,{\rm{kpc}}$ and vertical distances $| z| \lt 3\,{\rm{kpc}}$. Known structures such as Gaia Enceladus and the Helmi stream are identified. The largest previously unknown structure, Nyx, is a vast stream consisting of at least 200 stars in the region of interest. This study displays the power of the machine-learning approach by not only successfully identifying known features but also discovering new kinematic structures that may shed light on the merger history of the Milky Way.

Astronomy & Astrophysics↗
DOE OSTI2022

CTF Improved Drag Model and Flow Regime Transition Criteria

The demand for accurate prediction of two-phase flow behavior in a boiling water reactor (BWR) requires a comprehensive understanding of flow regime, void fraction, heat transfer, and pressure drop. The CTF subchannel code, which is used for the Thermal/Hydraulic (T/H) solution in the Consortium for Advanced Simulation of Light Water Reactors (CASL)-developed Virtual Environment for Reactor Application (VERA) core simulator, is being further developed for BWR applications. In support of this goal, the present work highlights some of the two-phase closure model developments towards improving the CTF void fraction prediction, especially for subcooled boiling. The drift-flux approach has been well-developed for upward dispersed two-phase flows and proven to be accurate in predicting void fraction in bubbly and slug flow regimes. In this work, these kinematic constitutive relations for the drift-flux velocity have been implemented into CTF to describe the interfacial drag of bubbly flow as an alternative to the existing model for better void fraction prediction. The success of these constitutive relations also relies on a good flow regime map that accounts for flow conditions and channel geometry. A more reliable flow regime transition criteria that account for the flow condition has also been implemented in this study for modeling the flow regime transition criteria. The newly implemented models are shown to give improved void fraction predictions in comparison to experimental data.

Hizoum, Belgacem↗
DOE OSTI2020

Phase distribution in Fischer-Tropsch mimicked slurry bubble column via computed tomography

This study aims to investigate the phases distribution and hydrodynamics in an air-Therminol LT-glass beads slurry bubble column reactor in bubbly and churn-turbulent flow regimes by implementing a gamma-ray computed tomography. It is found that gas holdup increases with increasing superficial gas velocity and operating pressure while increasing solids loading has an opposite effect on gas holdup measurements. It is noticed that the pressure effect is more significant at high solids loading and the solids loading is relatively more substantial at lower pressure. Solids loading has the most significant a stronger effect on the solids holdup profile than on the gas holdup profile. Current results were compared with the Sedimentation-Dispersion Model (SDM) predictions. Overall, the findings of this study advance the knowledge of the true behavior of phases distribution in FT mimicked slurry bubble columns which are pertinent in designing such reactors.

42 ENGINEERING↗
DOE OSTI2021

Chemistry and P-V-T equation of state of FeO 2 H x at the base of Earth’s lower mantle and their geophysical implications

The hydrogen-absorbing ability of a mantle mineral in its structure determines the capacity of the water reservoirs hosted by the mineral. Water reservoirs at the base of Earth’s mantle directly influence the fate of water brought down by slab subduction and the seismic heterogeneity such as ultralow-velocity zones (ULVZs) at the core-mantle boundary. Pyrite-FeO 2 H x (0 ≤ x ≤ 1) presents a possibility of such reservoirs in the deep mantle. Ever since the discovery of this mineral phase, however, its chemistry at the lower mantle conditions has been debated. We conducted kinetics experiments of pyrite-FeO 2 H x dehydrogenation at 110 GPa/2100 K, 110 GPa/2300 K, and 120 GPa/2300 K and P-V-T equation of state analysis using in situ synchrotron X-ray diffraction. We found that x approaches 0.80, 0.75, and 0.79, respectively, at the above conditions. The collective P-V-T data yield K 0 = 241(13) GPa, K' = 4.2(4), dK/dT = –0.028(1) GPa/K, α 0 = 4.32(13) × 10 –5 K –1 , and α 1 = 0.31(10) × 10 –8 K –2 for the composition of x = 0.75 ± 0.04. Our first-principles calculations indicate that FeO 2 H 0.75 with a slightly distorted pyrite structure is stable at 100 GPa. These results indicate that this mineral is likely present in the deep mantle with rather a partially dehydrogenated composition than FeO 2 or FeOOH. Furthermore, the results also clarify the difference between the ULVZs originated from pyrite-FeO 2 H x and those from partial melting in terms of shear and compressional wave seismic velocity reduction ratio δlnV S /δlnV P .

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

Carbon Tetrachloride Degradation Results for 200-ZP-1 Operable Unit

Carbon tetrachloride (CT) contamination in the 200-ZP-1 Operable Unit (OU) at the Hanford Site originated from large-volume discharges to the subsurface during plutonium production operations between 1955 and 1973. Contamination migrated through more than 70 meters of unsaturated sediment to reach the underlying unconfined aquifer, where it persists as a large and complex groundwater plume. The 200-ZP-1 OU Record of Decision (ROD) requires that groundwater CT concentrations be reduced to 3.4 µg/L within 125 years. Current groundwater modeling projections estimate that the existing pump-and-treat, even when combined with monitored natural attenuation (specifically hydrolysis), will not achieve this target within the designated timeframe. A fundamental contributor to this shortfall is the extremely slow rate of CT hydrolysis under Hanford aquifer conditions, which has been estimated to have a half-life of 630 years. If faster-acting biotic and abiotic degradation processes are operating within the aquifer, their contribution to CT mass reduction could have a meaningful impact. However, site-specific measurements of these processes and their rates have not previously been performed. This report documents the results of a two-phase laboratory investigation designed to characterize and quantify the capacity of site-specific 200-ZP-1 OU sediments and groundwater to support natural attenuation of CT through biotic and abiotic pathways. In this context, degradation capacity is defined as the intrinsic potential of the subsurface matrix to transform CT under optimized, controlled conditions. System capacity is evaluated in two ways: (1) as rate-limited capacity, which establishes the maximum kinetic velocity of CT transformation and is measured using half-lives and first order rate constants; and (2) as mass limited capacity, which defines the total contaminant mass the batch experimental system can degrade before reactants are exhausted, representing the maximum amount of contaminant the microbial community and reactive mineral phases can transform under the experimental conditions.

abiotic degradation↗
DOE OSTI2021

OSSOS. XXI. Collision Probabilities in the Edgeworth–Kuiper Belt

Here, we present results on the intrinsic collision probabilities, P {sub I}, and range of collision speeds, V {sub I}, as a function of the heliocentric distance, r, in the trans-Neptunian region. The collision speed is one of the parameters that serves as a proxy for a collisional outcome (e.g., disruption and scattering of fragments, or formation of a crater, as both processes are related to the impact energy). We utilize an improved and debiased model of the trans-Neptunian object (TNO) region from the “Outer Solar System Origins Survey” (OSSOS). It provides a well-defined model of TNO orbital distribution, based on multiple opposition observations of more than 1000 bodies. We compute collisional probabilities for the OSSOS models of the main classical, resonant, detached+outer, and scattering TNO populations. The intrinsic collision probabilities and collision speeds are computed using Öpik’s approach, as revised and modified by Wetherill for noncircular and inclined orbits. The calculations are carried out for each of the dynamical TNO groups, allowing for inter-population collisions as well as collisions within each TNO population, resulting in 28 combinations in total. Our results indicate that collisions in the trans-Neptunian region are possible over a wide range in (r, V {sub I}) phase space. Although collisions are calculated to happen within r ∼ 20–200 au and V {sub I} ∼ 0.1 km s{sup −1} to as high as V {sub I} ∼ 9 km s{sup −1}, most of the collisions are likely to happen at low relative velocities V {sub I} < 1 km s{sup −1} and are dominated by the main classical belt.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Delving into the depths of NGC 3783 with XRISM V. Broad-band modelling of ionized outflows

The Seyfert 1 galaxy NGC 3783 hosts a multi-phase warm absorber (WA) that has been extensively studied in the X-ray band. High-resolution spectra from 2000−2001 revealed a complex outflow with multiple ionization and velocity components. Two decades later, new XMM-Newton and XRISM observations allow us to investigate the long-term evolution of these outflows. We performed joint spectral modelling of the XMM-Newton /RGS and XRISM/Resolve time-averaged spectra using the pion photoionization code within SPEX. We derived the ionization parameter, column density, turbulent velocity, and outflow velocity for each absorption component, and investigated their thermal stability and absorption measure distribution (AMD) to characterize the physical and dynamical properties of the WA in NGC 3783 in 2024. We compared these results with the 2000−2001 epoch to assess long-term variability, stability, and possible changes in the absorber population. We identify eight WA components spanning log ξ = 1.08 − 3.38 and outflow velocities of 480−1230 km s −1 . The ranges of column densities and turbulent velocities remain broadly consistent with the WAs from 2000−2001, but the earlier data contained more low-ionization, high-velocity components. The total column density of all the outflows in 2024 is 1.5 times larger than in 2000−2001, which means that it has been replenished by fresh material. The dominant unresolved transition array (UTA) absorber (component B3) has increased its column density by a factor of three while maintaining a similar ionization parameter. The WAs in NGC 3783 have undergone significant structural and dynamical evolution over the past 24 years.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Coupled confined phase behavior and transport of methane in slit nanopores

The proximity of the order of magnitude molecular mean free path and pore sizes in nanopores leads to remarkable interactions between molecules and walls. In such systems, the thermodynamic property and the transport behavior of fluids deviate from those at bulk conditions. Molecular dynamics (MD) simulation may be used to investigate the effects of confinement on fluid physics in nanosize pores. However, translating these subscale observations into a larger scale of interest remains a challenging task. In this work, we propose a modified extension of Peng-Robinson equation of state (PR-EOS) motivated and guided by MD simulation results. The shift of critical properties, i.e., pressure and temperature, are evaluated independently. A temperature dependent parameter is introduced to account for the existence of capillary pressure in the two-phase region. Here, this formulation is capable of capturing the shift in critical properties as well as density phase diagrams under various confinement scenarios. We incorporate the proposed EOS in the lattice Boltzmann method (LBM) to study coupled confined phase behavior and transport of methane in nano-size slit pores. Adsorptive strength is determined such that the resulting density ratio matches that from MD simulation. Results indicate that the use of the proposed EOS leads to smaller fluid velocities compared to PR-EOS. This effect is due to stronger interactions between fluid particles under confinement. Also, results show that transport characteristics are impacted by the pressure of the system: in systems with a relatively low pressure, transport seems to be dominated by Knudsen diffusion; however, at a higher pressure, the contribution from viscous flow increases as the pore widens while the influence of Knudsen and surface diffusion diminishes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗