Search NASASearch

Engineering topics

Meelan Choudhari

Publications and source records attributed to Meelan Choudhari.

At least 19 records

Evaluation of Transport-Equation-Based Transition Models for High-Speed Boundary Layers Using OVERFLOW

Accurate modeling of laminar-turbulent transition is crucial for the design of hypersonic flight systems. However, the current transition models used in production CFD codes are insufficient for high-speed flows. Many extensions to low-speed models have been suggested; however, a thorough verification and validation effort is needed before these models can be used in design settings. Challenges include potentially missing details of the model implementation requirements and/or a complete specification of the input parameters needed to replicate the test findings. A meaningful assessment of the generalization capability of these models is also hindered by a lack of information regarding the specific flow configurations and associated grids employed for model calibration. As a key first step toward model verification, we present an independent assessment of two recently proposed models for high-speed transition, namely, a model within the SST-$\gamma$ framework and a model based on the SST-$\gamma-\nu_L$ equations. These models are implemented in the NASA OVERFLOW 2.3e solver and their performance in predicting first mode, second mode, and crossflow transition has been evaluated for several test cases in the supersonic and hypersonic regimes. Besides the test cases employed by the model developers, which could have also been used for model calibration, the present assessment includes supplementary configurations that contribute to an unbiased assessment of the models. The outcomes presented in this study indicate the potential for the models to be applied to high-speed flight configurations. Key steps toward future improvements to these models are also outlined.

High-speed flow

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

boundary layer transition

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

Boundary layer transition

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

boundary layer transition

CFD Integrated Transition Modeling for High-Speed Flows via Coupled OVERFLOW-LASTRAC Analysis

This work details ongoing efforts at the NASA Langley Research Center to develop and validate a general-use CFD tool that includes built-in predictions of boundary-layer transition in high-speed flows. Existing tools for the coupling of the NASA OVERFLOW structured overset RANS solver and LASTRAC stability analysis code have been extended to capture boundary-layer transition in high-speed flows driven by either Mack’s first-mode or secondmode instabilities and/or crossflow instabilities. The efficacy of this coupled approach has been demonstrated by examining a variety of supersonic and hypersonic test cases for which experimental validation data is available. The intermittency prescription parameters inherited from prior low-speed applications are found to be suitable for some, but not all high-speed flow scenarios. The method is found to be robust for straight cone configurations and improvements necessary to accurately capture transition on geometries with spatially complex amplification factor envelopes with intermediate regions of slow N-factor variation are examined. The automated, coupled analysis is demonstrated for a 3d supersonic test case and found to perform well within the limits of the linear stability analysis on which it relies. Finally, a preliminary investigation of the method’s robustness to non-ideal CFD meshes is conducted.

High-speed flow

Nonlinear Evolution of Instability in a Laminar Separation Bubble at Hypersonic Mach Number

The development of both convective stationary perturbation as well as global instabilities in the vicinity of a laminar separation bubble above an axisymmetric compression corner in a hypersonic flow is investigated using numerical simulations. The flow configuration of primary interest corresponds to the cone-cylinder-flare model used in experimental measurements in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. For a flare angle of 10 degrees and a unit Reynolds number of 11.5 x 10 6 m -1 , their surface flow visualizations identified the presence of streamwise elongated thermal streaks near the reattachment position and the dominant azimuthal spacing between the streaks was determined to be approximately 10 degrees (i.e., an azimuthal mode number of 36). Previous linear stability analyses predicted that the amplification characteristics of small amplitude, unsteady, convective instabilities within this flow were consistent with the surface pressure fluctuations measured in the experiment. However, their accompanying analysis of global instabilities had found the separation bubble to be weakly unstable at the 10 degrees flare angle, with the most unstable global mode corresponding to a stationary disturbance with an azimuthal wavenumber of m ≈ 5, which was well below the measured wavenumber of m = 36. Besides confirming the presence of the global instability at these flow conditions, the present numerical simulations quantify the details of the stationary equilibrium state associated with the supercritical bifurcation resulting from the nonlinear saturation of the unstable global mode. Although velocity perturbations associated with the saturated global mode are dominated by the fundamental spanwise wavelength associated with the linear global instability, the surface heat flux downstream of reattachment is dominated by m = 36 in agreement with experimental measurements. However, the shorter wavelength peak ( m = 36) in the azimuthal spectrum is not manifested as streaks with an azimuthal spacing of 10 degrees. An additional simulation of the combined evolution of both unsteady convective and global instabilities indicates that the nonlinear interactions between a broadband spectrum of oblique Mack’s first mode instabilities generate stationary vortex modes with a higher dominant wavenumber ( m = 60) as that observed in the experiments. Additional analyses are performed to predict the global mode behavior at other flare angles and how it differs from the convective centrifugal instability due to concave streamline curvature over the compression corner.

