Numerical study of gas-cyclone airflow: an investigation of turbulence modelling approaches

A numerical study of unsteady single-phase vortical flow inside a cyclone is presented. Two different geometric configurations have been considered, with the goal of assessing several different turbulence modelling approaches for this class of problem. The models investigated include three Reynolds-averaged Navier–Stokes models: a commonly used two-equation eddy-viscosity model, a differential Reynolds stress model (DRSM) and an eddy-viscosity model sensitised to rotational and curvature (RC) effects which was recently developed and implemented into a commercial CFD (computational fluid dynamics) code by the authors. Results were also obtained using large eddy simulation (LES). The computational results are analysed and compared with available experimental data. The RC-sensitised eddy-viscosity model shows significant improvement over the standard eddy-viscosity model. The RC-sensitised model, DRSM and LES model predictions of the mean flowfield are in good agreement with the experimental data. The results suggest that curvature- and rotation-sensitive eddy-viscosity models may provide a practical alternative to more computationally intensive approaches.     

Reynolds averaged simulation of flow and heat transfer in ribbed ducts

The accuracy of modern eddy-viscosity type turbulence models in predicting turbulent flows and heat transfer in complex passages is investigated. The particular geometries of interest here are those related to turbine blade cooling systems. This paper presents numerical data from the calculation of the turbulent flow field and heat transfer in two-dimensional (2D) cavities and three-dimensional (3D) ribbed ducts. It is found that heat transfer predictions obtained using the v²–f turbulence model for the 2D cavity are in good agreement with experimental data. However, there is only fair agreement with experimental data for the 3D ribbed duct. On the wall of the duct where ribs exist, predicted heat transfer agrees well with experimental data for all configurations (different streamwise rib spacing and the cavity depth) considered in this paper. But heat transfer predictions on the smooth-side wall do not concur with the experimental data. Evidence is provided that this is mainly due to the presence of strong secondary flow structures which might not be properly simulated with turbulence models based on eddy viscosity.     

Predicting Noninertial Effects with Linear and Nonlinear Eddy-Viscosity, and Algebraic Stress Models

Three types of turbulence models which account for rotational effects in noninertial frames of reference are evaluated for the case of incompressible, fully developed rotating turbulent channel flow. The different types of models are a Coroiolis-modified eddy-viscosity model, a realizable nonlinear eddy-viscosity model, and an algebraic stress model which accounts for dissipation rate anisotropies. A direct numerical simulation of a rotating channel flow is used for the validation of the turbulence models. This simulation differs from previous studies in that significantly higher rotation numbers are investigated. Flows at these higher rotation numbers are characterized by a relaminarization on the cyclonic or suction side of the channel, and a linear velocity profile on the anticyclonic or pressure side of the channel. The predictive performance of the three types of models are examined in detail, and formulation deficiencies are identified which cause poor predictive performance for some of the models. Criteria are identified which allow for accurate prediction of such flows by algebraic stress models and their corresponding Reynolds stress formulations.     

On stagnation pressure increases in calorically perfect,ideal gases

When stagnation pressure rises in a natural or numerically simulated flow it is frequently a cause for concern, as one usually expects viscosity and turbulence to cause stagnation pressure to decrease. In fact, if stagnation pressure increases, one may suspect measurement or numerical errors. However, this need not be the case, as the laws of nature do not require that stagnation pressure continually decreases. In order to help clarify matters, the objective of this work is to understand the conditions under which stagnation pressure will rise in the unsteady/steady flows of compressible, viscous, calorically perfect, ideal gases. Furthermore, at a more practical level, the goal is to understand the conditions under which stagnation pressure will increase in flows simulated with the Reynolds averaged Navier–Stokes equations and eddy-viscosity turbulence models. In order to provide an improved understanding of increases in stagnation pressure for both these scenarios, transport equations are derived that govern its behavior in the unaveraged and Reynolds averaged settings. These equations are utilized to precisely determine the relationship between changes in stagnation pressure and zeroth, first, and second derivatives of fundamental flow quantities. Furthermore, these equations are utilized to demonstrate the relationship between changes in stagnation pressure and fundamental non-dimensional quantities that govern the conductivity, viscosity, and compressibility of the flow. In addition, based on an analysis of the Reynolds averaged equation (for eddy-viscosity turbulence models), it is shown that stagnation pressure is particularly likely to experience a spurious rise at the outer edges of shear layers that are undergoing convex curvature. Thereafter, numerical experiments are performed which confirm the primary aspects of the theoretical analysis.     

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.     

