Numerical aspects of calculation of confined swirling flows with internal recirculation

Predictions were performed for two different confined swirling flows with internal recirculation zones. The convection terms in the elliptic governing equations were discretized using three different finite differencing schemes: hybrid, quadratic upwind interpolation and skew upwind differencing. For each flow case, calculations were carried out with these schemes and successively refined grids were employed. For the turbulent flow case the k-ε turbulence model was used. The predicted cases were a laminar swirling flow investigated by Bornstein and Escudier, and a turbulent low-swirl case studied by Roback and Johnson. In both cases an internal recirculation zone was present. The laminar case is well predicted when account is taken of the estimated radial velocity component at the chosen inlet plane. The quadratic upwind interpolation and skew upwind schemes predict the main features of the internal recirculation zone also with a coarse grid. The turbulent case is well predicted with the coarse as well as the finer grids, the skew upwind and quadratic upwind interpolation schemes yielding results very close to the measurements. It is concluded that the skew upwind scheme reaches grid independence slightly before the quadratic upwind scheme, both considerably earlier than the hybrid scheme.     

Development of a hybrid k-ε turbulence model for swirling recirculating flows under moderate to strong swirl intensities

A hybrid k-ε turbulence model, based on the concept that the modification of anisotropic effects should not be made in the flow regions inherent to small streamline curvatures, has been developed and examined with the swirling recirculating flows, with the swirl levels ranging from 0·6 to 1·23 in abrupt pipe expansion. A fairly satisfactory agreement of model predictions with the experimental data shows that this hybrid k-ε model can perform better simulation of swirling recirculating flows as compared to the standard k-ε model and the modified k-ε model proposed by Abujelala and Lilley.     

Computations of strongly swirling flows with second‐moment closures

The present study is concerned with simulating turbulent, strongly swirling flows by eddy viscosity model and Reynolds stress transport model variants adopting linear and quadratic form of the pressure–strain models. Flows with different inlet swirl numbers, 2.25 and 0.85, were investigated. Detailed comparisons of the predicted results and measurements were presented to assess the merits of model variants. For the swirl number 2.25 case, due to the inherent capability of the Reynolds stress models to capture the strong swirl and turbulence interaction, both the linear and quadratic form of the pressure–strain models predict the flow adequately. In strong contrast, the k–ϵ model predicts an excessively diffusive flow fields. For the swirl number 0.85 case, both the k–ϵ and Reynolds stress model with linear pressure–strain process, show an excessive diffusive transport of the flow fields. The quadratic pressure–strain model, on the other hand, mimics the correct flow development with the recirculating region being correctly predicted. Copyright © 1999 John Wiley & Sons, Ltd.     

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.     

Buoyancy effects on turbulent swirling flows

A three-parameter model of turbulence applicable to free boundary layers has been developed and applied for the prediction of axisymmetric turbulent swirling flows in uniform and stagnant surroundings under the action of buoyancy forces. The turbulent momentum and heat fluxes appearing in the time-averaged equations for the mean motion have been determined from algebraic expressions, derived by neglecting the convection and diffusion terms in the differential transport equations for these quantities, which relate the turbulent fluxes to the kinetic energy of turbulence, k, the dissipation length scale of turbulence, L, and the temperature covariance, T ′². Differential transport equations have been used to determine these latter quantities. The governing equations have been solved using fully implicit finite difference schemes. The turbulence model is capable of reproducing the gross features of pure jet flows, buoyant flows and swirling flows for weak and moderate swirl. The behaviour of a turbulent buoyant swirling jet has been found to depend solely on exit swirl and Froude numbers. The predicted results indicate that the incorporation of buoyancy can cause significant changes in the behaviour of a swirling jet, particularly when the buoyancy strength is high. The jet exhibits similarity behaviour in the initial region for weak swirl and weak buoyancy strengths only, and the asymptotic case of a swirling jet under the action of buoyancy forces is a pure plume in the far field. The predicted results have been found to be in satisfactory agreement with the available experimental data and in good qualitative agreement with other predicted results.     

Computation of strongly swirling confined flows with cubic eddy‐viscosity turbulence models

