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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Stockmayer fluid simulations for viscosity and glass transition temperature of ionic liquids

We develop a Stockmayer fluid model for molecular dynamics simulations of ionic liquids that captures molecular polarization, ionic conductivity, viscosity, and glass transition temperature, using ethylammonium nitrate (EAN) as an example. The ions in EAN are treated as spheres interacting via the Lennard-Jones potential with an embedded point charge and a permanent dipole moment. We show that our simulation results for EAN are consistent with experimental data and then explore the effects of the molecular parameters on the viscosity of ionic liquids. Our results indicate that viscosity monotonically increases with ionic charge and dipole moment but non-monotonically changes with ionic diameter (or molar volume). This non-monotonic trend arises from the competition among the electrostatic interactions, molecular packing, and size asymmetry between the cation and anion. In conclusion, our model also shows that long-lived ion pairs result in higher viscosities.

Coarse-grained simulations↗

DOE OSTI2020

A small-volume, high-throughput approach for surface tension and viscosity measurements of liquid fuels

An ongoing effort looking at bioblendstocks to function as drop-in replacements or blending components for gasoline has identified a large number of candidate fuels. This work documents an approach for rapid screening of candidate fuels using relatively small sample sizes ($\mathcal{O}(\mu\,{\rm l})$) that targets two key physical properties of liquid fuels—surface tension and viscosity. The approach utilizes shape oscillation dynamics of single droplets generated by a piezo-electric device and their decay over time. Strobed imaging of the oscillation process is used along with image processing, edge detection, and data analysis to capture the decay of the oscillation over time. The time constant of the decay process along with oscillation frequency are then used to estimate viscosity and surface tension using a theory for small amplitude droplet oscillations. Measurements are obtained for primary reference fuels (isooctane and n-heptane) as well as candidate fuels from four bio-derived functional groups of interest. Measurement results for surface tension and viscosity are correlated with literature data as well as measurements from standard reference instruments. The measurement results show that the droplet oscillation based approach is capable of reproducing surface tension and viscosity values for the tested fuels within deviations of 7% and 13% respectively from literature data. Results are obtained using an average of 5 µl per fuel within about 20 s, thus demonstrating a small-volume, high-throughput approach that can be used for screening of candidate bioblendstock fuels.

09 BIOMASS FUELS↗

DOE OSTI2024

Bulk viscosity of the rigid rotor one-component plasma

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Effect of Oil Viscosity and Impact of Ionic Liquid Additive on Electrically Induced Pitting in Rolling Contact

Electrically induced bearing damage (EIBD) is a growing concern in electric drivetrains, where electrical discharge currents generate localized pitting that accelerates wear and premature failure. While ionic liquids (ILs) have shown promise as lubricant additives due to their natural physical adsorption, strong capabilities of tribofilm formation, and high molecular tunability, their role in mitigating EIBD is little known. This gap is especially evident because lubricant viscosity and film thickness affect discharge behavior, and additive performance in these regimes has not been well explored. In this study, tribological tests were performed using a ball-on-disk pure rolling-contact system under an applied voltage on polyalphaolefin oils of three viscosities (4, 10, and 150 cSt measured at 100 °C), with and without the addition of a phosphonium-phosphate IL. Surface damage was analyzed using scanning electron microscopy (SEM) and stylus profilometry. Damage quantification was based on pit morphology analysis of SEM images. The results provide new insights into how oil viscosity governs the discharge behavior and how an IL additive could influence the pit formation. While viscosity seems to play a significant role, an IL additive shows potential to reduce EIBD and deserves further study.

electrically induced bearing damage (EIBD)↗

DOE OSTI2022

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

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

79 ASTRONOMY AND ASTROPHYSICS↗

DOE 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 OSTI2023

Correlating viscosity and die swell in 3D printing of polyphenylsulfone: A thermo-mechanical optimization modus operandi

Additive manufacturing (AM) of polyphenylsulfone (PPSU), an amorphous high-performance polymer (HPP), has shown promise in exhibiting structural integrity for many applications. With thermoplastic material extrusion 3D printing, parts consisting of layered melts exhibited voids, where the adhesion quality and void size depend on the processing parameters implemented during printing. However, the structure-process-property (SPP) relationships between the thermal properties and rheological and mechanical behavior of AM PPSU are poorly understood. Here, in this study, a comprehensive analysis of the thermal, thermo-mechanical, and layer-by-layer build-up of PPSU is investigated and correlated with viscosity and extrudate die swell behavior with void fraction volume as confirmed by high-resolution x-ray computed tomography (XCT). Gas chromatography–mass spectrometry (GC–MS) pyrolysis confirmed PPSU composition and degradation at 480 °C, consistent with thermo-gravimetric analysis TGA. A 141% difference in viscosity was observed when processing temperatures increased from 350 °C to 400 °C. Likewise, extrudate die swell decreased with increasing temperatures, resulting in a 17% difference. Tensile and compressive properties confirmed that higher extruding temperatures resulted in lower porosity and higher mechanical strength. In general, higher extruding temperatures lower the viscosity to allow for stronger cohesion between layers; however, increasing the temperature too high results in a below unity extrudate die swell, reducing mechanical properties. Therefore, to optimize the mechanical integrity of 3D printed amorphous PPSU parts, an increase in the viscosity, near unity of the extrudate die swell, and tempered extrusion temperature are needed.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Odd-viscosity-induced passivation of Moffatt vortices

