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

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

An Attempt to Combine Large Eddy Simulation with the k-ε Model in a Channel-Flow Calculation

Large eddy simulation (LES) is combined with the Reynolds-averaged Navier–Stokes (RANS) equation in a turbulent channel-flow calculation. A one-equation subgrid-scale model is solved in a three-dimensional grid in the near-wall region whereas the standard k–ε model is solved in a one-dimensional grid in the outer region away from the wall. The two grid systems are overlapped to connect the two models smoothly. A turbulent channel flow is calculated at Reynolds numbers higher than typical LES and several statistical quantities are examined. The mean velocity profile is in good agreement with the logarithmic law. The profile of the turbulent kinetic energy in the near-wall region is smoothly connected with that of the turbulent energy for the k–ε model in the outer region. Turbulence statistics show that the solution in the near-wall region is as accurate as a usual LES. The present approach is different from wall modeling in LES that uses a RANS model near the wall. The former is not as efficient as the latter for calculating high-Reynolds-number flows. Nevertheless, the present method of combining the two models is expected to pave the way for constructing a unified turbulence model that is useful for many purposes including wall modeling. Received 11 June 1999 and accepted 15 December 2000     

Numerical prediction of turbulent flow in a conical diffuser usingk-ε model

Fully developed incompressible turbulent flow in a conical diffuser having a total divergence angle of 8° and an area ratio of 4∶1 has been simulated by ak-ε turbulence model with high Reynolds number and adverse pressure gradient. The research has been done for pipe entry Reynolds numbers of 1.16×10⁵ and 2.93×10⁵. The mean flow velocity and turbulence energy are predicted successfully and the advantage of Boundary Fit Coordinates approach is discussed. Furthermore, thek-ε turbulence model is applied to a flow in a conical diffuser having a total divergence angle of 30° with a perforated screen. A simplified mathematical model, where only the pressure drop is considered, has been used for describing the effect of the perforated screen. The optimum combination of the resistance coefficient and the location of the perforated screen is predicted for high diffuser efficiency or the uniform velocity distribution.     

Three-Dimensional Calculations of the Flow Field around a Turbine Blade with Film Cooling Injection near the Leading Edge

Injection of coolant air from a showerhead injection system at the leading edge of a high pressure turbine blade is investigated using a fully implicit three-dimensional finite-volume method on multi-block grids. For various blowing rates, the calculation results for the velocity and pressure fields and turbulence intensity are compared with available experimental data. The present method yields excellent agreement with the experiments for the isentropic Mach number distributions on the blade surface. The standardk–ε turbulence model with wall functions is already capable of capturing the major details of the flow field including the injection-induced secondary-flow vortices, particularly so on the suction side. On the pressure side, however, the lateral jet spreading is under-predicted somewhat together with an exaggeration of the near-wall sink-flow vortices. On this side with convex walls, where turbulence anisotropy is appreciable according to the experiments, overall better predictions were obtained with the anisotropy correction of Bergeles et al. [23] promoting the Reynolds stress in the lateral direction. The correction has no beneficial effect on the suction side with concave walls where the turbulence anisotropy was observed to be much smaller. This revised version was published online in July 2006 with corrections to the Cover Date.     

A turbulent plane offset jet with small offset ratio

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

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.     

Prediction of Flow and Heat Transfer through Stationary and Rotating Ribbed Ducts Using a Non-linear <Emphasis Type="Italic">k</Emphasis>−<Emphasis Type="Italic">ε</Emphasis> Model

The present paper deals with the prediction of three-dimensional fluid flow and heat transfer in rib-roughened ducts of square cross-section, which are either stationary, or rotate in orthogonal mode. The main objective is to assess how a recently developed variant of a cubic non-linear k−ε model (proposed by Craft et al. Flow Turbul Combust 63:59–80, 1999) can predict three-dimensional flow and heat transfer characteristics through stationary and rotating ribbed ducts. The present paper discusses turbulent air flow and heat transfer through two different configurations, namely: (I) a stationary square duct with “in-line” normal and (II) a square duct with normal ribs in a “staggered” arrangement under stationary and rotating conditions, with the axis of rotation normal to the flow direction and parallel to the ribs. In this paper the flow and thermal predictions of the linear k−ε model (EVM) are also included, as a set of baseline predictions. The mean flow predictions show that both linear and non-linear k−ε models can successfully reproduce most of the measured data for stream-wise and cross-stream velocity components. Moreover, the non-linear model is able to produce better results for the turbulent stresses. The heat transfer predictions show that both EVM and NLEVM2, the more recent variant of the non-linear k−ε, with the algebraic length-scale correction term, overestimate the measured Nusselt numbers for both geometries examined. While the EVM with the differential length-scale correction term underestimates heat transfer levels, the Nusselt number predictions with the NLEVM2 and the ‘NYP’ term are in close agreements with the measured data. Comparisons with our earlier work, Iacovides and Raisee (Int J Heat Fluid Flow, 20:320–328, 1999), show that the NLEVM2 thermal predictions are of similar quality to those of a second-moment closure.     

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

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