An investigation on the predictive performance of four cubic eddy‐viscosity turbulence models for two strongly swirling confined flows is presented. Comparisons of the prediction with the experiments show clearly the superiority of cubic models over the linear k–εmodel. The linear k–εmodel does not contain any mechanism to describe the stabilizing effects of swirling motion and as a consequence it performs poorly. Cubic models return a lower level of Reynolds stresses and the combined forced‐free vortex profiles of tangential velocity close to the measurements in response to the interaction between swirl‐induced curvature and stresses. However, a fully developed rotating pipe flow is too simple to contain enough flow physics, so the calibration of cubic terms is still a topic of investigation. It is shown that explicit algebraic stress models require fewer calibrations and contain more flow physics. Copyright © 2003 John Wiley & Sons, Ltd.     

Computation of confined coflow jets with three turbulence models

A numerical study of confined jets in a cylindrical duct is carried out to examine the performance of two recently proposed turbulence models: an RNG-based K-? model and a realizable Reynolds stress algebraic equation model. The former is of the same form as the standard K-? model but has different model coefficients. The latter uses an explicit quadratic stress-strain relationship to model the turbulent stresses and is capable of ensuring the positivity of each turbulent normal stress. The flow considered involves recirculation with unfixed separation and reatachment points and severe adverse pressure gradients, thereby providing a valuable test of the predictive capability of the models for complex flows. Calculations are performed with a finite volume procedure. Numerical credibility of the solutions is ensured by using second-order-accurate differencing schemes and sufficiently fine grids. Calculations with the standard K-? model are also made for comparison. Detailed comparisons with experiments show that the realizable Reynolds stress algebraic equation model consistently works better than does the standard K-? model in capturing the essential flow features, while the RNG-based K-? model does not seem to give improvements over the standard K-? model under the flow conditions considered.     

环形通道内湍流旋流流动的数值模拟

本文应用一种考虑湍流－旋流相互作用及湍流脉动各向异性的新的代数Ｒｅｙｎｏｌｄｓ应力模型,对环形通道内的湍流旋流流动进行了数值模拟,研究了改主为旋流流数,进口轴向速度及半径比等参数对环形通道内湍流流动的影响,以及对强化环形通道内传热的作用。     

Confined and Agitated Swirling Flows with Applications in Chemical Engineering

Turbulent, swirling flows are encountered frequently in chemical engineering practice. In this article, experiments and simulations on two classes of swirling flows, viz. agitated flows (stirred tanks), and confined swirling flows are discussed. Results of large-eddy simulations of stirred tank flow are compared with experimental data, mainly phase-resolved LDA data of the flow in the vicinity of the impeller. Next to the average velocity field, also the turbulent kinetic energy, and the anisotropy of the Reynolds stress tensor have been assessed. An important application of confined swirling flow is the cyclone separator (hydrocyclones for the separation of liquids, gas cyclones for gas-solid separation). The flow in a swirl tube geometry exhibiting many of the typical features of swirl flows (e.g. vortex breakdown) is discussed. Furthermore, a large-eddy simulation of the gas flow in a high-efficiency Stairmand cyclone separator is presented. Two examples of process modeling based on flow simulations are briefly treated: orthokinetic agglomeration of crystals in a stirred tank, and particle separation in a cyclone. This revised version was published online in July 2006 with corrections to the Cover Date.     

Numerical study of hysteresis in annular swirling jets with a stepped‐conical nozzle

This study investigates the experimentally observed hysteresis in the mean flow field of an annular swirling jet with a stepped‐conical nozzle. The flow is simulated using the Reynolds‐averaged Navier–Stokes (RANS) approach for incompressible flow with a k–ε and a Reynolds stress transport (RSTM) turbulence model. Four different flow structures are observed depending on the swirl number: ‘closed jet flow’, ‘open jet flow low swirl’, ‘open jet flow high swirl’ and ‘coanda jet flow’. These flow patterns change with varying swirl number and hysteresis at low and intermediate swirl numbers is revealed when increasing and subsequently decreasing the swirl. The influence of the inlet velocity profile on the transitional swirl numbers is investigated. When comparing computational fluid dynamics with experiments, the results show that both turbulence models predict the four different flow structures and the associated hysteresis and multiple solutions at low and intermediate swirl numbers. Therefore, a good agreement exists between experiments and numerics. Copyright © 2006 John Wiley & Sons, Ltd.     

