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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Computational Fluid Dynamics Simulations of Glass Vitrification Refractory Coupon Tests

The Waste Treatment and Immobilization Plant (WTP) at the Hanford site is nearing the start of the Direct-Feed Low-Activity Waste (DFLAW) operations. DFLAW is destined to convert a pretreated low activity waste portion of the 56 million gallons of tank waste into a stable solid glass. In the subsequent decade completion of the high-level waste (HLW) facility is anticipated. Sustained operational missions of both LAW and HLW melter facilities are expected over multiple decades. In high-temperature glass melters, the refractory lining corrodes over time, which could potentially be an issue for longer term operations, this refractory corrosion is higher at the level of the glass-air interface due to surface tension driven flow. The glass viscosity, melt pool temperature, and glass chemical composition can impact the rate at which the refractory corrodes. This rate is important to quantify for the various waste glasses to be produced at the WTP since the integrity of the refractory should not be a limiting factor affecting the lifetime of the melter. To this end, a series of glasses representative of the first batches of waste glass produced by the WTP will be melted in small-scale crucibles with Monofrax® K-3 coupons inserted. The corrosion of the K-3 will be measured in the melt and at the meltline (or neckline). A model for the corrosion rate will be constructed and implemented into a previously developed framework for a computational fluid dynamics (CFD) model of the full-scale WTP. To assist with experimental design and validate the implementation of the model in the full-scale melter, CFD simulations of the small-scale crucible tests were performed. The bubbling that occurs in the small-scale crucible is initially validated here with a model that uses silicone oil at room temperature. The viscosity of the oil ranges from 1 to 100 Pa•s, which corresponds to operating glass pool temperatures near 1150 °C down to idling temperatures near 950 °C. The simulation results show good agreement with the bubble sizes that form during experiments. CFD modeling of the crucible setup was used to determine bubbling characteristics to match the range of near-wall velocities expected in the full-scale WTP. This study presents the initial CFD modeling results, corrosion testing plan, and some preliminary corrosion samples with an outline for the next steps for the development of the corrosion model.

Abboud, Alexander W. [Idaho National Lab]↗
DOE OSTI2026

Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification

The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2024

Numerical investigation of multiphase flow through self-affine rough fractures

Multiphase flow through fractures has great significance in subsurface energy recovery and gas storage applications. Different fracture and flow properties affect flow through a fracture which is difficult to control in laboratory experiments. Here, we perform lattice Boltzmann simulations in an ensemble of synthetically generated fractures. Drainage simulations are performed at different capillary numbers, wettability, and viscosity ratios. We track the invading front and quantify breakthrough saturations and show that roughness and wettability have a strong effect on fluid invasion through a complex fracture. Invading a more viscous fluid results in more stable displacement regardless of the capillary number while at very low capillary numbers, fluid migration is dependent on the inherent structure of the fracture. We develop a fluid displacement phase diagram in a single rough fracture and compare our results from that in the literature. Finally, we extend the phase diagrams across multiple fractures and demonstrate the importance of natural fracture features of roughness and wettability in identifying stable versus unstable displacement regimes during multiphase flow through rough fractures. Our work presents an end-to-end numerical pathway for testing on experimental data and expanding numerical data sets for testing combinations of different physical phenomenon and make valuable predictions on fluid flow through rough fractures.

02 PETROLEUM↗
DOE OSTI2024

Glass production rate in an electric melter: Melting rate correlation and primary foam stability

Waste loading and melting rate are major factors determining the lifecycle of nuclear waste disposal by vitrification in electric melters in which the heat to the cold cap is delivered from the melt pool. One of the crucial parameters that affect the rate of melting is the temperature at which the primary foam collapses at the cold cap bottom. Further, apart from the transient glass-forming melt viscosity, primary foam stability is affected by the shear stress imposed by the flow in the melt pool and the presence of residual solid particles.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The effect of hydroxyl spacing in diols on the solvation structure, dynamics, and transport properties of choline chloride-based deep eutectic solvents

Deep eutectic solvents (DESs) are a class of liquids that offer great potential in alleviating some of the challenges present in today's long-term energy storage methods because they have physical properties that are favorable for storable electrolyte solutions. In this work, a series of glycols (ethylene glycol, 1,3-propanediol, 1,4-butanediol, and 1,5-pentanediol) were studied as potential hydrogen bond donors (HBD) with a common choline chloride (ChCl) as the hydrogen bond acceptor (HBA). The solvation dynamics of the prepared systems were studied by measuring the solvent reorganization response using femtosecond transient absorption spectroscopy (fs-TA). Conductivity, viscosity, density, ET(30) polarity, and dynamics of the prepared DESs were analyzed, with a particular interest in determining the effect of HBD chain length on these parameters. Here, classical molecular dynamics simulations were employed to investigate how the local liquid structure, solvent dynamics, and bulk solvent properties vary with changes in glycol chain length.

