Numerical prediction of a turbulent curved wake and comparison with experimental data

Numerical studies of the curved wake of a NACA 0012 airfoil of chord length 0.150 m are presented. The airfoil is placed in air at 10 m/s in the straight section of a duct of 0.457 m × 0.457 m cross‐section followed by a 90° bend with a mean radius‐to‐height ratio of 1.17. The trailing edge is located at one chord length upstream of the bend entry plane. The authors' own measurements are used to define the boundary conditions and for comparison with the predicted results. The numerical models are based on the time‐averaged, three‐dimensional conservation equations of fluid flow, incorporating the k–ε, RNG k–ε, realizable k–ε and the Reynolds stress turbulence models. The results show that the models are capable of predicting the effects of curvature on the wake development. However, quantitative differences between prediction and experiment exist. The results obtained using the Reynolds stress model show better agreement with the experimental data, compared with the k–ε based models, but not consistently for all parameters. There are also better predictions by the RNG k–ε and realizable k–ε models compared with the standard k–ε model. The predicted results using the RNG k–ε are closer to experimental data than the realizable k–ε. Copyright © 2005 John Wiley & Sons, Ltd.     

Parameter estimation of engineering turbulence model

A parameter estimation algorithm is introduced and used to determine the parameters in the standardk-∈ two equation turbulence model (SKE). It can be found from the estimation results that although the parameter estimation method is an effective method to determine model parameters, it is difficult to obtain a set of parameters for SKE to suit all kinds of separated flow and a modification of the turbulence model structure should be considered. So, a new nonlineark-∈ two-equation model (NNKE) is put forward in this paper and the corresponding parameter estimation technique is applied to determine the model parameters. By implementing the NNKE to solve some engineering turbulent flows, it is shown that NNKE is more accurate and versatile than SKE. Thus, the success of NNKE implies that the parameter estimation technique may have a bright prospect in engineering turbulence model research.     

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.     

Evaluation of two RANS turbulence models in predicting developing mixed convection within a heated horizontal pipe

The developing weakly turbulent regime of mixed convection in a uniformly heated horizontal pipe was first studied experimentally, by means of heat transfer measurements in the following ranges of dimensionless numbers: 3.19 < Re × 10^? 3 < 6.39, 1.80 < Gr _h  × 10^? 8 < 4.20. The working fluid was FC-72?, with Pr = 12.4.

In order to gain a better insight into the thermo-fluid dynamics involved in the phenomenon and obtain the velocity and temperature fields at every point of the fluid domain, numerical simulations were performed by means of commercial software. Turbulence was modelled by using the Reynolds averaged Navier–Stokes equations (RANS) approach. Two closures of the governing equations were evaluated: realizable κ–? (RKE) model and renormalization-group κ–? (RNG) model.

Both models were capable of reproducing the observed physical trends. However, deviations from the experimental data lower than 20% were obtained only in the entry-zone with the RKE model, while the RNG model gave fair predictions only in developed or quasi-developed flow.     

OPERATOR–SPLITTING COMPUTATION OF TURBULENT FLOW IN AN AXISYMMETRIC 180° NARROWING BEND USING SEVERAL k—ϵ MODELS AND WALL FUNCTIONS

This paper describes a finite element implementation of an operator-splitting algorithm for solving transient/steady turbulent flows and presents solutions for the turbulent flow in an axisymmetric 180° narrowing bend, a benchmark problem dealt with at the 1994 WUA-CFD annual meeting. Three k–ϵ based models are used: the standard linear k–ϵ model, a non-linear k–ϵ model and an RNG k–ϵ model. Flow separation after the bend, as observed in the experiment, is predicted by the RNG model and by both the linear and non-linear k–ϵε models with van Driest mixing length wall functions. Good agreement with experimental data of pressure distribution on bending walls is obtained by the present numerical simulation. Results show that there is very little difference between the linear and non-linear k–ϵε models in terms of predicted velocity fields and that the non-linearities mainly affect the distribution of turbulent normal stress and pressure, in analogy to the effect of second-order viscoelastic fluid models on laminar flow. Both the linear and non-linear k–ϵε models fail to predict any flow separation if logarithmic wall functions are used.     