Control volume finite element solution of confined turbulent swirling flows

A control volume finite element method that uses a triangular grid has been applied for solving confined turbulent swirling flows. To treat the velocity-pressure coupling, the vorticity-streamfunction formulation has been used. For turbulence effects the k-? model has been adopted. Consistent with the use of wall functions in the near-wall regions, a boundary condition for the calculation of the vorticity at computational boundaries is proposed and used effectively. The discretized equations are obtained by making use of an exponential interpolation function. Its use has been beneficial in reducing numerical diffusion. Comparisons of the current predictions with available experimental and numerical data from the literature showed generally fair agreement.     

不同参数下环形通道内湍流旋流流动的模拟

应用一种合理考虑湍流一旋流相互作用及湍流脉动各向异性的新的代数ＲｅｙｎｏｌｄＳ应力模型,对环形通道内的湍流旋流流动进行了数值模拟．研究了旋流数、进口轴向速度和内外半径比等参数对环形通道内湍流旋流流动的影响,以及由此产生的流场变化对强化环形通道内传热的作用．     

Numerical and experimental study of an axially induced swirling pipe flow

In-line flow segregators based on axial induction of swirling flow have important applications in chemical, process and petroleum production industries. In the later, the segregation of gas bubbles and/or water droplets dispersed into viscous oil by swirling pipe flow may be beneficial by either providing a pre-separation mechanism (bubble and/or drop coalescer) or, in the case of water-in-oil dispersions, by causing a water-lubricated flow pattern to establish in the pipe (friction reduction). Works addressing these applications are rare in the literature. In this paper, the features and capabilities of swirling pipe flow axially induced by a vane-type swirl generator were investigated both numerically and experimentally. The numerical analysis has been carried out using a commercial CFD package for axial Reynolds numbers less than 2000. Pressure drop, tangential and axial velocity components as well as swirl intensity along a 5 cm i.d. size and 3 m long pipe were computed. Single phase flow experiments have been performed using a water–glycerin solution of 54 mPa s viscosity and 1210 kg/m³ density as working fluid. The numerical predictions of the pressure drop were compared with the experimental data and agreement could be observed within the range of experimental conditions. The experiments confirmed that swirl flow leads to much higher friction factors compared with theoretical values for non-swirl (i.e. purely axial) flow. Furthermore, the addition of a conical trailing edge reduces vortex breakdown. Visualization of the two-phase swirling flow pattern was achieved by adding different amounts of air to the water–glycerin solution upstream the swirl generator.     

A new second‐moment closure approach for turbulent swirling confined flows

An improved anisotropic model for the dissipation rate—ε—of the turbulent kinetic energy (k), to be used together with a non‐linear pressure‐strain correlations model, is proposed. Experimental data from the open literature for two confined turbulent swirling flows are used to assess the performance of the proposed model in comparison to the standard ε transport equation and to a linear approach to model the pressure‐strain term that appears in the exact equations for the Reynolds‐stress tensor. For the less strongly swirling flow the predictions show much more sensitivity to the εtransport equation than to the pressure‐strain model. In opposition, for the more strongly swirling flow, the results show that the predictions are much sensitive to the pressure‐strain model. Nevertheless, the improved εtransport equation together with the non‐linear pressure strain model yield predictions in good agreement with experiments in both studied cases. Copyright © 2003 John Wiley & Sons, Ltd.     

Hot flow analysis of swirling sudden-expansion dump combustor

Hot flow of a sudden-expansion dump combustor with swirling is analysed by employing an infinite chemical reaction rate. Turbulence properties are closed using one type of algebraic Reynolds stress model and two types of κ–? model. One of the κ–? models includes a swirling effect modification to the ε-equation. Computations have been performed by the SIMPLE-C algorithm with a power-law scheme. The calculated results of the momentum fields and turbulence quantities for swirling flow are compared with the available experimental data. It is shown that the standard κ–? model gives poor prediction of the mean velocity, particularly the tangential velocity. For the hot flow analysis of a sudden-expansion dump combustor with swirling flow it is suggested that it is necessary to use the modified κ–? model or algebraic Reynolds stress model.     

