Cubic eddy‐viscosity turbulence models for strongly swirling confined flows with variable density

An investigation on the predictive performance of cubic eddy‐viscosity turbulence models for strongly swirling confined flows with variable density is presented. Comparisons of the prediction with the experiments show some improvements of cubic models over the linear k–ε model. The linear k–ε model does not contain any mechanism to represent the interaction of swirl and density variation and as a consequence it performs poorly. With appropriate modelling, two‐equation cubic turbulence models can capture the subcritical nature of the flow, represent the azimuthal velocity profiles of combined forced‐free vortex motion, and predict the combined effects of swirl and density variation fairly well. However, the calibration of model coefficients is still a topic of investigation. Further amendments are also needed for the equations of k and ε to take into account the effects of swirl and density gradients correctly. Copyright © 2004 John Wiley & Sons, Ltd.     

Prediction of developing turbulent flow in 90°‐curved ducts using linear and non‐linear low‐Re k–ε models

This paper reports the outcome of applying two different low‐Reynolds‐number eddy‐viscosity models to resolve the complex three‐dimensional motion that arises in turbulent flows in ducts with 90^° bends. For the modelling of turbulence, the Launder and Sharma low‐Re k–ε model and a recently produced variant of the cubic non‐linear low‐Re k–ε model have been employed. In this paper, developing turbulent flow through two different 90^° bends is examined: a square bend, and a rectangular bend with an aspect ratio of 6. The numerical results indicate that for the bend of square cross‐section the curvature induces a strong secondary flow, while for the rectangular cross‐section the secondary motion is confined to the corner regions. For both curved ducts, the secondary motion persists downstream of the bend and eventually slowly disappears. For the bend of square cross‐section, comparisons indicate that both turbulence models can produce reasonable predictions. For the bend of rectangular cross‐section, for which a wider range of data is available, while both turbulence models produce satisfactory predictions of the mean flow field, the non‐linear k–ε model returns superior predictions of the turbulence field and also of the pressure and friction coefficients. Copyright © 2006 John Wiley & Sons, Ltd.     

Turbulent flow simulation of liquid jet emanating from pressure-swirl atomizer

The performances of three linear eddy viscosity models (LEVM) and one algebraic Reynolds stress model (ARSM) for the simulation of turbulent flow inside and outside pressure-swirl atomizer are evaluated by comparing the interface position with available experimental data and by comparing the turbulence intensity profiles at the atomizer exit. It is found that the turbulence models investigated exhibit zonal behaviors, i.e. none of the models investigated performs well throughout the entire flow field. The turbulence intensity has a significant influence on the global characteristics of the flow field. The turbulence models with better predictions of the turbulence intensity, such as Gatski-Speziale’s ARSM model, can yield better predictions of the global characteristics of the flow field, e.g. the reattachment lengths for the backward-facing step flow and the sudden expansion pipe flow, or the discharge coefficient, film thickness and the liquid sheet outer surface position for the atomizer flows. The standard k–ε model predicts stronger turbulence intensity as compared to the other models and therefore yields smaller film thickness and larger liquid sheet outer surface position. In average, the ARSM model gives both quantitatively and qualitatively better results as compared to the standard k–ε model and the low Reynolds number models.     

A two‐scale low‐Reynolds number turbulence model

In this study, a two‐scale low‐Reynolds number turbulence model is proposed. The Kolmogorov turbulence time scale, based on fluid kinematic viscosity and the dissipation rate of turbulent kinetic energy (ν, ε), is adopted to address the viscous effects and the rapid increasing of dissipation rate in the near‐wall region. As a wall is approached, the turbulence time scale transits smoothly from a turbulent kinetic energy based (k, ε) scale to a (ν, ε) scale. The damping functions of the low‐Reynolds number models can thus be simplified and the near‐wall turbulence characteristics, such as the ε distribution, are correctly reproduced. The proposed two‐scale low‐Reynolds number turbulence model is first examined in detail by predicting a two‐dimensional channel flow, and then it is applied to predict a backward‐facing step flow. Numerical results are compared with the direct numerical simulation (DNS) budgets, experimental data and the model results of Chien, and Lam and Bremhorst respectively. It is proved that the proposed two‐scale model indeed improves the predictions of the turbulent flows considered. Copyright © 2000 John Wiley & Sons, Ltd.     

