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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

On the applicability of the Redlich-Kister framework for viscosity estimation of molten halide salt mixtures

For molten halide salt mixtures already being utilized or under consideration for carbon-free energy production systems, it is crucial that their viscosity is well understood so that system thermal hydraulics can be reliably assessed. Because of the difficulty in accurately measuring molten halide viscosity and the sheer size of the matrix of possible higher order salt mixtures that may be of interest to the energy industry, there are several gaps in the quantified understanding of molten halide viscosity across this matrix. As such, both first-principles and semi-empirical modeling techniques may be crucial for rapidly assessing this broad, complex compositional domain. Herein, the Redlich-Kister framework is applied to assess the feasibility of broadly interpolating and estimating the viscosity of several pseudobinary and pseudoternary molten halide salt systems that may be of key interest to the energy industry. The framework is based on the assumption that an ideal component and a nonideal component collectively describe the viscosity as a function of composition and temperature for a given molten halide system. Three different ideal models were considered for the ideal component, including Grunburg-Nissan, Katti-Chaudhri, and Gambill methods. Regarding the pseudobinary interpolations, the Redlich-Kister models with either the Grunburg-Nissan or Katti-Chaudhri models as the ideal component resulted in either highly (average error less than 5%) or reasonably (average error less than 15%) accurate interpolations of pseudobinary halide viscosity; BeF 2 - or UF 4 -bearing salts tended to result in reasonably accurate interpolations, whereas other pseudobinary mixtures tended to show high accuracy. Regarding the pseudoternary extrapolations, the Redlich-Kister framework shows reasonable success at estimating the extent to which a pseudoternary system may indicate deviations from ideal Grunburg-Nissan mixing, where discrepancies with comparative experimental data generally stay within 30%. Finally, the primary reasons identified for such discrepancies are (1) inaccuracy in the underlying experimental data, (2) different complexation behavior in the higher order systems compared to the pseudobinary subsystems, and (3) extrapolation into temperatures too far out of the domain, which is valid for the underlying experimental data feeding the Redlich-Kister model.

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

Viscosity of transient glass-forming melt and its relation to foaming during batch-to-glass conversion

Primary foam, which affects the heat transfer into the glass batch and the final glass quality, occurs when a sufficient quantity of transient glass-forming melt evolves with viscosity low enough to close the open porosity of the reacting glass batch. To better understand how the fraction of transient melt and its viscosity affect the primary foam temperature range, we determined, with x-ray diffraction, the fraction and composition of the transient glass-forming melt in a heated waste glass melter feed as a function of temperature. Then we prepared a set of transient melts that occurred within the foaming temperature interval and measured their viscosities with spindle and falling sphere viscometers. Further, the results agree with the Adam-Gibbs and VFT viscosity-composition models, even for transient melt compositions outside of the models compositional validity range. As silica and other refractory particles dissolved in the predominantly borate transient melt while temperature increased, viscosity increased from the initial value of ~500 Pa s at the onset of foaming (~650°C) to a maximum of ~770 Pa s when silica dissolution was almost complete (~700°C). As temperature increased further, transient melt viscosity decreased to ~220 Pa s (~850°C) when the primary foam collapsed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Distinct Gas-Particle Partitioning and Viscosity Characteristics of Secondary Organic Aerosols Derived from α-Pinene versus Ocimene

