Meelan M. Choudhari

Secondary instability of crossflow vortices and swept-wing boundary-layer transition

Crossflow instability of a three-dimensional boundary layer is a common cause of transition in swept-wing flows. The boundary-layer flow modified by the presence of finite-amplitude crossflow modes is susceptible to high-frequency secondary instabilities, which are believed to harbinger the onset of transition. The role of secondary instability in transition prediction is theoretically examined for the recent swept-wing experimental data by Reibert et al . (1996). Exploiting the experimental observation that the underlying three-dimensional boundary layer is convectively unstable, non-linear parabolized stability equations are used to compute a new basic state for the secondary instability analysis based on a two-dimensional eigenvalue approach. The predicted evolution of stationary crossflow vortices is in close agreement with the experimental data. The suppression of naturally dominant crossflow modes by artificial roughness distribution at a subcritical spacing is also confirmed. The analysis reveals a number of secondary instability modes belonging to two basic families which, in some sense, are akin to the ‘horseshoe’ and ‘sinuous’ modes of the Gortler vortex problem. The frequency range of the secondary instability is consistent with that measured in earlier experiments by Kohama et al . (1991), as is the overall growth of the secondary instability mode prior to the onset of transition (e.g. Kohama et al . 1996). Results indicate that the N -factor correlation based on secondary instability growth rates may yield a more robust criterion for transition onset prediction in comparison with an absolute amplitude criterion that is based on primary instability alone.

Roughness-induced generation of crossflow vortices in three-dimensional boundary layers

The receptivity theory of Goldstein and Ruban is extended within the nonasymptotic (quasi-parallel) framework of Zavolskii et al. to predict the roughness-induced generation of stationary and nonstationary instability waves in three-dimensional, incompressible boundary layers. The influence of acoustic-wave orientation, as well as that of different types of roughness geometries, including isolated roughness elements, periodic arrays, and two-dimensional lattices of compact roughness shapes, as well as random, but spatially homogeneous roughness distributions, is examined. The parametric study for the Falkner-Skan-Cooke family of boundary layers supports our earlier conjecture that the initial amplitudes of roughness-induced stationary vortices are likely to be significantly larger than the amplitudes of similarly induced nonstationary vortices in the presence of acoustic disturbances in the free stream. Maximum unsteady receptivity occurs when the acoustic velocity fluctuation is aligned with the wave-number vector of the unsteady vortex mode. On the other hand, roughness arrays that are oriented somewhere close to the group velocity direction are likely to produce higher instability amplitudes. Limitations of the nonasymptotic theory are discussed, and future work is suggested.

Tandem Cylinder Noise Predictions

In an effort to better understand landing-gear noise sources, we have been examining a simplified configuration that still maintains some of the salient features of landing-gear flow fields. In particular, tandem cylinders have been studied because they model a variety of component level interactions. The present effort is directed at the case of two identical cylinders spatially separated in the streamwise direction by 3.7 diameters. Experimental measurements from the Basic Aerodynamic Research Tunnel (BART) and Quiet Flow Facility (QFF) at NASA Langley Research Center (LaRC) have provided steady surface pressures, detailed off-surface measurements of the flow field using Particle Image Velocimetry (PIV), hot-wire measurements in the wake of the rear cylinder, unsteady surface pressure data, and the radiated noise. The experiments were conducted at a Reynolds number of 166 105 based on the cylinder diameter. A trip was used on the upstream cylinder to insure a fully turbulent shedding process and simulate the effects of a high Reynolds number flow. The parallel computational effort uses the three-dimensional Navier-Stokes solver CFL3D with a hybrid, zonal turbulence model that turns off the turbulence production term everywhere except in a narrow ring surrounding solid surfaces. The current calculations further explore the influence of the grid resolution and spanwise extent on the flow and associated radiated noise. Extensive comparisons with the experimental data are used to assess the ability of the computations to simulate the details of the flow. The results show that the pressure fluctuations on the upstream cylinder, caused by vortex shedding, are smaller than those generated on the downstream cylinder by wake interaction. Consequently, the downstream cylinder dominates the noise radiation, producing an overall directivity pattern that is similar to that of an isolated cylinder. Only calculations based on the full length of the model span were able to capture the complete decay in the spanwise correlation, thereby producing reasonable noise radiation levels.

Transition Analysis for the HIFiRE-5 Vehicle