Modelling turbulent separated flow in the context of aerodynamic applications

In contrast to rapid advances in computing, numerical methods and visualisation, the predictive capabilities of statistical models of turbulence are limited and improve only slowly, despite much intensive research in the recent past. The intuitive nature of turbulence modelling, its strong reliance on calibration and validation, the extreme sensitivity of model performance to seemingly minor variations in modelling details and flow conditions, and the fact that the non-local dynamics of turbulence are not well captured by single-point closure, all conspire to make turbulence modelling an especially demanding component of CFD, but one that is crucially important for the correct prediction of complex flows. This applies in particular to separation from streamlined bodies, which is, from a computational point of view, the most challenging flow feature in aeronautical CFD.This paper reviews some aspects of the foundation and application of turbulence models to flows that relate to aeronautical practice, with particular emphasis being placed on turbulence-transport models at a closure level higher than that based on the Boussinesq-viscosity hypothesis. Following a review of basic modelling issues, including aspects of linear-eddy-viscosity two-equation modelling, some recent experience and current work on predicting separation from continuous surfaces with non-linear eddy-viscosity models and second-moment closure are reported. The predictive performance of several anisotropy-resolving models is illustrated by reference to computational solutions for a number of flows, both two- and three-dimensional, some compressible and others incompressible.     

内锥流量计流出系数预测方法研究

采用标准k-ε模型、RNG(Renormalization Group)k-ε模型、Realizable k-ε模型和Reynolds应力方程模型 RSM(Reynolds Stress Model) 对100 mm口径6种结构的内锥流量计内流场进行了数值模拟.在等效直径比β值为0.65的三种结构内锥流量计流出系数的仿真计算中,四种湍流模型计算结果与物理实验结果误差的平均值分别为4.19%,2.84%,2.88%和-0.822%;对β值为0.85的情况,各模型计算误差的平均值分别为11.8%,9.62%,9.30%和4.76%.研究结果表明,RSM模型在6种结构内锥流量计流出系数的预测中,计算精度较高,表现出了较好的性能,优于三种k-ε涡粘模型,更适于内锥流量计流场数值模拟与流出系数的预测.     

Development of Gas-Particle Euler-Euler LES Approach: A Priori Analysis of Particle Sub-Grid Models in Homogeneous Isotropic Turbulence

A new large eddy simulation (LES) approach for particle-laden turbulent flows in the framework of the Eulerian formalism for inertial particle statistical modelling is developed. Local instantaneous Eulerian equations for the particle cloud are first written using the mesoscopic Eulerian formalism (MEF) proposed by Février et al. (J Fluid Mech 533:1–46, 2005), which accounts for the contribution of an uncorrelated velocity component for inertial particles with relaxation time larger than the Kolmogorov time scale. Second, particle LES equations are obtained by volume filtering the mesoscopic Eulerian ones. In such an approach, the particulate flow at larger scales than the filter width is recovered while sub-grid effects need to be modelled. Particle eddy-viscosity, scale similarity and mixed sub-grid stress (SGS) models derived from fluid compressible turbulence SGS models are presented. Evaluation of such models is performed using three sets of particle Lagrangian results computed from discrete particle simulation (DPS) coupled with fluid direct numerical simulation (DNS) of homogeneous isotropic decaying turbulence. The two phase flow regime corresponds to the dilute one where two-way coupling and inter-particle collisions are not considered. The different particle Stokes number (based on Kolmogorov time scale) are initially equal to 1, 2.2 and 5.1. The mesoscopic field properties are analysed in detail by considering the particle velocity probability function (PDF), correlated velocity power spectra and random uncorrelated velocity moments. The mesoscopic fields measured from DPS+DNS are then filtered to obtain large scale fields. A priori evaluation of particle sub-grid stress models gives comparable agreement than for fluid compressible turbulence models. It has been found that the standard Smagorinsky eddy-viscosity model exhibits the smaller correlation coefficients, the scale similarity model shows very good correlation coefficient but strongly underestimates the sub-grid dissipation and the mixed model is on the whole superior to pure eddy-viscosity model.     

