Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence modeling”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Turbulence Modeling for Shock Wave/Turbulent Boundary Layer Interactions

Accurate aerodynamic computational predictions are essential for the safety of space vehicles, but these computations are of limited accuracy when large pressure gradients are present in the flow. The goal of the current project is to improve the state of compressible turbulence modeling for high speed flows with shock wave / turbulent boundary layer interactions (SWTBLI). Emphasis will be placed on models that can accurately predict the separated region caused by the SWTBLI. These flows are classified as nonequilibrium boundary layers because of the very large and variable adverse pressure gradients caused by the shock waves. The lag model was designed to model these nonequilibrium flows by incorporating history effects. Standard one- and two-equation models (Spalart Allmaras and SST) and the lag model will be run and compared to a new lag model. This new model, the Reynolds stress tensor lag model (lagRST), will be assessed against multiple wind tunnel tests and correlations. The basis of the lag and lagRST models are to preserve the accuracy of the standard turbulence models in equilibrium turbulence, when the Reynolds stresses are linearly related to the mean strain rates, but create a lag between mean strain rate effects and turbulence when nonequilibrium effects become important, such as in large pressure gradients. The affect this lag has on the results for SWBLI and massively separated flows will be determined. These computations will be done with a modified version of the OVERFLOW code. This code solves the RANS equations on overset grids. It was used for this study for its ability to input very complex geometries into the flow solver, such as the Space Shuttle in the full stack configuration. The model was successfully implemented within two versions of the OVERFLOW code. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWBLI assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Lillard, Randolph P.↗
NASA NTRS1991

An analysis of RNG based turbulence models for homogeneous shear flow

In a recent paper, the authors compared the performance of a variety of turbulence models including the k-epsilon model and the second-order closure model based on Renormalization Group (RNG) Methods. The performance of these RNG models in homogeneous turbulent shear flow was found to be quite poor, apparently due to the value of the constant C(sub epsilon1) in the modeled dissipation rate equation which was substantially lower than its traditional value. However, recently a correction has been made in the RNG based calculation of C(sub epsilon1). It is shown that with the new value of C(sub epsilon1), the performance of the RNG k-epsilon model is substantially improved. On the other hand, while the predictions of the revised RNG second-order closure model are better, some lingering problems still remain which can be easily remedied by the addition of higher order terms.

Speziale, Charles G.↗
NASA NTRS1991

An analysis of RNG-based turbulence models for homogeneous shear flow

In a recent paper, the authors compared the performance of a variety of turbulence models including the k-epsilon model and the second-order closure model based on Renormalization Group (RNG) Methods. The performance of these RNG models in homogeneous turbulent shear flow was found to be quite poor, apparently due to the value of the constant C(sub epsilon 1) in the modeled dissipation rate equation which was substantially lower than its traditional value. However, recently a correction has been made in the RNG-based calculation of C(sub epsilon 1). It is shown that with the new value of C(sub epsilon 1), the performance of the RNG k-epsilon model is substantially improved. On the other hand, while the predictions of the revised RNG second-order closure model are better, some lingering problems still remain which can be easily remedied by the addition of higher order terms.

Speziale, Charles G.↗
NASA NTRS2006

Solver and Turbulence Model Upgrades to OVERFLOW 2 for Unsteady and High-Speed Applications

An implicit unfactored SSOR algorithm has been added to the overset Navier-Stokes CFD code OVERFLOW 2 for unsteady and moving body applications. The HLLEM and HLLC third-order spatial upwind convective flux models have been added for high-speed flow applications. A generalized upwind transport equation has been added for solution of the two-equation turbulence models and the species equations. The generalized transport equation is solved using an unfactored SSOR implicit algorithm. Three hybrid RANS/DES turbulence models have been added for unsteady flow applications. Wall function boundary conditions that include compressibility and heat transfer effects have been also been added to OVERFLOW 2.

Nichols, Robert H.↗
DOE OSTI2023

RANS Simulation of Variable Density Turbulent Round Jets with Coflow using xRAGE Hydrodynamic Code and BHR Turbulence Models

