Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pressure gradient”

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 1,261 records · Page 70

Compressible Laminar Boundary Layer over a Yawed Infinite Cylinder with Heat Transfer and Arbitrary Prandtl Number

The equations are presented for the development of the compressible laminar boundary layer over a yawed infinite cylinder. For compressible flow with a pressure gradient the chordwise and spanwise flows are not independent. Using the Stewartson transformation and a linear viscosity-temperature relation yields a set of three simultaneous ordinary differential equations in a form yielding similar solutions. These equations are solved for stagnation-line flow for surface temperatures from zero to twice the free-stream stagnation temperature and for a wide range of yaw angle and free-stream Mach number. The results indicate that the effect of yaw on the heat-transfer coefficient at the stagnation line depends markedly on the free-stream Mach number. An unusual result of the solutions is that for large yaw angles and stream Mach numbers the chordwise velocity within the boundary layer exceeds the local external chordwise velocity, even for a highly cooled wall.

Reshotko, Eli↗

Investigation at supersonic speeds of the variation with Reynolds number and Mach number of the total, base, and skin-friction drag of seven boattail bodies of revolution designed for minimum wave drag

An investigation has been conducted in the Langley 9-inch supersonic tunnel to determine the effect of the variation with Reynolds number and Mach number of the total, base, and skin-friction drag at zero lift of seven boattail bodies of revolution designed for minimum wave drag according to NACA TN 2550. The tests covered a Reynolds number range from approximately 1.0 x 10(6) to 10.0 x 10(6) at Mach numbers of 1.62, 1.93, and 2.41, respectively. The results show that base drag and, in general, the total drag increase with increasing values of the ratio of base area to maximum area B/S(sub max), although the results reported in NACA TN 3054 showed that the wave drag decreased with increasing values of B/S(sub max). The laminar skin-friction drag is in agreement with the theoretical predictions used, and, within the Mach number range of these tests, the simple Blasius incompressible theory gives a satisfactory prediction. Except for values of B/S (sub max) near 1, the Reynolds number of transition increases with increasing Mach number and as this ratio approaches 1 this variation is seen to reverse. These variations in Reynolds number of transition with Mach number appear to be associated with the changes in pressure gradient over the rear of the bodies.

Bromm, August F , Jr↗

Turbulent boundary layer of an airfoil

A need has arisen for a new determination of the velocity profiles in the boundary layer. Assuming that the character of the velocity distribution depends to a large extent on the character of the shear distribution across the boundary layer, we shall consider the nature of the shear distribution for a boundary layer with a pressure gradient.

Fediaevsky, K↗

A Laboratory Model of a Hydrogen/Oxygen Engine for Combustion and Nozzle Studies

A small laboratory diagnostic thruster was developed to augment present low thrust chemical rocket optical and heat flux diagnostics at the NASA Lewis Research Center. The objective of this work was to evaluate approaches for the use of temperature and pressure sensors for the investigation of low thrust rocket flow fields. The nominal engine thrust was 110 N. Tests were performed at chamber pressures of about 255 kPa, 370 kPa, and 500 kPa with oxidizer to fuel mixture ratios between 4.0 and 8.0. Two gaseous hydrogen/gaseous oxygen injector designs were tested with 60 percent and 75 percent fuel film cooling. The thruster and instrumentation designs were proven to be effective via hot fire testing. The thruster diagnostics provided inner wall temperature and static pressure measurements which were compared to the thruster global performance data. For several operating conditions, the performance data exhibited unexpected trends which were correlated with changes in the axial wall temperature distribution. Azimuthal temperature distributions were found to be a function of operating conditions and hardware configuration. The static pressure profiles showed that no severe pressure gradients were present in the rocket. The results indicated that small differences in injector design can result in dramatically different thruster performance and wall temperature behavior, but that these injector effects may be overshadowed by operating at a high fuel film cooling rate.

Morren, Sybil Huang↗

Large eddy simulations of time-dependent and buoyancy-driven channel flows

The primary goal of this work has been to assess the performance of the dynamic SGS model in the large eddy simulation (LES) of channel flows in a variety of situations, viz., in temporal development of channel flow turned by a transverse pressure gradient and especially in buoyancy-driven turbulent flows such as Rayleigh-Benard and internally heated channel convection. For buoyancy-driven flows, there are additional buoyant terms that are possible in the base models, and one objective has been to determine if the dynamic SGS model results are sensitive to such terms. The ultimate goal is to determine the minimal base model needed in the dynamic SGS model to provide accurate results in flows with more complicated physical features. In addition, a program of direct numerical simulation (DNS) of fully compressible channel convection has been undertaken to determine stratification and compressibility effects. These simulations are intended to provide a comparative base for performing the LES of compressible (or highly stratified, pseudo-compressible) convection at high Reynolds number in the future.

