SEARCH · Engineering Papers

Results for “Statistical hydrodynamics”

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.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2026

ELG×LRG Distribution through Dark Matter Halo Dynamics

We investigate the clustering and halo occupation distribution (HOD) of DESI Y1 emission-line (ELGs) and luminous red (LRGs) galaxies at 0.8 < z < 1.1, including their cross-correlation (ELG×LRG), using the A BACUS S UMMIT suite and a new Halo Occupation Model (H OME ) for galaxy multitracers. This integrates intrahalo dynamics, halo exclusion, and quenching, bridging insights from hydrodynamical, HOD, abundance-matching, and semianalytic studies. Leveraging full phase-space information from the Uchuu N-body simulation, and sampling satellites from dark-matter particle positions via physically motivated prescriptions, Home reproduces the anisotropic clustering down to s = 200 h −1 kpc with unprecedented accuracy. Model parameters are inferred solely from two-point statistics using a two-level Bayesian framework, yielding high-fidelity ELG, LRG, and cross-reference catalogs. We find that satellite ELGs behave as incoherent flows within their parent halos, dominating the clustering below 4 h −1 Mpc. The HOD from the best-fit Home has the following properties: (i) 90.50% (85.91%) of ELGs (LRGs) are central galaxies without satellites, residing in halos of M vir ∼ 6.6 × 10 11 (1.2 × 10 13 ) h −1 M ⊙ ; (ii) the ELG×LRG cross-correlation is governed by central-central pairs and shaped by halo exclusion on 2–5 h −1 Mpc scales; (iii) 9.50% (14.09%) of ELGs (LRGs) are satellites, of which 1.09% (3.52%) inhabit halos with a central galaxy of the same species in a maximally conformal configuration, 7.02% (0.005%) orbit complementary hosts in a minimally conformal state, and 0.58% (10.57%) are orphans. The high sensitivity of Home precisely captures the dynamics of satellites in different host environments, opening a promising avenue for understanding systematics and the dynamical nature of dark matter, potentially distinguishing gravity models.