Evaluation of one‐ and two‐equation low‐Re turbulence models. Part I—Axisymmetric separating and swirling flows

This first segment of the two‐part paper systematically examines several turbulence models in the context of three flows, namely a simple flat‐plate turbulent boundary layer, an axisymmetric separating flow, and a swirling flow. The test cases are chosen on the basis of availability of high‐quality and detailed experimental data. The tested turbulence models are integrated to solid surfaces and consist of: Rodi's two‐layer k–ε model, Chien's low‐Reynolds number k–ε model, Wilcox's k–ω model, Menter's two‐equation shear‐stress‐transport model, and the one‐equation model of Spalart and Allmaras. The objective of the study is to establish the prediction accuracy of these turbulence models with respect to axisymmetric separating flows, and flows of high streamline curvature. At the same time, the study establishes the minimum spatial resolution requirements for each of these turbulence closures, and identifies the proper low‐Mach‐number preconditioning and artificial diffusion settings of a Reynolds‐averaged Navier–Stokes algorithm for optimum rate of convergence and minimum adverse impact on prediction accuracy. Copyright © 2003 John Wiley & Sons, Ltd.     

Near‐wall turbulence modelling with enhanced dissipation

A low‐Reynolds number k–ε turbulence model is proposed that incorporates diffusion terms and modified C_ε(1,2) coefficients to amplify the level of dissipation in non‐equilibrium flow regions, thus reducing the kinetic energy and length scale magnitudes to improve prediction of adverse pressure gradient flows, involving flow separation and reattachment. Unlike the conventional k–ε model, it requires no wall function/distance parameter that bridges the near‐wall integration. The model is validated against a few flow cases, yielding predictions in good agreement with the direct numerical simulation (DNS) and experimental data. Copyright © 2003 John Wiley & Sons, Ltd.     

SOME VALIDATION OF STANDARD,MODIFIED AND NON-LINEAR k–ε TURBULENCE MODELS

Standard, modified and non-linear k–ε: turbulence models are validated against three axisymmetric flow problems—flow through a pipe expansion, flow through a pipe constriction and an impinging jet problem—to underpin knowledge about the solution quality obtained from two-equation turbulence models. The extended models improve the prediction of turbulence as a flow approaches a stagnation point and the non-linear model allows for the prediction of anisotropic turbulence. Significantly different values for the non-linear model coefficients are proposed in comparison with values found in the literature. Nevertheless, current turbulence models are still unable to accurately predict the spreading rate of shear layers. © 1997 by John Wiley & Sons, Ltd. Int. j. numer. methods fluids, 24: 965–986, 1997.     

Numerical simulations of non‐equilibrium turbulent boundary layer flowing over a bump

Large‐eddy simulation (LES) and Reynolds‐averaged Navier–Stokes simulation (RANS) with different turbulence models (including the standard k?ε, the standard k?ω, the shear stress transport k?ω (SST k?ω), and Spalart–Allmaras (S–A) turbulence models) have been employed to compute the turbulent flow of a two‐dimensional turbulent boundary layer over an unswept bump. The predictions of the simulations were compared with available experimental measurements in the literature. The comparisons of the LES and the SST k?ω model including the mean flow and turbulence stresses are in satisfied agreements with the available measurements. Although the flow experiences a strong adverse pressure gradient along the rear surface, the boundary layer is unique in that intermittent detachment occurring near the wall. The numerical results indicate that the boundary layer is not followed by mean‐flow separation or incipient separation as shown from the numerical results. The resolved turbulent shear stress is in a reasonable agreement with the experimental data, though the computational result of LES shows that its peak is overpredicted near the trailing edge of the bump, while the other used turbulence models, except the standard k?ε, underpredicts it. Analysis of the numerical results from LES confirms the experimental data, in which the existence of internal layers over the bump surface upstream of the summit and along the downstream flat plate. It also demonstrates that the quasi‐step increase in skin friction is due to perturbations in pressure gradient. The surface curvature enhances the near‐wall shear production of turbulent stresses, and is responsible for the formation of the internal layers. The aim of the present work is to examine the response and prediction capability of LES with the dynamic eddy viscosity model as a sub‐grid scale to the complex turbulence structure with the presence of streamline curvature generated by a bumpy surface. Aiming to reduce the computational costs with focus on the mean behavior of the non‐equilibrium turbulent boundary layer of flow over the bump surface, the present investigation also explains the best capability of one of the used RANS turbulence models to capture the driving mechanism for the surprisingly rapid return to equilibrium over the trailing flat plate found in the measurements. Copyright © 2010 John Wiley & Sons, Ltd.     

