Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinematic 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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Modeling the effect of water on mantle rheology

To study the thermal history of the Earth we use a parameterized model of mantle convection. This model includes a mathematical description of de- and regassing processes of water from the Earth's mantle. The rates of this processes are considered to be directly proportional to the seafloor spreading rate. The kinematic viscosity of the mantle depends on the temperature/pressure as well as on the volatile content. Dissolved volatiles such as water weaken the minerals by reducing their activation energy for solid state creep. Karato and Toriumi showed a power law dependence between creep rate and water fugacity derived from experimental results. Therefore, we use such flow parameters of diffusion creep in olivine under wet and dry conditions to calculate the mantle viscosity as a function of the water content. Because the creep rate is proportional to the concentration of water-related point deflects we assume that the water fugacity is proportional to the water weight fraction. An equation for the steady-state strain rate under wet conditions is established. To assess the unknown constant K in this equation, we use flow law parameters given by Karato and Wu as well as the results of McGovern and Schubert.

Bounama, CH.↗
NASA NTRS1996

Experimental study of vortex breakdown in a cylindrical, swirling flow

The stability of a steady, vortical flow in a cylindrical container with one rotating endwall has been experimentally examined to gain insight into the process of vortex breakdowwn. The dynamics of the flow are governed by the Reynolds number (Re) and the aspect ratio of the cylinder. Re is given by Omega R(sup 2)/nu, where Omega is the speed of rotation of the endwall, R is the cylinder radius, and nu is the kinematic viscosity of the fluid filling the cylinder. The aspect ratio is H/R, where H is the height of the cylinder. Numerical simulation studies disagree whether or not the steady breakdown is stable beyond a critical Reynolds number, Re(sub c). Previous experimental researches have considered the steady and unsteady flows near Re(sub c), but have not explored the stability of the steady breakdown structures beyond this value. In this investigation, laser induced fluorescence was utilized to observe both steady and unsteady vortex breakdown at a fixed H/R of 2.5 with Re varying around Re(sub c). When the Re of a steady flow was slowly increased beyond Re(sub c), the breakdown structure remained steady even though unsteadiness was possible. In addition, a number of hysteresis events involving the oscillation periods of the unsteady flow were noted. The results show that both steady and unsteady vortex breakdown occur for a limited range of Re above Re(sub c). Also, with increasing Re, complex flow transformations take place that alter the period at which the unsteady flow oscillates.

Stevens, J. L.↗
NASA NTRS1996

Turbulence Scales, Rise Times, Caustics, and the Simulation of Sonic Boom Propagation

The general topic of atmospheric turbulence effects on sonic boom propagation is addressed with especial emphasis on taking proper and efficient account of the contributions of the portion oi the turbulence that is associated with extremely high wavenumber components. The recent work reported by Bart Lipkens in his doctoral thesis is reexamined to determine whether the good agreement between his measured rise times with the 1971 theory of the author is fortuitous. It is argued that Lipken's estimate of the distance to the first caustic was a gross overestimate because of the use of a sound speed correlation function shaped like a gaussian curve. In particular, it is argued that the expected distance to the first caustic varies with the kinematic viscosity nu and the energy epsilon dissipated per unit mass per unit time, and the sound speed c as : d(sub first caustic) = nu(exp 7/12) c(exp 2/3)/ epsilon(exp 5/12)(nu x epsilon/c(exp 4))(exp a), where the exponent a is greater than -7/12 and can be argued to be either O or 1/24. In any event, the surprising aspect of the relationship is that it actually goes to zero as the viscosity goes to zero with s held constant. It is argued that the apparent overabundance of caustics can be grossly reduced by a general computational and analytical perspective that partitions the turbulence into two parts, divided by a wavenumber k(sub c). Wavenumbers higher than kc correspond to small-scale turbulence, and the associated turbulence can be taken into account by a renormalization of the ambient sound speed so that the result has a small frequency dependence that results from a spatial averaging over of the smaller-scale turbulent fluctuations. Selection of k(sub c). can be made so large that only a very small number of caustics are encountered if one adopts the premise that the frequency dispersion of pulses is caused by that part of the turbulence spectrum which lies in the inertial range originally predicted by Kolmogoroff. The acoustic propagating wave's dispersion relation has the acoustic wavenumber being of the form k = (omega/c) + F(omega), where c is a spatially averaged sound speed and where, for mechanical turbulence, the extra term F(omega) must depend on only the angular frequency omega, the sound speed c, and the turbulent energy dissipation epsilon per unit fluid mass and per unit time. If the turbulence is weak, then the quantity F(omega) has to be of second order in the portions of the turbulent fluid velocity in the inertial range, so, following Kolmogoroff's reasoning, it must vary with epsilon as epsilon(exp 2/3). Simple dimensional analysis then reveals that F(omega) is K epsilon(exp 2/3) c(exp -7/3) omega(exp l/3), K being a universal dimensionless complex constant.

