Search NASA⌕ Search

SEARCH · Search NASA

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

Effective viscosities in a hydrodynamically expanding boost-invariant QCD plasma

Background: The near-equilibrium properties of a QCD plasma can be encoded into transport coefficients such as bulk and shear viscosity. In QCD, the ratio of these transport coefficients to entropy density, ζ/$\textit{s}$ and η/$\textit{s}$, depends nontrivially on the plasma's temperature. This is unlike in conformal systems where they take constant values such as η/$\textit{s}$ = 1/(4π). Purpose: In this work, we show that in a (0 + 1)D boost-invariant fluid with no transverse expansion, a temperature-dependent ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ can be described by an equivalent effective viscosity $\langleζ/s\rangle_{\text{eff}}$ or $\langleη/s\rangle_{\text{eff}}$. This effective viscosity combines the actual temperature-dependent ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ with the temperature profile of the fluid. Further, we extend the concept of effective viscosity in systems with transverse expansion and discuss how effective viscosities can be used to identify families of ζ/$\textit{s(T)}$ and η /$\textit{s(T)}$ that lead to similar hydrodynamic evolution. Methods: The Navier-Stokes relativistic hydrodynamic equations are used to provide a first definition of effective viscosity, in (0 + 1)D and (1 + 1)D. In the (0 + 1)D case, the analysis is extended to Israel-Stewart-type second-order hydrodynamics to clarify the effect of higher-order hydrodynamics corrections on the effective viscosity. Results: In a boost-invariant fluid with no transverse expansion [(0 + 1)D], the effective viscosity is expressed as a simple integral of ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ over temperature, with a weight determined by the speed of sound of the fluid. The result is general for any equation of state with a moderate temperature dependence of the speed of sound, including the QCD equation of state. This definition of effective viscosity can be used to identify infinite families of ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ that produce essentially indistinguishable temperature profiles. In a boost-invariant cylindrical system [(1 + 1)D], a similar definition of effective viscosity is obtained in terms of characteristic trajectories in time and transverse direction. This leads to an infinite number of constraints on an infinite functional space for ζ/$\textit{s(T)}$ and η/$\textit{s(T)}$. Realistic examples are presented by using a finite number of constraints on a finite functional space. Conclusions: The definition of effective viscosity in a (0 + 1)D system clarifies how infinite families of ζ/$\textit{s(T)}$ and η/$\textit{s(T)}$ can result in nearly identical hydrodynamic temperature profiles. By extending the study to a boost-invariant cylindrical [(1 + 1)D] fluid, we identify an approximate but more general definition of effective viscosity that highlights the potential and limits of the concept of effective viscosity in fluids with limited symmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Predictions of the glass transition temperature and viscosity of organic aerosols from volatility distributions

Abstract. Volatility and viscosity are important properties of organic aerosols (OA), affecting aerosol processes such as formation, evolution, and partitioning of OA. Volatility distributions of ambient OA particles have often been measured, while viscosity measurements are scarce. We have previously developed a method to estimate the glass transition temperature (Tg) of an organic compound containing carbon, hydrogen, and oxygen. Based on analysis of over 2400 organic compounds including oxygenated organic compounds, as well as nitrogen- and sulfur-containing organic compounds, we extend this method to include nitrogen- and sulfur-containing compounds based on elemental composition. In addition, parameterizations are developed to predict Tg as a function of volatility and the atomic oxygen-to-carbon ratio based on a negative correlation between Tg and volatility. This prediction method of Tg is applied to ambient observations of volatility distributions at 11 field sites. The predicted Tg values of OA under dry conditions vary mainly from 290 to 339 K and the predicted viscosities are consistent with the results of ambient particle-phase-state measurements in the southeastern US and the Amazonian rain forest. Reducing the uncertainties in measured volatility distributions would improve predictions of viscosity, especially at low relative humidity. We also predict the Tg of OA components identified via positive matrix factorization of aerosol mass spectrometer (AMS) data. The predicted viscosity of oxidized OA is consistent with previously reported viscosity of secondary organic aerosols (SOA) derived from α-pinene, toluene, isoprene epoxydiol (IEPOX), and diesel fuel. Comparison of the predicted viscosity based on the observed volatility distributions with the viscosity simulated by a chemical transport model implies that missing low volatility compounds in a global model can lead to underestimation of OA viscosity at some sites. The relation between volatility and viscosity can be applied in the molecular corridor or volatility basis set approaches to improve OA simulations in chemical transport models by consideration of effects of particle viscosity in OA formation and evolution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Source of Processable Vitrimer Viscosities: Swap Frequencies and Steric Factors