Odd viscosity is associated with a number of striking physical phenomena: a disk rotating in a viscous liquid experiences a normal compressive or tensile stress depending on its sense of rotation; chiral active matter interfacial deformations level off in a manner akin to surface tension; and unstable thermocapillary films become stabilized. In this article we show that odd viscosity decreases the size and intensity of Moffatt vortices in a corner formed by a solid substrate and a free interface. The effect is present because odd viscosity creates a normal-stress-induced shear stress at the liquid–gas interface. The new shear stress is out of phase and oppositely directed with respect to the interfacial liquid velocity. It squeezes the vortices and damps their rotation. As a consequence, Moffatt vortices can now be observed with greater facility in an experiment. Odd viscosity also increases the range of corner angles over which Moffatt vortices can be observed, and, in general, it shows a tendency to moderate singular behaviour.

Mechanics↗

DOE OSTI2023

Ultra-Low Viscosity and High Magnetic Susceptibility Magnetic Ionic Liquids Featuring Functionalized Diglycolic Acid Ester Rare-Earth and Transition Metal Chelates

Magnetic ionic liquids (MILs) comprise a subcategory of ionic liquids (ILs) and contain a paramagnetic metal center allowing them to be readily manipulated by an external magnetic field. While MILs are popularly employed as solvents in catalysis, separations, and organic synthesis, most low viscosity combinations possess a hydrophilic character that limits their use in aqueous matrices. To date, no study has reported the synthesis and characterization of hydrophobic MILs with viscosities similar to those of hydrophilic MILs and organic solvents while simultaneously exhibiting enhanced magnetic and thermal properties. In this study, diglycolic acid esters are employed as ligands to chelate with paramagnetic metals to produce cations that are paired with metal chelates composed of hexafluoroacetylacetonate ligands to form MILs incorporating multiple metal centers in the cation and anion. Viscosity values below 31.6 cP were obtained for these solvents, the lowest ever reported for hydrophobic MILs. Solubilities in nonpolar solvents such as benzene were observed to be as high as 50% (w/v) MIL-to-solvent ratio while being insoluble in water at concentrations as low as 0.01% (w/v). Effective paramagnetic moment values for these solvents ranged from 5.33 to 15.56 Bohr magnetons (μB), with mixed metal MILs containing multiple lanthanides in the anion generally offering higher magnetic susceptibilities. MILs composed of ligands containing octyl substituents were found to possess thermal stabilities up to 190 °C. The synthetic strategies explored in this study exploit the highly tunable nature of the employed cation and anion pairs to design versatile ultra-low viscosity magnetoactive solvents that possess tremendous potential and applicability in liquid–liquid separation systems, catalysis, and microfluidics where the mechanical movement of the solvent can be easily facilitated using electromagnets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Predicting the influence of particle size on the glass transition temperature and viscosity of secondary organic material

Abstract Atmospheric aerosols can assume liquid, amorphous semi-solid or glassy, and crystalline phase states. Particle phase state plays a critical role in understanding and predicting aerosol impacts on human health, visibility, cloud formation, and climate. Melting point depression increases with decreasing particle diameter and is predicted by the Gibbs–Thompson relationship. This work reviews existing data on the melting point depression to constrain a simple parameterization of the process. The parameter $$\xi $$ ξ describes the degree to which particle size lowers the melting point and is found to vary between 300 and 1800 K nm for a wide range of particle compositions. The parameterization is used together with existing frameworks for modeling the temperature and RH dependence of viscosity to predict the influence of particle size on the glass transition temperature and viscosity of secondary organic aerosol formed from the oxidation of $$\alpha $$ α -pinene. Literature data are broadly consistent with the predictions. The model predicts a sharp decrease in viscosity for particles less than 100 nm in diameter. It is computationally efficient and suitable for inclusion in models to evaluate the potential influence of the phase change on atmospheric processes. New experimental data of the size-dependence of particle viscosity for atmospheric aerosol mimics are needed to thoroughly validate the predictions.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Effects of electron viscosity on resonant layer responses to non-axisymmetric magnetic perturbations

