Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Effect of Viscosity on the Crystallization of Undercooled Liquids

There have been numerous studies of glasses indicating that low-gravity processing enhances glass formation. NASA PI s are investigating the effect of low-g processing on the nucleation and crystal growth rates. Dr. Ethridge is investigating a potential mechanism for glass crystallization involving shear thinning of liquids in 1-g. For shear thinning liquids, low-g (low convection) processing will enhance glass formation. The study of the viscosity of glass forming substances at low shear rates is important to understand these new crystallization mechanisms. The temperature dependence of the viscosity of undercooled liquids is also very important for NASA s containerless processing studies. In general, the viscosity of undercooled liquids is not known, yet knowledge of viscosity is required for crystallization calculations. Many researchers have used the Turnbull equation in error. Subsequent nucleation and crystallization calculations can be in error by many orders of magnitude. This demonstrates the requirement for better methods for interpolating and extrapolating the viscosity of undercooled liquids. This is also true for undercooled water. Since amorphous water ice is the predominant form of water in the universe, astrophysicists have modeled the crystallization of amorphous water ice with viscosity relations that may be in error by five orders-of-magnitude.

Source record↗
DOE OSTI2023

Phonon Hall Viscosity of Ionic Crystals

When time-reversal symmetry is broken, the low-energy description of acoustic lattice dynamics allows for a dissipationless component of the viscosity tensor, the phonon Hall viscosity, which captures how phonon chirality grows with the wave vector. In this work, we show that, in ionic crystals, a phonon Hall viscosity contribution is produced by the Lorentz forces on moving ions. We calculate typical values of the Lorentz force contribution to the Hall viscosity using a simple square lattice toy model, and we compare it with literature estimates of the strengths of other Hall-viscosity mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of artificial viscosity on shocked particle-laden flows for staggered grid Lagrangian methods

Abstract Shocked particle-laden flows are important to many natural and industrial processes. When simulating these systems, artificial viscosity is often required to prevent numerical artifacts, such as ringing, from arising in the pressure and density fields. The linear and quadratic coefficients of the artificial viscosity determine the amount of smoothing that occurs in these fields. For particle-laden flows, however, many of the fluid–particle interaction forces, for example, the pressure gradient force and unsteady forces, depend on gradients in the fluid fields. Furthermore, while the shock passes over a particle, these forces can be more dominant than drag. This means that the artificial viscosity coefficients affect how a particle and fluid interact when simulating shocked particle systems. Here this effect is investigated for isolated particles and for a particle curtain using a staggered grid Lagrangian approach. The artificial viscosity coefficients have a significant impact on the maximum force that a fluid imparts to a particle, which is important for determining whether a particle will break up in response to the shock. Furthermore, it is found that the density ratio between the particle and the fluid is important in determining whether the artificial viscosity coefficients have a significant impact on the particle’s motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Effect of CMAS viscosity on the infiltration depth in thermal barrier coatings of different microstructures

Calcium magnesium aluminum silicate (CMAS) is one of the leading concerns for the gas turbine industry. The effects of CMAS viscosity and the coating microstructure on CMAS infiltration depth were explored by conducting a time dependent interaction study. Three CMAS compositions were used from literature, and their viscosities predicted through FactSage viscosity module were drastically different. The interaction was carried out on three different TBCs synthesized using the solution precursor plasma spray process (SPPS): two of the TBCs were made of yttrium aluminum garnet (YAG) having different microstructures that promote different modes of CMAS infiltration, and one TBC was made of gadolinium zirconate (GZO). All samples had stress relieving vertical cracks and different intensities of horizontally banded porosity, (inter pass boundaries IPBs). A concentration of 100 mg/cm 2 of CMAS was applied on the TBCs which were then subjected to a 5-minute interaction at 1300 °C. Samples were analyzed using scanning electron microscopy (SEM), electron dispersive X-ray spectroscopy (EDXS), and transmission electron microscopy (TEM). Low viscosity CMAS readily penetrated the TBCs while more viscous CMAS showed less penetration. The depth of CMAS infiltration depended on the coating microstructure. In the YAG with IPBs, the CMAS spread horizontally in the IPBs before infiltrating deeper, resulting in reduced infiltration depth compared to other samples in spite of having wider vertical cracks. TEM and EDXS analysis were performed to investigate the phases present in the CMAS-TBC interaction region in YAG. Two regions were chosen, the top TBC surface in direct contact with the sea of CMAS, and the region at the CMAS penetration ended within the coating. Finally, the results showed that no secondary phases like apatite were observed in YAG, thus it can be concluded that the arrest of CMAS happened solely because of CMAS viscosity and the short infiltration time.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The Sensitivity of Joint Inversions of Seismic and Geodynamic Data to Mantle Viscosity

