Search NASA⌕ Search

SEARCH · Search NASA

Results for “wall-modeled large-eddy simulation”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

WMLES of K-Type and Bypass Boundary-Layer Transition

A new wall model for transitional flows is applied to the wall-modeled large-eddy simulation of the Klebanoff (fundamental) and bypass transition mechanisms, continuing an effort to improve computational efficiency by mitigating resolution requirements for the disturbances upstream of transition. Comparison with direct-numerical simulation results shows good agreement in the region prior to transition, but agreement during and after transition is not as good as previous wall-modeled large-eddy simulations of the Herbert (subharmonic) transition mechanism using this wall model. This is apparently due to weaknesses in the modelling of the transition process itself.

wall-modeled large-eddy simulation↗
NASA NTRSDate not supplied

Turbulence Simulations of Transonic Flows over NACA-0012 and OAT15A Airfoils

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012, OAT15A airfoils at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 and 3million, respectively. The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

Turbulence Simulations of Transonic Flows over NACA-0012 and OAT15A Airfoils

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012, OAT15A airfoils at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 and 3million, respectively. The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

Turbulence Simulations of Transonic Flows over an NACA-0012 Airfoil

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012 airfoil at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 million. These conditions are the same as the wind-tunnel experimental conditions of McDevitt and Okuno (1985). The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Transonic Buffet↗
NASA NTRSDate not supplied

Turbulence Simulations of Transonic Flows over an NACA-0012 Airfoil

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012 airfoil at different angles of attack covering pre-and post-buffet-onset regimes. The freestream Mach number is 0.75,and the Reynolds number based on the chord length is 10million. These conditions are the same as the wind-tunnel experimental conditions of McDevitt and Okuno (1985). The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre-and post-buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the span wise direction, WMLES predictions show buffeting consistent with the experiment.

Transonic Buffet↗
NASA NTRS2023

Turbulence Simulations of Transonic Flows over an NACA-0012 Airfoil

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES),and wall-modeled large-eddy simulation (WMLES)are employed to simulate transonic flow over an NACA-0012 airfoil at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75,and the Reynolds number based on the chord length is 10million. These conditions are the same as the wind-tunnel experimental conditions of McDevitt and Okuno (1985).1The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model.The unsteady RANS simulations,only with the compressibility correction,predict the pre- and post-buffet characteristics,which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment.The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment

Transonic Buffet↗
NASA NTRS2023

Numerical Study on the Aerodynamics of an Iced Airfoil With Scale-Resolving Simulations

Scale-resolving simulations of the NACA 23012 airfoil with horn ice accretion on the leading edge are conducted using the hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (hybrid RANS/LES) and wall-modeled large-eddy simulation (WMLES) approaches implemented in the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. Aerodynamic results at the Reynolds number of 1.8 million show good comparison with the experimental measurements at different angles of attack from pre-stall to post-stall regimes. The pressure plateaus caused by the flow separation and the recovery of pressure inside the separation bubble around the iced leading edge are well predicted with the scale-resolving simulations when sufficient grid resolution is used around the accreted ice. The unsteadiness of the turbulent flows around the iced airfoil is also examined through the turbulent kinetic energy with the Reynolds normal stress anisotropy. Kelvin–Helmholtz instability (KHI) arises at the shear layer triggered by the upper ice horn and leads to rapid laminar-to-turbulent transition over a large range of angle of attack. With the increase of the angle of attack, the region with high turbulence intensity induced by the unstable shear layer spreads quickly over the entire upper surface of the airfoil. The coherent KHI modes from the upper and lower ice horns are extracted using the spectral proper orthogonal decomposition (SPOD) technique. The SPOD modes extracted from the upper shear layer have large-scale variations in the spanwise direction and low-rank behavior where the energy of the leading SPOD mode at each Strouhal number of the KHI largely represents the total energy when the mode number in the spanwise direction is small.

ARMD↗
NASA NTRS2023

Numerical Study on the Aerodynamics of an Iced Airfoil With Scale-Resolving Simulations

