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 433 records · Page 24

Laminar Boundary Layer Behind a Strong Shock Moving into Air

The laminar wall boundary layer behind a strong shock advancing into stationary air has been determined. Numerical results have been obtained for shock Mach numbers up to 14 using real gas values for density and viscosity and assuming Prandtl and Lewis numbers of 0.72 and 1, respectively. The numerical results for shear and heat transfer agree, within 4 percent, with a previously presented approximate analytical expression for these quantities. A slight modification of this expression results in agreement with the numerical data to within 2.5 percent. Analytical expressions for boundary-layer thickness and displacement thickness, correct to within 4 percent for the present data, have also been obtained.

Mirels, Harold↗

Integrated High Payoff Rocket Propulsion Technology (IHPRPT) SiC Recession Model

SiC stability and recession rates were modeled in hydrogen/oxygen combustion environments for the Integrated High Payoff Rocket Propulsion Technology (IHPRPT) program. The IHPRPT program is a government and industry program to improve U.S. rocket propulsion systems. Within this program SiC-based ceramic matrix composites are being considered for transpiration cooled injector faceplates or rocket engine thrust chamber liners. Material testing under conditions representative of these environments was conducted at the NASA Glenn Research Center, Cell 22. For the study described herein, SiC degradation was modeled under these Cell 22 test conditions for comparison to actual test results: molar mixture ratio, MR (O2:H2) = 6, material temperatures to 1700 C, combustion gas pressures between 0.34 and 2.10 atm, and gas velocities between 8,000 and 12,000 fps. Recession was calculated assuming rates were controlled by volatility of thermally grown silica limited by gas boundary layer transport. Assumptions for use of this model were explored, including the presence of silica on the SiC surface, laminar gas boundary layer limited volatility, and accuracy of thermochemical data for volatile Si-O-H species. Recession rates were calculated as a function of temperature. It was found that at 1700 C, the highest temperature considered, the calculated recession rates were negligible, about 200 m/h, relative to the expected lifetime of the material. Results compared favorably to testing observations. Other mechanisms contributing to SiC recession are briefly described including consumption of underlying carbon and pitting. A simple expression for liquid flow on the material surface was developed from a one-dimensional treatment of the Navier-Stokes Equation. This relationship is useful to determine under which conditions glassy coatings or thermally grown silica would flow on the material surface, removing protective layers by shear forces. The velocity of liquid flow was found to depend on the gas velocity, the viscosity of gas and liquid, as well as the thickness of the gas boundary layer and the liquid layer. Calculated flow rates of a borosilicate glass coating compared well to flow rates observed for this coating tested on a SiC panel in Cell 22.

Opila, E. J.↗

Model-Invariant Hybrid Computations of Separated Flows for RCA Standard Test Cases

NASA's Revolutionary Computational Aerosciences (RCA) subproject has identified several smooth-body separated flows as standard test cases to emphasize the challenge these flows present for computational methods and their importance to the aerospace community. Results of computations of two of these test cases, the NASA hump and the FAITH experiment, are presented. The computations were performed with the model-invariant hybrid LES-RANS formulation, implemented in the NASA code VULCAN-CFD. The model- invariant formulation employs gradual LES-RANS transitions and compensation for model variation to provide more accurate and efficient hybrid computations. Comparisons revealed that the LES-RANS transitions employed in these computations were sufficiently gradual that the compensating terms were unnecessary. Agreement with experiment was achieved only after reducing the turbulent viscosity to mitigate the effect of numerical dissipation. The stream-wise evolution of peak Reynolds shear stress was employed as a measure of turbulence dynamics in separated flows useful for evaluating computations.

Woodruff, Stephen↗

A Rotary Flow Channel for Shear Stress Sensor Calibration

