Search NASA⌕ Search

SEARCH · Search NASA

Results for “PRESSURE COEFFICIENT”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

Comparison of Ares I-X Wind-Tunnel Derived Buffet Environment with Flight Data

The Ares I-X Flight Test Vehicle (FTV), launched in October 2009, carried with it over 243 buffet verification pressure sensors and was one of the most heavily instrumented launch vehicle flight tests. This flight test represented a unique opportunity for NASA and its partners to compare the wind-tunnel derived buffet environment with that measured during the flight of Ares I-X. It is necessary to define the launch vehicle buffet loads to ensure that structural components and vehicle subsystems possess adequate strength, stress, and fatigue margins when the vehicle structural dynamic response to buffet forcing functions are considered. Ares I-X buffet forcing functions were obtained via wind-tunnel testing of a rigid buffet model (RBM) instrumented with hundreds of unsteady pressure transducers designed to measure the buffet environment across the desired frequency range. This paper discusses the comparison of RBM and FTV buffet environments, including fluctuating pressure coefficient and normalized sectional buffet forcing function root-mean-square magnitudes, frequency content of power-spectral density functions, and force magnitudes of an alternating flow phenomena. Comparison of wind-tunnel model and flight test vehicle buffet environments show very good agreement with root-mean-square magnitudes of buffet forcing functions at the majority of vehicle stations. Spectra proved a challenge to compare because of different wind-tunnel and flight test conditions and data acquisition rates. However, meaningful and promising comparisons of buffet spectra are presented. Lastly, the buffet loads resulting from the transition of subsonic separated flow to supersonic attached flow were significantly over-predicted by wind-tunnel results.

Piatak, David J.↗

Preliminary Computational Analysis of the (HIRENASD) Configuration in Preparation for the Aeroelastic Prediction Workshop

This paper presents preliminary computational aeroelastic analysis results generated in preparation for the first Aeroelastic Prediction Workshop (AePW). These results were produced using FUN3D software developed at NASA Langley and are compared against the experimental data generated during the HIgh REynolds Number Aero- Structural Dynamics (HIRENASD) Project. The HIRENASD wind-tunnel model was tested in the European Transonic Windtunnel in 2006 by Aachen University0s Department of Mechanics with funding from the German Research Foundation. The computational effort discussed here was performed (1) to obtain a preliminary assessment of the ability of the FUN3D code to accurately compute physical quantities experimentally measured on the HIRENASD model and (2) to translate the lessons learned from the FUN3D analysis of HIRENASD into a set of initial guidelines for the first AePW, which includes test cases for the HIRENASD model and its experimental data set. This paper compares the computational and experimental results obtained at Mach 0.8 for a Reynolds number of 7 million based on chord, corresponding to the HIRENASD test conditions No. 132 and No. 159. Aerodynamic loads and static aeroelastic displacements are compared at two levels of the grid resolution. Harmonic perturbation numerical results are compared with the experimental data using the magnitude and phase relationship between pressure coefficients and displacement. A dynamic aeroelastic numerical calculation is presented at one wind-tunnel condition in the form of the time history of the generalized displacements. Additional FUN3D validation results are also presented for the AGARD 445.6 wing data set. This wing was tested in the Transonic Dynamics Tunnel and is commonly used in the preliminary benchmarking of computational aeroelastic software.

Chwalowski, Pawel↗

Development, Verification and Use of Gust Modeling in the NASA Computational Fluid Dynamics Code FUN3D

This paper presents the implementation of gust modeling capability in the CFD code FUN3D. The gust capability is verified by computing the response of an airfoil to a sharp edged gust. This result is compared with the theoretical result. The present simulations will be compared with other CFD gust simulations. This paper also serves as a users manual for FUN3D gust analyses using a variety of gust profiles. Finally, the development of an Auto-Regressive Moving-Average (ARMA) reduced order gust model using a gust with a Gaussian profile in the FUN3D code is presented. ARMA simulated results of a sequence of one-minus-cosine gusts is shown to compare well with the same gust profile computed with FUN3D. Proper Orthogonal Decomposition (POD) is combined with the ARMA modeling technique to predict the time varying pressure coefficient increment distribution due to a novel gust profile. The aeroelastic response of a pitch/plunge airfoil to a gust environment is computed with a reduced order model, and compared with a direct simulation of the system in the FUN3D code. The two results are found to agree very well.