Scale-resolving simulations of the NACA 23012 airfoil with horn ice accretion on the leading edge are conducted using the hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (hybrid RANS/LES) and wall-modeled large-eddy simulation (WMLES) approaches implemented in the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. Aerodynamic results at the Reynolds number of 1.8 million show good comparison with the experimental measurements at different angles of attack from pre-stall to post-stall regimes. The pressure plateaus caused by the flow separation and the recovery of pressure inside the separation bubble around the iced leading edge are well predicted with the scale-resolving simulations when sufficient grid resolution is used around the accreted ice. The unsteadiness of the turbulent flows around the iced airfoil is also examined through the turbulent kinetic energy with the Reynolds normal stress anisotropy. Kelvin–Helmholtz instability (KHI) arises at the shear layer triggered by the upper ice horn and leads to rapid laminar-to-turbulent transition over a large range of angle of attack. With the increase of the angle of attack, the region with high turbulence intensity induced by the unstable shear layer spreads quickly over the entire upper surface of the airfoil. The coherent KHI modes from the upper and lower ice horns are extracted using the spectral proper orthogonal decomposition (SPOD) technique. The SPOD modes extracted from the upper shear layer have large-scale variations in the spanwise direction and low-rank behavior where the energy of the leading SPOD mode at each Strouhal number of the KHI largely represents the total energy when the mode number in the spanwise direction is small.

ARMD↗
NASA NTRSDate not supplied

Summary of the 5th Propulsion Aerodynamics Workshop Nozzle Test Case: Heated Nozzle Exhaust Passing Over A Film-Cooled Plate

This paper summarizes findings from the fifth AIAA Propulsion Aerodynamics Workshop nozzle test case. The experimental configuration examined by workshop participants was a subsonic nozzle with a square exit blowing a heated exhaust over a film-cooled plate. Computational fluid dynamics solutions were obtained and compared to experimental measurements of flow velocities and temperatures, as well as plate surface temperatures. The heated nozzle operated with a Mach 0.3 exit flow at a static temperature ratio of 2.7. Cooling air blowing ratios of 0, 1, and 2 were considered. Computational meshes for the flow domain were provided for participants. Workshop participants from eight separate organizations represented government, industry, and academia. A variety of flow solutions were obtained: Most of the flow solutions employed a Reynolds-averaged Navier-Stokes (RANS) approach, but scale-resolving simulations were used in some cases, including wall-modeled large-eddy simulation (LES), and a few hybrid RANS-LES approaches. One Lattice-Boltzmann solver was employed. For heat transfer to the test article, many participants used a conjugate heat transfer approach. Effects of mesh sensitivity, flow solution approach, and wall heat transfer approach are considered. In general, a fully three-dimensional conjugate heat transfer approach enabled better prediction of surface temperatures than simpler wall temperature boundary treatments. Also, the scale resolving approaches provided better prediction of static temperatures in the flow immediately above the plate consisting of the hot jet exhaust boundary layer interacting with the cooling film. Comparisons of plate surface friction drag and heat transfer obtained from the computations are also presented.

nozzle↗
NASA NTRSDate not supplied

Summary of the 5th Propulsion Aerodynamics Workshop Nozzle Test Case: Heated Nozzle Exhaust Passing Over A Film-Cooled Plate

This paper summarizes findings from the fifth AIAA Propulsion Aerodynamics Workshop nozzle test case. The experimental configuration examined by workshop participants was a subsonic nozzle with a square exit blowing a heated exhaust over a film-cooled plate. Computational fluid dynamics solutions were obtained and compared to experimental measurements of flow velocities and temperatures, as well as plate surface temperatures. The heated nozzle operated with a Mach 0.3 exit flow at a static temperature ratio of 2.7. Cooling air blowing ratios of 0, 1, and 2 were considered. Computational meshes for the flow domain were provided for participants. Workshop participants from eight separate organizations represented government, industry, and academia. A variety of flow solutions were obtained: Most of the flow solutions employed a Reynolds-averaged Navier-Stokes (RANS) approach, but scale-resolving simulations were used in some cases, including wall-modeled large-eddy simulation (LES), and a few hybrid RANS-LES approaches. One Lattice-Boltzmann solver was employed. For heat transfer to the test article, many participants used a conjugate heat transfer approach. Effects of mesh sensitivity, flow solution approach, and wall heat transfer approach are considered. In general, a fully three-dimensional conjugate heat transfer approach enabled better prediction of surface temperatures than simpler wall temperature boundary treatments. Also, the scale resolving approaches provided better prediction of static temperatures in the flow immediately above the plate consisting of the hot jet exhaust boundary layer interacting with the cooling film. Comparisons of plate surface friction drag and heat transfer obtained from the computations are also presented.

nozzle↗
NASA NTRSDate not supplied

WMLES of Boundary-Layer Transition

The potential for the accurate prediction of laminar-turbulent transition by Wall-Modeled Large-Eddy Simulation (WMLES) is investigated through analysis of Direct Numerical Simulation (DNS) and WMLES results. The analysis indicates satisfactory prediction of transition by WMLES is possible, provided a) a sufficiently accurate wall model for the initial stage of transition is employed, b) a satisfactory transition is made between this initial-stage wall model and the wall model for the downstream fully-developed turbulence, and c) limitations on temporal and spatial resolution, and on the thickness of the wall-modeled region, imposed by the need to accurately resolve the disturbance waves and the transition region are observed. This is illustrated with a WMLES computation that yields reasonable agreement with DNS results; however, the wall-modeled region thickness, and thus the cost savings, is less than is typical for WMLES of fully developed turbulence.