Favole, Ginevra [Universidad de La Laguna (Spain);↗

DOE OSTI2025

Emergent hydrodynamic mode on SU(2) plaquette chains and quantum simulation

We search for emergent hydrodynamic modes in real-time Hamiltonian dynamics of 2+1-dimensional SU(2) lattice gauge theory on a quasi-one-dimensional plaquette chain, by numerically computing symmetric correlation functions of energy densities on lattice sizes of about 20 with the local Hilbert space truncated at 𝑗 max = $\frac{1}{2}$. Because of the Umklapp processes, we only find a mode for energy diffusion. The symmetric correlator exhibits transport peak near zero frequency with a width approximately proportional to momentum squared at small momentum, when the system is fully quantum ergodic, as indicated by the eigenenergy level statistics. This transport peak leads to a power-law 𝑡 −$\frac{1}{2}$ decay of the symmetric correlator at late time, also known as the long-time tail, as well as diffusionlike spreading in position space. We also introduce a quantum algorithm for computing the symmetric correlator on a quantum computer and find it gives results consistent with exact diagonalization when tested on the IBM emulator. Finally we discuss the future prospect of searching for the sound modes.

Hamiltonian systems↗

DOE OSTI2025

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

Forecasting constraints on the high-z IGM thermal state from the Lyman-α forest flux autocorrelation function

ABSTRACT The autocorrelation function of the Lyman-$\alpha$ (Ly $\alpha$) forest flux from high-z quasars probes the small-scale structure of the intergalactic medium (IGM). The thermal state of the IGM, determined by the physics of reionization, sets the small-scale power observed in the Ly $\alpha$ forest. To explore the sensitivity of the autocorrelation function to the IGM’s thermal state, we compute the autocorrelation function from a cosmological hydrodynamical simulation with an instantaneous reionization model and 135 post-processed thermal states. Using mock data sets of 20 quasars, we forecast constraints on $T_0$ and $\gamma$, which characterize the post-processed IGM thermal state, at $5.4 \le z \le 6$. While this model simplifies the IGM’s thermal state, it serves as a key first step in assessing future observational prospects. We also perform an inference test on mocks and re-weight out posterior distributions to guarantee that they exhibit statistically correct behaviour. At $z = 5.4$, we find that an idealized data set constrains $T_0$ to 59 per cent and $\gamma$ to 16 per cent at the 1$\sigma$ equivalent confidence level. To explore more realistic, non-instantaneous reionization scenarios, we analyse four models combining temperature and ultraviolet background (UVB) fluctuations at $z = 5.8$. We find that mock data generated from a model with both temperature and UVB fluctuations can rule out a model with only temperature fluctuations at the $> 1\sigma$ level 73.9 per cent of the time.

Wolfson, Molly↗

DOE OSTI2025

Exploring scenarios for enhanced fuel compression and performance on the National Ignition Facility with machine-learning-aided design techniques

Recent fusion experiments on the National Ignition Facility (NIF) have achieved ignition, producing multi-MJ fusion yields for input laser energies of roughly 2 MJ [Abu-Shawareb et al., Phys. Rev. Lett. 132, 065102 (2024)]. Building on the success of the target designs that have achieved ignition, we explore new implosion scenarios predicted to generate significantly more compression of the dense DT ice layer and correspondingly higher yields while preserving many of the key physics characteristics of present-day ignition designs. Our main result is a novel 3-shock implosion scheme that effectively minimizes the shock-induced entropy in the dense, accelerating DT shell and maximizes the resulting fuel compression subject to a fixed leading shock strength consistent with present-day ignition experiments, which is necessary to melt the crystalline high-density carbon ablator. Compared to the first NIF experiment to fulfill Lawson's ignition criterion, shot N210808 [Abu-Shawareb et al., Phys. Rev. Lett. 129, 075001 (2022)], our design exhibits a 40% increase in simulated peak areal density (ρR) and a 5× increase in 1D fusion yield using a 4% lighter ablator and identical DT payloads. We also present a complete integrated 2D hohlraum design and laser pulse specifications capable of generating the desired 3-shock drive and maintaining control of the low-mode capsule implosion symmetry, where the increase in simulated 2D yield relative to N210808 is > 10×. This new implosion regime was discovered with help from a machine-learning-enabled capsule design optimization framework. We outline the workflow this automated tool uses to identify improved design candidates by running several rounds of capsule simulations, constructing a surrogate model mapping input variations to key physics output quantities, and querying the resulting statistical model to propose adjustments to the x-ray drive and capsule to reach a set of physics objectives prescribed by the designer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Uncertainty Quantification Enabled by Automatic Differentiation for Hydrodynamic Simulation of Shock‐to‐Detonation Transition in High Explosives

Quantifying the effects of uncertainty in a reactive burn model on the run-to-detonation time in high explosives (HEs) provides a robust methodology for assessing the probability of an HE failing the IHE qualification standard. Moreover, uncertainty quantification helps evaluate whether the model calibration accurately represents data outside the calibration set. This study uses a specialized hydrodynamic simulation code for modeling detonation to determine the run-to-detonation time of the HE PBX 9502 for various impact velocities. To quickly approximate uncertainties in the model, a surrogate was constructed using a Taylor series expansion centered at the mean of the input parameters. To obtain the sensitivities required for constructing the Taylor series, HYP-percomplex Automatic Differentiation (HYPAD) was implemented. HYPAD is a methodology for infusing existing codes with automatic differentiation capabilities by augmenting variables with one or more imaginary units to compute step-size independent partial derivatives. These derivatives are accurate to machine precision with respect to the implemented numerical algorithm, meaning their accuracy reflects that of the underlying method (e.g., integration or discretization schemes). Using reduced order modeling techniques, the mean and standard deviation of the run-to-detonation time of a shock within PBX 9502 were computed for a number of initial impact velocities. A weighted least squares regression was then performed to obtain a best fit curve and prediction interval for the computed statistics. Historical data points from explosively driven wedge tests were utilized to validate the prediction interval, ensuring its reliability in predicting future outcomes. With this prediction interval and a known safety constraint curve, the most probable point of failure and the probability of failure for the HE PBX 9502 were determined.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Dwarf Galaxies at Cosmic Noon: New JWST Constraints on Satellite Models and Subhalo Tidal Evolution

The advent of JWST has revolutionized the study of faint satellite galaxies at z ≳ 1, enabling statistical constraints on galaxy evolution and the galaxy–halo connection in a previously unexplored mass and redshift regime. We compare satellite abundances at 1 < z < 3.5 from recent JWST observations with predictions from cosmological dark-matter-only zoom-in simulations. We identify and quantify several sources of biases that can impact theoretical satellite counts, finding that assumptions about subhalo tidal evolution introduce the largest uncertainty in predictions for the satellite mass function. Using a flexible galaxy disruption model, we explore a range of disruption scenarios, spanning hydrodynamically motivated and idealized prescriptions, to bracket plausible physical outcomes. We show that varying galaxy durability can change the predicted satellite mass functions by a factor of ∼3.5. The JWST data and our fiducial model are consistent within 1σ–2σ across the full redshift (1 < z < 3.5) and stellar mass (M⋆ > 10 7 M ⊙ ) range probed. We find evidence that subhalos are at least as long-lived as predicted by hydrodynamic simulations. Our framework will enable robust constraints on the tidal evolution of subhalos with future observations. This work presents the first direct comparison between cosmological models and observations of the high-redshift satellite population in this low-mass regime. These results showcase JWST’s emerging power to test structure formation in the first half of the Universe in a new domain and to constrain the physical processes driving the evolution of low-mass galaxies across cosmic time.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Measurement of (anti)alpha production in central Pb–Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV

In this letter, measurements of (anti)alpha production in central (0–10%) Pb–Pb collisions at a center-of-mass energy per nucleon–nucleon pair of $\sqrt{s_{NN}} = 5.02$ TeV are presented, including the first measurement of an antialpha transverse-momentum spectrum. Owing to its large mass, the production of (anti)alpha is expected to be sensitive to different particle production models. The production yields and transverse-momentum spectra of nuclei are of particular interest because they provide a stringent test of these models. The averaged antialpha and alpha spectrum is compared to the spectra of lighter particles, by including it into a common blast-wave fit capturing the hydrodynamic-like flow of all particles. This fit is indicating that the (anti)alpha also participates in the collective expansion of the medium created in the collision. A blast-wave fit including only protons, (anti)alpha, and other light nuclei results in a similar flow velocity as the fit that includes all particles. A similar flow velocity, but a significantly larger kinetic freeze-out temperature is obtained when only protons and light nuclei are included in the fit. The coalescence parameter B 4 is well described by calculations from a statistical hadronization model but significantly underestimated by calculations assuming nucleus formation via coalescence of nucleons. Similarly, the (anti)alpha-to-proton ratio is well described by the statistical hadronization model. On the other hand, coalescence calculations including approaches with different implementations of the (anti)alpha substructure tend to underestimate the data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Mapping the Perseus galaxy cluster with XRISM: Gas kinematic features and their implications for turbulence

We present extended gas kinematic maps of the Perseus cluster based on a combination of five new XRISM/Resolve pointings observed in 2025 with four performance verification datasets from 2024, totaling a net exposure of 745 ks. To date, Perseus remains the only cluster that has been extensively mapped out to ≃0.7 r 2500 by XRISM/Resolve, while simultaneously offering sufficient spatial resolution to resolve gaseous substructures driven by mergers and active galactic nucleus (AGN) feedback. Our observations cover multiple radial directions and a broad range of dynamical scales, enabling us to characterize the kinematic properties of the intracluster medium up to a scale of ∼500 kpc. In the measurements, we detected high-velocity dispersions (≃300km s −1 ) in the eastern region of the cluster that are spatially coincident with the extended X-ray surface brightness excess and correspond to a nonthermal pressure fraction of ≃7 − 13%. The velocity field outside the AGN-dominant region can be effectively described by a single, large-scale kinematic driver based on the velocity structure function, which statistically favors an energy injection scale of at least a few hundred kpc. The estimated turbulent dissipation energy is comparable to the gravitational potential energy released by a recent merger, implying a significant role of turbulent cascade in the merger energy conversion. In the bulk velocity field, we observed a dipole-like pattern along the east-west direction with an amplitude of ≃ ± 200 − 300 km s −1 , indicating rotational motions induced by the recent merger event. This feature constrains the viewing direction to ≃30° −50° relative to the normal of the merger plane. Our hydrodynamic simulations suggest that Perseus has experienced at least two energetic mergers since redshift z ∼ 1, the most recent of which is associated with the radio galaxy IC310, in agreement with recent SRG/eROSITA findings. This study showcases exciting scientific opportunities for future missions with high-resolution spectroscopic capabilities (e.g., HUBS, LEM, and NewAthena).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

r -process Nucleosynthesis from Hyperaccreting Neutron Stars in Common Envelopes

We investigate nuclear reactions and feedback in hyperaccreting neutron star environments, considering accretion rates in the range 0.3–3 × 10 4 M ⊙ yr −1 , typical of short-period compact-object binaries in common envelopes. Our models account for weak reactions, neutrino energy loss, nuclear energy release, pair production, degenerate equations of state, and general relativistic hydrodynamics. Depending on the accretion rates, these systems can develop both proton- and neutron-rich atmospheres with strong convective instabilities linking the neutrino sphere to the outgoing accretion shock inside the radiation trapping zone. Convection drives nucleons through multiple heating and cooling cycles, with photodisintegration dominating during the heating phase and heavy element synthesis during the cooling phase, ejecting material with abundances that depend on the accretion rate and depth of the final decompression trajectory. The turbulent nature of convective currents is conducive to creating a wide range of nuclear products through a variety of effects, including nuclear statistical equilibrium freeze-out and the r-, p-, and γ-processes. We also observe a novel multistep process in reheated trajectories, consisting of proton-capture and photodissociation reactions operating on r-process seeds, producing overall neutron-deficient isotopes. A significant amount of infalling gas experiences high entropy and short (millisecond) freeze-out timescales capable of making r-process elements with high overabundances through a disequilibrium effect between neutrons and α-particles that does not require an excess of neutrons.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Temperature measurement of Quark-Gluon plasma at different stages

In a Quark-Gluon Plasma (QGP), the fundamental building blocks of matter, quarks and gluons, are under extreme conditions of temperature and density. A QGP could exist in the early stages of the Universe, and in various objects and events in the cosmos. The thermodynamic and hydrodynamic properties of the QGP are described by Quantum Chromodynamics (QCD) and can be studied in heavy-ion collisions. Despite being a key thermodynamic parameter, the QGP temperature is still poorly known. Thermal lepton pairs (e + e − and μ + μ − ) are ideal penetrating probes of the true temperature of the emitting source, since their invariant-mass spectra suffer neither from strong final-state interactions nor from blue-shift effects due to rapid expansion. Here we measure the QGP temperature using thermal e+e− production at the Relativistic Heavy Ion Collider (RHIC). The average temperature from the low-mass region (in-medium ρ 0 vector-meson dominant) is (2.01 ± 0.23) × 10 12 K, consistent with the chemical freeze-out temperature from statistical models and the phase transition temperature from Lattice QCD. The average temperature from the intermediate mass region (above the ρ0 mass, QGP dominant) is significantly higher at (3.25 ± 0.60) × 10 12 K. This work provides essential experimental thermodynamic measurements to map out the QCD phase diagram and understand the properties of matter under extreme conditions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Subexponential Decay of Local Correlations from Diffusion-Limited Dephasing

Chaotic quantum systems at finite entropy density are expected to act as their own heat baths, rapidly dephasing local quantum superpositions. Here, we argue that in fact this dephasing is generically subexponential in one-dimensional systems with conservation laws: all local correlation functions decay as exp⁡[−𝒪⁡(𝑡 𝛼 )] with 0 ≤ 𝛼 ≤ 2/3, even when the operators are orthogonal to all hydrodynamic modes. The mechanism is diffusion-limited dephasing, in which rare low-entropy regions (“voids”) protect quantum coherences. This intrinsically quantum effect lies beyond standard hydrodynamics and disappears under extrinsic dephasing. In random charge-conserving circuits we find 𝛼 = 1/2, while in generic translation-invariant Floquet systems we bound 𝛼 ≤ 2/3. Our arguments are general, subject principally to the assumption that thermal fluctuations can create regions of zero entropy density. In systems with energy conservation, this assumption is automatically satisfied because of the third law of thermodynamics.

information scrambling↗

DOE OSTI2026

Strangeness enhancement at its extremes: multiple (multi-)strange hadron production in pp collisions at $\sqrt{s}=5.02$ TeV

The probability to observe a specific number of strange and multi-strange hadrons (nS), denoted as P(nS), is measured by ALICE at midrapidity (|y| < 0.5) in $$\sqrt{s}=5.02$$ TeV proton-proton (pp) collisions, dividing events into several multiplicity-density classes. Exploiting, for the first time, a technique based on counting the number of strange-particle candidates event-by-event, this measurement allows one to extend the study of strangeness production beyond the mean of the distribution. This constitutes a new test bench for production mechanisms, probing events with a large imbalance between strange and non-strange content. The analysis of a large-statistics data sample makes it possible to extract P(nS) up to a maximum nS of 7 for $${\text{K}}_{\text{S}}^{0}$$, 5 for Λ and $$\overline{\Lambda }$$, 4 for Ξ− and $${\overline{\Xi } }^{+}$$, and 2 for Ω− and $${\overline{\Omega } }^{+}$$. From this, the probability of producing strange hadron multiplets per event is calculated, thereby enabling the extension of the study of strangeness enhancement to extreme situations where several strange quarks hadronize in a single event at midrapidity. Moreover, comparing hadron combinations with different u and d quark compositions and equal overall s quark content, the contribution to the enhancement pattern coming from non-strangeness related mechanisms is isolated. The results are compared with state-of-the-art phenomenological models implemented in commonly used Monte Carlo event generators, including PYTHIA 8 Monash 2013, PYTHIA 8 with QCD-based Color Reconnection and Rope Hadronization (QCD-CR + Ropes), and EPOS LHC, which incorporates both partonic interactions and hydrodynamic evolution. These comparisons show that the new approach dramatically enhances the sensitivity to the different underlying physics mechanisms modeled by each generator.

Abualrob, I J↗

DOE OSTI2024

The 3D Lyman- α forest power spectrum from eBOSS DR16

We measure the three-dimensional power spectrum (P3D) of the transmitted flux in the Lyman-α (Ly α) forest using the complete extended Baryon Oscillation Spectroscopic Survey data release 16 (eBOSS DR16). This sample consists of ~205 000 quasar spectra in the redshift range 2 ≤ z ≤ 4 at an effective redshift z = 2.334. We propose a pair-count spectral estimator in configuration space, weighting each pair by exp( i k ∙ r), for wave vector k and pixel pair separation r, effectively measuring the anisotropic power spectrum without the need for fast Fourier transforms. This accounts for the window matrix in a tractable way, avoiding artefacts found in Fourier-transform based power spectrum estimators due to the sparse sampling transverse to the line of sight of Ly α skewers. We extensively test our pipeline on two sets of mocks: (i) idealized Gaussian random fields with a sparse sampling of Ly α skewers, and (ii) log-normal LyaCoLoRe mocks including realistic noise levels, the eBOSS survey geometry and contaminants. On eBOSS DR16 data, the Kaiser formula with a non-linear correction term obtained from hydrodynamic simulations yields a good fit to the power spectrum data in the range $(0.02 ≤ k ≤ 0.35)$ h Mpc -1 at the 1–2σ level with a covariance matrix derived from LyaCoLoRe mocks. We demonstrate a promising new approach for full-shape cosmological analyses of Ly α forest data from cosmological surveys such as eBOSS, the currently observing Dark Energy Spectroscopic Instrument and future surveys such as the Prime Focus Spectrograph, WEAVE-QSO, and 4MOST.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Merger Response of Halo Anisotropy Properties

Anisotropy properties—halo spin, shape, position offset, velocity offset, and orientation—are an important family of dark matter halo properties that indicate the level of directional variation of the internal structures of halos. These properties reflect the dynamical state of halos, which in turn depends on the mass assembly history. In this work, we study the evolution of anisotropy properties in response to merger activity using the IllustrisTNG simulations. We find that the response trajectories of the anisotropy properties significantly deviate from secular evolution. These trajectories have the same qualitative features and timescales across a wide range of merger and host properties. We propose explanations for the behavior of these properties and connect their evolution to the relevant stages of merger dynamics. We measure the relevant dynamical timescales. We also explore the dependence of the strength of the response on time of merger, merger ratio, and mass of the main halo. These results provide insight into the physics of halo mergers and their effects on the statistical behavior of halo properties. This study paves the way toward a physical understanding of scaling relations, particularly to how systematics in their scatter are connected to the mass assembly histories of halos.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Galaxy cluster profiles: a Gaussian mixture model approach to halo miscentering

Measurements of the galaxy density and weak-lensing profiles of galaxy clusters typically rely on an assumed cluster center, which is taken to be the brightest cluster galaxy or other proxies for the true halo center defined as the minimum in the potential well. Departure of the assumed cluster center from the true halo center bias the resultant profile measurements, an effect known as miscentering bias. Currently, miscentering is typically modeled in stacked profiles of clusters with a two parameter model. We use an alternate approach in which the profiles of individual clusters are used with the corresponding likelihood computed using a Gaussian mixture model. We test the approach using halos and the corresponding subhalo profiles from the IllustrisTNG hydrodynamic simulations. We obtain significantly improved estimates of the miscentering parameters for both 3D and projected 2D profiles relevant for imaging surveys. We discuss applications to upcoming cosmological surveys. Our Python package for the Gaussian mixture model is publicly available at https://github.com/KyleMiller1/Halo-Miscentering-Mixture-Model.

Bayesian reasoning↗

DOE OSTI2025

Gas and dust dynamics in γ Cephei-type disks

Context. Giant planets are observed orbiting the primary stars of close binary systems. Such planets may have formed in compact circumprimary disks, which once surrounded these stars, under conditions much different than those encountered around single stars. Aims. In order to quantify the effects of the strong gravitational perturbations exerted on circumprimary disk material, the three-dimensional (3D) dynamics of gas and dust in orbit around the primary star of a compact and eccentric binary system was modeled by applying the stellar and orbital parameters of γ Cephei, a well-known system that can be representative of a class of close binaries. Methods. Circumprimary gas was approximated as an Eulerian viscous and compressible fluid and modeled by means of 3D hydrodynamical simulations, assuming locally isothermal conditions in the medium around the primary star. Dust grains were modeled as Lagrangean particles, subjected to gravity and aerodynamic drag forces. Models that include a giant planet were also considered. Results. Models indicate that spiral density waves excited around pericenter passage propagate toward the inner boundary of the disk, through at least a few pressure scale-heights from the mid-plane, inducing radial and vertical mixing in the gas. However, perturbations imparted to gas, both in terms of eccentricity and precession, are far weaker than previously estimated by two-dimensional (2D) simulations. Models predict small eccentricities, ≲0.03, and slow retrograde precession. The addition of a giant planet does not change the low eccentricity state of the disk. The parameters applied to the disk would lead to the formation of a massive planet, many times the mass of Jupiter, in agreement with some observations. Micron to mm-size dust grains are well coupled to the gas, resulting in similar dynamics and statistically similar distributions of orbital elements. The planet only affects the dust distributions locally. In agreement with outcomes of recent 2D models, the lifetime of an isolated circumprimary disk would be brief, ~10 5 years, because of its compact nature, requiring a long-term external supply of mass to allow for the in situ formation of a giant planet.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Monitored Fluctuating Hydrodynamics

We introduce a hydrodynamic framework for describing monitored classical stochastic processes. We study the conditional ensembles for these monitored processes—i.e., we compute spacetime correlation functions conditioned on a fixed, typical measurement record. In the presence of global symmetries we show that these conditional ensembles can undergo measurement-induced “sharpening” phase transitions as a function of the monitoring rate; moreover, even weak monitoring can give rise to novel critical phases, derived entirely from a classical perspective. We give a simple hydrodynamic derivation of the known “charge-fuzzy phase” for weakly monitored diffusive many-body quantum systems. We show that although the unmonitored symmetric and asymmetric exclusion processes are in different universality classes of transport, the fluctuations in their conditional ensembles flow to the same fixed point with emergent relativistic invariance under monitoring. On the other hand, weakly monitored systems with non-Abelian symmetries enter a novel strongly coupled fixed point with nontrivial dynamical exponent, which we characterize. Our formalism naturally accounts for monitoring general observables, such as currents or density gradients, and allows for a direct calculation of information-theoretic diagnostics of sharpening transitions, including the Shannon entropy of the measurement record.

classical statistical mechanics↗