Bartels, Robert E.↗

CFD Simulations for the Effect of Unsteady Wakes on the Boundary Layer of a Highly Loaded Low-Pressure Turbine Airfoil (L1A)

This report is the third part of a three-part final report of research performed under an NRA cooperative Agreement contract. The first part was published as NASA/CR-2012-217415. The second part was published as NASA/CR-2012-217416. The study of the very high lift low-pressure turbine airfoil L1A in the presence of unsteady wakes was performed computationally and compared against experimental results. The experiments were conducted in a low speed wind tunnel under high (4.9%) and then low (0.6%) freestream turbulence intensity for Reynolds number equal to 25,000 and 50,000. The experimental and computational data have shown that in cases without wakes, the boundary layer separated without reattachment. The CFD was done with LES and URANS utilizing the finite-volume code ANSYS Fluent (ANSYS, Inc.) under the same freestream turbulence and Reynolds number conditions as the experiment but only at a rod to blade spacing of 1. With wakes, separation was largely suppressed, particularly if the wake passing frequency was sufficiently high. This was validated in the 3D CFD efforts by comparing the experimental results for the pressure coefficients and velocity profiles, which were reasonable for all cases examined. The 2D CFD efforts failed to capture the three dimensionality effects of the wake and thus were less consistent with the experimental data. The effect of the freestream turbulence intensity levels also showed a little more consistency with the experimental data at higher intensities when compared with the low intensity cases. Additional cases with higher wake passing frequencies which were not run experimentally were simulated. The results showed that an initial 25% increase from the experimental wake passing greatly reduced the size of the separation bubble, nearly completely suppressing it.

Vinci, Samuel, J.↗

Supersonic Retropropulsion CFD Validation with Ames Unitary Plan Wind Tunnel Test Data

A validation study of Computational Fluid Dynamics (CFD) for Supersonic Retropropulsion (SRP) was conducted using three Navier-Stokes flow solvers (DPLR, FUN3D, and OVERFLOW). The study compared results from the CFD codes to each other and also to wind tunnel test data obtained in the NASA Ames Research Center 90 70 Unitary PlanWind Tunnel. Comparisons include surface pressure coefficient as well as unsteady plume effects, and cover a range of Mach numbers, levels of thrust, and angles of orientation. The comparisons show promising capability of CFD to simulate SRP, and best agreement with the tunnel data exists for the steadier cases of the 1-nozzle and high thrust 3-nozzle configurations.

Schauerhamer, Daniel G.↗

Surface-Pressure and Flow-Visualization Data at Mach Number of 1.60 for Three 65 deg Delta Wings Varying in Leading-Edge Radius and Camber

An experimental investigation of the effect of leading-edge radius, camber, Reynolds number, and boundary-layer state on the incipient separation of a delta wing at supersonic speeds was conducted at the Langley Unitary Plan Wind Tunnel at Mach number of 1.60 over a free-stream Reynolds number range of 1 x 106 to 5 x 106 ft-1. The three delta wing models examined had a 65 deg swept leading edge and varied in cross-sectional shape: a sharp wedge, a 20:1 ellipse, and a 20:1 ellipse with a -9.750 circular camber imposed across the span. The wings were tested with and without transition grit applied. Surface-pressure coefficient data and flow-visualization data are electronically stored on the CD-ROM. The data indicated that by rounding the wing leading edge or cambering the wing in the spanwise direction, the onset of leading-edge separation on a delta wing can be raised to a higher angle of attack than that observed on a sharp-edged delta wing. The data also showed that the onset of leading-edge separation can be raised to a higher angle of attack by forcing boundary-layer transition to occur closer to the wing leading edge by the application of grit or the increase in free-stream Reynolds number.