Conductivity↗
DOE OSTI2026

Beam-energy dependence of correlations between mean transverse momentum and anisotropic flow of charged particles in Au+Au collisions at RHIC

The correlation between the mean transverse momentum, [ p T ], and the squared anisotropic flow, $v^2_n$, on an event-by-event basis has been suggested to be influenced by the initial conditions in heavy-ion collisions. We present measurements of the variances and covariance of [ p T ] and $v^2_n$, along with their dimensionless ratio, for Au+Au collisions at various beam energies: $\sqrt{s_{NN}}$ = 14.6, 19.6, 27, 54.4, and 200 GeV. Our measurements reveal a distinct energy-dependent behavior in the variances and covariances. In addition, the dimensionless ratio displays a similar behavior across different beam energies. We compare our measurements with hydrodynamic models and similar measurements from Pb+Pb collisions at the Large Hadron Collider (LHC). These findings provide valuable insights into the beam energy dependence of the specific shear viscosity ( η/s ) and initial-state effects, allowing for differentiating between different initial-state models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Temperature-Dependent Molecular Diffusional Properties in Deep Eutectic Solvents and Eutectogels

Eutectogels (ETGs) prepared from deep eutectic solvents (DESs), a gelator, and water have many uses in separations, catalysis, and energy storage systems. In these applications, temperature-dependent molecular diffusional properties and intermolecular interactions play a critical role in their function. Diffusional properties of Alexa Fluor 633 and ATTO 647N were measured across a range of temperatures in choline chloride:2glycerol DESs comprised of one molar equivalent of choline chloride and two molar equivalents of glycerol (also known as glyceline) as well as ETGs made from this DES, a xanthan gum gelator, and 10% w/w water or 20% w/w water. Fluorescence recovery after photobleaching (FRAP) was employed to evaluate potential changes in molecular diffusion within DESs and ETGs from 20 to 100 °C. Surprisingly, the ETGs have a larger sample viscosity but also exhibit faster molecular diffusion when compared to the dry DES. This is attributed to macroscopic properties of the ETGs (e.g., pores, three-dimensional structure). The FRAP data also show an irreversible temperature-dependent decrease in the diffusion coefficient of both fluorophores in the ETGs. This is consistent with differential scanning calorimetry data of the ETG, which shows a shift to a more negative glass transition temperature (−36 °C to −47 °C) after the first cooling/heating cycle. Raman data reveal no detectable changes in the intermolecular interactions in either the DESs or ETGs as a function of temperature or water content. Finally, the findings on temperature- and water-dependent diffusional properties of DESs and ETGs provide a foundation for optimizing the use of these materials across various applications, particularly where repeated heating/cooling cycles may be used.

alcohols↗
DOE OSTI2025

Automation-Accelerated Electrolyte Design Mitigates Solubility Competition between Redox-Active Molecules and Supporting Salts

In nonaqueous redox-flow batteries (NRFBs), redox-active organic molecules (ROMs) and supporting salts compete for solvation sites, limiting achievable energy density. We combine automated high-throughput experimentation (HTE) with camera-based saturation monitoring and quantitative NMR to measure paired (ROM, salt) solubilities across single and mixed organic solvents. Using 2,1,3-benzothiadiazole (BTZ) with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) as a model system, we find that a binary m-xylene/acetonitrile mixture dissolves ≈3 M of both BTZ and LiTFSI─surpassing the previously reported 2 M ceiling for neat acetonitrile─by leveraging complementary solvation (MX is BTZ-philic and salt-phobic; ACN stabilizes LiTFSI). A random-forest model (RMSE ≈ 0.24) trained on solvent descriptors highlights log P and salt concentration as dominant predictors and predicts MX/ACN ≈0.3/0.7 (v/v) to be near-optimal. These formulations retain practical viscosity and ∼5 mS·cm –1 conductivity at high loading. In conclusion, the workflow provides a reproducible, data-centric route to NRFB electrolyte design and motivates an open, standardized dual-solute solubility resource for accelerated electrolyte discovery.

Electrolytes↗
DOE OSTI2026

How to Minimize Faradaic Efficiency Error in Electrochemical CO 2 Reduction for Gas Products

Faradaic efficiency (FE) is an important metric for evaluating electrochemical processes, such as carbon dioxide reduction, that is often used for performance comparisons. Although the equation to calculate FE is well-known, details needed to yield accurate values are often overlooked, potentially leading to errors exceeding 100%. Avoiding any errors is crucial for drawing correct conclusions, especially during a catalysis optimization process. Factors with high potential for FE errors include incorrect mass flow rates and temperature values, imprecise gas chromatography calibrations, and uncorrected gas viscosities, while factors with less potential are also identified. In conclusion, this manuscript presents FE calculation guidelines intended for both newcomers and experts in the field of electrochemical fuel production.

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

