Search NASASearch

SEARCH · Search NASA

Results for “Wall Model”

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 163 records · Page 9

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

On the calculation of turbulent heat and mass transport downstream from an abrupt pipe expansion

A numerical study is reported of heat/mass transfer in the separated flow region created by an abrupt pipe expansion. Computations have employed a hybrid method of central and upwind finite differencing to solve the full Navier-Stokes equations with turbulent model (k approximately equal to epsilon). The study has given its main attention to the simulation of the region in the immediate vicinity of the wall, by formulating near-wall model for the evaluation of the mean generation and destruction rate of the epsilon equation. The computed results were compared with the experimental data and they showed generally encouraging agreement with the measurements.

Amano, R. S.

Wall-temperature effects on the aerodynamics of a hydrogen-fueled transport concept in Mach 8 blowdown and shock tunnels

Results are presented from two separate tests on the same blended wing-body hydrogen fueled transport model at a Mach number of about 8 and a range of Reynolds numbers (based on theoretical body length) of 0.597 x 10 to the 6th power to about 156.22 x 10 to the 6th power. Tests were made in conventional hypersonic blowdown tunnel and a hypersonic shock tunnel at angles of attack of -2 deg to about 8 deg, with an extensive study made at a constant angle of attack of 3 deg. The model boundary-layer flow varied from laminar at the lower Reynolds numbers to predominantly turbulent at the higher Reynolds numbers. Model wall temperatures and stream static temperatures varied widely between the two tests, particularly at the lower Reynolds numbers. These temperature differences resulted in marked variations of the axial-force coefficients between the two tests, due in part to the effects of induced pressure and viscous interaction variations. The normal-force coefficient was essentially independent of Reynolds number. Analysis of results utilized current theoretical computer programs and basic boundary-layer theory.

Penland, J. A.

Sensitivity to Snapshot Frequency in the POD-Based Reduced-Order Modeling of Flow over a Gaussian Bump

The flow over a Gaussian bump at a Reynolds number of 10^6, based on the bump length and freestream conditions, and a Mach number of 0.2, was simulated using wall-modeled large eddy simulation (WMLES). The incoming turbulent flow exhibits relaminarization upstream of the bump apex, due to the favourable pressure gradient acceleration induced by the bump curvature, as well as separation, reattachment and a detached shear layer downstream of the apex. A Galerkin reduced-order model (ROM) was then constructed using a proper orthogonal decomposition (POD) modal basis. The complex, three-dimensional flow features are captured by spatial POD modes differently based on the frequency of the high-fidelity snapshots originally provided. POD modes constructed from differently sampled snapshot matrices were then truncated to retain 90% of the turbulence kinetic energy in the domain of interest. A parametric study is conducted to evaluate the effect of the frequency of snapshots used to construct the POD-ROM.

turbulence

Comparison Study of Conventional and Adaptive Mesh Refinement in Organic Material Decomposition Models

This study compares conventional mesh refinement techniques, specifically Uniform Mesh Refinement (UMR), with a new Adaptive Mesh Refinement (AMR) method, applied to Organic Material Decomposition (OMD) models. The proposed benefit of AMR is that only areas that require refinement, based on minimizing a specific field gradient, are refined thus decreasing model wall time compared to conventional UMR methods. This work specifically focuses on comparing UMR and AMR methods on decomposing (both No-Flow and Porous-Flow material models) Polymeric Methylene Diisocyanate (PMDI) polyurethane foam. Throughout the work, the geometry increased in complexity to assess the refinement methods performance at varying levels geometric intricacy. While AMR has been shown to work well in a variety of applications, the UMR approach proved to be computationally faster, for many of the geometries and foam decomposition models, than AMR. However, it was observed that at higher levels of refinement, greater than 3 UMR, AMR begins to be computationally better. Additionally, the settings used to perform AMR greatly impact its performance, and lessons learned, in terms of OMD models, are shared. Due to physics involved in material decomposition, specifically the evolution of state variables, these problems don’t fully benefit from the advantages of AMR.