高阶矩湍流模型研究进展及挑战

高阶矩模型是湍流模式理论研究中的难点和前沿.自周培源先生首次建立一般湍流的雷诺应力输运方程起,为了更精确的预测复杂流动,人们从未间断过对高阶矩模型的研究.尤其进入新世纪以来,随着计算机硬件水平的飞跃和高精度数值算法的突破,湍流模拟方法正由RANS向LES转变.而无论对于RANS框架、LES框架还是两者混合,高阶矩模式都...     

2D numerical flow modeling in a macro‐rough channel

A 2D numerical flow model, developed at the Research unit of Hydrology, Applied Hydrodynamics and Hydraulic Constructions at ULg, has been applied to flows in a macro‐rough channel. The model solves the shallow water equations (SWE) with a two length scale, depth‐integrated k‐type approach for turbulence modeling. Data for the comparison have been provided by experiments conducted at the Laboratory of Hydraulic Constructions at EPFL. In the experiments with different non‐prismatic channel configurations, namely large‐scale cavities at the side walls, three different 2D flow characteristics could be observed in cavities. With the used numerical model features, especially regarding turbulence and friction modeling, a single set of bottom and side wall roughness could be found for a large range of discharges investigated in a prismatic channel. For the macro rough configurations, the numerical model gives an excellent agreement between experimental and numerical results regarding backwater curves and flow patterns if the side wall cavities have low aspect ratios. For configurations with high aspect ratios, the head loss generated by the preservation of important recirculation gyres in the cavities is slightly underestimated. The results of the computations reveal clearly that the separation of turbulence sources in the mathematical model is of great importance. Indeed, the turbulence related to 2D transverse shear effects and the 3D turbulence, generated by bed friction, can have very different amplitude. When separating these two effects in the numerical models, most of the flow features observed experimentally can be reproduced accurately. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

Introduction of turbulent model in a mixed finite volume/finite element method

The aim of this work is to present a new numerical method to compute turbulent flows in complex configurations. With this in view, a k-? model with wall functions has been introduced in a mixed finite volume/finite element method. The numerical method has been developed to deal with compressible flows but is also able to compute nearly incompressible flows. The physical model and the numerical method are first described, then validation results for an incompressible flow over a backward-facing step and for a supersonic flow over a compression ramp are presented. Comparisons are performed with experimental data and with other numerical results. These simulations show the ability of the present method to predict turbulent flows, and this method will be applied to simulate complex industrial flows (flow inside the combustion chamber of gas turbine engines). The main goal of this paper is not to test turbulence models, but to show that this numerical method is a solid base to introduce more sophisticated turbulence model.     

Separation-induced boundary layer transition: Modeling with a non-linear eddy-viscosity model coupled with the laminar kinetic energy equation

We present an effort to model the separation-induced transition on a flat plate with a semi-circular leading edge, using a cubic non-linear eddy-viscosity model combined with the laminar kinetic energy. A non-linear model, compared to a linear one, has the advantage to resolve the anisotropic behavior of the Reynolds-stresses in the near-wall region and it provides a more accurate expression for the generation of turbulence in the transport equation of the turbulence kinetic energy. Although in its original formulation the model is not able to accurately predict the separation-induced transition, the inclusion of the laminar kinetic energy increases its accuracy. The adoption of the laminar kinetic energy by the non-linear model is presented in detail, together with some additional modifications required for the adaption of the laminar kinetic energy into the basic concepts of the non-linear eddy-viscosity model. The computational results using the proposed combined model are shown together with the ones obtained using an isotropic linear eddy-viscosity model, which adopts also the laminar kinetic energy concept and in comparison with the existing experimental data.     

无壁面参数低雷诺数非线性涡黏性模式研究

建立了一个低雷诺数的非线性涡黏性湍流模式,该模式的一个显著特性是它不包含壁面参数（如y^ ,n等）,因而特别适用于复杂几何流场的计算,本模式在几种包括回流、分离、激波等典型流动中进行了验证,结果令人满意。     