This work for the fiscal year 2023 (FY23) is a continuation of previous efforts to evaluate the BHR turbulence models for their ability to accurately simulate variable density turbulent round jets with coflow. As before, RANS simulations are carried out using the xRAGE hydrodynamic code. The following are some of the previous findings. Israel showed that i) the symmetry boundary conditions for the BHR models in axisymmetric simulations were in error, ii) three grids of different resolutions did not lead to converging solutions, and iii) BHR 3.1 simulation exhibited instabilities and did not reach a steady state. Saenz and Rauenzahn derived and implemented into xRAGE the correct BHR boundary conditions at the symmetry axis. Cline conducted sensitivity studies with various parameters including the BHR models (versions 2, 3.1 and 4), gravity, material pressure, specific heat, initial turbulent kinetic energy and initial turbulent length scale, and found that the largest factor impacting on simulation results was the BHR model version, followed by the initial turbulent length scale. In addition, freeze boundary conditions at the exit and the side wall of the computational domain were used to remove anomalous flow behavior. Cline adjusted the inlet jet velocity, initial turbulent kinetic energy and initial turbulent length scale to obtain the best reasonable match with the experimental data of Charonko and Prestridge. It was found that the BHR 2 and 3.1 models performed in a similar manner, but the BHR 2 model produced much lower levels of density-specific-volume covariance and turbulent kinetic energy. The main focus for the FY23 is to study the effects of computational parameters related to the boundary conditions, mesh, domain size and timestep size. The reasoning behind this is that, unless simulation results can be shown to be reasonably independent from the aforementioned computational parameters, it would be difficult to attribute any discrepancies between simulation and experimental results to turbulence models. This important aspect has largely been overlooked in the previous years, and therefore will be studied comprehensively here. Additionally, effects of varying the initial turbulent length scale will be examined because it was previously identified as a major factor affecting the flow fields.

42 ENGINEERING↗
NASA NTRS2006

Apparent Transition Behavior of Widely-Used Turbulence Models

The Spalart-Allmaras and the Menter SST kappa-omega turbulence models are shown to have the undesirable characteristic that, for fully turbulent computations, a transition region can occur whose extent varies with grid density. Extremely fine two-dimensional grids over the front portion of an airfoil are used to demonstrate the effect. As the grid density is increased, the laminar region near the nose becomes larger. In the Spalart-Allmaras model this behavior is due to convergence to a laminar-behavior fixed point that occurs in practice when freestream turbulence is below some threshold. It is the result of a feature purposefully added to the original model in conjunction with a special trip function. This degenerate fixed point can also cause nonuniqueness regarding where transition initiates on a given grid. Consistent fully turbulent results can easily be achieved by either using a higher freestream turbulence level or by making a simple change to one of the model constants. Two-equation kappa-omega models, including the SST model, exhibit strong sensitivity to numerical resolution near the area where turbulence initiates. Thus, inconsistent apparent transition behavior with grid refinement in this case does not appear to stem from the presence of a degenerate fixed point. Rather, it is a fundamental property of the kappa-omega model itself, and is not easily remedied.

Rumsey, Christopher L.↗
NASA NTRS2007

Apparent Transition Behavior of Widely-Used Turbulence Models

The Spalart-Allmaras and the Menter SST k-omega turbulence models are shown to have the undesirable characteristic that, for fully turbulent computations, a transition region can occur whose extent varies with grid density. Extremely fine two-dimensional grids over the front portion of an airfoil are used to demonstrate the effect. As the grid density is increased, the laminar region near the nose becomes larger. In the Spalart-Allmaras model this behavior is due to convergence to a laminar-behavior fixed point that occurs in practice when freestream turbulence is below some threshold. It is the result of a feature purposefully added to the original model in conjunction with a special trip function. This degenerate fixed point can also cause non-uniqueness regarding where transition initiates on a given grid. Consistent fully turbulent results can easily be achieved by either using a higher freestream turbulence level or by making a simple change to one of the model constants. Two-equation k-omega models, including the SST model, exhibit strong sensitivity to numerical resolution near the area where turbulence initiates. Thus, inconsistent apparent transition behavior with grid refinement in this case does not appear to stem from the presence of a degenerate fixed point. Rather, it is a fundamental property of the k-omega model itself, and is not easily remedied.

Rumsey, Christopher L.↗
DOE OSTI2020

Assessment of Low-Altitude Atmospheric Turbulence Models for Aircraft Aeroelasticity

We investigate the dynamic aeroelastic response of large but slow aircraft in low-altitude atmospheric turbulence. To this end, three turbulence models of increasing fidelity, namely, the one-dimensional von Kármán model, the two-dimensional Kaimal model and full three-dimensional wind fields extracted from large-eddy simulations (LES) are used to simulate ambient turbulence near the ground. Load calculations and flight trajectory predictions are conducted for a representative high-aspect-ratio wing aircraft, using a fully coupled nonlinear flight dynamics/aeroelastic model, when it operates in background atmospheric turbulence generated by the aforementioned models. Comparison of load envelopes and spectral content, on vehicles of varying flexibility, shows strong dependency between the selected turbulence model and aircraft aeroelastic response (e.g. 58% difference in the predicted magnitude of the wing root bending moment between LES and von Kármán models). This is mainly due to the presence of large flow structures at low altitudes that have comparable dimensions to the vehicle, and which despite the relatively small wind speeds within the Earth boundary layer, result in overall high load events for slow-moving vehicles. Results show that one-dimensional models that do not capture those effects provide fairly non-conservative load estimates and are unsuitable for very flexible airframe design.