An immersed boundary/wall modeling method for RANS simulation of compressible turbulent flows

In this paper, an immersed boundary (IB) method is developed to simulate compressible turbulent flows governed by the Reynolds‐averaged Navier‐Stokes equations. The flow variables at the IB nodes (interior nodes in the immediate vicinity of the solid wall) are evaluated via linear interpolation in the normal direction to close the discrete form of the governing equations. An adaptive wall function and a 2‐layer wall model are introduced to reduce the near‐wall mesh density required by the high resolution of the turbulent boundary layers. The wall shear stress modified by the wall modeling technique and the no‐penetration condition are enforced to evaluate the velocity at an IB node. The pressure and temperature at an IB node are obtained via the local simplified momentum equation and the Crocco‐Busemann relation, respectively. The SST k ? ω and S‐A turbulence models are adopted in the framework of the present IB approach. For the Shear‐Stress Transport (SST) k ? ω model, analytical solutions in near‐wall region are utilized to enforce the boundary conditions of the turbulence equations and evaluate the turbulence variables at an IB node. For the S‐A model, the turbulence variable at an IB node is calculated by using the near‐wall profile of the eddy viscosity. In order to validate the present IB approach, numerical experiments for compressible turbulent flows over stationary and moving bodies have been performed. The predictions show good agreements with the referenced experimental data and numerical results.     

A simple sliding‐mesh interface procedure and its application to the CFD simulation of a tidal‐stream turbine

An effective way of using computational fluid dynamics (CFD) to simulate flow about a rotating device—for example, a wind or marine turbine—is to embed a rotating region of cells inside a larger, stationary domain, with a sliding interface between. This paper describes a simple but effective method for implementing this as an internal Dirichlet boundary condition, with interfacial values obtained by interpolation from halo nodes. The method is tested in two finite‐volume codes: one using block‐structured meshes and the other unstructured meshes. Validation is performed for flow around simple, isolated, rotating shapes (cylinder, sphere and cube), comparing, where possible, with experiment and the alternative CFD approach of fixed grid with moving walls. Flow variables are shown to vary smoothly across the sliding interface. Simulations of a tidal‐stream turbine, including both rotor and support, are then performed and compared with towing‐tank experiments. Comparison between CFD and experiment is made for thrust and power coefficients as a function of tip‐speed ratio (TSR) using Reynolds‐averaged Navier–Stokes turbulence models and large‐eddy simulation (LES). Performance of most models is good near the optimal TSR, but simulations underestimate mean thrust and power coefficients in off‐design conditions, with the standard k– ? turbulence model performing noticeably worse than shear stress transport k– ω and Reynolds‐stress‐transport closures. LES gave good predictions of mean load coefficients and vital information about wake structures but at substantial computational cost. Grid‐sensitivity studies suggest that Reynolds‐averaged Navier–Stokes models give acceptable predictions of mean power and thrust coefficients on a single device using a mesh of about 4 million cells. Copyright © 2013 John Wiley & Sons, Ltd.     

Evaluation of erosion in a two‐way coupled fluid–particle system

In this study, the effects of flow turbulence intensity, temperature, particle sizes and impinging velocity on erosion by particle impact are demonstrated numerically. Underlying turbulent flow on an Eulerian frame is described by the compressible Reynolds averaged Navier–Stokes equations with a RNG k–ε turbulence model. The particle trajectories and particle–wall interactions are evaluated by a Eulerian–Lagrangian approach in a two‐way coupling system. An erosion model considering material weight removal from surfaces is used to predict erosive wear. Computational validation against measured data is demonstrated satisfactorily. The analysis of erosion shows that the prevention of erosion is enhanced by increasing the effects of flow temperature and turbulence intensity and reducing particle inertial momentum. Copyright © 2001 John Wiley & Sons, Ltd.     