High speed flow

Stability Analysis of Streaks Induced By Optimized Vortex Generators

Numerical computations are performed to investigate the potential for transition control in an axisymmetric boundary layer via fully realizable, streamwise stationary streaks induced by an azimuthally periodic array of surface mounted vortex generators (VGs). Previous work has shown that suitable streaks of this type can significantly reduce the growth of Mack’s second mode instabilities, but large streak amplitudes can make the flow susceptible to previously absent streak instabilities that can become the leading cause of transition. Here, we use the adjoint capabilities of the SU2 flow solver to optimize the VG shape to maximize the reduction in the growth of second-mode disturbances while also preventing the streak amplitudes from reaching large enough values to precipitate an earlier onset of transition via streak instabilities. The geometry and the freestream flow conditions are selected to match a relevant trajectory lo-cation from the HIFiRE-1 flight experiment. Results show that the optimized VGs can increase the mean streak amplitude by 117% with respect to a manually developed baseline design. The stability of this optimized basic state is analyzed via the plane-marching parabolized stability equations, predicting a fully laminar flow over the entire cone, or equivalently, yielding transition delay of 130% versus the 17% for the baseline VGs.

Boundary layer transition

Nonlinear Evolution of Instabilities in a Laminar Separation Bubble at a Hypersonic Mach Number

The development of both convective stationary perturbation as well as global instabilities in the vicinity of a laminar separation bubble above an axisymmetric compression corner in a hypersonic flow is investigated using numerical simulations. The flow configuration of primary interest corresponds to the cone-cylinder-flare model used in experimental measurements in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. For a flare angle of 10 degrees and a unit Reynolds number of 11.5 x 10 6 m -1 , their surface flow visualizations identified the presence of streamwise elongated thermal streaks near the reattachment position and the dominant azimuthal spacing between the streaks was determined to be approximately 10 degrees (i.e., an azimuthal mode number of 36). Previous linear stability analyses predicted that the amplification characteristics of small amplitude, unsteady, convective instabilities within this flow were consistent with the surface pressure fluctuations measured in the experiment. However, their accompanying analysis of global instabilities had found the separation bubble to be weakly unstable at the 10 degrees flare angle, with the most unstable global mode corresponding to a stationary disturbance with an azimuthal wavenumber of m ≈ 5, which was well below the measured wavenumber of m = 36. Besides confirming the presence of the global instability at these flow conditions, the present numerical simulations quantify the details of the stationary equilibrium state associated with the supercritical bifurcation resulting from the nonlinear saturation of the unstable global mode. Although velocity perturbations associated with the saturated global mode are dominated by the fundamental spanwise wavelength associated with the linear global instability, the surface heat flux downstream of reattachment is dominated by m = 36 in agreement with experimental measurements. However, the shorter wavelength peak ( m = 36) in the azimuthal spectrum is not manifested as streaks with an azimuthal spacing of 10 degrees. An additional simulation of the combined evolution of both unsteady convective and global instabilities indicates that the nonlinear interactions between a broadband spectrum of oblique Mack’s first mode instabilities generate stationary vortex modes with a higher dominant wavenumber ( m = 60) as that observed in the experiments. Additional analyses are performed to predict the global mode behavior at other flare angles and how it differs from the convective centrifugal instability due to concave streamline curvature over the compression corner.

