Note on a correlation of boundary-layer transition results on highly cooled blunt bodies
Boundary layer transition from hypersonic laminar to turbulent flow around blunt bodies
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.
Boundary layer transition from hypersonic laminar to turbulent flow around blunt bodies
An experiment is conducted on the transition from quasi-periodic to weakly turbulent flow of a fluid contained between concentric cylinders with the inner cylinder rotating and the outer cylinder at rest. Power spectra, phase-space portraits, and circle maps obtained from velocity time-series data indicate that the nonperiodic behavior observed is deterministic, that is, it is described by strange attractors. Various problems that arise in computing the dimension of strange attractors constructed from experimental data are discussed and it is shown that these problems impose severe requirements on the quantity and accuracy of data necessary for determining dimensions greater than about 5. In the present experiment the attractor dimension increases from 2 at the onset of turbulence to about 4 at a Reynolds number 50-percent above the onset of turbulence.
Hydrodynamic and semi-kinetic treatments of plasma flow along closed geomagnetic field lines are compared. The hydrodynamic treatment is based on a simplified 16-moment set of transport equations as the equations for the heat flows are not solved; the heat flows are treated heuristically. The semi-kinetic treatment is based on a particle code. The comparison deals with the distributions of the plasma density, flow velocity, and parallel and perpendicular temperatures as obtained from the two treatments during the various stages of the flow. In the kinetic treatment, the appropriate boundary condition is the prescription of the velocity distribution functions for the particles entering the flux tubes at the ionospheric boundaries; those particles leaving the system are determined by the processes occurring in the flux tube. The prescribed distributions are half-Maxwellian with temperature T(sub 0) and density n(sub 0). In the hydrodynamic model, the prescribed boundary conditions are on density (n(sub 0)), flow velocity (V(sub 0)) and temperature (T(sub 0). It was found that results from the hydrodynamic treatment critically depend on V(sub 0); for early stages of the flow this treatment yields results in good agreement with those from the kinetic treatment, when V(sub 0) = square root of (kT(sub 0)/2 (pi)m), which is the average velocity of particles moving in a given direction for a Maxwellian distribution. During this early stage, the flows developing form the conjugate ionospheres show some distinct transitions. For the first hour or so, the flows are highly supersonic and penetrate deep into the opposite hemispheres, and both hydrodynamics and kinetic treatments yield almost similar features. It is found that during this period heatflow effects are negligibly small. When a flow penetrates deep into the opposite hemisphere, the kinetic treatment predicts reflection and setting up of counterstreaming. In contrast, the hydrodynamic treatment yields a shock in the flow. The reasons for this difference in the two treatments is discussed, showing that in view of the relatively warm ions, the coupling of ion beams and the consequent shock formation in the offequatorial region are not likely due to the enhancements in the beam temperatures. The counterstreaming in the kinetic treatment and the shock in the hydrodynamic treatment first advance upward to the equator and then downward to the ionospheric boundary from where the flow originated. The transit time for this advancement is found to be about 1 hour for the respective models. After 2 hours or so, both models predict that the flows from the ionospheric boundaries are generally subsonic with respect to the local ion-sound speed. At late stages of the flow, when a substantial fraction of ions entering the flux tube begin to return back in the kinetic treatment, the hydrodynamic treatment with the boundary condition V(sub 0) = square root of (kT(sub 0)/2(pi)m) yields an over-refilling, and the choice of V(sub 0) becomes uncertain.
An overview of the NASA Glenn Low-Pressure Turbine (LPT) Flow Physics Program will be presented. The flow in the LPT is unique for the gas turbine. It is characterized by low Reynolds number and high freestream turbulence intensity and is dominated by interplay of three basic mechanisms: transition, separation and wake interaction. The flow of most interest is on the suction surface, where large losses are generated due to separation. The LPT is a large, multistage, heavy, jet engine component that suffers efficiency degradation between takeoff and cruise conditions due to decrease in Reynolds number with altitude. The performance penalty is around 2 points for large commercial bypass engines and as much as 7 points for small, high cruise altitude, military engines. The gas-turbine industry is very interested in improving the performance of the LPT and in reducing its weight, part count and cost. Many improvements can be accomplished by improved airfoil design, mainly by increasing the airfoil loading that can yield reduction of airfoils and improved performance. In addition, there is a strong interest in reducing the design cycle time and cost. Key enablers of the needed improvements are computational tools that can accurately predict LPT flows. Current CFD tools in use cannot yet satisfactorily predict the unsteady, transitional and separated flow in the LPT. The main reasons are inadequate transition & turbulence models and incomplete understanding of the LPT flow physics. NASA Glenn has established its LPT program to answer these needs. The main goal of the program is to develop and assess models for unsteady CFD of LPT flows. An approach that consists of complementing and augmenting experimental and computational work elements has been adopted. The work is performed in-house and by several academic institutions, in cooperation and interaction with industry. The program was reviewed at the Minnowbrook II meeting in 1997. This review will summarize the progress that was made since and will introduce newly started projects. The LPT program is focused on three areas: acquisition of experimental and numerical databases and on modeling and computation. Priority was initially given to experiments. There are three classes of experiments: simulated LPT passages, linear cascade, both with and without wakes, and low-speed rotating rig. They are being conducted as follows: At NASA GRC on a flat surface with blade pressure distribution, at the US Naval Academy on a curved surface. The addition of wakes is studied at the University of Minnesota in a curved passage with a retractable wake generator, and at Texas A&M University in a linear cascade with continuously running wake generator. The pressure distribution of the Pratt & Whitney blade 'Pak B' is used in all these experiments. Experiments have been performed also in the GEAE Low-Speed Rotating Turbine (LSRT) rig with GE-designed airfoils. Work on numerically generated database is in progress at the University of Kentucky, using the DNS/LES code LESTool developed there. Turbulence/transition model assessment and development is performed also at the University of Kentucky, where a new intermittency transport model was developed and many experimental test cases have been numerically computed. Assessments of models using simulations of multistage LPT experiments were performed at Virginia Commonwealth University using the Corsair code. Work on suction surface separation delay, using passive and active flow-control, has also been initiated. Following the overview, Principal Investigators attending the workshop will present in detail several of the projects supported by NASA.
Report presents the results of experimental and theoretical research conducted on flow separation associated with steps, bases, compression corners, curved surfaces, shock-wave boundary-layer reflections, and configurations producing leading-edge separation. Results were obtained from pressure-distribution measurements, shadowgraph observations, high-speed motion pictures, and oil-film studies. The maximum scope of measurement encompassed Mach numbers between 0.4 and 3.6, and length Reynolds numbers between 4,000 and 5,000,000.
Experiments were conducted on a 76 degree swept cylinder to establish the behavior of the attachment line transition process in a low-disturbance level, Mach number 1.6 flow. For a near adiabatic wall condition, the attachment-line boundary layer remained laminar up to the highest attainable Reynolds number. The attachment-line boundary layer transition under the influence of trip wires depended on wind tunnel disturbance level, and a transition onset condition for this flow is established. Internal heating raised the surface temperature of the attachment line to induce boundary layer instabilities. This was demonstrated experimentally for the first time and the frequencies of the most amplified disturbances were determined over a range of temperature settings. Results were in excellent agreement to those predicted by a linear stability code, and provide the first experimental verification of theory. Transition onset along the heated attachment line at an R-bar of 800 under quiet tunnel conditions was found to correlate with an N factor of 13.2. Increased tunnel disturbance levels caused the transition onset to occur at lower cylinder surface temperatures and was found to correlate with an approximate N factor of 1 1.9, so demonstrating that the attachment-line boundary layer is receptive to increases in the tunnel disturbance level.
The profile of turbulent intensity versus Reynolds number from bursts found in an ignition front agrees with the Orr-Sommerfeld solutions for unstable boundary-layer flow. This result provides the first evidence of a formal connection between bursting transition to turbulence in flows which share the boundary-layer approximation but which are otherwise unrelated.
Report describes computer-simulation of instability and transition to turbulence in plane Poisseuille flow in rotating channel. Motions of fluid governed by Navier-Stokes equations of incompressible flow, using periodic boundary conditions in streamwise and spanwise directions. Principal parameters affecting flow are Reynolds and rotation numbers.
IR thermography was used to characterize the transition front on a S-NLF test article at chord Reynolds numbers in excess of 30 million Changes in transition due to Mach number, Reynolds number, and surface roughness were investigated - Regions of laminar flow in excess of 80% chord at chord Reynolds numbers greater than 14 million IR thermography clearly showed the transition front and other flow features such as shock waves impinging upon the surface A series of parallel oblique shocks, of yet unknown origin, were found to cause premature transition at higher Reynolds numbers. NASA has a current goal to eliminate barriers to the development of practical supersonic transport aircraft Drag reduction through the use of supersonic natural laminar flow (S-NLF) is currently being explored as a means of increasing aerodynamic efficiency - Tradeoffs work best for business jet class at M<2 Conventional high-speed designs minimize inviscid drag at the expense of viscous drag - Existence of strong spanwise pressure gradient leads to crossflow (CF) while adverse chordwise pressure gradients amplifies and Tollmien-Schlichting (TS) instabilities Aerion Corporation has patented a S-NLF wing design (US Patent No. 5322242) - Low sweep to control CF - dp/dx < 0 on both wing surfaces to stabilize TS - Thin wing with sharp leading edge to minimize wave drag increase due to reduction in sweep NASA and Aerion have partnered to study S-NLF since 1999 Series of S-NLF experiments flown on the NASA F-15B research test bed airplane Infrared (IR) thermography used to characterize transition - Non-intrusive, global, good spatial resolution - Captures significant flow features well
Experimental results suggest that the boundary-layer mean profiles in hypersonic flows may be highly transitional in the outer part of the boundary layer before the transition process is detected at the surface. Means of determining the initial location where the transition process begins in the outer part of the boundary layer are considered.
The evolution of three-dimensional temporally evolving plane mixing layers through as many as three pairings has been simulated numerically. All simulations were begun from a few low-wavenumber disturbances, usually derived from linear stability theory, in addition to the mean velocity. Three-dimensional perturbations were used with amplitudes ranging from infinitesimal to large enough to trigger a rapid transition to turbulence. Pairing is found to inhibit the growth of infinitesimal three-dimensional disturbances, and to trigger the transition to turbulence in highly three-dimensional flows. The mechanisms responsible for the growth of three-dimensionality and onset of transition to turbulence are described. The transition to turbulence is accompanied by the formation of thin sheets of spanwise vorticity, which undergo secondary rollups. The post-transitional simulated flow fields exhibit many properties characteristic of turbulent flows.
This report documents a study of the use of liquid crystals to visualize boundary layer development on a turbomachine blade. A turbine blade model in a linear cascade of blades was used for the tests involved. Details of the boundary layer development on the suction surface of the turbine blade model were known from previous research. Temperature sensitive and shear sensitive liquid crystals were tried as visual agents. The temperature sensitive crystals were very effective in their ability to display the location of boundary layer flow separation and reattachment. Visualization of natural transition from laminar to turbulent boundary layer flow with the temperature sensitive crystals was possible but subtle. The visualization of separated flow reattachment with the shear sensitive crystals was easily accomplished when the crystals were allowed to make a transition from the focal-conic to a Grandjean texture. Visualization of flow reattachment based on the selective reflection properties of shear sensitive crystals was achieved only marginally because of the larger surface shear stress and shear stress gradient levels required for more dramatic color differences.
A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.
Quantitative observations of transitional boundary layers in regions of strong flow deceleration on an axial compressor stator blade are reported. Measurements are obtained at a fixed chordwise position, and the blade incidence was varied by changing the compressor throughflow so as to move the transition region relative to the stationary probe. It was thus possible to observe typical boundary layer behavior at various stages of transition in the turbomachine environment. The range of observations covers separating laminar flow at transition onset, and reattachment of intermittently turbulent periodically separated shear layers.
Recent progress in computational vortex flow aerodynamics at the Langley Research Center is reviewed. Emphasis is placed on Navier-Stokes methodology, both for compressible and incompressible flows, and results are presented from central and upwind-biased schemes for fully laminar, transitional, and fully turbulent flows. In addition, results are presented from selected potential-based methods to address the hierarchy of formulations presently available for the computational analysis of aerodynamic vortex flows. Some comparisons among this hierarchy of methods are shown.
A model is proposed for calculating laminar-turbulent transition in the boundary layer of flows in chemical shock tubes. It is shown that the time to turbulent transition may be treated as the time required by the shocked gas to move over one shock-tube diameter. The model is shown to provide a means of predicting the turbulent transition time in any shock tube experiment, provided the condition of limiting flow is met. Knowing the transition time, it becomes possible to obtain the density and temperature profiles and, hence, to determine precisely the effect of shock tube boundary layers on the measurement of chemical rate constants.
Aerobraking has been proposed as an efficient means of decelerating spacecraft for planetary missions. Most current aerobrake designs feature a blunt forebody shielding the payload from the intense heat generated during atmospheric entry. Although this forebody will absorb the largest portion of the heat pulse, accurate prediction of heating in the near wake is of great importance, since large local heating values can occur at points of shear-layer impingement. In order to address the various issues associated with these blunt-body wake flowfields, the Advisory Group for Aerospace Research and Development (AGARD) formed Working Group 18 in 1992. One of the objectives of this activity was to examine real-gas effects in high-speed flow fields around a 70 deg. blunted cone. To date, many researchers have conducted experiments using this geometry in various facilities, such as the Large Energy National Shock (LENS) tunnel at Cubric/Calspan and the HEG shock tunnel at DLR-Goettingen. Several computational studies have also been conducted in concert with these tests. Many of the experimental results have indicated the possible presence of a transitional shear layer through a large increase in heat transfer downstream of the reattachment point. The presence of transition could in fact lead to much higher peak heating than if the separated flow is entirely laminar or turbulent. In the shock-tunnel tests, however, it is difficult to separate such viscous-flow phenomena from real-gas effects. In order to help make this distinction, Horvath et al. recently conducted a set of experiments in the NASA Langley 20-Inch Mach 6 Tunnel, and compared the results to laminar Navier-Stokes calculations. They found heat-transfer distributions similar to those obtained in the high-enthalpy facilities, with the measured peak heating along the sting support markedly greater than that predicted by the laminar computations. These trends point to the need to find transitional and turbulent computational solutions for these flowfields.
This paper reports results from pressure tests on an ogive-cylinder in the low-turbulence 12-foot pressure wind tunnel. The results consist of pressure distributions over a wide range of Reynolds numbers and angles of attack. The tests encompassed a complete coverage of different roll orientations shown to be essential in order to fully define all the possible flow conditions. When the various roll-angle results are combined, it is possible to interpret the effects of changing angle of attack or Reynolds number. Two basic mechanisms for producing asymmetric flow are identified. One mechanism operates in both the laminar and the fully turbulent separation regimes; this mechanism is the one qualitatively described by the impulsive flow analogy. The other mechanism occurs only in the transitional separation regime. This asymmetric flow has the same form as that in the two-dimensional crossflow on a circular cylinder in the transitional flow regime. Finally, these results make it possible to draw up critical Reynolds number boundaries between the laminar, transitional, and fully turbulent separation regimes throughout the angle-of-attack range from 20 to 90 deg.