Cabot, William H.↗

Progress in modeling hypersonic turbulent boundary layers

A good knowledge of the turbulence structure, wall heat transfer, and friction in turbulent boundary layers (TBL) at high speeds is required for the design of hypersonic air breathing airplanes and reentry space vehicles. This work reports on recent progress in the modeling of high speed TBL flows. The specific research goal described here is the development of a second order closure model for zero pressure gradient TBL's for the range of Mach numbers up to hypersonic speeds with arbitrary wall cooling requirements.

Zeman, Otto↗

A transport equation for eddy viscosity

A transport equation for eddy viscosity is proposed for wall bounded turbulent flows. The proposed model reduces to a quasi-homogeneous form far from surfaces. Near to a surface, the nonhomogeneous effect of the wall is modeled by an elliptic relaxation model. All the model terms are expressed in local variables and are coordinate independent; the model is intended to be used in complex flows. Turbulent channel flow and turbulent boundary layer flows with/without pressure gradient are calculated using the present model. Comparisons between model calculations and direct numerical simulation or experimental data show good agreement.

Durbin, P. A.↗

Computation of confined coflow jets with three turbulence models

A numerical study of confined jets in a cylindrical duct is carried out to examine the performance of two recently proposed turbulence models: an RNG-based K-epsilon model and a realizable Reynolds stress algebraic equation model. The former is of the same form as the standard K-epsilon model but has different model coefficients. The latter uses an explicit quadratic stress-strain relationship to model the turbulent stresses and is capable of ensuring the positivity of each turbulent normal stress. The flow considered involves recirculation with unfixed separation and reattachment points and severe adverse pressure gradients, thereby providing a valuable test of the predictive capability of the models for complex flows. Calculations are performed with a finite-volume procedure. Numerical credibility of the solutions is ensured by using second-order accurate differencing schemes and sufficiently fine grids. Calculations with the standard K-epsilon model are also made for comparison. Detailed comparisons with experiments show that the realizable Reynolds stress algebraic equation model consistently works better than does the standard K-epsilon model in capturing the essential flow features, while the RNG-based K-epsilon model does not seem to give improvements over the standard K-epsilon model under the flow conditions considered.

Zhu, J.↗

Research on two equation models

The k-epsilon model is the most widely used turbulence model in engineering calculations. However, the model has several deficiencies that need to be fixed. This document presents improvements to the capabilities of the k-epsilon model in the following areas: a Galilean and tensorial invariant k-epsilon model for near wall turbulence; a new set of wall functions for attached flows; a new model equation for the dissipation rate, which has a better theoretical basis, contains the contribution of flow inhomogeneity, and captures the effect of the pressure gradient accurately; and a better model for bypass transition due to freestream turbulence.

Yang, Z.↗

Moist wind relationships

The equations describing the temporal and spatial behavior of the kinematic moisture and heat flux are described. In these nonlinear equations, the contribution by diabatic processes to the large-scale flux is composed of two parts. One part is associated with a Rayleigh damping term while the other arises from temporal and spatial changes in the pressure gradient term. It was found that the influence of diabatic processes on large-scale moisture fluxes depends greatly on the degree of balance between forcing and damping terms in the governing kinematic flux equations. The existence of a near balance requires a reduction in the large-scale horizontal geostrophic wind speed. Based on an examination of the moisture flux equations, it is argued that reductions in the large-scale horizontal wind speed observed within major cumulus cloud systems help conserve large-scale moisture fluxes. The deviation of the wind from geostrophic conditions is easily estimated for a near balanced state. This wind modification induces secondary vertical circulations that contribute to convergence, creating or supporting long-lived mesoscale flows. We believe this process to be a major supporter of the mesoscale circulations observed in severe storms and squall lines. In the tropics the wind modification has an antitriptic relationship. These diagnostic findings suggest possible modifications to the wind field in the application of a cumulus parameterization, and may be important in diabatic initialization of numerical weather prediction models.

Raymond, William H.↗

Modeling of the coupled magnetospheric and neutral wind dynamos