wall-modeled large-eddy simulation, boundary-layer↗
NASA NTRSDate not supplied

WMLES of Boundary-Layer Transition

An analysis of the potential for the accurate prediction of laminar-turbulent transition by Wall-Modeled Large-Eddy Simulation (WMLES) indicates that satisfactory prediction of transition by WMLES is possible, provided a) a sufficiently accurate wall model for the initial stage of transition is employed, b) a satisfactory transition is made between this initial-stage wall model and the wall model for the downstream fully-developed turbulence, and c) limitations on temporal and spatial resolution, and on the thickness of the wall-modeled region, imposed by the need to accurately resolve the disturbance waves and the transition region are observed. This is illustrated with WMLES computations of subharmonic transition in a flat-plate boundary layer that yield reasonable agreement with DNS results. However, the wall-modeled region thickness, and thus the cost savings, is less than is typical for WMLES of fully developed turbulence.

wall-modeled large-eddy simulation↗
NASA NTRSDate not supplied

WMLES of Boundary-Layer Transition

An analysis of the potential for the accurate prediction of laminar-turbulent transition by Wall-Modeled Large-Eddy Simulation (WMLES) indicates that satisfactory prediction of transition by WMLES is possible, provided a) a sufficiently accurate wall model for the initial stage of transition is employed, b) a satisfactory transition is made between this initial-stage wall model and the wall model for the downstream fully-developed turbulence, and c) limitations on temporal and spatial resolution, and on the thickness of the wall-modeled region, imposed by the need to accurately resolve the disturbance waves and the transition region are observed. This is illustrated with WMLES computations of subharmonic transition in a flat-plate boundary layer that yield reasonable agreement with DNS results. However, the wall-modeled region thickness, and thus the cost savings, is less than is typical for WMLES of fully developed turbulence.

Wall-modeled large-eddy simulation↗
NASA NTRSDate not supplied

Wall Modeled Large Eddy Simulations for NASA’s Jet Noise Consensus Database of Single-Flow, Round, Convergent Jets

A campaign of wall-modeled large-eddy simulations (WMLES) using structured curvilinear overlapping grids has been performed with the Launch Ascent and Vehicle Aerodynamics(LAVA) computational fluid dynamics (CFD) software to predict jet noise for single-stream axisymmetric round jets. The simulations address the new Prediction Uncertainty Reduction(PUR) technical challenge within the context of NASA’s Commercial Supersonic Technology(CST) project. The goal of PUR is to quantify and reduce uncertainties from scale-resolving simulations to assess noise characteristics of next generation quiet supersonic commercial jets during takeoff and landing conditions where the noise from the exhaust jet dominates. The focus of this effort is to generate a simulation database for single-stream axisymmetric round nozzles at several conditions both for static (no ambient co-flow), which is the focus of this article, and in-flight (flight stream co-flow) conditions, which are beyond the current scope. Nine different flow conditions ranging in jet exit Mach number from 0.38 to 1.0 with nozzle temperature ratios (NTR) from 0.84 to 2.7 have been conducted. Details of the structured overset grids, numerical discretization and wall-model are provided. Near-field comparisons to PIV data show great agreement for both velocity and normal stresses, however a systematic TKE overshoot at the nozzle exit is seen in the lip line shear-layer. A permeable Ffowcs Williams Hawkings (FWH) surface, enclosing the jet, is used to predict far-field noise from the simulated flow-field. Comparison of CFD predictions to microphone array measurements demonstrate excellent agreement within the resolved frequency range. A systematic under-prediction of far-aft observer angles larger than 150 degrees has been observed across all simulations. We achieved a cost reduction of an order of magnitude for these WMLES compared to an earlier study of this configuration due to algorithmic and software improvements. The accuracy of the results and short turnaround time demonstrate that WMLES within the LAVA framework is a cost-effective approach for jet noise predictions that could soon be incorporated into the design cycle of jet noise reduction technologies.

CST↗
NASA NTRSDate not supplied

Wall-Modeled Large Eddy Simulation Method for Unstructured-Grid Navier-Stokes Solvers