Verification testing in computational fluid dynamics: an example using Reynolds‐averaged Navier–Stokes methods for two‐dimensional flow in the near wake of a circular cylinder

Verification testing was performed for various Reynolds‐averaged Navier–Stokes methods for uniform flow past a circular cylinder at Re= 5232. The standard and renormalized group (RNG) versions of the k–ε method were examined, along with the Boussinesq, Speziale and Launder constitutive relationships. Wind tunnel experiments for flow past a circular cylinder were also performed to obtain a comparative data set. Preliminary studies demonstrate poor convergence for the Speziale relationship. Verification testing with the standard and RNG k–ε models suggests that the simulations exhibit global monotonic convergence for the Boussinesq models. However, the global order of accuracy of the methods was much lower than the expected order of accuracy of 2. For this reason, pointwise convergence ratios and orders of accuracy were computed to show that not all sampling locations had converged (standard k–ε model: 19% failed to converge; RNG k–ε model: 14% failed to converge). When the non‐convergent points were removed from consideration, the average orders of accuracy are closer to the expected value (standard k–ε model: 1.41; RNG k–ε model: 1.27). Poor iterative and global grid convergence was found for the RNG k–ε/Launder model. The standard and RNG k–ε models with the Boussinesq relationship were compared with experimental data and yielded results significantly different from the experiments. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

Suitability of the k–ω turbulence model for scramjet flowfield simulations

The suitability of Wilcox's 2006 k– ω turbulence model for scramjet flowfield simulations is demonstrated by validation against five test cases that have flowfields representative of those to be expected in scramjets. The five test cases include a 2D flat plate, an axisymmetric cylinder, a backward‐facing step, the mixing of a pair of coaxial jets and the interaction between a shock wave and turbulent boundary layer. A generally good agreement between the numerical and experimental results is obtained for all test cases. These tests reveal that despite the turbulence model's sensitivity to freestream turbulence properties, the numerically predicted skin friction agrees with experimental data and theoretical correlations to their degree of uncertainty. The tests also confirm the importance of using a y⁺ value of less than 1 in getting accurate surface heat transfer distributions. In the coaxial jets case, the importance of matching the turbulence intensities at the inflow plane in improving the predictions of the turbulent mixing phenomena is also shown. A review of guidelines with regard to the setting up of grids and specification of freestream turbulence properties for turbulent Reynolds‐averaged Navier–Stokes CFD simulations is also included in this paper. Copyright © 2011 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.     

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.     

Calculation of turbulent fluid flow and heat transfer in ducts by a full Reynolds stress model

A computational method has been developed to predict the turbulent Reynolds stresses and turbulent heat fluxes in ducts by different turbulence models. The turbulent Reynolds stresses and other turbulent flow quantities are predicted with a full Reynolds stress model (RSM). The turbulent heat fluxes are modelled by a SED concept, the GGDH and the WET methods. Two wall functions are used, one for the velocity field and one for the temperature field. All the models are implemented for an arbitrary three‐dimensional channel. Fully developed condition is achieved by imposing cyclic boundary conditions in the main flow direction. The numerical approach is based on the finite volume technique with a non‐staggered grid arrangement. The pressure–velocity coupling is handled by using the SIMPLEC‐algorithm. The convective terms are treated by the van Leer scheme while the diffusive terms are handled by the central‐difference scheme. The hybrid scheme is used for solving the ε equation. The secondary flow generation using the RSM model is compared with a non‐linear k–ε model (non‐linear eddy viscosity model). The overall comparison between the models is presented in terms of the friction factor and Nusselt number. Copyright © 2003 John Wiley & Sons, Ltd.     

Least‐squares finite element processes in h,p, k mathematical and computational framework for a non‐linear conservation law