36 MATERIALS SCIENCE

Annual Research Briefs - 1996, Center for Turbulence Research

Partial contents include: (1) Simulation of jet impingement heat transfer with the k-e-V2 model; (2) A new methodology for turbulence modelers using Direct Numerical Simulation (DNS) database analysis; (3) A particle representation model for the deformation of homogeneous turbulence; (4) Modeling of inhomogeneous compressible turbulence using a two-scale statistical theory; (5) Lewis number effects on partially premixed flames; (6) Conditional moment closure of mixing and reaction in turbulent non-premixed combustion; (7) Dynamic models for LES of turbulent front propagation with a spectral method; (8) Jet noise using large eddy simulation; (9) Towards numerical simulations of trailing-edge aeroacoustics; (10) A model for the onset of vortex breakdown; (11) Compressing an elliptic vortex: transition to turbulence by tumble breakdown; (12) A new method for the adaptive control of vortex-wall interactions; (13) Symmetries in turbulent boundary layer flows; (14) Near-wall models in large eddy simulations of flow behind a backward-facing step; (15) Assessment of the subgrid-scale models at low and high Reynolds numbers; (16) Large-eddy simulations of flow around a NACA 4412 airfoil using unstructured grids; (17) Progress on LES of flow past a circular cylinder; (18) Experimental investigation of flow through an asymmetric plane diffuser; (19) Progress in the large-eddy simulation of an asymmetric plane diffuser; (20) A numerical study of self-similarity in a turbulent plane wake using large-eddy simulation; (21) On the generation of turbulent inflow conditions for boundary layer simulations, and (22) A new class of finite difference schemes.

COMPUTATIONAL FLUID DYNAMICS

Turbulent separated flow past a backward-facing step: A critical evaluation of two-equation turbulence models

The ability of two-equation models to accurately predict separated flows is analyzed from a combined theoretical and computational standpoint. Turbulent flow past a backward facing step is chosen as a test case in an effort to resolve the variety of conflicting results that were published during the past decade concerning the performance of two-equation models. It is found that the errors in the reported predictions of the k-epsilon model have two major origins: (1) numerical problems arising from inadequate resolution, and (2) inaccurate predictions for normal Reynolds stress differences arising from the use of an isotropic eddy viscosity. Inadequacies in near wall modelling play a substantially smaller role. Detailed calculations are presented which strongly indicate the standard k-epsilon model - when modified with an independently calibrated anisotropic eddy viscosity - can yield surprisingly good predictions for the backstep problem.

Thangam, S.

Turbulent flow past a backward-facing step - A critical evaluation of two-equation models

The ability of two-equation models to accurately predict separated flows is analyzed from a combined theoretical and computational standpoint. Turbulent flow past a backward facing step is chosen as a test case in an effort to resolve the variety of conflicting results that were published during the past decade concerning the performance of two-equation models. It is found that the errors in the reported predictions of the k-epsilon model have two major origins: (1) numerical problems arising from inadequate resolution, and (2) inaccurate predictions for normal Reynolds stress differences arising from the use of an isotropic eddy viscosity. Inadequacies in near wall modelling play a substantially smaller role. Detailed calculations are presented which strongly indicate the standard k-epsilon model - when modified with an independently calibrated anisotropic eddy viscosity - can yield surprisingly good predictions for the backstep problem.

Thangam, S.

A near-wall two-equation model for compressible turbulent flows