Vitrimers exhibit high, processable viscosities, where other polymers do not, and are among the most promising polymers for closed-loop material circularity. We sought to investigate the underlying chemical kinetic factors that result in high viscosities for vitrimers, which are crucial to designing vitrimers with tunable viscosity. To interrogate these factors, we achieved the first simulated predictions of real vitrimer viscosities, using a novel kinetic Monte Carlo molecular dynamics method, overcoming the time and length scale gaps to predict experimental bulk viscosities. The vitrimer architecture investigated is based on poly(dimethylsiloxane) chains and vinylogous urethane bond swaps. We probed the effects of the extent of free swapping groups, %F, the activation energy, E A , and the steric factor, ρ. The steric factor is related to the intrinsic reaction probability for molecules with sufficient energy. All three factors were found to be significant, but the role of ρ was found to be the biggest and also the most underappreciated. The results show that the inclusion of accurate ρ is of critical importance for viscosity predictions, with the evidence suggesting that the typical assumption of ρ = 1 is not valid for vitrimers and that, indeed, very low steric factors are present in bond-swap vitrimers such that values of ρ < 10 –10 may be typical. This greatly influences the bond exchange rates and, ultimately, the viscosities. Recognition of this result is necessary for the prediction of vitrimer viscosities from molecular simulations and to make vitrimers by design from molecular dynamics. We also investigated the effects that E A , ρ, and the number of free swapping groups have upon vitreous range temperatures, TV, with respect to achieving a specific viscosity (η V = 1 × 10 8 Pa·s), as well as for a commonly reported higher viscosity extrapolation (η V = 1 × 10 12 Pa·s). The evidence suggests that vitrimers may follow universal curves for E A vs T V , as a function of ρ. Finally, this study achieves the first of these comparisons of molecular simulations to experiments and reveals critical insights toward creating vitrimers by design, while providing a route for the prediction of T V from kinetic Monte Carlo molecular dynamics simulations.

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

Thermophysical modeling of raw glaze liquidus temperature and viscosity

Raw glazes, which do not contain pre-melted glass frits, are widely used due to their low cost. The natural variability of the mineral compositions used in raw glazes affects their viscosity during firing. Molten glaze viscosity plays an important role in determining the final surface quality; therefore, it is important to determine how raw glaze composition and firing conditions affect viscosity. Thermophysical modeling provides a way to evaluate the liquidus temperature and viscosity of molten glazes as a function of composition and oxygen activity. Equilibrium phase volume fraction and composition of 242 raw glaze recipes were evaluated under oxidizing and reducing conditions. Undissolved or precipitation of solid phases explain the secondary flux transition from anti-fluxing to fluxing near the liquidus temperature. The liquidus temperature and viscosity of iron glazes decrease as a function of increasing iron content under a reducing atmosphere. The empirical probability distribution of molten glaze viscosity follows a lognormal distribution, with lead glazes having significantly lower viscosity compared with glazes without lead. Here, the peak firing temperature viscosity of lead-free glazes is near the working point of glass. Multiple linear regression analysis shows that the peak firing temperature viscosity and liquidus temperature significantly predict the viscosity factor.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Universal Solvent Viscosity Reduction via Hydrogen Bonding Disruptors