Secondary organic aerosols (SOA) have complex, multicomponent composition that controls particle viscosity and gas-particle partitioning, key factors to their atmospheric evolution. This study investigates the chemical composition, volatility and viscosity of SOA formed by ozonolysis of cyclic α-pinene (PSOA) and acyclic ocimene (OSOA) monoterpenes. Using Temperature-Programmed Desorption combined with Direct Analysis in Real-Time ionization and High-Resolution Mass Spectrometry, we determined the molecular composition and saturation mass concentration of individual SOA constituents. These data enabled gas-particle partitioning and viscosity estimates under varied atmospheric conditions. PSOA, composed of higher molecular weight and less oxidized species, shows higher condensability and viscosity under high total organic mass (tOM) loadings. Here, in contrast, OSOA, consisting of more oxidized, lower molecular weight species, exhibits greater sensitivity to tOM, with viscosity increasing significantly upon dilution. Poke-flow experiments support this trend, indicating that OSOA undergoes more dynamic compositional and phase changes during atmospheric aging. These observations reveal distinct dynamic trends in the atmospheric transformations and reactivity of SOA from cyclic and acyclic monoterpenes, with the latter showing greater compositional changes during aging that alter viscosity and diffusion. This highlights the importance of incorporating such dynamic transformations into atmospheric models to improve predictions of SOA atmospheric loadings, lifetimes, and impacts.

cyclic and acyclic monoterpenes↗
DOE OSTI2024

High-Throughput Screening and Accurate Prediction of Ionic Liquid Viscosities Using Interpretable Machine Learning

Ionic liquids (ILs) are a novel group of green solvents with great promise for various industrial applications, including carbon capture and lignocellulosic biomass deconstruction. However, the use of ILs at the industrial scale remains challenging due to their high viscosities at ambient temperatures. To develop ILs with lower viscosities, a systematic study of their quantitative structure–property relationship (QSPR) is desirable. Here, we developed four machine learning (ML) models to predict viscosity at various temperature and pressure ranges, trained over a wide range of ILs consisting of various cationic and anionic families. ML methods including two-factor polynomial regression (two-factor PR), support vector regression (SVR), feed-forward neural networks (FFNN), and categorical boosting (CATBoost) were developed based on features that have proven useful in previous ML studies: COSMO-RS (conductor-like screening model for real solvents)-derived surface screening charge densities (sigma profiles). FFNN and CATBoost were the most accurate in predicting IL viscosities with lower average absolute relative deviation and higher R2 values on the test set. Tanimoto similarity scores were calculated to characterize the chemical space and structural similarity of the investigated ions. Furthermore, SHapley Additive exPlanation (SHAP) analysis was employed to interpret the ML results. Temperature, the polar area of ILs, and the nonpolar regions of ions are key features that influence the viscosity predictions. Importantly, the IL viscosity prediction here is the most accurate reported to date.

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

High-throughput small-angle X-ray scattering reveals effective structure factor transitions linked to high-concentration antibody viscosity

High-concentration monoclonal antibody (mAb) formulations are often constrained by elevated viscosity, largely driven by protein–protein interactions, which complicates manufacturing and limits subcutaneous delivery. Early viscosity risk assessment is essential during discovery, yet traditional measurements require large sample volumes, and lack high-throughput capability. Here, we develop a high-throughput small-angle X-ray scattering (SAXS) protocol to detect mAb self-association at dilute concentrations, enabling early predictive insights into high-concentration viscosity. Synchrotron SAXS measurements were conducted for 21 mAbs formulated in a histidine buffer at pH 6.0. An initial subset of 10 mAbs analyzed across 1–150 mg/mL revealed that effective structure factor transitions in the low-q region, indicative of interparticle interactions, consistently emerged below 25 mg/mL. Subsequently, 11 additional mAbs were analyzed at 1–25 mg/mL using automated liquid handling and flow cells to enable high-throughput screening. High-viscosity mAbs exhibited detectable low-q upturns at concentrations ≤10 mg/mL, whereas low-viscosity mAbs showed downturns. A classification criterion based on effective structure factor transitions accurately classified all high- and low-viscosity mAbs at 150 mg/mL, offering a scalable, sample-efficient alternative to conventional methods. These results extend recent findings on the concentration-dependent sensitivity of SAXS to short-range attractions, demonstrating that they can emerge at lower concentrations than previously reported. This study presents the most comprehensive and diverse SAXS dataset for mAbs reported to date within a single formulation, providing a valuable resource for developing and validating coarse-grained models that can more accurately capture intermolecular interactions governing high-concentration solution behavior, thereby enabling rational antibody engineering and improved developability.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Bulk viscosity of nuclear matter with pions in the neutrino-trapped regime