A near-wall two-equation turbulence model of the K - epsilon type is developed for the description of high-speed compressible flows. The Favre-averaged equations of motion are solved in conjunction with modeled transport equations for the turbulent kinetic energy and solenoidal dissipation wherein a variable density extension of the asymptotically consistent near-wall model of So and co-workers is supplemented with new dilatational models. The resulting compressible two-equation model is tested in the supersonic flat plate boundary layer - with an adiabatic wall and with wall cooling - for Mach numbers as large as 10. Direct comparisons of the predictions of the new model with raw experimental data and with results from the K - omega model indicate that it performs well for a wide range of Mach numbers. The surprising finding is that the Morkovin hypothesis, where turbulent dilatational terms are neglected, works well at high Mach numbers, provided that the near wall model is asymptotically consistent. Instances where the model predictions deviate from the experiments appear to be attributable to the assumption of constant turbulent Prandtl number - a deficiency that will be addressed in a future paper.

Zhang, H. S.

A near-wall two-equation model for compressible turbulent flows

A near-wall two-equation turbulence model of the k-epsilon type is developed for the description of high-speed compressible flows. The Favre-averaged equations of motion are solved in conjunction with modeled transport equations for the turbulent kinetic energy and solenoidal dissipation wherein a variable density extension of the asymptotically consistent near-wall model of So and co-workers is supplemented with new dilatational models. The resulting compressible two-equation model is tested in the supersonic flat plate boundary layer - with an adiabatic wall and with wall cooling - for Mach numbers as large as 10. Direct comparisons of the predictions of the new model with raw experimental data and with results from the K-omega model indicate that it performs well for a wide range of Mach numbers. The surprising finding is that the Morkovin hypothesis, where turbulent dilatational terms are neglected, works well at high Mach numbers, provided that the near wall model is asymptotically consistent. Instances where the model predictions deviate from the experiments appear to be attributable to the assumption of constant turbulent Prandtl number - a deficiency that will be addressed in a future paper.

Zhang, H. S.

Dynamic Smagorinsky Modeled Large-Eddy Simulations of Turbulence Using Tetrahedral Meshes

Eddy-resolving numerical computations of turbulent flows are emerging as viable alternatives to Reynolds Averaged Navier-Stokes (RANS) calculations for flows with an intrinsically steady mean state due to the advances in large-scale parallel computing. In these computations, medium to large turbulent eddies are resolved by the numerics while the smaller or subgrid scales are either modeled or taken care of by the inherent numerical dissipation. To advance the state of the art of unstructured-mesh turbulence simulation capabilities, large eddy simulations (LES) using the dynamic Smagorinsky model (DSM) on tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method. In contrast to what has been reported in the literature, the present implementation of dynamic models allows for active backscattering without any ad-hoc limiting of the eddy viscosity calculated from the subgrid-scale model. For the benchmark problems involving compressible isotropic turbulence decay as well as the shock/turbulent boundary layer interaction benchmark problems, no numerical instability associated with kinetic energy growth is observed and the volume percentage of the backscattering portion accounts for about 38-40% of the simulation domain. A slip-wall model in conjunction with the implemented DSM is used to simulate a relatively high Reynolds number Mach 2.85 turbulent boundary layer over a 30° ramp with several tetrahedral meshes and a wall-normal spacing of either Δγ& = 10 or Δγ& = 20. The computed mean wall pressure distribution, separation region size, mean velocity profiles, and Reynolds stress agree reasonably well with experimental data.

CESE method

Large-Eddy Simulation Study of Flow and Combustion Dynamics in a Full-Scale Hydrogen–Air Rotating Detonation Combustor-Stator Integrated System

