Two new two-equation eddy-viscosity turbulence models will be presented. They combine different elements of existing models that are considered superior to their alternatives. The first model, referred to as the baseline (BSL) model, utilizes the original k-ω model of Wilcox in the inner region of the boundary layer and switches to the standard k-e model in the outer region and in free shear flows. It has a performance similar to the Wilcox model, but avoids that model's strong freestream sensitivity

/pdf/two-equation-eddy-viscosity-turbulence-models-for-54j7shl2v6.pdf

Two-equation eddy-viscosity turbulence models for engineering applications

This paper presents a reformulated version of the author'sk-ω model of turbulence. Revisions include the addition of just one new closure coefficient and an adjustment to the dependence of eddy viscosity on turbulence properties. The result is a significantly improved model that applies to both boundary layers and free shear flows and that has very little sensitivity to finite freestream boundary conditions on turbulence properties. The improvements to the k-ω model facilitate a significant expansion of its range of applicability. The new model, like preceding versions, provides accurate solutions for mildly separated flows and simple geometries such as that of a backward-facing step. The model's improvement over earlier versions lies in its accuracy for even more complicated separated flows. This paper demonstrates the enhanced capability for supersonic flow into compression corners and a hypersonic shock-wave/ boundary-layer interaction. The excellent agreement is achieved without introducing any compressibility modifications to the turbulence model.

Formulation of the k-w Turbulence Model Revisited

Reynolds stress models and traditional large-eddy simulations are reexamined with a view toward developing a combined methodology for the computation of complex turbulent flows. More specifically, an entirely new approach to time-dependent Reynolds-averaged Navier-Stokes (RANS) computations and very large-eddy simulations (VLES) is presented in which subgrid scale models are proposed that allow a direct numerical simulation (DNS) to go continuously to a RANS computation in the coarse mesh/infinite Reynolds number limit. In between these two limits, we have a large eddy simulation (LES) or VLES, depending on the level of resolution. The Reynolds stress model that is ultimately recovered in the coarse mesh/infinite Reynolds number limit has built in nonequilibrium features that make it suitable for time-dependent RANS. The fundamental technical issues associated with this new approach, which has the capability of bridging the gap between DNS, LES and RANS, are discussed in detail. Illustrative calculations are presented along with a discussion of the future implications of these results for the simulation of the turbulent flows of technological importance.

Turbulence modeling for time-dependent RANS and VLES : a review

The Wilcox (1988, 1991) k-omega model for eddy-viscosity turbulence does not require damping functions in the viscous sublayer, and its equations are less stiff near the wall. It has been designed to predict the requisite wake length in equilibrium, adverse pressure-gradient boundary-layer flows. When applied to free shear layers, however, a strong dependency of its results on the freestream value of omega has been noted. This feature is presently investigated via the self-similar equations for incompressible equilibrium boundary layers.

Influence of freestream values on k-omega turbulence model predictions

The interaction of a spatially developing adiabatic boundary layer flow at M∞=2.25 and Reθ=3725 with an impinging oblique shock wave (β=33.2°) is analyzed by means of direct numerical simulation of the compressible Navier-Stokes equations. Under the selected flow conditions the incoming boundary layer undergoes mild separation due to the adverse pressure gradient. Coherent structures are shed near the average separation point and the flow field exhibits large-scale low-frequency unsteadiness. The formation of the mixing layer is primarily responsible for the amplification of turbulence, which relaxes to an equilibrium state past the interaction. Complete equilibrium is attained in the inner part of the boundary layer, while in the outer region the relaxation process is incomplete. Far from the interaction zone, turbulence exhibits a universal behavior and it shows similarities with the incompressible case. The interaction of the coherent structures with the incident shock produces acoustic waves that prop...

Direct numerical simulation of impinging shock wave/turbulent boundary layer interaction at M=2.25

A variety of two-equation turbulence models,including several versions of the K-epsilon model as well as the K-omega model, are analyzed critically for near wall turbulent flows from a theoretical and computational standpoint. It is shown that the K-epsilon model has two major problems associated with it: the lack of natural boundary conditions for the dissipation rate and the appearance of higher-order correlations in the balance of terms for the dissipation rate at the wall. In so far as the former problem is concerned, either physically inconsistent boundary conditions have been used or the boundary conditions for the dissipation rate have been tied to higher-order derivatives of the turbulent kinetic energy which leads to numerical stiffness. The K-omega model can alleviate these problems since the asymptotic behavior of omega is known in more detail and since its near wall balance involves only exact viscous terms. However, the modeled form of the omega equation that is used in the literature is incomplete-an exact viscous term is missing which causes the model to behave in an asymptotically inconsistent manner. By including this viscous term and by introducing new wall damping functions with improved asymptotic behavior, a new K-tau model (where tau is identical with 1/omega is turbulent time scale) is developed. It is demonstrated that this new model is computationally robust and yields improved predictions for turbulent boundary layers.

/pdf/critical-evaluation-of-two-equation-models-for-near-wall-24uja9ue6i.pdf

Critical Evaluation of Two-Equation Models for Near-Wall Turbulence

Evaluation of two-equation turbulence models for predicting transitional flows

/pdf/a-critical-evaluation-of-two-equation-models-for-near-wall-n6e20avt00.pdf

A critical evaluation of two-equation models for near wall turbulence

An explicit algebraic stress equation, developed by Gatski and Speziale, is used in the framework of the K-UJ and K-e formulations to predict two-dimensional separated turbulent flows. The nonequilibrium effects are modeled through coefficients that depend nonlinearly on both rotational and irrotational strains. The proposed models were implemented in the Courant-Friedrichs-Lewy three-dimensional Navier-Stokes code. Comparisons with the experimental data are presented, which clearly demonstrate that explicit algebraic stress models improve the response of standard two-equation models to nonequilibrium effects. I. Introduction

Prediction of Nonequilibrium Turbulent Flows with Explicit Algebraic Stress Models

Turbulent channel flow and homogeneous shear flow have served as basic building block flows for the testing and calibration of Reynolds stress models. A direct theoretical connection is made between homogeneous shear flow in equilibrium and the log-layer of fully-developed turbulent channel flow. It is shown that if a second-order closure model is calibrated to yield good equilibrium values for homogeneous shear flow it will also yield good results for the log-layer of channel flow provided that the Rotta coefficient is not too far removed from one. Most of the commonly used second-order closure models introduce an ad hoc wall reflection term in order to mask deficient predictions for the log-layer of channel flow that arise either from an inaccurate calibration of homogeneous shear flow or from the use of a Rotta coefficient that is too large. Illustrative model calculations are presented to demonstrate this point which has important implications for turbulence modeling.

/pdf/predicting-equilibrium-states-with-reynolds-stress-closures-24hwzdng0b.pdf

Ridha Abid

Papers

Critical Evaluation of Two-Equation Models for Near-Wall Turbulence

Evaluation of two-equation turbulence models for predicting transitional flows

A critical evaluation of two-equation models for near wall turbulence

Prediction of Nonequilibrium Turbulent Flows with Explicit Algebraic Stress Models

Predicting equilibrium states with Reynolds stress closures in channel flow and homogeneous shear flow