McMIllin, S. Naomi↗

Computational Aerodynamic Analysis of Three-Dimensional Ice Shapes on a NACA 23012 Airfoil

The present study identifies a process for performing computational fluid dynamic calculations of the flow over full three-dimensional (3D) representations of complex ice shapes deposited on aircraft surfaces. Rime and glaze icing geometries formed on a NACA23012 airfoil were obtained during testing in the NASA Glenn Research Center's Icing Research Tunnel (IRT). The ice shape geometries were scanned as a cloud of data points using a 3D laser scanner. The data point clouds were meshed using Geomagic software to create highly accurate models of the ice surface. The surface data was imported into Pointwise grid generation software to create the CFD surface and volume grids. It was determined that generating grids in Pointwise for complex 3D icing geometries was possible using various techniques that depended on the ice shape. Computations of the flow fields over these ice shapes were performed using the NASA National Combustion Code (NCC). Results for a rime ice shape for angle of attack conditions ranging from 0 to 10 degrees and for freestream Mach numbers of 0.10 and 0.18 are presented. For validation of the computational results, comparisons were made to test results from rapid-prototype models of the selected ice accretion shapes, obtained from a separate study in a subsonic wind tunnel at the University of Illinois at Urbana-Champaign. The computational and experimental results were compared for values of pressure coefficient and lift. Initial results show fairly good agreement for rime ice accretion simulations across the range of conditions examined. The glaze ice results are promising but require some further examination.

Aerodynamics↗

Computational Aerodynamic Analysis of Three-Dimensional Ice Shapes on a NACA 23012 Airfoil

The present study identifies a process for performing computational fluid dynamic calculations of the flow over full three-dimensional (3D) representations of complex ice shapes deposited on aircraft surfaces. Rime and glaze icing geometries formed on a NACA23012 airfoil were obtained during testing in the NASA Glenn Research Centers Icing Research Tunnel (IRT). The ice shape geometries were scanned as a cloud of data points using a 3D laser scanner. The data point clouds were meshed using Geomagic software to create highly accurate models of the ice surface. The surface data was imported into Pointwise grid generation software to create the CFD surface and volume grids. It was determined that generating grids in Pointwise for complex 3D icing geometries was possible using various techniques that depended on the ice shape. Computations of the flow fields over these ice shapes were performed using the NASA National Combustion Code (NCC). Results for a rime ice shape for angle of attack conditions ranging from 0 to 10 degrees and for freestream Mach numbers of 0.10 and 0.18 are presented. For validation of the computational results, comparisons were made to test results from rapid-prototype models of the selected ice accretion shapes, obtained from a separate study in a subsonic wind tunnel at the University of Illinois at Urbana-Champaign. The computational and experimental results were compared for values of pressure coefficient and lift. Initial results show fairly good agreement for rime ice accretion simulations across the range of conditions examined. The glaze ice results are promising but require some further examination.

Aerodynamics↗

Spatial Convergence of Three Dimensional Turbulent Flows

Finite-volume and finite-element schemes, both implemented within the FUN3D flow solver, are evaluated for several test cases described on the Turbulence-Modeling Resource (TMR) web site. The cases include subsonic flow over a hemisphere cylinder, subsonic flow over a swept bump configuration, and supersonic flow in a square duct. The finite- volume and finite-element schemes are both used to obtain solutions for the first two cases, whereas only the finite-volume scheme is used for the supersonic duct. For the hemisphere cylinder, finite-element solutions obtained on tetrahedral meshes are compared with finite- volume solutions on mixed-element meshes. For the swept bump, finite-volume solutions have been obtained for both hexahedral and tetrahedral meshes and are compared with finite-element solutions obtained on tetrahedral meshes. For the hemisphere cylinder and the swept bump, solutions are obtained on a series of meshes with varying grid density and comparisons are made between drag coefficients, pressure distributions, velocity profiles, and profiles of the turbulence working variable. The square duct shows small variation due to element type or the spatial accuracy of turbulence model convection. It is demonstrated that the finite-element scheme on tetrahedral meshes yields similar accuracy as the finite- volume scheme on mixed-element and hexahedral grids, and demonstrates less sensitivity to the mesh topology (biased tetrahedral grids) than the finite-volume scheme.

