Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “viscosity, shear”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Calculation of shock-separated turbulent boundary layers

Numerical solutions of the complete, time-averaged conservation equations using several eddy-viscosity models for the Reynolds shear stress to close the equations are compared with experimental measurements in a compressible, turbulent separated flow. An efficient time-splitting, explicit difference scheme was used to solve the two-dimensional conservation equations. The experiment used for comparison was a turbulent boundary layer that was separated by an incident shock wave in a Mach 2.93 flow with a unit Reynolds number of 5.7 x 10 to the seventh power m. Comparisons of predicted and experimental values of surface pressure, shear stress along the wall, and velocity profiles are shown. One of the tested eddy-viscosity models which allows the shear stress to be out of equilibrium with the mean flow produces substantially better agreement with the experimental measurements than the simpler models. A tool is thereby provided for inferring additional information about the flow, such as static pressures in the stream, which might not be directly obtainable from experiments.

Baldwin, B. S.↗
DOE OSTI2021

Toward Optimizing Electrospun Nanofiber Fuel Cell Catalyst Layers: Polymer–Particle Interactions and Spinnability

We investigate the effect of the poly(acrylic acid) (PAA) carrier polymer concentration on the microstructure and rheological properties of catalyst inks for electrospun polymer–electrolyte membrane fuel-cell catalyst layers. Characterization of an ink microstructure using oscillatory shear rheology showed that the catalyst particles (platinum on carbon) are significantly agglomerated in the absence of PAA or an ionomer. Both the ionomer and PAA promoted the stability of the particles against agglomeration via electrosteric stabilization by adsorbing onto the particle surface. Increasing the PAA concentration increased the stability of the particles (or reduced the agglomerated structure) due to increasing PAA coverage onto the free surface area of the particles. However, beyond a certain increase in concentration, PAA was found to predominantly remain as an excess free polymer in the ink due to an insufficient free/available surface area on the particles for further PAA coverage. Extensional rheology measurements demonstrated that PAA enhances the extensional viscosities of the inks. Consequently, increasing the PAA concentration in the ink promoted the evolution of uniform nanofibers. However, beyond a certain concentration, a significant increase in the shear viscosities of the inks led to defective fiber morphologies because of the onset of flow instabilities. Electrochemical performance comparisons between catalyst layers with different PAA concentrations showed maximum performance at the PAA concentration that led to the least agglomerated structure of the catalyst, most uniform fiber morphologies, and low concentrations of free (non-adsorbing) PAA in the electrode. These results provide a rationale for optimization of electrospun catalyst nanofibers for both spinnability and electrochemical performance.

25 ENERGY STORAGE↗
NASA NTRS1992

Search for subgrid scale parameterization by projection pursuit regression

The dependence of subgrid-scale stresses on variables of the resolved field is studied using direct numerical simulations of isotropic turbulence, homogeneous shear flow, and channel flow. The projection pursuit algorithm, a promising new regression tool for high-dimensional data, is used to systematically search through a large collection of resolved variables, such as components of the strain rate, vorticity, velocity gradients at neighboring grid points, etc. For the case of isotropic turbulence, the search algorithm recovers the linear dependence on the rate of strain (which is necessary to transfer energy to subgrid scales) but is unable to determine any other more complex relationship. For shear flows, however, new systematic relations beyond eddy viscosity are found. For the homogeneous shear flow, the results suggest that products of the mean rotation rate tensor with both the fluctuating strain rate and fluctuating rotation rate tensors are important quantities in parameterizing the subgrid-scale stresses. A model incorporating these terms is proposed. When evaluated with direct numerical simulation data, this model significantly increases the correlation between the modeled and exact stresses, as compared with the Smagorinsky model. In the case of channel flow, the stresses are found to correlate with products of the fluctuating strain and rotation rate tensors. The mean rates of rotation or strain do not appear to be important in this case, and the model determined for homogeneous shear flow does not perform well when tested with channel flow data. Many questions remain about the physical mechanisms underlying these findings, about possible Reynolds number dependence, and, given the low level of correlations, about their impact on modeling. Nevertheless, demonstration of the existence of causal relations between sgs stresses and large-scale characteristics of turbulent shear flows, in addition to those necessary for energy transfer, provides important insight into the relation between scales in turbulent flows.

Meneveau, C.↗
NASA NTRS1978

Rheology of composite solid propellants during motor casting