This paper reports on the implementation and assessment of a Wall-Modeled Large-Eddy Simulation (WMLES) methodology in an unstructured-grid, node-centered flow solver, FUN3D that is developed and supported at the NASA Langley Research Center. Finite-volume (FV) and finite-element (FE) discretization schemes considered in the study provide formal second-order spatial accuracy. Large-Eddy Simulations (LES) resolve large-scale turbulent-flow features and filter out small-scale effects using the Vreman subgrid-scale model. At solid-wall boundaries, a shear-stress model is employed to provide a proper boundary-flux closure. The nonlinear equations are integrated in time using either an optimized backward difference formula or an implicit multistage Runge-Kutta temporal scheme. The implicit equations at each time step are solved by strong nonlinear iteration schemes. WMLES demonstrations are shown for two high-lift configurations, namely, the McDonnell Douglas 30P30N multielement airfoil and a NASA High-Lift Common Research Model. Results show that the WMLES approaches implemented in the FV and FE discretization methods produce consistent solutions and are capable of capturing key aerodynamic characteristics and flow structures for high-lift configurations at a wide range of angles of attack including maximum-lift conditions. In the 30P30N example, correct trends in the variations of integrated aerodynamic forces and moments, surface pressure distributions, and boundary-layer profiles are captured as the Reynolds number is increased.

CFD; turbulence modeling; High-Lift flow simulatio↗
NASA NTRS2018

Large-Eddy Simulation of Axisymmetric Compression Corner Flow

The Wall-Modeled Large Eddy Simulation (WMLES) approach is used to study the interaction of a shock wave with a high Reynolds number turbulent boundary layer. Since the near wall region is modeled, high Reynolds number turbulent flows can be simulated at a moderate computational cost. The case considered is that of an axisymmetric Mach 2.85 turbulent boundary layer over a 30 compression corner. The Reynolds number of the boundary layer upstream of the interaction based on momentum thickness (Re theta = u sub infinity theta/v sub infinity) is ~12,000. The geometry and flow conditions match the experiments of Dunagan et al. (NASA TM 88227, 1986). The simulations were performed using equilibrium and non-equilibrium wall models. The agreement with experiment is encouraging for the finest grid with respect to the separation bubble length, unsteady shock structure and wall pressure distribution. Sensitivity ofWMLES results to grid, wall model, and blockage effects in the tunnel are reported.

Iyer, Prahladh S.↗
NASA NTRSDate not supplied

Noise Prediction of Multi-stream Internally Mixed Jets with External and Internal Plugs using Large-Eddy Simulation

This paper presents results from ongoing research on jet noise prediction using wall-modeled large eddy simulations (WMLES) conducted within the Launch Ascent and Vehicle Aerodynamics (LAVA) computational framework. The study primarily focuses on the aeroacoustic implications of multi-stream nozzle configurations with internal mixing and both internal and external plugs, at a Reynolds number of 1×10 6 based on the nozzle exit diameter. While internal mixing nozzles have long been considered for jet noise reduction, the complexities of their impact on overall noise levels remain insufficiently understood. This research applies established best practices for WMLES to these complex nozzle configurations, aiming to assess their efficacy and identify limitations in accurately predicting noise behaviors. Through detailed comparisons with experimental data obtained from NASA’s Glenn Research Center, initial findings underscore the challenges inherent in current simulation practices when confronted with intricate geometric and operational conditions. Responding to these challenges, the study explores innovative computational approaches to rectify discrepancies noted between experimental outcomes and CFD predictions. Once example is the introduction of a resonating sound source in the simulation environment to mimic potential unobserved acoustic phenomena.A second example is the strategic modifications to the geometry of external plugs to account for real-world deformations caused by heating and gravity. These novel strategies aim to enhance the accuracy and reliability of noise predictions from supersonic jets, advancing our understanding and capability to effectively reduce jet noise in commercial supersonic aircraft

CST↗
DOE OSTI2024

An enrichment wall modeling framework for spectral element methods

In the present work, a first-of-its-kind enrichment wall-model is developed within the spectral element method (SEM) framework for large-eddy simulations (LES) of wall-bounded turbulent flows. The method augments the polynomial solution in the wall-adjacent elements with an analytical law-of-the-wall enrichment function representing the mean velocity near the wall. In the solution representation, this enrichment function captures the large gradients in the boundary layer, which allows the polynomial modes to represent the turbulent fluctuations. The enriched solution is able to resolve the shear stress at the wall without any modification to the no-slip wall boundary conditions, which allows for greater accuracy in the near-wall region compared to traditional methods. The enrichment wall modeling approach is implemented in a high-order SEM computational fluid dynamics solver, Nek5000, and its performance is assessed in turbulent channel flow wall-modeled LES for a range of Reynolds numbers. It is demonstrated that the enrichment wall-model improves solution accuracy on under-resolved near-wall grids as compared to traditional shear stress wall-models.

42 ENGINEERING↗