Park, Michael A.↗

Wall-Modeled Lattice Boltzmann and Navier-Stokes Approaches for Separated Flows

Lattice Boltzmann (LB) and hybrid Reynolds-averaged Navier-Stokes/large eddy simulation (RANS/LES) methods within the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are applied to NASA's Revolutionary Computational Aerosciences (RCA) standard test cases for separated flows. A detailed comparison between the performance and accuracy of the two emerging numerical methodologies for turbulence resolving simulations, i.e. the LB and hybrid RANS/LES methods will be presented. This contribution addresses the RCA technical challenge to identify and down-select critical turbulence, transition, and numerical method technologies for 40% reduction in predictive error for standard turbulence separated flow test cases. Results for the 2D NASA wall-mounted hump and the axisymmetric transonic bump including time-averaged pressure coefficient, skin friction, and velocity pro les, as well as resolved and modeled Reynolds stresses for both numerical approaches will be presented and differences between LB and hybrid RANS/LES will be discussed.

Lattice Boltzmann↗

Measurement and Analysis of Terminal Shock Oscillation and Buffet Forcing Functions on a Launch Vehicle Payload Fairing

The buffet loads on a launch vehicle payload shroud can be impacted by the unsteadiness associated with a terminal shock at high subsonic speeds. At these conditions, flow accelerates to supersonic speeds on the nose of the payload fairing and is terminated by a normal shock on the cylindrical section downstream of the nose cone/cylinder shoulder. The location of the terminal shock and associated separated boundary layer is affected by the freestream Mach number, Reynolds number, and the pitch/yaw of the launch vehicle. Furthermore, even when the freestream conditions and vehicle attitude are constant, this terminal shock oscillates on the surface of the vehicle. The time-varying surface pressure associated with the terminal shock results in unsteady aerodynamic loads that may interact with vehicle structural dynamic modes and the guidance and control of the vehicle. Buffet testing of a 3-percent scale rigid buffet model of a launch vehicle cargo configuration with a tangent-ogive payload shroud was conducted in 2012 and in 2016. Initial buffet forcing functions (BFFs) utilized a coarse pressure sensor distribution on the vehicle surface in which a single longitudinal station with eight sensors observed the terminal shock environment at Mach 0.90. An examination of these circumferential pressures reveal large impulse-like pressure fluctuations and an asymmetry in pressure when the vehicle is at a nonzeroangle of attack that result in high BFFs. Revisions to the shock integration region were made based on computational fluid dynamics and shadowgraph video of shock motion to better represent the BFFs and reduce the high loads resulting from this environment. To more clearly understand this terminal shock environment, a second wind tunnel test was conducted with a dense distribution of 256 sensors at the terminal shock location. These sensor arrays presents a unique opportunity to observe the unsteady terminal shock environment and to characterize the impact of various integration schemes on the BFFs. This paper presents a summary of the development of BFFs for this terminal shock and a detailed analyses of shock region pressure coefficients, coherence, BFFs, shock location time histories, and power spectral density to help guide development of BFFs for other launch vehicle test and analysis programs.

Piatak, David J.↗

Influence of Airfoil Angle of Attack on Ice Accretion Roughness

