Langley Full-scale-tunnel Investigation of the Fuselage Boundary Layer on a Typical Fighter Airplane with a Single Liquid-cooled Engine
Explore the source record for details and available documents.
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.
Explore the source record for details and available documents.
For each of two inlet-boundary-layer thicknesses, performance and boundary-layer characteristics have been determined for a 12 degree, 10-inch-inlet-diameter diffuser, a 12 degree, 21-inch-inlet-diameter diffuser, and a 23 degree, 21-inch-inlet-diameter diffuser. The investigation covered an inlet Mach number range from about 0.10 to coking. The corresponding inlet Reynolds number, based on inlet diameter, varied from about 0.5 x 10(6) to 7.5 x 10(6).
The Navier Stokes equation is solved computationally by the Dispersion-Relation-Preserving (DRP) scheme for the flow and acoustic fields associated with a laminar boundary layer flow over an automobile door cavity. In this work, the flow Reynolds number is restricted to R(sub delta*) < 3400; the range of Reynolds number for which laminar flow may be maintained. This investigation focuses on two aspects of the problem, namely, the effect of boundary layer thickness on the cavity tone frequency and intensity and the effect of the size of the computation domain on the accuracy of the numerical simulation. It is found that the tone frequency decreases with an increase in boundary layer thickness. When the boundary layer is thicker than a certain critical value, depending on the flow speed, no tone is emitted by the cavity. Computationally, solutions of aeroacoustics problems are known to be sensitive to the size of the computation domain. Numerical experiments indicate that the use of a small domain could result in normal mode type acoustic oscillations in the entire computation domain leading to an increase in tone frequency and intensity. When the computation domain is expanded so that the boundaries are at least one wavelength away from the noise source, the computed tone frequency and intensity are found to be computation domain size independent.
An extensive parametric study of vortices shed from airfoil vortex generators has been conducted to determine the dependence of initial vortex circulation and peak vorticity on elements of the airfoil geometry and impinging flow conditions. These elements include the airfoil angle of attack, chord length, span, aspect ratio, local boundary layer thickness, and free stream Mach number. In addition, the influence of airfoil-to-airfoil spacing on the circulation and peak vorticity has been examined for pairs of co-rotating and counter-rotating vortices. The vortex generators were symmetric airfoils having a NACA-0012 cross-sectional profile. These airfoils were mounted either in isolation, or in pairs, on the surface of a straight pipe. The turbulent boundary layer thickness to pipe radius ratio was about 17 percent. The circulation and peak vorticity data were derived from cross-plane velocity measurements acquired with a seven-hole probe at one chord-length downstream of the airfoil trailing edge location. The circulation is observed to be proportional to the free-stream Mach number, the angle-of-attack, and the span-to-boundary layer thickness ratio. With these parameters held constant, the circulation is observed to fall off in monotonic fashion with increasing airfoil aspect ratio. The peak vorticity is also observed to be proportional to the free-stream Mach number, the airfoil angle-of-attack, and the span-to-boundary layer thickness ratio. Unlike circulation, however, the peak vorticity is observed to increase with increasing aspect ratio, reaching a peak value at an aspect ratio of about 2.0 before falling off again at higher values of aspect ratio. Co-rotating vortices shed from closely spaced pairs of airfoils have values of circulation and peak vorticity under those values found for vortices shed from isolated airfoils of the same geometry. Conversely, counter-rotating vortices show enhanced values of circulation and peak vorticity when compared to values obtained in isolation. The circulation may be accurately modeled with an expression based on Prandtl's relationship between finite airfoil circulation and airfoil geometry. A correlation for the peak vorticity has been derived from a conservation relationship equating the moment at the airfoil tip to the rate of angular momentum production of the shed vortex, modeled as a Lamb (ideal viscous) vortex. This technique provides excellent qualitative agreement to the observed behavior of peak vorticity for low aspect ratio airfoils typically used as vortex generators.
To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, three-dimensional direct numerical simulations (DNS) were conducted to investigate the interaction of a tunnel-like acoustic disturbance field with a Mach 8 laminar boundary layer over a 7 deg half-angle blunt cone with a nose radius of 𝑅 𝑛 = 5.2 mm. The flow conditions and the axisymmetric cone geometry matched those measured in the Sandia Hypersonic Wind Tunnel at Mach 8 (Sandia HWT-8), and the incident freestream acoustic disturbances in the DNS were generated by acoustic radiation from the nozzle-wall turbulent boundary layer in the same tunnel. Consistent with the predictions of the nonmodal instability analysis, the study shows the emergence of a second spectral peak for moderately oblique waves between 150 kHz and 200 kHz in the aft part of the cone. The spectral energy of temperature and streamwise velocity fluctuations is initially concentrated within the entropy layer. These fluctiations are eventually “swallowed” by the boundary layer as the entropy-layer thickness 𝛿 𝑆 becomes similar to the boundary-layer thickness 𝛿h𝑡 at downstream locations. The numerical schlieren contours show the same inclined structures commonly observed in blunt cone experiments, and these structures are found to be three dimensional and advect with a constant speed approximately equal to the mean flow speed near the edge of the boundary layer.
To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, three-dimensional direct numerical simulations (DNS) were conducted to investigate the interaction of a tunnel-like acoustic disturbance field with a Mach 8 laminar boundary layer over a 7 deg half-angle blunt cone with a nose radius of 𝑅 𝑛 = 5.2 mm. The flow conditions and the axisymmetric cone geometry matched those measured in the Sandia Hypersonic Wind Tunnel at Mach 8 (Sandia HWT-8), and the incident freestream acoustic disturbances in the DNS were generated by acoustic radiation from the nozzle-wall turbulent boundary layer in the same tunnel. Consistent with the predictions of the nonmodal instability analysis, the study shows the emergence of a second spectral peak for moderately oblique waves between 150 kHz and 200 kHz in the aft part of the cone. The spectral energy of temperature and streamwise velocity fluctuations is initially concentrated within the entropy layer. These fluctiations are eventually “swallowed” by the boundary layer as the entropy-layer thickness 𝛿 𝑆 becomes similar to the boundary-layer thickness 𝛿h𝑡 at downstream locations. The numerical schlieren contours show the same inclined structures commonly observed in blunt cone experiments, and these structures are found to be three dimensional and advect with a constant speed approximately equal to the mean flow speed near the edge of the boundary layer.
A rational foundation is presented for the application of the linear shear flows to transition prediction, and an explicit method is given for carrying out the necessary calculations. The expansions used are shown to be complete. Sample calculations show that a typical boundary layer is very sensitive to vorticity disturbances in the inner boundary layer, near the critical layer. Vorticity disturbances three or four boundary layer thicknesses above the boundary are nearly uncoupled from the boundary layer in that the amplitudes of the discrete Tollmien-Schlicting waves are an extremely small fraction of the amplitude of the disturbance.
A rational foundation is provided for the application of the linear stability theory of parallel shear flows to transition prediction. An explicit method is given for carrying out the necessary calculations. The expansions are shown to be complete. Sample calculations show that a typical boundary layer is very sensitive to vorticity disturbance in the inner boundary layer near the critical layer. Vorticity disturbances three or four boundary layer thicknesses above the boundary are nearly uncoupled from the boundary layer, in that the amplitudes of the discrete Tollmein-Schlichting waves are an extremely small fraction of the amplitude of the disturbance.
An experimental study is performed to determine the mechanisms by which a forward-facing step impacts the growth and breakdown to turbulence of the stationary cross-flow instability. Particle image velocimetry measurements are obtained in the boundary layer of a 30◦swept flat plate with a pressure body. Step heights range from 53 % to 71 % of the boundary-layer thickness. The critical step height is approximately 60 % of the boundary-layer thickness for the current study, although it is also shown that the critical step height depends on the initial amplitude of the stationary cross-flow vortices. For the critical cases, the stationary cross-flow amplitude grows sharply downstream of the step, decays for a short region and then grows again. The initial growth region is linear, and can be explained primarily through the impact of the step on the mean flow. Namely, the step causes abrupt changes to the mean flow, resulting in large values of wall-normal shear, as well as highly inflectional profiles, due to either cross-flow reversal, separation or both. These inflectional profiles are highly unstable for the stationary cross-flow. Additionally,the reversed flow regions are significantly modulated by the stationary cross-flow vortices. The second region of growth occurs due to the stationary-cross-flow-induced modulation of the shear layer, which leads to multiple smaller wavelength stream wise vortices. High-frequency fluctuations indicate that the unsteady transition mechanism for the critical cases relates to the shedding of vortices downstream of reattachment of the modulated separated regions.
An experimental investigation has been made of turbulent boundary-layer separation associated with compression corners, curved surfaces of various radii, and incident shock waves. The purpose of the investigation was to provide design information, and to define significant physical trends, which would aid in the prediction of turbulent separation for various aerodynamic devices, such as compressor blades, flaps, spoilers, and diffusers. A characteristic change in the longitudinal static-pressure distribution (i.e., a change from a curve with one inflection point to a curve with three inflection points) was employed to detect the occurrence of separation. The effects of Reynolds number (10(exp 6) to 10(exp 7) per foot or l.5 x 10(exp 4) to 7.5 x 10(exp 4) based upon boundary-layer thickness) and Mach number (1.6 to 4.2) on the onset of turbulent boundary-layer separation were investigated. The pressure gradient of the boundary-layer flow ahead of the interaction region was essentially zero. The results show a considerable effect of Mach number on the pressure rise for incipient separation for all configurations. For a curved-surface model, the static pressure-rise ratio required to cause separation varied from about 2.5, at a Mach number of 2 to about 16, at a Mach number of 3.5. A substantial effect of Reynolds number on the pressure rise for incipient separation was observed in the upper Mach number range and in the lower Reynolds number range; namely, the pressure rise required for separation decreased with increasing Reynolds number. For low Mach numbers and high Reynolds numbers, there appeared to be no Reynolds number effect. The effects of Mach number and of Reynolds number were similar for all models. Model shape was also found to be an important variable affecting the onset of separation. Large gains were realized in the pressure-rise ratio with no separation when the radius of curvature of the model surface was increased. At a Mach number of 3.4, for instance, the pressure-rise ratio with no separation increased from about 5 to 15 as a result of an increase in the radius of curvature from approximately 0 to 30 boundary-layer thicknesses.
The experimental database established by this investigation of the flow in a large rectangular turning duct is of benchmark quality. The experimental Reynolds numbers, Deans numbers and boundary layer characteristics are significantly different from previous benchmark curved-duct experimental parameters. This investigation extends the experimental database to higher Reynolds number and thinner entrance boundary layers. The 5% to 10% thick boundary layers, based on duct half-width, results in a large region of near-potential flow in the duct core surrounded by developing boundary layers with large crossflows. The turbulent entrance boundary layer case at R sub ed = 328,000 provides an incompressible flowfield which approaches real turbine blade cascade characteristics. The results of this investigation provide a challenging benchmark database for computational fluid dynamics code development.
Displacement thickness for two-dimensional or axisymmetric compressible boundary layer with blowing
Acoustic propagation through a lined two-dimensional duct is examined in order to assess the influence of the shape of the mean-velocity profile on the attenuation rate. Five mean-velocity profiles are considered: linear, parabolic, Pohlhausen, 1/7th power law with a linear sub-layer, and linear with slip at the wall. It is shown that when the attenuation rate is tabulated as a function of the boundary-layer thickness, as is usually done, substantially different results are obtained from the several mean profiles. However, when the displacement thickness is used, a considerable collapse is achieved in the attenuation curves that are obtained from the various profiles. For downstream propagation, all profiles produce essentially the same results over a reasonable range of boundary-layer thickness. However, for many cases of upstream propagation, the results from the 'turbulent' boundary-layer profiles differ significantly from the results of the other profiles even when compared on the basis of displacement thickness.
Preliminary compliant wall skin friction test results obtained in a low-turbulence pressure tunnel are reported. Compliant surface skins consisted of 0.0025 cm thick mylar, stretched under tension and area-bonded or longitudinally strip-bonded with silicone rubber adhesive to polyurethane foam. Mean velocity and fluctuating survey data were obtained with a single slanted hot wire. Mean velocity profiles for the area-bonded mylar surface skins indicated up to a 20% reduction in boundary-layer thickness (and lower momentum thickness) over rigid surfaces. This reduction in boundary-layer thickness indicates that a drag reduction occurred. In addition, a 16% reduction in wall shear was evident for the mylar/compressor foam compliant surface.
The transformation permits a uniform mesh to be used in the computational coordinate which extends across the layer. This coordinate transformation uses the local value of the skin friction coefficient to scale the thickness of the wall layer region, and the local maximum value of turbulent viscosity to scale the boundary-layer thickness. Results are presented for two dimensional boundary layers in both positive and negative pressure gradients and comparisons are made with experimental data and conventional variable-grid results for low speed turbulent boundary-layers. The cases chosen illustrate the capability of this new transformation to capture the boundary layer growth over the full extent of laminar, transitional, and turbulent flow with no grid adjustment as well as its ability to consistently enlarge the wall layer region for accurate shear stress representation. Results of mesh refinement studies using the new coordinate transformation are presented.
The flow conditions for which heating distributions were measured on corrugated surfaces and wavy walls in turbulent boundary layers are shown, along with the ratio of the displacement thickness to the roughness height versus the local edge Mach number for an equivalent smooth surface. The present data are seen to greatly extend the range of data available on corrugated surfaces in turbulent boundary layers. These data were obtained by testing fullscale corrugation roughened panels in the wall boundary layer of a supersonic and hypersonic wind tunnel. The experimental program used to obtain the data is described. The data are analyzed and correlated in terms of the pertinent flow and geometric parameters. The developed correlations are compared with the available thin boundary layer data, as well as with previously published correlation techniques.
Further improvements made to advance the current Boundary Layer Integral Matrix Procedure - Version J (BLIMPJ) containing previously modeled simplified calculation methods by accounting for condensed phase, thick boundary layer and free stream turbulence effects are discussed. The condensed phase effects were included through species composition effect considered via input to the code and through particle damping effect considered via a turbulence model. The thrust loss calculation procedure for thick boundary layer effects was improved and the optimization of net thrust with respect to nozzle length was performed. The effects of free stream turbulence were approximately modeled in the turbulence model.
Measurements have been made of the heat transfer through a turbulent boundary layer on a convexly curved isothermal wall and on a flat plate following the curved section. Data were taken for one free-stream velocity and two different ratios of boundary layer thickness to radius of curvature delta/R = 0.051 and delta/R = 0.077. Only small differences were observed in the distribution of heat transfer rates for the two boundary layer thicknesses tested, although differences were noted in the temperature distributions within the boundary layer