Search NASA⌕ Search

SEARCH · Search NASA

Results for “viscosity, shear”

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 OSTI2022

Predicted Melt Curve and Liquid Shear Viscosity of RDX up to 30 GPa

Recent grain-scale simulations of HMX and TATB have shown that predictions for hot spot formation in high explosives are particularly sensitive to accurate determinations of the pressure-dependent melt curve and the shear viscosity of the liquid phase. These physics terms are poorly constrained beyond ambient pressure for the explosive RDX. Here we adopt an all-atom modeling approach using molecular dynamics (MD) simulations to predict the melt curve of RDX near to detonation conditions (30 GPa) and determine the shear viscosity of the liquid as a function of temperature and pressure above the melt curve. Phase-coexistence simulations were used to determine the melt curve, which is predicted to vary by almost 1100 K as the pressure increases from 0 GPa to 30 GPa. Equilibrium MD simulations and the Green-Kubo formalism were used to obtain the pressure-temperature-dependent shear viscosity. The shear viscosity of RDX is predicted to be of similar magnitude to the viscosity of TATB at low GPa-range pressures, and to be roughly an order of magnitude lower than the viscosity of HMX. The temperature dependence of the shear viscosity is Arrhenius at a given pressure, and the exponential pre-factor and activation term exhibit a strong, yet complicated, pressure dependence. An empirical pressure-temperature dependent function for RDX shear viscosity is developed that simultaneously captures a wide range of MD predictions while taking an analytic form that extrapolates smoothly beyond the fitted regime. The relative strength of the pressure and temperature dependencies of these two physics terms is found to be of similar magnitude for RDX, HMX, and TATB, which motivates incorporating these results in future RDX grain-scale modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Early-Stage Shear Viscosity far from Equilibrium via Holography

We report shear viscosity is a crucial property of QCD matter which determines the collective behavior of the the quark-gluon plasma (QGP) in ultrarelativistic heavy-ion collisions. Extending the near-equilibrium, high-precision investigations in theory and experiment, we take into account the fact that, in a collision, the QGP is generated far from equilibrium. We use the AdS/CFT correspondence to study a strongly coupled plasma and find a significant impact on the ratio of shear viscosity to entropy density, η/s. In particular, we investigate the initial heating phase and find a decrease reaching down to below 60% followed by an overshoot to 110% of the near-equilibrium value. This finding might be highly relevant for the extraction of transport coefficients from anisotropic flow measurements at RHIC and LHC.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Shear viscosity from perturbative quantum chromodynamics to the hadron resonance gas at finite baryon, strangeness, and electric charge densities

Through model-to-data comparisons from heavy-ion collisions, it has been shown that the quark gluon plasma has an extremely small shear viscosity at vanishing densities. At large baryon densities, significantly less is known about the nature of the shear viscosity from quantum chromodynamics (QCD). Within heavy-ion collisions, there are three conserved charges: baryon number (B), strangeness (S), and electric charge (Q). Here we calculate the shear viscosity in two limits using perturbative QCD (pQCD) and an excluded-volume hadron resonance gas at finite BSQ densities. We then develop a framework that interpolates between these two limits such that shear viscosity is possible to calculate across a wide range of finite BSQ densities. We find that the pQCD and hadron resonance gas calculations have different BSQ density dependencies such that a rather nontrivial shear viscosity appears at finite densities.

Phenomenology↗
DOE OSTI2024

Classical and quantum computing of shear viscosity for ( 2 + 1 ) D SU(2) gauge theory

We perform a nonperturbative calculation of the shear viscosity for ( 2 + 1 )-dimensional SU(2) gauge theory by using the lattice Hamiltonian formulation. The retarded Green’s function of the stress-energy tensor is calculated from real time evolution via exact diagonalization of the lattice Hamiltonian with a local Hilbert space truncation, and the shear viscosity is obtained via the Kubo formula. When taking the continuum limit, we account for the renormalization group flow of the coupling but no additional operator renormalization. We find the ratio of the shear viscosity and the entropy density η s is consistent with a well-known holographic result 1 4 π at several temperatures on a 4 × 4 honeycomb lattice with the local electric representation truncated at j max = 1 2 . We also find the ratio of the spectral function and frequency ρ x y ( ω ) ω exhibits a peak structure when the frequency is small. Both the exact diagonalization method and simple matrix product state classical simulation method beyond j max = 1 2 on bigger lattices require exponentially growing resources. So we develop a quantum computing method to calculate the retarded Green’s function and analyze various systematics of the calculation including j max truncation and finite size effects, Trotter errors and the thermal state preparation efficiency. Our thermal state preparation method still requires resources that grow exponentially with the lattice size, but with a very small prefactor at high temperature. We test our quantum circuit on both the Quantinuum emulator and the IBM simulator for a small lattice and obtain results consistent with the classical computing ones. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX

A multiscale modeling strategy is used to quantify factors governing the temperature rise in hot spots formed by pore collapse from supported and unsupported shock waves in the high explosive HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Two physical aspects are examined in detail, namely the melting temperature and liquid shear viscosity. All-atom molecular dynamics simulations of phase coexistence are used to predict the pressure-dependent melting temperature up to 5 GPa. Equilibrium simulations and the Green–Kubo formalism are used to obtain the temperature- and pressure-dependent liquid shear viscosity. Starting from a simplified continuum-based grain-scale model of HMX, in this study we systematically increase the complexity of treatments for the solid–liquid phase transition and liquid shear viscosity in simulations of pore collapse. Using a realistic pressure-dependent melting temperature completely suppresses melting for supported shocks, which is otherwise predicted when treating it as a constant determined at atmospheric pressure. Alternatively, melt pools form around collapsed pores when the pressure (and melting temperature) are reduced during the release stage of unsupported shocks. Capturing the pressure dependence of the shear viscosity increases the peak temperature of melt pools by hundreds of Kelvin through viscous work. The complicated interplay of the solid-phase plastic work, solid–liquid phase transition, and liquid-phase viscous work identified here motivate taking a systematic approach to building increasingly complex grain-scale models.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Extraction of the specific shear viscosity of quark-gluon plasma from two-particle transverse momentum correlations

The specific shear viscosity, η/s, of the quark-gluon plasma formed in ultrarelativistic heavy-ion collisions at RHIC and LHC is estimated based on the progressive longitudinal broadening of transverse momentum two-particle correlators, G 2 , reported as a function of collision centrality by the STAR and ALICE experiments. Estimates are computed as a function of collision centrality using the Gavin ansatz which relates the G 2 longitudinal broadening to the specific shear viscosity. Freeze out times required for the use of the ansatz are computed using a linear fit of freeze out times reported as a function of the cubic root of the charged particle pseudorapidity density (dN ch /dη) 1/3 . Estimates of η/s based on ALICE data exhibit little to no dependence on collision centrality at LHC energy, while estimates obtained from STAR data hint that η/s might be a function of collision centrality at top RHIC energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Building a testable shear viscosity across the QCD phase diagram

Current experiments at the Relativistic Heavy Ion Collider (RHIC) are probing finite baryon densities where the shear viscosity to enthalpy ratio ηT/w of the Quark Gluon Plasma remains unknown. Here we use the Hadron Resonance Gas (HRG) model with the most up-to-date hadron list to calculate ηT/w at low temperatures and at finite baryon densities ρ B . We then match ηT/w to a QCD-based shear viscosity calculation within the deconfined phase to create a table across {T,μ B } for different cross-over and critical point scenarios at a specified location. We find that these new ηT/w(T,μ B ) values would require initial conditions at significantly larger ρB, compared to ideal hydrodynamic trajectories, in order to reach the same freeze-out point.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Impact of fragment formation on shear viscosity in the nuclear liquid-gas phase transition region

Within the improved quantum molecular dynamics (ImQMD) model we follow the evolution of nuclear matter for planar Couette flow in a periodic box. We focus on the region of liquid-gas phase transition and extract the shear viscosity coefficient from the local stress tensor, directly following the viscosity definition. Here, by switching on and off the mean field and thus inducing the phase transition, we are able to observe the impact of clumping in the phase transition region on the viscosity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Scaling of Shear Rheology of Concentrated Charged Colloidal Suspensions across Glass Transition

Electrostatic interparticle interactions are a key component in controlling and designing rheological characteristics of concentrated charged colloidal suspensions. Herein, we investigate electroviscous effects on shear rheology using highly charged silica particles. By fixing the volume fraction but varying the salinity, the system undergoes a glass transition as evidenced by the evolution of the yield stress and zero-shear viscosity. We show that the steady shear viscosities obey a critical scaling relation that scales the flow curves into a super- and a sub-critical branch with glass transition salinity serving as the bifurcation point; we also demonstrate an isoviscosity scaling that collapses all isoviscosity lines into a single master curve that exhibits no singularity. Based on each scaling relation, in conjunction with common modeling equations, the quantitative relationships between the shear viscosity, stress, and salinity are established. This study demonstrates a new framework to model the steady shear rheology of concentrated charged colloids.

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

The effect of blood viscosity on shear‐induced hemolysis using a magnetically levitated shearing device