The influence of airfoil angle attack on roughness evolution and spatial variation was investigated in the Icing Research Tunnel (IRT) at the NASA Glenn Research Center. Two airfoil models were used for the study: a 53.34-cm (21-in.) NACA 0012 model and a 152.4-cm (60-in.) HAARP-II model. For the NACA 0012, the angle of attack was varied from 0° to 3°, and the ice accretion roughness was characterized while keeping the other icing parameters such as freestream velocity, freestream stagnation temperature, accumulation parameter, and median volumetric diameter (MVD) constant. For the HAARP-II, the angle of attack was varied from -1.9°, which is the non-lifting angle of attack for the airfoil, to 3°. A series of accumulation time progression cases for the HAARP-II at 2.3° angle of attack was also performed in Appendix C and Appendix O conditions and compared to the measurements presented by McClain et al. (2018) for the non-lifting angle of attack condition. The results demonstrate the strong influence of the changing local static pressure along the surface influencing the roughness characteristics. A correlation approach for the maximum equivalent sand-grain roughness along the surface of an airfoil is presented and considers the scaled accumulation time, the stagnation point freezing fraction, and the surface pressure relative to the static and stagnation pressures of the flow. Finally, an attempt to relate the maximum roughness quantities to local pressure coefficients and local icing quantities is presented and discussed.

Icing↗

Comparison of Computational Predictions of the Mach 0.80 Transonic Truss-Braced Wing Configuration with Experimental Data

The NASA Langley and Ames Research Centers have teamed together for comparisons of computational predictions of the Boeing Mach 0.80Transonic Truss-Braced Wing (TTBW) configuration with a high-speed experimental dataset. The Mach 0.80 TTBW vehicle is a high wing, high aspect ratio configuration, designed for a high lift-to-drag ratio. System studies have predicted significant fuel burn and emissions benefits with the TTBW technology moving toward meeting NASA Subsonic Transport Systems-Level-Metrics. A 4.5% scale Mach 0.80 design TTBW model was recently tested at the NASA Ames Research Center 11-by 11-Foot Transonic Wind Tunnel(11-Ft TWT),providing a valuable dataset to validate computational tools and investigate best practices as risk reduction efforts continue for the development of the advanced TTBW vehicle. The NASA Computational Fluid Dynamics (CFD) team has computed free-air flow solutions on the Mach 0.80 design flight configuration and two wind tunnel configuration variants using the USM3D and LAVA flow solvers. Accurate modeling of the configuration tested in the wind tunnel environment is critical to validating the CFD tools, thus the team has included the internal cavity region and sting in their modeling of the configuration, similar to that tested in the 11-Ft TWT. Overall, the CFD simulations compared well and show similar trends as the corrected experimental data for lift and drag polars. The CFD predicted lift curve is shifted in angle of attack from what was observed in the experiment. The shift in lift also was seen in the pitching moment comparison plots. CFD solutions were computed at constant CL test point values and showed overall very good agreement when comparing constant spanwise cuts of pressure coefficient data on the wing and strut with experimental data. CFD cavity corrections were also investigated using the 11-Ft TWT cavity correction method, similar to that used to correct the wind tunnel data. Results showed some improvement in pitching moment coefficient predictions, and an increase in drag, shifting the data to the right in drag polars, further from the experimental data at lower lift conditions, good agreement near the design CL, and slight improvement at the higher lift conditions.

Transonic Truss-Braced Wing↗

A Model for Ice Accretion Roughness Evolution and Spatial Variations

