Search NASA⌕ Search

SEARCH · Search NASA

Results for “bulk viscosity”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Intrinsic bulk viscosity of the one-component plasma

The intrinsic bulk viscosity of the one-component plasma (OCP) is computed and analyzed using equilibrium molecular dynamics simulations and the Green-Kubo formalism. It is found that bulk viscosity exhibits a maximum at Γ ≈ 1, corresponding to the condition that the average kinetic energy of particles equals the potential energy at the average interparticle spacing. The weakly coupled and strongly coupled limits are analyzed and used to construct a model that captures the full range of coupling strengths simulated: Γ ≈ 10 –2 –10 2 . Simulations are also run of the Yukawa one-component plasma (YOCP) in order to understand the impact of electron screening. It is found that electron screening leads to a smaller bulk viscosity due to a reduction in the excess heat capacity of the system. Bulk viscosity is shown to be at least an order of magnitude smaller than shear viscosity in both the OCP and YOCP. The generalized frequency-dependent bulk viscosity coefficient is also analyzed. This is found to exhibit a peak near twice the plasma frequency in strongly coupled conditions, which is associated with the oscillatory decay observed in the bulk viscosity autocorrelation function. As a result, the generalized shear and bulk viscosity coefficients are found to have a similar magnitude for ω ≳ 2⁢ω p at strongly coupled conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Projecting the likely importance of weak-interaction-driven bulk viscosity in neutron star mergers

ABSTRACT In this work, we estimate how much bulk viscosity driven by Urca processes is likely to affect the gravitational wave signal of a neutron star coalescence. In the late inspiral, we show that bulk viscosity affects the binding energy at fourth post-Newtonian order. Even though this effect is enhanced by the square of the gravitational compactness, the coefficient of bulk viscosity is likely too small to lead to observable effects in the waveform during the late inspiral, when only considering the orbital motion itself. In the post-merger, however, the characteristic time-scales and spatial scales are different, potentially leading to the opposite conclusion. We post-process data from a state-of-the-art equal-mass binary neutron star merger simulation to estimate the effects of bulk viscosity (which was not included in the simulation itself). In that scenario, we find that bulk viscosity can reach high values in regions of the merger. We compute several estimates of how much it might directly affect the global dynamics of the considered merger scenario, and find that it could become significant. Even larger effects could arise in different merger scenarios or in simulations that include non-linear effects. This assessment is reinforced by a quantitative comparison with relativistic heavy-ion collisions where such effects have been explored extensively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Bulk viscosity of the rigid rotor one-component plasma

Bulk viscosity of a plasma consisting of strongly coupled diatomic ions is computed using molecular dynamics simulations. The simulations are based on the rigid rotor one-component plasma, which is introduced as a model system that adds two degrees of molecular rotation to the traditional one-component plasma. It is characterized by two parameters: the Coulomb coupling parameter, Γ, and the bond length parameter, Ω. Results show that the long-range nature of the Coulomb potential can lead to long rotational relaxation times, which in turn yield large values for bulk viscosity. Here, the bulk-to-shear viscosity ratio is found to span from small to large values depending on the values of Γ and Ω. Although bulk viscosity is often neglected in plasma modeling, these results motivate that it can be large in molecular plasmas with rotational degrees of freedom.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Bulk viscosity from Urca processes: n p e μ matter in the neutrino-transparent regime

We study the bulk viscosity of moderately hot and dense, neutrino-transparent relativistic npeμ matter arising from weak-interaction direct Urca processes. This work parallels our recent study of the bulk viscosity of npeμ matter with a trapped neutrino component. The nuclear matter is modeled in a relativistic density functional approach with two different parametrizations—DDME2 (which does not allow for the low-temperature direct-Urca process at any density) and NL3 (which allows for low-temperature direct-Urca process above a low-density threshold). Here, we compute the equilibration rates of Urca processes of neutron decay and lepton capture, as well as the rate of the muon decay, and find that the muon decay process is subdominant to the Urca processes at temperatures T ≥ 3 MeV in the case of DDME2 model and T ≥ 1 MeV in the case of NL3 model. Thus, the Urca-process-driven bulk viscosity is computed with the assumption that pure leptonic reactions are frozen. As a result the electronic and muonic Urca channels contribute to the bulk viscosity independently and at certain densities the bulk viscosity of npeμ matter shows a double-peak structure as a function of temperature instead of the standard one-peak (resonant) form. In the final step, we estimate the damping timescales of density oscillations by the bulk viscosity. We find that, e.g., at a typical oscillation frequency f = 1~kHz, the damping of oscillation is most efficient at temperatures 3 ≤ T ≤ 5~MeV and densities n B ≤ 2n 0 where they can affect the evolution of the post-merger object.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Bulk Viscosity of Relativistic npeμ Matter in Neutron-Star Mergers