The resonant field penetration to magnetic islands is the central MHD mechanism of non-axisymmetric plasma responses in a tokamak such as disruptive locking or favorable ELM stabilization. The resonant field penetration can be induced by any non-ideal processes as manifested in the delicate balance under the generalized Ohm's law. Here, we show that the viscous effects by electrons are not ignorable in the field penetration unlike previous presumption, even if the electron viscosity is as small as the square root of its mass compared to the ions. It is clear that its effects become only bigger if the electron viscosity becomes anomalously large. The work strictly follows the three-field model in the linear regime targeting the prediction of the onset of the field penetration and successfully extending it with electron viscosity and identifying new regimes. The results also indicate that the error field thresholds become more strongly dependent on plasma density than ones predicted in the linear regimes without the electron viscosity, which is consistent with experimental observations and thus a significant implication.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Reference Correlations for the Density and Viscosity of Molten Alkali and Alkaline Earth Fluoride Salts

While there is a significant body of literature pertaining to thermophysical property measurements of molten salts, there is often a wide degree of variability among independent measurements of the same compounds. As such, the scientific community benefits greatly from an unbiased, independent assessment of duplicate datasets, so that reference correlations which describe these thermophysical properties as functions of temperature can be determined and then commonly used by researchers, scientists, and engineers. With regard to molten fluoride compounds, a significant time has elapsed since density and viscosity reference correlations have been determined; Janz conducted the most recent effort, in 1988, to provide reference correlations for the densities and viscosities of molten fluoride compounds via the National Standard Reference Data System coordinated by the National Bureau of Standards. Since then, new data have been published for molten fluoride compounds, and a new precedent has surfaced for putting forth reference correlations that involve fitting to multiple primary datasets. In this work, reference correlations are put forth for molten alkali and alkaline earth fluoride compounds in an effort to provide updated, improved correlations for general use. For molten alkali fluoride densities, estimated uncertainties with a 95% confidence interval are summarized as follows: LiF (0.63%), NaF (0.48%), KF (0.76%), RbF (0.93%), and CsF (0.75%). For molten alkaline earth fluoride densities, an estimated uncertainty was not able to be quantified for BeF 2 because of limited data; however, estimated uncertainties with a 95% confidence interval are summarized as follows for the remaining alkaline earth fluorides: MgF 2 (1.5%), CaF 2 (0.92%), SrF 2 (1.6%), and BaF 2 (0.23%). For molten alkali fluoride viscosities, uncertainty was not able to be quantified for RbF and CsF because of limited data; however, estimated uncertainties with a 95% confidence interval are summarized as follows for the remaining alkali fluorides: LiF (4.4%), NaF (3.0%), and KF (4.0%). For molten alkaline earth fluoride viscosities, limited consistent data resulted in the recommendation of single datasets (from literature) that are deemed to be the most trustworthy based on the quality of the underlying experimental studies.