Secondary Organic Aerosol from OH Oxidation of Acyclic Terpenes Is More Viscous and Less Volatile than That of Their Cyclic Analogs

Biogenic volatile organic compounds (BVOCs), a dominant source of secondary organic aerosol (SOA) globally, exhibit emission rates and composition that are plant species-specific and vary with environmental stressors. A common outcome of plant stress is increased emissions of acyclic terpenes. The paucity of information about acyclic terpene SOA chemistry contributes to uncertainties in predictions of SOA global loadings and impacts on Earth’s radiative budget, particularly in a changing climate where acyclic terpene emissions could become more prominent. This study compared properties of SOA derived from OH oxidation of acyclic and cyclic monoterpenes (ß-ocimene, a-pinene) and sesquiterpenes (ß-farnesene, ß-caryophyllene). Single particle mass spectrometry was used for assessing shape, density, and evaporation kinetics of size-selected SOA particles, and nanospray desorption electrospray ionization high-resolution mass spectrometry (nano-DESI-HRMS) was used to measure molecular composition of SOA. Acyclic terpene SOA exhibited higher viscosity and lower volatility compared to cyclic terpene SOA, and had a greater volume fraction remaining (VFR) after ~24 hours of evaporation - approximately 1.3-1.6 times higher VFR than that of cyclic terpene SOA. Additionally, HRMS analysis revealed greater chemical diversity and higher fractions of extremely low volatility compounds (56-62% ELVOC/LVOC) in acyclic terpene SOA compared to cyclic counterparts (25-37% ELVOC/LVOC). Our findings highlight the potential importance of accounting for acyclic terpene aerosol chemistry under conditions of plant stress to improve predictions of SOA loadings and impacts.

SOA physical properties↗
DOE OSTI2025

Prolonging All-Optical Molecular Electron Spin Coherence in the Tissue Transparency Window

Coherent electron spin states within paramagnetic molecules hold significant potential for microscopic quantum sensing. However, all-optical coherence measurements amenable to high spatial and temporal resolution under ambient conditions remain a significant challenge. Here we conduct room-temperature, picosecond time-resolved Faraday ellipticity/rotation (TRFE/R) measurements of the electron spin decoherence time T 2 * in [IrBr 6 ] 2- . Decoherence is strongly sensitive to solution phase viscosity, pointing to molecular tumbling as an important decoherence mechanism. Accordingly, immobilization of [IrBr 6 ] 2- molecules in thin polymer films results in an order-of-magnitude increase in coherence lifetime and significantly greater magnetic field sensitivity. Here, by tuning energies of ligand-to-metal charge transfer (LMCT) states, TRFE/R enables spin initialization and readout in the tissue transparency window, paving the way toward all-optical, ultrafast molecular electron spin coherence imaging in biological systems.

Electron paramagnetic resonance spectroscopy↗
DOE OSTI2024

Review of the second charged-particle transport coefficient code comparison workshop

We report the results of the second charged-particle transport coefficient code comparison workshop, which was held in Livermore, California on 24–27 July 2023. This workshop gathered theoretical, computational, and experimental scientists to assess the state of computational and experimental techniques for understanding charged-particle transport coefficients relevant to high-energy-density plasma science. Data for electronic and ionic transport coefficients, namely, the direct current electrical conductivity, electron thermal conductivity, ion shear viscosity, and ion thermal conductivity were computed and compared for multiple plasma conditions. Additional comparisons were carried out for electron–ion properties such as the electron–ion equilibration time and alpha particle stopping power. Overall, 39 participants submitted calculated results from 18 independent approaches, spanning methods from parameterized semi-empirical models to time-dependent density functional theory. In the cases studied here, we find significant differences—several orders of magnitude—between approaches, particularly at lower temperatures, and smaller differences—roughly a factor of five—among first-principles models. We investigate the origins of these differences through comparisons of underlying predictions of ionic and electronic structure. The results of this workshop help to identify plasma conditions where computationally inexpensive approaches are accurate, where computationally expensive models are required, and where experimental measurements will have high impact.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Collision integrals within the Chapman–Enskog theory for a generalized Lennard-Jones potential

Here we report the values of the collision integrals, needed for the calculation of the macroscopic transport properties such as viscosity (η) and diffusion coefficient (D) of gases within the Chapman–Enskog kinetic gas theory, for a generalized Lennard-Jones potential (gLJ), a more general potential with an adjustable long range 1/r dependence that can describe a wide range of intermolecular interactions.

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