Seismic tomography has revealed the existence of large-scale velocity heterogeneities in the mantle. The interpretation of seismic velocity anomalies in terms of temperature and chemical composition is nonunique. We use geodynamic observations including gravity, plate motions, dynamic topography, and excess ellipticity of the core-mantle boundary combined with seismic observations to investigate the thermo-chemical structure of the mantle through joint inversions. An outstanding issue, however, is the physical connection between mantle density anomalies and the surface geodynamic observations, which requires knowledge of the mantle viscosity structure. Here we perform joint inversions assuming different viscosity profiles and examine the dependence of the results on the viscosity. We first assume that mantle heterogeneity is due to thermal variations, which places a constraint on the relation between seismic velocity and density, and we subsequently relax the constraint to allow for potential nonthermal effects. In all of our joint inversions, a nonthermal origin of density anomalies is required to explain the geodynamic data, though the amount varies with the assumed viscosity structure. A common observation is a high-density chemical signal in the center of the large low-shear-velocity provinces at the base of the mantle resulting in a near neutral or slightly dense overall buoyancy there. Using the derived density models and their corresponding viscosity profiles, we also calculate instantaneous mantle flow fields. The predicted flow fields derived from joint inversions are generally similar but are quite different from flow fields using density models derived from a posteriori scaling of pure seismic tomography models.

58 GEOSCIENCES↗
DOE OSTI2023

Estimating viscosity of individual substrate-deposited particles from measurements of their height-to-width ratios

Airborne particles alter the radiative forcing of climate and have further consequences on air visibility, atmospheric chemistry, and human health. Recent studies reported the existence of highly viscous semisolid and even solid amorphous organic aerosol (OA) particles. Particle viscosity has an impact on the heterogeneous chemistry, gas-particle partitioning, and ice nucleation properties. Consequently, variations in particle viscosity must be considered when predicting the atmospheric impact of OA. Here, in this work, we use scanning electron microscopy (SEM) and scanning transmission X-ray microscopy (STXM) to estimate the viscosity of individual particles deposited on substrates based on their characteristic height-to-width ratios, which are affected by changes in morphology upon deposition. The height-to-width ratios obtained from SEM and STXM exhibit a strong correlation, demonstrating that both imaging approaches can be applied separately for viscosity assessment of the substrate-deposited particles. While these metrics are largely qualitative, this method enables rapid assessment of particle viscosity ranges, distinguishing between semisolid (>10 10 Pa·s), viscous (10 4– 10 8 Pa·s), and liquid (10°–10 1 Pa·s) particles within ensembles of ambient particles collected for microscopy studies.

54 ENVIRONMENTAL SCIENCES↗
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 OSTI2020

The Viscosity and Atomic Structure of Volatile-Bearing Melilititic Melts at High Pressure and Temperature and the Transport of Deep Carbon