Liquid Ion Solutions LLC (DBA RoCo Global) in partnership with Carnegie Mellon University and Carbon Capture Scientific LLC, has performed lab-scale development and evaluation of novel additives that lower the viscosity of water-lean amine solvents for post-combustion carbon dioxide capture. This project focuses on developing additives that minimize the formation of long-range electrostatic and hydrogen bonding (HB) networks, decreasing the solvent viscosity, improving diffusion, and improving the process economics. The project objectives included: 1) performing computer simulation to understand the molecular interactions of the additive molecules in water-lean CO 2 capture solvents, 2) design and synthesis of HB disruptors additives, 3) performance testing with additive molecules on model amine solvents, and 4) demonstration of the effectiveness of the optimized additives in the presence of synthetic flue gas. To meet the abovementioned objectives, the project team utilized a holistic approach that combines molecular simulation, experimental testing, and economic analysis studies. The project team developed ab initio molecular model and then perform computer simulation to develop relationship between hydrogen bonding, viscosity, and performed quantitative analysis of additive on the viscosity of the solvent. The team completed computational comparative study on a range of organic functional groups such as ethers, esters, cyclic carbonates, alkanes, and ammonium salts for their effect on viscosity gaining key insights into molecular interactions and the impact of various functional groups and molecular shapes on viscosity. Assisted with molecular simulation insights, the project team conducted additive synthesis and testing, including a proof-of-concept study, additive screening, optimization, and synthetic flue gas testing. The experimental proof-of-concept study proved that the hydrogen bonding acceptors result in significant decrease of viscosities. Detailed additive screening (exploring various functionalities and molecular structures) has been performed. Several promising additives showed excellent reduction in viscosity (30-41%) at 5% additive loading, and over 50% viscosity reduction at 10% additive loading for the model solvents. The team also performed complex screening studies on additive loadings and mixing effect among additives using the design of experiments. Based on multiple screening experiments, one additive-solvent candidate was down-selected for synthetic flue gas testing. A 100-hour continuous absorption/desorption study was conducted under simulated flue gas using a lab-scale continuous capture and separation system. No degradation (for both solvent and additive) was observed based on the GC results of the solvent samples collected from the continuous study. The team conducted preliminary engineering analyses and cost-benefit analyses to quantify the potential economic benefits of the additive approach for solvent viscosity reduction. Based on the experimental data, CO 2 capture cost savings from the capital and operating cost savings are estimated at $\$$4.7/tonne and $\$$0.3/tonne CO 2 captured, respectively. Considering the additive cost, the net benefit is estimated to be between $\$$4.32~$\$$4.86/tonne CO 2 captured.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

Reduced viscosity steadily weakens oceanic currents

The viscosity of both air and water is temperature dependent. A rising temperature leads to an increased viscosity for air but a decreased viscosity for water. As climate becomes warmer, this increased air viscosity can partly inhibit the reduction of wind stress over the ocean, and the reduced water viscosity causes less downward momentum and heat transport. As these opposing effects of warming on air and water viscosity are not included in the state-of-the-art climate models, the understanding of their potential impacts on the response of the climate system to the anthropogenic warming is lacking. Here, via analyzing the Simple Ocean Data Assimilation oceanic reanalysis dataset, we show that the ocean heat content increases at a rate of ~1.3 × 10 22 J/yr over 35 years, which leads to a continuous reduction of oceanic viscosity. As a result, the ocean vertical shear enhances with a shoaling of the mixed layer depth and a reduced vertical linkage in the ocean. Our calculations show a reduction of the oceanic kinetic energy at a rate of ~2.4 × 10 16 J/yr. Potentially, this could generate far-reaching impacts on the energy storage of the climate system and, hence, could pace the global warming. Thus, it is important to include the temperature-dependent viscosity in our climate models. Freshwater discharged from polar ice sheets and mountain glaciers also contributes to the reduction in oceanic viscosity but, at present, to a lesser extent than that in oceanic warming. Reduced oceanic viscosity, therefore, is an important, but hitherto overlooked, response to a warming climate and contributes to many recent weather extremes including heavier rainfall rates in hurricanes, slackening of the polar vortex, and oceanic heat waves.

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

Arrhenian to non-Arrhenian crossover in glass melt viscosity

