Analysis and modelling of the relation between the shear rate and Reynolds stress tensors in transitional boundary layers

A non-linear eddy-viscosity transition model is presented, tuned by a large experimental data set describing transitional boundary layers. Data have been acquired by TR-PIV on a flat plate placed in a 2D converging-diverging channel with variable opening angle, allowing variation of the adverse pressure gradient, the free-stream turbulence intensity and the flow Reynolds number. Overall, 48 different combinations of these flow parameters encompass different modes of transition from bypass to separated-flow mechanisms, thus allowing fine tuning of the model, spanning significantly different conditions. The model is tuned locally as a function of the turbulent kinetic energy, a Reynolds number based on the wall distance and the ℓ₂-norm of the shear rate tensor. A first correlation determines the rotation for alignment of the principal axes of the shear and stress tensors. By a second correlation, the eigenvalues of the stress tensor are obtained. The non-linear eddy-viscosity relation reproduces the anisotropy of the turbulence field observed for both bypass and separated-flow transitional cases. The relation has been applied to another experimental data set that did not participate to the fitting of the model and that is characterized by a different range of Reynolds number and turbulence intensity and a significantly stronger adverse pressure gradient with respect to the tuning dataset. Such application further strengthens the capability of the proposed correlations, that can easily be implemented in existing CFD solvers.     

Turbulence properties in the cylinder wake at high Reynolds numbers

The present contribution analyses the turbulence properties in unsteady flows around bluff body wakes and provides a database for improvement and validation of turbulence models, concerning the present class of nonequilibrium flows. The flow around a circular cylinder with a low aspect ratio and a high blockage coefficient is investigated. This confined environment is used in order to allow direct comparisons with realisable 3-D Navier–Stokes computations avoiding ‘infinite’ conditions. The flow is investigated in the beginning of the critical regime at Reynolds number 140 000. The analysis is carried out by means of 2-D PIV, of 3-C PIV and of high-frequency 2-D PIV. The experimental analysis contributes to confirm the validity of advanced statistical turbulence modelling for unsteady flows around bodies.     

A dynamic global-coefficient mixed subgrid-scale model for large-eddy simulation of turbulent flows

A dynamic global-coefficient mixed subgrid-scale eddy-viscosity model for large-eddy simulation of turbulent flows in complex geometries is developed. In the present model, the subgrid-scale stress is decomposed into the modified Leonard stress, cross stress, and subgrid-scale Reynolds stress. The modified Leonard stress is explicitly computed assuming a scale similarity, while the cross stress and the subgrid-scale Reynolds stress are modeled using the global-coefficient eddy-viscosity model. The model coefficient is determined by a dynamic procedure based on the global-equilibrium between the subgrid-scale dissipation and the viscous dissipation. The new model relieves some of the difficulties associated with an eddy-viscosity closure, such as the nonalignment of the principal axes of the subgrid-scale stress tensor and the strain rate tensor and the anisotropy of turbulent flow fields, while, like other dynamic global-coefficient models, it does not require averaging or clipping of the model coefficient for numerical stabilization. The combination of the global-coefficient eddy-viscosity model and a scale-similarity model is demonstrated to produce improved predictions in a number of turbulent flow simulations.     

A dynamic subgrid-scale tensorial Eddy viscosity model

In the Navier-Stokes equations the removal of the turbulent fluctuating velocities with a frequency above a certain fixed threshold, employed in the Large Eddy Simulation (LES), causes the appearance of a turbulent stress tensor that requires a number of closure assumptions. In this paper insufficiencies are demonstrated for those closure models which are based on a scalar eddy viscosity coefficient. A new model, based on a tensorial eddy viscosity, is therefore proposed; it employs the Germano identity [1] and allows dynamical evaluation of the single required input coefficient. The tensorial expression for the eddy viscosity is deduced by removing the widely used scalar assumption of the high-frequency viscous dissipation and replacing it by its tensorial counterpart arising in the balance of the Reynolds stress tensor. The numerical simulations performed for a lid driven cavity flow show that the proposed model allows to overcome the drawbacks encountered by the scalar eddy viscosity models. Received November 25, 1997     

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.     