Results of casting studies are reviewed so as to define the viscosity criteria insuring the fabrication of defect-free grains. The rheology of uncured propellants is analyzed showing that a realistic assessment of a propellant's flow properties must include measurement of viscosity as a function of shear stress and time after curing agent. Methods for measuring propellant viscosity are discussed, with particular attention given to the Haake-Rotovisko rotational viscometer. The effects of propellant compositional and processing variables on apparent viscosity are examined, as are results relating rheological behavior to grain defect formation during casting.

Klager, K.↗
DOE OSTI2025

Sedimentation and shear-induced dynamics of spheroids in fluids with spatial viscosity variations

A generalized reciprocal theorem is used to relate the force and torque induced on a particle in an inertia-less fluid with small variation in viscosity to integrals involving Stokes flow fields and the spatial dependence of viscosity. These resistivity expressions are analytically evaluated using spheroidal harmonics and then used to obtain the mobility of the spheroid during sedimentation, and in linear flows, of a fluid with linear viscosity stratification. The coupling between the rotational and translational motion induced by stratification rotates the spheroid’s centerline, creating a variety of rotational and translational dynamics dependent upon the particle’s aspect ratio, κ, and the component of the stratification unit vector in the gravity direction, d g . Spheroids with 0.55 ⪅ κ ⪅ 2.0 exhibit the largest variety of settling behaviors. Interestingly, this range covers most microplastics and typical microorganisms. One of the modes include a stable orientation dependent only on κ and d g , but independent of initial orientation, thus allowing for the potential control of settling angles and sedimentation rates. In a simple shear flow, cross-streamline migration occurs due to the stratification-induced force generated on the particle. Similarly, a particle no longer stays at the stagnation point of a uniaxial extensional flow. While fully analytical results are obtained for spheroids, numerical simulations provide a source of validation. These simulations also provide additional insights into the stratification-induced force- and torque-producing mechanisms through the stratification-induced stress, which is not accessed in the reciprocal theorem-based analytical calculations.

Geophysical and Geological Flows: Stratified flows↗
DOE OSTI2026

Insights into mixing of non-isothermal multi-polymer melts for complex plastics recycling

Catalytic recycling or upcycling of plastics is often limited not by catalyst performance, but by transport, arising from highly viscous, non-Newtonian polymer melts. In this work, we develop a reactor-scale framework that integrates rheological measurements, constitutive modeling, computational fluid dynamics (CFD), and experiments to quantify mixing, heat transfer, and dispersion in surrogate hydrocarbon melts representing mixed plastics systems. Temperature- and shear rate-dependent viscosity of low-density polyethylene (LDPE) and high-density polyethylene (HDPE) is measured to create two surrogate polymers (PLD and PHD) that capture the dominant shear-thinning flow behavior while neglecting strong elastic effects, enabling tractable simulation of non-isothermal, polymer-melt mixing using a Carreau-Arrhenius generalized Newtonian framework. Three-dimensional CFD simulations are employed to evaluate impeller performance in PLD using mixing time, cavern volume, thermal uniformity, and interfacial area for regimes in which viscoelastic effects are not dominant. We show that magnetic stir bars commonly used in lab-scale studies produce large thermal gradients (~60 °C) and poor mixing, even under idealized power delivery and polymer flow conditions. In contrast, close-clearance anchor impellers achieve near-isothermal operation, reduce mixing times by up to 5×, and provide >90% active circulation volume. We further demonstrate that, at low pseudo-Deborah number (De*), motor power requirements can be predicted directly from shear rate-dependent rheology using the Carreau-Arrhenius framework, enabling rational selection of operating conditions. Extension to surrogate immiscible multi-polymer systems based on PLD and PHD shows that interfacial area is highly sensitive to operating conditions and impeller design, with coaxial anchor-turbine configurations enhancing dispersion by up to 4 × .

Close-clearance impellers↗
DOE OSTI2020

Molecular origins of bulk viscosity in liquid water

The shear viscous response of water is closely associated with changes in network connectivity on the sub ps timescale. The bulk viscous response is shown here to be associated with local density fluctuations and rotational motion around 1–3 ps.

Yahya, Ahmad↗
DOE OSTI2020

Viscous flow properties and hydrodynamic diameter of phenothiazine-based redox-active molecules in different supporting salt environments