The activation energy for glass melt viscosity is a function of temperature, i.e., non-Arrhenian, within a middle range of viscosity values (typically, between 10 2 and 10 10 Pa s) and constant, i.e., Arrhenian, outside this range. At the low-viscosity end, the Arrhenian to non-Arrhenian transition (termed the crossover) falls within the glass processing range: glass melting and fining occur within the Arrhenian range while glass forming occurs within the non-Arrhenian range. By the Adam-Gibbs equation, the configuration entropy is nearly constant at high temperatures (low viscosities), where the glass melt turns into a simple liquid, and decreases with decreasing temperature (increasing viscosity) as the glass is increasingly polymerized. In the present study, the power-law function and hyperbolic tangent function are used to express the configuration entropy as a function of temperature. The viscosity crossover is defined as the point at which onset of structural grouping is causing the configuration entropy to deviate from its high-temperature value. In conclusion, it appears that the limited range of the configuration entropy of the fully depolymerized melt imposes a restriction on the crossover viscosity; for well-defined glass families, such as float glasses or nuclear waste glasses, the crossover point can be defined using a fixed, i.e., composition-independent, viscosity value.

36 MATERIALS SCIENCE↗
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 OSTI2022

Error quantification of the Arrhenius blending rule for viscosity of hydrocarbon mixtures

Six hundred and seventy-five measurements of dynamic viscosity and density have been used to assess the prediction error of the Arrhenius blending rule for kinematic viscosity of hydrocarbon mixtures. Major trends within the data show that mixture complexity–binary to hundreds of components—and temperature are more important determinants of prediction error than differences in molecular size or hydrogen saturation between the components of the mixtures. Over the range evaluated, no correlation between prediction error and mole fractions was observed, suggesting the log of viscosity truly is linear in mole fraction, as indicated by the Arrhenius blending rule. Mixture complexity and temperature also impact molar volume and its prediction. However, a linear regression between the two model errors explains less than 20% of the observed variation, indicating that mixture viscosity and/or molar volume are not linear with respect to temperature and/or mixture complexity. Extensive discussion of the intermolecular forces and the geometric arrangement of molecules and vacancies in liquids, which ultimately determines its viscosity, is brought into context with the implicit approximations behind the Arrhenius blending rule. The complexity of this physics is not compatible with a simple algebraic correction to the model. However, sufficient data is now available to determine confidence intervals around the prediction of fuel viscosity based on its component mole fractions and viscosities. At -40°C, when all identified components are pure molecules the modeling error is 13.2% of the predicted (nominal) viscosity times the root mean square of the component mole fractions.

10 SYNTHETIC FUELS↗
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 OSTI2021

Estimation of secondary organic aerosol viscosity from explicit modeling of gas-phase oxidation of isoprene and α -pinene

Abstract. Secondary organic aerosols (SOA) are major components of atmospheric fine particulate matter, affecting climate and air quality. Mounting evidence exists that SOA can adopt glassy and viscous semisolid states, impacting formation and partitioning of SOA. In this study, we apply the GECKO-A (Generator of Explicit Chemistry and Kinetics of Organics in the Atmosphere) model to conduct explicit chemical modeling of isoprene photooxidation and α-pinene ozonolysis and their subsequent SOA formation. The detailed gas-phase chemical schemes from GECKO-A are implemented into a box model and coupled to our recently developed glass transition temperature parameterizations, allowing us to predict SOA viscosity. The effects of chemical composition, relative humidity, mass loadings and mass accommodation on particle viscosity are investigated in comparison with measurements of SOA viscosity. The simulated viscosity of isoprene SOA agrees well with viscosity measurements as a function of relative humidity, while the model underestimates viscosity of α-pinene SOA by a few orders of magnitude. This difference may be due to missing processes in the model, including autoxidation and particle-phase reactions, leading to the formation of high-molar-mass compounds that would increase particle viscosity. Additional simulations imply that kinetic limitations of bulk diffusion and reduction in mass accommodation coefficient may play a role in enhancing particle viscosity by suppressing condensation of semi-volatile compounds. The developed model is a useful tool for analysis and investigation of the interplay among gas-phase reactions, particle chemical composition and SOA phase state.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Toward closure between predicted and observed particle viscosity over a wide range of temperatures and relative humidity