Three-dimensional LDA measurements in the impeller region of a turbulently stirred tank

Three-dimensional, angle-resolved LDA measurements of the turbulent flow field (Re = 2.9 × 10⁴) in the vicinity of a Rushton turbine in a baffled mixing tank have been performed. For this goal, a procedure for accurate beam alignment, based on a submerged micro-mirror system, has been developed. Results on the average flow field as well as on the complete set of Reynolds stresses are presented. The anisotropy of the turbulence has been characterized by the invariants of the anisotropy tensor. The trailing vortex structure, which is characteristic for the flow induced by a Rushton turbine, is demonstrated to be associated with strong, anisotropic turbulent activity. Received: 14 September 1998/Accepted: 22 February 1999     

Large eddy simulations on the effect of the irregular roughness shape on turbulent channel flows

Large Eddy Simulations (LES) are carried out to investigate on the mean flow in turbulent channel flows over irregular rough surfaces. Here the attention is focused to selectively investigate on the effect induced by crests or cavities of the roughness. The irregular shape is generated through the super-imposition of sinusoidal functions having random amplitude and four different wave-lengths. The irregular roughness profile is reproduced along the spanwise direction in order to obtain a 2D rough shape. The analysis of the mean velocity profiles shows that roughness crests induce higher effect in the outer-region whereas roughness cavities cause the highest effects in the inner-region with a reduced effect in the external region. The numerical simulations have been carried out at friction Reynolds number Re_τ=395. Similar results have been found for the higher order statistics: turbulence intensities or shear stresses. The analysis of the Reynolds stress anisotropy tensor confirms the existence of specific roles of cavities and crests in the turbulence modulation.     

Turbulence modelling of the flow past a pitching NACA0012 airfoil at and Reynolds numbers

This paper provides a study of the NACA0012 dynamic stall at Reynolds numbers 10⁵ and 10⁶ by means of two- and three-dimensional numerical simulations. The turbulence effect on the dynamic stall is studied by statistical modelling. The results are compared with experiments concerning each test case. Standard URANS turbulence modelling have shown a quite dissipative character that attenuates the instabilities and the vortex structures related to the dynamic stall. The URANS approach Organised Eddy Simulation (OES) has shown an improved behaviour at the high Reynolds number range. Emphasis is given to the physical analysis of the three-dimensional dynamic stall structure, for which there exist few numerical results in the literature, as far as the Reynolds number range is concerned. This study has shown that the downstroke phases of the pitching motion are subjected to strong three-dimensional turbulence effects along the span, whereas the flow is practically two-dimensional during the upstroke motion.     

Computation of high Reynolds number flow around a circular cylinder with surface roughness

Incompressible high-Reynolds-number flows around a circular cylinder are analyzed by direct integration of the Navier-Stokes equations using finite-difference method. A generalized coordinate system is used so that a sufficient number of grid points are distributed in the boundary layer and the wake. A numerical scheme which suppresses non-linear instability for calculations of high-Reynolds-number flows is developed. The computation of an impulsively started flow at Re = 1200 is compared with corresponding experimental observations, and excellent agreements are obtained.A series of computations are carried out on the flow around a circular cylinder with surface roughness. The height of the roughness in these computations is 0.5% of the diameter. The range of Reynolds numbers is from 10³ to 10⁵; no turbulence model is employed. Sharp reduction of drag coefficient is observed near Re = 2 × 10⁴, which indicates that the critical Reynolds number is captured in the present computation.     

From single-scale turbulence models to multiple-scale and subgrid-scale models by Fourier transform