On the construction of manufactured solutions for one and two‐equation eddy‐viscosity models

This paper presents manufactured solutions (MSs) for some well‐known eddy‐viscosity turbulence models, viz. the Spalart & Allmaras one‐equation model and the TNT and BSL versions of the two‐equation k–ω model. The manufactured flow solutions apply to two‐dimensional, steady, wall‐bounded, incompressible, turbulent flows. The two velocity components and the pressure are identical for all MSs, but various alternatives are considered for specifying the eddy‐viscosity and other turbulence quantities in the turbulence models. The results obtained for the proposed MSs with a second‐order accurate numerical method show that the MSs for turbulence quantities must be constructed carefully to avoid instabilities in the numerical solutions. This behaviour is model dependent: the performance of the Spalart & Allmaras and k–ω models is significantly affected by the type of MS. In one of the MSs tested, even the two versions of the k–ω model exhibit significant differences in the convergence properties. Copyright © 2006 John Wiley & Sons, Ltd.     

A non‐linear k–ε model with realizability for prediction of flows around bluff bodies

The incompressible flow around bluff bodies (a square cylinder and a cube) is investigated numerically using turbulence models. A non‐linear k–ε model, which can take into account the anisotropy of turbulence with less CPU time and computer memory then RSM or LES, is adopted as a turbulence model. In tuning of the model coefficients of the non‐linear terms are adjusted through the examination of previous experimental studies in simple shear flows. For the tuning of the coefficient in the eddy viscosity (=C_μ), the realizability constraints are derived in three types of basic 2D flow patterns, namely, a simple shear flow, flow around a saddle and a focal point. C_μ is then determined as a function of the strain and rotation parameters to satisfy the realizability. The turbulence model is first applied to a 2D flow around a square cylinder and the model performance for unsteady flows is examined focussing on the period and the amplitude of the flow oscillation induced by Karman vortex shedding. The applicability of the model to 3D flows is examined through the computation of the flow around a surface‐mounted cubic obstacle. The numerical results show that the present model performs satisfactorily to reproduce complex turbulent flows around bluff bodies. Copyright © 2003 John Wiley & Sons, Ltd.     

Time‐dependent turbulent cavitating flow computations with interfacial transport and filter‐based models

Turbulent cavitating flow computations need to address both cavitation and turbulence modelling issues. A recently developed interfacial dynamics‐based cavitation model (IDCM) incorporates the interfacial transport into the computational modelling of cavitation dynamics. For time‐dependent flows, it is known that the engineering turbulence closure such as the original k–ε model often over‐predicts the eddy viscosity values reducing the unsteadiness. A recently proposed filter‐based modification has shown that it can effectively modulate the eddy viscosity, rendering better simulation capabilities for time‐dependent flow computations in term of the unsteady characteristics. In the present study, the IDCM along with the filter‐based k–ε turbulence model is adopted to simulate 2‐D cavitating flows over the Clark‐Y airfoil. The chord Reynolds number is Re=7.0 × 10⁵. Two angles‐of‐attack of 5 and 8° associated with several cavitation numbers covering different flow regimes are conducted. The simulation results are assessed with the experimental data including lift, drag and velocity profiles. The interplay between cavitation and turbulence models reveals substantial differences in time‐dependent flow results even though the time‐averaged characteristics are similar. Copyright © 2005 John Wiley & Sons, Ltd.     

