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↗
NASA NTRS1995

Concentration Dependence of Solution Shear Viscosity and Solute Mass Diffusivity in Crystal Growth from Solutions

The physical properties of a supersaturated binary solution such as its density rho, shear viscosity eta, and solute mass diffusivity D are dependent on the solute concentration c: rho = rho(c), eta = eta(c), and D = D(c). The diffusion boundary layer equations related to crystal growth from solution are derived for the case of natural convection with a solution density, a shear viscosity, and a solute diffusivity that are all depen- dent on solute concentration. The solution of these equations has demonstrated the following. (1) At the vicinity of the saturation concentration c(sub s) the solution shear viscosity eta depends on rho as eta(sub s) = eta(rho(sub s))varies as square root of rho(c(sub s)). This theoretically derived result has been verified in experiments with several aqueous solutions of inorganic and organic salts. (2) The maximum solute mass transfer towards the growing crystal surface can be achieved for values of c where the ratio of d ln(D(c)/dc) to d ln(eta(c)/dc) is a maximum.

Izmailov, Alexander F.↗
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↗
NASA NTRS2018

Molecular Dynamics of ULTEM 9085 for 3D Manufacturing: Spectra, Thermodynamic Properties, and Shear Viscosity.

We present results of a molecular dynamic analysis of welding at the polymer-polymer interface. The analysis is performed for polyetherimide/ polycarbonate polymer blends. The work is motivated by the applications to 3D manufacturing in space. In the first part of the report, we discuss bulk and spectral characteristics of the amorphous polymer blends. The vibrational and infra-red spectra obtained using auto-correlation functions calculations in molecular dynamics are compared with the experimental spectra. The mechanical and thermal properties of the samples including heat capacity, bulk modulus, and thermal expansion coefficients are estimated and compared with experimental values. In the second part of the report, we discuss the result of molecular dynamical modeling of shear viscosity in a fully atomistic model of amorphous polymer blends with flat interface. The key result of the research is the demonstration of shear thinning behavior of the shear viscosity as a function of shear rate which is in good agreement with experimental data.

Luchinsky, Dmitry G.↗
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↗
NASA NTRS1996

Relationship Between Solution Shear Viscosity and Density at the Saturation Point

Properties of supersaturated solutions such as the density, viscosity, and solute diffusivity are dependent on the solute concentration. The diffusion-boundary-layer equations are derived and solved for the natural convection case with the viscosity and density dependent on the solute concentration. The solution obtained demonstrates that, at the vicinity of the saturation concentration c(sub s), there is a non-trivial dependence of the solution viscosity eta on its density rho: eta(c(sub s)) = eta(rho(sub s)) varies as rho(sub s)(exp 1/2), where rho(sub s) = rho(c(sub s)). This result has been verified in experiments with aqueous solutions of inorganic and organic salts.

Izmailov, Alexander F.↗
NASA NTRS1999

Comparison of observed rheological properties of hard wheat flour dough with predictions of the Giesekus-Leonov, White-Metzner and Phan-Thien Tanner models

The measured rheological behavior of hard wheat flour dough was predicted using three nonlinear differential viscoelastic models. The Phan-Thien Tanner model gave good zero shear viscosity prediction, but overpredicted the shear viscosity at higher shear rates and the transient and extensional properties. The Giesekus-Leonov model gave similar predictions to the Phan-Thien Tanner model, but the extensional viscosity prediction showed extension thickening. Using high values of the mobility factor, extension thinning behavior was observed but the predictions were not satisfactory. The White-Metzner model gave good predictions of the steady shear viscosity and the first normal stress coefficient but it was unable to predict the uniaxial extensional viscosity as it exhibited asymptotic behavior in the tested extensional rates. It also predicted the transient shear properties with moderate accuracy in the transient phase, but very well at higher times, compared to the Phan-Thien Tanner model and the Giesekus-Leonov model. None of the models predicted all observed data consistently well. Overall the White-Metzner model appeared to make the best predictions of all the observed data.

Non-NASA Center↗
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↗