SEARCH · Search NASA
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.
A calculation was made for the flow and heat transfer in the turbulent free-convection boundary layer on a vertical flat plate. Formulas for the heat-transfer coefficient, boundary layer thickness, and the maximum velocity in the boundary layer were obtained.
With the use of Karman's integrated momentum equation for the boundary layer and data on the wall-shearing stress and heat transfer in forced-convection flow, a calculation was carried out for the flow and heat transfer in the turbulent free-convection boundary layer on a vertical flat plate. The calculation is for a fluid with a Prandtl number that is close to 1. A formula was derived for the heat-transfer coefficient that was in good agreement with experimental data in the range of Grashof numbers from 10sup10 to 10sup12. Because of the good agreement between the theoretical formula and the experimental data, the formula may be used to obtain data for high Grashof numbers. The calculation also yielded formulas for the maximum velocity in the boundary layer and for boundary-layer thickness.
The present report consists of two parts. The first part deals with the two-dimensional stationary flow in the presence of local supersonic zones. A numerical method of integration of the equation of gas dynamics is developed. Proceeding from solutions at great distance from the body the flow pattern is calculated step by step. Accordingly the related body form is obtained at the end of the calculation. The second part treats the relationship between the displacement thickness of laminar and turbulent boundary layers and the pressure distribution at high speeds. The stability of the boundary layer is investigated, resulting in basic differences in the behavior of subsonic and supersonic flows. Lastly, the decisive importance of the boundary layer for the pressure distribution, particularly for thin profiles, is demonstrated.
Explore the source record for details and available documents.
Two-dimensional eigenvalue analysis is used on a massive scale to study spatial instabilities of compressible shear flows with two inhomogeneous directions. The main focus of the study is crossflow dominated swept-wing boundary layers although the methodology can also be applied to study other type of flows, such as the attachment-line flow. Certain unique aspects of formulating a spatial, two-dimensional eigenvalue problem for the secondary instability of finite amplitude crossflow vortices are discussed, namely, fixing the spatial growth direction unambiguously through a non-orthogonal formulation of the linearized disturbance equations. A primary test case used for parameter study corresponds to the low-speed, NLF-0415(b) airfoil configuration as tested in the ASU Unsteady Wind Tunnel, wherein a spanwise periodic array of roughness elements was placed near the leading edge in order to excite stationary crossflow modes with a specified fundamental wavelength. The two classes of flow conditions selected for this analysis include those for which the roughness array spacing corresponds to either the naturally dominant crossflow wavelength, or a subcritical wavelength that serves to reduce the growth of the naturally excited dominant crossflow modes. Numerical predictions are compared with the measured database, both as indirect validation for the spatial instability analysis and to provide a basis for comparison with a higher Reynolds number, supersonic swept-wing configuration. Application of the eigenvalue analysis to the supersonic configuration reveals that a broad spectrum of stationary crossflow modes can sustain sufficiently strong secondary instabilities as to potentially cause transition over this configuration. Implications of this finding for transition control in swept wing boundary layers are examined.
An asymptotic description is proposed for supersonic laminar flow over a wedge or a backward-facing step, for large Reynolds number and for a base or step height which is small compared with the boundary-layer length. The analysis is carried out for adiabatic wall conditions and a viscosity coefficient proportional to temperature. In a particular limit corresponding to a very thick boundary layer, a similarity law is obtained for the base pressure. For a thinner boundary layer an asymptotic form for the base pressure is obtained which shows the dependence on the parameters explicitly and which permits good agreement with experiment. This latter result is based on an inviscid-flow approximation for the corner expansion and for reattachment with viscous forces important primarily in a thin sublayer about the dividing streamline. A prediction of the pressure distribution at reattachment is given and the result is compared with experimental pressure distributions.
A set of higher-order boundary-layer equations is derived valid for three-dimensional compressible flows. The equations are written in a generalized curvilinear coordinate system, in which the surface coordinates are nonorthogonal; the third axis is restricted to be normal to the surface. Also, higher-order viscous terms which are retained depend on the surface curvature of the body. Thus, the equations are suitable for the calculation of the boundary layer about arbitrary vehicles. As a starting point, the Navier-Stokes equations are derived in a tensorian notation. Then by means of an order-of-magnitude analysis, the boundary-layer equations are developed. To provide an interface between the analytical partial differentiation notation and the compact tensor notation, a brief review of the most essential theorems of the tensor analysis related to the equations of the fluid dynamics is given. Many useful quantities, such as the contravariant and the covariant metrics and the physical velocity components, are written in both notations.
Langley Stability and Transition Analysis Code (LASTRAC) is a general-purpose, physics-based transition prediction code released by NASA for laminar flow control studies and transition research. This paper describes the LASTRAC extension to general three-dimensional (3D) boundary layers such as finite swept wings, cones, or bodies at an angle of attack. The stability problem is formulated by using a body-fitted nonorthogonal curvilinear coordinate system constructed on the body surface. The nonorthogonal coordinate system offers a variety of marching paths and spanwise waveforms. In the extreme case of an infinite swept wing boundary layer, marching with a nonorthogonal coordinate produces identical solutions to those obtained with an orthogonal coordinate system using the earlier release of LASTRAC. Several methods to formulate the 3D parabolized stability equations (PSE) are discussed. A surface-marching procedure akin to that for 3D boundary layer equations may be used to solve the 3D parabolized disturbance equations. On the other hand, the local line-marching PSE method, formulated as an easy extension from its 2D counterpart and capable of handling the spanwise mean flow and disturbance variation, offers an alternative. A linear stability theory or parabolized stability equations based N-factor analysis carried out along the streamline direction with a fixed wavelength and downstream-varying spanwise direction constitutes an efficient engineering approach to study instability wave evolution in a 3D boundary layer. The surface-marching PSE method enables a consistent treatment of the disturbance evolution along both streamwise and spanwise directions but requires more stringent initial conditions. Both PSE methods and the traditional LST approach are implemented in the LASTRAC.3d code. Several test cases for tapered or finite swept wings and cones at an angle of attack are discussed.
It was previously observed that an incident acoustic field on a wing with laminar flow can cause transition to turbulent flow if the fluctuating acoustic velocities are of sufficient amplitude and in the critical frequency range for an unstable laminar boundary layer. A section of a wing was modified with a natural laminar flow (NLF) glove to allow direct measurement of the effect of varying engine noise on the extent of laminar flow. The flight test program was completed in June, 1985. At each flight condition, the engine power was varied from about 2600 r/min (idle) to about 4500 r/min (maximum continuous power). The spectral data provides considerable insight into the influences of the various sound sources that contribute to the overall noise levels. Additional analysis will be required to assess the impact of these sources on boundary layer transition. These results demonstrate that substantial laminar flow on the wing of a transport configuration with wing-mounted engines can be obtained.
Shock generated plasma boundary layers, discussing equilibrium ionization and quantitative schlieren technique
Papers are presented on the application of stability theory to laminar flow control, secondary instabilities in boundary layers, a Floquet analysis of secondary instability in shear flows, and the generation of Tollmien-Schlichting waves by long wavelength free stream disturbances. Also considered are numerical experiments on boundary-layer receptivity, short-scale inviscid instabilities in the flow past surface-mounted obstacles, wave phenomena in a high Reynolds number compressible boundary layer, and instability of time-periodic flows. Other topics include high frequency Rayleigh instability of Stokes layers, stability and resonance in grooved-channel flows, finite length Taylor Couette flow, and vortical structures in the breakdown stage of transition.
Boundary layer stability was analyzed for the HIFiRE-1 flight vehicle geometry for ground tests conducted at the CUBRC LENS I hypersonic shock test facility and the Langley Research Center (LaRC) 20- inch Mach 6 Tunnel. Boundary layer stability results were compared to transition onset location obtained from discrete heat transfer measurements from thin film gauges during the CUBRC test and spatially continuous heat transfer measurements from thermal phosphor paint data during the LaRC test. The focus of this analysis was on conditions at non-zero angles of attack as stability analysis has already been performed at zero degrees angle of attack. Also, the transition onset data obtained during flight testing was at nonzero angles of attack, so this analysis could be expanded in the future to include the results of the flight test data. Stability analysis was performed using the 2D parabolized stability software suite STABL (Stability and Transition Analysis for Hypersonic Boundary Layers) developed at the University of Minnesota and the mean flow solutions were computed using the DPLR finite volume Navier-Stokes computational fluid dynamics (CFD) solver. A center line slice of the 3D mean flow solution was used for the stability analysis to incorporate the angle of attack effects while still taking advantage of the 2D STABL software suite. The N-factors at transition onset and the value of Re(sub theta)/M(sub e), commonly used to predict boundary layer transition onset, were compared for all conditions analyzed. Ground test data was analyzed at Mach 7.2 and Mach 6.0 and angles of attack of 1deg, 3deg and 5deg. At these conditions, the flow was found to be second mode dominant for the HIFiRE-1 slender cone geometry. On the leeward side of the vehicle, a strong trend of transition onset location with angle of attack was observed as the boundary layer on the leeward side of the vehicle developed inflection points at streamwise positions on the vehicle that correlated to angle of attack. Inflection points are a strong instability mechanism that lead to rapid breakdown and transition to turbulence. The transition onset location on the windward side of the vehicle displayed no trend with angle of attack or freestream Reynolds number and transition was observed farther down the vehicle than observed on the leeward side of the vehicle. In analysis of both windward and leeward sides of the vehicle, use of the N factor methodology to develop trends to predict boundary layer transition onset showed improvements over the Re(sub theta)/M(sub e) empirical correlation methodology. Stronger correlations and less scatter in the data were observed when using the N factor method for these cases.
Flight test results are reported from an experiment designed to study the detailed growth of disturbances in the laminar boundary layer. A gloved wing section incorporating closely-spaced flush-mounted streamwise-located instrumentation for measuring instability frequencies and amplitude growths as well as pressure distributions was used. The growth of Tollmien-Schlichting (T-S) and crossflow instabilities is predicted by the linear e exp n method and compared to the measured boundary-layer disturbance frequencies. The predictions showed good agreement with the measured data. The results exhibited fair agreement with previous n(T-S) and n(CF) flight correlations for several of the conditions analyzed. It is inferred from the high n(T-S) values for these data that moderately swept wings at compressible speeds can withstand higher combinations of n(T-S) and n(CF) values and still remain laminar than previously thought.
The present critical-layer asymptotic analysis for the nonlinear interaction of frequency-detuned modes in boundary-layer transition indicates that the interaction between a plane mode at the fundamental frequency and a pair of symmetrical oblique waves at the near-subharmonic frequency amplifies another pair of symmetrical oblique waves at the 'mirror frequency'. This type of interaction is stronger in the frequency-detuned case than the resonant triad case, and leads to a sharp drop in the oblique waves' peak with small detuning.
There is a growing body of observational and theoretical evidence which suggests that local climate characteristics are associated with variations in the earth's surface. The link between surface variability and local-scale processes must be made if we are to improve our understanding of the feedback mechanisms involved in surface-atmosphere dynamics. However, to understand these interactions, the surface-atmosphere interface must be studied as a large-scale spatial system. Lidars are ideal tools to study the spatial properties of the atmosphere. The described techniques were developed for use with the Los Alamos Water Raman-Lidar, but are applicable to many other types of lidar. The methodology of the analysis of lidar data is summarized in order to determine meteorological parameters in the atmospheric boundary layer. The techniques are not exhaustive but are intended to show the depth and breadth of the information which can be obtained from lidars. Two methods for the computation of water-vapor fluxes were developed. The first uses the fact that the water vapor concentration in the vertical direction follows a logarithmic profile when corrected for atmospheric stability. The second method involves using inertial dissipation techniques in which lidar-derived spatial and temporal power spectra are used to determine the flux.
The viscous subsonic flow past two-dimensional and infinite-span swept multi-component airfoils is studied theoretically and experimentally. The computerized analysis is based on iteratively coupled boundary-layer and potential-flow analysis. The method, which is restricted to flows with only slight separation, gives surface pressure distribution, chordwise and spanwise boundary-layer characteristics, lift, drag, and pitching moment for airfoil configurations with up to four elements. Merging confluent boundary layers are treated. Theoretical predictions are compared with an exact theoretical potential flow solution and with experimental measures made in the Ames 40- by 80-Foot Wind Tunnel for both two-dimensional and infinite-span swept wing configurations. Section lift characteristics are accurately predicted for zero and moderate sweep angles where flow separation effects are negligible.
A dual potential decomposition of the velocity field into a scalar and a vector potential function is extended to three dimensions and used in the finite-difference simulation of steady three-dimensional inviscid rotational flows and viscous flow. The finite-difference procedure was used to simulate the flow through the 80 by 120 ft wind tunnel at NASA Ames Research Center. Rotational flow produced by the stagnation pressure drop across vanes and screens which are located at the entrance of the inlet is modeled using actuator disk theory. Results are presented for two different inlet vane and screen configurations. The numerical predictions are in good agreement with experimental data. The dual potential procedure was also applied to calculate the viscous flow along two and three dimensional troughs. Viscous effects are simulated by injecting vorticity which is computed from a boundary layer algorithm. For attached flow over a three dimensional trough, the present calculations are in good agreement with other numerical predictions. For separated flow, it is shown from a two dimensional analysis that the boundary layer approximation provides an accurate measure of the vorticity in regions close to the wall; whereas further away from the wall, caution has to be exercised in using the boundary-layer equations to supply vorticity to the dual potential formulation.
The purpose of these experiments and analysis was to augment the heating database and tools used for assessment of impact-induced shallow-cavity damage to the thermal protection system of the Space Shuttle Orbiter. The effect of length and depth on the local heating disturbance of rectangular cavities tested at hypersonic freestream conditions has been globally assessed using the two-color phosphor thermography method. These rapid-response experiments were conducted in the Langley 31-Inch Mach 10 Tunnel and were initiated immediately prior to the launch of STS-114, the initial flight in the Space Shuttle Return-To-Flight Program, and continued during the first week of the mission. Previously-designed and numerically-characterized blunted-nose baseline flat plates were used as the test surfaces. Three-dimensional computational predictions of the entire model geometry were used as a check on the design process and the two-dimensional flow assumptions used for the data analysis. The experimental boundary layer state conditions were inferred using the measured heating distributions on a no-cavity test article. Two test plates were developed, each containing 4 equally-spaced spanwise-distributed cavities. The first test plate contained cavities with a constant length-to-depth ratio of 8 with design point depth-to-boundary-layer-thickness ratios of 0.1, 0.2, 0.35, and 0.5. The second test plate contained cavities with a constant design point depth-to-boundary-layer-thickness ratio of 0.35 with length-to-depth ratios of 8, 12, 16, and 20. Cavity design parameters and the test condition matrix were established using the computational predictions. Preliminary results indicate that the floor-averaged Bump Factor (local heating rate nondimensionalized by upstream reference) at the tested conditions is approximately 0.3 with a standard deviation of 0.04 for laminar-in/laminar-out conditions when the cavity length-to-boundary-layer thickness is between 2.5 and 10 and for cavities in the depth-to-boundary-layer-thickness range of 0.3 to 0.8. Over this same range of conditions and parameters, preliminary results also indicate that the maximum Bump Factor on the cavity centerline falls between 2.0 and 2.75, as long as the cavity-exit conditions remain laminar. Cavities with length-to-boundary-layer-thickness ratio less than 2.5 can not be easily classified with this approach and require further analysis.