Recent work has shown that a population of thermal pions could modify the equation of state and transport properties of hot and dense neutron-rich matter and introduce new reaction pathways to change the proton fraction. In this article we study their impact on the bulk viscosity of dense matter, focusing on the neutrino-trapped regime that would be realized in neutron star mergers and supernovae. We find that the presence of a thermal population of pions alters the bulk viscosity by modifying the EoS (via the susceptibilities) and by providing new reaction pathways to achieve beta equilibrium. In neutron star merger conditions, the bulk viscosity in neutrino-trapped npeμ matter (without pions) has its peak at temperatures of at most a couple MeV and is quite small at temperatures of tens of MeV. We find that thermal pions enhance the low-temperature peak of the bulk viscosity by a factor of a few and shift it to slightly lower temperatures. At higher temperatures, where the pion abundance is large but the bulk viscosity is traditionally small, pions can increase the bulk viscosity by an order of magnitude or more, although it is still orders of magnitude smaller than its peak value.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Evolution of eddy viscosity in the wake of a wind turbine

The eddy viscosity hypothesis is a popular method in wind turbine wake modeling for estimating turbulent Reynolds stresses. We document the downstream evolution of eddy viscosity in the wake of a wind turbine from experimental and large-eddy-simulation data. Wake eddy viscosity is isolated from its surroundings by subtracting the inflow profile, and the driving forces are identified in each wake region. Eddy viscosity varies in response to changes in turbine geometry and nacelle misalignment with larger turbines generating stronger velocity gradients and shear stresses. We propose a model for eddy viscosity based on a Rayleigh distribution. Model parameters are obtained from scaling the eddy viscosity hypothesis and demonstrate satisfactory agreement with the reference data. The model is implemented in the curled wake formulation in the FLOw Redirection and Induction in Steady State (FLORIS) framework and assessed through comparisons with the previous formulation. Our approach produced more accurate flow field estimates with lower total error for the majority of cases.

17 WIND ENERGY↗
DOE OSTI2022

Optimization of matrix viscosity improves polypropylene fiber dispersion and properties of engineered cementitious composites

Engineered cementitious composites (ECC) is a durable cementitious material with high tensile ductility and strain-hardening characteristics. Although considered as a cost-effective fiber for ECC, polypropylene (PP) fiber is reportedly difficult to disperse in mortar matrix due to its high aspect ratio and hydrophobicity. In this study, the matrix viscosity of a low-carbon ECC based on limestone calcined clay cement was tailored as a variable to improve PP fiber dispersion under a pre-determined mixing protocol. The effect of matrix viscosity on the composite fresh and hardened properties was investigated experimentally. Results suggested an optimal range of matrix viscosity (10.3–11.5 Pa . s) favors the composite tensile strength and strain capacity at 28 days. At the optimal state with a 0.1 % viscosity modifying admixture (VMA)-to-binder mass ratio, PP-ECC achieved 7.0 % tensile strain capacity and 3.5 MPa ultimate tensile strength. When matrix viscosity falls outside the desired range, both ultimate tensile strength and strain capacity were diminished. By tailoring the VMA dosage, the matrix viscosity can be adjusted for desired fiber dispersion, workability, and mechanical properties. Finally, the findings of this study provide a technical reference for the practical design and application of PP-ECC.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Impact of melt viscosity on filler dispersion in elastomeric nanocomposites