Abstract. Atmospheric aerosols can exist in amorphous semi-solid or glassy phase states whose viscosity varies with atmospheric temperature and relative humidity. The temperature and humidity dependence of viscosity has been hypothesized to be predictable from the combination of a water–organic binary mixing rule of the glass transition temperature, a glass-transition-temperature-scaled viscosity fragility parameterization, and a water uptake parameterization. This work presents a closure study between predicted and observed viscosity for sucrose and citric acid. Viscosity and glass transition temperature as a function of water content are compiled from literature data and used to constrain the fragility parameterization. New measurements characterizing viscosity of sub-100 nm particles using the dimer relaxation method are presented. These measurements extend the available data of temperature- and humidity-dependent viscosity to −28 ∘C. Predicted relationships agree well with observations at room temperature and with measured isopleths of constant viscosity at ∼107 Pa s at temperatures warmer than −28 ∘C. Discrepancies at colder temperatures are observed for sucrose particles. Simulations with the kinetic multi-layer model of gas–particle interactions suggest that the observed deviations at colder temperature for sucrose can be attributed to kinetic limitations associated with water uptake at the timescales of the dimer relaxation experiments. Using the available information, updated equilibrium phase-state diagrams (-80∘C

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Anomalously strong viscosity behavior in mixed oxy-sulfide Na 4 P 2 S 7-x O x invert glasses

In this paper, the first observations have been made of anomalously strong viscosity behavior and correspondingly strong bulk glass formation and resistance to crystallization in a melt quenched (MQ) prepared series of invert, modifier content greater than the glass former content, mixed oxy-sulfide Na 4 P 2 S 7-x O x , 0 ≤ x ≤ 5, (NaPSO) glasses. We find that these glasses even with 67 mole % modifier, Na 2 S, and only 33 mole% glass former P 2 S 5-x O x , exhibit a calorimetric fragility index, m cal , as small as 29 compared to typical values of ~ 60 to 80 for alkali silicate glasses with > 50 mole% SiO 2 and the theoretical minimum m cal of ~ 15. The m cal values were converted to kinetic viscosity fragility, m vis , values and used with the measured glass transition temperatures, T g s, to produce viscosity curves for the glasses using the Mauro-Yue-Ellison-Gupta-Allan (MYEGA) model. The viscosity of the glass with a minimum m cal value of x = 2.5 was then shown to exhibit a viscosity that was as strong as, and for some compositions stronger than, the viscosity of alkali aluminosilicate glasses which have three dimensional network structures compared to the one dimensional chain structures observed in these MOS invert glasses. Consistent with these small values of fragility indices and therefore predicted higher viscosity above the glass transition temperature, T g , in the supercooled liquid regime, these glasses are also particularly resistant to crystallization. For example, the x = 2.5 glass composition exhibits a crystallization temperature, T C , more than 150K above its T g . We find that the calorimetric fragility of these glasses varies systematically with the extent of network forming bridging oxygens (BOs) in the glasses. The structural change caused by incrementally replacing sulfur with oxygen has been found to have consistent, large effects on the thermal properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accurate Machine Learning for Predicting the Viscosities of Deep Eutectic Solvents

Deep eutectic solvents (DESs) are emerging as environmentally friendly designer solvents for mass transport and heat transfer processes in industrial applications; however, the lack of accurate tools to predict and thus control their viscosities under both a range of environmental factors and formulations hinders their general application. While DESs may serve as designer solvents, with nearly unlimited combinations, this unfortunately makes it experimentally infeasible to comprehensively measure the viscosities of all DESs of potential industrial interest. To assist in the design of DESs, we have developed several new machine learning (ML) models that accurately and rapidly predict the viscosities of a diverse group of DESs at different temperatures and molar ratios using, to date, one of the most comprehensive data sets containing the properties of over 670 DESs over a wide range of temperatures (278.15–385.25 K). Three ML models, including support vector regression (SVR), feed forward neural networks (FFNNs), and categorical boosting (CatBoost), were developed to predict DES viscosity as a function of temperature and molar ratio and contrasted with multilinear and two-factor polynomial regression baselines. Further, quantum chemistry-based, COSMO-RS-derived sigma profile (σ-profile) features were used as inputs for the ML models. The CatBoost model is excellent at externally predicting DES viscosity, as indicated by high R 2 (0.99) and low root-mean-square-error (RMSE) and average absolute relative deviations (AARD) (5.22%) values for the testing data sets, and 98% of the data points lie within the 15% of AARD deviations. Furthermore, SHapley additive explanation (SHAP) analysis was employed to interpret the ML results and rationalize the viscosity predictions. The result is an ML approach that accurately predicts viscosity and will aid in accelerating the design of appropriate DESs for industrial applications.

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

Hanford low-activity waste glass composition-temperature-melt viscosity relationships

This study developed a model for predicting viscosity of alkali-alumino-borosilicate glass melts as functions of composition and temperature. The model is based on a total of 3935 viscosity-temperature data from 574 glasses with viscosity values ranging from 2.53 to 7260 Poise (P) in the temperature range of 900–1260°C. Several different model forms were surveyed, including those based on Arrhenius, Vogel-Fulcher-Tammann, Avramov-Milchev, and Mauro-Yue-Ellison-Gupta-Allen. For each of these models, combinations of the temperature-independent parameters were fitted to composition. It was found that generally fitting more than one temperature-independent parameter as functions of composition resulted in overfitting. The Avramov-Milchev-based model was found to best represent the Hanford low-activity waste glass melt viscosity data based on model fit and validation statistics. A 21-term partial quadratic mixture model was recommended for use. This model predicts melt viscosity with a root-mean square error of .1736 ln(P), which is similar to the error in viscosity measurements from replicate glass analyses of .1383 ln(P). Viscosity was found to be most increased by SiO 2 > Al 2 O 3 > ZrO 2 > SnO 2 and most decreased by Li 2 O > Na 2 O > B 2 O 3 > CaO > K 2 O > MgO, at temperatures from 900 to 1260°C.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Drag viscosity of metals and its connection to Coulomb drag

Recent years have seen a surge of interest in studies of hydrodynamic transport in electronic systems. We investigate the electron viscosity of metals and find a component that is closely related to Coulomb drag. By using linear-response theory, viscosity, which is a transport coefficient for momentum, can be extracted from the retarded correlation function of the momentum flux, i.e., the stress tensor. There exists a previously overlooked contribution to the shear viscosity from the interacting part of the stress tensor which accounts for the momentum flow induced by interactions. This contribution, which we dub drag viscosity, is caused by the frictional drag force due to long-range interactions. It is therefore linked to Coulomb drag which also originates from the interaction-induced drag force. Starting from the Kubo formula and using the Keldysh technique, we compute the drag viscosity of two- and three-dimensional metals along with the drag resistivity of double-layer two-dimensional electronic systems. Both the drag resistivity and drag viscosity exhibit a crossover from quadratic-in-$\textit{T}$ behavior at low temperatures to a linear behavior at higher temperatures. Although the drag viscosity appears relatively small compared with the normal Drude component for the clean metals, it may dominate hydrodynamic transport in some systems, which are discussed in the conclusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Insight into the Relationship Between Viscosity and Structure of CaO-SiO{sub 2}-MgO-Al{sub 2}O{sub 3} Molten Slags

This article elucidates the quantitative relationship between viscosity and structure in a basic slag system of CaO-SiO{sub 2}-MgO-Al{sub 2}O{sub 3} and focuses on the role of Al{sub 2}O{sub 3}. Slag viscosity was measured by the rotating cylinder method, and structural information was obtained using Fourier transformation infrared, Raman and magic angular spinning nuclear magnetic resonance (MAS-NMR) techniques. The results show that, as the Al{sub 2}O{sub 3} content increased, slag viscosity increased initially and decreased afterwards, directly indicating that Al{sub 2}O{sub 3} had an amphoteric effect on slag viscosity. The Raman spectra verified that with increasing Al{sub 2}O{sub 3} content, the concentrations of Q{sup 0}(Si) and Q{sup 2}(Si) decreased first and then increased, while that of Q{sup 1}(Si) kept increasing and that of Q{sup 3}(Si) increased first and then decreased. The {sup 27}Al MAS-NMR spectra proved that the mole ratios of AlO{sub 5} and AlO{sub 6} to AlO{sub 4} kept increasing with the increase of Al{sub 2}O{sub 3} content, and, overall, Al{sub 2}O{sub 3} changed from a network former to a network modifier. The relationship between the viscosity and structure of the molten slags was further analyzed quantitatively based on the modified (NBO/T), denoted as (NBO/T)′, and we found a fine linear correlation between the logarithm of viscosity and (NBO/T)′. Moreover, the variations of thermodynamic properties of this system also indirectly supported the present experimental results.

36 MATERIALS SCIENCE↗