We discuss the bulk viscosity of hot and dense npeμ matter arising from weak-interaction direct Urca processes. We consider two regimes of interest: (a) the neutrino-transparent regime with T ≤ T tr (T tr ≃ 5 ÷ 10 MeV is the neutrino-trapping temperature); and (b) the neutrino-trapped regime with T ≥ T tr . Nuclear matter is modeled in relativistic density functional approach with density-dependent parametrization DDME2. The maximum of the bulk viscosity is achieved at temperatures T ≃ 5 ÷ 6 MeV in the neutrino-transparent regime, then it drops rapidly at higher temperatures where neutrino-trapping occurs. As an astrophysical application, we estimate the damping timescales of density oscillations by the bulk viscosity in neutron star mergers and find that, e.g., at the oscillation frequency f = 10 kHz, the damping will be very efficient at temperatures 4 ≤ T ≤ 7 MeV where the bulk viscosity might affect the evolution of the post-merger object.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Bulk viscosity from Urca processes: npeμ matter in the neutrino-trapped regime

In this work, we extend our previous study of the bulk viscosity of hot and dense npe matter induced by the Urca process in the neutrino trapped regime to npeμ matter by adding the muonic Urca processes as well as the purely leptonic electroweak processes involving electron-muon transition. Here, the nuclear matter is modeled in a relativistic density functional approach with two different parametrizations which predict neutrino dominated matter (DDME2 model) and anti-neutrino dominated matter (NL3 model) at temperatures for which neutrinos/anti-neutrinos are trapped. In the case of neutrino-dominated matter, the main equilibration mechanism is lepton capture, whereas in the case of antineutrino-dominated matter this is due to neutron decay. We find that the equilibration rates of Urca processes are higher than that of the pure leptonic processes, which implies that the Urca bulk viscosity can be computed with the leptonic reactions assumed to be frozen. We find that the bulk viscosity decreases with temperature as ζ~T –2 at moderate temperatures. At high temperatures this scaling breaks down by sharp drops of the bulk viscosity close to the temperature where the proton fraction is density-independent and the matter becomes scale-invariant. This occurs also when the matter undergoes a transition from the antineutrino-dominated regime to the neutrino-dominated regime where the bulk viscosity attains a local maximum. We also estimate the bulk viscous dissipation timescales and find that these are in the range ≳ 1 s for temperatures above the neutrino trapping temperature. These timescales would be relevant only for long-lived objects formed in binary neutron star mergers and hot proto-neutron stars formed in core-collapse supernovas.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Bulk viscosity of nuclear matter with pions in the neutrino-trapped regime

Recent work has shown that a population of thermal pions could modify the equation of state and transport properties of hot and dense neutron-rich matter and introduce new reaction pathways to change the proton fraction. In this article we study their impact on the bulk viscosity of dense matter, focusing on the neutrino-trapped regime that would be realized in neutron star mergers and supernovae. We find that the presence of a thermal population of pions alters the bulk viscosity by modifying the EoS (via the susceptibilities) and by providing new reaction pathways to achieve beta equilibrium. In neutron star merger conditions, the bulk viscosity in neutrino-trapped npeμ matter (without pions) has its peak at temperatures of at most a couple MeV and is quite small at temperatures of tens of MeV. We find that thermal pions enhance the low-temperature peak of the bulk viscosity by a factor of a few and shift it to slightly lower temperatures. At higher temperatures, where the pion abundance is large but the bulk viscosity is traditionally small, pions can increase the bulk viscosity by an order of magnitude or more, although it is still orders of magnitude smaller than its peak value.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Symmetry-Energy Dependence of the Bulk Viscosity of Nuclear Matter

We clarify how the weak-interaction-driven bulk viscosity ζ and the bulk relaxation time τ Π of neutrino-transparent n p e matter depend on the nuclear symmetry energy. We show that, at saturation density, the equation-of-state dependence of these transport quantities is fully determined by the experimentally constrained nuclear symmetry energy S and its slope L . Variations of L can change the bulk viscosity by orders of magnitude, which can affect both the dissipative and the conservative tidal response of neutron stars. Furthermore, this suggests that both conservative and dissipative effects encoded in the gravitational-wave signatures of binary neutron star inspirals may help constrain nuclear symmetry-energy properties.

Neutron stars & pulsars↗
DOE OSTI2025

Bulk viscosity of two-flavor color superconducting quark matter in neutron star mergers