Numerical Simulation of Delft-Jet-in-Hot-Coflow (DJHC) Flames Using the Eddy Dissipation Concept Model for Turbulence–Chemistry Interaction

In this paper, we report results of a numerical investigation of turbulent natural gas combustion for a jet in a coflow of lean combustion products in the Delft-Jet-in-Hot-Coflow (DJHC) burner which emulates MILD (Moderate and Intense Low Oxygen Dilution) combustion behavior. The focus is on assessing the performance of the Eddy Dissipation Concept (EDC) model in combination with two-equation turbulence models and chemical kinetic schemes for about 20 species (Correa mechanism and DRM19 mechanism) by comparing predictions with experimental measurements. We study two different flame conditions corresponding to two different oxygen levels (7.6% and 10.9% by mass) in the hot coflow, and for two jet Reynolds number (Re = 4,100 and Re = 8,800). The mean velocity and turbulent kinetic energy predicted by different turbulence models are in good agreement with data without exhibiting large differences among the model predictions. The realizable k-ε model exhibits better performance in the prediction of entrainment. The EDC combustion model predicts too early ignition leading to a peak in the radial mean temperature profile at too low axial distance. However the model correctly predicts the experimentally observed decreasing trend of lift-off height with jet Reynolds number. A detailed analysis of the mean reaction rate of the EDC model is made and as possible cause for the deviations between model predictions and experiments a low turbulent Reynolds number effect is identified. Using modified EDC model constants prediction of too early ignition can be avoided. The results are weakly sensitive to the sub-model for laminar viscosity and laminar diffusion fluxes.     

Turbulence Modelling of Unsteady Turbulent Flows Using the Stress Strain Lag Model

This paper reports the application of a recently developed turbulence modelling scheme known as the C _as model. This model was specifically developed to capture the effects of stress-strain misalignment observed in turbulent flows with mean unsteadiness. Earlier work has reported the approach applied within a linear k-ε modelling framework, and some initial testing of it within the k-ω SST model of Menter (AIAA J 32:1598–1605, 1994). The resulting k-ε-C _as and SST-C _as models have been shown to result in some of the advantages of a full Reynolds Stress transport Model (RSM), whilst retaining the computational efficiency and stability benefits of a eddy viscosity model (EVM). Here, the development of the the high-Reynolds-number version of the C _as model is outlined, with some example applications to steady and unsteady homogeneous shear flows. The SST-C _as form of the model is then applied to further, more challenging cases of 2-D flow around a NACA0012 aerofoil beyond stall and the 3-D flow around a circular cylinder in a square duct, both being flows which exhibit large, unsteady, separated flow regions. The predictions returned by a range of other common turbulence modelling schemes are included for comparison and the SST-C _as scheme is shown to return generally good results, comparable in some respects to those obtainable from far more complex schemes, for only moderate computing resource requirements.     

Testing of Model Equations for the Mean Dissipation using Kolmogorov Flows

Direct Numerical Simulations (DNS) of Kolmogorov flows are performed at three different Reynolds numbers Re _λ between 110 and 190 by imposing a mean velocity profile in y-direction of the form U(y) = F sin(y) in a periodic box of volume (2π)³. After a few integral times the turbulent flow turns out to be statistically steady. Profiles of mean quantities are then obtained by averaging over planes at constant y. Based on these profiles two different model equations for the mean dissipation ε in the context of two-equation RANS (Reynolds Averaged Navier–Stokes) modelling of turbulence are compared to each other. The high Reynolds number version of the k-ε-model (Jones and Launder, Int J Heat Mass Transfer 15:301–314, 1972), to be called the standard model and a new model by Menter et al. (2006), to be called the Menter–Egorov model, are tested against the DNS results. Both models are solved numerically and it is found that the standard model does not provide a steady solution for the present case, while the Menter–Egorov model does. In addition a fairly good quantitative agreement of the model solution and the DNS data is found for the averaged profiles of the kinetic energy k and the dissipation ε. Furthermore, an analysis based on flow-inherent geometries, called dissipation elements (Wang and Peters, J Fluid Mech 608:113–138, 2008), is used to examine the Menter–Egorov ε model equation. An expression for the evolution of ε is derived by taking appropriate moments of the equation for the evolution of the probability density function (pdf) of the length of dissipation elements. A term-by-term comparison with the model equation allows a prediction of the constants, which with increasing Reynolds number approach the empirical values.     