A theoretical method based on mathematical physics formalism that allows transposition of turbulence modeling methods from URANS (unsteady Reynolds averaged Navier–Stokes) models, to multiple-scale models and large eddy simulations (LES) is presented. The method is based on the spectral Fourier transform of the dynamic equation of the two-point fluctuating velocity correlations with an extension to the case of non-homogenous turbulence. The resulting equation describes the evolution of the spectral velocity correlation tensor in wave vector space. Then, we show that the full wave number integration of the spectral equation allows one to recover usual one-point statistical closure whereas the partial integration based on spectrum splitting gives rise to partial integrated transport models (PITM). This latter approach, depending on the type of spectral partitioning used, can yield either a statistical multiple-scale model or subfilter transport models used in LES or hybrid methods, providing some appropriate approximations are made. Closure hypotheses underlying these models are then discussed by reference to physical considerations with emphasis on identification of tensorial fluxes that represent turbulent energy transfer or dissipation. Some experiments such as the homogeneous axisymmetric contraction, the decay of isotropic turbulence, the pulsed turbulent channel flow and a wall injection induced flow are then considered as typical possible applications for illustrating the potentials of these models.      

Anisotropy measurements in the boundary layer over a flat plate with suction

Laser Doppler velocity measurements are carried out in a turbulent boundary layer subjected to concentrated wall suction (through a porous strip). The measurements are taken over a longitudinal distance of 9× the incoming boundary layer thickness ahead of the suction strip. The mean and rms velocity profiles are affected substantially by suction. Two-point measurements show that the streamwise and wall-normal autocorrelations of the streamwise velocity are reduced by suction. It is found that suction alters the redistribution of the turbulent kinetic energy k between its components. Relative to the no-suction case, the longitudinal Reynolds stress contributes more to k than the other two normal Reynolds stresses; in the outer region, its contribution is reduced which suggests structural changes in the boundary layer. This is observed in the anisotropy of the Reynolds stresses, which depart from the non-disturbed boundary layer. With suction, the anisotropy level in the near-wall region appears to be stronger than that of the undisturbed layer. It is argued that the mean shear induced by suction on the flow is responsible for the alteration of the anisotropy. The variation of the anisotropy of the layer will make the development of a turbulence model quite difficult for the flow behind suction. In that respect, a turbulence model will need to reproduce well the effects of suction on the boundary layer, if the model is to capture the effect of suction on the anisotropy of the Reynolds stresses.     

Applications of URANS on predicting unsteady turbulent separated flows

Accurate prediction of unsteady separated turbulent flows remains one of the toughest tasks and a practi cal challenge for turbulence modeling. In this paper, a 2D flow past a circular cylinder at Reynolds number 3,900 is numerically investigated by using the technique of unsteady RANS （URANS）. Some typical linear and nonlinear eddy viscosity turbulence models （LEVM and NLEVM） and a quadratic explicit algebraic stress model （EASM） are evaluated. Numerical results have shown that a high-performance cubic NLEVM, such as CLS, are superior to the others in simulating turbulent separated flows with unsteady vortex shedding.     

Invariant analysis of the Reynolds stress tensor for a nuclear fuel assembly with spacer grid and split type vanes

Invariant analysis of the Reynolds stress tensor anisotropy can give an accurate and deep intuitive understanding of the turbulent structure of a turbulent flow. Lumley's triangle has proven to be a powerful representation of the invariant analysis of the second-order statistics collection provided by the Reynolds stress tensor. In the present work the spectral element code Nek5000 has been used to investigate the turbulent structure of the flow across a pressurized water reactor spacer grid with split type mixing vanes. Wall-resolved large eddy simulation of the flow in a prototypical rod bundle geometry at Re = 14,000 and P/D = 1.32 are performed and validated against particle image velocimetry data. The results are then used to perform an in-depth invariant analysis. The results show a reorganization of the Reynolds stresses components in the downstream region of the spacer grid. The mixing vanes orientation produces a symmetric behavior between sub-channels. The turbulent structure in the fully developed region has the typical behavior of fully-developed channel flow turbulence. When averaging the state across regions of the sub-channels, we observed a transition from disk-like turbulence in the mixing vanes region to rod-like turbulence in the fully developed region.     