Over the past decade, multiple investigations of ice accretion roughness and spatial variations have been performed in the Icing Research Tunnel (IRT) at the NASA Glenn Research Center. The early investigations used models of NACA 0012 airfoils with different chord sizes and focused on temporal scaling and primary cloud scaling parameters such as stagnation point collection efficiency. Subsequent investigations included the effects of model sweep, airfoil shape, airfoil lifting condition, and freestream static temperature. To develop a predictive model for roughness evolution in generalized icing situations, the maximum roughness values for the airfoils and conditions used in the angle of attack and freestream temperature investigations were scaled to eliminate the temporal variations. LEWICE simulations were employed to identify the local collection efficiency at the locations of maximum roughness, and a two-dimensional panel-method code was used to identify the local static pressure at the location of the maximum roughness. Because of the stochastic nature of ice accretion roughness, a physics-directed approach was employed to develop a multi-dimensional correlation based on the local pressure coefficient, the local total temperature, the cloud properties, and the cloud exposure time. The resulting correlation predictions are compared to ice shapes measured in the IRT for both a 21-in. NACA 0012 airfoil model and a 60-in. HAARP-II model. The resulting predictions indicate that the local freezing fraction must be included in the roughness predictive model. Implications regarding icing heat transfer predictions using the resulting roughness model are also discussed.

Icing↗

Angle-of-Attack Sweep with Mesh Adaptation for High-Lift Configurations

This paper reports an angle-of-attack (alpha) sweep study with mesh adaptation for the high-lift configuration - JAXA Standard Model (JSM) from the 3rd AIAA High-Lift Prediction Workshop. The method of angle-of-attack (alpha) continuation, where the flow-field at a certain angle of attack is initialized from the solution obtained at the previous angle of attack, is applied to adapted unstructured meshes. Since the meshes are continuously adapted, flow-field initialization on a particular mesh for a certain angle is done through linear interpolation of the solution from the previous mesh. Interpolated initial conditions usually prevent massively separated solutions that may otherwise occur with impulsive initial conditions. Alpha continuations are done in both forward (increasing alpha) and reverse (decreasing alpha) directions, resulting in the formation of hysteresis loops. Skin frictions are compared with the wind-tunnel oil flow images and pressure coefficients are compared with the wind-tunnel values, at selected angles of attack. Final paper will include results also for the high-lift version of the Common Research Model (CRM-HL) from the 4th AIAA High-Lift Prediction Workshop. A thorough comparison of the results will also be made with the available wind-tunnel data.

Mesh Adaptation↗

High-Lift Common Research Model: RANS, HRLES and WMLES Perspectives for CLmax Prediction Using LAVA

A unified assessment of three turbulence treatments: Reynolds Averaged Navier-Stokes(RANS), Hybrid RANS/LES (HRLES) and Equilibrium Wall-Modelled Large Eddy Simulation(WMLES) is presented for the High-Lift Common Research Model (CRM-HL). For the free-air configuration, steady-state RANS simulations show very accurate drag polar predictions in the low-𝛼linear regime. However, strong grid sensitivity is reported near the maximum lift-state(𝐶𝐿max), with finer-grids showing larger errors and predicting erroneous flow topologies on the wing. Our RANS simulations show that several corrections for the Spalart-Allmaras (SA)turbulence model widely used in the community lead to more erroneous results compared to the baseline closure, without exception. Both scale-resolving methods (HRLES and WMLES)address these drawbacks and predict an outboard separation pattern on the main element that is in good agreement with the oil flow photographs taken from the QinetiQ wind tunnel experiments, when LES-appropriate grids and numerical discretizations are used. While RANS simulations with the baseline SA closure do not show any wing-root separation post𝐶𝐿max, both HRLES and WMLES show onset of corner flow separation with varying degrees of progression, along with a weak pitch break in the wing-contribution of the overall pitching moment. This post-𝐶𝐿max pitch break seen in the free-air simulations is weaker than the break observed in experiments, with a weaker break reported in WMLES for each iteration of grid-refinement. In-tunnel simulations using both SA-baseline RANS and WMLES show a much stronger post-𝐶𝐿max break with the WMLES predictions showing excellent agreement with the experiment in terms of both the flow-topology observed and the pressure-coefficients at various spanwise stations. Sensitivity to the tunnel wall boundary layer is characterized via comparisons between viscous and inviscid treatments for the tunnel walls. WMLES predictions show moderate sensitivity at the predicted inboard flow-state at 𝐶𝐿max along with the progression towards a post-𝐶𝐿max stall; however, this stalled state at 𝛼≈20◦(inside the tunnel) obtained with both tunnel wall treatments appears to be largely identical.