This paper considers numerical simulation of time‐dependent non‐linear partial differential equation resulting from a single non‐linear conservation law in h, p, k mathematical and computational framework in which k=(k₁, k₂) are the orders of the approximation spaces in space and time yielding global differentiability of orders (k₁?1) and (k₂?1) in space and time (hence k‐version of finite element method) using space–time marching process. Time‐dependent viscous Burgers equation is used as a specific model problem that has physical mechanism for viscous dissipation and its theoretical solutions are analytic. The inviscid form, on the other hand, assumes zero viscosity and as a consequence its solutions are non‐analytic as well as non‐unique (Russ. Math. Surv. 1962; 17 (3):145–146; Russ. Math. Surv. 1960; 15 (6):53–111). In references (Russ. Math. Surv. 1962; 17 (3):145–146; Russ. Math. Surv. 1960; 15 (6):53–111) authors demonstrated that the solutions of inviscid Burgers equations can only be approached within a limiting process in which viscosity approaches zero. Many approaches based on artificial viscosity have been published to accomplish this including more recent work on H(Div) least‐squares approach (Commun. Pure Appl. Math. 1965; 18 :697–715) in which artificial viscosity is a function of spatial discretization, which diminishes with progressively refined discretizations. The thrust of the present work is to point out that: (1) viscous form of the Burgers equation already has the essential mechanism of viscosity (which is physical), (2) with progressively increasing Reynolds (Re) number (thereby progressively reduced viscosity) the solutions approach that of the inviscid form, (3) it is possible to compute numerical solutions for any Re number (finite) within hpk framework and space–time least‐squares processes, (4) the space–time residual functional converges monotonically and that it is possible to achieve the desired accuracy, (5) space–time, time marching processes utilizing a single space–time strip are computationally efficient. It is shown that viscous form of the Burgers equation without linearizing provides a physical and viablemechanism for approaching the solutions of inviscid form with progressively increasing Re. Numerical studies are presented and the computed solutions are compared with published work. Copyright © 2008 John Wiley & Sons, Ltd.     

Calculation of three‐dimensional turbulent flow with a finite volume multigrid method

A numerical method for the efficient calculation of three‐dimensional incompressible turbulent flow in curvilinear co‐ordinates is presented. The mathematical model consists of the Reynolds averaged Navier–Stokes equations and the k–ε turbulence model. The numerical method is based on the SIMPLE pressure‐correction algorithm with finite volume discretization in curvilinear co‐ordinates. To accelerate the convergence of the solution method a full approximation scheme‐full multigrid (FAS‐FMG) method is utilized. The solution of the k–ε transport equations is embedded in the multigrid iteration. The improved convergence characteristic of the multigrid method is demonstrated by means of several calculations of three‐dimensional flow cases. Copyright © 1999 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.     

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.     

Impact of turbulence modelling on external supersonic flow field simulations in rocket aerodynamics

This paper contains a three-dimensional study on the influence of different turbulence models for external supersonic flow field simulations, aiming at reaching the best accuracy in rocket aerodynamics. A well-studied test case -a slender body- has been used for the validation process, which involves the major turbulence models available. The SST k–ω model has been selected as the most suitable one for this kind of flows. Good agreements between numerical, theoretical and experimental results are obtained, which are used to set up some guidelines regarding the configuration of Reynolds-averaged Navier–Stokes Computational Fluid Dynamics (CFD) supersonic models for these flight regimes.     

Experimental validation of two depth‐averaged turbulence models

A finite volume turbulence model for the resolution of the two‐dimensional shallow water equations with turbulent term is presented. After making a finite volume discretization of the depth‐averaged k–ε equations in conservative form, the q–r equations, that give stability to the process, are obtained. Wall and inlet boundary conditions for the turbulent equations and wall conditions for the hydrodynamic equations are discussed. A comparison between the k–ε and q–r models and some experimental results is made. Copyright © 2008 John Wiley & Sons, Ltd.