We report viscous flow properties of a redox-active organic molecule, N-(2-(2-methoxyethoxy)ethyl)phenothiazine (MEEPT), a candidate for non-aqueous redox flow batteries, and two of its radical cation salts. A microfluidic viscometer enabled the use of small sample volumes in determining viscosity as a function of shear rate and concentration in the non-aqueous solvent, acetonitrile, both with and without supporting salts. All solutions tested show Newtonian behavior over shear rates of up to 30 000 s−1, which was rationalized by scaling arguments for the diffusion-based relaxation time of a single MEEPT molecule without aggregation. Neat MEEPT is flowable but with a large viscosity (412 mPa⋅s at room temperature), which is ∼1000 times larger than that of acetonitrile. MEEPT solutions in acetonitrile have low viscosities; at concentrations up to 0.5 M, the viscosity increases by less than a factor of two. From concentration-dependent viscosity measurements, molecular information was inferred from intrinsic viscosity (hydrodynamic diameter) and the Huggins coefficient (interactions). Model fit credibility was assessed using the Bayesian Information Criterion. It is found that the MEEPT and its charged cations are “flowable” and do not flocculate at concentrations up to 0.5 M. MEEPT has a hydrodynamic diameter of around 8.5 Å, which is almost insensitive to supporting salt and state of charge. This size is comparable to molecular dimensions of single molecules obtained from optimized structures using density functional theory calculations. The results suggest that MEEPT is a promising candidate for redox flow batteries in terms of its viscous flow properties.

Wang, Yilin (ORCID:0000000316102995)↗
NASA NTRS1973

Investigations of lubricant rheology as applied to elastohydrodynamic lubrication

The pressure viscometer was modified to permit the measurement of viscosity at elevated pressures and shear stresses up to 5 x 10 to the 6th power N/sq m (720 psi). This shear stress is within a factor of three of the shear stress occurring in a sliding ehd point contact such as occurs in the ehd simulator. Viscosity data were taken on five lubricant samples, and it was found that viscous heating effects on the viscosity were predominant and not non-Newtonian behavior at the high shear stresses. The development of the infrared temperature measuring technique for the ehd simulator was completed, and temperature data for a set of operating conditions and one lubricant are reported. The numerical analysis of the behavior of nonlinear lubricants in the lubrication of rollers is reported.

Carlson, S.↗
NASA NTRS1994

Two-layer viscous instability in a rotating couette device

A novel experiment to study the interfacial shear instability between two liquids is described. Density-matched immiscible liquids are confined between concentric cylinders such that the interface is parallel to the cylinder walls. Interfacial waves that develop because of viscosity differences between the shearing fluids are studied as a function of rotation rate and depth ratio using optical techniques. Conditions neutral stability and the most unstable wavenumber agree reasonably well with predictions from linear stability analysis of the Navier-Stokes equations. Illumination using laser sheets allows precise measurement of the interface shape. Future experiments will verify the correctness of weakly nonlinear theories that describe energy transfer and saturation of wave growth by nonlinear effects. Measurements of solitary wave shapes, that occur far above neutral stability, will be compared to similar measurements for systems that have gravity as an important force to determine how gravity effects large disturbances. These results will be used to interpret slug and annular flow data that have been obtained in other mu g studies.

Gallagher, Christopher↗
NASA NTRS2006

Development and Implementation of Non-Newtonian Rheology Into the Generalized Fluid System Simulation Program (GFSSP)

One of the goals of this program was to develop the experimental and analytical/computational tools required to predict the flow of non-Newtonian fluids through the various system components of a propulsion system: pipes, valves, pumps etc. To achieve this goal we selected to augment the capabilities of NASA's Generalized Fluid System Simulation Program (GFSSP) software. GFSSP is a general-purpose computer program designed to calculate steady state and transient pressure and flow distributions in a complex fluid network. While the current version of the GFSSP code is able to handle various systems components the implicit assumption in the code is that the fluids in the system are Newtonian. To extend the capability of the code to non-Newtonian fluids, such as silica gelled fuels and oxidizers, modifications to the momentum equations of the code have been performed. We have successfully implemented in GFSSP flow equations for fluids with power law behavior. The implementation of the power law fluid behavior into the GFSSP code depends on knowledge of the two fluid coefficients, n and K. The determination of these parameters for the silica gels used in this program was performed experimentally. The n and K parameters for silica water gels were determined experimentally at CFDRC's Special Projects Laboratory, with a constant shear rate capillary viscometer. Batches of 8:1 (by weight) water-silica gel were mixed using CFDRC s 10-gallon gelled propellant mixer. Prior to testing the gel was allowed to rest in the rheometer tank for at least twelve hours to ensure that the delicate structure of the gel had sufficient time to reform. During the tests silica gel was pressure fed and discharged through stainless steel pipes ranging from 1", to 36", in length and three diameters; 0.0237", 0.032", and 0.047". The data collected in these tests included pressure at tube entrance and volumetric flowrate. From these data the uncorrected shear rate, shear stress, residence time, and viscosity were evaluated using formulae for non-Newtonian, power law fluids. The maximum shear rates (corrected for entrance effects) obtained in the rheometer with the current setup were in the 150,000 to 170,000sec- range. GFSSP simulations were performed with a flow circuit simulating the capillary rheometer and using Power Law gel viscosity coefficients from the experimental data. The agreement between the experimental data and the simulated flow curves was within +/-4% given quality entrance effect data.