A proposed shear sensor calibrator consists of a rotating wheel with the sensor mounted tangential to the rim and positioned in close proximity to the rim. The shear stress generated by the flow at the sensor position is simply tau(sub omega) = (mu)r(omega)/h, where mu is the viscosity of the ambient gas, r the wheel radius, omega the angular velocity of the wheel, and h the width of the gap between the wheel rim and the sensor. With numerical values of mu = 31 (mu)Pa s (neon at room temperature), r = 0.5 m, omega = 754 /s (7200 rpm), and h = 50.8 m, a shear stress of tau(sub omega) = 231 Pa can be generated. An analysis based on one-dimensional flow, with the flow velocity having only an angular component as a function of the axial and radial coordinates, yields corrections to the above simple formula for the curvature of the wheel, flatness of the sensor, and finite width of the wheel. It is assumed that the sensor mount contains a trough (sidewalls) to render a velocity release boundary condition at the edges of the rim. The Taylor number under maximum flow conditions is found to be 62.3, sufficiently low to obviate flow instability. The fact that the parameters entering into the evaluation of the shear stress can be measured to high accuracy with well-defined uncertainties makes the proposed calibrator suitable for a physical standard for shear stress calibration.

Zuckerwar, Allan J.↗

Avoiding surface instability and slurry jamming in simultaneous multilayer coating of structured particulate films

Simultaneous multilayer coating techniques are widely known, but their industrial application remains limited to narrow market sectors. One barrier to adoption may be the mismatch between industries that are familiar with such processes but have no need, and industries that have need but are not familiar. Moreover, there are application-specific technical challenges to developing multilayer-coating processes. In this article, we describe our resolution of two specific issues in simultaneous multilayer coating of all aqueous highly loaded slurries for new and emerging energy applications. The first issue is particle jamming (associated with Shear-thickening) of the highly loaded slurries in the die internals, which we alleviated by adding small amounts of viscosity modifiers without reducing the solid loading. The second issue is a Marangoni-driven surface instability that resembles top layer de-wetting, which we solved by carefully selecting surfactants to tune the dynamic surface tensions of each slurry. Both issues were resolved early in a step-wise development, saving significant development cost which in our case was driven by expensive materials.

08 HYDROGEN↗

Correlating inter-particle forces and particle shape to shear-induced aggregation/fragmentation and rheology for dilute anisotropic particle suspensions: A complementary study via capillary rheometry and in-situ small and ultra-small angle X-ray scattering

Predicting the stability of suspensions composed of anisotropic particles is difficult due to the complexity and interplay of hydrodynamic and colloidal forces. For example, it is not well understood how non-spherical particle shapes affect suspension-rheological responses. This study presents evidence showing how the stacking of boehmite platelets dispersed in aqueous solution leads to moderate shear thinning, even at low boehmite mass fractions. Wide angle X-ray, small angle X-ray, and ultra-small angle X-ray scattering were used to simultaneously monitor changes to the size and fractal dimensions of boehmite aggregates from 6 – 10,000 Å as the sample was recirculated through an in situ capillary rheometer. The latter also provided simultaneous measurements of the viscosity of the suspensions. Computational fluid dynamics modeling of the apparatus was used to provide a more rigorous analysis of the fluid flow. In addition to shear-induced aggregation/fragmentation correlated to a complicated balance between hydrodynamic and colloidal forces, the results indicate that orientationally-dependent interactions, arising from non-spherical boehmite aggregate sub-units, may explain the relatively large experimental viscosity when the hydrodynamic force is small compared to colloidal forces. Ultimately, the results from this study can provide a foundation for future work investigating particle shape effects on colloidal and hydrodynamics forces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pressure-viscosity measurements for several lubricants to 550 meganewtons per square meter /80,000 PSI/ and 149 C /300 F/

The viscosities of a number of liquid lubricants and lubricant formulations, determined as function of pressure, temperature, and shear stress by means of a high-pressure capillary viscometer, are reviewed. Where possible, these results are compared with those obtained by other techniques (optical elastohydrodynamics, oscillating crystal, and low shear capillary viscometry).