In the present work, a first-of-its-kind three-dimensional (3D) large-eddy simulation (LES) study is conducted to numerically investigate the combustion dynamics as well as aero-thermal phenomena in a full-scale nonpremixed hydrogen–air rotating detonation engine (RDE) (with a diverging-shaped lower-end wall), when integrated with nozzle guide vanes (NGV) acting as the turbine stator. The wall-modeled LES framework incorporates hydrogen–air detailed chemical kinetics and adaptive mesh refinement (AMR). A comparative analysis is carried out for two operating conditions with different fuel/air mass flow rates but global equivalence ratio of unity, and considering RDE configurations without and with stator. The LES model is validated against available experimental data for the low mass flux condition with respect to detonation wave speed/height, wave dynamics, and axial static pressure distribution. Numerical results indicate significant deflagrative combustion occurring in the fill region near the inner wall due to formation of recirculation zones in the injection near-field driven by the backward facing step. The leading detonation wave is found to be trailed by an azimuthal reflected-shock combustion (ARSC) wave, consistent with experimental observations, which consumes unburned vitiated reactants that leak through the main detonation wave. The main detonation wave characteristics, such as detonation wave speed/height and combustion efficiency, do not change appreciably with the presence of NGV. A novel combustion diagnostic technique based on chemical explosive mode analysis (CEMA) is employed to quantify the fraction of heat release occurring in the detonative mode versus deflagrative mode for the simulated conditions. The exit flow is found to be nearly fully subsonic and supersonic for the low and high mass flux conditions, respectively. Further analysis of the exit flow profiles shows that the presence of NGV renders the flow more axial and significantly impacts the exit Mach number and total pressure, while the total temperature shows negligible change. In addition, the low mass flux operating point, despite exhibiting more deflagrative losses within the combustor, yields overall lower pressure drop from plenum to exhaust, which is mainly attributed to lower pressure drop across the injectors. Lastly, the rotating detonation engine-nozzle guide vanes (RDE-NGV) configuration exhibits higher total pressure loss compared to rotating detonation engine (RDE) without stator across both the mass flux conditions. In conclusion, this study extends the state-of-the-art in numerical modeling of pressure gain combustion (PGC) systems by demonstrating high-fidelity 3D reacting LES of full-scale RDE-NGV systems relevant to RDE-turbine integration for stationary power generation.

Combustion

Numerical Investigation of C L,max Prediction on the NASA High-Lift Common Research Model

An assessment of a Hybrid RANS/LES (HRLES) approach for C L,max prediction is presented for the NASA High-Lift Common Research Model (CRM-HL). Both free air and wind tunnel configurations of the CRM-HL are investigated and the results are compared to the QinetiQ wind tunnel experiments and two other numerical approaches: Reynolds Averaged Navier-Stokes (RANS) and Wall-Modeled Large Eddy Simulations (WMLES). For the free-air configuration, HRLES is shown to address some of the known shortcomings with RANS and prevent inboard and outboard flow separation particularly in the region of C L,max and post-stall. To achieve these improvements over RANS, LES-appropriate grids and numerical discretizations are required. It was also found that when applying HRLES to a RANS best practice grid and numerics that the HRLES method significantly underperformed RANS. For the in-tunnel configuration, HRLES showed good agreement with the loads, surface pressure and oil-flow photographs obtained in the experiment. HRLES was able to improve upon the RANS simulations, which showed a sharp loss of lift at the two highest angles-of-attack due to large scale inboard and outboard separation on the wing, by correctly predicting the corner flow separation and showing remarkably close agreement in the flow topologies with the experiment.

TTT

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

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

A Hybrid RANS-LES Perspective for the High Lift Common Research Model Using LAVA

An assessment of a Hybrid RANS/LES (HRLES) approach for𝐶𝐿,max prediction is presented for the NASA High-Lift Common Research Model (CRM-HL). Both the free air and the wind tunnel configuration of the CRM-HL are investigated and the results are compared to the QinetiQ wind tunnel experiments and to two other numerical approaches: Reynolds Averaged Navier-Stokes (RANS) and Wall-Modeled Large Eddy Simulations (WMLES). For the free-air configuration, HRLES was shown to address some of the known shortcomings in RANS methods and prevent inboard and outboard flow separation particularly in the region of𝐶𝐿,max and post-stall. To achieve these improvements over RANS, LES-appropriate grids and numerical discretizations are required. HRLES predicts a weak pitch break at the highest angle-of-attack due to onset of wing-root corner flow separation whereas the free-air corrected experiment values indicate an occurrence of a much stronger pitch break. The improvements of HRLES over a URANS approach has been objectively shown by computing a set of solutions with the same grid, same numerics and time-step size and comparing the solutions. It was also found that when applying HRLES to a RANS best practice grid and numerics that the HRLES method significantly under performed RANS. For the in tunnel configuration, HRLES showed good agreement with the loads, surface pressure and oil-flow photographs obtained in the experiment. HRLES was able to improve upon the RANS simulations, which showed a sharp loss of lift at the two highest angles-of-attack due to large scale inboard and outboard separation on the wing, by correctly predicting the corner flow separation and showing remarkably close agreement in the flow topologies with the experiment.