Understanding the viscosity of mantle-derived magmas is needed to model their migration mechanisms and ascent rate from the source rock to the surface. High pressure–temperature experimental data are now available on the viscosity of synthetic melts, pure carbonatitic to carbonate–silicate compositions, anhydrous basalts, dacites and rhyolites. However, the viscosity of volatile-bearing melilititic melts, among the most plausible carriers of deep carbon, has not been investigated. In this study, we experimentally determined the viscosity of synthetic liquids with ~31 and ~39 wt% SiO 2 , 1.60 and 1.42 wt% CO 2 and 5.7 and 1 wt% H 2 O, respectively, at pressures from 1 to 4.7 GPa and temperatures between 1265 and 1755 °C, using the falling-sphere technique combined with in situ X-ray radiography. Our results show viscosities between 0.1044 and 2.1221 Pa·s, with a clear dependence on temperature and SiO 2 content. The atomic structure of both melt compositions was also determined at high pressure and temperature, using in situ multi-angle energy-dispersive X-ray diffraction supported by ex situ microFTIR and microRaman spectroscopic measurements. Our results yield evidence that the T–T and T–O (T = Si,Al) interatomic distances of ultrabasic melts are higher than those for basaltic melts known from similar recent studies. Based on our experimental data, melilititic melts are expected to migrate at a rate ~from 2 to 57 km·yr -1 in the present-day or the Archaean mantle, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Viscosity-temperature relation based on the evolution of medium-range structures of silica

The viscosity of glass is the most important technological property to glass manufactures and various applications. Practically, finding an accurate equation to express the glass viscosity behavior in the entire glass transition temperature range is tremendously challenge because it spans more than ten orders of magnitude. After a brief review of existing empirical viscosity equations, this work focuses on the correlating silica viscosity behavior with the evolution of glass medium-range structure, based on the recently proposed nanoflake model. From this new model, a new equation is constructed, which correctly describes the Arrhenius-type behavior of silica viscosity η above the melting temperature Tm, and non-Arrhenius-type behavior from Tm to a critical temperature Tc. At temperature lower than Tc, the equation predicts Arrhenius-type behavior again for the viscous flow with increased activation energy. The new equation agrees with experimental data in the entire temperature spanned from extremely high to extremely low. Furthermore, the application of the new equation is shown to extend to sodium silicate glasses as well.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Revisiting the validity of eddy viscosity models for predicting airflow over water waves

In this study, we revisit the validity of eddy viscosity models for predicting wave-induced airflow disturbances over ocean surface waves. We first derive a turbulence curvilinear model for the phase-averaged Navier–Stokes equations, extending the work of Cao, Deng & Shen (2020 J. Fluid Mech. 901, A27), by incorporating turbulence stress terms previously neglected in the linearised viscous curvilinear model. To verify our formulation, we perform a priori tests by numerically solving the model using mean wind and turbulence stress profiles from large-eddy simulations (LES) of airflow over waves across various wave ages. Results show that including turbulence stress terms improves wave-induced airflow predictions compared with the previous viscous curvilinear model. We further show that using a standard mixing-length eddy viscosity yields inaccurate predictions at certain wave ages, as it fails to capture wave-induced turbulence, which fundamentally differs from mean shear-driven turbulence. The LES data show that accurate representations of wave-induced stresses require a complex-valued eddy viscosity. The maximum magnitude of this eddy viscosity scales as ∼𝑢 𝜏 ⁢𝜁 𝑖𝑛𝑛𝑒𝑟 , where 𝑢 𝜏 is the friction velocity and 𝜁 𝑖𝑛𝑛𝑒𝑟 is the inner-layer thickness, the height at which the eddy-turnover time matches the wave advection time scale. This scaling aligns with the prediction by Belcher & Hunt (1993 J. Fluid Mech. 251, 109–148). Overall, the findings demonstrate that traditional eddy viscosity models are inadequate for capturing wave-induced turbulence. More sophisticated turbulence models are essential for the accurate prediction of airflow disturbances and form drag in wind–wave interaction models.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Numerical minimization of neoclassical poloidal viscosity for supersonic equilibria in tokamak geometry

