BOLT: Streak Instability Analysis (Including AoA and Yaw Effects)
Explore the source record for details and available documents.
boundary layer transition↗
Engineering topics
Meelan Choudhari
Publications and source records attributed to Meelan Choudhari.
Assessment of RANS-based Transition Models based on Experimental Data of the Common Research Model with Natural Laminar Flow
Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework often rely upon transition correlations that were derived from a limited number of low-speed experiments and these models often fail to account for all of the relevant transition mechanisms and/or the variation in those mechanisms with respect to changes in the significant flow parameters. Available data from a recent experiment on the Common Research Model with Natural Laminar Flow (CRM-NLF) in the National Transonic Facility at the NASA Langley Research Center are used to assess the current transition modeling capability in NASA's OVERFLOW 2.3b solver for a swept wing configuration with nonzero taper and transonic cruise conditions. Specifically, the OVERFLOW solutions are used to evaluate the accuracy and robustness of the transport-equation-based transition models. Results highlight that the Spalart-Allmaras-based amplification factor transport (AFT-2017b) equation model and Menter’s shear-stress transport equation (SST2003)-based Langtry-Menter transition models (either with or without the modeling of crossflow transition) significantly underpredict the reported extent of laminar flow region over the entire span of the wing, irrespective of which instability mechanism(s) is expected to dominate the onset of the transition process. We show that the transition correlations underlying these models fail to account for the stabilizing effect of compressibility on the Tollmien-Schlichting transition, which is likely to be a major contributor to the underprediction of the laminar flow region on the CRM-NLF. The SST-2003-based Langtry-Menter model also appears to inaccurately predict the chordwise pressure variation along the majority of the wing span at all the flow conditions studied herein, due to how the turbulence intensity levels were enforced in the computations and how that was interfering with the functioning of the underlying turbulence model within the boundary layer. The AFT and Langtry-Menter models appear to be sensitive to the level of the freestream turbulence intensity, but the degree of sensitivity varies across the models.
CFD modeling↗
A Computational Analysis of Boundary Layer Instability over the BOLT Configuration (Effect of Nonzero AoA and Yaw)
The complex boundary layer flow over the BOLT flight configuration is known to exhibit multiple and potentially interacting instability mechanisms. This paper represents a continuation of our numerical investigation of the flow instabilities over the main test surface of the BOLT configuration by using state-of-the-art tools in multidimensional instability analysis. Specifically, the paper extends our previous computations by considering the separate effects of a nonzero angle of attack and a nonzero yaw on the modal instability characteristics of the boundary layer streaks adjacent to the minor-axis symmetry plane, specifically near both ends of the azimuthal region of a thick boundary layer in the middle, where this region rapidly changes to a thinner boundary layer on either side. At the t = 28.8767 s condition from the ascent part of the anticipated flight trajectory with a flight Mach number of M ∞ = 5.53 and unit Reynolds number of 4.25 x 10 6 /m, either type of departure from the design condition is shown to have a considerable impact on the structure of the basic state rollup within the region of interest and, hence, also on the amplification characteristics of instability waves within the resulting streaks. Yet, for a yaw angle of β = 4 degrees, the computations indicate only a slight reduction with respect to the peak N-factor of nearly 11 at the design condition of zero degrees yaw and zero degrees angle of attack. In contrast, an angle of attack equal to α = 4 degrees, the peak N-factor decreases to nearly 6 on the leeward side and increases above 16 on the windward side, making the onset of transition highly likely on the windward side. Computations also highlight the role of streak instabilities that originate as Mack mode disturbances and also demonstrate the potential pitfalls in using surface pressure sensors alone to gauge the magnitude of instability amplification.
Boundary layer transition↗
Hypersonic Second-Mode Instability Response to Shaped Roughness
An experimental campaign was conducted on a 7-degree half-angle cone in the NASA Langley Research Center 20-Inch Mach 6 Wind Tunnel to examine the influence of arrays of regularly spaced roughness elements on instability growth and transition. The primary element shape was a pair of elliptical planform ramps that were inclined at equal and opposite angles with respect to the local streamwise direction. The element shapes were designed to induce transient growth disturbances that would lead to sustained azimuthal modulation of the boundary layer flow while limiting the nearfield disturbances to avoid an immediate, i.e., effective tripping of the boundary layer. The bulk of the run matrix consisted of testing different element height sat free stream unit Reynolds numbers ranging from 9.8 to 13.1 million per meter. Other element shapes previously designed for tripping hypersonic boundary layers were also implemented. The model was instrumented with surface mounted Kulite ® and PCB ® pressure transducers and thermocouples. Spectra from the PCBs ® indicated clear suppression of the second-mode instability; however, neither the PCB ® spectra nor the heat transfer data presented strong evidence for delayed turbulent flow. Complementary stability computations likewise demonstrated second-mode reduction, particularly just downstream of the roughness, but also revealed a rise in first mode (streak-instability) amplitudes from the baseline that was likely responsible for the earlier transition observed for taller roughness cases.
Hypersonics↗
Assessment and Improvement of RANS-based Transition Models based on Experimental Data of the Common Research Model with Natural Laminar Flow
Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework often rely upon the correlations that were derived from a limited number of low-speed experiments and do not account for all of the transition mechanisms and/or their variation with the significant flow parameters. Available data from a recent experiment in the National Transonic Facility at the NASA Langley Research Center are used to assess the current transition modeling capability in NASA's OVERFLOW 2.2o code for a swept wing configuration at transonic cruise conditions. Specifically, the OVERFLOW solutions are used together with detailed stability analysis of the boundary layer flow over the new Common Research Model with Natural Laminar Flow (CRM-NLF) to evaluate the accuracy and the robustness of the transport-equation-based transition models, with the goal of proposing improvements that would help to strengthen the physical basis of these models for the important class of flows involving the combined effects of crossflow and flow compressibility. Results highlight the significant underprediction of the laminar flow extent within the inboard region of the wing, wherein the onset of transition may be attributed to a gradual amplification of Tollmien-Schlichting instabilities.
Boundary layer transition↗
A Computational Analysis of Boundary Layer Instability over the BOLT Configuration
The complex boundary layer flow over the BOLT flight configuration is known to exhibit multiple and potentially interacting instability mechanisms. This paper represents a continuation of our numerical investigation of the flow instabilities over the main test surface of the BOLT configuration by using state-of-the-art tools in multidimensional instability analysis. Specifically, the paper extends our previous computations by considering the separate effects of a nonzero angle-of-attack and a nonzero yaw on the modal instability characteristics of the boundary layer streaks adjacent to the minor-axis symmetry plane. At a selected condition from the ascent part of the anticipated flight trajectory, either type of departure from the design condition is shown to have a considerable impact on the structure of the laminar basic state and, hence, also on the amplification characteristics of the instability waves. A broader set of computations, including N-factor predictions for the relevant instabilities, will be presented in the final paper.
Boundary layer transition↗
A Computational Analysis of Boundary Layer Instability over the BOLT Configuration
The complex boundary layer flow over the BOLT flight configuration is known to exhibit multiple and potentially interacting instability mechanisms. This paper represents a continuation of our numerical investigation of the flow instabilities over the main test surface of the BOLT configuration by using state-of-the-art tools in multidimensional instability analysis. Specifically, the paper extends our previous computations by considering the separate effects of a nonzero angle of attack and a nonzero yaw on the modal instability characteristics of the boundary layer streaks adjacent to the minor-axis symmetry plane, specifically near both ends of the azimuthal region of a thick boundary layer in the middle, where this region rapidly changes to a thinner boundary layer on either side. At the t = 28.8767 s condition from the ascent part of the anticipated flight trajectory with a flight Mach number of M ∞ = 5.53 and unit Reynolds number of 4.25 x 10 6 /m, either type of departure from the design condition is shown to have a considerable impact on the structure of the basic state rollup within the region of interest and, hence, also on the amplification characteristics of instability waves within the resulting streaks. Yet, for a yaw angle of β = 4 degrees, the computations indicate only a slight reduction with respect to the peak N-factor of nearly 11 at the design condition of zero degrees yaw and zero degrees angle of attack. In contrast, an angle of attack equal to α = 4 degrees, the peak N-factor decreases to nearly 6 on the leeward side and increases above 16 on the windward side, making the onset of transition highly likely on the windward side. Computations also highlight the role of streak instabilities that originate as Mack mode disturbances and also demonstrate the potential pitfalls in using surface pressure sensors alone to gauge the magnitude of instability amplification.
Boundary layer transition↗
Streak Instability Analysis for BOLT Configuration
The BOLT configuration designed for an upcoming flight experiment exemplifies the complexities of hypersonic boundary layer transition in the presence of multiple and potentially interacting instability mechanisms. A numerical study of the boundary layer instabilities over this configuration has been performed with the goal of applying state-of-the-art tools in multidimensional stability analysis that complement the other ongoing research related to the BOLT flight campaign. Specifically, the paper outlines modal instability characteristics of the boundary layer streaks adjacent to the minor-axis symmetry plane of a 1/3 rd scale BOLT configuration that was tested in the low-disturbance hypersonic wind tunnel at Purdue University and both quiet and conventional tunnels at Texas A&M University (TAMU). The basic state computations agree well with the heat flux measurements in the Purdue quiet tunnel as well as with the off-body velocity field described in a previous DNS study at the University of Minnesota. The encouraging agreement between the predicted frequency spectra of the streak instability in the region of increased boundary layer thickness and the surface pressure spectra measured during the TAMU experiments suggests that these streak instabilities may play an important role during the transition process near the symmetry plane. The modest value of the predicted peak N-factor (N max ≈ 5) for the streak instabilities provides a potential explanation for why transition is not observed under quiet flow conditions, but may occur in a noisy wind tunnel. Similar instability analyses are performed for the flight configuration at two selected flow conditions from the ascent portion of the planned flight trajectory. The predicted values of the peak N-factor suggest that transition onset in the vicinity of the center plane is possible at the ascent time of t = 26.98 seconds (N max ≈ 12.5 at X/L = 0.77), but the streak instabilities weaken sufficiently at the later time of t = 28.88 seconds (N max ≈ 11.5 at X/L ≈ 1.00) and are less likely to cause transition at that time.
Laminar turbulent transition↗
Streak Instabilities on HIFiRE-5 Elliptic Cone
The Hypersonic International Flight Research and Experimentation (HIFiRE) 5 flight experiment by the Air Force Research Laboratories and Australian Defense Science and Technology Organization yielded inflight, hypersonic transition data for a canonical, 3D configuration in the form of an elliptic cone at various combinations of angle-of-attack and yaw. Previous analysis has shown that boundary-layer stability calculations for the as-flown, zero-yaw flight conditions are able to capture the second mode dominated transition along the attachment line as well as the crossflow induced transition in between the attachment line and the minor axis symmetry plane. However, there is still a need to understand the transition in the vicinity of the minor axis plane of symmetry, where the boundary layer thickness exhibits an azimuthally compact bulge due to the collision of the secondary flow approaching from the attachment lines on both sides. This paper reports a computational analysis of the stability of the boundary layer flow along the bulge regions on both the leeward and windward sides of the vehicle surface. At a representative flight condition corresponding to a time instant of t = 514.813 seconds, computations are performed to identify several different unstable modes within the bulge regions, with the eventual goal of developing a linear stability correlation for the transition onset near the minor axis symmetry plane. The most amplified flow instabilities in the pretransitional region correspond to a “stem” mode that is concentrated within a pair of nearly vertical shear layers below the cap of the mushroom shaped bulge in the mean velocity contours. An excellent agreement is found between the predictions of instability amplification based on the plane marching parabolized stability equations (PSE) and the direct numerical solutions (DNS). Given the large number of unstable modes within the bulge region, the observed agreement demonstrates, for first time, that the PSE can provide accurate predictions for the development of flow instabilities in streak-like flow features over three-dimensional bodies. The peak N-factor of approximately 15 at the transition location estimated from the thermocouple measurements along the leeward centerline is close to the transition N-factor for Mack mode transition during the HIFiRE-1 circular cone flight experiment with nearly zero degrees angle of attack.
Meelan Choudhari↗
Assessment of Transition Modeling Capability in OVERFLOW with Emphasis on Swept-Wing Configurations
In preparation for comparisons with data obtained from the recently concluded experiments in the National Transonic Facility at the NASA Langley Research Center on the new common research model with natural laminar flow (CRM-NLF), an assessment of the current transition modeling capability in the NASA OVERFLOW 2.2o code has been carried out. A combination of the available experimental data and linear stability analysis is used to evaluate the accuracy and robustness of these models for selected swept-wing type configurations, with significant crossflow. An additional goal for this work involves providing a comparative assessment of the relevant transition models in the context of a single flow solver and identifying model limitations as well as the potential for future improvements that would help strengthen the physical basis of such transition models. Included in this investigation is an assessment of the sensitivities of the underlying transition models to grid resolution (wall-normal, as well as streamwise and spanwise) and the values of extra input parameters such as the level of surface roughness and freestream turbulence variables. The flow configurations targeted in this assessment include the NASA NLF(2)-0415 swept-wing configuration, the sickle-shaped wing introduced by the Technical University of Braunschweig, and the wing-body configuration of the CRM model from the fourth and fifth AIAA CFD Drag Prediction Workshops.
Freestream velocity↗
Effects of 3D Roughness Patch on Transition in High-Speed Boundary Layers
No abstract available
Meelan Choudhari↗
Boundary Layer Transition Induced by Distributed Roughness Array
The effects of a finite, spanwise-periodic array of cylindrical roughness elements on boundary layer transition over a NACA 0012 airfoil are investigated at a chord-based Reynolds number of 1.44×105 by using hotwire anemometry and infrared thermography. Both the number and the spanwise spacing of roughness elements in the array are varied in order to study their effect on the wake flow topology. Spanwise interaction between the roughness elements has an effect on the connection and the merging of neighbouring low-speed regions, which results in the formation of merged low-speed blobs (MLSs) that modify the spatial distribution and the amplitudes of the velocity streaks. When the spanwise distance between adjacent roughness elements equals 1.5 times the cylinder diameter, the transition location moves rapidly upstream. In this case, the two neighbouring low-speed regions overlap with each other in the near wake of the roughness, leading to the maximum growth in the velocity streak amplitude and the velocity fluctuations. The number of roughness elements affects the total number of MLSs within the boundary layer. For a single MLS behind a pair of cylinders, the Kelvin-Helmholtz instability dominates the growth of velocity fluctuations around the three-dimensional shear layers. When three cylinders are placed in the array, two MLSs appear in the near wake, which coalesce in to one low-speed blob downstream before the onset of transition, revealing the importance of Kelvin-Helmholtz instability.
Qingqing Ye↗
Nonlinear Görtler Vortices and Their Secondary Instability in a Hypersonic Boundary Layer
Nonlinear development of the Görtler instability over a concave surface gives rise to a highly distorted inflectional flow field in the boundary layer that leads to both wall-normal and spanwise gradients in the flow. Such nonlinear structures are susceptible to strong, high-frequency secondary instabilities that may lead to the onset of laminar-turbulent transition. The present numerical study uses direct numerical simulations and linear secondary instability theory to investigate finite amplitude Görtler vortices and their secondary instability characteristics, respectively, in the hypersonic flow over an axisymmetric cone with a concave aft body. To complement previous studies in the literature wherein the Görtler instability was usually studied for a flat plate and initiated at some upstream location by imposing an eigenfunction as the inflow condition or by blowing and suction at the wall, the present investigation is focused on fully realizable Görtler instability that is excited by an azimuthally periodic array of surface protuberances. Furthermore, while the previous work had mostly focused on the secondary instability of Görtler vortices with cross-plane velocity contours that resembled bell-shaped structures, the present results confirm that fully developed mushroom structures also exist in the hypersonic regime when the Görtler vortex amplitude is sufficiently large. Computations further reveal that the dominant modes of secondary instability in these mushroom-shaped structures correspond to an antisymmetic (i.e., sinuous) “stem” mode that concentrates within the strong, nearly wall-normal internal shear layers surrounding the stem regions underneath the caps of the mushroom structures. Additionally, there exist a multitude of other significantly unstable secondary instability modes of both symmetric and antisymmetric types. Analogous to the secondary instability of crossflow vortices in hypersonic flows, secondary instability modes of both symmetric and antisymmetric types. Analogous to the secondary instability of crossflow vortices in hypersonic flows, secondary instability modes originating from the Mack mode instability play an important role during the nonlinear breakdown process.
Laminar turbulent transition↗
Effect of 3D Roughness Patch on Instability Amplification in a Supersonic Boundary Layer
Surface roughness is known to have a substantial impact on the aerothermodynamic loading of high-speed vehicles, particularly via its influence on the laminar-turbulent transition process within the boundary layer. Numerical simulations are performed to investigate the effects of a distributed region of densely packed, sinusoidal shape roughness elements on a Mach 3.5 flat plate boundary layer for flow conditions corresponding to the planned conditions of an upcoming experiment in the Mach 3.5 Supersonic Low Disturbance Tunnel at the NASA Langley Research Center. Analysis of convective instabilities in the wake of the roughness patch was reported in a previous paper and the current work extends that analysis to instability amplification across the length of the roughness patch. Quasi-parallel stability analysis of the modified boundary layer flow over the patch indicates two dominant families of unstable disturbances, namely, a group of high frequency modes that amplify in localized regions along the roughness patch and another group of low frequency modes that have smaller peak amplification rates but amplify steadily both above the roughness patch and in the wake region behind it. The results suggest that the amplification factors associated with the high-frequency modes are sufficiently low, at least for the roughness patches considered in this paper, so that these modes are unlikely to have a major influence on the transition process. The amplification of the low-frequency modes within the region of the roughness patch is further quantified via direct numerical simulations. Results confirm the strongly destabilizing influence of the roughness patch on the first mode instabilities, yielding an N-factor increment of ΔN ≈ 3.6 for a roughness patch length of eight wavelengths.
Laminar turbulent transition↗
Görtler Instability and Its Control via Surface Suction over an Axisymmetric Cone at Mach 6
The characteristics of Görtler instability over an axisymmetric cone with an aft concave section are studied via linear and nonlinear instability analysis and direct numerical simulations. Several options for the cone geometry have been investigated numerically, subject to a fixed forecone section and constraints on the maximum cone diameter, overall cone length, and minimum N-factor for the most amplified Görtler modes. Computations show that it is possible to design a cone with a peak N-factor of Nmax > 8 at the target Reynolds number of 12.1×10(exp 6) per meter, corresponding to the maximum quiet Reynolds number in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. Direct numerical simulations show that an array of roughness elements corresponding to a peak roughness height of 0.1006 mm at the center can excite Görtler vortices that evolve into sufficiently strong streamwise streaks that may break down via high-frequency secondary instability. Thus, the selected axisymmetric configuration of interest should provide an acceptable baseline to investigate the feasibility of several aspects of laminar-flow control via boundary-layer suction. The apparatus that is being used for measurements in the Boeing/AFOSR Quiet Tunnel is described, along with some preliminary experimental results.
Gortler vortices↗
Effect of Distributed Patch of Smooth Roughness Elements on Transition in a High-Speed Boundary Layer
Surface roughness is known to have a substantial impact on the aerothermodynamic loading of hypersonic vehicles, particularly via its influence on the laminar-turbulent transition process within the boundary layer. Numerical simulations are performed to investigate the effects of a distributed region of densely packed, smooth-shaped roughness elements on the laminar boundary layer over a 7-degree half-angle, circular cone for flow conditions corresponding to a selected trajectory point from the ascent phase of the HIFiRE-1 flight experiment. For peak-to-valley roughness heights of 50 percent or less in comparison with the thickness of the unperturbed boundary layer, the computations converge to a stationary flow, suggesting that the flow is globally stable. Analysis of convective instabilities in the wake of the roughness patch indicates two dominant families of unstable disturbances, namely, a high frequency mode that corresponds to Mack mode waves modified by the wake and a lower frequency mode that corresponds to shear layer instabilities associated with the streaks in the roughness wake. Even though the peak growth rate of the later mode is more than 35 percent greater than the peak growth rates of the Mack modes, the latter modes achieve higher amplification ratios, and hence, are likely to dominate the onset of transition, which is estimated to occur slightly later than that in the unperturbed, i.e., smooth surface boundary layer. Additional computations are performed to investigate the effects of various roughness patch configurations on a Mach 3.5 flat plate boundary layer, to help guide an upcoming experiment in the Mach 3.5 Supersonic Low Disturbance Tunnel at NASA Langley Research Center. In this case, the cumulative reinforcement of basic state distortion over the length of the roughness patch is predicted to yield a significantly earlier transition than that over a smooth plate or a plate with a shorter length roughness patch.
Meelan Choudhari↗
Nonlinear Evolution and Breakdown of Azimuthally Compact Crossflow Vortex Pattern over a Yawed Cone
Hypersonic boundary-layer flows over a circular cone at moderate angle of incidence can support strong crossflow instability in between the windward and leeward rays on the plane of symmetry. Due to the more efficient excitation of stationary crossflow vortices by surface roughness, a possible path to transition in such flows corresponds to rapid amplification of the high-frequency secondary instabilities of finite amplitude stationary crossflow vortices. In the present paper, the previous analyses of crossflow instability over a 7- degree half-angle, yawed circular cone in a Mach 6 free stream have been extended to the nonlinear evolution of azimuthally localized crossflow vortex packets and the amplification characteristics and nonlinear breakdown of high-frequency secondary instabilities associated with those packets. A comparison between plane marching PSE and direct Navier-Stokes simulations (DNS) reveals favorable agreement in regard to mode shapes, most amplified disturbance frequencies, and N-factor evolution. In contrast, the quasi-parallel predictions are found to result in severe underprediction of the N-factors. The direct numerical simulations also indicate that the breakdown of secondary instabilities in a 3D hypersonic boundary layer shares certain common features with the previous computations of crossflow transition over subsonic swept wings.
Vortex structure↗
TPSAS-NF1676L-26887-DND
Explore the source record for details and available documents.
Meelan Choudhari↗