Molecular dynamics simulation of hydrodynamic transport coefficients in plasmas

Molecular dynamics simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green–Kubo formalism over a broad range of Coulomb coupling strength, 0.01 ≤ Γ ≤ 140. Emphasis is placed on testing standard results of the Chapman–Enskog solution in the weakly coupled regime (Γ ≪ 1) using these first-principles simulations. As expected, the results show good agreement for Γ ≲ 0.1. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green–Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal, and thermoelectric coefficients computed in the Chapman–Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime (Γ ≫ 1). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.

Electrical conductivity↗
DOE OSTI2026

Plasma hydrodynamics from mean force kinetic theory

Mean force kinetic theory is used to evaluate the electrical conductivity, thermal conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity of a two-component plasma. The results are compared with molecular dynamics simulations. These simulations are made possible by assuming a repulsive Coulomb force for all interactions. Good agreement is found for all coefficients up to a Coulomb coupling strength of Γ ≈ 20. This is over 100-times larger than the coupling strength at which traditional theories break down. It is concluded that mean force kinetic theory is a promising framework for extending hydrodynamics to dense plasmas.

Electrical conductivity↗
DOE OSTI2025

Quantum Ornstein-Zernike theory for two-temperature two-component plasmas

Laboratory plasma production almost always preferentially heats either the ions or electrons, leading to a two-temperature state. In this state, density functional theory molecular dynamic simulation is the state of the art for modeling bulk material properties. We construct a statistical mechanics model for the two temperature limit that is theoretically consistent with the molecular dynamics method. We proceed to derive the electron-ion multi-temperature quantum Ornstein-Zernike equations for the first time. This allows the construction of a two-temperature two-component plasma model using the average atom from which we can compute bulk material properties at a fraction of the computation time of the two-temperature density functional theory simulation. The accuracy of the model is benchmarked against ion pair correlation and self-diffusion results from ab initio simulation. Here, we proceed to compute the viscosity and ion thermal conductivity as a function of both ion and electron temperature.

Ab initio molecular dynamics↗
DOE OSTI2024

Thermoplastic Elastomers and Their Physical Gels Electrospun into Tunable Microfibrous Nonwoven Mats: Structure Formation and Property Enhancement

Abstract Thermoplastic elastomers (TPEs) based on styrenic block copolymers constitute excellent examples of self‐networking macromolecules that are employed in a wide range of contemporary technologies as molded parts. In such applications, these TPEs exist as dense (nonporous) films or other shapes. Here, it is first demonstrated that a series of commercial TPEs possessing comparable compositions can be electrospun from solution to form microfibers that are arranged into nonwoven mats that are breathable. An important consideration for microfiber formation is the copolymer molecular weight, which regulates i) the viscosity of the parent solution prior to electrospinning, ii) the ability of these copolymers to self‐assemble during electrospinning, iii) the microfiber morphology, and iv) the mechanical properties of the resultant microfibers. The addition of a midblock‐selective aliphatic oil to these TPEs yields thermoplastic elastomer gels (TPEGs), wherein the copolymer morphology and mechanical properties become highly composition‐tunable. Electrospinning TPEGs from a binary oil+solvent solution introduces a micelle inversion mechanism that begins with an oil‐rich micellar core and ends with a styrene‐rich micellar core, required for network stabilization, as the solvent dries during microfiber solidification. This work has implications for the production of controllably low‐modulus microfibrous materials possessing modestly improved toughness but exceptional extensibility and enhanced optical transparency.

Shamsi, Mohammad↗
DOE OSTI2024

Additive Manufacturing with Cellulose-Based Composites: Materials, Modeling, and Applications

Recent advances in large-scale additive manufacturing (AM) with polymer-based composites have enabled efficient production of high-performance materials. Cellulose nanomaterials (CNMs) have emerged as bio-based feedstocks due to their exceptional strength and sustainability. However, challenges such as hornification and poor dispersion in polymer matrices still limit large-scale CNM–polymer composite manufacturing, requiring novel strategies. Here, this review outlines an approach starting with atomic-level simulations to link molecular composition to key parameters like bulk density, viscosity, and modulus. These simulations provide data for finite element analysis (FEA), which informs large-scale experiments and reduces the need for extensive trials. The strategy explores how atomic interactions impact the morphology, adhesion, and mechanical properties of CNM-based composites in AM processes. The review also discusses current developments in AM, along with predictions of mechanical and thermal properties for structural applications, packaging, flexible electronics, and hydrogel scaffolds. By integrating experimental findings with molecular dynamics (MD) simulations and finite element modeling (FEM), valuable insights for material design, process optimization, and performance enhancement in CNM-based AM are provided to address ongoing challenges.

36 MATERIALS SCIENCE↗