Coupling of the FUN3D Unstructured Flow Solver and the LASTRAC Stability Code to Model Transition
Explore the source record for details and available documents.
Transition Airfoil Spheroid↗
Engineering topics
Pedro Paredes
Publications and source records attributed to Pedro Paredes.
Transition Modeling Based on the Dual N-factor Method for the CRM-NLF Wind Tunnel Configuration
The dual N-factor method is used to model the boundary-layer transition over the common research model with natural laminar flow (CRM-NLF) aircraft configuration. The flow conditions match selected test conditions from a wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. The paper presents a systematic methodology for transition prediction in the presence of a dual shock system and extends the prior capability for iteratively coupled computational fluid dynamics (CFD) predictions to incorporate three-dimensional, transonic wings. The method employs stability computations based on the linear parabolized stability equations (PSE), along with a dual N-factor criterion. The iterative process begins with the fully turbulent Reynolds-averaged-Navier-Stokes (RANS) mean flow solution. For the first iteration, a mean flow solution is calculated with an imposed transition front that aligns with the shock front from the fully turbulent solution. Subsequently, stability computations are performed along a set of streamlines across the wing to calculate the amplification of planar Tollmien-Schlichting (TS) and stationary crossflow (CF) modes. The transition criterion based on the dual N-factor method is used to infer the updated transition front and the process is successively repeated until convergence of the solution. Within three iterations, the predicted fronts for angles of attack of 1.45, 1.98, 2.46 and 2.94 degrees and a mean-aerodynamic-chord Reynolds number equal to 15 million, approach visual convergence in most regions of the studied cases, and the resulting predictions are in good agreement with the transition fronts deduced from measurements of temperature-sensitive paint. Even though surface pressure measurements based on fully-turbulent flow agree well with the measured pressure coefficient distributions, strong viscous-inviscid interaction effects cause significant shifts in the shock locations based on the imposed transition front, underscoring the intrusive nature of static pressure measurements using surface mounted ports on the CRM-NLF configuration.
Boundary Layer Transition↗
Nonlinear Nonmodal Analysis of Hypersonic Flow over Blunt Cones
The linear amplification of modal disturbances that lead to boundary-layer transition in two-dimensional/axisymmetric hypersonic configurations is strongly reduced by the presence of a blunt nosetip, and the mechanisms underlying the observed onset of transition over the cone frustum are currently unknown. Linear nonmodal analysis has shown that both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification for moderate to large nosetip bluntness. The present study extends the previous linear analysis by including the nonlinear effects. Specifically, the perturbation form of the 2D, harmonic Navier-Stokes equations (HNSE) are solved with a fully implicit formulation and the Newton-Raphson method. The increased number of degrees of freedom for the nonlinear system presents difficulties for solution strategies based on direct solution of the linearized system. Such difficulties are overcome by using the GMRES iterative method with a preconditioner corresponding to a simplified Jacobian without the cross derivative terms. The HNSE solver is verified by comparing with nonlinear parabolized stability equation (NPSE) results for the nonlinear evolution of planar waves in an incompressible Blasius boundary layer and in a Mach 6 flow over a blunt cone. Finally, nonlinear nonmodal results are presented for planar traveling disturbances over the blunt cone. The nonmodal analysis demonstrates that entropy-layer disturbances generated close to the nose tip can seed the amplification of higher frequency Mack’s second-mode instabilities further downstream.
boundary layer transition↗
Eigenvalue Sensitivity Computations for Linear Stability Theory
To realize the drag reduction benefit of boundary-layer transition control strategies, it is crucial to integrate transition prediction into the vehicle design through an optimization process. The integration of transition prediction based on linear stability analysis into adjoint d design optimization requires coupling an adjoint enabled computational fluid dynamics (CFD) solver with an adjoint enabled linear stability code. In particular, the boundary-layer transition location is often predicted using the N-factor method based on linear stability theory (LST). Thus, sensitivity of the linear-stability eigenvalues constitute an essential building block for optimizing the laminar flow performance. The present paper describes an implementation of LST eigenvalue sensitivity analysis that can be easily coupled with a CFD solver. Specifically, we describe a discrete adjoint formulation for the transition location prediction based on the N-factor method. The verification of this formulation is carried out by comparing the adjoint-based sensitivity of the local growth rate of a given instability mode with respect to the disturbance frequency, and the adjoint-based sensitivity of the transition location with respect to spanwise wavenumber with those sensitivities computed using a finite-difference approximation. Finally, the adjoint LST formulation is applied to flat-plate boundary-layer flows at transonic, supersonic, and hypersonic conditions, to determine the behavior and sensitivities of the transition location with respect to a range of disturbance spanwise wavenumbers.
Boundary Layer Transition↗
Grid Refinement Techniques for the 𝜸-Re𝜽𝒕 Transition Model in FUN3D
There has been an increased focus on the overall accuracy and grid convergence of Reynolds-averaged Navier-Stokes (RANS)-based transition models from the recent AIAA and NATO-AVT workshops. Even though satisfactory grid convergence could be achieved for simple two-dimensional flow configurations, it required mesh counts that are substantially larger than those used in typical applications. In this paper, we focus our efforts on understanding how the grid resolution and topology influences the accuracy and convergence of results by studying the Schubauer and Skramstad flat-plate configuration and the NLF-0416 airfoil at an angle of attack equal to five degrees using FUN3D, a second-order finite-volume code. By focusing on these cases, we can analyze both natural and separation-induced transition scenarios. Multiple grid refinement strategies are investigated. First, we determine the relative effectiveness of zonal streamwise refinement in the transition region within structured grids as an alternative to the costly option of globally uniform refinement of a baseline grid. We also complement this zonal technique by globally varying the wall-normal resolution keeping the streamwise resolution fixed. The zonal streamwise refinement can accurately model natural transition in a flat-plate boundary layer and separation-induced transition, but struggles to accurately model natural transition in airfoil flows. A series of unstructured prismatic grids that have similar node counts and viscous wall spacings as the structured hexahedral grids are also tested, and they do not achieve grid convergence until an extremely fine resolution. Last, we employ adjoint-based unstructured grid adaptation in FUN3D to natural and separation-induced transition on the NLF-0416 airfoil. The adjoint-based refinement process converges to the same solution as the baseline family of structured grids, but leads to smaller errors on coarser grids.
Grid Adaptation↗
Modeling Boundary-Layer Transition in Subsonic Flow over a Swept Wing
Predicting the onset of boundary-layer transition is often more accurate using physics-based models that directly compute disturbance growth rather than phenomenological models often implemented into industrial CFD codes. The aim of this ongoing study is to calibrate linear, physics-based computations of transition in subsonic flows over swept wings against a large set of experimental data. Advancing the calibration of linear models of transition contributes to the CFD-Vision-2030 goal of automated boundary-layer transition prediction. This progress report uses the dual N-factor method to model transition over the swept NACA 64-2-015A wing. The flow conditions match selected test conditions from an extensive experimental dataset acquired from the NASA Ames 12-ft Pressure Tunnel. The OVERFLOW 2.4b flow solver is used to obtain laminar basic states based on an infinite-span assumption. Stability analyses are performed on 365 distinct configurations with linear stability theory (LST) and parabolized stability equations (PSE) from the Langley Stability and Transition Analysis Codes (LASTRAC), modeling the growth of Tollmien-Schlichting (TS) and stationary crossflow (SCF) disturbances. From a total of 67 data points for unswept, i.e., TS-dominant configurations, the critical N-factor based on PSE is found to be N_TS = 9. The SCF critical N-factor is found to be near 8 for the highly swept, SCF-dominant configurations. Dual N-factor curves for both LST and PSE computations demonstrate a high level of interaction between TS and SCF. It may be worthwhile to investigate an alternate metric to visualize maximal SCF amplification upstream of the transition location to account for the growth of SCF modes near the leading edge, which is not considered in the conventional applications of the dual N-factor criterion.
boundary-layer transition↗
Separation and Transition on the ROTEX-T Cone-Flare
As part of NATO STO AVT-346 “Predicting Hypersonic Boundary-Layer Transition on Complex Geometries,” coordinated experimental and computational studies were conducted on the ROTEX-T, a cone-flare geometry used in a successful flight-test experiment. At the as-flown conditions, a separation bubble existed at the compression corner. Separation, reattachment, and the multifaceted linear instability paths leading this bubble to transition to turbulence are challenging to predict, but have significant impact on surface pressure and heat flux. High-resolution background-oriented schlieren and infrared thermography measurements were made in the AFOSR–Notre Dame Large Mach-6 Quiet Tunnel at freestream unit Reynolds numbers from 5.8×10(exp 6) to 12.2×10(exp 6) m-1and nominally zero angle of attack. High-speed self-aligned focusing schlieren, infrared thermography, and focused laser differential interferometry measurements were made in the AFRL Mach-6 Ludwieg Tube from 2.2×10(exp 6) to 24.7×10(exp 6) m-1and nominally zero degrees angle of attack. The surface heat-flux and Stanton number distributions were computed. Separation and reattachment locations, as well as the flow state at each, were determined from the combination of surface and off-wall measurements. The convective and global boundary-layer instabilities of the axisymmetric laminar flow at the experimental conditions were investigated computationally. Amplification of Mack’s first and second modes were observed to have logarithmic amplification factors between 5 to 7.5 at the separation location, depending on conditions. The flow was found to be globally unstable to stationary three-dimensional disturbances concentrated in the reattachment region. Previous analysis of the ROTEX-T flight data had not assessed reattachment location or the flow state upon reattachment. Thanks to the insights gained from the coordinated, on- and off-wall ground-test measurements, these evaluations have now been made. The separation location indicated by laminar simulations is consistently numerically predicted to be upstream of the experimentally observed location for a transitional separation bubble. The cause of this difference is understood to lie in the steady-state and axisymmetric assumptions made by both solvers employed to compute the basic states analyzed, as flow topology considerations assert that unsteady two-dimensional or axisymmetric separation bubbles are structurally unstable and will become three-dimensional. Computed laminar heating rates prior to separation agreed well with experiment; transitional heating rates after reattachment were between laminar and turbulent computations.
CFD Validation↗
Separation and Transition on a CCF: Experimental Campaigns
Several experimental campaigns have been conducted across a number of Mach-6 facilities on the hypersonic flow around a CCF (CCF) geometry with a 5° half-angle cone and a 12° half-angle flare. These experiments were conducted as part of the NATO STO Research Task Group AVT-346, which is focused on predicting hypersonic boundary-layer transition on complex geometries. Two conventional wind tunnels (AFRL M6LT and ONERA R2Ch) and one quiet tunnel (Purdue BAM6QT) were used to test the same CCF geometry and to study the instabilities in both the boundary layer in the attached parts of the flow and the shear layer above the axisymmetric separation bubble near the cylinder-cone junction. Two nosetip radii (one nominally sharp and one blunt with a 5 mm radius) were tested. For the sharp nosetip case, there was a great deal of agreement between the measurements of both second-mode and shear-layer instabilities across the two conventional facilities. However, the measured spectra for the blunt nosetip case showed more significant differences between the two tunnels, potentially due to an alternate dominant instability mechanism coupled with the variations in the freestream noise spectra. The quiet facility resulted in a flow that remained laminar to much higher freestream unit Reynolds numbers, as well as in instability measurements that had more distinct spectral peaks for the sharp tip case and broadband rises for the blunt one. The instability mechanisms at play in the sharp quiet case were found to be the same as those in the conventional facilities.
Boundary Layer Transition↗
Hypersonic Boundary-Layer Transition over a Blunt Circular Cone in a Mach 8 Digital Wind Tunnel
To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, three-dimensional direct numerical simulations (DNS) were conducted to investigate the interaction of a tunnel-like acoustic disturbance field with a Mach 8 laminar boundary layer over a 7 deg half-angle blunt cone with a nose radius of 𝑅 𝑛 = 5.2 mm. The flow conditions and the axisymmetric cone geometry matched those measured in the Sandia Hypersonic Wind Tunnel at Mach 8 (Sandia HWT-8), and the incident freestream acoustic disturbances in the DNS were generated by acoustic radiation from the nozzle-wall turbulent boundary layer in the same tunnel. Consistent with the predictions of the nonmodal instability analysis, the study shows the emergence of a second spectral peak for moderately oblique waves between 150 kHz and 200 kHz in the aft part of the cone. The spectral energy of temperature and streamwise velocity fluctuations is initially concentrated within the entropy layer. These fluctiations are eventually “swallowed” by the boundary layer as the entropy-layer thickness 𝛿 𝑆 becomes similar to the boundary-layer thickness 𝛿h𝑡 at downstream locations. The numerical schlieren contours show the same inclined structures commonly observed in blunt cone experiments, and these structures are found to be three dimensional and advect with a constant speed approximately equal to the mean flow speed near the edge of the boundary layer.
Hypersonic Flow↗
Ground Tests on the BOLT Geometry at Mach 7: Cross-Facility Comparison and Stability Analysis
Transition measurements have been obtained through two experimental campaigns conducted independently by the German Aerospace Center and the French Aerospace Lab with the French Alternative Energies and Atomic Energy Commission on subscale models of the BOLT-1 flight experiment geometry. This paper details a cross-facility comparison of measurements obtained at Mach 7, as well as subsequent computational analysis. Infrared (IR) thermography measurements obtained by both campaigns have facilitated a global comparison of the transition front across facilities at analagous conditions, which are found to be in good agreement. High-frequency surface pressure fluctuation data demonstrate significant amplification of instabilities with Mack-mode characteristics in the outboard regions of the acreage. These measurements are compared to stability analyses of varying fidelity. The computational methods employed to characterize the boundary-layer transition phenomena include the traditional line-marching implementation of the Parabolized Stability Equations (PSE), 2D eigenvalue analysis coupled with PSE, and Adaptive Mesh Refinement Wavepacket Tracking. While the 2D eigenvalue analysis is found to predict instabilities, which correlate in terms of frequency and acreage location to the experimental measurements, the predicted amplification for these instabilities is lower than would typically be expected for transition. Line-marching results for traveling crossflow produce the best match to the experimental transition front, with a consistent transition N factor of approximately 3-3.5.
Boundary Layer Transition↗
Separation and Transition on a Cone-Cylinder-Flare: Computational Investigations
Base flow computation and stability analysis were conducted for hypersonic flow over a cone-cylinder-flare (CCF) geometry for conditions that correspond to the experimental runs carried out in three wind tunnels. Owing to the presence of an attached boundary layer, a separation bubble induced by a shock-boundary layer interaction, and a reattachment region, the chosen flow configuration is physically rich. The complexity of this flowfield encompasses a combination of convective instabilities developing on the cone, global instabilities in the separation bubble, and shear-layer modes and streaks in the reattachment region. Thus, the selected CCF configuration provides the opportunity for a comprehensive comparison of the currently available methodologies for analyzing boundary layer instabilities. Various tools are used for the analysis, including global stability codes as well as convective instability analyses based on a local theory, a weakly non-parallel analysis, and tools that are applicable to strongly non-parallel flows. The paper presents a comparison of the convective instability characteristics based on different methodologies, such as linear stability theory (LST), the harmonic form of linearized Navier-Stokes equations (HLNSE), and resolvent analysis. The CCF configuration provided an effective framework for conducting a detailed cross-validation of this type, which had not yet been addressed in existing literature. This document is accompanied by a companion paper that is focused on the experimental aspects of the CCF configuration. Both papers are being presented in a dedicated session that highlights the research activities of the NATO STO Research Task Group AVT-346.
Boundary Layer Transition↗
Grid Refinement Techniques for the 𝜸-Re 𝜽 𝒕 Transition Model in FUN3D
There has been an increased focus on the overall accuracy and grid convergence of Reynolds-averaged Navier-Stokes (RANS)-based transition models from the recent AIAA and NATO-AVT workshops. Even though satisfactory grid convergence could be achieved for simple two-dimensional flow configurations, it required mesh counts that are substantially larger than those used in typical applications. In this paper, we focus our efforts on understanding how the grid resolution and topology influences the accuracy and convergence of results by studying the Schubauer and Skramstad flat-plate configuration and the NLF-0416 airfoil at an angle of attack equal to five degrees using FUN3D, a second-order finite-volume code. By focusing on these cases, we can analyze both natural and separation-induced transition scenarios. Multiple grid refinement strategies are investigated. First, we determine the relative effectiveness of zonal streamwise refinement in the transition region within structured grids as an alternative to the costly option of globally uniform refinement of a baseline grid. We also complement this zonal technique by globally varying the wall-normal resolution keeping the streamwise resolution fixed. The zonal streamwise refinement can accurately model natural transition in a flat-plate boundary layer and separation-induced transition, but struggles to accurately model natural transition in airfoil flows. A series of unstructured prismatic grids that have similar node counts and viscous wall spacings as the structured hexahedral grids are also tested, and they do not achieve grid convergence until an extremely fine resolution. Last, we employ adjoint-based unstructured grid adaptation in FUN3D to natural and separation-induced transition on the NLF-0416 airfoil. The adjoint-based refinement process converges to the same solution as the baseline family of structured grids, but leads to smaller errors on coarser grids.
Grid Adaptation↗
Modeling Boundary-Layer Transition in Subsonic Flow over a Swept Wing
Predicting the onset of boundary-layer transition is often more accurate using physics-based models that directly compute disturbance growth rather than phenomenological models often implemented into industrial CFD codes. The aim of this ongoing study is to calibrate linear, physics-based computations of transition in subsonic flows over swept wings against a large set of experimental data. Advancing the calibration of linear models of transition contributes to the CFD-Vision-2030 goal of automated boundary-layer transition prediction. This progress report uses the dual N-factor method to model transition over the swept NACA 64-2-015A wing. The flow conditions match selected test conditions from an extensive experimental dataset acquired from the NASA Ames 12-ft Pressure Tunnel. The OVERFLOW 2.4b flow solver is used to obtain laminar basic states based on an infinite-span assumption. Stability analyses are performed on 365 distinct configurations with linear stability theory (LST) and parabolized stability equations (PSE) from the Langley Stability and Transition Analysis Codes (LASTRAC), modeling the growth of Tollmien-Schlichting (TS) and stationary crossflow (SCF) disturbances. From a total of 67 data points for unswept, i.e., TS-dominant configurations, the critical N-factor based on PSE is found to be N_TS = 9. The SCF critical N-factor is found to be near 8 for the highly swept, SCF-dominant configurations. Dual N-factor curves for both LST and PSE computations demonstrate a high level of interaction between TS and SCF. It may be worthwhile to investigate an alternate metric to visualize maximal SCF amplification upstream of the transition location to account for the growth of SCF modes near the leading edge, which is not considered in the conventional applications of the dual N-factor criterion.
computational modeling↗
Hypersonic Boundary-Layer Transition over a Blunt Circular Cone in a Mach 8 Digital Wind Tunnel
To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, three-dimensional direct numerical simulations (DNS) were conducted to investigate the interaction of a tunnel-like acoustic disturbance field with a Mach 8 laminar boundary layer over a 7 deg half-angle blunt cone with a nose radius of 𝑅 𝑛 = 5.2 mm. The flow conditions and the axisymmetric cone geometry matched those measured in the Sandia Hypersonic Wind Tunnel at Mach 8 (Sandia HWT-8), and the incident freestream acoustic disturbances in the DNS were generated by acoustic radiation from the nozzle-wall turbulent boundary layer in the same tunnel. Consistent with the predictions of the nonmodal instability analysis, the study shows the emergence of a second spectral peak for moderately oblique waves between 150 kHz and 200 kHz in the aft part of the cone. The spectral energy of temperature and streamwise velocity fluctuations is initially concentrated within the entropy layer. These fluctiations are eventually “swallowed” by the boundary layer as the entropy-layer thickness 𝛿 𝑆 becomes similar to the boundary-layer thickness 𝛿h𝑡 at downstream locations. The numerical schlieren contours show the same inclined structures commonly observed in blunt cone experiments, and these structures are found to be three dimensional and advect with a constant speed approximately equal to the mean flow speed near the edge of the boundary layer.
Hypersonic Flow↗
Aerodynamic Design Optimization for Natural Laminar Flow Airfoils
Natural laminar flow technology is a passive laminar flow control (LFC) strategy that seeks to delay the onset of boundary-layer transition (BLT) through shape optimization to reduce the drag of the aerodynamic vehicle. Adjoint-based design optimization for LFC is proposed in an integrated multidisciplinary framework, which includes the computational fluid dynamics (CFD), geometry and grid deformation, and linear stability analysis (LSA) for transition prediction. In particular, the BLT location is predicted using the dual N-factor method that is based on a linear stability theory (LST) eigenvalue problem. The dual N-factor criterion accounts for the amplification of planar Tollmien-Schlichting (TS) and stationary crossflow (CF) boundary-layer instabilities to predict the transition location in three-dimensional boundary-layer flows. The adjoint-based shape optimization procedure is based on an iteratively coupled CFD and LSA methodology to converge the transition location and flow solutions, as well as to calculate the sensitivities of the aerodynamic metrics of interest with respect to the flow and shape design parameters. The RAE 2822 airfoil at 0 and 30 degrees yaw angles, an angle of attack of 0.72 degrees, and subsonic conditions (M∞ = 0.19, Rec = 5.6 × 106 ) are used as baseline configurations for design optimization. The angle of attack and the vertical displacement of free-form-deformation control points are used as design variables to reduce the drag coefficient while reaching a specified lift coefficient. The optimized unswept airfoil designs achieve a 30% drag reduction accompanied by a downstream shift of the transition locations over both suction and pressure sides of the airfoil. The initial design iterations for the swept case also show a favorable trend in the drag reduction with transition delay over both sides.
Transition↗
The harmonic linearized Navier–Stokes equations for transition prediction in three-dimensional flows
Explore the source record for details and available documents.
Pedro Paredes↗
Zonal and adaptive mesh refinement for RANS-based transition models in FUN3D
Explore the source record for details and available documents.
Nathaniel Hildebrand↗