DiSalvo, Roberto↗
NASA NTRS1974

Dissipative heating effects and end corrections for viscous Newtonian flow in high shear capillary tube viscometry

The effect of dissipation heating on the apparent viscosity measured in capillary tube viscometry is described in this paper. Conditions of low Reynolds number and high shear are assumed. End corrections to the tube flow are incorporated. The flow curves show decreasing apparent viscosity when the shear stress increases. The configuration of the flow curves plotted in logarithmic presentation are found to be identical for fluids with Newtonian behavior. Convection is the predominant mechanism in removal of the heat in short capillary tube. The estimated upper bound for the shear stress obtainable in short length capillary tubes appears to be of the order of magnitude of 10 MPa limited primarily by the pressure drop associated with the constant end correction from the flat ended inlet and exit of the tube.

Jakobsen, J.↗
DOE OSTI2022

With a Little Help from My Cell Wall: Structural Modifications in Pectin May Play a Role to Overcome Both Dehydration Stress and Fungal Pathogens

Cell wall structural modifications through pectin cross-linkages between calcium ions and/or boric acid may be key to mitigating dehydration stress and fungal pathogens. Water loss was profiled in a pure pectin system and in vivo. While calcium and boron reduced water loss in pure pectin standards, the impact on Allium species was insignificant (p > 0.05). Nevertheless, synchrotron X-ray microscopy showed the localization of exogenously applied calcium to the apoplast in the epidermal cells of Allium fistulosum. Exogenous calcium application increased viscosity and resistance to shear force in Allium fistulosum, suggesting the formation of calcium cross-linkages (“egg-box” structures). Moreover, Allium fistulosum (freezing tolerant) was also more tolerant to dehydration stress compared to Allium cepa (freezing sensitive). Furthermore, the addition of boric acid (H 3 BO 3 ) to pure pectin reduced water loss and increased viscosity, which indicates the formation of RG-II dimers. The Arabidopsis boron transport mutant, bor1, expressed greater water loss and, based on the lesion area of leaf tissue, a greater susceptibility to Colletotrichum higginsianum and Botrytis cinerea. While pectin modifications in the cell wall are likely not the sole solution to dehydration and biotic stress resistance, they appear to play an important role against multiple stresses.

59 BASIC BIOLOGICAL SCIENCES↗
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↗
NASA NTRS1988

Shear rheological characterization of motor oils

Measurements of high pressure viscosity, traction coefficient, and EHD film thickness were performed on twelve commercial automotive engine oils, a reference oil, two unformulated base oils and two unformated base oil and polymer blends. An effective high shear rate inlet viscosity was calculated from film thickness and pressure viscosity coefficient. The difference between measured and effective viscosity is a function of the polymer type and concentration. Traction measurements did not discriminate mileage formulated oils from those not so designated.

Bair, Scott↗
DOE OSTI2020

Active- and transfer-learning applied to microscale-macroscale coupling to simulate viscoelastic flows