Hypersonic flow

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling

Wake Instability Behind Isolated Trip Near the Leading Edge of the BOLT-II Configuration

The BOLT-II (Holden Mission) configuration is an extended version of the BOLT flight article and depicts hypersonic boundary-layer transition in the presence of multiple and potentially interacting instability mechanisms. Several numerical studies of the boundary-layer instabilities over these configurations have been reported in the recent literature, including our previous studies of the modal instability characteristics of boundary-layer streaks adjacent to the minor-axis symmetry plane of the BOLT configuration and the wake instabilities behind a diamond planform (“pizza-box”) trip along the symmetry plane on the secondary side of the BOLT-II configuration. The present work extends the latter study to a scaled version of the same trip that is located in the region of nonzero crossflow in the vicinity of the leading edge at X/L = 0.5. Collectively, the trips along the symmetry plane and near the leading edge are the focus of the NASA roughness experiment on the secondary side of the BOLT-II configuration. The laminar basic state computation at the nominal flight design condition of Re ∞ = 5.44 x 10 6 /m and Re ∞ = 2.5 x 10 6 /m shows that the leading-edge trip with k/δ ≈ 0.70 and planform-side-length-to-height ratio of b/k = 3.0 induces multiple asymmetric, longitudinal streaks within the trip wake. The most prominent streak among these resembles a finite amplitude crossflow vortex, and it supports the amplification of multiple families of unstable modes. The application of multi-dimensional instability analysis to the wake flow reveals that the most amplified unstable mode can achieve a peak N-factor of up to 20 by X/L = 0.80, indicating that the onset of transition is more than likely to occur within eighty percent of the model length. To the best of our knowledge, the present study represents the first analysis including nonparallel and curvature effects on hypersonic tripwake instabilities in the presence of boundary-layer crossflow over a three-dimensional configuration.

Boundary layer transition

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling

Hypersonic Boundary-Layer Transition on Blunted Cones at Angle of Attack

Experimental studies of cones at several high-speed facilities have demonstrated that, for small nosetip bluntness, transition onset over a circular cone moves upstream along the leeward side and downstream along the windward side, but this trend may be reversed at large bluntness values, where transition onset moves downstream along the leeward side and upstream along the windward side. A theoretical and numerical investigation is performed to characterize the effects of nose bluntness on disturbance amplification over the circular cone for several angles of attack, with the goal of understanding the potential physical mechanisms behind the experimental observations. The three-dimensional laminar basic states over a 1.5 m long, 7-degree half-angle cone with 9.525 mm nosetip radius are computed for selected angles of attack values and freestream conditions that are selected to match the Mach 10 experiments conducted within the Hypervelocity Wind Tunnel 9 at the Arnold Engineering Development Complex (AEDC). The solutions at a freestream unit Reynolds number of 17.1 million per meter are used to perform detailed instability analyses for angles of attack equal to 0, 1, 3, and 5 degrees. Results indicate that the linear amplification of stationary crossflow waves along inflection lines may begin to influence transition along the acreage of the cone for angles of attack equal to or larger than 5 degrees. The measured trend in transition front with respect to increasing angle of attack is found to be consistent with the predicted increase in the amplification factors for Mack mode disturbances along the streamline trajectories. The increase in Mack mode amplification along the windward ray for higher angles of attack is shown to be the result of a progressively earlier entropy-layer swallowing. Computations also indicate that the transition amplification factor along the windward ray is not constant and increases with the angle of attack and that the transition #-factors along the leeward ray are rather small. The nonmodal analysis for zero degrees angle of attack shows that entropy-layer disturbances with appreciably strong energy growth can coexist with Mack mode instabilities at the measured transition location.

boundary layer transition

Pretest Computational Assessment of Boundary Layer Transition in the NASA Juncture Flow Model with an NACA 0015-Based Wing