Simulation of cross-flow-induced vibration of cylinder arrays by surface vorticity method

The surface vorticity method (SVM), which is a fast and practical grid-free two-dimensional (2-D) method, and a fluid–structure interaction model incorporating the effects of cylinder motions and displacements is used to simulate the vortex-induced vibration of cylinder arrays at sub-critical Reynolds number Re=2.67×10⁴. The SVM is found to be most suitable for simulating a 2-D cylinder row with large-amplitude vibrations where the vorticity field and the fluid forces of the cylinder row change drastically, and the effect of the stream on the transverse direction vibration is very significant. The fluidelastic instability of a flexible cylinder row at small pitch ratio is also investigated, and the critical reduced velocity of the cylinder row at a reduced damping parameter S_G=1.29 is calculated, which is in good agreement with experimental and analytical results of the unsteady model. Vortex-induced vibration of a staggered cylinder array is simulated using different structural parameters. When the cylinders are relatively more flexible, the flow pattern changes dramatically and the fluid–structure interaction has a dominant impact on the flow field. Compared with grid-based methods, the grid-free SVM is a fast and practical method for the simulation of the FIV of cylinder arrays due to vortex shedding at sub-critical Reynolds numbers.     

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.     

Simulation of turbulent flows around a circular cylinder using nonlinear eddy-viscosity modelling: steady and oscillatory ambient flows

Steady incident flow past a circular cylinder for sub- to supercritical Reynolds number has been simulated as an unsteady Reynolds-averaged Navier–Stokes (RANS) equation problem using nonlinear eddy-viscosity modelling assuming two-dimensional flow. The model of Craft et al. (Int. J. Heat Fluid Flow 17 (1996) 108), with adjustment of the coefficients of the ‘cubic’ terms, predicts the drag crisis at a Reynolds number of about 2×10⁵ due to the onset of turbulence upstream of separation and associated changes in Strouhal number and separation positions. Slightly above this value, at critical Reynolds numbers, drag is overestimated because attached separation bubbles are not simulated. These do not occur at supercritical Reynolds numbers and drag coefficient, Strouhal number and separation positions are in approximate agreement with experimental measurements (which show considerable scatter). Fluctuating lift predictions are similar to sectional values measured experimentally for subcritical Reynolds numbers but corresponding measurements have not been made at supercritical Reynolds numbers. For oscillatory ambient flow, in-line forces, as defined by drag and inertia coefficients, have been compared with the experimental values of Sarpkaya (J. Fluid Mech. 165 (1986) 61) for values of the frequency parameter, β=D²/νT, equal to 1035 and 11240 and Keulegan–Carpenter numbers, KC=U₀T/D, between 0.2 and 15 (D is cylinder diameter, ν is kinematic viscosity, T is oscillation period, and U₀ is the amplitude of oscillating velocity). Variations with KC are qualitatively reproduced and magnitudes show best agreement when there is separation with a large-scale wake, for which the turbulence model is intended. Lift coefficients, frequency and transverse vortex shedding patterns for β=1035 are consistent with available experimental information for β≈250−500. For β=11240, it is predicted that separation is delayed due to more prominent turbulence effects, reducing drag and lift coefficients and causing the wake to be more in line with the flow direction than transverse to it. While these oscillatory flows are highly complex, attached separation bubbles are unlikely and the flows probably two dimensional.     

Comparison of flow structures in the downstream region of a cylinder and sphere

