Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics modeling”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Harnessing complexity: Nonlinear optical phenomena in L-shapes, nanocrescents, and split-ring resonators

Here, we conduct systematic studies of the optical characteristics of plasmonic nanoparticles that exhibit C 2v symmetry. In particular, we analyze three distinct geometric configurations: an L-type shape, a crescent, and a split-ring resonator shaped like the Greek letter π. Optical properties are examined using the finite-difference time-domain method. It is demonstrated that all three shapes exhibit two prominent plasmon modes associated with the two axes of symmetry. This is in addition to a wide range of resonances observed at high frequencies corresponding to quadrupole modes and peaks due to sharp corners. Next, to facilitate nonlinear analysis, we employ a semiclassical hydrodynamic model, where the electron pressure term is explicitly accounted for. This model goes beyond the standard Drude description and enables capturing nonlocal and nonlinear effects. Employing this model enables us to rigorously examine the second-order angular resolved nonlinear optical response of these nanoparticles in each of the three configurations. Two pumping regimes are considered, namely, continuous wave (CW) and pulsed excitations. For CW pumping, we explore the properties of the second harmonic generation (SHG). Polarization and angle-resolved SHG spectra are obtained, revealing strong dependence on the nanoparticle geometry and incident wave polarization. The C 2v symmetry is shown to play a key role in determining the polarization states and selection rules of the SHG signal. For pulsed excitations, we discuss the phenomenon of broadband terahertz (THz) generation induced by the difference-frequency generation . It is shown that the THz emission spectra exhibit unique features attributed to the plasmonic resonances and symmetry of the nanoparticles. The polarization of the generated THz waves is also examined, revealing interesting patterns tied to the nanoparticle geometry. To gain deeper insight, we propose an analytical theory that agrees very well with the numerical experiments. The theory shows that the physical origin of the THz radiation is the mixing of various frequency components of the fundamental pulse by the second-order nonlinear susceptibility. An expression for the far-field THz intensity is derived in terms of the incident pulse parameters and the nonlinear response tensor of the nanoparticle. The results presented in this work offer new insights into the linear and nonlinear optical properties of nanoparticles with C 2v symmetry. The demonstrated strong SHG response and efficient broadband THz generation hold great promise for applications in nonlinear spectroscopy, nanophotonics, and optoelectronics. The proposed theoretical framework also provides a valuable tool for understanding and predicting the nonlinear behavior of other related nanostructures.

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

Experimental and computational study of phase space dynamics in strongly coupled plasmas with steep density gradients

Understanding how plasmas thermalize when density gradients are steep remains a fundamental challenge in plasma physics, with direct implications for fusion experiments and astrophysical phenomena. Standard hydrodynamic models break down in these regimes, and kinetic theories make predictions that have never been directly tested. Here, we present the first detailed phase-space measurements of a strongly coupled plasma as it evolves from sharp density gradients to thermal equilibrium. Using laser-induced fluorescence imaging of an ultracold calcium plasma, we track the complete ion distribution function f(x,v,t). We discover that commonly used kinetic models (Bhatnagar–Gross–Krook and Lenard–Bernstein) overpredict thermalization rates, even while correctly capturing the initial counterstreaming plasma formation. Our measurements reveal that the initial ion acceleration response scales linearly with electron temperature, and that the simulations underpredict the initial ion response. In our geometry we demonstrate the formation of well-controlled counterpropagating plasma beams. This experimental platform enables precision tests of kinetic theories and opens new possibilities for studying plasma stopping power and flow-induced instabilities in strongly coupled systems.

Bergeson, Scott (ORCID:0000000231249226)↗
DOE OSTI2025

Determining the nanostructure of polymer foams using 3D ptycho-tomography for inertial fusion energy applications