Pierce, Allan D.↗
NASA NTRS1996

Buoyancy-Driven Heat Transfer During Application of a Thermal Gradient for the Study of Vapor Deposition at Low Pressure Using and Ideal Gas

A mathematical model has been developed to determine heat transfer during vapor deposition of source materials under a variety of orientations relative to gravitational accelerations. The model demonstrates that convection can occur at total pressures as low as 10-2 mm Hg. Through numerical computation, using physical material parameters of air, a series of time steps demonstrates the development of flow and temperature profiles during the course of vapor deposition. These computations show that in unit gravity vapor deposition occurs by transport through a fairly complicated circulating flow pattern when applying heat to the bottom of the vessel with parallel orientation with respect to the gravity vector. The model material parameters for air predict the effect of kinematic viscosity to be of the same order as thermal diffusivity, which is the case for Prandtl number approx. 1 fluids. Qualitative agreement between experiment and the model indicates that 6-(2-methyl-4-nitroanilino)-2,4-hexadiyn-l-ol (DAMNA) at these pressures indeed approximates an ideal gas at the experiment temperatures, and may validate the use of air physical constants. It is apparent that complicated nonuniform temperature distribution in the vapor could dramatically affect the homogeneity, orientation, and quality of deposited films. The experimental test i's a qualitative comparison of film thickness using ultraviolet-visible spectroscopy on films generated in appropriately oriented vapor deposition cells. In the case where heating of the reaction vessel occurs from the top, deposition of vapor does not normally occur by convection due to a stable stratified medium. When vapor deposition occurs in vessels heated at the bottom, but oriented relative to the gravity vector between these two extremes, horizontal thermal gradients induce a complex flow pattern. In the plane parallel to the tilt axis, the flow pattern is symmetrical and opposite in direction from that where the vessel is positioned vertically. The ground-based experiments are sufficient preliminary tests of theory and should be of significant interest regarding vapor deposited films in microgravity.

Frazier, D. O.↗
NASA NTRS1998

Oscillatory Thermocapillary Flow Experiment-2 (OTFE-2)

Oscillatory thermocapillary flow experiments were performed in the Glovebox aboard the USML-2 Spacelab which was launched on October 20, 1995. Open cylindrical containers of 1.2 and 2.0 cm in diameter were used. The ratio of container depth to radius (aspect ratio) was set at 0.5 and 2. Silicone oil of 2 centistokes kinematic viscosity was the test fluid. The fluid was heated by a cylindrical heater placed along the center axis of the container. The fluid motion was studied by flow visualization. The effect of aspect ratio on the onset of oscillations and on the oscillatory flow was investigated. It was found that the onset of oscillations was delayed when the container was made more shallow.

Kamotani, Yasuhiro↗
NASA NTRS1998

Experiment 6: Surface Tension Driven Convection Experiment-2 (STDCE-2)