Compounding of commercial nanocomposites usually involves the addition of viscosity enhancers such as binder resins in ink jet inks, and paints. Contrary to this, plasticizers such as process oils are added to reduce the melt viscosity and ease processability of reinforced elastomers. Nanofillers such as silica and carbon black are typically added to reinforce rubber and enhance performance of automotive tire treads. Filler dispersion has traditionally been qualitatively (indirectly) assessed by measuring the properties of reinforced elastomers. While dispersion can be quantified by examining filler agglomeration through surface roughness measurements and microscopy, the size-scale dependence for these hierarchical fillers has usually been ignored. Here, we have recently devised a method to quantify nano-scale dispersion of fillers using Ultra small-angle X-ray scattering (USAXS) techniques. This method is advantageous since it directly links the controllable processing/compounding parameters such as the mixing speed, mixer geometry, residence time (or mixing duration), melt density, flow gap distance, and melt viscosity to nano-scale dispersion. While our previous studies have explored the impact of different processing parameters, this study specifically investigates the impact of melt viscosity on nano-scale filler dispersion in elastomer compounds. Commercially available polybutadienes with different Mooney viscosities were used in conjunction with different grades and amounts of process oils to modify the melt viscosity.

Carbon black↗
DOE OSTI2020

Hall viscosity of the composite-fermion Fermi seas for fermions and bosons

The Hall viscosity has been proposed as a topological property of incompressible fractional quantum Hall states and can be evaluated as Berry curvature. This paper reports on the Hall viscosities of composite-fermion Fermi seas at ν = 1/m, where m is even for fermions and odd for bosons. A well-defined value for the Hall viscosity is not obtained by viewing the 1/m composite-fermion Fermi seas as the n → ∞ limit of the Jain ν = n/(nm ± 1) states, whose Hall viscosities (±n + m)ℏρ/4 (ρ is the two-dimensional density) approach ±∞ in the limit n → ∞. A direct calculation shows that the Hall viscosities of the composite-fermion Fermi sea states are finite and also relatively stable with system size variation, although they are not topologically quantized in the entire τ space. Here, I find that the ν = 1/2 composite-fermion Fermi sea wave function for a square torus yields a Hall viscosity that is expected from particle-hole symmetry and is also consistent with the orbital spin of 1/2 for Dirac composite fermions. I compare my numerical results with some theoretical conjectures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Bulk viscosity from Urca processes: n p e μ matter in the neutrino-transparent regime

We study the bulk viscosity of moderately hot and dense, neutrino-transparent relativistic npeμ matter arising from weak-interaction direct Urca processes. This work parallels our recent study of the bulk viscosity of npeμ matter with a trapped neutrino component. The nuclear matter is modeled in a relativistic density functional approach with two different parametrizations—DDME2 (which does not allow for the low-temperature direct-Urca process at any density) and NL3 (which allows for low-temperature direct-Urca process above a low-density threshold). Here, we compute the equilibration rates of Urca processes of neutron decay and lepton capture, as well as the rate of the muon decay, and find that the muon decay process is subdominant to the Urca processes at temperatures T ≥ 3 MeV in the case of DDME2 model and T ≥ 1 MeV in the case of NL3 model. Thus, the Urca-process-driven bulk viscosity is computed with the assumption that pure leptonic reactions are frozen. As a result the electronic and muonic Urca channels contribute to the bulk viscosity independently and at certain densities the bulk viscosity of npeμ matter shows a double-peak structure as a function of temperature instead of the standard one-peak (resonant) form. In the final step, we estimate the damping timescales of density oscillations by the bulk viscosity. We find that, e.g., at a typical oscillation frequency f = 1~kHz, the damping of oscillation is most efficient at temperatures 3 ≤ T ≤ 5~MeV and densities n B ≤ 2n 0 where they can affect the evolution of the post-merger object.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Viscosity of the magnetized strongly coupled one-component plasma

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

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

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

Krisher, James A.↗
DOE OSTI2023

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