TTT↗

Exploring Airfoil Table Generation using XFOIL and OVERFLOW

The rotorcraft design process is a continuously evolving field of research that incorporates a number of software programs. An accurate airfoil table is critical in the design and testing process for rotorcraft. With multiple flow solvers available and flow conditions of multirotor UAM vehicles potentially covering a wide range of Reynolds and Mach numbers, a documented approach for developing airfoil tables is needed. Using benchmark data from legacy airfoil tables and wind tunnel tests for comparison, simulations for a comprehensive test matrix could guide rotorcraft design engineers in generating their own airfoil tables using the XFOIL and OVERFLOW solvers. The motivation for this study is to investigate flow solver features to develop a best practices document for airfoil table generation. The study uses the OVERFLOW and XFOIL flow solvers, coupled with the airfoil table generator AFTGen, to analyze three airfoils for a specific Reynolds numbers flow regime and provide details on how well each flow solver performs within a specific angle of attack range, Mach number range, Reynolds number range, and in different flow conditions, such as turbulent and transitional flow. OVERFLOW analyses in AFTGen for fully turbulent and transition flow are compared with XFOIL results and experimental test data for the section lift, section drag, and pressure coefficients. XFOIL ultimately yields results that are accurate within the linear angle of attack range and below a Mach number of 0.4 but tends to overpredict lift and underpredict drag unless the flow is in the compressible regime. XFOIL cannot accurately model stall and post-stall conditions due to the nature of the solver. This is evident in nearly every case run with XFOIL, where the linear range is usually predicted acceptably and the lift coefficient is overpredicted as the stall angle of attack is approached (with the exception being the generally poor correlation with most of the SSC-A09 cases). OVERFLOW is limited at low Mach numbers, and appears to perform best at Mach numbers of 0.4 and above. The exploration of airfoil table generation using XFOIL and OVERFLOW yielded moderately successful results for the NACA 0012 airfoil, reasonably good results for the RC(4)-10 airfoil, and less accurate results for the SSC-A09 airfoil.

Airfoil Table Generation↗

Implicit Large-Eddy Simulation of Axisymmetric and Asymmetric Turbulent Shock Wave Boundary Layer Interaction at Mach 2.5

High-speed inlet flows feature turbulent shock wave boundary layer interactions at supersonic Mach numbers. A canonical Mach 2.5 shock wave boundary layer interaction experiment at the NASA Glenn Research Center was designed for axisymmetric and asymmetric interactions depending on the placement of the shock generator. Four different wall-resolved implicit large-eddy simulations of the experiment were carried out at one-tenth of the Reynolds number of the experiment: (1) Two axisymmetric interactions with shock generator radii of 9.2% and 14.7% of the test section diameter and (2) two asymmetric interactions with 14.7% shock generator radius and shock generator centerline offsets of one-third and one-fourth of the test section radius. For the two axisymmetric interactions, various sub-domains upstream of the interaction, within the interaction, and downstream of the interaction were analyzed with the proper orthogonal decomposition. For all three regions, the dominant modes are streamwise coherent structures with Strouhal numbers based on separation length of 0.15 and below. Fourier spectra of the wall-pressure coefficient for the interaction region have peaks at Strouhal numbers based on separation length of 0.15 and 0.04. For the two asymmetric interactions, the streamwise extent of the separation varies considerably in the azimuthal direction. For both shock generator offsets, the interaction is strongest for the azimuthal location that is closest to the shock generator. Compared to the axisymmetric interaction, the streamwise extent of the separated flow region is noticeably reduced and a substantial crossflow is observed. Fourier analyses of the unsteady flow indicate low-frequency content f or the upstream part of the interaction which is more two-dimensional. Away from the upstream region, with increasing sweep angle and crossflow, the low-frequency content is diminished.

Computational Fluid Dynamics↗