Results are reported from the Second Surface Tension Driven Convection Experiment (STDCE-2) performed aboard the Second United States Microgravity Laboratory (USML-2), which flew as Space Shuttle mission STS-73 and launched on October 20, 1995. Oscillatory thermocapillary flows were investigated in open cylindrical containers filled with 2 centistokes kinematic viscosity (Pr=27 at 25 C) silicone oil. Two different heating modes were investigated in detail. The main objectives of the experiments were to determine the onset of oscillatory thermocapillary flow under highly reduced buoyancy and gravity conditions and to study the important features of the oscillatory flow. The onset conditions were determined in three different size containers (1.2-, 2-, and 3-cm diameter) and for various free surface shapes (flat and curved). Numerical and scaling analyses were also performed to understand the basic steady flows. The analysis shows that the main flow is viscous-dominated near the onset of oscillations. The onset conditions determined in the present experiments together with our earlier ground-based data show that one could not describe the oscillation phenomenon if the fluid free surface is assumed to be rigid. Therefore, a parameter representing free surface deformation was derived, and the data are shown to be correlated well by that parameter. The oscillation patterns and frequencies are also presented.

Ostrach, Simon↗
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 NTRS2003

Suppression of Buoyancy in Gaseous Media at High Temperatures

Consider a rectangular box filled with a fluid having a heated bottom and a cold top surface, and insulated side-walls (Benard problem). As the temperature difference between the horizontal top and bottom surfaces increases, a critical condition, defined quantitatively by the Rayleigh number, is reached beyond which density stratification can no longer be sustained by conduction and the fluid disrupts from its stable, quiescent state into an unstable, convective mode in which lighter and heavier gas mix. This paper suggests that such a statement is not necessarily true for gaseous media under normalized temperature differences that are much larger than justifiable for the Boussinesq approximation! In fact, there may be situations where a system cannot ever be made unstable with respect to the onset on buoyant convection no matter how large the temperature (density) difference becomes at a given pressure even under normal gravity! This unexpected behavior is primarily attributed to highly temperature-sensitive kinematic viscosity which counteracts the tendency toward instability and dampens convection by making the gas more viscous at higher temperatures. This compensation of the buoyant force by the viscous force exhibits itself by the formation of a peak hot-surface temperature beyond which a system will tend to be more stable as the hot-surface temperature increases.

Gokoglu, Suleyman A.↗
NASA NTRS2004

PIV Investigations of the Flow Field in the Volute of a Rotary Blood Pump

A full-size acrylic model of a rotary blood pump was developed in order to utilize Particle Image Velocimetry (PIV) to make measurements of the fluid velocities and turbulent stresses throughout the device. The development of an understanding of the hemodynamics within the blood pump is critical to the development and validation of computational models. A blood analog solution, consisting of sodium iodide solution and glycerin, was developed to match physiological kinematic viscosity. The refractive indecies of the fluid, the pump casing and the impeller were matched to facilitate the use of PIV to make velocity measurements. Velocity measurements made in the volute exit/diffuser region are presented for pumps speeds of 3000-3850 rpm. At each speed data were obtained at a physiological pressure of 90 mmHg and at a maximum flow condition. Four hundred data pairs were used for each resultant mean velocity vector value, representing greater than an order of magnitude more data pairs than reported previously in the literature on similar devices and resulting in velocity uncertainty levels of approximately 2.9%.

Sankovic, John M.↗
NASA NTRS2005

Experimental Study of Large-Amplitude Faraday Waves in Rectangular Cylinders

Experiment on single-mode Faraday waves having two, thee, and four wavelengths across a rectangular cylinder of high aspect ratio is the subject of discussion. Previous experiments recently done by Henderson & Miles (1989) and by Lei Jiang et. a1 (1996) focused on Faraday waves with one and two wavelengths across rectangular cylinders. In this experimental study the waves steepness ranges from small at threshold levels to a large amplitude which according to Penny & Price theory (1952) approaches the maximum sustainable amplitude for a standing wave. The waves characteristics for small amplitudes are evaluated against an existing well known linear theory by Benjamin & Ursell (l954) and against a weakly nonlinear theory by J. Miles (1984) which includes the effect of viscous damping. The evaluation includes the wave neutral stability and damping rate. In addition, a wave amplitude differential equation of a linear theory including viscous effect by Cerda & Tirapegui (1998) is solved numerically to yield prediction of temporal profiles of both wave damping and wave formation at the threshold. An interesting finding from this exercise is that the fluid kinematic viscosity needs to increase ten times in order to obtain good agreement between the theoretical prediction and the experimental data for both wave damping and wave starting. For large amplitude waves, the experimental data are evaluated against the theory of Penny & Price which predicts wave characteristics of any amplitude up to the point at which the wave reaches its maximum amplitude attainable for a standing wave. The theory yields two criteria to show the maximum wave steepness, the vertical acceleration at the wave crest of half the earth gravity field acceleration and the including angle at the crest of 90 degrees. Comparison with experimental data shows close agreement for the wave crest acceleration but a large discrepancy for the including angle. Additional information is included in the original extended abstract.