As presented in this report, thermal conductivity and viscosity measurements were performed on key chloride pseudo-binary molten salt systems of relevance to molten salt reactor developers. Thermal conductivity measurements were conducted with a variable gap technique, in which a known heat flux is driven across a molten salt specimen and the temperature difference is measured, allowing calculation of the thermal conductivity. This is achieved by establishing a small gap between the bottom of a cylindrical inner containment, which houses electrical heating elements, and an outer containment, which houses cooling channels; the gap size can be varied by compression of a formed bellows. A new calibration scheme was developed herein, involving a correction factor to the heat flux based on He measurements at various temperatures. Furthermore, the data processing methodology was improved to minimize the impact of radiative heat transfer in the calculation of salt specimen thermal conductivity from the temperature difference measurements. Viscosity measurements were conducted with a rolling ball viscometer, in which a ball rolls some known distance in an angled tubular crucible, and the terminal velocity can be used to calculate the viscosity of the salt. The measurement can be performed in a quartz crucible, with which a standard camera can be used to track the ball, or in a metal crucible, with which x-ray radiography is required to track the ball. A new custom x-ray system was made and dedicated to the rolling ball viscometer to enable high throughput automated measurements with salts which require containment with metal. Both systems have been integrated with a ventilation stack which allows for off-gassing of radioactive material, enabling future measurements with U-bearing salts. The thermal conductivity measurements performed herein were with NaCl-KCl (44 mol% NaCl). This salt system was measured in the previous fiscal year, however the thermal conductivity values obtained were comparatively low, and so the system was remeasured with the aforementioned calibration scheme and improved post-processing techniques. The newly obtained thermal conductivity values for NaCl-KCl (44 mol% NaCl) indicate good agreement with kinetic theory and ab-initio models (within 5–10 %). The viscosity measurements performed herein were with three different compositions of NaCl-KCl: 75, 50, and 25 mol% NaCl; a new calibration scheme was applied to account for variable flow effects in the laminar regime. The results show reasonable agreement with literature (5–20 %, depending on the temperature, composition, and study); however, literature values are likely higher than true values based on pure end-member measurements performed in the comparative studies. The measurements conducted herein do show a trend such that viscosity increases with increasing NaCl concentration, which agrees well with one of two comparative studies. In general, the measurements conducted herein gives confidence in the capability to use these systems to accurately measure thermal conductivity and viscosity of actinide-bearing salts within the next fiscal year.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

A theoretical/computational framework to measure SiO2 and MgO viscosity at high pressure

The convection of the mantle of Earth and super-Earths is important for many terrestrial phenomena, from plate tectonics to outgassing. Rheological properties, such as viscosity, regulate the transport of thermal energy and mass. However, the viscosity of mantle-relevant materials, such as MgO, at relevant pressures (>120 GPa) are not well constrained. The objective of this work was to develop a computational platform to simulate novel experiments aiming to measure the viscosity of MgO at high pressures. Experiments performed by our collaborators at Johns Hopkins University and Lawrence Livermore National Laboratory use the OMEGA EP laser facility to shock a corrugated MgO interface to 170 GPa, with resulting velocity evolution governed by the viscous Richtmyer-Meshkov instability. We used an in-house hydrocode to simulate this process and thus provide a bound on viscosity by comparing our simulations results to these experiments, as well as to examine the physical processes at play. We simulated the experiments with different values of MgO viscosity (from inviscid to 10,000 Pa∙s) while taking into account the unsteady laser pulse, the material's equation of state, and the rate-dependent constitutive relation for MgO, and the materials' equations of state. Our results suggest that MgO at these conditions has a viscosity within an order of magnitude of 5000 Pa∙s.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A spline-based method to obtain spatially dependent viscosity in confined flows

Coupling chemical physics to continuum theories is a critical step to understanding multi-scale phenomena. This paper will connect non-equilibrium molecular dynamics simulations to a continuum-based Navier-Stokes equation that has relaxed the assumption of spatial uniformity in viscosity. Using a form for viscosity based on spline interpolation, viscosity as a function of position is obtained from the least squares fit of the velocity profile measured from molecular simulations of flow in a nanochannel. Viscosity can vary widely, particularly near the channel boundaries, indicating that uniform viscosity is no longer appropriate. Variations of the viscosity near the channel surfaces imply that considering solution and surface chemistry could be necessary to rigorously understand molecular-scale flows in nanochannels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