Birri, A. [Oak Ridge National Laboratory (ORNL), O↗

DOE OSTI2024

Viscosity bounds in liquids with different structure and bonding types

Recently, it was realized that liquid viscosity has a lower bound which is nearly constant for all liquids and is governed by fundamental physical constants. This was supported by experimental data in noble and molecular liquids. Here, we perform large-scale molecular dynamics simulations to ascertain this bound in two other important liquid types: the ionic molten salt system LiF and metallic Pb. We find that these ionic and metallic systems similarly have lower viscosity bounds corresponding to the minimum of kinematic viscosity of ∼ 10 − 7 m 2 / s . We show that this agrees with experimental data in other systems with different structures and bonding types, including noble, molecular, metallic, and covalent liquids. This expands the universality of viscosity bounds into the main system types known. Published by the American Physical Society 2024

Materials Science↗

DOE OSTI2019

Nanoscale viscosity of confined polyethylene oxide

Complex fluids near interfaces or confined within nanoscale volumes can exhibit substantial shifts in physical properties compared to bulk, including glass transition temperature, phase separation, and crystallization. Because studies of these effects typically use thin film samples with one dimension of confinement, it is generally unclear how more extreme spatial confinement may influence these properties. In this work, we used x-ray photon correlation spectroscopy and gold nanoprobes to characterize polyethylene oxide confined by nanostructured gratings (<100nm width) and measured the viscosity in this nanoconfinement regime to be ~500 times the bulk viscosity. Here, this enhanced viscosity occurs even when the scale of confinement is several times the polymer's radius of gyration, consistent with previous reports of polymer viscosity near flat interfaces.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Stokesian dynamics with odd viscosity

Stokesian dynamics is a well-established computational method for simulating dynamics of many particles suspended in a conventional passive fluid medium. Active fluids composed of self-driving particles with broken time-reversal symmetry permit the emergence of a so-called odd viscosity. Here, in this work, we extended the conventional Stokesian dynamics formalism to incorporate the additional hydrodynamic effects due to odd viscosity, which enables simulating collective behavior of many particles suspended in an active fluid medium with both even viscosity and odd viscosity.

42 ENGINEERING↗

DOE OSTI2020

Viscosity of glass-forming melt at the bottom of high-level waste melter feed cold caps: Effects of temperature and incorporation of solid components

At the final stages of conversion of melter feed (glass batch) to molten glass, the transient glass-forming melt becomes a continuous liquid phase encapsulating dissolving solid particles and gas bubbles that produce primary foam at the bottom of the cold cap (the reacting melter feed in an electric glass melting furnace). The glass-forming melt viscosity plays a dominant role in primary foam formation, stability, and eventual collapse, thus affecting the rate of melting (the glass production rate per cold cap area). For several melter feeds designed for nuclear waste vitrification, we have traced the glass-forming melt viscosity during the final stages of feed-to-glass conversion as it changes in response to changing temperature and composition (resulting from dissolving solid particles). Starting with a relatively low values at the moment when the melt connects, melt viscosity reaches maximum within the primary foam layer and then decreases to its final melter-operating temperature value. We paid a particular attention to the cold-cap bottom—the boundary between the primary foam layer and the thermal boundary layer—where the melt viscosity affects the rate of melting predominantly through its effect on the temperature at which primary foam is collapsing.

Lee, Seung Min↗

DOE OSTI2022

FY2022 Progress Report on Viscosity and Thermal Conductivity Measurements of Molten Salts

This report describes the development of and experimental measurements made with the thermal conductivity and viscosity measurement systems for molten salts at Oak Ridge National Laboratory. The thermal conductivity system is based on a steady-state measurement technique in which a heat flux is driven across a variable gap, and the changes in temperature difference across the gap are measured as the gap size is varied. The system comprises an inner containment, which possesses the heating element, and an outer containment, which possesses the salt specimen and cooling channels. The gap is formed between the outer bottom of the inner containment and the inner bottom of the outer containment. The viscosity system is based on the falling-ball technique in which a salt-filled tubular crucible is fixed at an angle and a ball rolls through the salt. The terminal velocity of the ball can be converted into viscosity via a calibration factor calculated for crucible-ball combinations which account for thermal expansion and flow characteristics. Modifications were made to both systems this fiscal year. Modifications to the thermal conductivity system included adding guard heaters to prevent axial heating losses, and modifications to the viscosity system included changing the crucible design to prevent air pocket formation and providing the ability to conduct measurements with a glass or steel crucible.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

High-Fidelity Simulations of Gravity Currents Using a High-Order Finite-Difference Spectral Vanishing Viscosity Approach

This numerical work investigates the potential of a high-order finite-difference spectral vanishing viscosity approach to simulate gravity currents at high Reynolds numbers. The method introduces targeted numerical dissipation at small scales through altering the discretisation of the second derivatives of the viscous terms in the incompressible Navier-Stokes equations to mimic the spectral vanishing viscosity (SVV) operator, originally designed for the regularisation of spectral element method (SEM) solutions of pure advection problems. Using a sixth-order accurate finite-difference scheme, the adoption of the SVV method is straightforward and comes with a negligible additional computational cost. In order to assess the ability of this high-order finite-difference spectral vanishing viscosity approach, we performed large-eddy simulations (LES) of a gravity current in a channelised lock-exchange set-up with our SVV model and with the well-known explicit static and dynamic Smagorinsky sub-grid scale (SGS) models. The obtained data are compared with a direct numerical simulation (DNS) based on more than 800 million mesh nodes, and with experimental measurements. A framework for the energy budget is introduced to investigate the behaviour of the gravity current. First, it is found that the DNS is in good agreement with the experimental data for the evolution of the front location and velocity field as well as for the stirring and mixing inside the gravity current. Secondly, the LES performed with less than 0.4% of the total number of mesh nodes compared to the DNS, can reproduce the main features of the gravity currents, with the SVV model yielding slightly more accurate results. It is also found that the dynamic Smagorinsky model performs better than its static version. For the present study, the static and dynamic Smagorinsky models are 1.8 and 2.5 times more expensive than the SVV model, because the latter does not require the calculation of explicit SGS terms in the Navier-Stokes equations nor spatial filtering operations.

direct numerical simulation↗