The first two phases of the NASA Juncture Flow experiment were carried out on a DLR-F6 swept-wing model and were designed to provide “CFD validation-quality” data toward the assessment and improvement of existing CFD turbulence models in predicting onset and extent of three-dimensional separated flow near the wing-juncture trailing-edge region. The next phase of experiments will involve an NACA 0015-based swept wing, as prior risk reduction experiments had indicated that this wing shape resulted in reduced separation near the juncture region than the DLR-F6 wing, thus providing a better option to evaluate the ability of CFD models to predict incipient turbulent separation. The NACA 0015 measurements will also include IR thermography to infer the variation of transition front with respect to an increasing angle of attack. The primary objective of this work is to computationally make a preliminary assessment of the transition front on both surfaces of the NACA 0015 wing at a crank-chord-based Reynolds number of 2.4 x 106 for four different angles of attack, (0°, 2.5°, 5°, and 7.5°) and to determine the dominant mechanisms responsible for transition. This assessment includes both RANS-based transition models from NASA’s OVERFLOW 2.3b flow solver and linear parabolized stability equations (PSE) stability analysis based on the Langley Stability and Transition Analysis code, LASTRAC. Linear PSE results indicate that the upper surface of the wing is dominated by Tollmien- Schlichting (TS) instabilities, and that the laminar flow region shrinks from about 50% chord to a very small region just downstream of the attachment line as the angle of attack is increased from 0° to 7.5°. Consequently, the transition fronts predicted by the Spalart- Allmaras-based amplification factor transport (AFT-2017b) equation model (which accounts for the TS instabilities alone) and the Menter’s shear-stress transport equation (SST2003)- based Langtry-Menter transition model with ability to account for both TS and crossflow effects (LM2015) compare well with those predicted using linear PSE. On the lower surface of the wing, stationary crossflow (CF) instabilities begin to appear on the inboard portion of the wing in addition to the TS-instabilities for the larger angles of attack (5° and 7.5°), further reducing the laminar flow extent within the inboard region. The LM2015 model that accounts for CF effects is able to replicate this trend but appears to predict a slightly earlier transition. The outcome of this effort will inform the experiment and, when the actual experimental data become available, provide further opportunity to assess and improve the various transition models.

CFD modeling

Transition Analysis for Isolated Trips on BOLT-II Wind-Tunnel and Flight Configuration

The BOLT-II 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. The present paper is related to the NASA-led roughness experiment on the secondary side of the BOLT-II configuration. First, computations are performed to aid mean-flow based correlations for the transition measurements acquired during an experimental campaign in the Ludwieg Tube facility at the United States Air Force Academy. To that end, both laminar and turbulent flow computations are performed for selected conditions of the experiment with a 1:3 scale model of the slightly shorter BOLT configuration. The variation in well-known correlation parameters for transition due to isolated roughness elements with respect to element position, height, and flow Reynolds number is delineated. Furthermore, the wake instability theory is used to examine the amplification of unstable disturbances behind the centerline trip at the design conditions of the BOLT-II flight experiment and for its dynamically scaled version on the 1:3 scale model of the wind tunnel experiment. To our knowledge, this study represents the first analysis of the wake instabilities in the presence of a prominent entropy layer outside the boundary layer. At the nominal flight design condition of Re ∞ ) = 5.44 and Re ∞ = 2.5 x 10 6 /m, a diamond planform (“pizza-box”) trip with k/δ = 0.70 and planform-halfwidth-to-height ratio of b/k = 3.0 at X/L = 0.254 induces a strong distortion of the boundary-layer flow in the vicinity of the symmetry plane, such that the wake structure includes a centerline streak produced by the flow retardation behind the trip along with multiple additional, prominent vorticity structures representing the necklace vortex system. Stability results reveal that each of these streaks can support a strong amplification of wake instabilities within a frequency band of 95 kHz–255 kHz, with a peak N-factor of greater than 15 beyond X/L ≈ 0.42. The dynamically scaled wind tunnel configuration with Re ∞ = 20.25 x 10 6 /m is also found to support N > 15 at similar axial stations.

BOLT

Pretest Computational Assessment of Boundary Layer Transition in the NASA Juncture Flow Model with an NACA 0015-Based Wing

