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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Hydration and transport properties of cesium hydroxide and mixed cesium hydroxide–sodium nitrite aqueous solutions

Here, this study explores the hydration and transport properties of aqueous cesium hydroxide (CsOH) solution, with or without 1 molar (M) sodium nitrite (NaNO 2 ). Historic studies of electrolyte solutions indicate that Cs + ions decrease viscosity and increase diffusion rates, whereas OH − ions have the opposite effect. Here, the influence of OH − was dominant in CsOH solutions, leading to increased viscosity and reduced diffusion rates. There was a linear relationship between diffusion coefficients and water activity, emphasizing the significant role of ion–water interactions in determining transport properties. This may be because the interaction between Cs + and the anions is weak even when they are in direct contact with each other. The weakness of the ion-pairing was established through thermodynamic analysis. The findings suggest that ion-pairing is not the only important interaction controlling transport properties when ion-pairing is weak. Nonetheless, ion-pairing or obstructions did result in more sluggish transport properties as electrolyte concentrations increased. Overall, the research enhances the understanding of the complexities underlying ion interactions in multicomponent solutions.

Reynolds, Jacob G. [Hanford Site (HNF), Richland, ↗
DOE OSTI2024

Quantifying electron transport in aggregated colloidal suspensions in the strong flow regime

Electron transport in complex fluids, biology, and soft matter is a valuable characteristic in processes ranging from redox reactions to electrochemical energy storage. These processes often employ conductor–insulator composites in which electron transport properties are fundamentally linked to the microstructure and dynamics of the conductive phase. While microstructure and dynamics are well recognized as key determinants of the electrical properties, a unified description of their effect has yet to be determined, especially under flowing conditions. In this work, the conductivity and shear viscosity are measured for conductive colloidal suspensions to build a unified description by exploiting both recent quantification of the effect of flow-induced dynamics on electron transport and well-established relationships between electrical properties, microstructure, and flow. These model suspensions consist of conductive carbon black (CB) particles dispersed in fluids of varying viscosities and dielectric constants. In a stable, well-characterized shear rate regime where all suspensions undergo self-similar agglomerate breakup, competing relationships between conductivity and shear rate were observed. To account for the role of variable agglomerate size, equivalent microstructural states were identified using a dimensionless fluid Mason number, Mn f , which allowed for isolation of the role of dynamics on the flow-induced electron transport rate. At equivalent microstructural states, shear-enhanced particle–particle collisions are found to dominate the electron transport rate. This work rationalizes seemingly contradictory experimental observations in literature concerning the shear-dependent electrical properties of CB suspensions and can be extended to other flowing composite systems.

Science & Technology - Other Topics↗
DOE OSTI2025

High-Fidelity CFD Modeling of Cryogenic Hydrogen Isotope Extrusion for Fusion Reactor Pellet Fueling

This study investigates the extrusion processes of deuterium and protium using ANSYS-Polyflow. The geometries and computational fluid dynamics (CFD) settings closely replicate the experimental setups and data acquired from the extruder experiments at Oak Ridge National Laboratory (ORNL) for validation purposes. We explore the impacts of (1) slip versus non-slip boundary conditions and (2) the use of constant, temperature-, and shear rate–dependent viscosities, concluding that the implementation of non-slip wall boundary conditions combined with shear rate–dependent viscosity produced more accurate predictions. The simulations achieved excellent agreement with the experimental data, with relative differences of only 5% for deuterium, and 3% to 6% for protium. This is the first time that experimental extrusion data at ORNL have been accurately predicted through high-fidelity CFD modeling. In conclusion, the advancements offer valuable insights and a foundational modeling tool for optimizing pellet injectors for ITER and other future reactor-scale devices.

ANSYS-Polyflow↗
DOE OSTI2024

Simulation of Li-pellet triggered ELMs in EAST with an impurity model implemented under BOUT++ framework

A simple impurity model has been developed under the BOUT++ framework to investigate the Li-pellet triggered edge localized mode (ELM) in EAST configuration. The present impurity model decouples the ion pressure enhancement effect (IPEE) (which affects radial E X B flow shear and gyro-viscosity of deuterium ions and is implicitly included in our previous work) to make the simulation of Li pellet injection more physically reasonable. In addition, the present impurity model also includes the impurity equilibrium effect (IEE) (which induces modifications on vorticity and gyro-viscosity). The simulation results show that without IPEE a turbulent ELM induced by Li pellet is triggered by multiple peeling-ballooning modes (PBMs) rather than by a single dominant mode, and the simultaneous growth of multiple modes is conducive to reducing the time for pedestal entering the energy loss state. When a turbulent ELM occurs, the nonlinearly dominant modes undergo a secondary fast growth during fast crash phase. The finding explains the secondary increase of D α emission observed in the DIII-D Li injection experiments. During the evolution of a turbulent ELM, the continuous pedestal collapse dominates the pedestal energy loss in turbulent transport phase; however, in saturation phase, the stabilizing effects on PBMs (n ≠ 0) by n = 0 mode and radial E X B flow shear within the steepest gradient region are more prominent. It is also found that with IEE and without IPEE, the size threshold of ELM triggering and the magnitude of ELM sizes both show a good agreement with EAST and DIII-D experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Impurity radiation seeding of neoclassical tearing mode growth

Abstract The physics of neoclassical tearing modes (NTMs) is of great concern to tokamak plasma stability and performance, especially in the burning plasma regime. Whereas in many situations the seed events can be clearly identified, such as sawteeth and edge localized modes, the potential seeding mechanism of NTMs due to the resistive tearing instability driven by the impurity radiation cooling still needs more study. Recent NIMROD simulations have demonstrated that local impurity radiation cooling can drive the seed island growth and trigger the subsequent onset of NTM instability. The seed island is mainly driven by the local helical perturbation of the diamagnetic current induced by the perturbed pressure gradient as a result of the impurity radiative cooling on the rational surface. A heuristic closure for the neoclassical viscosity is adopted, and the seed island is further driven by the perturbed bootstrap current induced from the neoclassical electron viscous stress in the extended Ohm’s law. The growth rate of the NTM in simulations is found proportional to the electron neoclassical viscosity, and a theoretical neoclassical driving term is adopted to account for the nonlinear neoclassical island growth in the simulations.

Zeng, Shiyong (ORCID:000000017098818X)↗
DOE OSTI2024

Droplet bioprinting of acellular and cell-laden structures at high-resolutions

Advances in digital light projection(DLP) based (bio) printers have made printing of intricate structures at high resolution possible using a wide range of photosensitive bioinks. A typical setup of a DLP bioprinter includes a vat or reservoir filled with liquid bioink, which presents challenges in terms of cost associated with bioink synthesis, high waste, and gravity-induced cell settling, contaminations, or variation in bioink viscosity during the printing process. Here, we report a vat-free, low-volume, waste-free droplet bioprinting method capable of rapidly printing 3D soft structures at high resolution using model bioinks and model cells. A multiphase many-body dissipative particle dynamics model was developed to simulate the dynamic process of droplet-based DLP printing and elucidate the roles of surface wettability and bioink viscosity. Process variables such as light intensity, photo-initiator concentration, and bioink formulations were optimized to print 3D soft structures (∼0.4–3 kPa) with a typical layer thickness of 50 µm, an XY resolution of 38 ± 1.5 μm and Z resolution of 237 ± 5.4 µm. To demonstrate its versatility, droplet bioprinting was used to print a range of acellular 3D structures such as a lattice cube, a Mayan pyramid, a heart-shaped structure, and a microfluidic chip with endothelialized channels. Droplet bioprinting, performed using model C3H/10T1/2 cells, exhibited high viability (90%) and cell spreading. Additionally, microfluidic devices with internal channel networks lined with endothelial cells showed robust monolayer formation while osteoblast-laden constructs showed mineral deposition upon osteogenic induction. Overall, droplet bioprinting could be a low-cost, no-waste, easy-to-use, method to make customized bioprinted constructs for a range of biomedical applications.

DLP↗
DOE OSTI2026

Radial oscillations of viscous neutron stars: Zero diffusion case

In this article, the spectrum of radial oscillations of neutron stars is systematically studied within two frameworks of viscous relativistic hydrodynamics: the relativistic Navier-Stokes and Israel-Stewart theories. A correspondence is established between the discrete stellar eigenmodes and the continuous dispersion relation of perturbations around a homogeneous fluid, providing a basis for interpreting our numerical results. We analyze the Newtonian limit and assess the impact of relativistic corrections, such as the gravitational redshifting of microscopic relaxation timescales. We show that bulk viscosity can significantly affect the behavior of both hydrodynamic and nonhydrodynamic fundamental modes, and that, depending on the magnitude of the viscous effects, it is the nonhydrodynamic mode that becomes unstable beyond the turning point in a sequence of equilibrium configurations. These results provide a useful step toward systematic studies of neutron star quasinormal modes in the presence of viscosity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Real-space atomic dynamics in metallic liquids investigated by inelastic neutron scattering

Understanding the dynamics of liquids at the atomic level remains a major challenge. Even though viscosity is one of the most fundamental properties of liquids, its atomistic origin is not fully elucidated. Through inelastic neutron scattering experiment on levitated metallic liquid droplets, the time-dependent pair correlation function, the Van Hove function, was determined for Z⁢r 50 ⁢C⁢u 50 and Z⁢r 80 ⁢P⁢t 20 liquids at various temperatures. The time for change in local atomic connectivity, τ LC , which is the timescale of atomic bond cutting and forming, is estimated based on the exponential decay of the nearest neighbor peak of the Van Hove function. At high temperatures above the crossover temperature T A , τ LC is equal to the Maxwell relaxation time, τ M =η/G ∞ , where η is the macroscopic shear viscosity and G ∞ is the high-frequency shear modulus. Below T A the ratio of τ M /τ LC increases with decreasing temperature, indicating increased atomic cooperativity as predicted by molecular dynamics simulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Simulations of Stochastic Fluid Dynamics near a Critical Point in the Phase Diagram

Here, we present simulations of stochastic fluid dynamics in the vicinity of a critical endpoint belonging to the universality class of the Ising model. This study is motivated by the challenge of modeling the dynamics of critical fluctuations near a conjectured critical endpoint in the phase diagram of quantum chromodynamics (QCD). We focus on the interaction of shear modes with a conserved scalar density, which is known as model H. We show that the observed dynamical scaling behavior depends on the correlation length and the shear viscosity of the fluid. As the correlation length is increased or the viscosity is decreased we observe a crossover from the dynamical exponent of critical diffusion, z≃4, to the expected scaling exponent of model H, z≃3. We use our method to investigate the time-dependent correlation function of non-Gaussian moments M n (t) of the order parameter. We find that the relaxation time depends in a nontrivial manner on the power n.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Dynamics and Structure of Molten Neodymium Titanate

Binary rare earth titanates can be vitrified using levitation melt quenching, and their glasses are promising for optical applications. To better understand their thermophysical properties and dynamics, we have used molecular dynamics (MD) and empirical potential structure refinement (EPSR) to simulate molten neodymium titanate, 83TiO 2 –17Nd 2 O 3 . The models are compared with high-energy synchrotron X-ray diffraction of the melt structure over 1060–1960 K and thermophysical property measurements of electrostatically levitated droplets over 1500–2200 K. As temperature increases, the network periodicity increases and the correlation length decreases, both in monotonic trends that suggest continuous structural evolution through the melt and supercooled liquid regimes. At 1900 K, the mean Ti–O and Nd–O coordination numbers are 5.8(5) and 6.7(5), compared to the glass, 5.72(6) and 7.70(26). The equilibrium melting point is measured as 1860(20) K, over 100 K higher than previous estimates. Viscosity data from MD are in good agreement with experimental measurements. The dynamic propensity of Nd is similar to that of the network-forming Ti, suggesting that Nd does not participate as a conventional modifier. The viscosity of this fragile liquid is likely controlled by the breaking of cation–oxygen bonds and rearrangements of Ti–O polyhedra.

molecular dynamics↗
DOE OSTI2025

Impact of the LiPF6 Concentration on the Interfacial Charge Transfer and Fast-charging Capabilities of Lithium-Ion Batteries

Fast-charging lithium-ion batteries (LIB) demand optimized electrolyte formulations to balance ionic conductivity, viscosity, and interfacial charge transfer kinetics. This study examines how LiPF 6 concentration shapes solvation structure, desolvation energy, charge transfer activation energy, and solid electrolyte interphase (SEI) properties, which are critical for fast-charging performance. Using Raman spectroscopy, electrochemical cycling, X-ray photoelectron spectroscopy, and atomistic modeling, we analyze how varying LiPF 6 concentrations impact interfacial and bulk transport properties. Our findings show that increasing LiPF 6 concentration alters lithium solvation structures, reduces desolvation energy, and accelerates charge transfer at the electrode interface. Higher concentrations lower the activation energy for charge transfer and suppress excessive SEI growth, improving interfacial kinetics. However, concentrations above a certain threshold increase viscosity and reduce ionic conductivity, limiting transport efficiency. These results offer insights into electrolyte solvation and interfacial charge transfer mechanisms, providing guidelines for designing next-generation fast-charging LIB electrolytes with enhanced efficiency, stability, and longevity.

Son, Seoung-Bum [Argonne National Laboratory (ANL)↗
DOE OSTI2024

An Overview of the Molten Salt Thermal Properties Database--Thermophysical, Version 3.1 (MSTDB-TP v.3.1)

This report presents the current status of the Molten Salt Thermal Properties Database–Thermophysical (MSTDB-TP). Information regarding version 3.1 is provided herein, which contains 820 individual salt entries (data from 180+ independent studies); the thermophysical properties contained in the database include density, viscosity, thermal conductivity, and heat capacity. The major updates to the database include a significant expansion of pseudobinary and higher-order chloride salt mixtures, many of which bearing actinides, and an incorporation of more recent literature data (i.e., that within the past 5 years). Also, modifications have been made to the pure compound data in the database as a consequence of an external quality assessment of duplicate datasets. The user-facing API for the MSTDB-TP, Saline, has been updated to include viscosity estimation capabilities based on the Redlich–Kister formalism; this is an advancement with respect to the existing density estimation capabilities. The graphical user interface was also updated to include a density estimation capability, backed by Saline. Finally, additional preliminary efforts to include surface tension into the database, as well as an investigation on formalisms that would be appropriate for thermal conductivity estimation, are reported herein.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development of a multiphase PIC model for slurry flow modeling

MFIX-Exa is a recently released multiphase CFD code originally developed for the simulation of particle-laden gas-solid flows. Due to its high-performance computing capabilities, MFIX-Exa is an ideal candidate for scale-up studies of slurry reactors, specifically the coarse-grained particle-in-cell (PIC) model with its statistical treatment of the particle phase. Unfortunately, several physical models that were neglected during original development because they are not relevant for high-density ratio gas-solid flows are important in slurry flows where the particle-to-fluid density ratio is near unity. In this preliminary work we focus on the effective (suspension) viscosity. The models of Brinkman (1952), Krieger and Dougherty (1956), and Cheng and Law (2003) are considered. The impact of the effective viscosity model is studied on horizontal pipe flow. The experimental data of Gillies et al. (2002) is used to assess the pressure drop predictions.

Fullmer, William D.↗
DOE OSTI2025

Investigations into the Ternary NaF-KF-UF4 Salt System – Phase A

A knowledge gap exists in the data and understanding of fresh fuel salt and irradiated multicomponent fuel salt systems thermophysical properties. Quantifying these properties is necessary for the design and construction of test reactors, as well as the licensing of future commercial molten-salt reactors. To facilitate thermal property determination on a proposed fuel salt composition for Seaborg Technologies, several samples containing depleted uranium tetrafluoride (UF4), sodium fluoride (NaF), and potassium fluoride (KF) were blended, and a melt temperature analysis was performed. From the melting temperature analysis, it was determined that sample Seaborg-7, a ternary salt composition of 26.4UF4-24.7KF-48.9NaF (mol%), was very near a ternary eutectic point. Therefore, thermal properties such as melting temperature, salt stability, density, heat capacity, thermal diffusivity, and viscosity were experimentally determined on the Seaborg-7 salt. These measurements document the baseline properties of fresh fuel salt as a function of temperature, where future experiments on irradiated fuel salt will provide a holistic perspective on the change of thermophysical properties during reactor operations. Several precision instruments were used to collect property data, and instrument calibrations and data collection were performed and documented in a standardized and reproducible manner with meticulous detail. This process ensured that the measurement procedures and resulting data can readily be duplicated elsewhere. The Seaborg-7 salt was shown to be stable at temperatures up to 900°C, as no mass change was observed upon repeated heating and cooling. The peak melting temperature was determined to be 547°C (557°C endset). The enthalpy of fusion (??H?_fus^o) was determined to be 167.5 ± 2.7 J/g while the enthalpy of crystallization (??H?_c^o) was determined to be -147.8 ± 13.3 J/g. In addition to the eutectic melting peak, upon heating, several pre eutectic peaks were observed, occurring at 470°C (onset) and 499°C (peak). Specific heat capacity measurements showed a slightly increasing trend with respect to temperature in the solid phase, while the liquid-specific heat capacity showed a somewhat flat trend with an average value of 106.1 ± 1.24 J/mol·K between 600 to 800°C. Three independent trials using the Seaborg-7 salt determined the density to be ?(T) = 4.908 – 0.000363·T(°C), validated between 32 to 200°C, and ?(T) = 4.808 – 0.00113·T(°C), validated between ~575 to 850°C. Thermal diffusivity was determined for the liquid state and is represented by the linear equation y = 0.1581 + 0.000207·T(°C) between 550 to 850°C. The viscosity was determined from 600 to 800°C and is represented by the exponential fit equation, ? (mPa·s) = 736.58e^(-0.006·T(°C)). This report documents the conclusion of fuel salt thermophysical property measurements for the Seaborg SPP, Phase A project.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Quality Ranking of Unary Chloride Salt Property Data Included in MSTDB-TP

Molten salt reactor developers rely on thermal property data to design, license and operate the reactors. The Molten Salt Thermal Database-Thermophysical Properties (MSTDB-TP) was established under the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program and is managed by Oak Ridge National Laboratory to serve as a single source of thermophysical property values measured for a wide variety of molten salt systems for use by researchers, molten salt reactor developers, and regulators. These properties include density, viscosity and thermal diffusivity and conductivity. Published measurements of molten salt properties are lacking for many salts of interest and the data that are available are often inconsistent. This creates a challenge for MSR developers when determining which property values to use when designing their reactors. It is the purpose of this work to apply a consistent ranking system to all data entries that indicates the quality of property values listed in the database. These rankings will be the technical basis for down-selections by the database developers and alert users about the quality of the available property values. MSTDB-TP collects all available property data and indicates preferred data sets or correlations. However, all available data sets are included in the database. Quality assessments and rankings are being applied to data in MSTDB-TP to provide an indication of the quality of each data set independent of consistency with other data. Previous reports detailed the ranking system that was followed and assessments of unary fluoride data sets. Documentation of the quality of data in MSTDB-TP was continued by reviewing and assessing all available sources of density, viscosity and thermal diffusivity or conductivity values for unary chloride salts in MSTDB-TP V3.0 using the same criteria.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Predicting Selective Laser Printing Print Quality of Polymer Powders through Melt Flow Index

Selective Laser Sintering (SLS) uses a precisely controlled laser to fuse polymer powder to build complex 3D shapes. While SLS covers a wide application space, the processing knowledge of polymer powder is limited, restricting the number of commercial powders available. This study serves to expand the processing knowledge of polypropylene and polyethylene in parallel with a “mature” SLS feedstock, nylon, through melt flow index (MFI) characterization. Differential Scanning Calorimetry (DSC) was used to explore the sintering window of the polymers. Polyethylene and polypropylene exhibited a relatively narrow sintering window between 4 – 5 °C, whereas the sintering window of nylon was much wider, between 23 – 24 °C. This suggests that the print quality between polyethylene and polypropylene would be similar; however, X-ray computed tomography revealed a higher volume of print defects, voids, and de lamination in polyethylene than in polypropylene. MFI analysis provided additional insight into the difference in print quality, as the MFI of polyethylene was 7.16 g/10 min, 9.95 g/10 min for polypropylene, and 17.23 g/10 min for nylon. MFI is inversely correlated to melt viscosity, and a low melt viscosity is desired for proper coalescing between the layers and particles. These results suggest that MFI is a promising tool, in conjunction with traditional thermal analysis, for screening candidate powder feedstocks for SLS and optimization of print parameters for novel powders.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Exploring the Feasibility of INCONEL® ALLOY 740H® for Power Plant Headers: Integrating Machine Learning with Computational Fluid Dynamics (CFD)

This keynote presentation explores the behavior of headers—essential components of pipeline systems—using ANSYS simulation software and machine learning techniques. The study aims to predict the thermal and mechanical performance of headers under diverse conditions through both steady-state and transient simulations. We investigate critical parameters such as heat transfer coefficient, fluid velocity, and temperature to optimize header design. Conducted as part of a DOE project led by NCAT in collaboration with UNC Charlotte, this research encompasses multiple key topics. The initial section focuses on the behavior of header systems under steady-state conditions using ANSYS simulation. It underscores the importance of headers in industrial infrastructure, especially in the energy sector, and examines the implications of material selection and flow direction on heat transfer dynamics. Methodologically, we employ Computational Fluid Dynamics (CFD) analysis through ANSYS, detailing the development of models, material properties, geometry specifications, boundary conditions, and meshing strategies. Our simulations explore various operational parameters, including temperature and mass flow rates, crucial for predicting heat transfer coefficients and enhancing header design. Results from the study include parametric investigations into mesh sensitivity, viscosity model evaluations, and the effects of heat transfer locations, all validated against theoretical calculations. We conclude with insights on mesh optimization, the suitability of viscosity models, and recommendations for future research aimed at improving header system efficiency and sustainability in industrial applications.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Spectrally Stabilized Interface Capturing Formulation and Implementation in Nek5000/NekRS

This report documents the formulation of a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation. The two-phase framework conceptualized herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al, is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The two-phase formulation is wholly implemented in Nek5000 and the SVV stabilization method is implemented in NekRS, which is the essential precursor to the two-phase framework, undergoing active development.

97 MATHEMATICS AND COMPUTING↗