Polymer foams play a critical role in contemporary inertial fusion energy (IFE) target designs by enhancing energy yield and optimizing implosion dynamics. However, the lack of high-resolution characterization of the nanostructure of these foams restricts progress in fusion science. In this work, we demonstrate the first high-resolution three-dimensional (3D) reconstruction of a low-density, Si-doped polymer foam fabricated via two-photon polymerization, using ptychographic x-ray computed tomography (PXCT) at an x-ray free electron laser (XFEL). This imaging method reconstructs two-dimensional (2D) attenuation and phase information at multiple sample angles that are combined into a 3D density map used to extract local mass density and determine structural dimensions. We achieve a 2D spatial resolution of 19 ± 3 nm on a high-contrast Ronchi pattern target and 78.7 ± 3 nm for low-contrast polymer foams, marking a significant advancement for XFEL-based ptychography of low-density materials. Furthermore, our experimental results reveal an average foam strut thickness of 1.17 ± 0.4 μm, consistent with fabrication expectations, and a reconstructed average mass density of 0.35 g/cc, aligning closely with the predicted density of 0.29 g/cc. These findings provide important insights for improving foam design and refining radiation hydrodynamics modeling in future IFE experiments. Our study establishes PXCT at an XFEL as a powerful tool for high-resolution characterization of fusion-relevant materials, paving the way for enhanced target performance in IFE research.