The Hypersonic International Flight Research and Experimentation (HIFiRE) 5 flight experiment by Air Force Research Laboratories and Australian Defense Science and Technology Organization is designed to provide in-flight boundary-layer transition data for a canonical 3D configuration at hypersonic Mach numbers. This paper outlines the progress, to date, on boundary layer stability analysis for the HIFiRE-5 flight configuration, as well as for selected test conditions from the wind tunnel experiments supporting the flight test. At flow conditions corresponding to the end of the test window, rather large values of linear amplification factor are predicted for both second mode (N>40) and crossflow (N>20) instabilities, strongly supporting the feasibility of first in-flight measurements of natural transition on a fully three-dimensional hypersonic configuration. Additional results highlight the rich mixture of instability mechanisms relevant to a large segment of the flight trajectory, as well as the effects of angle of attack and yaw angle on the predicted transition fronts for ground facility experiments at Mach 6. 1. Background

Characterization of Unsteady Flow Structures Near Leading-Edge Slat. Part 1; PIV Measurements

A comprehensive computational and experimental study has been performed at the NASA Langley Research Center as part of the Quiet Aircraft Technology (QAT) Program to investigate the unsteady flow near a leading-edge slat of a two-dimensional, high-lift system. This paper focuses on the experimental effort conducted in the NASA Langley Basic Aerodynamics Research Tunnel (BART) where Particle Image Velocimetry (PIV) data was acquired in the slat cove and at the slat trailing edge of a three-element, high-lift model at 4, 6, and 8 degrees angle of attack and a freestream Mach Number of 0.17. Instantaneous velocities obtained from PIV images are used to obtain mean and fluctuating components of velocity and vorticity. The data show the recirculation in the cove, reattachment of the shear layer on the slat lower surface, and discrete vortical structures within the shear layer emanating from the slat cusp and slat trailing edge. Detailed measurements are used to examine the shear layer formation at the slat cusp, vortex shedding at the slat trailing edge, and convection of vortical structures through the slat gap. Selected results are discussed and compared with unsteady, Reynolds-Averaged Navier-Stokes (URANS) computations for the same configuration in a companion paper by Khorrami, Choudhari, and Jenkins (2004). The experimental dataset provides essential flow-field information for the validation of near-field inputs to noise prediction tools.

Effect of Three-Dimensional Shear-Layer Structures on Slat Cove Unsteadiness

Numerical simulations are used to investigate the local and global dynamics of large-scale, three-dimensional vorticity structures within the free shear layer originating from the slat cusp of a multielement airfoil configuration. Results indicate that accounting for the local three-dimensionality of flow fluctuations leads to substantially improved agreement between the computed unsteady near-field solution and wind tunnel measurements based on particle image velocimetry. Analysis of simulation data indicates the potential significance of high intensity turbulent fluctuations near the reattachment location along the slat lower surface toward the generation of broadband slat noise. The computed acoustic characteristics, in terms of the frequency spectrum and spatial distribution within short distances from the slat, resemble the previously reported, subscale measurements of slat noise.

Noise Radiation from a Leading-Edge Slat

This paper extends our previous computations of unsteady flow within the slat cove region of a multi-element high-lift airfoil configuration, which showed that both statistical and structural aspects of the experimentally observed unsteady flow behavior can be captured via 3D simulations over a computational domain of narrow spanwise extent. Although such narrow domain simulation can account for the spanwise decorrelation of the slat cove fluctuations, the resulting database cannot be applied towards acoustic predictions of the slat without invoking additional approximations to synthesize the fluctuation field over the rest of the span. This deficiency is partially alleviated in the present work by increasing the spanwise extent of the computational domain from 37.3% of the slat chord to nearly 226% (i.e., 15% of the model span). The simulation database is used to verify consistency with previous computational results and, then, to develop predictions of the far-field noise radiation in conjunction with a frequency-domain Ffowcs Williams-Hawkings solver.

Unsteady flowfield around tandem cylinders as prototype component interaction in airframe noise