Low Reynolds number k–ε model for near‐wall flow

A wall‐distance free k–ε turbulence model is developed that accounts for the near‐wall and low Reynolds number effects emanating from the physical requirements. The model coefficients/functions depend non‐linearly on both the strain rate and vorticity invariants. Included diffusion terms and modified C_ε(1,2) coefficients amplify the level of dissipation in non‐equilibrium flow regions, thus reducing the kinetic energy and length scale magnitudes to improve prediction of adverse pressure gradient flows, involving flow separation and reattachment. The model is validated against a few flow cases, yielding predictions in good agreement with the direct numerical simulation (DNS) and experimental data. Copyright © 2005 John Wiley & Sons, Ltd.     

Time‐dependent simulation of cavitating flow with k − ℓ turbulence models

A compressible, multiphase, one‐fluid Reynolds‐averaged Navier–Stokes solver has been developed to study turbulent cavitating flows. The interplay between turbulence and cavitation regarding the unsteadiness and structure of the flow is complex and not well understood. This constitutes a critical point to accurately simulate the dynamic behavior of sheet cavities. In the present study, different formulations based on a k ? ? transport‐equation model are investigated and a scale‐adaptive formulation is proposed. Numerical results are given for a Venturi geometry and comparisons are made with experimental data. The scale‐adaptive model shows several improvements compared with standard turbulence models. Copyright © 2011 John Wiley & Sons, Ltd.     

Turbulent flow investigation inside and outside plain-orifice atomizers with rounded orifice inlets

The performances of three linear eddy viscosity models (LEVM) and one algebraic Reynolds stress model (ARSM) on the simulation of the internal and external flows in the plain-orifice atomizers with rounded orifice inlets are evaluated. The validity of the computational model is first assessed through the testing of a backward facing step flow, a sudden expansion pipe flow and a liquid column collapsing problem. Then the atomizer internal and external flows are analyzed by comparing the computed discharge coefficients with available experimental data and by comparing the turbulence intensity profiles at the orifice exit. The results are also illustrated by the fluid/air interface plot. It is found that the turbulence models investigated exhibit zonal behaviors, i.e., none of the models investigated performs well throughout the entire flow field. It is worthwhile to note that the standard k-model is not necessarily the worst among the models investigated. In average, the ARSM model gives better results as compared to the standard k-model and the low Reynolds number models. The turbulence strength has a significant influence on the global characteristics of the flow field. The models with better predictions of the turbulence kinetic energy, such as Gatski–Speziales ARSM model and Nagano–Hishidas low Reynolds number model, can yield better predictions of the global characteristics of the flow field, e.g., the reattachment lengths for the backward-facing step flow and the sudden expansion pipe flow, and the discharge coefficient for the atomizer flow.     

Finite element computation of a turbulent flow over a two-dimensional backward-facing step

The present paper is devoted to the computation of turbulent flows by a Galerkin finite element method. Effects of turbulence on the mean field are taken into account by means of a (k-ε) turbulence model. The wall region is treated through wall laws and, more specifically, Reichardt's law. An inlet profile for ε is proposed as a numerical treatment for physically meaningless values of k and ε. Results obtained for a recirculating flow in a two-dimensional channel with a sudden expansion in width are presented and compared with experimental values.     

Two‐dimensional prediction of time dependent,turbulent flow around a square cylinder confined in a channel