On a Generalized Nonlinear K–ε Model and the Use of Extended Thermodynamics in Turbulence

A resent extension of the nonlinear K–ε model is critically discussed from a basic theoretical standpoint. While it was said in the paper that this model was formulated to incorporate relaxation effects, it will be shown that the model is incapable of describing one of the most basic such turbulent flows as is obvious but is described for clarity. It will be shown in detail that this generalized nonlinear K–ε model yields erroneous results for the Reynolds stress tensor when the mean strains are set to zero in a turbulent flow – the return-to-isotropy problem which is one of the most elementary relaxational turbulent flows. It is clear that K–ε type models cannot describe relaxation effects. While their general formalism can describe relaxation effects, the nonlinear K–ε model – which the paper is centered on – cannot. The deviatoric part of the Reynolds stress tensor is predicted to be zero when it actually only gradually relaxes to zero. Since this model was formulated by using the extended thermodynamics, it too will be critically assessed. It will be argued that there is an unsubstantial physical basis for the use of extended thermodynamics in turbulence. The role of Material Frame-Indifference and the implications for future research in turbulence modeling are also discussed. Received 19 February 1998 and accepted 23 October 1998     

Influence of geometry on separation efficiency in a hydrocyclone

A numerical study of the gas–liquid–solid multiphase flow in hydrocyclones is presented. Three models of turbulence, the RNG k–ε model, the Reynolds stress model and Large eddy simulation with the volume of fluid model (VOF) multiphase model for simulating air core are compared in order to predict axial and tangential velocity distributions. This presentation is mainly aimed at identifying an optimal method, used to study effective parameters, based on which, eventually, effect of inlet flow rate variations and body dimension variations such as underflow diameter, overflow diameter and cone angle on the separation performance and pressure drop are investigated. The results are then used in the simulation of particle flow described by the stochastic Lagrangian model. The results suggest that the predicted size classifications are approximately similar to those of RSM and LES methods. Predictions using the RSM model are found in agreement with experimental results with a marginal error within the range of 4 to 8%. Proceeding model validation, parametric studies have been carried out concerning the influence of velocity inlet, particle size and body dimension such as underflow and overflow diameter and cone angle. The predictions demonstrate that the flow fields in the hydrocyclones with different sizes and lengths are different, which yields different performances.     

Turbulence kinetic energy budget in bubbly flows in a vertical duct

Understanding turbulence kinetic energy (TKE) budget in gas–liquid two-phase bubbly flows is indispensable to develop and improve turbulence models for the bubbly flows. In this study, a molecular tagging velocimetry based on photobleaching reaction was applied to turbulent bubbly flows with sub-millimeter bubbles in a vertical square duct to examine the applicability of the k–ε models to the bubbly flows. Effects of bubbles on TKE budget are discussed and a priori tests of the standard and low Reynolds number k–ε models are carried out to examine the applicability of these models to the bubbly flows. The conclusions obtained are as follows: (1) The photobleaching molecular tagging velocimetry is of use for validating turbulence models. (2) The bubbles increase the liquid velocity gradient in the near wall region, and therefore, enhance the production and dissipation rates of TKE. (3) The k–ε models can reasonably evaluate the production rate of TKE in the bubbly flows. (4) The modulations of diffusion due to the bubbles have different characteristics from the diffusion enhancement due to shear-induced turbulence. Hence, the k–ε models fail in evaluating the diffusion rate in the near wall region in the bubbly flows. (5) The k–ε models represent the trends of the production, dissipation, and diffusion rates of ε in the bubbly flow, although more accurate experimental data are required for quantitative validation of the ε equation.     

Longitudinal vortex pair imbedded in a turbulent boundary layer: Comparative performances of several turbulence models

Four turbulence models, namely, the basic and nonlinear stress-transport models and the basic and anisotropick-ε models, have been tested in the case of interaction between a longitudinal vortex pair and a flat-plate boundary layer. The results of their predictions were compared with Mehta and Bradshaw's measurements. In this paper, part of the results involving those of the nonlinear stress-transport model and anisotropick-ε model are presented and discussed. The project supported by the National Natural Science Foundation of China under Contract No. 19132012     

A Turbulence Model Sensitivity Study for CH4/H2 Bluff-Body Stabilized Flames

The objective of this work is to assess the performances of different turbulence models in predicting turbulent diffusion flames in conjunction with the flamelet model.The k– model, the Explicit Algebraic Stress Model (EASM) and the k– model withvaried anisotropy parameter C (LEA k– model)are first applied to the inert turbulent flow over a backward-facing step, demonstrating the quality of the turbulence models. Following this, theyare used to simulate the CH₄/H₂ bluff-body flame studied by the University of Sydney/Sandia.The numerical results are compared to experimental values of the mixture fraction, velocity field, temperature and constituent mass fractions.The comparisons show that the overall result depends on the turbulence model used, and indicate that theEASM and the LEA k– models perform better than the k– model and mimic most of the significant flow features.     