This work investigates the bulk viscosity of warm, dense, neutrino-transparent, color-superconducting quark matter, where damping of density oscillations in the kHz frequency range arises from weak-interaction-driven direct Urca processes involving quarks. We study the two-flavor red-green paired color-superconducting (2SC) phase, while allowing for the presence of unpaired strange quarks and blue color light quarks of all flavors. Our calculations are based on the SU(3) Nambu–Jona–Lasinio model, extended to include both vector interactions and the ‘t Hooft determinant term. The primary focus is on how variations in the NJL Lagrangian parameters—specifically, the diquark and vector coupling strengths—affect both the static properties of quark matter, such as its equation of state and composition, and its dynamical behavior, including bulk viscosity and associated damping timescales. We find that the bulk viscosity and corresponding damping timescale can change by more than an order of magnitude upon varying the vector coupling by a factor of two at high densities and by a lesser degree at lower densities. This sensitivity primarily arises from the susceptibility of 2SC matter, with a smaller contribution from modifications to the weak interaction rates. In comparison, changes in the diquark coupling have a more limited impact. The damping of density oscillations in 2SC matter is similar quantitatively to nucleonic matter and can be a leading mechanism of dissipation in merging hybrid stars containing color superconducting cores.

Gravitational Waves↗
DOE OSTI2025

Taking control of compressible modes: bulk viscosity and the turbulent dynamo

Many polyatomic astrophysical plasmas are compressible and out of chemical and thermal equilibrium, introducing a bulk viscosity into the plasma via the internal degrees of freedom of the molecular composition, directly impacting the decay of compressible modes, $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$. This is especially important for small-scale, turbulent dynamo processes in the interstellar medium (ISM), which are known to be sensitive to the effects of compression. To control the viscous properties of $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$, we perform trans-sonic, visco-resistive dynamo simulations with additional bulk viscosity $\nu _{\text{bulk}}$, deriving a new $\nu _{\text{bulk}}$ Reynolds number $\text{Re}_{\text{bulk}}$, and viscous Prandtl number $\text{P}\nu \equiv \text{Re}_{\text{bulk}}/ \text{Re}_{\text{shear}}$, where $\text{Re}_{\text{shear}}$ is the shear viscosity Reynolds number. We derive a framework for decomposing $E_{\rm mag}$ growth rates into incompressible and compressible terms via orthogonal tensor decompositions of $\boldsymbol {\nabla }\otimes \mathrm{{\boldsymbol {\mathit {v}}}}$, where $\mathrm{{\boldsymbol {\mathit {v}}}}$ is the fluid velocity. We find that $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ play a dual role, growing and decaying $E_{\rm mag}$, and that field-line stretching is the main driver of growth, even in compressible dynamos. In the absence of $\nu _{\text{bulk}}$ ($\text{P}\nu \rightarrow \infty$), $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ pile up on small-scales, creating a spectral bottleneck, which disappears for $\text{P}\nu \approx 1$. As $\text{P}\nu$ decreases, $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ are dissipated at increasingly larger scales, in turn suppressing incompressible modes through a coupling between high-k modes. We emphasize the importance of further understanding the role of $\nu _{\text{bulk}}$ in compressible astrophysical plasmas, which we estimate could be as strong as the shear viscosity in the cold ISM, and highlight that compressible direct numerical simulations without bulk viscosity have unresolved compressible mode dissipation scales.

MHD↗
DOE OSTI2021

Strangeness-changing rates and hyperonic bulk viscosity in neutron star mergers

In this paper we present a computation of the rates of strangeness-changing processes and the resultant bulk viscosity in matter at the densities and temperatures typical of neutron star mergers. To deal with the high temperature in this environment we go beyond the Fermi surface approximation in our rate calculations and numerically evaluate the full phase space integral. Here, we include processes where quarks move between baryons via meson exchange: these have generally been omitted in previous analyses but provide the dominant contribution to the rates of strangeness-changing processes and the bulk viscosity. The calculation of these rates is an essential step towards any calculation of dissipation mechanisms in hyperonic matter in mergers. As one application, we calculate the dissipation times for density oscillations at the frequencies seen in merger simulations. We find that hyperon bulk viscosity for temperatures in the MeV regime can probably be neglected in this context, but becomes highly relevant for keV-range temperatures.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Emergence of Microphysical Bulk Viscosity in Binary Neutron Star Postmerger Dynamics

Abstract In nuclear matter in isolated neutron stars, the flavor content (e.g., proton fraction) is subject to weak interactions, establishing flavor ( β -)equilibrium. However, there can be deviations from this equilibrium during the merger of two neutron stars. We study the resulting out-of-equilibrium dynamics during the collision by incorporating direct and modified Urca processes (in the neutrino-transparent regime) into general-relativistic hydrodynamics simulations with a simplified neutrino transport scheme. We demonstrate how weak-interaction-driven bulk viscosity in postmerger simulations can emerge and assess the bulk viscous dynamics of the resulting flow. We further place limits on the impact of the postmerger gravitational-wave strain. Our results show that weak-interaction-driven bulk viscosity can potentially lead to a phase shift of the postmerger gravitational-wave spectrum, although the effect is currently on the same level as the numerical errors of our simulation.