Synergistic application of experiments and numerical simulations is crucial to understanding the underlying physics of airframe noise sources. The current effort is aimed at characterizing the details of the flow interaction between two cylinders in a tandem configuration. This setup is viewed to be representative of several component-level flow interactions that occur when air flows over the main landing gear of large civil transports. Interactions of this type are likely to have a significant impact on the noise radiation associated with the aircraft undercarriage. The paper is focused on two-dimensional, time-accurate flow simulations for the tandem cylinder configuration. Results of the unsteady Reynolds Averaged Navier-Stokes (URANS) computations with a two-equation turbulence model, at a Reynolds number of 0.166 million and a Mach number of 0.166, are presented. The experimental measurements of the same flow field are discussed in a separate paper by Jenkins, Khorrami, Choudhari, and McGinley (2005). Two distinct flow regimes of interest, associated with short and intermediate separation distances between the two cylinders, are considered. Emphasis is placed on understanding both time averaged and unsteady flow features between the two cylinders and in the wake of the rear cylinder. Predicted mean flow quantities and vortex shedding frequencies show reasonable agreement with the measured data for both cylinder spacings. Computations for short separation distance indicate decay of flow unsteadiness with time, which is not unphysical; however, the predicted sensitivity of mean lift coefficient to small angles of attack explains the asymmetric flowfield observed during the experiments.

Characterization of Unsteady Flow Structures Around Tandem Cylinders for Component Interaction Studies in Airframe Noise

A joint computational and experimental study has been performed at NASA Langley Research Center to investigate the unsteady flow generated by the components of an aircraft landing gear system. Because the flow field surrounding a full landing gear is so complex, the study was conducted on a simplified geometry consisting of two cylinders in tandem arrangement to isolate and characterize the pertinent flow phenomena. This paper focuses on the experimental effort where surface pressures, 2-D Particle Image Velocimetry, and hot-wire anemometry were used to document the flow interaction around the two cylinders at a Reynolds Number of 1.66 x 10(exp 5), based on cylinder diameter, and cylinder spacing-todiameter ratios, L/D, of 1.435 and 3.70. Transition strips were applied to the forward cylinder to produce a turbulent boundary layer upstream of the flow separation. For these flow conditions and L/D ratios, surface pressures on both the forward and rear cylinders show the effects of L/D on flow symmetry, base pressure, and the location of flow separation and attachment. Mean velocities and instantaneous vorticity obtained from the PIV data are used to examine the flow structure between and aft of the cylinders. Shedding frequencies and spectra obtained using hot-wire anemometry are presented. These results are compared with unsteady, Reynolds-Averaged Navier-Stokes (URANS) computations for the same configuration in a companion paper by Khorrami, Choudhari, Jenkins, and McGinley (2005). The experimental dataset produced in this study provides information to better understand the mechanisms associated with component interaction noise, develop and validate time-accurate computer methods used to calculate the unsteady flow field, and assist in modeling of the radiated noise from landing gears.

Numerical simulations of boundary-layer bypass transition due to high-amplitude free-stream turbulence

Direct numerical simulations (DNS) of bypass transition due to high-amplitude free-stream turbulence (FST) are carried out for a flat-plate boundary layer. The computational domain begins upstream of the plate leading edge and extends into the fully turbulent region. Thus, there is no ad hoc treatment to account for the initial ingestion of FST into the laminar boundary layer. We study the effects of both the FST length scale and the disturbance behaviour near the plate leading edge on the details of bypass transition farther downstream. In one set of simulations, the FST parameters are chosen to match the ERCOFTAC benchmark case T3B. The inferred FST integral length scale L 11 is significantly larger ( R L = UL 11 /ν = 6580) than that employed in previous simulations of bypass transition ( R L ≃ 1000). An additional set of simulations was performed at R L = 1081 to compare the transition behaviour in the T3B case with that of a smaller value of FST length scale. The FST length scale is found to have a profound impact on the mechanism of transition. While streamwise streaks (Klebanoff modes) are observed at both values of the FST length scale, they appear to have clear dynamical significance only at the smaller value of R L , where transition is concomitant with streak breakdown. For the T3B case, turbulent spots form upstream of the region where streaks could be detected. Spot precursors are traced to quasi-periodic spanwise structures, first observed as short wavepackets in the wall-normal velocity component inside the boundary layer. These structures are reoriented to become horseshoe vortices, which break down into young turbulent spots. Two of the four spots examined for this case had a downstream-pointing shape, similar to those found in experimental studies of transitional boundary layers. Additionally, our simulations indicate the importance of leading-edge receptivity for the onset of transition. Specifically, higher fluctuations of the vertical velocity at the leading edge of the plate result in higher levels of streamwise Reynolds stress inside the developing boundary layer, facilitating breakdown to turbulence.