TTT

Improved Models for the Porous Surface with Passive Control

A computational investigation of afterbody flow using a passive control method is conducted. The passive control method consists of a porous surface placed over a plenum. The purpose of the passive control method is to exploit the adverse pressure gradient present in afterbody flow in an attempt to reduce boundary layer separation and afterbody drag. Four different porous wall models are used to model the transpiration velocity in the region of passive control. A three-dimensional, time-dependent, Reynolds-averaged, simplified Navier-Stokes solver, PAB3D, is used to simulate afterbody flow with and without passive control. Three afterbody configurations with boat-tail angles of 10, 20, and 30 deg. are used to obtain two-dimensional solutions with a freestream Mach number of 0.6 and nozzle pressure ratio of 6. The region of passive control was initially placed from 20-60% of the nozzle length. The effect of the porous placement and porous extent is also studied. Baseline (no porosity) two-dimensional solutions are qualitatively similar to experimental data but under-predict the magnitude of the pressure recovery. Results for the subsonic solutions show losses in the pressure recovery for some cases with passive control. Three-dimensional effects are also investigated and seen to be very significant. Three-dimensional baseline solutions, for both sub- and super-critical freestream Mach numbers, compare very favorably with the experimental data in comparison to the two-dimensional solution. Future work is required to examine three-dimensional afterbody flows with passive porosity.

Kraushaar, Sandra L.

Hybrid Reynolds-Averaged/Large Eddy Simulation of a Cavity Flameholder; Assessment of Modeling Sensitivities

Steady-state and scale-resolving simulations have been performed for flow in and around a model scramjet combustor flameholder. The cases simulated corresponded to those used to examine this flowfield experimentally using particle image velocimetry. A variety of turbulence models were used for the steady-state Reynolds-averaged simulations which included both linear and non-linear eddy viscosity models. The scale-resolving simulations used a hybrid Reynolds-averaged / large eddy simulation strategy that is designed to be a large eddy simulation everywhere except in the inner portion (log layer and below) of the boundary layer. Hence, this formulation can be regarded as a wall-modeled large eddy simulation. This effort was undertaken to formally assess the performance of the hybrid Reynolds-averaged / large eddy simulation modeling approach in a flowfield of interest to the scramjet research community. The numerical errors were quantified for both the steady-state and scale-resolving simulations prior to making any claims of predictive accuracy relative to the measurements. The steady-state Reynolds-averaged results showed a high degree of variability when comparing the predictions obtained from each turbulence model, with the non-linear eddy viscosity model (an explicit algebraic stress model) providing the most accurate prediction of the measured values. The hybrid Reynolds-averaged/large eddy simulation results were carefully scrutinized to ensure that even the coarsest grid had an acceptable level of resolution for large eddy simulation, and that the time-averaged statistics were acceptably accurate. The autocorrelation and its Fourier transform were the primary tools used for this assessment. The statistics extracted from the hybrid simulation strategy proved to be more accurate than the Reynolds-averaged results obtained using the linear eddy viscosity models. However, there was no predictive improvement noted over the results obtained from the explicit Reynolds stress model. Fortunately, the numerical error assessment at most of the axial stations used to compare with measurements clearly indicated that the scale-resolving simulations were improving (i.e. approaching the measured values) as the grid was refined. Hence, unlike a Reynolds-averaged simulation, the hybrid approach provides a mechanism to the end-user for reducing model-form errors.

Baurle, R. A.