Scalar turbulence model investigation with variable turbulent Prandtl number in heated jets and diffusion flames

Traditional turbulence models using constant turbulent Prandtl number fail to predict the experimentally observed anisotropies in the thermal eddy diffusivity and thermal turbulent intensity fields. Accurate predictions depend strongly on the turbulence model employed. Consequently, the objective of this paper is to assess the performance of turbulence model with variable turbulent Prandtl number in predicting of thermal and scalar fields quantities. The model is applied to axisymmetric turbulent round jet with variable density and in turbulent hydrogen diffusion flames using the flamelet concept. The k − ɛ turbulence model is used in conjunction with thermal field; the model involves solving supplemental scalar equations for the temperature variance and its dissipation rate. The model predictions are compared with available experimental data for the purpose of validating model. In reacting cases, velocity and scalar (including temperature and mass fractions) predictions agree relatively well in the near field of the investigated diluted hydrogen flames.     

Numerical simulation of turbulent flow through series stenoses

The flow fields in the neighbourhoods of series vascular stenoses are studied numerically for the Reynolds numbers from 100 to 4000, diameter constriction ratios of 0.2–0.6 and spacing ratios of 1, 2, 3, 4 and ∞. In this study, it has been further verified that in the laminar flow region, the numerical predictions by k–ω turbulence model matched those by the laminar‐flow modelling very well. This suggests that the k–ω turbulence model is capable of the prediction of the laminar flow as well as the prediction of the turbulent stenotic flow with good accuracy. The extent of the spreading of the recirculation region from the first stenosis and its effects on the flow field downstream of the second stenosis depend on the stenosis spacing ratio, constriction ratio and the Reynolds number. For c₁ = 0.5 with c₂ ≤ c₁, the peak value of wall vorticity generated by the second stenosis is always less than that generated by the first stenosis. However, the maximum centreline velocity and turbulence intensity at the second stenosis are higher than those at the first stenosis. In contrast, for c₁ = 0.5 with c₂ = 0.6, the maximum values at the second stenosis are much higher than those at the first stenosis whether for centreline velocity and turbulence intensity or for wall vorticity. The peak values of the wall vorticity and the centreline disturbance intensity both grow up with the Reynolds number increasing. The present study shows that the more stenoses can result in a lower critical Reynolds number that means an earlier occurrence of turbulence for the stenotic flows. Copyright © 2003 John Wiley & Sons, Ltd.     

Simulation of Wind Flow Around a Building with a k–ε Model

The three-dimensional numerical simulation of airflow around a building using a k–ε two-equation turbulence model is presented in this paper. Several cases of numerical simulation of airflow around a building are carried out to estimate the influence of mesh spacing on simulated results. The accuracy of simulations is examined by comparing the predicted results with wind-tunnel experiments. It is confirmed that numerical simulations by means of the k–ε model reproduce the velocity fields well when using fine mesh resolution. In the latter part of the paper, the simulation method is applied to predict the flow field around a building with different width-to-height ratios, under light wind conditions. Received 16 June 1999 and accepted 20 July 2000     

Prediction of a Turbulent Non-Premixed Natural Gas Flame in a Semi-Industrial Scale Furnace using a Radiative Flamelet Combustion Model

A mixedness-reactedness flamelet combustion model coupled with a comprehensive radiation heat transfer model based on the discrete transfer method of solution of the radiative transport equation is applied for the simulation of a 3 MW non-swirling turbulent non-premixed natural gas flame in the experimental furnace at the International Flame Research Foundation. In the calculation, turbulence is represented by the standard k − ε and a differential Reynolds-stress model. Predictions are compared with measurements of mean gas velocity, temperature, major species concentrations and incident radiation wall flux. The radiative mixedness-reactedness flamelet combustion model, irrespective of the model for turbulence, is able to reproduce the basic structure of the experimental flame, which is stabilised downstream of the burner nozzle. In the near burner region, encompassing the non-reacting lift-off zone, good quality predictions are obtained using both the turbulence models, whereas further downstream, within the combusting zone of the jet, the Reynolds-stress turbulence model generates better predictions at and about the furnace axis. The nitric oxide (NO) formation via the thermal- and prompt-NO routes was also calculated and compared with in-flame and flue-gas NO data. The measured NO level at the furnace exit is well reproduced in the calculation, however discrepancies exist near the burner where NO concentrations around the furnace axis are overpredicted.