Iek, Chanthy↗
NASA NTRS2004

PIV Investigations of the Flow Field in the Volute of a Rotary Blood Pump

A full-size acrylic model of a rotary blood pump was developed in order to utilize Particle Image Velocimetry (PIV) to make measurements of the fluid velocities and turbulent stresses throughout the device. The development of an understanding of the hemodynamics within the blood pump is critical to the development and validation of computational models. A blood analog solution, consisting of sodium iodide solution and glycerin, was developed to match physiological kinematic viscosity. The refractive indices of the fluid, the pump casing, and the impeller were matched to facilitate the use of PIV to make velocity measurements. Velocity measurements made in the volute exit/diffuser region are presented for pumps speeds of 3000-3850 rpm. At each speed data were obtained at a physiological pressure of 12 kPa and at a maximum flow condition. Four hundred data pairs were used for each resultant mean velocity vector value, representing greater than an order of magnitude more data pairs than reported previously in the literature on similar devices and resulting in velocity uncertainty levels of approximately 22.9%.

Sankovic, John M.↗
NASA NTRS2016

Origin of the Two Scales of Wind Ripples on Mars

Earth's sandy deserts host two main types of bedforms - decimeter-scale ripples and larger dunes. Years of orbital observations on Mars also confirmed the existence of two modes of active eolian bedforms - meter-scale ripples, and dunes. By analogy to terrestrial ripples, which are thought to form from a grain mechanism, it was hypothesized that large martian ripples also formed from grain impacts, but spaced further apart due to elongated saltation trajectories from the lower martian gravity and different atmospheric properties. However, the Curiosity rover recently documented the coexistence of three scales of bedforms in Gale crater. Because a grain impact mechanism cannot readily explain two distinct and coeval ripple modes in similar sand sizes, a new mechanism seems to be required to explain one of the scales of ripples. Small ripples are most similar to Earth's impact ripples, with straight crests and subdued profiles. In contrast, large martian ripples are sinuous and asymmetric, with lee slopes dominated by grain flows and grainfall deposits. Thus, large martian ripples resemble current ripples formed underwater on Earth, suggesting that they may form from a fluid-drag mechanism. To test this hypothesis, we develop a scaling relation to predict the spacing of fluid-drag ripples from an extensive flume data compilation. The size of large martian ripples is predicted by our scaling relation when adjusted for martian atmospheric properties. Specifically, we propose that the wavelength of martian wind-drag ripples arises from the high kinematic viscosity of the low-density atmosphere. Because fluid density controls drag-ripple size, our scaling relation can help constrain paleoatmospheric density from wind-drag ripple stratification.

Lapotre, Mathieu G. A.↗
NASA NTRS2015

Comparison of Cryogenic Flow Boiling in Liquid Nitrogen and Liquid Hydrogen Chilldown Experiments

This paper presents a comparison between experimental results from recent liquid hydrogen (LH 2 ) transfer line chilldown experiments at high Reynolds (Re) numbers versus liquid nitrogen (LN 2 ) experiments conducted at low Re numbers. Parasitic heat leak, inner wall temperatures, inner wall heat fluxes, and heat transfer coefficients are computed to compare between the two systems. Analysis of temperature traces and flow visualization indicates that the chilldown process evolves much more rapidly at higher Re numbers due to a quick transition from vapor flow to annular liquid flow and near immediate liquid contact along the pipe walls. The lower kinematic viscosity and surface tension of LH 2 , along with reduced parasitic heat leak and higher Re numbers relative to the LN 2 experiments, causes chilldown to proceed almost immediately into the nucleate boiling regime, in comparison to low Re flows where >75% of the chilldown is spent in vapor film boiling.