This paper presents two‐dimensional and unsteady RANS computations of time dependent, periodic, turbulent flow around a square block. Two turbulence models are used: the Launder–Sharma low‐Reynolds number k–ε model and a non‐linear extension sensitive to the anisotropy of turbulence. The Reynolds number based on the free stream velocity and obstacle side is Re=2.2×10⁴. The present numerical results have been obtained using a finite volume code that solves the governing equations in a vertical plane, located at the lateral mid‐point of the channel. The pressure field is obtained with the SIMPLE algorithm. A bounded version of the third‐order QUICK scheme is used for the convective terms. Comparisons of the numerical results with the experimental data indicate that a preliminary steady solution of the governing equations using the linear k–ε does not lead to correct flow field predictions in the wake region downstream of the square cylinder. Consequently, the time derivatives of dependent variables are included in the transport equations and are discretized using the second‐order Crank–Nicolson scheme. The unsteady computations using the linear and non‐linear k–ε models significantly improve the velocity field predictions. However, the linear k–ε shows a number of predictive deficiencies, even in unsteady flow computations, especially in the prediction of the turbulence field. The introduction of a non‐linear k–ε model brings the two‐dimensional unsteady predictions of the time‐averaged velocity and turbulence fields and also the predicted values of the global parameters such as the Strouhal number and the drag coefficient to close agreement with the data. Copyright © 2009 John Wiley & Sons, Ltd.     

An eddy viscosity model with near‐wall modifications

An extended version of the isotropic k–ε model is proposed that accounts for the distinct effects of low‐Reynolds number (LRN) and wall proximity. It incorporates a near‐wall correction term to amplify the level of dissipation in nonequilibrium flow regions, thus reducing the kinetic energy and length scale magnitudes to improve prediction of adverse pressure gradient flows, involving flow separation and reattachment. The eddy viscosity formulation maintains the positivity of normal Reynolds stresses and the Schwarz' inequality for turbulent shear stresses. The model coefficients/functions preserve the anisotropic characteristics of turbulence. The model is validated against a few flow cases, yielding predictions in good agreement with the direct numerical simulation (DNS) and experimental data. Comparisons indicate that the present model is a significant improvement over the standard eddy viscosity formulation. Copyright © 2005 John Wiley & Sons, Ltd.     

Flows through plane sudden‐expansions

A calculation method has been developed and used to represent flows downstream of plane symmetric expansions with dimensions and velocities encompassing laminar and turbulent flows. Except for very low Reynolds numbers, the flows are time‐dependent and asymmetric and the calculated results are appraised first in relation to published measurements of laminar flows and then to new measurements obtained at a Reynolds number of 26 500. The time‐dependent laminar simulations indicate that the critical Reynolds numbers are predicted with excellent accuracy for different expansion ratios and the details of the asymmetric velocity profiles are in good agreement with experimental measurements. The laminar flow calculations also show that increasing the thickness of the separating boundary layer leads to longer regions of separation and no dominant frequency for Reynolds numbers up to those at which the third separation region was observed. The turbulent flow simulations made use of the k–ε turbulence model and provided a satisfactory representation of measurements, except in regions close to the wall and within the recirculation regions. Also, the longer reattachment length was underestimated. Limitations are discussed in relation to these and higher‐order assumptions. Copyright © 2000 John Wiley & Sons, Ltd.     

Application of the two-equation k-ϵ turbulence model to a two-dimensional,steady, free surface flow problem with separation

An application of the depth-integrated k-? turbulence model is presented for separated flow in a wide, shallow, rectangular channel with an abrupt expansion in width. The well-known numerical problems associated with the use of upwind and central finite differences for convection are overcome by the adoption of the spatially third-order accurate QUICK finite difference technique. Results show that modification of the depth-integrated k-? turbulence closure model for streamline curvature leads to significant improvement in the agreement between model predictions and experimental measurements.     

A turbulent plane offset jet with small offset ratio

 Mean velocities and turbulence characteristics of a turbulent plane offset jet with a small offset ratio of 2.125 have been studied using laser Doppler anemometry (LDA). Static pressure measurements highlight the importance of side plates in enhancing two-dimensionality of the jet. The spatial distributions of turbulence intensities and Reynolds shear stress show a high turbulence recirculating flow region close to the nozzle plate between the jet and the offset plate. The LDA results have been used to examine the capability of three different turbulence models (i.e. k–ɛ, RNG and Reynolds stress) in predicting the velocity field of this jet. While all three models are able to predict qualitatively the recirculation, converging and reattachment regions observed experimentally, the standard k–ɛ turbulence model predicts a reattachment length that best agrees with the experimentally determined value. Received: 11 September 1996/Accepted: 30 May 1997