atmospheric turbulence↗
NASA NTRS2000

Implementation of Advanced Two Equation Turbulence Models in the USM3D Unstructured Flow Solver

USM3D is a widely-used unstructured flow solver for simulating inviscid and viscous flows over complex geometries. The current version (version 5.0) of USM3D, however, does not have advanced turbulence models to accurately simulate complicated flow. We have implemented two modified versions of the original Jones and Launder k-epsilon "two-equation" turbulence model and the Girimaji algebraic Reynolds stress model in USM3D. Tests have been conducted for three flat plate boundary layer cases, a RAE2822 airfoil and an ONERA M6 wing. The results are compared with those from direct numerical simulation, empirical formulae, theoretical results, and the existing Spalart-Allmaras one-equation model.

Wang, Qun-Zhen↗
NASA NTRS1999

Solving Navier-Stokes Equations with Advanced Turbulence Models on Three-Dimensional Unstructured Grids

USM3D is a widely-used unstructured flow solver for simulating inviscid and viscous flows over complex geometries. The current version (version 5.0) of USM3D, however, does not have advanced turbulence models to accurately simulate complicated flows. We have implemented two modified versions of the original Jones and Launder k-epsilon two-equation turbulence model and the Girimaji algebraic Reynolds stress model in USM3D. Tests have been conducted for two flat plate boundary layer cases, a RAE2822 airfoil and an ONERA M6 wing. The results are compared with those of empirical formulae, theoretical results and the existing Spalart-Allmaras one-equation model.

Wang, Qun-Zhen↗
NASA NTRS2001

Turbulence Modeling Effects on the Prediction of Equilibrium States of Buoyant Shear Flows