Active- and transfer-learning are applied to microscale dynamics of polymer flows for the multiscale discovery of effective constitutive approximations required in viscoelastic flow simulation. The result is macroscopic rheology directly connected to a microstructural model. Micro and macroscale simulations are adaptively coupled by means of Gaussian process regression (GPR) to run the expensive microscale computations only as necessary. This multiscale method is demonstrated with flows of a polymer solution as a model system. At the microscale level dissipative particle dynamics (DPD) is employed to model the fluid as a suspension of bead-spring micro-structures subjected to steady shear flow. The results yield the non-Newtonian viscosity and the first normal stress difference at strain rates as training data used in a GPR model. DPD parameters are calibrated with respect to experimental data for a real polymer solution. Compliance with these data requires adjustment of the DPD model's cutoff radius, which then becomes a function of the second invariant of the strain rate tensor. The FENE-P model is chosen for the macroscale description using the spectral element method (SEM) to simulate channel flow and flow past a circular cylinder. The DPD results at the lowest possible shear strain rate yield an estimate of the zero-shear rate viscosity, which allows the initiation of the macroscale flow by SEM as a Newtonian fluid. The resulting strain-rate field is surveyed to determine additional shear strain rate sampling points for the DPD system. This new information allows an initial fitting of parameters of the constitutive equation followed by new SEM simulations at the macroscale. Additionally, guided by active-learning GPR to select new sampling points, this process continues until convergence is achieved. The effectiveness of this new simulation paradigm for viscoelastic flows is tested with different macroscale operating conditions. The effective closure learned in the channel simulation is then transferred directly to the flow past a circular cylinder at low Reynolds number, where the results show that only two additional DPD simulations are required to achieve a satisfactory constitutive model. With an increase of the Reynolds number, the active-learning scheme automatically detects the inaccuracy of the learned constitutive model, and initiates additional DPD simulations for the extra data needed to once again close the microscale-macroscale coupled system. This new paradigm of active- and transfer-learning for multiscale modeling is readily applicable to other microscale-macroscale coupled simulations of complex fluids and other materials. Furthermore, the coupling between microscale and macroscale solvers can be seamlessly implemented with our open source multiscale universal interface (MUI) library.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A constitutive model for sheared dense suspensions of rough particles

In a concentrated suspension, particles come into contact due to the presence of asperities on their surfaces. As a result, the contact forces and interparticle friction become one of the important factors governing the rheology of rigid particle suspensions at high concentrations. We show that a load-dependent friction model can be used to reproduce the experimentally observed shear thickening [ST—continuous and discontinuous (DST)] behavior with an increase in the asperity size. Increasing the particle surface roughness size leads to a decrease in the critical shear rate above which shear-thickening takes place, a reduction in the critical volume fraction for DST and an increase in the viscosity jump across non-shear-thickened to shear-thickened regime. In this paper, we propose a constitutive model to quantify the effect of increasing the roughness size on the rheology of dense ST suspensions as well as on the critical shear rate for ST and the critical volume fraction for DST. We fit this model to our simulation data for stress controlled shear flow of dense rough particle suspensions. Once the fitting is complete, these equations are used to predict exact volume fractions and shear stress values for transitions between three regimes on the shear stress-shear rate flow state diagram for different roughness values. Furthermore, the results of this study can be used to tune the particle surface roughness for manipulating the dense suspension rheology according to different applications.

42 ENGINEERING↗
DOE OSTI2026

Reduced-order CFD modeling of cryogenic hydrogen isotope extrusion for pellet fueling

This study presents a reduced-order model (ROM) for computational fluid dynamics (CFD) simulations of cryogenic hydrogen isotope extrusions, focusing on protium (H₂) and deuterium (D₂) piston extruders. Using a 2D axisymmetric ROM in ANSYS-Polyflow, significant computational savings were achieved (runtime reduced from 9∼24 h to 3∼5 min), with extrusion force discrepancies between the 2D ROM and 3D models being on the order of 1%. Parametric studies identified optimal cutoff shear rates in the viscosity model (0.01/s for H₂ and 0.001/s for D₂), providing recommendations for future simulations. Finally, a comprehensive comparison of ROM results with experimental data was performed across varying geometries, cryogenic materials, temperatures, extrusion lengths, and piston velocities. Predictions at low extrusion temperatures met the objective of providing quick and efficient solutions with an acceptable extrusion force error of approximately 10% or less, validating the effectiveness of the 2D ROM approach. However, at high temperatures closer to the triple point, extrusion force error grows, which necessitates developing an improved model that accounts for temperature effects, e.g. melting. Nevertheless, the findings still represent a significant improvement in efficiency of CFD modeling of cryogenic hydrogenic extrusion. The ROM framework can also be extended to tritium (T2) and screw extruders, which will ultimately provide a fast and effective tool for optimizing pellet injector design for ITER and future reactor systems.

Fan, Joy [ORNL] (ORCID:0000000229751735)↗