SEARCH · Engineering Papers
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.
Results of experimental study of regularities of a natural transition of a laminar boundary layer to a turbulent layer at low subsonic air flow velocities are presented, analyzed and compared with theory and model experiments.
The viscous shock layer equations applicable to hypersonic laminar, transitional, and turbulent flows of a perfect gas over two-dimensional plane or axially symmetric blunt bodies are presented. The equations are solved by means of an implicit finite difference scheme, and the results are compared with a turbulent boundary layer analysis. The agreement between the two solution procedures is satisfactory for the region of flow where streamline swallowing effects are negligible. For the downstream regions, where streamline swallowing effects are present, the expected differences in the two solution procedures are evident.
The results of a finite-difference calculation of hypersonic transitional boundary-layer flow are presented. Precursor-transition and low Reynolds number turbulence effects are included. These calculation results are compared with data from a source where profiles are available at the beginning, middle and end of the hypersonic transitional flow. The results indicate that the low Reynolds number effect is the most important.
This work presents the experimental results of research on the influence of gravity on flow pattern transitions, pressure drop and flow characteristics for cocurrent gas-liquid two-phase flow through packed columns. The flow pattern transition data indicates that the pulse flow regime exists over a wider range of gas and liquid flow rates under reduced gravity conditions compared to normal gravity cocurrent down-flow. This is illustrated by comparing the flow regime transitions found in reduced gravity with the transitions predicted by Talmor. Next, the effect of gravity on the total pressure drop in a packed column is shown to depend on the flow regime. The difference is roughly equivalent to the liquid static head for bubbly flow but begins to decrease at the onset of pulse flow. As the spray flow regime is approached by increasing the gas to liquid ratio, the effect of gravity on pressure drop becomes negligible. Finally, gravity tends to suppress the amplitude of each pressure pulse. An example of this phenomenon is presented.
The propagation of power flow through a dynamically loaded beam model with 90 degree bends is investigated using NASTRAN and McPOW. The transitioning of power flow types (axial, torsional, and flexural) is observed throughout the structure. To get accurate calculations of the torsional response of beams using NASTRAN, torsional inertia effects had to be added to the mass matrix calculation section of the program. Also, mass effects were included in the calculation of BAR forces to improve the continuity of power flow between elements. The importance of including all types of power flow in an analysis, rather than only flexural power, is indicated by the example. Trying to interpret power flow results that only consider flexural components in even a moderately complex problem will result in incorrect conclusions concerning the total power flow field.
The work in this report was conducted at NASA Ames Research Center during the period from August 1993 to January 1995 deals with the direct numerical simulation of transitional and turbulent flow at low Mach numbers using high-order-accurate finite-difference techniques. A computation of transition to turbulence of the spatially-evolving boundary layer on a heated flat plate in the presence of relatively high freestream turbulence was performed. The geometry and flow conditions were chosen to match earlier experiments. The development of the momentum and thermal boundary layers was documented. Velocity and temperature profiles, as well as distributions of skin friction, surface heat transfer rate, Reynolds shear stress, and turbulent heat flux were shown to compare well with experiment. The numerical method used here can be applied to complex geometries in a straightforward manner.
Recent experiments have shown that low-density gas jets injected into a high-density ambient gas undergo an instability mode, leading to highly-periodic oscillations in the flow-field for certain conditions. The transition from laminar to turbulent flow in these jets is abrupt, without the gradual change in scales. Even the fine scale turbulent structure repeats itself with extreme regularity from cycle to cycle. Similar observations were obtained in buoyancy-dominated and momentum-dominated jets characterized by the Richardson numbers, Ri = [gD(rho(sub a)-rho(sub j))/rho(sub j)U(sub j)(exp 2) ] where g is the gravitational acceleration, D is the jet diameter, rho(sub a) and rho(sub a) are, respectively, the free-stream and jet densities, and U(sub j) is the mean jet exit velocity. At high Richardson numbers, the instability is presumably caused by buoyancy since the flow-oscillation frequency (f) or the Strouhal number, St = [fD/U(sub j)] scales with Ri. In momentum-dominated jets, however, the Strouhal number of the oscillating flow is relatively independent of the Ri. In this case, a local absolute instability is predicted in the potential core of low-density jets with S [= rho(sub j)/rho(sub a)] < 0.7, which agrees qualitatively with experiments. Although the instability in gas jets of high Richardson numbers is attributed to buoyancy, direct physical evidence has not been acquired in experiments. If the instability is indeed caused by buoyancy, the near-field flow structure of the jet will change significantly when the buoyancy is removed, for example, in the microgravity environment. Thus, quantitative data on the spatial and temporal evolutions of the instability, length and time scale of the oscillating mode and its effects on the mean flow and breakdown of the potential core are needed in normal and microgravity to delineate gravitational effects in buoyant jets. In momentum dominated low-density jets, the instability is speculated to originate in the potential core. However, experiments have not succeeded in identifying the direct physical cause of the instability. For example, the theory predicts an oscillating mode for S<0.62 in the limit of zero momentum thickness, which contradicts with the experimental findings of Kyle and Sreenivasan. The analyses of momentum-dominated jets neglect buoyancy effects because of the small Richardson number. Although this assumption is appropriate in the potential core, the gravitational effects are important in the annular region surrounding the jet, where the density and velocity gradients are large. This reasoning provides basis for the hypothesis that the instability in low Richardosn number jets studied by Kyle and Sreenivasan and Monkewitz et al. is caused by buoyancy. The striking similarity in characteristics of the instability and virtually the identical conclusions reached by Subbarao and Cantwell in buoyant (Ri>0.5) helium jets on one hand and by Kyle and Sreenivasan in momentum-dominated (Ri<1x10(exp -3)) helium jets on the other support this hypothesis. However, quantitative experiments in normal and microgravity are necessary to obtain direct physical evidence of buoyancy effects on the flow instability and structure of momentum-dominated low-density jets. The primary objective of this new research project is to quantify how buoyancy affects the flow instability and structure in the near field of low-density jets. The flow will be described by the spatial and temporal evolutions of the instability, length and time scales of the oscillating mode, and the mean and fluctuating concentration fields. To meet this objective, concentration measurements will be obtained across the whole field using quantitative Rainbow Schlieren Deflectometry, providing spatial resolution of 0.1mm and temporal resolution of 0.017s to 1ms. The experimental effort will be supplemented with linear stability analysis of low-density jets by considering buoyancy. The first objective of this research is to investigate the effects of gravity on the flow instability and structure of low-density jets. The flow instability in these jets has been attributed to buoyancy. By removing buoyancy in our experiments, we seek to obtain the direct physical evidence of the instability mechanism. In the absence of the instability, the flow structure will undergo a significant change. We seek to quantify these changes by mapping the flow field (in terms of the concentration profiles) of these jets at non-buoyant conditions. Such information is presently lacking in the existing literature. The second objective of this research is to determine if the instability in momentum-driven, low-density jets is caused by buoyancy. At these conditions, the buoyancy effects are commonly ignored because of the small Richardson based on global parameters. By eliminating buoyancy in our experiments, globally as well as locally, we seek to examine the possibility that the instability mechanism in self-excited, buoyant or momentum-driven jets is the same. To meet this objective, we would quantify the jet flow in normal and microgravity, while systematically decreasing the Richardson number from buoyancy-driven to momentum driven flow regime. The third objective of this research is to perform a linear stability analysis of low-density gas jets by including the gravitational effects. The flow oscillations in these jets are attributed to an absolute instability, whereby the disturbance grows exponentially at the site to ultimately contaminate the entire flow field. We seek to study the characteristics of both convective and absolute instabilities and demarcate the boundary between them.
The transition process characteristics of flows over swept wings were computationally modelled. The crossflow instability and crossflow/T-S wave interaction are analyzed through the numerical solution of the full three dimensional Navier-Stokes equations including unsteadiness, curvature, and sweep. The leading-edge region of a swept wing is considered in a three-dimensional spatial simulation with random disturbances as the initial conditions.
Experimental results are presented to document hydrodynamic and thermal development of flat-plate boundary layers undergoing natural transition. Local heat transfer coefficients, skin friction coefficients and profiles of velocity, temperature and Reynolds normal and shear stresses are presented. A case with no transition and transitional cases with 0.68% and 2.0% free-stream disturbance intensities were investigated. The locations of transition are consistent with earlier data. A late-laminar state with significant levels of turbulence is documented. In late-transitional and early-turbulent flows, turbulent Prandtl number and conduction layer thickness values exceed, and the Reynolds analogy factor is less than, values previously measured in fully turbulent flows.
Experimental results are presented to document hydrodynamic and thermal development of flat-plate boundary layers undergoing natural transition. Local heat transfer coefficients, skin friction coefficients and profiles of velocity, temperature and Reynolds normal and shear stresses are presented. A case with no transition and transitional cases with 0.68 percent and 2.0 percent free-stream disturbance intensities were investigated. The locations of transition are consistent with earlier data. A late-laminar state with significant levels of turbulence is documented. In late-transitional and early-turbulent flows, turbulent Prandtl number and conduction layer thickness values exceed, and the Reynolds analogy factor is less than, values previously measured in fully turbulent flows.
A computer program is described which solves the compressible laminar, transitional, or turbulent boundary-layer equations for planar or axisymmetric flows. Three-point implicit difference relations are used to reduce the momentum and energy equations to finite-difference form. These equations are solved simultaneously without iteration. Turbulent flow is treated by the inclusion of either a two-layer eddy-viscosity model or a mixing-length formulation. The eddy conductivity is related to the eddy viscosity through a static turbulent Prandtl number which may be an arbitrary function of the distance from the wall boundary. The transitional boundary layer is treated by the inclusion of an intermittency function which modifies the fully turbulent model. The laminar-boundary-layer equations are recovered when the intermittency is zero, and the fully turbulent equations are solved when the intermittency is unity.
An algebraic eddy-viscosity model is derived from the renormalization group (RNG) theory of turbulence. A new length scale, based on boundary layer characteristics (displacement thickness, shape factor), is proposed. The model was applied to transitional boundary layer flow over a flat plate. Integral characteristics, such as the total skin friction coefficient, and mean velocity profile across the boundary layer, are found to be in good agreement with experimental data.
Various secondary and tertiary instabilities in plane channael flow are explored via time-dependent numerical simulations using the incompressible Navier-Stokes equations. Comparisons are made between transitional flows at Reynolds numbers 1500, 5000, and 8000. The lambda vortex, detached shear layer, and inverted vortex regions are identified, and the origin of the latter is explained. The laminar breakdown of the Re = 1500 flow is computed with high resolution, and the nature of its ensuing hairpin eddies is clarified by numerical particle paths. The potential of center-mode, rather than wall-mode transitions is proposed, and the resulting flow structure is described.
A computer program was developed to do the calculations for two-dimensional or axisymmetric configurations from low speeds to hypersonic speeds with arbitrary streamwise pressure, temperature, and Mach number distributions. Options are provided for obtaining initial conditions either from experimental information or from a theoretical similarity solution. The transition region can be described either by an arbitrary distribution of intermittency or by a function based on Emmons' probability theory. Correlations were developed for use in estimating the parameters of the theoretical intermittency function. Correlations obtained from other sources are used for estimating the transition point. Comparisons were made between calculated and measured boundary layer quantities for laminar, transitional, and turbulent flows on flat plates, cones, cone flares, and a waisted body of revolution. Excellent agreement was obtained between the present theory and two other theories based on the method of finite differences. The intermittency required to reproduce some experimental heat transfer results in hypersonic flow was found to be quite different from the theoretical function. It is suggested that the simple probability theory of Emmons may not be valid for representing the intermittency of hypersonic transitional boundary layers and that the program could be useful as a tool for detailed study of the intermittency of the transition region.
Results of an experimental study of flow separation and transition in either attached boundary layers or separated shear layers over the suction surface of a simulation of a low-pressure turbine airfoil flow are presented. Detailed velocity profiles were measured with the hot-wire technique. Static pressure distributions are also presented. Flow transition is documented using measured intermittency distributions in the boundary layer and the separated shear layer. Cases for Reynolds numbers of 50,000, 100,000, 200,000 and 300,000 are reported. These Reynolds numbers are based on suction surface length and exit velocity. Three Free Stream Turbulence Intensity values, 0.5%, 2.5% and 10%, are represented. Flow separation is observed for all the low-FSTI cases. Of these, the lowest Reynolds number case was not able to complete transition of the shear layer and the separation bubble persisted over the entire blade surface. For the other low-FSTI cases, transition is observed in the shear layer over the separation bubble. This transition proceeded quickly, spreading rapidly toward the wall. Elevated FSTI drives an earlier transition than in the low-FSTI cases and the separation bubbles are smaller. For the highest Reynolds number cases with 2.5% and 10% FSTI, transition is of the attached boundary layer and no separation exists. Flow separation with shear flow transition is observed for the lower-Re cases. Models for intermittency and transition length and location from the modern literature are assessed.
An investigation was conducted on the limiting case of a gas flow when the constant pressure in the surrounding medium is exactly equal to the critical pressure for the given initial state of the gas.
This study investigated the Navier-Stokes computations of the surface heat transfer coefficients of a transition duct flow. A transition duct from an axisymmetric cross section to a non-axisymmetric cross section, is usually used to connect the turbine exit to the nozzle. As the gas turbine inlet temperature increases, the transition duct is subjected to the high temperature at the gas turbine exit. The transition duct flow has combined development of hydraulic and thermal entry length. The design of the transition duct required accurate surface heat transfer coefficients. The Navier-Stokes computational method could be used to predict the surface heat transfer coefficients of a transition duct flow. The Proteus three-dimensional Navier-Stokes numerical computational code was used in this study. The code was first studied for the computations of the turbulent developing flow properties within a circular duct and a square duct. The code was then used to compute the turbulent flow properties of a transition duct flow. The computational results of the surface pressure, the skin friction factor, and the surface heat transfer coefficient were described and compared with their values obtained from theoretical analyses or experiments. The comparison showed that the Navier-Stokes computation could predict approximately the surface heat transfer coefficients of a transition duct flow.
We present experimental data on flow pattern transitions, pressure drop and flow characteristics for cocurrent gas-liquid flow through packed columns in microgravity. The flow pattern transition data indicates that the pulse flow regime exists over a wider range of gas and liquid flow rates under microgravity conditions compared to 1-g and the widely used Talmor map in 1-g is not applicable for predicting the transition boundaries. A new transition criterion between bubble and pulse flow in microgravity is proposed and tested using the data. Since there is no static head in microgravity, the pressure drop measured is the true frictional pressure drop. The pressure drop data, which has much smaller scatter than most reported 1-g data clearly shows that capillary effects can enhance the pressure drop (especially in the bubble flow regime) as much as 200% compared to that predicted by the single phase Ergun equation. The pressure drop data are correlated in terms of a two-phase friction factor and its dependence on the gas and liquid Reynolds numbers and the Suratman number. The influence of gravity on the pulse amplitude and frequency is also discussed and compared to that under normal gravity conditions. Experimental work is planned to determine the gas-liquid mass transfer coefficients. Because of enhanced interfacial effects, we expect the gas-liquid transfer coefficients k(L)a and k(G)a (where a is the gas-liquid interfacial area) to be higher in microgravity than in normal gravity at the same flow conditions. This will be verified by gas absorption experiments, with and without reaction in the liquid phase, using oxygen, carbon dioxide, water and dilute aqueous amine solutions. The liquid-solid mass transfer coefficient will also be determined in the bubble as well as the pulse flow regimes using solid benzoic acid particles in the packing and measuring their rate of dissolution. The mass transfer coefficients in microgravity will be compared to those in normal gravity cocurrent flow to determine the mass transfer enhancement and propose new mass transfer correlations for two-phase gas-liquid flows through packed beds in microgravity.