The effects of turbulence modeling on the prediction of equilibrium states of turbulent buoyant shear flows were investigated. The velocity field models used include a two-equation closure, a Reynolds-stress closure assuming two different pressure-strain models and three different dissipation rate tensor models. As for the thermal field closure models, two different pressure-scrambling models and nine different temperature variance dissipation rate, Epsilon(0) equations were considered. The emphasis of this paper is focused on the effects of the Epsilon(0)-equation, of the dissipation rate models, of the pressure-strain models and of the pressure-scrambling models on the prediction of the approach to equilibrium turbulence. Equilibrium turbulence is defined by the time rate (if change of the scaled Reynolds stress anisotropic tensor and heat flux vector becoming zero. These conditions lead to the equilibrium state parameters. Calculations show that the Epsilon(0)-equation has a significant effect on the prediction of the approach to equilibrium turbulence. For a particular Epsilon(0)-equation, all velocity closure models considered give an equilibrium state if anisotropic dissipation is accounted for in one form or another in the dissipation rate tensor or in the Epsilon(0)-equation. It is further found that the models considered for the pressure-strain tensor and the pressure-scrambling vector have little or no effect on the prediction of the approach to equilibrium turbulence.

Zhao, C. Y.↗
NASA NTRS1995

Evaluation of Turbulence Models for Unsteady Flows of an Oscillating Airfoil

Unsteady flowfields of a two-dimensional oscillating airfoil are calculated using an implicit, finite-difference, Navier Stokes numerical scheme. Five widely used turbulence models are used with the numerical scheme to assess the accuracy and suitability of the models for simulating the retreating blade stall of helicopter rotor in forward flight. Three unsteady flow conditions corresponding to an essentially attached flow, light-stall, and deep-stall cases of an oscillating NACA 0015 wing experiment were chosen as test cases for computations. Results of unsteady airloads hysteresis curves, harmonics of unsteady pressures, and instantaneous flowfield patterns are presented. Some effects of grid density, time-step size, and numerical dissipation on the unsteady solutions relevant to the evaluation of turbulence models are examined. Comparison of unsteady airloads with experimental data show that all models tested are deficient in some sense and no single model predicts airloads consistently and in agreement with experiment for the three flow regimes. The chief findings are that the simple algebraic model based on the renormalization group theory (RNG) offers some improvement over the Baldwin Lomax model in all flow regimes with nearly same computational cost. The one-equation models provide significant improvement over the algebraic and the half-equation models but have their own limitations. The Baldwin-Barth model overpredicts separation and underpredicts reattachment. In contrast, the Spalart-Allmaras model underpredicts separation and overpredicts reattachment.

Srinivasan, G. R.↗
NASA NTRS1983

Turbulence modeling methods for the compressible Navier-Stokes equations

Turbulence modeling methods for the compressible Navier-Stokes equations, including several zero- and two-equation eddy-viscosity models, are described and applied. Advantages and disadvantages of the models are discussed with respect to mathematical simplicity, conformity with physical theory, and numerical compatibility with methods. A new two-equation model is introduced which shows advantages over other two-equation models with regard to numerical compatibility and the ability to predict low-Reynolds-number transitional phenomena. Calculations of various transonic airfoil flows are compared with experimental results. A new implicit upwind-differencing method is used which enhances numerical stability and accuracy, and leads to rapidly convergent steady-state solutions.

Coakley, T. J.↗
NASA NTRS1990

Improvements to a nonequilibrium algebraic turbulence model

It has been noted that while the nonequilibrium turbulence model of Johnson and King (1985, 1987) performed significantly better than alternative methods, differences between predicted and observed shock locations for certain weak interactions are produced due to a defficiency in the model's inner eddy viscosity formulation. A novel formulation for the model is presented which removes this deficiency, while satisfying the law of the wall for adverse pressure-gradient conditions better than either the original formulation or mixing-length theory.

Johnson, D. A.↗
NASA NTRS1993

An algebraic turbulence model for three-dimensional viscous flows

An algebraic turbulence model is proposed for use with three-dimensional Navier-Stokes analyses. It incorporates features of both the Baldwin-Lomax and Cebeci-Smith models. The Baldwin-Lomax model uses the maximum of a function f(y) to determine length and velocity scales. An analysis of the Baldwin-Lomax model shows that f(y) can have a spurious maximum close to the wall, causing numerical problems and non-physical results. The proposed model uses integral relations to determine delta(*) u(sub e) and delta used in the Cebeci-Smith mode. It eliminates a constant in the Baldwin-Lomax model and determines the two remaining constants by comparison to the Cebeci-Smith formulation. Pressure gradient effects, a new wake model, and the implementation of these features in a three-dimensional Navier-Stokes code are also described. Results are shown for a flat plate boundary layer, an annular turbine cascade, and endwall heat transfer in a linear turbine cascade. The heat transfer results agree well with experimental data which shows large variations in endwall Stanton number contours with Reynolds number.

Chima, R. V.↗
NASA NTRS1993

An algebraic turbulence model for three-dimensional viscous flows

An algebraic turbulence model is proposed for use with three-dimensional Navier-Stokes analyses. It incorporates features of both the Baldwin-Lomax and Cebeci-Smith models. The Baldwin-Lomax model uses the maximum of a function f(v) to determine length and velocity scales. An analysis of the Baldwin-Lomax model shows that f(v) can have a spurious maximum close to the wall, causing numerical problems and non-physical results. The proposed model uses integral relations to determine delta(*) u(sub e) and delta used in the Cebeci-Smith mode. It eliminates a constant in the Baldwin-Lomax model and determines the two remaining constants by comparison to the Cebeci-Smith formulation. Pressure gradient effects, a new wake model, and the implementation of these features in a three-dimensional Navier-Stokes code are also described. Results are shown for a flat plate boundary layer, an annular turbine cascade, and endwall heat transfer in a linear turbine cascade. The heat transfer results agree well with experimental data which shows large variations in endwall Stanton number contours with Reynolds number.

Chima, R. V.↗
NASA NTRS2017

Evaluation of Full Reynolds Stress Turbulence Models in FUN3D

Full seven-equation Reynolds stress turbulence models are a relatively new and promising tool for todays aerospace technology challenges. This paper uses two stress-omega full Reynolds stress models to evaluate challenging flows including shock-wave boundary layer interactions, separation and mixing layers. The Wilcox and the SSGLRR full second-moment Reynolds stress models are evaluated for four problems: a transonic two-dimensional diffuser, a supersonic axisymmetric compression corner, a compressible planar shear layer, and a subsonic axisymmetric jet. Simulation results are compared with experimental data and results using the more commonly used Spalart-Allmaras (SA) one-equation and the Menter Shear Stress Transport (SST) two-equation models.

computational fluid dynamics↗
NASA NTRS1993

An algebraic turbulence model for turbomachinery

This paper presents a description and verification of RVC3D (rotor viscous code 3-D) which provides a Euler or Navier-Stokes analysis for steady three dimensional flows in turbomachinery. A motivation for this analysis is the calculation of turbine endwall heat transfer. Features of the turbulence model code include thin-layer formulation, Baldwin-Lomax or Cebeci-Smith turbulence models, node-centered finite difference formulation, and explicit four-stage Runge-Kutta time marching scheme. Results for flat plate, annular turbine cascade, turbine endwall heat transfer, and supersonic compressor blade test cases are presented.

Chima, Rodrick V.↗