The solar wind interaction with the earth's magnetosphere generates electric fields and currents that flow from the magnetosphere to the ionosphere at high latitudes. Consequently, the neutral atmosphere is subject to the dissipation and conversion of this electrical energy to thermal and mechanical energy through Joule heating and Lorentz forcing. As a result of the mechanical energy stored within the neutral wind (caused in part by Lorentz--and pressure gradient--forces set up by the magnetospheric flux of electrical energy), electric currents and fields can be generated in the ionosphere through the neutral wind dynamo mechanism. At high latitudes this source of electrical energy has been largely ignored in past studies, owing to the assumed dominance of the solar wind/magnetospheric dynamo as an electrical energy source to the ionosphere. However, other researchers have demonstrated that the available electrical energy provided by the neutral wind is significant at high latitudes, particularly in the midnight sector of the polar cap and in the region of the magnetospheric convection reversal. As a result, the conclusions of a number of broad ranging high-latitude investigations may be modified if the neutral-wind contribution to high-latitude electrodynamics is properly accounted for. These include the following: studies assessing solar wind-magnetospheric coupling by comparing the cross polar cap potential with solar wind parameters; research based on the alignment of particle precipitation with convection or field aligned current boundaries; and synoptic investigations attributing seasonal variations in the observed electric field and current patterns to external sources. These research topics have been initiated by satellite and ground-based observations and have been attributed to magnetospheric causes. However, the contribution of the neutral wind to the high-latitude electric field and current systems and their seasonal and local time dependence has yet to be quantitatively evaluated. In this program, we are evaluating the coupled magnetospheric and neutral wind dynamos at high latitudes under various conditions. In addition to examining the impact of seasonal variations, we are investigating the consequences of the separate dynamos having pure current-source or voltage-source behaviors.

Thayer, Jeff P.↗

The use of transition region characteristics to improve the numerical simulation of heat transfer in bypass transitional flows

A method is presented for improving the numerical prediction of bypass transition heat transfer on a flat plate in a high-disturbance environment with zero or favorable pressure gradient. The method utilizes low Reynolds number k-epsilon turbulence models in combination with the characteristic parameters of the transition region. The parameters representing the characteristics of the transition region used are the intermittency, transition length and turbulent spot properties. An analysis is made of the transition length in terms of turbulent spot variables. The nondimensional spot formation rate, required for the prediction of the transition length, is shown by the analysis to be a function of the spot spreading angle, the dimensionless spot velocity ratio and the dimensionless spot area ratio. The intermittency form of the k-epsilon equations were derived from conditionally averaged equations which have been shown to be an improvement over global-time-averaged equations for the numerical calculation of the transition region. The numerical predictions are in general good agreement with the experimental data and indicate the potential use of the method in accelerating flows. Turbulence models of the k-epsilon type are known to underpredict the transition length. The present work demonstrates how incorporating transition region characteristics improves the ability of two-equation turbulence models to simulate bypass transition for flat plates with potential application to turbine vanes and blades.

Simon, Frederick F.↗

Features of a SINDA/FLUINT model of a liquid oxygen supply line

The modeling features used in a steady-state heat transfer problem using SINDA/FLUINT are described. The problem modeled is a 125 feet long, 3 inch diameter pipe, filled with liquid oxygen flow driven by a given pressure gradient. The pipe is fully insulated in five sections. Three sections of 1 inch thick spray-on foam and two sections of vacuum jacket. The model evaluates friction, turns losses and convection heat transfer between the fluid and the pipe wall. There is conduction through the foam insulation with temperature dependent thermal conductivity. The vacuum space is modeled with radiation and gas molecular conduction, if present, in the annular gap. Heat is transferred between the outer surface and surrounding ambient by natural convection and radiation; and, by axial conduction along the pipe and through the vacuum jacket spacers and welded seal flanges. The model makes extensive use of SINDA/FLUINT basic capabilities such as the GEN option for nodes and conductors (to generate groups of nodes or conductors), the SIV option (to generate single, temperature varying conductors), the SIM option (for multiple, temperature varying conductors) and the M HX macros for fluids (to generate strings of lumps, paths, and ties representing a diabatic duct). It calls subroutine CONTRN (returns the relative location in the G-array of a network conductor, given an actual conductor number) enabling an extensive manipulation of conductor (calculates an assignment of their values) with DO loops. Models like this illustrate to the new and even to the old SINDA/FLUINT user, features of the program that are not so obvious or known, and that are extremely handy when trying to take advantage of both, the automation of the DATA headers and make surgical modifications to specific parameters of the thermal or fluid elements in the OPERATIONS portion of the model.

Simmonds, Boris G.↗

Skin friction measurement in complex flows using thin oil film techniques

