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.

SourceResearch connectionFrom yearThrough year

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

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Buoyancy Driven Shear Flows of Bubble Suspensions

In this work the gas volume fraction and the root-mean-squared fluid velocity are measured in buoyancy driven shear flows of bubble suspensions in a tall, inclined, rectangular channel. The experiments are performed under conditions where We << 1a nd Re >> 1, for which comparisons are made with kinetic theory and numerical simulations. Here Re = gamma(a(exp 2)/nu is the Reynolds number and We = rho(gamma(exp 2))a(exp 3)/sigma is the Weber number; gamma is the shear rate, a is the bubble radius, nu is the kinematic viscosity of the liquid, rho is the density of the liquid, and sigma is the surface tension of the gas/liquid interface. Kang et al. calculated the bubble phase pressure and velocity variance of sheared bubble suspensions under conditions where the bubbles are spherical and the liquid phase velocity field can be approximated using potential flow theory, i.e. We= 0 and Re >> 1. Such conditions can be achieved in an experiment using gas bubbles, with a radius of O(0.5mm), in water. The theory requires that there be no average relative motion of the gas and liquid phases, hence the motivation for an experimental program in microgravity. The necessity of performing preliminary, Earth based experiments, however, requires performing experiments where the gas phase rises in the liquid, which significantly complicates the comparison of experiments with theory. Rather than comparing experimental results with theory for a uniform, homogeneous shear flow, experiments can be compared directly with solutions of the averaged equations of motion for bubble suspensions. This requires accounting for the significant lift force acting on the gas phase when the bubbles rise parallel to the average velocity of the sheared suspension. Shear flows can be produced in which the bubble phase pressure gradient, arising from shear induced collisions amongst the bubbles, balances a body force (centrifugal or gravitational) on the gas phase. A steady, non-uniform gas volume fraction can be measured, from which the bubble phase pressure gradient can be obtained and compared to theory and numerical simulations. The presence of bounding walls further complicates the experiments, since the detailed interactions of the bubbles with bounding walls is not well understood, especially in the presence of gravity, where the momentum and energy exchange depends on the inclination of the wall.

Koch, D. L.↗

NASA NTRS1993

Can weak crust explain the correlation of geoid and topography on Venus?

The effect on geoid and topography of low viscosity crust overlying a steady-state convecting mantle is estimated under the assumption that the shear between crust and mantle does not alter the mantle flow. The weak crustal layer can change the sign of the geoid to topography ratio (admittance). The positive long wavelength admittance for Venus is consistent with a weak crust overlying a mantle with a viscosity that increases strongly with depth. The accepted interpretation of the strong positive correlation of geoid and topography on Venus, is that the convecting mantle of Venus has a constant viscosity with depth. Topography results from vertical normal stresses caused by mantle convection and highlands occur where mantle upwells. For topography to be supported by normal stress, the time scale for crustal flow must be long compared to the time scale for changes in the pattern of mantle flow. Because the high surface temperature of Venus may cause the crust to have a low viscosity, this assumption may be false. Topography should then be dominated by shear coupling between the crust and mantle. In the absence of a crustal layer, convection in a constant viscosity layer gives rise to a geoid anomaly that correlates positively with surface topography. When the viscosity in the layer increases with depth by several orders of magnitude, the surface topography and geoid anomaly become anti-correlated.

Buck, W. Roger↗

NASA NTRS1994

Nonlinear Instability of a Uni-directional Transversely Sheared Mean Flow

It is well known that the presence of a weak cross flow in an otherwise two-dimensional shear flow results in a spanwise variation in the mean streamwise velocity profile that can lead to an amplification of certain three-dimensional disturbances through a kind of resonant-interaction mechanism (Goldstein and Wundrow 1994). The spatial evolution of an initially linear, finite-growth-rate, instability wave in such a spanwise-varying shear flow is considered, The base flow, which is governed by the three-dimensional parabolized Navier-Stokes equations, is initiated by imposing a spanwise- periodic cross-flow velocity on an otherwise two-dimensional shear flow at some fixed streamwise location. The resulting mean-flow distortion initially grows with increasing streamwise distance, reaches a maximum and eventually decays through the action of viscosity. This decay, which coincides with the viscous spread of of the shear layer, means that the local growth rate of the instability wave will eventually decrease as the wave propagates downstream. Nonlinear effects can then become important within a thin spanwise-modulated critical layer once the local instability-wave amplitude and growth rate become sufficiently large and small, respectively. The amplitude equation that describes this stage of evolution is shown to be a generalization of the one obtained by Goldstein and Choi (1989) who considered the related problem of the interaction of two oblique modes in a two-dimensional shear layer.

Wundrow, David W.↗

DOE OSTI2024

Porous Polymer Structures with Tunable Mechanical Properties Using a Water Emulsion Ink

Recently, the manufacturing of porous polydimethylsiloxane (PDMS) with engineered porosity has gained considerable interest due to its tunable material properties and diverse applications. An innovative approach to control the porosity of PDMS is to use transient liquid phase water to improve its mechanical properties, which has been explored in this work. Adjusting the ratios of deionized water to the PDMS precursor during blending and subsequent curing processes allows for controlled porosity, yielding water emulsion foam with tailored properties. The PDMS-to-water weight ratios were engineered ranging from 100:0 to 10:90, with the 65:35 specimen exhibiting the best mechanical properties with a Young’s Modulus of 1.17 MPa, energy absorption of 0.33 MPa, and compressive strength of 3.50 MPa. This led to a porous sample exhibiting a 31.46% increase in the modulus of elasticity over a bulk PDMS sample. Dowsil SE 1700 was then added, improving the storage capabilities of the precursor. The optimal storage temperature was probed, with -60 °C resulting in great pore stability throughout a three-week duration. The possibility of using these water emulsion foams for paste extrusion additive manufacturing (AM) was also analyzed by implementing a rheological modifier, fumed silica. Fumed silica’s impact on viscosity was examined, revealing that 9 wt% of silica demonstrates optimal rheological behaviors for AM, bearing a viscosity of 10,290 Pa·s while demonstrating shear-thinning and thixotropic behavior. This study suggests that water can be used as pore-formers for PDMS in conjunction with AM to produce engineered materials and structures for aerospace, medical, and defense industries as sensors, microfluidic devices, and lightweight structures.

36 MATERIALS SCIENCE↗

NASA NTRS1975

Calculation of three-dimensional compressible laminar and turbulent boundary layers. Calculation of three-dimensional compressible boundary layers on arbitrary wings

A very general method for calculating compressible three-dimensional laminar and turbulent boundary layers on arbitrary wings is described. The method utilizes a nonorthogonal coordinate system for the boundary-layer calculations and includes a geometry package that represents the wing analytically. In the calculations all the geometric parameters of the coordinate system are accounted for. The Reynolds shear-stress terms are modeled by an eddy-viscosity formulation developed by Cebeci. The governing equations are solved by a very efficient two-point finite-difference method used earlier by Keller and Cebeci for two-dimensional flows and later by Cebeci for three-dimensional flows.

Cebeci, T.↗

NASA NTRS1985

On the importance of viscous dissipation in Io

A model of Io is presented that consists of an elastic inner core, a low strength asthenosphere, and a thin elastic outer shell. The middle layer is ssumed to posses negligible shear strength and to be characterized by a Newtonian viscosity. The fluid in the viscous layer is forced to circulate mainly by the tidal distortion in the outer shell, modeled here as a variation in the distortion amplitude. As a result, heat is generated in the fluid by viscous dissipation. There are three important unconstrained parameters in the model: the fluid viscosity, the thickness of the fluid layer, and the degree to which the distortion of the outer shell is affected by the fluid viscosity. For a wide range of these model parameters viscous heating can generate just as much or even more heat than does elastic dissipation is the outer shell. The model suggests that much of Io's heat flow may be generated below the outer shell and could provide a source of energy for any silicate volcanism on the satellite.

Ross, M. N.↗

NASA NTRS1988

Turbulent reconnection

The reconnection rate in the presence of tearing mode turbulence is calculated. Tearing modes can produce turbulence in three-dimensional sheared magnetic fields. They produce a turbulent anomalous electron viscosity or hyperresistivity. Using quasi-linear and direct interaction approximations, the scaling of the hyperresistivity with magnetic perturbation strength and other parameters is estimated. This yields a reconnection rate, or inflow Mach number M, scaling as delta B, the ratio of rms magnetic fluctuations to mean magnetic field. This can greatly exceed the Sweet-Parker rate. The rate of heating of the current sheet is also large and can be many orders of magnitude larger than in a Sweet-Parker current sheet. Hyperresistivity produced by tearing mode turbulence may explain both fast reconnection and the heating of the solar corona.

Strauss, H. R.↗

NASA NTRS1953

On the theory of the turbulent boundary layer

As a rule, a division of the turbulent boundary layer is admissible: a division into a part near the wall, where the flow is governed only by the wall effects, and into an outer part, where the wall roughness and the viscosity of the flow medium affects only the wall shearing stress occurring as boundary condition but does not exert any other influence on the flow. Both parts may be investigated to a large extent independently. Under certain presuppositions there result for the outer part "similar" solutions. The theoretical considerations give a cue how to set up, by appropriate experiments and their evaluation, generally valid connections which are required for the approximate calculation of the turbulent boundary layer according to the momentum and energy theorem.

TURBULENCE↗

NASA NTRS1994

A vorticity dynamics based model for the turbulent dissipation: Model development and validation

A new model dissipation rate equation is proposed based on the dynamic equation of the mean-square vorticity fluctuation for large Reynolds number turbulence. The advantage of working with the vorticity fluctuation equation is that the physical meanings of the terms in this equation are more clear than those in the dissipation rate equation. Hence, the model development based on the vorticity fluctuation equation is more straightforward. The resulting form of the model equation is consistent with the spectral energy cascade analysis introduced by Lumley. The proposed model dissipation rate equation is numerically well behaved and can be applied to any level of turbulence modeling. It is applied to a realizable eddy viscosity model. Flows that are examined include: rotating homogeneous shear flows; free shear flows; a channel flow and flat plate boundary layers with and without pressure gradients; and backward facing step separated flows. In most cases, the present model predictions show considerable improvement over the standard kappa-epsilon model.

Shih, Tsan-Hsing↗

NASA NTRS2004

The Effects of Gravity on the Crystallization Behavior of Heavy Metal Fluoride Glasses

Heavy metal fluoride glasses are used in such applications as fiber lasers and laser amplifiers. ZrF4-BaF2-LaF3-AlF3-NaF (ZBLAN) is one of the more commonly used heavy metal fluoride glasses. ZBLAN is an infrared transmitter and has a theoretical attenuation coefficient of 0.002 db/km. However, due to impurities and small crystallites this attenuation coefficient has not been achieved to date. ZBLAN is a fragile glass which can lead to rapid crystallization, if the glass is not cooled rapidly to below the glass transition temperature or if the glass is reheated near the crystallization temperature for any period of time. Studies carried on at Marshall Space Flight Center and the University of Alabama in Huntsville since 1993 have shown that heating ZBLAN glass at the crystallization temperature in reduced gravity results in a suppression of crystallization when compared to ZBLAN processed in unit gravity. These studies utilized NASA's KC-135 aircraft and the Conquest sounding rocket. In the first series of experiments, short lengths of ZBLAN fiber were heated to the crystallization temperature in reduced gravity on board the KC- 135 and the Conquest sounding rocket and compared with fibers heated in unit gravity. The fibers processed in reduced gravity showed no evidence of crystallization when studied with x-ray diffraction and scanning electron microscopy. However, the fibers processed in unit gravity were completely crystallized. Subsequent experiments included heating small pieces of ZBLAN glass at the crystallization temperature while viewing with a video camera to follow the crystallization phenomenon. In this experiment crystallization was observed in reduced gravity, however, it was suppressed when compared to heating in unit gravity. In the most recent experiment on board the KC-135, rapid thermal analysis of ZBLAN was performed. A mechanism to explain the observations has been proposed. This mechanism is based on shear thinning whereby, the glass undergoes a reduction in viscosity in the 10(exp 5) - 10(exp 7) poise range, allowing more rapid diffusion and thus crystallization. It is proposed that this mechanism is suppressed in reduced gravity. An experiment is presently being conducted to test this theory. With increased knowledge of ZBLAN behavior in reduced gravity, three low earth orbit tiber drawing facilities have been designed. One would be suitable for use on the International Space Station, another while aboard the Space Shuttle and the third system is a fully automated facility which would operate independently of the ISS or Shuttle in a free float mode. The primary benefits of free floating a facility in LEO includes a higher quality of microgravity and reduced safety concerns since it is not in a manned environment.

Tucker, Dennis S.↗

NASA NTRSDate not supplied

Assessment of Numerical and Modeling Errors of RANS based Transition Models for Low-Reynolds Numbers 2-D Flows

In this paper we report the outcome of selected workshops organized as part of the NATO Applied Vehicle Technology (AVT)-313 activity Incompressible Laminar-to-Turbulent Flow Transition Study that focused on assessing the numerical and modeling accuracy of the γ−Reθ and γ transition models coupled to the k−ω Shear-Stress Transport (SST) two-equation eddy-viscosity model. Three different test cases involving nominally 2D flow configurations were selected: flow over a flat plate with two different levels of turbulence intensity at the inlet; flow around the Eppler 387 foil at a Reynolds number of 3×10^5 and angles of attack of 1 deg. and 7 deg. flow around the NACA 0015 foil at a Reynolds number of 1.8×10^5 and angles of attack of 5 deg. and10 deg. The flat plate flow conditions correspond to natural and by-pass transition, whereas the other two test cases include laminar separation bubbles that lead to separation-induced transition. For each test case, the selected quantities of interest include both integral and local flow quantities. Geometrically similar grids with a wide range of grid refinement ratios were generated for each of the test cases to allow the estimation of numerical uncertainties for all quantities of interest selected for this study. Several RANS flow solvers were used, employing common grids with the same boundary conditions and mathematical models. Therefore, it is possible to analyze the consistency of the results, i.e., to check if the intervals defined by the different numerical solutions with their respective uncertainties overlap with each other. Modeling errors can also be addressed for the selected flow quantities that have experimental data available. However, the experimental information available in these cases is not sufficient to guarantee that experiments and simulations are performed with the same settings. Nonetheless, the available experimental data is sufficient to guarantee that modeling errors are significantly reduced with the use of the transition models when compared to simulations performed using only the k−ω SST model.

CFD Modeling↗

DOE OSTI2021

Relativistic spin hydrodynamics with torsion and linear response theory for spin relaxation

Using the second law of local thermodynamics and the first-order Palatini formalism, we formulate relativistic spin hydrodynamics for quantum field theories with Dirac fermions, such as QED and QCD, in a torsionful curved background. We work in a regime where spin density, which is assumed to relax much slower than other non-hydrodynamic modes, is treated as an independent degree of freedom in an extended hydrodynamic description. Spin hydrodynamics in our approach contains only three non-hydrodynamic modes corresponding to a spin vector, whose relaxation time is controlled by a new transport coefficient: the rotational viscosity. We study linear response theory and observe an interesting mode mixing phenomenon between the transverse shear and the spin density modes. We propose several field-theoretical ways to compute the spin relaxation time and the rotational viscosity, via the Green-Kubo formula based on retarded correlation functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1993

Experiments on near-wall structure of three-dimensional boundary layers

Investigations of three-dimensional turbulent boundary layers have shown basic differences between two- and three-dimensional flows. These differences can significantly impact the modeling of three-dimensional flows since many flow models are based on results from two-dimensional boundary layers. In many cases the shear stress vector direction has been shown to lag relative to the direction of the mean velocity gradient as the cross flow grows downstream. Coincidence of these vectors is necessary for a scalar eddy viscosity modeling assumption. A second effect is a reduction in magnitude of the shear stress and/or the shear stress to turbulence energy ratio, a(sub 1). This reduction has been observed in several experiments. Recent numerical simulations also indicate wall-layer structural differences between two- and three-dimensional boundary layers. The differences in structure between two- and three-dimensional boundary layers was also addressed in the experiment of Littell & Eaton. The experiment used two-point correlations to investigate the vortical structures in a three dimensional boundary layer on a spinning disk. It was found that each sign of longitudinal vortex is equally likely to exist, but one sign of vorticity is associated with a structure which is better at producing ejections. The goal of the current investigation is to study the structure of the inner layers. Among other questions, the differences between the effects deduced from the three-dimensional flow simulations and the effects seen in experiments can be examined. The research concentrates on the structure of the wall-layer through flow visualization and direct turbulence measurements down to y(+) = 5.

Flack, Karen A.↗

NASA NTRS1993

A Realizable Reynolds Stress Algebraic Equation Model

The invariance theory in continuum mechanics is applied to analyze Reynolds stresses in high Reynolds number turbulent flows. The analysis leads to a turbulent constitutive relation that relates the Reynolds stresses to the mean velocity gradients in a more general form in which the classical isotropic eddy viscosity model is just the linear approximation of the general form. On the basis of realizability analysis, a set of model coefficients are obtained which are functions of the time scale ratios of the turbulence to the mean strain rate and the mean rotation rate. The coefficients will ensure the positivity of each component of the mean rotation rate. These coefficients will ensure the positivity of each component of the turbulent kinetic energy - realizability that most existing turbulence models fail to satisfy. Separated flows over backward-facing step configurations are taken as applications. The calculations are performed with a conservative finite-volume method. Grid-independent and numerical diffusion-free solutions are obtained by using differencing schemes of second-order accuracy on sufficiently fine grids. The calculated results are compared in detail with the experimental data for both mean and turbulent quantities. The comparison shows that the present proposal significantly improves the predictive capability of K-epsilon based two equation models. In addition, the proposed model is able to simulate rotational homogeneous shear flows with large rotation rates which all conventional eddy viscosity models fail to simulate.

Shih, Tsan-Hsing↗

DOE OSTI2021

Direct Numerical Simulation of Flow Over a Wall-Mounted Cube with the Nek5000 Spectral Element Code: DNS at Re h = 3900

The direct numerical simulation (DNS) of the canonical wall-mounted unit cube subjected to two distinct incident velocity profiles is performed at a Reynolds number where the bulk flow characteristics are known to become relatively Reynolds number insensitive. The aim of this work is to highlight the sensitivity of such bluff-body flows to mean shear and to provide a representative set of data for such flow scenarios where common turbulence modeling techniques often fail. Simple boundary conditions and a small domain are selected to allow for cost-effective and easy comparisons for model development purposes. In addition to mean velocity profiles, select turbulence statistics are presented in detail. Further, the basic efficacy of eddy viscosity models is evaluated and found to be adequate only for the shear stress components for the bluff-body flows examined here. Other failure mechanisms for RANS are proposed.

42 ENGINEERING↗

NASA NTRS1984

Shearing stability of lubricants

Shearing stabilities of lubricating oils containing a high mol. wt. polymer as a viscosity index improver were studied by use of ultrasound. The oils were degraded by cavitation and the degradation generally followed first order kinetics with the rate of degradation increasing with the intensity of the ultrasonic irradiation and the cumulative energy applied. The shear stability was mainly affected by the mol. wt. of the polymer additive and could be determined in a short time by mechanical shearing with ultrasound.

Shiba, Y.↗

NASA NTRS1999

Dynamic Nucleation of Supercooled Melts and Measurement of the Surface Tension and Viscosity

We investigate the phenomenon of acoustic pressure-induced nucleation by using a novel approach involving the large amplitude resonant radial oscillations and collapse of a single bubble intentionally injected into a supercooled liquid. Using a combination of previously developed and proven techniques, the bubble is suspended in a fluid host by an ultrasonic field which supplies both the levitation capability as well as the forcing of the radial oscillations. We observe the effects of an increase in pressure (due to bubble collapse) in a region no larger than 100 microns within the supercooled melt to rigorously probe the hypothesis of pressure-induced nucleation of the solid phase. The use of single bubbles operating in narrow temporal and spatial scales will allow the direct and unambiguous correlation between the origin and location of the generation of the disturbance and the location and timing of the nucleation event. In a companion research effort, we are developing novel techniques for the non-contact measurements of the surface tension and viscosity of highly viscous supercooled liquids. Currently used non-invasive methods of surface tension measurement for the case of undercooled liquids generally rely of the quantitative determination of the resonance frequencies of drop shape oscillations, of the dynamics of surface capillary waves, or of the velocity of streaming flows. These methods become quickly ineffective when the liquid viscosity rises to a significant value. An alternate and accurate method which would be applicable to liquids of significant viscosity is therefore needed. We plan to develop such a capability by measuring the equilibrium shape of levitated undercooled melt droplets as they undergo solid-body rotation. The experimental measurement of the characteristic point of transition (bifurcation point) between axisymmetric and two-lobed shapes will be used to calculate the surface tension of the liquid. Such an approach has already been validated through the experimental verification of numerical modeling results. The experimental approach involves levitation, melting, and solidification of undercooled droplets using a hybrid ultrasonic-electrostatic technique in both a gaseous as well as a vacuum environment. A shape relaxation method will be investigated in order to derive a reliable method to measure the viscosity of undercooled melts. The analysis of the monotonic relaxation to equilibrium shape of a drastically deformed and super-critically damped free drop has been used to derive interfacial tension of immiscible liquid combinations where one of the component has high viscosity. A standard approach uses the initial elongation of a droplet through shear flows, but an equivalent method could involve the initial deformation of a drop levitated in a gas by ultrasonic radiation pressure, electric stresses, or even solid body rotation. The dynamic behavior of the free drop relaxing back to equilibrium shape will be modeled, and its characteristic time dependence should provide a quantitative means to evaluate the liquid viscosity.

Trinh, E. H.↗

NASA NTRS1989

Prediction of separated transonic wing flows with a non-equilibrium algebraic model

A nonequilibrium algebraic turbulence model, which is based on the turbulence closure scheme of Johnson and King (1985), is proposed to predict separated transonic wing flows. The influence of history effects are modeled by solving a partial differential equation for the maximum total Reynolds shear stress, which is then used to scale the eddy viscosity of an algebraic model. The turbulence model is implemented in a three-dimensional, Reynolds-averaged Navier-Stokes code. Comparisons with experimental data are presented which show clearly that the nonequilibrium type of turbulence model is essential for accurate prediction of transonic separated flows.

Abid, Ridha↗