The first two phases of the NASA Juncture Flow experiment were carried out on a DLR-F6 swept-wing model and were designed to provide “CFD validation-quality” data toward the assessment and improvement of existing CFD turbulence models in predicting onset and extent of three-dimensional separated flow near the wing-juncture trailing-edge region. The next phase of experiments will involve an NACA 0015-based swept wing, as prior risk reduction experiments had indicated that this wing shape resulted in reduced separation near the juncture region than the DLR-F6 wing, thus providing a better option to evaluate the ability of CFD models to predict incipient turbulent separation. The NACA 0015 measurements will also include IR thermography to infer the variation of transition front with respect to an increasing angle of attack. The primary objective of this work is to computationally make a preliminary assessment of the transition front on both surfaces of the NACA 0015 wing at a crank-chord-based Reynolds number of 2.4 x 106 for four different angles of attack, (0°, 2.5°, 5°, and 7.5°) and to determine the dominant mechanisms responsible for transition. This assessment includes both RANS-based transition models from NASA’s OVERFLOW 2.3b flow solver and linear parabolized stability equations (PSE) stability analysis based on the Langley Stability and Transition Analysis code, LASTRAC. Linear PSE results indicate that the upper surface of the wing is dominated by Tollmien- Schlichting (TS) instabilities, and that the laminar flow region shrinks from about 50% chord to a very small region just downstream of the attachment line as the angle of attack is increased from 0° to 7.5°. Consequently, the transition fronts predicted by the Spalart- Allmaras-based amplification factor transport (AFT-2017b) equation model (which accounts for the TS instabilities alone) and the Menter’s shear-stress transport equation (SST2003)- based Langtry-Menter transition model with ability to account for both TS and crossflow effects (LM2015) compare well with those predicted using linear PSE. On the lower surface of the wing, stationary crossflow (CF) instabilities begin to appear on the inboard portion of the wing in addition to the TS-instabilities for the larger angles of attack (5° and 7.5°), further reducing the laminar flow extent within the inboard region. The LM2015 model that accounts for CF effects is able to replicate this trend but appears to predict a slightly earlier transition. The outcome of this effort will inform the experiment and, when the actual experimental data become available, provide further opportunity to assess and improve the various transition models.

CFD modeling

Transition Analysis for Isolated Trips on BOLT-II Wind-Tunnel and Flight Configuration

The BOLT-II 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. The present paper is related to the NASA-led roughness experiment on the secondary side of the BOLT-II configuration. First, computations are performed to aid mean-flow based correlations for the transition measurements acquired during an experimental campaign in the Ludwieg Tube facility at the United States Air Force Academy. To that end, both laminar and turbulent flow computations are performed for selected conditions of the experiment with a 1:3 scale model of the slightly shorter BOLT configuration. The variation in well-known correlation parameters for transition due to isolated roughness elements with respect to element position, height, and flow Reynolds number is delineated. Furthermore, the wake instability theory is used to examine the amplification of unstable disturbances behind the centerline trip at the design conditions of the BOLT-II flight experiment and for its dynamically scaled version on the 1:3 scale model of the wind tunnel experiment. To our knowledge, this study represents the first analysis of the wake instabilities in the presence of a prominent entropy layer outside the boundary layer. At the nominal flight design condition of Re ∞ ) = 5.44 and Re ∞ = 2.5 x 10 6 /m, a diamond planform (“pizza-box”) trip with k/δ = 0.70 and planform-halfwidth-to-height ratio of b/k = 3.0 at X/L = 0.254 induces a strong distortion of the boundary-layer flow in the vicinity of the symmetry plane, such that the wake structure includes a centerline streak produced by the flow retardation behind the trip along with multiple additional, prominent vorticity structures representing the necklace vortex system. Stability results reveal that each of these streaks can support a strong amplification of wake instabilities within a frequency band of 95 kHz–255 kHz, with a peak N-factor of greater than 15 beyond X/L ≈ 0.42. The dynamically scaled wind tunnel configuration with Re ∞ = 20.25 x 10 6 /m is also found to support N > 15 at similar axial stations.

BOLT