BOUNDARY TREATMENT AND AN EFFICIENT PRESSURE ALGORITHM FOR INTERNAL TURBULENT FLOWS USING THE PDF METHOD

The generalized Langevin model, which is used to model the motion of stochastic particles in the velocity–composition joint probability density function (PDF) method for reacting turbulent flows, has been extended to incorporate solid wall effects. Anisotropy of Reynolds stresses in the near-wall region has been addressed. Numerical experiments have been performed to demonstrate that the forces in the near-wall region of a turbulent flow cause the stochastic particles approachi ng a solid wall to reverse their direction of motion normal to the wall and thereby, leave the near-wall layer. This new boundary treatment has subsequently been implemented in a full-scale problem to prove its validity. The test problem considered here is that of an isothermal, non-reacting turbulent flow in a two-dimensional channel with plug inflow and a fixed back-pressure. An efficient pressure correction method, developed in the spirit of the PISO algorithm, has been implemented. The pressure correction strategy is easy to implement and is completely consistent with the time- marching scheme used for the solution of the Lagrangian momentum equations. The results show remarkable agreement with both k–ϵ and algebraic Reynolds stress model calculations for the primary velocity. The secondary flow velocity and the turbulent moments are in better agreement with the algebraic Reynolds stress model predictions than the k– ϵ predictions. © 1997 by John Wiley & Sons, Ltd.     

RAS one-equation turbulence model with non-singular eddy-viscosity coefficient

A simplified consistency formulation for P_k/ε (production to dissipation ratio) is devised to obtain a non-singular C_μ (coefficient of eddy-viscosity) in the explicit algebraic Reynolds stress model of Gatski and Speziale. The coefficient C_μ depends non-linearly on both rotational/irrotational strains and is used in the framework of an improved RAS (Rahman–Agarwal–Siikonen) one-equation turbulence model to calculate a few well-documented turbulent flows, yielding predictions in good agreement with the direct numerical simulation and experimental data. The strain-dependent C_μ assists the RAS model in constructing the coefficients and functions such as to benefit complex flows with non-equilibrium turbulence. Comparisons with the Spalart–Allmaras one-equation model and the shear stress transport k-ω model demonstrate that C_μ improves the response of RAS model to non-equilibrium effects.     

On Algebraic Second Moment Models

Equilibrium and bifurcation analysis is used to explore algebraic second moment models. It is shown that the three-dimensional, explicit algebraic stress solution for the anisotropy tensor precludes rotational stabilization unless two invariants of the mean velocity gradient vanish. If these vanish the irrotational part of the flow must be a plane strain: essentially the model can only bifurcate and stabilize in two-dimensional mean flow. However, it is also shown that those same two invariants must vanish if the mean flow is steady. The full equilibrium analysis described herein provides a consistent picture of a model with equilibria that respond appropriately to rotation. However, if the algebraic stress approximation is used as a constitutive equation, without imposing full equilibrium, the bifurcation criterion $$\mathop {\lim }\limits_{\left| S \right| \to \infty } {\text{ }}b_{ij} S_{ij} = O(1)$$ will not be met in three-dimensional flow. Hence the model cannot bifurcate to the stable solution branch. Similarly, ad hoc non-linear constitutive formulas that do not satisfy the bifurcation criterion preclude rotational stabilization. The bifurcation criterion is a simple and powerful guidance to turbulence model formulations.     

Reynolds stress modelling of jet and swirl interaction inside a gas turbine combustor

Influences of the inlet swirl levels on the interaction between the dilution air jets and the swirling cross‐flow to the interior flow field inside a gas turbine combustor were investigated numerically by Reynolds stress transport model (RSTM). Due to the intense swirl and jet interaction, a high level of swirl momentum is transported to the centreline and hence, an intense vortex core is formed. The strength of the centreline vortex core was found to depend on the inlet swirl levels. For the higher swirling inlet, the decay of the swirling motion causes strong streamline variation of pressure; and consequently leads to an elevated level of deceleration of its axial velocity. Predictions contrasted with measurements indicate that the stress model reproduces the flow correctly and is able to reflect the influences of inlet swirl levels on the interior flow structure. Copyright © 1999 John Wiley & Sons, Ltd.