Extensive experimental evidence has shown that the presence of poloidal flow in tokamaks can dramatically improve transport properties. However, theory indicates that poloidal flows are damped by poloidal viscosity, thus necessitating external drivers, such as neutral beam injection or radio frequency heating. In this work, ideal magnetohydrodynamic equilibria are calculated via the FORTRAN code FLOW [Guazzotto et al., Phys. Plasmas 11, 604 (2004)] and a postprocessor is used to estimate the neoclassical poloidal viscosity. The equilibrium inputs, which correspond to intuitive physical quantities, are then numerically optimized to reduce a viscosity figure of merit. We present supersonic equilibria in tokamak geometry with minimized neoclassical poloidal viscosities for various velocity free function inputs, plasma aspect ratios, and collisionality regimes. Benchmarks are made against an analytic theory as well as a classical expression of poloidal viscosity. Numerical confirmation of the analytic theory is obtained in the high aspect ratio and high collisionality limit. Good agreement is also seen near the plasma core and edge, with discrepancies arising in the intermediate region. Outside of these limits, rotation input function profiles are found that provide ∼order of magnitude improvements over the analytic theory, with additional progress being made toward predictions for tokamak-relevant equilibria.

Physics↗
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 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 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 OSTI2022

Verification study of xRAGE’s multi-ion viscosity model

This report on the adequacy of the plasma viscosity model implemented in LANL’s ASC program code, xRAGE is organized as follows. Section I provides a brief introduction on the kinetic theory which underpins plasma viscosity, and Section II describes how plasma kinetic theory can be used to calculate plasma viscosity coefficients accurate within the hydrodynamic limit. Section III describes the existing model for calculating the plasma viscosity coefficient in xRAGE. Next, Section IV shows tests of the xRAGE viscosity model against established hydrodynamic theory in the case of binary and trinary plasma mixtures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Optimization of artificial viscosity in production codes based on Gaussian Regression surrogate models

To accurately model flows with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity has to be introduced to convert kinetic energy into internal energy, thereby increasing the entropy across shocks. Determining the appropriate strength of the artificial viscosity is an art and strongly depends on the particular problem and experience of the researcher. The objective of this study is to pose the problem of finding the appropriate strength of artificial viscosity as an optimization problem and solve this problem using machine learning (ML) tools, specifically using surrogate models based on Gaussian Process regression and Bayesian analysis. We describe the optimization method and discuss various practical details of its implementation. The shock-containing problems for which we apply this method all have been implemented in the LANL code FLAG. First, we apply ML to find optimal values to isolated shock problems of different strengths. Second, we apply ML to optimize viscosity for a 1D propagating detonation problem based on Zel’dovich-von Neumann-Doring (ZND) detonation theory using a reactive burn model. We compare results for default (currently used values in FLAG) and optimized values of artificial viscosity for these problems demonstrating the potential for significant improvement in the accuracy of computations.

42 ENGINEERING↗
NASA NTRS1971

The viscosity of magmatic silicate liquids: A model for calculation

A simple model has been designed to allow reasonably accurate calculations of viscosity as a function of temperature and composition. The problem of predicting viscosities of anhydrous silicate liquids has been investigated since such viscosity numbers are applicable to many extrusive melts and to nearly dry magmatic liquids in general. The fluidizing action of water dissolved in silicate melts is well recognized and it is now possible to predict the effect of water content on viscosity in a semiquantitative way. Water was not incorporated directly into the model. Viscosities of anhydrous compositions were calculated, and, where necessary, the effect of added water and estimated. The model can be easily modified to incorporate the effect of water whenever sufficient additional data are accumulated.

Bottinga, Y.↗
NASA NTRS1979

Development of high viscosity coatings for advanced Space Shuttle applications

Laboratory studies for increasing the thermal resistance of high viscosity coatings for silica reusable surface insulation are presented. The coatings are intended for the reentry temperature associated with advanced Space Shuttle applications which will involve aerodynamic shear forces during entry from earth orbits. Coating viscosity was increased by (1) reduction in the concentration of the low viscosity additive B2O3; (2) reduction in the particle size of the constituent powders in coatings; and (3) addition of a high viscosity glass former (GeO2). A coating system was produced by combining the three methods which showed apparent higher viscosity than the current coating, while satisfying all the current Shuttle Orbiter coating requirements.

Garofalini, S. H.↗