Astronomy & Astrophysics↗
DOE OSTI2024

Bulk Viscosity in Dense Nuclear Matter

NEUTRON star mergers are an extreme environment where, within milliseconds, the cold, dense matter in neutron stars is dramatically heated to some of the hottest temperatures encountered in astrophysics. Matter in neutron star mergers is located solidly in the interior of the QCD phase diagram, and therefore, the observation and simulation of these mergers provides an opportunity to learn about the strong interaction at a set of densities and temperatures that compliment those obtained in relativistic heavy-ion collisions. Matter during the merger of two neutron stars is far from hydrostatic equilibrium, allowing for the study of not only the equation of state of dense matter, but of transport properties as well. Perhaps the most likely transport property to play a role in neutron star mergers is the bulk viscosity, which acts to resist changes in density. Furthermore, the bulk viscosity of matter at high densities has long been anticipated to damp oscillations in isolated neutron stars, and the study of its role in neutron star mergers is just beginning.

Harris, Steven P.↗
DOE OSTI2020

Exploring the influence of bulk viscosity of QCD on dilepton tomography

The collective behavior of hadrons and of electromagnetic radiation in heavy-ion collisions has been widely used to study the properties of the quark-gluon plasma (QGP). Indeed this collectivity, as measured by anisotropic flow coefficients, can be used to constrain the transport properties of QGP. The goal of this contribution is to investigate the influence of the specific bulk viscosity (ζ/s) on dilepton production, both at Relativistic Heavy-Ion Collider and Large Hadron Collider energies. We explore the sensitivity of dileptons to dynamical features that bulk viscosity induces on the evolution of a strongly interacting medium, and highlight what makes them a valuable probe in the pursuit to also constrain ζ/s.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Bulk viscosity and cavitation in heavy ion collisions

Relativistic heavy ion collisions generate nuclear-sized droplets of quark-gluon plasma (QGP) that exhibit nearly inviscid hydrodynamic expansion. Smaller collision systems such as p + Au, d + Au, and 3 He + Au at the Relativistic Heavy Ion Collider, as well as p + Pb and high-multiplicity p + p at the Large Hadron Collider may create even smaller droplets of QGP. If so, the standard time evolution paradigm of heavy ion collisions may be extended to these smaller systems. These small systems present a unique opportunity to examine pre-hydrodynamic physics and extract properties of the QGP, such as the bulk viscosity, where the short lifetimes of the small droplets makes them more sensitive to these contributions. Here we focus on the influence of bulk viscosity, its temperature dependence, and cavitation effects on the dynamics in small and large systems using the publicly available hydrodynamic codes SONIC and MUSIC. In conclusion, we also compare pre-hydrodynamic physics in different frameworks including AdS/CFT strong coupling, IP-GLASMA weak coupling, and free streaming or no coupling.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Shear and bulk viscosity for a pure glue theory using an effective matrix model

At nonzero temperatures, the deconfining phase transition can be analyzed using an effective matrix model to characterize the change in holonomy. The model includes gluons and two-dimensional ghost fields in the adjoint representation, or “teens.” As ghosts, the teen fields are responsible for the decrease of the pressure as 𝑇 →𝑇 𝑑 , with 𝑇 𝑑 the transition temperature for deconfinement. Using the solution of this matrix model for a large number of colors, the parameters of the teen fields are adjusted so that the expectation value of the Polyakov loop is close to the values from the lattice. The shear, 𝜂, and bulk, 𝜁, viscosities are computed at nonzero holonomy to leading logarithmic order in weak coupling. In the pure glue theory, the value of the Polyakov loop is relatively large in the deconfined phase, ≈1/2 at 𝑇 𝑑 . Consequently, if 𝑠 is the entropy density, while 𝜂/𝑠 decreases as 𝑇 →𝑇 𝑑 , it is still well above the conformal bound. In contrast, 𝜁/𝑠 is largest at 𝑇 𝑑 , comparable to 𝜂/𝑠, then falls off rapidly with increasing temperature and is negligible by ∼2⁢𝑇 𝑑 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Molecular origins of bulk viscosity in liquid water

The shear viscous response of water is closely associated with changes in network connectivity on the sub ps timescale. The bulk viscous response is shown here to be associated with local density fluctuations and rotational motion around 1–3 ps.

Yahya, Ahmad↗