Jones, W. R., Jr.↗

Calculations of transonic potential flow over cascades

Transonic flow through a cascade was studied by using the full potential equation and the finite volume method of Jameson and Caughey. The C-type computational grid is generated by an electrostatic analogy and simple shearing transformation. The solution algorithm includes an option of using either an artificial density or an artificial viscosity formulation of the dissipative term. Using the developed code, flows through a cascade of NACA 0012 airfoils and flows through a cascade of shockless blades were computed. It is found that the designed flow through the shockless blade is accurately predicted, the artificial density formulation shows more tolerance to the mesh irregularity, and the C-type mesh does not extend very far upstream for a small pitch-cord ratio.

Jou, W. H.↗

The isostatic state of the lunar Apennines and regional surroundings

High-resolution gravity and topography data taken over the Apennine Mountains have been used to compute their isostatic state. Results show that the Apennines are uncompensated; thus this state implies that the lunar crust and upper mantle have been strong enough over 3.9 b.y. to support the load exerted by this topographic excess. The Apennines produce a maximum shear stress of 60 bars at a depth of 60 km. A lower bound on the lunar crustal viscosity of 10 to the 27th power P is calculated on the basis of the assumption of a 10% relaxation over 3.9 b.y. Studies of a broad negative regional anomaly located between Maria Serenitatis and Imbrium necessitate a locally thicker crust to satisfy the observed data. This anomaly may have been produced by a lateral transport of crustal material from beneath the giant impact basins as a result of rebound at the crust-mantle interface.

Ferrari, A. J.↗

The existence of Hadley convective regimes of atmospheric motion

The solutions of the equations describing deep global convection on a rotating planet are discussed. The existence of generalized steady axisymmetric solutions is established. It is then shown that these are classical solutions when the heat source is sufficiently smooth. The solutions are shown to be unique when the heating is sufficiently weak and asymptotically stable when the shear is sufficiently small. Finally, the application of these results to earth's atmosphere is discussed, with eddy viscosity replacing molecular viscosity.

Dutton, J. A.↗

Local magnetohydrodynamic instabilities and the wave-driven dynamo in accretion disks

We consider the consequences of magnetic buoyancy and the magnetic shearing instability (MSI) on the strength and organization of the magnetic field in a thin accretion disk. We discuss a model in which the wave-driven dynamo growth rate is balanced by the dissipative effects of the MSI. As in earlier work, the net helicity is due to small advective motions driven by nonlinear interactions between internal waves. Assuming a simple model of the internal wave spectrum generated from the primary m = 1 internal waves, we find that the magnetic energy density saturates at about (H/r) exp 4/3 times the local pressure (where H is the disk thickness and r is its radius). On very small scales the shearing instability will produce an isotropic fluctuating field. For a stationary disk this is equivalent to a dimensionless 'viscosity' of about (H/r) exp 4/3. The vertical and radial diffusion coefficients will be comparable to each other. Magnetic buoyancy will be largely suppressed by the turbulence due to the MSI. We present a rough estimate of its effects and find that it removes magnetic flux from the disk at a rate comparable to that caused by turbulent diffusion.

Vishniac, Ethan T.↗

Studies of Two-Phase Gas-Liquid Flow in Microgravity