The NASA Grant NAG2-261 was initiated to support a program of research to study complex flows that occur in flight and laboratory experiments by building, testing and optimizing an on-board technique for direct measurement of surface shear stress using thin oil film techniques. The program of research has proceeded under the supervision of the NASA Ames Research Center and with further cooperation from the NASA Ames-Dryden and NASA Langley Research Centers. In accordance with the original statement of work, the following research milestones were accomplished: (1) design and testing of an internally mounted one-directional skin friction meter to demonstrate the feasibility of the concept; (2) design and construction of a compact instrument capable of measuring skin friction in two directions; (3) study of transitional and fully turbulent boundary layers over a flat plate with and without longitudinal pressure gradients utilizing the compact two-directional skin friction meter; (4) study of the interaction between a turbulent boundary layer and a shock wave generated by a compression corner using the two-directional meter; and (5) flight qualification of the compact meter and accompanying electronic and pneumatic systems, preliminary installation into flight test fixture.

Source record↗

Computation of turbulent flows over backward and forward-facing steps using a near-wall Reynolds stress model

Separation and reattachment of turbulent shear layers is observed in many important engineering applications, yet it is poorly understood. This has motivated many studies on understanding and predicting the processes of separation and reattachment of turbulent shear layers. Both of the situations in which separation is induced by adverse pressure gradient, or by discontinuities of geometry, have attracted attention of turbulence model developers. Formulation of turbulence closure models to describe the essential features of separated turbulent flows accurately is still a formidable task. Computations of separated flows associated with sharp-edged bluff bodies are described. For the past two decades, the backward-facing step flow, the simplest separated flow, has been a popular test case for turbulence models. Detailed studies on the performance of many turbulence models, including two equation turbulence models and Reynolds stress models, for flows over steps can be found in the papers by Thangam & Speziale and Lasher & Taulbee). These studies indicate that almost all the existing turbulence models fail to accurately predict many important features of back step flow such as reattachment length, recovery rate of the redeveloping boundary layers downstream of the reattachment point, streamlines near the reattachment point, and the skin friction coefficient. The main objectives are to calculate flows over backward and forward-facing steps using the NRSM and to make use of the newest DNS data for detailed comparison. This will give insights for possible improvements of the turbulence model.

Ko, Sung HO↗

Direct simulation of isothermal-wall supersonic channel flow

The motivation for this work is the fact that in turbulent flows where compressibility effects are important, they are often poorly understood. A few examples of such flows are those associated with astrophysical phenomena and those found in combustion chambers, supersonic diffusers and nozzles, and over high-speed airfoils. For this project, we are primarily interested in compressibility effects near solid surfaces. Our main objective is an improved understanding of the fundamentals of compressible wall-bounded turbulence, which can in turn be used to cast light upon modeling concepts such as the Morkovin hypothesis and the Van Driest transformation. To this end, we have performed a direct numerical simulation (DNS) study of supersonic turbulent flow in a plane channel with constant-temperature walls. All of the relevant spatial and temporal scales are resolved so that no sub grid scale or turbulence model is necessary. The channel geometry was chosen so that finite Mach number effects can be isolated by comparing the present results to well established incompressible channel data. Here the fluid is assumed to be an ideal gas with constant specific heats, constant Prandtl number, and power-law temperature-dependent viscosity. Isothermal-wall boundary conditions are imposed so that a statistically stationary state may be obtained. The flow is driven by a uniform (in space) body force (rather than a mean pressure gradient) to preserve stream wise homogeneity, with the body force defined so that the total mass flux is constant.

Coleman, Gary N.↗

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↗

Computational fluid mechanics

Two papers are included in this progress report. In the first, the compressible Navier-Stokes equations have been used to compute leading edge receptivity of boundary layers over parabolic cylinders. Natural receptivity at the leading edge was simulated and Tollmien-Schlichting waves were observed to develop in response to an acoustic disturbance, applied through the farfield boundary conditions. To facilitate comparison with previous work, all computations were carried out at a free stream Mach number of 0.3. The spatial and temporal behavior of the flowfields are calculated through the use of finite volume algorithms and Runge-Kutta integration. The results are dominated by strong decay of the Tollmien-Schlichting wave due to the presence of the mean flow favorable pressure gradient. The effects of numerical dissipation, forcing frequency, and nose radius are studied. The Strouhal number is shown to have the greatest effect on the unsteady results. In the second paper, a transition model for low-speed flows, previously developed by Young et al., which incorporates first-mode (Tollmien-Schlichting) disturbance information from linear stability theory has been extended to high-speed flow by incorporating the effects of second mode disturbances. The transition model is incorporated into a Reynolds-averaged Navier-Stokes solver with a one-equation turbulence model. Results using a variable turbulent Prandtl number approach demonstrate that the current model accurately reproduces available experimental data for first and second-mode dominated transitional flows. The performance of the present model shows significant improvement over previous transition modeling attempts.

Hassan, H. A.↗