Abstract Introduction Blood contacting medical devices, including rotary blood pumps, can cause shear‐induced blood damage that may lead to adverse effects in patients. Due in part to an inadequate understanding of how cell‐scale fluid mechanics impact red blood cell membrane deformation and damage, there is currently not a uniformly accepted engineering model for predicting blood damage caused by complex flow fields within ventricular assist devices (VADs). Methods We empirically investigated hemolysis in a magnetically levitated axial Couette flow device typical of a rotary VAD. The device is able to accurately control the shear rate and exposure time experienced by blood and to minimize the effects of other uncharacterized stresses. Using this device, we explored the effects of both hematocrit and plasma viscosity on shear‐induced hemolysis to characterize blood damage based on the viscosity‐independent shear rate, rather than on shear stress. Results Over a shear rate range of 20 000 – 80 000 1/s, the Index of Hemolysis (IH) was found to be dependent upon and well‐predicted by the shear rate alone. IH was independent of hematocrit, bulk viscosity, or the suspension media viscosity and less correlated to shear stress (MSE = 0.46–0.75) than to shear rate (MSE = 0.06–0.09). Conclusion This study recommends that future investigations of shear‐induced blood damage report findings with respect to the viscosity‐neutral term of shear rate, in addition to the bulk whole blood viscosity measured at an appropriate shear rate relevant to the flow conditions of the device.

Krisher, James A.↗
DOE OSTI2023

Predicted viscosity of liquid HMX up to 40 GPa

Viscous flow serves as a significant heating mechanism during the formation of hot spots, but the shear viscosity which determines this response is poorly characterized for most high explosives. Most recently, a model was proposed for the shear viscosity of liquid HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocane) that was fit to pressures reaching 5 GPa, but this work uncovered uncertainties in the viscosity at 0 GPa and remains untested at higher pressures. We use molecular dynamics (MD) simulations and the Green-Kubo formalism to predict the temperature- and pressure-dependent shear viscosity of HMX over the pressure interval 0 GPa≤P≤40 GPa. Reassessment of the viscosity at 0 GPa rules out several potential explanations for discrepancies between earlier reports; we attribute these differences to details of MD trajectory integration impacting molecular flexibility. The shear viscosity of HMX exhibits an Arrhenius temperature dependence at each pressure considered, with exponential prefactor and activation energy terms that are also strong functions of pressure. An analytic form for the viscosity is developed based on an extension of the well-known Andrade equation that simultaneously captures the temperature and pressure dependencies in the MD data up to 40 GPa. Comparison against a recently developed model for the viscosity of liquid RDX (1,3,5-trinitro-1,3,5-triazinane) shows that both materials exhibit similar functional dependencies with the viscosity of HMX being higher by roughly an order of magnitude at a given temperature-pressure state.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Viscosity of the magnetized strongly coupled one-component plasma

The viscosity tensor of the magnetized one-component plasma, consisting of five independent shear viscosity coefficients, a bulk viscosity coefficient, and a cross coefficient, is computed using equilibrium molecular dynamics simulations and the Green-Kubo relations. A broad range of Coulomb coupling and magnetization strength conditions are studied. Magnetization is found to strongly influence the shear viscosity coefficients when the gyrofrequency exceeds the Coulomb collision frequency. Three regimes are identified as the Coulomb coupling strength and magnetization strength are varied. Here, the Green-Kubo relations are used to separate kinetic and potential energy contributions to each viscosity coefficient, showing how each contribution depends upon the magnetization strength. The shear viscosity coefficient associated with the component of the pressure tensor parallel to the magnetic field, and the two coefficients associated with the component perpendicular to the magnetic field, are all found to merge to a common value at strong Coulomb coupling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Nonlinear Shear Rheology of Entangled Polymer Rings

Steady-state shear viscosity η($\bar{γ}$) of unconcatenated ring polymer melts as a function of the shear rate $\bar{γ}$ is studied by a combination of experiments, simulations, and theory. Experiments using polystyrenes with Z ≈ 5 and Z ≈ 11 entanglements indicate weaker shear thinning for rings compared to linear polymers exhibiting power law scaling of shear viscosity η ~ $\bar{γ}$ –0.56 ± 0.02 , independent of chain length, for Weissenberg numbers up to about 10 2 . Nonequilibrium molecular dynamics simulations using the bead-spring model reveal a similar behavior with η ~ $\bar{γ}$ -0.57 ± 0.08 for 4 ≤ Z ≤ 57. Viscosity decreases with chain length for high $\bar{γ}$. In our experiments, we see the onset of this regime, and in simulations, which we extended to Wi ~ 10 4 , the nonuniversality is fully developed. In addition to a naive scaling theory yielding for the universal regime η ~ $\bar{γ}$ –0.57 , we developed a novel shear slit model explaining many details of observed conformations and dynamics as well as the chain length-dependent behavior of viscosity at large $\bar{γ}$. The signature feature of the model is the presence of two distinct length scales: the size of tension blobs and much larger thickness of a shear slit in which rings are self-consistently confined in the velocity gradient direction and which is dictated by the size of a chain section with relaxation time 1/$\bar{γ}$. These two length scales control the two normal stress differences. In this model, the chain length-dependent onset of nonuniversal behavior is set by tension blobs becoming as small as about one Kuhn segment. This model explains the approximate applicability of the Cox–Merz rule for ring polymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