Two-phase gas-liquid flows are expected to occur in many future space operations. Due to a lack of buoyancy in the microgravity environment, two-phase flows are known to behave differently than those in earth gravity. Despite these concerns, little research has been conducted on microgravity two-phase flow and the current understanding is poor. This dissertation describes an experimental and modeling study of the characteristics of two-phase flows in microgravity. An experiment was operated onboard NASA aircraft capable of producing short periods of microgravity. In addition to high speed photographs of the flows, electronic measurements of void fraction, liquid film thickness, bubble and wave velocity, pressure drop and wall shear stress were made for a wide range of liquid and gas flow rates. The effects of liquid viscosity, surface tension and tube diameter on the behavior of these flows were also assessed. From the data collected, maps showing the occurrence of various flow patterns as a function of gas and liquid flow rates were constructed. Earth gravity two-phase flow models were compared to the results of the microgravity experiments and in some cases modified. Models were developed to predict the transitions on the flow pattern maps. Three flow patterns, bubble, slug and annular flow, were observed in microgravity. These patterns were found to occur in distinct regions of the gas-liquid flow rate parameter space. The effect of liquid viscosity, surface tension and tube diameter on the location of the boundaries of these regions was small. Void fraction and Weber number transition criteria both produced reasonable transition models. Void fraction and bubble velocity for bubble and slug flows were found to be well described by the Drift-Flux model used to describe such flows in earth gravity. Pressure drop modeling by the homogeneous flow model was inconclusive for bubble and slug flows. Annular flows were found to be complex systems of ring-like waves and a substrate film. Pressure drop was best fitted with the Lockhart- Martinelli model. Force balances suggest that droplet entrainment may be a large component of the total pressure drop.

Bousman, William Scott↗

Heat Transfer on a Flat Plate with Uniform and Step Temperature Distributions

Heat transfer associated with turbulent flow on a step-heated or cooled section of a flat plate at zero angle of attack with an insulated starting section was computationally modeled using the GASP Navier-Stokes code. The algebraic eddy viscosity model of Baldwin-Lomax and the turbulent two-equation models, the K- model and the Shear Stress Turbulent model (SST), were employed. The variations from uniformity of the imposed experimental temperature profile were incorporated in the computations. The computations yielded satisfactory agreement with the experimental results for all three models. The Baldwin- Lomax model showed the closest agreement in heat transfer, whereas the SST model was higher and the K-omega model was yet higher than the experiments. In addition to the step temperature distribution case, computations were also carried out for a uniformly heated or cooled plate. The SST model showed the closest agreement with the Von Karman analogy, whereas the K-omega model was higher and the Baldwin-Lomax was lower.

Bahrami, Parviz A.↗

Navier-Stokes Computations With One-Equation Turbulence Model for Flows Along Concave Wall Surfaces

This report presents the use of a time-marching three-dimensional compressible Navier-Stokes equation numerical solver with a one-equation turbulence model to simulate the flow fields developed along concave wall surfaces without and with a downstream extension flat wall surface. The 3-D Navier- Stokes numerical solver came from the NASA Glenn-HT code. The one-equation turbulence model was derived from the Spalart and Allmaras model. The computational approach was first calibrated with the computations of the velocity and Reynolds shear stress profiles of a steady flat plate boundary layer flow. The computational approach was then used to simulate developing boundary layer flows along concave wall surfaces without and with a downstream extension wall. The author investigated the computational results of surface friction factors, near surface velocity components, near wall temperatures, and a turbulent shear stress component in terms of turbulence modeling, computational mesh configurations, inlet turbulence level, and time iteration step. The computational results were compared with existing measurements of skin friction factors, velocity components, and shear stresses of the developing boundary layer flows. With a fine computational mesh and a one-equation model, the computational approach could predict accurately the skin friction factors, near surface velocity and temperature, and shear stress within the flows. The computed velocity components and shear stresses also showed the vortices effect on the velocity variations over a concave wall. The computed eddy viscosities at the near wall locations were also compared with the results from a two equation turbulence modeling technique. The inlet turbulence length scale was found to have little effect on the eddy viscosities at locations near the concave wall surface. The eddy viscosities, from the one-equation and two-equation modeling, were comparable at most stream-wise stations. The present one-equation turbulence model is an effective approach for turbulence modeling in the near solid wall surface region of flow over a concave wall.

Wang, Chi R.↗

On the Asymptotic Regimes and the Strongly Stratified Limit of Rotating Boussinesq Equations