Cryogenic Heat Transfer Coefficient↗
NASA NTRS1989

Navier-Stokes solutions about the F/A-18 forebody-LEX configuration

Three-dimensional viscous How computations are presented for the F/A-18 fore body-LEX geometry. Solutions are obtained from an algorithm for the compressible Navier-Stokes equations which incorporates an upwind-biased, flux-difference-splitting approach along with longitudinally-patched grids. Results are presented for both laminar and fully turbulent flow assumptions and include correlations with wind tunnel as well as flight-test results. A good quantitative agreement for the forebody surface pressure distribution is achieved between the turbulent computations and wind tunnel measurements at M ∞ = 0.6. The computed turbulent surface flow patterns on the forebody qualitatively agree well with in-flight surface flow patterns obtained on an F/A-18 aircraft at M ∞ = 0.34.

Transitional flow↗
NASA NTRS1993

Advancements in the LEWICE Ice Accretion Model

Recent evidence has shown that the NASA/Lewis Ice Accretion Model, LEWICE, does not predict accurate ice shapes for certain glaze ice conditions. This paper will present the methodology used to make a first attempt at improving the ice accretion prediction in these regimes. Importance is given to the correlations for heat transfer coefficient and ice density, as well as runback flow, selection of the transition point, flow field resolution, and droplet trajectory models. Further improvements and refinement of these modules will be performed once tests in NASA's Icing Research Tunnel, scheduled for 1993, are completed.

Icing↗
NASA NTRS1989

Navier-Stokes Solutions about the F/A-18 Forebody-LEX Configuration

Three-dimensional viscous flow computations are presented for the F/A-18 forebody-LEX geometry. Solutions are obtained from an algorithm for the compressible Navier-Stokes equations which incorporates an upwind-biased, flux-difference-splitting approach along with longitudinally-patched grids. Results are presented for both laminar and fully turbulent flow assumptions and include correlations with wind tunnel as well as flight-test results. A good quantitative agreement for the forebody surface pressure distribution is achieved between the turbulent computations and wind tunnel measurements at Mach number of 0.6. The computed turbulent surface flow patterns on the forebody qualitatively agree well with in-flight surface flow patterns obtained on an F/A-fS aircraft at Mach number of 0.34.

Transitional flow↗
NASA NTRS1976

Heating of the solar chromosphere and corona. I - Generalized inhomogeneous wave equation for magnetoacoustic motions

The generalized inhomogeneous wave equation that governs magnetoacoustic, vortical, and thermal motions in compressible fluids and that is applicable to the problem of heating of the solar chromosphere and corona is obtained. The effects of kinematic and bulk viscosity, heat conduction, Joule dissipation, and magnetic diffusivity are included. Under the usual assumptions, the generalized wave equation reduces to the well-known equations of Lighthill, Kulsrud, Phillips, and others. The major problems encountered in applying Lighthill's (1952) mechanism to sound generation in turbulent media are reviewed for both the subsonic and supersonic cases.

Anand, S. P. S.↗
NASA NTRS1974

A finite element solution algorithm for the Navier-Stokes equations

A finite element solution algorithm is established for the two-dimensional Navier-Stokes equations governing the steady-state kinematics and thermodynamics of a variable viscosity, compressible multiple-species fluid. For an incompressible fluid, the motion may be transient as well. The primitive dependent variables are replaced by a vorticity-streamfunction description valid in domains spanned by rectangular, cylindrical and spherical coordinate systems. Use of derived variables provides a uniformly elliptic partial differential equation description for the Navier-Stokes system, and for which the finite element algorithm is established. Explicit non-linearity is accepted by the theory, since no psuedo-variational principles are employed, and there is no requirement for either computational mesh or solution domain closure regularity. Boundary condition constraints on the normal flux and tangential distribution of all computational variables, as well as velocity, are routinely piecewise enforceable on domain closure segments arbitrarily oriented with respect to a global reference frame.

Baker, A. J.↗