Application of RANS-based Langtry-Menter Transition Model in OVERFLOW 2.3b to a 6:1 Prolate Spheroid
Explore the source record for details and available documents.
CFD Modeling↗
Engineering topics
Meelan M Choudhari
Publications and source records attributed to Meelan M Choudhari.
Aeroacoustic Simulations of the High-Lift Common Research Model and Validation with Experiment
Aerodynamic and acoustic data from unsteady numerical simulations of the High-Lift Common Research Model are compared with experimental measurements from an open-jet wind tunnel test of a 10\-scale model. Time-averaged surface pressure data is in good agreement for all of the configurations evaluated during the test. Unsteady surface pressure spectra from the leading-edge slat and the nose of the main element are presented, demonstrating that the highest levels are associated with the wakes of slat brackets. At these locations the simulations and experiment match reasonably well at low enough frequencies (f < 10 kHz based on the model scale) before excessive numerical dissipation leads to a rapid roll off in the predicted frequency spectra. The agreement is inconsistent at other locations on the slat and main-element surfaces, where the pressure fluctuation levels are lower. Noise predictions using a synthetic microphone array are compared with equivalent experimental results, and the changes in levels predicted between three configurations are similar to those measured. However, some details of the noise predictions are inconsistent with the experiment, and the use of solid- instead of porous-surface data for the acoustic processing of the simulations likely contributes to the discrepancies. All of the simulations were completed before the test, and, despite some shortcomings, they provided valuable insight into the aeroacoustic performance of the model that greatly aided test planning and execution. Hence, the post-test comparisons presented here allow an assessment of the predictive methodology in the context of a realistic high-lift configuration.
High-lift↗
Effects of Porous Gap Fillers on 30P30N Leading-Edge SlatNoise. Part I: Surface Pressure and Acoustics
The leading edge slat of a high-lift system is one of the main contributors to airframe noise during approach. In a previous experimental study, we assessed the performance of an impermeable slat gap filler as a passive flow control device to reduce the slat noise associated with a two-dimensional, three-element high-lift airfoil. The present paper, the first of two parts,represents a follow-on investigation to assess the relative efficacy of permeable gap fillers thatallow successively higher amounts of flow to pass through the gap. To evaluate the influence of this passive flow control device on the acoustics generated by the unsteady flow near the slat, experiments are conducted in an anechoic wind tunnel by mounting a gap filler to the slat element of the two-dimensional 30P30N high-lift configuration. Measurements are performed at a single geometric angle of attack (𝛼𝑘=8◦) and three different flow speeds that correspond to Reynolds numbers of𝑅𝑒𝑐=1.2𝑒6,1.5𝑒6, and1.71𝑒6, respectively. Steady surface pressure measurements are used to gauge the influence of the permeable gap filler treatments on the overall lift. The effect of each treatment on the radiated noise is analyzed via acoustic array measurements, followed by delay-and-sum beamforming to locate the slat noise sources and to provide the integrated acoustic spectra. The porous gap fillers are found to eliminate the narrowband peaks in the acoustic spectra and, also, to yield a 10 dB reduction in the broadband noise in comparison with the baseline case with no gap filler. The porous gap filler with thelowest permeability acts similar to the impermeable gap filler examined previously. However, the aerodynamics and noise reduction both degrade with increasing permeability. An accompanying abstract describes the particle image velocimetry measurements of the flowfield within the slat cove and on either side of the permeable gap f
slat noise↗
Effects of Porous Gap Fillers on Leading-Edge Slat Noise of 30P30N. Part II: PIV Measurements
The leading edge slat of a high-lift system is one of the main contributors to the airframe noise during approach conditions. This paper, the second of two parts, continues our previous studies on the slat gap filler as a passive noise control device on the two-dimensional, 30P30N multielement airfoil. Whereas the earlier study was focused on the effects of an impermeable gap filler that completely blocks the flow through the gap, this follow-on assessment is devoted to permeable slat gap fillers that allow limited amounts of flow to pass through the gap. Part I of this two-part abstract described the aerodynamic and acoustic effects of the permeable gap fillers, as inferred from both the measurements of static and unsteady surface pressures and the microphone array data for the radiated noise. To understand the physical mechanism responsible for the noise reduction documented in Part I, as well as for the accompanying aerodynamic penalty due to the porous gap fillers, Particle Image Velocimetry (PIV) is used in this paper to measure the flow details in the slat-cove region of the 30P30N model. A single angle of attack (𝛼 = 5.5◦) and a chord based Reynolds number of 1.71𝑒6 are selected as the test conditions. The PIV results show that the slat flow features are significantly altered with the presence of the porous gap filler, resulting in a more stable slat-cove shear layer and, thus, reduced velocity fluctuations with a successive decrease in the permeability. The porous gap filler with the lowest permeability acts similar to the solid gap filler. However, flow separation is observed on the upper side of the porous interface, which leads to an aerodynamic performance penalty via a reduction in lift on the main wing.
slat noise↗
Toward Verification of the γ-Reθt Transition Model in OVERFLOW and FUN3D
The results of an ongoing assessment of the transition modeling capability in NASA's OVERFLOW and FUN3D programs are presented, with a focus on the Langtry-Menter γ-Reθt transition model in combination with the shear-stress transport (SST) turbulence model. While the effect of numerics and boundary conditions on the accuracy and iterative convergence of RANS solutions for fully turbulent flows has been well documented, especially in the context of canonical flow configurations, the same cannot be said for transport-equation-based transition models coupled with RANS-based turbulence models. Given the criticality of transition modeling for new aircraft design and optimization, there has been a renewed focus on the accuracy of such transition models and their inconsistent implementation across different flow solvers as seen from the AIAA and NATO-AVT workshops. In this work, we aim to establish the verification of such models via these two well established CFD codes with different numerics. The goal is to produce high-quality data, such as grids, solutions, and other auxiliary data, that may be utilized for code verification by others in the CFD community. As a first step, the work reported here is focused on the Langtry-Menter transition model as applied to two simple 2D configurations, namely the flat plate, and the NLF-0416, respectively. The paper also highlights how the boundary conditions and baseline turbulence model can affect the solutions from the SST-based Langtry-Menter transition model. A preliminary evaluation of the automatic mesh adaption capabilities of these solvers and potential benefits for flow configurations involving a mix of laminar, transitional, and fully turbulent flows is also reported.
CFD Modeling↗
Effect of the Reynolds Number on the Freestream Disturbance Environment in a Mach 6 Nozzle
To understand the impact of unit Reynolds number on the acoustic disturbance field inside a high-speed wind tunnel, we use Direct Numerical Simulations (DNS) to model the turbulent boundary layers along the walls of a quasi-two-dimensional nozzle configuration. Intended as a stepping stone to fully three-dimensional simulations of freestream noise inside the NASA 20-Inch Mach 6 Wind Tunnel, the present simulations are based on periodic boundary conditions across the spanwise width of the computational domain that corresponds to about one third of the actual tunnel width. These simulations are performed at four different unit Reynolds numbers, ranging from 3.56e6 to 14.0e6 per meter. The predominantly hydrodynamic fluctuations inside the boundary layer are shown to be nearly unaffected by the presence of freestream forcing associated with the impinging acoustic radiation from the opposite wall. Thus, the Reynolds number trends associated with boundary-layer quantities are consistent with previously published DNS of flat-plate boundary layers at similar Mach numbers and wall temperature ratios. The unsteady disturbance environment within the nozzle core region is found to be approximately spatially homogeneous and purely acoustic in nature. The numerical results are used to make comparisons with tunnel noise measurements by Chou et al. [1,2]. Unlike previous comparisons involving static-pressure fluctuations based on the DNS and pitot-pressure fluctuations measured in the wind tunnel, direct comparisons involving the fluctuations in the same physical quantity, namely, the streamwise mass flux, have been reported for the first time. The predicted decrease in the root-mean-square fluctuations in pressure and mass flux with an increasing unit Reynolds number is in agreement with the measurements in the NASA 20-Inch Mach 6 Wind Tunnel. Additional details of the acoustic radiation field are quantified and should be useful toward a digital synthesis of the tunnel disturbance environment that would enable realistic simulations of the natural transition process.
Acoustics↗
Toward Verification of the 𝜸−𝑹𝒆_𝜽𝒕 Transition Model in OVERFLOW and FUN3D
Explore the source record for details and available documents.
boundary layer transition↗
DeepONet-Assisted Optimization of Surface Topography for Transition Delay in a Mach 4.5 Boundary Layer
We use deep learning, an ensemble variational technique (EnVar), and direct numerical simulations(DNS) to design an optimal topography for a two-dimensional roughness element that delays the on-set of laminar-turbulent transition in a Mach 4.5 flat-plate boundary layer. Deep operator networks (DeepONets), which have the known ability to learn complex nonlinear operators within dynamical systems, are used for machine learning. For the baseline configuration of a smooth flat plate, the second-mode waves at the DNS inflow cause a quick nonlinear breakdown of the high-speed boundary layer within the computational domain. Results reported in the present study validate the ability of DeepONets to model the transition delay via a given topography of the roughness element. The computing cost to optimize the rough-ness element for minimal skin-friction drag is substantially lowered by the DeepONets-based reduced-order model. In comparison to the baseline method of EnVar optimization based on DNS alone, the DeepONets-based EnVar optimizer is able to delay transition past the outflow boundary of the computational domain while utilizing almost 5–6 times fewer DNS.
Machine Learning↗
Recent Progress on RANS-Based Transition Model Verification
The current efforts to assess and improve the Reynolds-averaged Navier-Stokes (RANS)-coupled transition models in the NASA FUN3D and OVERFLOW codes are summarized in this study. The first AIAA Transition Modeling Workshop and the NATO AVT-313 Transition Workshop both emphasized the need for code verification for transport equations based transition models as a top priority. We discuss the methods used for the model verification, the resulting grid families, the flow solutions, and other supporting information collected with at least two established NASA flow solvers, namely, FUN3D and OVERFLOW. These results, which will be uploaded onto the NASA Turbulence Modeling Resource, should assist other members of the computational fluid dynamics (CFD) community in verifying their own implementations of various transition models, such as the Langtry-Menter (LM2009) model, the one-equation γ model, and Coder’s amplification factor transport (AFT) model. Grid convergence is assessed using both global and local flow metrics of interest such as lift and drag as well as local skin-friction coefficients. We also explore the anisotropic unstructured metric-based adaptive mesh refinement library known as refine with the NASA FUN3D solver to determine if this capability can achieve the same accuracy as handcrafted structured grids with a significantly smaller node count and to learn the characteristics of the resulting grid distribution, especially in the vicinity of the transition zone.
RANS↗
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↗
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↗
Enhancements to Linear Stability-Based, CFD-integrated Transition Prediction for High-Speed Flows
Combining linear stability calculations with computational fluid dynamics (CFD) simulations has great potential for the automated modeling of high-speed flows, especially when adequate information about the configuration and the disturbance environment is available. However, a significant impediment to the applicability of this technique is the lack of an efficient method to calculate the crucial amplification ratio corresponding to the onset of transition in hypersonic flows. This ratio, also known as the "transition N-factor," is dependent upon the freestream disturbance environment as well as the surface properties of the test article. In response to the need for an engineering solution to predict the transition N-factor within conventional hypersonic wind tunnels, this paper presents a data-driven correlation that expands the existing correlations from straight circular cones with a narrow range of half angles to a broader array of axisymmetric configurations. Furthermore, when tested against a chosen dataset that was not used in its calibration, the suggested correlation shows good predictive accuracy with an RMS error of only 6.9%. Although similar accuracy may also be achieved via existing correlations based on similar datasets, predictions based on the proposed correlation have the advantage of not requiring an extensive amount of configuration-specific data. Practical applications often have access to the input parameters needed for this correlation, such as the freestream disturbance intensity, Mach number, and body-based slenderness Reynolds number. Additionally, this correlation outperforms the traditional assumption of a constant N-factor, particularly for configurations with blunted nose geometries. The development of this correlation is grounded in an extensive dataset encompassing conical models with body half-angles varying between 5 degrees and 16 degrees, Mach numbers ranging from 5 to 14, and nosetip-based Reynolds numbers approaching the transition reversal limit for blunt-nosed cones.
CFD↗
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↗
Computational Study of Noise Characteristics for 30P30N Leading-edge Slat and Krueger Flap
In this study, numerical simulations are performed to investigate the flowfield, the noise sources, and their detailed characteristics around high-lift configurations featuring either a conventional leading-edge slat or a Krueger flap. While the flowfields and potential noise sources are similar for the two types of configurations, the noise directivity patterns exhibit some differences. The influence of the cavity where the Krueger flap is stowed on the lower surface is also examined via numerical simulations with and without the cavity. Results show that the cavity does not have a significant influence on the flowfield around the Krueger flap and the gap region. However, when the cavity is present, additional noise sources are observed around its rear part, although they are not as significant as the noise sources on the lower side of the Krueger flap. These additional noise sources contribute to the directivity of high noise levels, especially in the upstream direction.
Aeroacoustics↗
Transition Prediction for the 30P30N Airfoil and the CRM-HL using FUN3D
The ability to accurately model transition over high-lift configurations can lead to improved agreement between stationary Reynolds-Averaged Navier-Stokes (RANS) computations and experimental measurements. In this work, selected transport-equation-based RANS models for transition prediction in the NASA FUN3D flow solver are used to simulate flow over the 30P30N three-element airfoil and the High-Lift Common Research Model (CRM-HL). The primary RANS-based transition model is the Langtry-Menter (LM) 𝜸-Re𝜽𝒕 model coupled with the Menter shear-stress transport model. Additionally, the existing capability for transition prediction based on linear stability correlations is used to allow the modeling of transition over multiple elements of the 30P30N configuration. This approach pairs an algebraic transition model with the Spalart-Allmaras (SA) turbulence model that relies on an imposed transition location in FUN3D and the parabolized stability equations as implemented in the LASTRAC code. Predictions based on these different transition prediction approaches and models are compared with each other, experimental data, and previous simulations by evaluating the lift and drag coefficients, the surface-pressure/skin-friction distributions, and the locations of transition onset. Results are presented for multiple angles of attack for the 30P30N airfoil and CRM-HL. For the CRM-HL, a grid resolution study is performed for the LM transition model and the SA turbulence model.
High Lift↗