Hancock, Levi [Colorado State Univ., Fort Collins,↗
DOE OSTI2025

The origin of the Stokes–Einstein relation in simple dense liquids

Here, we investigate the origin of the universal relation between structural relaxation and diffusion in simple dense liquids, known as the Stokes–Einstein (SE) relation. The fact that this relation, originally derived from a hydrodynamic model of a macroscopic particle in a viscous medium, can describe the microscopic-scale liquid dynamics still eludes understanding. We introduce a new universal measure of structural relaxation in a system of N identical particles based on an explicit decomposition of the configuration space into N! congruent convex polyhedra. This measure makes it possible to quantify the correlation between two distinct particle configurations in terms of their minimal Euclidean distance, optimized with respect to particle permutations. Using this measure alongside a model of independent random walkers under the single-occupancy constraint, we derive a master equation that quantifies the SE relation. It allows us to demonstrate that the universal relation between structural relaxation and diffusion in simple dense liquids is caused by two conditions: (a) the confinement of the dominant density fluctuations to the first coordination shell, manifested by de Gennes narrowing, and (b) Gaussianity of the diffusion process; the former is shown to be violated in low-density fluids, and the latter is known to be violated in supercooled liquids.

Physics - Condensed matter physics↗
DOE OSTI2022

Plasma breakdown in bubbles passing between two pin electrodes

The ignition of plasmas in liquids has applications from medical instrumentation to manipulation of liquid chemistry. Formation of plasmas directly in a liquid often requires prohibitively large voltages to initiate breakdown. Producing plasma streamers in bubbles submerged in a liquid with higher permittivity can significantly lower the voltage needed to initiate a discharge by reducing the electric field required to produce breakdown. The proximity of the bubble to the electrodes and the shape of the bubbles play critical roles in the manner in which the plasma is produced in, and propagates through, the bubble. In this paper, we discuss results from a three-dimensional direct numerical simulation (DNS) used to investigate the shapes of bubbles formed by injection of air into water. Comparisons are made to results from a companion experiment. A two-dimensional plasma hydrodynamics model was then used to capture the plasma streamer propagation in the bubble using a static bubble geometry generated by the DNS The simulations showed two different modes for streamer formation depending on the bubble shape. In an elliptical bubble, a short electron avalanche triggered a surface ionization wave (SIWs) resulting in plasma propagating along the surface of the bubble. In a circular bubble, an electron avalanche first traveled through the middle of the bubble before two SIWs began to propagate from the point closest to the grounded electrode where a volumetric streamer intersected the surface. In an elliptical bubble approaching a powered electrode in a pin-to-pin configuration, we experimentally observed streamer behavior that qualitatively corresponds with computational results. Optical emission captured over the lifetime of the streamer curve along the path of deformed bubbles, suggesting propagation of the streamer along the liquid/gas boundary interface. Plasma generation supported by the local field enhancement of the deformed bubble surface boundaries is a mechanism that is likely responsible for initiating streamer formation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Atmospheric pressure plasmas interacting with wet and dry microchannels: reverse surface ionization waves

Atmospheric pressure plasma jets (APPJs) are increasingly being used to functionalize polymers and dielectric materials for biomedical and biotechnology applications. Once such application is microfluidic labs-on-a-chip consisting of dielectric slabs with microchannel grooves hundreds of microns in width and depth. The periodic channels, an example of a complex surface, present challenges in terms of directly and uniformly exposing the surface to the plasma. In this paper, we discuss results from computational and experimental investigations of negative APPJs sustained in Ar/N 2 mixtures flowing into ambient air and incident onto a series of microchannels. Results from two-dimensional plasma hydrodynamics modeling are compared to experimental measurements of electric field and fast-camera imaging. The propagation of the plasma across dry microchannels largely consists of a sequence of surface ionization waves (SIWs) on the top ridges of the channels and bulk ionization waves (IWs) crossing over the channels. The IWs are directed into electric field enhanced vertices of the next ridge. The charging of these ridges produce reverse IWs responsible for the majority of the ionization. The propagation of the plasma across water filled microchannels evolve into hopping SIWs between the leading edges of the water channels, regions of electric enhancement due to polarization of the water. Positive, reverse IWs follow the pre-ionized path of the initial negative waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Imaging nuclear shape through anisotropic and radial flow in high-energy heavy-ion collisions

Most atomic nuclei exhibit ellipsoidal shapes characterized by quadrupole deformation β 2 and triaxiality γ, and sometimes even a pear-like octupole deformation β 3 . The STAR experiment introduced a new ‘imaging-by-smashing’ technique to image the nuclear global shape by colliding nuclei at ultra-relativistic speeds and analyzing outgoing debris. Features of nuclear shape manifest in collective observables like anisotropic flow vn and radial flow via mean transverse momentum [p T ]. We present new measurements of the variances of v n (n = 2, 3, and 4) and [p T ], and the covariance of $v$$^{2}_{n}$ with [p T ], in collisions of highly deformed 238 U and nearly spherical 197 Au. Ratios of these observables between the two systems effectively suppress common final-state effects, isolating the strong impact of uranium’s deformation. By comparing results with state-of-the-art hydrodynamic model calculations, we extract β 2U and γ U values consistent with those deduced from low-energy nuclear structure measurements. Measurements of v 3 and its correlation with [p T ] also provide the first experimental suggestion of a possible octupole deformation for 238 U. These findings provide significant support for using high-energy collisions to explore nuclear shapes on femtosecond timescales, with implications for both nuclear structure and quark-gluon plasma studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Mixing and dilution controls on marine CO 2 removal using alkalinity enhancement

Marine CO 2 removal (CDR) using enhanced-alkalinity seawater discharge was simulated in the estuarine waters of the Salish Sea, Washington, US. The high-alkalinity seawater would be generated using bipolar membrane electrodialysis technology to remove acid and the alkaline stream returned to the sea. Response of the receiving waters was evaluated using a shoreline resolving hydrodynamic model with biogeochemistry, and carbonate chemistry. Two sites, and two deployment scales, each with enhanced TA of 2997 mmol m -3 and a pH of 9 were simulated. The effects on air-sea CO 2 flux and pH in the near-field as well as over the larger estuary wide domain were assessed. The large-scale deployment (addition of 164 Mmoles TA yr -1 ) in a small embayment (Sequim Bay, 12.5 km 2 ) resulted in removal of 2066 T of CO 2 (45% of total simulated) at rate of 3756 mmol m -2 yr -1 , higher than the 63 mmol m -2 yr -1 required globally to remove 1.0 GT CO 2 yr -1 . It also reduced acidity in the bay, ΔpH ≈ +0.1 pH units, an amount comparable to the historic impacts of anthropogenic acidification in the Salish Sea. The mixing and dilution of added TA with distance from the source results in reduced CDR rates such that comparable amount 2176 T CO 2 yr -1 was removed over >1000 fold larger area of the rest of the model domain. There is the potential for more removal occurring beyond the region modeled. The CDR from reduction of outgassing between October and May accounts for as much as 90% of total CDR simulated. Of the total, only 375 T CO 2 yr -1 (8%) was from the open shelf portion of the model domain. With shallow depths limiting vertical mixing, nearshore estuarine waters may provide a more rapid removal of CO 2 using alkalinity enhancement relative to deeper oceanic sites.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

DES16C3cje: A low-luminosity, long-lived supernova

ABSTRACT We present DES16C3cje, a low-luminosity, long-lived type II supernova (SN II) at redshift 0.0618, detected by the Dark Energy Survey (DES). DES16C3cje is a unique SN. The spectra are characterized by extremely narrow photospheric lines corresponding to very low expansion velocities of ≲1500 km s−1, and the light curve shows an initial peak that fades after 50 d before slowly rebrightening over a further 100 d to reach an absolute brightness of Mr ∼ −15.5 mag. The decline rate of the late-time light curve is then slower than that expected from the powering by radioactive decay of 56Co, but is comparable to that expected from accretion power. Comparing the bolometric light curve with hydrodynamical models, we find that DES16C3cje can be explained by either (i) a low explosion energy (0.11 foe) and relatively large 56Ni production of 0.075 M⊙ from an ∼15 M⊙ red supergiant progenitor typical of other SNe II, or (ii) a relatively compact ∼40 M⊙ star, explosion energy of 1 foe, and 0.08 M⊙ of 56Ni. Both scenarios require additional energy input to explain the late-time light curve, which is consistent with fallback accretion at a rate of ∼0.5 × 10−8 M⊙ s−1.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Second-generation dust in planetary systems: the case of HD 163296

Observations indicate that large, dust-laden protoplanetary discs are common. Some features, like gaps, rings, and spirals, suggest they may host young planets, which can excite the orbits of nearby leftover planetesimals. Energetic collisions among these bodies can lead to the production of second-generation dust. Grains produced by collisions may have a dynamical behaviour different from that of first-generation, primordial dust out of which planetesimals and planets formed. Here, we aim to study these differences for the HD 163296 system and determine whether dynamical signatures in the mixture of the two dust populations can help separate their contributions. We use three-dimensional (3D) hydrodynamic models to describe the gaseous disc with three, Saturn- to Jupiter-mass, embedded planets. Dust grains, of sizes $1\, \mu \mathrm{m}$–$1\, \mathrm{mm}$, are treated as Lagrangian particles with resolved thermodynamics and mass-loss. Initial disc and planet configurations are derived from observation-based work, which indicates low gas viscosity. The 3D approach also allows us to detect the formation of vortices induced by Rossby waves, where dust becomes concentrated and may contribute to planetesimal formation. We find that the main differences in the dynamical behaviour of first- and second-generation dust occur in the vertical distribution. The two populations have similar distributions around the disc mid-plane, although second-generation dust shows longer residence times close to the radial locations of the planets’ gas gaps. Sedimentation rates of $\mu$m-sized grains are comparable to or lower than the production rates by planetesimals’ collisions, making this population potentially observable. These outcomes can be extended to similar systems harbouring giant planets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

SN 2021foa: the bridge between SN IIn and Ibn

ABSTRACT We present the long-term photometric and spectroscopic analysis of a transitioning SN IIn/Ibn from –10.8 d to 150.7 d post V-band maximum. SN 2021foa shows prominent He i lines comparable in strength to the H $\alpha$ line around peak, placing SN 2021foa between the SN IIn and SN Ibn populations. The spectral comparison shows that it resembles the SN IIn population at pre-maximum, becomes intermediate between SNe IIn/Ibn, and at post-maximum matches with SN IIn 1996al. The photometric evolution shows a precursor at –50 d and a light curve shoulder around 17 d. The peak luminosity and colour evolution of SN 2021foa are consistent with most SNe IIn and Ibn in our comparison sample. SN 2021foa shows the unique case of an SN IIn where the narrow P-Cygni in H $\alpha$ becomes prominent at 7.2 d. The H $\alpha$ profile consists of a narrow (500–1200 km s$^{-1}$) component, intermediate width (3000–8000 km s$^{-1}$) and broad component in absorption. Temporal evolution of the H $\alpha$ profile favours a disc-like CSM geometry. Hydrodynamical modelling of the light curve well reproduces a two-component CSM structure with different densities ($\rho \propto$ r$^{-2}$–$\rho \propto$ r$^{-5}$), mass-loss rates (10$^{-3}$–10$^{-1}$ M$_{\odot }$ yr$^{-1}$) assuming a wind velocity of 1000 km s$^{-1}$ and having a CSM mass of 0.18 M$_{\odot }$. The overall evolution indicates that SN 2021foa most likely originated from an LBV star transitioning to a WR star with the mass-loss rate increasing in the period from 5 to 0.5 yr before the explosion or it could be due to a binary interaction.

Gangopadhyay, Anjasha (ORCID:0000000238845637)↗
DOE OSTI2020

Fast electron transport dynamics and energy deposition in magnetized, imploded cylindrical plasma

Inertial confinement fusion approaches involve the creation of high-energy-density states through compression. High gain scenarios may be enabled by the beneficial heating from fast electrons produced with an intense laser and by energy containment with a high-strength magnetic field. Here, we report experimental measurements from a configuration integrating a magnetized, imploded cylindrical plasma and intense laser-driven electrons as well as multi-stage simulations that show fast electrons transport pathways at different times during the implosion and quantify their energy deposition contribution. The experiment consisted of a CH foam cylinder, inside an external coaxial magnetic field of 5 T, that was imploded using 36 OMEGA laser beams. Two-dimensional (2D) hydrodynamic modelling predicts the CH density reaches 9.0 g cm − 3 , the temperature reaches 920 eV and the external B-field is amplified at maximum compression to 580 T. At pre-determined times during the compression, the intense OMEGA EP laser irradiated one end of the cylinder to accelerate relativistic electrons into the dense imploded plasma providing additional heating. The relativistic electron beam generation was simulated using a 2D particle-in-cell (PIC) code. Finally, three-dimensional hybrid-PIC simulations calculated the electron propagation and energy deposition inside the target and revealed the roles the compressed and self-generated B-fields play in transport. During a time window before the maximum compression time, the self-generated B-field on the compression front confines the injected electrons inside the target, increasing the temperature through Joule heating. For a stronger B-field seed of 20 T, the electrons are predicted to be guided into the compressed target and provide additional collisional heating. This article is part of a discussion meeting issue ‘Prospects for high gain inertial fusion energy (part 2)’.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Observation of Long-Range Collective Flow in O + O and Ne + Ne Collisions and Implications for Nuclear Structure Studies

The long-range collective flow of particles produced in oxygen-oxygen ( O + O ) and neon-neon ( Ne + Ne ) collisions is measured with the CMS detector at the CERN LHC. The data samples were collected at a center-of-mass energy per nucleon pair of 5.36 TeV, with integrated luminosities of 7 nb - 1 and 0.8 nb - 1 for O + O and Ne + Ne collisions, respectively. Two- and four-particle azimuthal correlations are measured over nearly five units of pseudorapidity. Significant elliptic ( v 2 ) and triangular ( v 3 ) flow harmonics are observed in both systems. The ratios of v n coefficients between Ne + Ne and O + O collisions suggest sensitivity to the intrinsic nuclear structure of the respective Ne and O ions. Hydrodynamic models with ab initio nuclear structure inputs qualitatively reproduce the collision-overlap dependence of both the v n values and the Ne + Ne to O + O ratios. These measurements provide new constraints when modeling the interplay of nuclear structure and collective dynamics in collisions of O 16 and Ne 20 ions, respectively.

Hayrapetyan, Aram [Yerevan Phys. Inst.]↗
DOE OSTI2025

First measurement of symmetric cumulants of hexagonal flow harmonics in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

Correlations between event-by-event fluctuations of anisotropic flow harmonics are measured in Pb-Pb collisions at a center-of-mass energy per nucleon pair of 5.02 TeV, as recorded by the ALICE detector at the LHC. This study presents correlations up to the hexagonal flow harmonic 𝑣 6 , which was measured for the first time. The magnitudes of these higher-order correlations are found to vary as a function of collision centrality and harmonic order. These measurements are compared to viscous hydrodynamic model calculations with EKRT initial conditions and to the iEBE-VISHNU model with T R ⁢ENTo initial conditions. The observed discrepancies between the data and the model calculations vary depending on the harmonic combinations. Due to the sensitivity of model parameters estimated with Bayesian analyses to these higher-order observables, the results presented in this work provide new and independent constraints on the initial conditions and transport properties in theoretical models used to describe the system created in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Observation of Λ Hyperon Local Polarization in 𝑝-Pb Collisions at $\sqrt{s_{NN}}$ = 8.16 TeV

The polarization of the Λ and $\bar{Λ}$ hyperons along the beam direction has been measured in proton-lead (𝑝-Pb) collisions at a center-of-mass energy per nucleon pair of 8.16 TeV. The data were obtained with the CMS detector at the LHC and correspond to an integrated luminosity of 186.0 ± 6.5 nb −1 . A significant azimuthal dependence of the hyperon polarization, characterized by the second-order Fourier sine coefficient 𝑃 𝑧,s⁢2 , is observed. The 𝑃 𝑧,s⁢2 values decrease as a function of charged particle multiplicity, but increase with transverse momentum. A hydrodynamic model that describes the observed 𝑃 𝑧,s⁢2 values in nucleus-nucleus collisions by introducing vorticity effects does not reproduce either the sign or the magnitude of the 𝑝-Pb results. These observations pose a challenge to the current theoretical implementation of spin polarization in heavy ion collisions and offer new insights into the origin of spin polarization in hadronic collisions at LHC energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Characterizing radial flow fluctuations in relativistic heavy-ion collisions at top RHIC and LHC energies

This study presents a systematic investigation of the transverse-momentum differential radial flow fluctuations observable 𝑣 0 ⁡ (𝑝 𝑇 ) in relativistic heavy-ion collisions at top Relativistic Heavy Ion Collider ($\sqrt{s_{NN}}$ = 200 GeV) and Large Hadron Collider ($\sqrt{s_{NN}}$ = 2.76 and 5.02 TeV) energies. Using a multistage hydrodynamic model, this study assesses the sensitivity of 𝑣 0 ⁡(𝑝 𝑇 ) to a wide range of physical effects, including bulk and shear viscosities, off-equilibrium corrections at particlization, the presence of a hadronic afterburner, and the nucleon size in the initial conditions. By employing complementary rescaling strategies, this study demonstrates how different physical effects leave distinct imprints on the shape of 𝑣 0 ⁡(𝑝 𝑇 ). A combined double-rescaling of 𝑣 0 ⁡(𝑝 𝑇 )/𝑣 0 versus 𝑝 𝑇 /⟨𝑝 𝑇 ⟩ reveals a universality across a wide range of energies and model assumptions in the low-𝑝 𝑇 regime, a robust signature of collective behavior. This allows us to disentangle the universal dynamics of the bulk medium from model-specific features that emerge at higher 𝑝 𝑇 . These results establish 𝑣 0⁡ (𝑝 𝑇 ) as a powerful and complementary observable for constraining quark-gluon plasma transport properties and initial-state granularity, offering a unique probe of the created QCD medium.

Du, Lipei [University of California, Berkeley, CA ↗
DOE OSTI2026

Evidence for the Collective Nature of Radial Flow in Pb+Pb Collisions with the ATLAS Detector

Anisotropic flow and radial flow are two key probes of the expansion dynamics and properties of the quark-gluon plasma (QGP). While anisotropic flow has been extensively studied, radial flow, which governs the system’s radial expansion, has received less attention. Notably, direct experimental evidence for the global and collective nature of radial flow fluctuations has been lacking. This Letter presents the first measurement of transverse momentum (𝑝 T ) dependence of radial flow fluctuations (𝑣 0 ⁡(𝑝 T )) over 0.5 < 𝑝 T < 10 GeV and demonstrates its collective nature using a two-particle correlation method in Pb+Pb collisions at $\sqrt{𝑠_{NN}}$ = 5.02 TeV. The data reveal three key features supporting the collective nature of radial flow: long-range correlation in pseudorapidity, factorization in 𝑝 T , and centrality-independent shape in 𝑝 T . The comparison with a hydrodynamic model demonstrates the sensitivity of 𝑣 0 ⁡(𝑝 T ) to bulk viscosity, a crucial transport property of the QGP. These findings establish a new, powerful tool for probing collective dynamics and properties of the QGP.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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↗