Asymptotic regimes of geophysical dynamics are described for different Burger number limits. Rotating Boussinesq equations are analyzed in the asymptotic limit, of strong stratification in the Burger number of order one situation as well as in the asymptotic regime of strong stratification and weak rotation. It is shown that in both regimes horizontally averaged buoyancy variable is an adiabatic invariant for the full Boussinesq system. Spectral phase shift corrections to the buoyancy time scale associated with vertical shearing of this invariant are deduced. Statistical dephasing effects induced by turbulent processes on inertial-gravity waves are evidenced. The 'split' of the energy transfer of the vortical and the wave components is established in the Craya-Herring cyclic basis. As the Burger number increases from zero to infinity, we demonstrate gradual unfreezing of energy cascades for ageostrophic dynamics. The energy spectrum and the anisotropic spectral eddy viscosity are deduced with an explicit dependence on the anisotropic rotation/stratification time scale which depends on the vertical aspect ratio parameter. Intermediate asymptotic regime corresponding to strong stratification and weak rotation is analyzed where the effects of weak rotation are accounted for by an asymptotic expansion with full control (saturation) of vertical shearing. The regularizing effect of weak rotation differs from regularizations based on vertical viscosity. Two scalar prognostic equations for ageostrophic components (divergent velocity potential and geostrophic departure ) are obtained.

Babin, A.↗

Boundary-layer computational model for predicting the flow and heat transfer in sudden expansions

Fully developed turbulent and laminar flows through symmetric planar and axisymmetric expansions with heat transfer were modeled using a finite-difference discretization of the boundary-layer equations. By using the boundary-layer equations to model separated flow in place of the Navier-Stokes equations, computational effort was reduced permitting turbulence modelling studies to be economically carried out. For laminar flow, the reattachment length was well predicted for Reynolds numbers as low as 20 and the details of the trapped eddy were well predicted for Reynolds numbers above 200. For turbulent flows, the Boussinesq assumption was used to express the Reynolds stresses in terms of a turbulent viscosity. Near-wall algebraic turbulence models based on Prandtl's-mixing-length model and the maximum Reynolds shear stress were compared.

Lewis, J. P.↗

Improved two-equation k-omega turbulence models for aerodynamic flows

Two new versions of the k-omega two-equation turbulence model will be presented. The new Baseline (BSL) model is designed to give results similar to those of the original k-omega model of Wilcox, but without its strong dependency on arbitrary freestream values. The BSL model is identical to the Wilcox model in the inner 50 percent of the boundary-layer but changes gradually to the high Reynolds number Jones-Launder k-epsilon model (in a k-omega formulation) towards the boundary-layer edge. The new model is also virtually identical to the Jones-Lauder model for free shear layers. The second version of the model is called Shear-Stress Transport (SST) model. It is based on the BSL model, but has the additional ability to account for the transport of the principal shear stress in adverse pressure gradient boundary-layers. The model is based on Bradshaw's assumption that the principal shear stress is proportional to the turbulent kinetic energy, which is introduced into the definition of the eddy-viscosity. Both models are tested for a large number of different flowfields. The results of the BSL model are similar to those of the original k-omega model, but without the undesirable freestream dependency. The predictions of the SST model are also independent of the freestream values and show excellent agreement with experimental data for adverse pressure gradient boundary-layer flows.

Menter, Florian R.↗

Characterization of annular two-phase gas-liquid flows in microgravity

A series of two-phase gas-liquid flow experiments were developed to study annular flows in microgravity using the NASA Lewis Learjet. A test section was built to measure the liquid film thickness around the perimeter of the tube permitting the three dimensional nature of the gas-liquid interface to be observed. A second test section was used to measure the film thickness, pressure drop and wall shear stress in annular microgravity two-phase flows. Three liquids were studied to determine the effects of liquid viscosity and surface tension. The result of this study provide insight into the wave characteristics, pressure drop and droplet entrainment in microgravity annular flows.

Bousman, W. Scott↗