An experimental investigation of flow structures downstream of a circular cylinder and sphere immersed in a free-stream flow is performed for Re = 5000 and 10,000 using qualitative and quantitative flow visualization techniques. The obtained results are presented in terms of time-averaged velocity vectors, patterns of streamlines, vorticity, Reynolds stress correlations and turbulent kinetic energy distributions. Flow data reveal that the size of wake flow region, the location of singular and double points, the peak values of turbulence quantities, such as Reynolds stress correlations, vorticity fluctuations and turbulent kinetic energy vary as a function of models’ geometry and Reynolds Numbers. The concentration of small scale vortices is more dominant in the wake of the sphere than that of the cylinder. The maximum value of turbulent kinetic energy (TKE) occurs close to the saddle point for the cylinder case while two maximum values of TKE occur along shear layers for the sphere one because of the 3-D flow behavior.     

Comparison of flow structures in the downstream region of a cylinder and sphere

An experimental investigation of flow structures downstream of a circular cylinder and sphere immersed in a free-stream flow is performed for Re = 5000 and 10,000 using qualitative and quantitative flow visualization techniques. The obtained results are presented in terms of time-averaged velocity vectors, patterns of streamlines, vorticity, Reynolds stress correlations and turbulent kinetic energy distributions. Flow data reveal that the size of wake flow region, the location of singular and double points, the peak values of turbulence quantities, such as Reynolds stress correlations, vorticity fluctuations and turbulent kinetic energy vary as a function of models’ geometry and Reynolds Numbers. The concentration of small scale vortices is more dominant in the wake of the sphere than that of the cylinder. The maximum value of turbulent kinetic energy (TKE) occurs close to the saddle point for the cylinder case while two maximum values of TKE occur along shear layers for the sphere one because of the 3-D flow behavior.     

Balance equation of the generalised sub-grid scale (SGS) turbulent kinetic energy in a new tensorial dynamic mixed SGS model

A new dynamic model is proposed in which the eddy viscosity is defined as a symmetric second rank tensor, proportional to the product of a turbulent length scale with an ellipsoid of turbulent velocity scales. The employed definition of the eddy viscosity allows to remove the local balance assumption of the SGS turbulent kinetic energy formulated in all the dynamic Smagorinsky-type SGS models. Furthermore, because of the tensorial structure of the eddy viscosity the alignment assumption between the principal axes of the SGS turbulent stress tensor and the resolved strain-rate tensor is equally removed, an assumption which is employed in the scalar eddy viscosity SGS models. The proposed model is tested for a turbulent channel flow. Comparison with the results obtained with other dynamic SGS models (Dynamic Smagorinsky Model, Dynamic Mixed Model and Dynamic K-equation Model) shows that the tensorial definition of the eddy viscosity and the removal of the local balance assumption of the SGS turbulent kinetic energy considerably improves the agreement between results obtained with Large Eddy simulation (LES) and Direct Numerical Simulations (DNS), respectevely. Received August 26, 1999     

Performance assessment of a non-linear eddy-viscosity turbulence model applied to the anisotropic wake flow of a low-pressure turbine blade

The wake flow produced by a low-pressure turbine blade is modeled using a non-linear eddy-viscosity turbulence model. The theoretical benefit of using a non-linear eddy-viscosity model is strongly related to the capability of resolving highly anisotropic flows in contrast to the linear turbulence models, which are unable to correctly predict anisotropy. The main aim of the present work is to practically assess the performance of the model, by examining its ability to capture the anisotropic behavior of the wake-flow, mainly focusing on the measured velocity and Reynolds-stress distributions and to provide accurate results for the turbulent kinetic energy balance terms. Additionally, the contribution of each term of its non-linear constitutive expression for the Reynolds stresses is also investigated, in order to examine their direct effect on the modeling of the wake flow. The assessment is based on the experimental measurements that have been carried-out by the same group in Thessaloniki, Sideridis et al. (2011). The computational results show that the non-linear eddy viscosity model is capable to predict, with a good accuracy, all the flow and turbulence parameters while it is easy to program it in a computer code thus meeting the expectations of its originators.