Performance of various RANS eddy-viscosity models for turbulent natural convection in tall vertical cavities

The present study is dedicated to the identification of turbulence models that are accurate and numerically economic for computing the natural air-flow and heat transfer by convection in tall cavities with differentially heated vertical walls. The eddy-viscosity models (EVM) are among the simplest to implement and the most economical to treat this problem. This study evaluated the dynamic, thermal and computational performances of twenty EVM turbulence models with one, two or three-equation closure. All the models were first implemented in several in-house codes using the finite volume method. The predictions of the retained models in terms of profiles of velocity, temperature and vertical velocity fluctuations in the cavity have been compared with those of experimental or numerical studies. The obtained results were used to identify the turbulence models that are accurate and numerically economic in predicting natural convection in vertical cavities with a high aspect ratio. The EVM models with three-equation (v2-f and ζ-f) provide the most accurate mean and fluctuating quantities, followed by the k-ε RNG (ReNormalization Group) and k-ω SST (Shear Stress Transport) models. The computing time of these four models is higher than that of the 2L (two-layer) and q-ω models, which provide fairly accurate results especially for the mean heat transfer between the vertical active walls. The other one-equation (Spalart and Allmaras model) and two-equation (k-ε, k-ω and hybrid models) turbulence models tested in this work, have a high computing time and/or predictions that are not sufficiently precise simultaneously for both velocity and temperature fields.     

Theoretical Investigation of an Eddy-Viscosity-Type Expression of the Reynolds Stress with Low-Reynolds-Number Effect

In low-Reynolds-number turbulent flows, the influence of the molecular viscosity is important. The turbulence models which are applied to those flows should include the low-Reynolds-number effect. In this study, turbulent flow with the molecular viscosity effect is analyzed theoretically with the aid of a two-scale direct-interaction approximation (TSDIA) and the energy spectrum and a new low-Reynolds-number-type eddy-viscosity representation are derived. An priori test for the model expression on the basis of the result of direct numerical simulation (DNS) for turbulent Couette flows is performed. Received 5 July 2002 and accepted 8 January 2003 Published online 25 March 2003 Communicated by T.B. Gatski     

Efficiency of Scale-Similarity Model for Study of Forced Compressible Magnetohydrodynamic Turbulence

Efficiency of scale-similarity model for study of forced compressible magnetohydrodynamic turbulence is studied. The scale-similarity model has several important advantages in contrast to the eddy-viscosity subgrid closures: good reproduction of the correlation between actual and model turbulent stress tensor even when the flow is highly anisotropic, and absence of special model constants. These advantages may be very essential for study of forced magnetohydrodynamic turbulence. Numerical computations under various similarity parameters are carried out and the obtained results are analyzed by means of comparison with results of direct numerical simulation and Smagorinsky closure for magnetohydrodynamics. Linear forcing algorithm is applied to keep the characteristics of turbulence stationary in time. Influence of discrete filter shapes on the scale-similarity model is studied as well. It is shown that the scale-similarity model provides good accuracy and the results agree well with the direct numerical simulation results. The present results show that the scale-similarity model might be a useful subgrid closure for study of scale-invariance properties of forced compressible magnetohydrodynamic turbulence in the inertial range and in contrast to decaying case the scale-similarity model can serve as a stand alone subgrid model.     

非线性湍流模式研究及进展

现代湍流模式研究已经超出了经典的Ｂｏｕｓｓｉｎｅｓｑ涡粘性概念和线性的雷诺应力输运范畴,湍流运动过程中的非线性本质已成为模式研究人员所关心的中心问题。其目的在于使湍流模式能更加真实地再现湍流运动的复杂性,提高模式的适用范围,使复杂湍流能够得到合理的模拟,非线性湍流模式在解决复杂湍流运动的计算中已经取得可喜进展,正逐步应用于工程湍流的计算。同时,工程中的湍流问题计算也已走出了简单剪切流动类型及传统的ｋ－ε（及其它形式的）二方程模式框架,二阶矩封闭模式在先进的工程计算中已被用来解决诸如可压缩的空气动力学、发动机气缸及三维复杂几何场内等具有重要应用背景的流动问题,并逐步进入计算流体力学商业软件包。