K-ε-T model of dense liquid-solid two-phase turbulent flow and its application to the pipe flow

To predict the characteristics of dense liquid-solid two-phase flow, K-ε-T model is established, in which the turbulent flow of fluid phase is described with fluid turbulent kinetic energy K_f and its dissipation rate ε_f, and the particles random motion is described with particle turbulent energy K_p and its dissipation rate ε_p and pseudothermal temperature T_p. The governing equations are also derived. With K-ε-T model, numerical study of dense liquid-solid two-phase turbulent up-flow in a pipe is performed. The calculated results are in good agreement with experimental data of Alajbegovic et al. (1994), and some flow features are captured.     

Finite element study of free surface turbulent flow in rotating channels

Under certain conditions of liquid flow through rotating channels, the Coriolis force can induce a free surface to be formed. This problem is of practical importance in a Coriolis wear tester, which is used for determining the sliding wear coefficient of wear materials in slurry handling equipment. A deforming Galerkin finite element method is presented for predicting two‐dimensional turbulent free surface mean flow in rotating channels. Reynolds‐averaged Navier–Stokes (RANS) equations are cast into weak(algebraic) form using primitive variables (velocity and pressure). Eddy viscosity is determined via a mixing length model. Velocity is interpolated biquadratically, while pressure is interpolated bilinearly. The kinematic condition is used to form the Galerkin residual for the free surface. The free surface is represented by Hermite polynomials of zeroeth order for continuity of position and slope. Combined Newton's iteration is used to simultaneously solve for the free surface and the field variables. Results of velocity and pressure fields, as well as the free surface are shown to converge with mesh‐size refinement. There is excellent respect for mass conservation. Results are presented for various values of Rossby number (Ro) and height‐based Reynolds number (Re_H). Parameter continuation in Ro and Re_H space is used to compute solutions at higher values of flow rate and angular velocity. Copyright © 2002 John Wiley & Sons, Ltd.     

Turbulent free surface flow simulation using a multilayer model

Based on the steady hydrodynamic equations, a multilayer (ML) model has been formulated for simulating turbulent flow in open channels. The model is imposed on a general curvilinear co-ordinate system with non-staggered finite volume discretization. The turbulent quantities in the model are described by the layer-averaged K-ε turbulence model with standard coefficients. Assuming a vertical hydrostatic pressure distribution, a depth correction scheme, originating in the Rhie and Chow approach for confined flows, is incorporated into the SIMPLE procedure to compute the water surface. Using the multilayer model, flows in a 180° channel bend, near a groin, and in straight open channels are computed. The results are compared with experimental data and with calculations of a depth-averaged model (DAV) having three-dimensional effect corrections. The comparisons show that the predictions of the ML model on mean flow values are in good agreement with the available data and are better than those of the DAV model. The vertical distribution of the turbulent energy dissipation rate is also shown to agree well with the open-channel measurements.     

SECOND-ORDER MOMENT MODEL FOR DENSE TWO-PHASE TURBULENT FLOW OF BINGHAM FLUID WITH PARTICLES

The USM-θmodel of Bingham fluid for dense two-phase turbulent flow was developed, which combines the second-order moment model for two-phase turbulence with the particle kinetic theory for the inter-particle collision. In this model, phases interaction and the extra term of Bingham fluid yield stress are taken into account. An algorithm for USM-θmodel in dense two-phase flow was proposed, in which the influence of particle volume fraction is accounted for. This model was used to simulate turbulent flow of Bingham fluid single-phase and dense liquid-particle two-phase in pipe. It is shown USM-θmodel has better prediction result than the five-equation model, in which the particle-particle collision is modeled by the particle kinetic theory, while the turbulence of both phase is simulated by the two-equation turbulence model. The USM-θmodel was then used to simulate the dense two-phase turbulent up flow of Bingham fluid with particles. With the increasing of the yield stress, the velocities of Bingham and particle decrease near the pipe centre. Comparing the two-phase flow of Bingham-particle with that of liquid-particle, it is found the source term of yield stress has significant effect on flow.     

Predictive Capabilities of an Improved Cubic k–ε Model for Inert Steady Flows

Through an improved ε transport equation, a major quality enhancement of the cubic k–ε model, earlier developed in[13], is obtained. The ε-equation of [13],yielding good results for wall-bounded and rotating flows, is combined with the one derived by Shih et al. [20], which produces good results for free shear flows (e.g. the plane jet–round jet anomaly is resolved).Results are presented for the following flows: fully developed stationary and rotating channel and pipe, backward-facing step, sudden pipe expansion, smooth channel expansion and contraction, plane and round jet. Heat transfer predictions in turbulent impinging jets are also discussed. Accurate results are obtained for the mean flow quantities for all test cases, without case dependent model tuning. This revised version was published online in July 2006 with corrections to the Cover Date.     

Numerical simulation of turbulent free surface flow with two‐equation k–ε eddy‐viscosity models

This paper presents a finite difference technique for solving incompressible turbulent free surface fluid flow problems. The closure of the time‐averaged Navier–Stokes equations is achieved by using the two‐equation eddy‐viscosity model: the high‐Reynolds k–ε (standard) model, with a time scale proposed by Durbin; and a low‐Reynolds number form of the standard k–ε model, similar to that proposed by Yang and Shih. In order to achieve an accurate discretization of the non‐linear terms, a second/third‐order upwinding technique is adopted. The computational method is validated by applying it to the flat plate boundary layer problem and to impinging jet flows. The method is then applied to a turbulent planar jet flow beneath and parallel to a free surface. Computations show that the high‐Reynolds k–ε model yields favourable predictions both of the zero‐pressure‐gradient turbulent boundary layer on a flat plate and jet impingement flows. However, the results using the low‐Reynolds number form of the k–ε model are somewhat unsatisfactory. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical analysis of turbulent structure in compound meandering open channel by algebraic Reynolds stress model

The turbulent flow in a compound meandering channel with a rectangular cross section is one of the most complicated turbulent flows, because the flow behaviour is influenced by several kinds of forces, including centrifugal forces, pressure‐driven forces and shear stresses generated by momentum transfer between the main channel and the flood plain. Numerical analysis has been performed for the fully developed turbulent flow in a compound meandering open‐channel flow using an algebraic Reynolds stress model. The boundary‐fitted coordinate system is introduced as a method for coordinate transformation in order to set the boundary conditions along the complicated shape of the meandering open channel. The turbulence model consists of transport equations for turbulent energy and dissipation, in conjunction with an algebraic stress model based on the Reynolds stress transport equations. With reference to the pressure–strain term, we have made use of a modified pressure–strain term. The boundary condition of the fluctuating vertical velocity is set to zero not only for the free surface, but also for computational grid points next to the free surface, because experimental results have shown that the fluctuating vertical velocity approaches zero near the free surface. In order to examine the validity of the present numerical method and the turbulent model, the calculated results are compared with experimental data measured by laser Doppler anemometer. In addition, the compound meandering open channel is clarified somewhat based on the calculated results. As a result of the analysis, the present algebraic Reynolds stress model is shown to be able to reasonably predict the turbulent flow in a compound meandering open channel. Copyright © 2005 John Wiley & Sons, Ltd.     

Parameter estimation of engineering turbulence model

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

A general one-equation turbulence model for free shear and wall-bounded flows

The purpose of thiswork is to introduce a complete and general one-equation model capable of correctly predicting a wide class of fundamental turbulent flows like boundary layer, wake, jet, and vortical flows. The starting point is the mature and validated two-equation k−ω turbulence model of Wilcox. The newly derived one-equation model has several advantages and yields better predictions than the Spalart-Allmaras model for jet and vortical flows while retaining the same efficiency and quality of the results for near-wall turbulent flows without using a wall distance. The derivation and validation of the new model using findings computed by the Spalart-Allmaras and the k−ω models are presented and discussed for several free shear and wall-bounded flows.     

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.     

Buffeting in transonic flow prediction using time‐dependent turbulence model

In transonic flow conditions, the shock wave/turbulent boundary layer interaction and flow separations on wing upper surface induce flow instabilities, ‘buffet’, and then the buffeting (structure vibrations). This phenomenon can greatly influence the aerodynamic performance. These flow excitations are self‐sustained and lead to a surface effort due to pressure fluctuations. They can produce enough energy to excite the structure. The objective of the present work is to predict this unsteady phenomenon correctly by using unsteady Navier–Stokes‐averaged equations with a time‐dependent turbulence model based on the suitable (k–ε) turbulent eddy viscosity model. The model used is based on the turbulent viscosity concept where the turbulent viscosity coefficient (Cμ) is related to local deformation and rotation rates. To validate this model, flow over a flat plate at Mach number of 0.6 is first computed, then the flow around a NACA0012 airfoil. The comparison with the analytical and experimental results shows a good agreement. The ONERA OAT15A transonic airfoil was chosen to describe buffeting phenomena. Numerical simulations are done by using a Navier–Stokes SUPG (streamline upwind Petrov–Galerkin) finite‐element solver. Computational results show the ability of the present model to predict physical phenomena of the flow oscillations. The unsteady shock wave/boundary layer interaction is described. Copyright © 2005 John Wiley & Sons, Ltd.     

Development of an Analytical <Emphasis Type="Italic">β</Emphasis>-Function PDF Integration Algorithm for Simulation of Non-premixed Turbulent Combustion

Among many presumed-shape pdf approaches for modeling non-premixed turbulent combustion, the presumed β-function pdf is widely used in the literature. However, numerical integration of the β-function pdf may encounter singularity difficulties at mixture fraction values of Z = 0 or 1. To date, this issue has been addressed by few publications. The present study proposes the Piecewise Integration Method (PIM), an efficient, robust and accurate algorithm to overcome these numerical difficulties with the added benefit of improving computational efficiency. Comparison of this method to the existing numerical integration methods shows that the PIM exhibits better accuracy and greatly increases computational efficiency. The PIM treatment of the β-function pdf integration is first applied to the Burke–Schumann solution in conjunction with the k − ε turbulence model to simulate a CH₄/H₂ bluff-body turbulent flame. The proposed new method is then applied to the same flow using a more complex combustion model, the laminar flamelet model. Numerical predictions obtained by using the proposed β-function pdf integration method are compared to experimental values of the velocity field, temperature and species mass fractions to illustrate the efficiency and accuracy of the present method.     

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.     

Large eddy simulation of vertical turbulent jets under JONSWAP waves

The effect of random waves on vertical plane turbulent jets is studied numerically and the mechanism behind the interaction of the jet and waves is analyzed. The large eddy simulation method is used and the σ-coordinate system is adopted. Turbulence is modeled by a dynamic coherent eddy model. The σ-coordinate transformation is introduced to map the irregular physical domain with a wavy free surface and an uneven bottom onto a regular computational domain. The fractional step method is used to solve the fil...     

Measurement of the gradient field of a turbulent free surface

We study the free surface above a turbulent channel flow. We describe a laser scanning technique that can be used to measure the space–time turbulent surface gradient field along a line. A harmonically swiveling laser beam is focused on the surface and its angle of refraction is measured using a position sensing device. The registered signals can be converted easily to the desired gradient field, and spectra and correlations can be measured. Examples of measured spectra and correlation functions of the surface above a turbulent channel flow (Reynolds number R _λ ≈ 250) demonstrate the viability of the technique. We further assess the validity of Taylor’s frozen turbulence hypothesis that implies that time-dependent signals measured along a line that is oriented perpendicularly to the mean channel velocity can be interpreted as 2D measurements of the surface slope. While Taylor’s hypothesis works for a turbulent velocity field, it does not work for its free surface.

Interaction of turbulent natural convection and surface thermal radiation in inclined square enclosures

Two-dimensional numerical studies of flow and temperature fields for turbulent natural convection and surface radiation in inclined differentially heated enclosures are performed. Investigations are carried out over a wide range of Rayleigh numbers from 10⁸ to 10¹², with the angle of inclination varying between 0° and 90°. Turbulence is modeled with a novel variant of the k–ε closure model. The predicted results are validated against experimental and numerical results reported in literature. The effect of the inclination of the enclosure on pure turbulent natural convection and the latter’s interaction with surface radiation are brought out. Profiles of turbulent kinetic energy and effective viscosity are studied to observe the net effect on the intensity of turbulence caused by the interaction of natural convection and surface radiation. The variations of local Nusselt number and average Nusselt number are presented for various inclination angles. Marked change in the convective Nusselt number is found with the orientation of enclosure. Also analyzed is the influence of change in emissivity on the flow and heat transfer. A correlation relevant to practical applications in the form of average Nusselt number, as a function of Rayleigh number, Ra, radiation convection parameter, N _RC and inclination angle of the enclosure, φ is proposed.     

Numerical simulation of turbulent free-surface flow in curved channel

This paper presents first results of numerical simulation of turbulent free-surface flow. Simple implementation of surface capturing method is based on the variable density approach. The flow is treated as if there is only one fluid, but with variable material properties (density, viscosity). The switch in these values is done by a function resulting from the mass conservation principle. This approach simplifies the implementation of turbulence model. In this case the SST k−ω model was chosen in modification given by Hellsten.Numerical solution was carried out by finite-volume method with explicit Runge-Kutta time-integration. The artificial compressibility method was used for time-marching search for steady state solution. The whole model was tested on horizontally placed square-sectioned 90^∘ bend, which was partially filled by the water. The main goal of this study was to demonstrate the applicability of this model and solution method for capturing the water-air interface as well as for predicting the turbulent effects in both fluids.     

One Equation Model for Turbulent Channel Flow with Second Order Viscoelastic Corrections

A modified second order viscoelastic constitutive equation is used to derive a k–l type turbulence closure to qualitatively assess the effects of elastic stresses on fully-developed channel flow. Specifically, the second order correction to the Newtonian constitutive equation gives rise to a new term in the momentum equation involving the time-averaged elastic shear stress and in the turbulent kinetic energy transport equation quantifying the interaction between the fluctuating elastic stress and rate of strain tensors, denoted by P _w , for which a closure is developed and tested. This closure is based on arguments of isotropic turbulence and equilibrium in boundary layer flows and a priori P _w could be either positive or negative. When P _w is positive, it acts to reduce the production of turbulent kinetic energy and the turbulence model predictions qualitatively agree with direct numerical simulation (DNS) results obtained for more realistic viscoelastic fluid models with memory which exhibit drag reduction. In contrast, P _w  < 0 leads to a drag increase and numerical breakdown of the model occurs at very low values of the Deborah number, which signifies the ratio of elastic to viscous stresses. Limitations of the turbulence model primarily stem from the inadequacy of the k–l formulation rather than from the closure for P _w . An alternative closure for P _w , mimicking the viscoelastic stress work predicted by DNS using the Finitely Extensible Nonlinear Elastic-Peterlin fluid model, which is mostly characterized by P _w  > 0 but has also a small region of negative P _w in the buffer layer, was also successfully tested. This second model for P _w leads to predictions of drag reduction, in spite of the enhancement of turbulence production very close to the wall, but the equilibrium conditions in the inertial sub-layer were not strictly maintained.     

Numerical modelling of complex turbulent free‐surface flows with the SPH method: an overview

The gridless smoothed particle hydrodynamics (SPH) method is now commonly used in computational fluid dynamics (CFD) and appears to be promising in predicting complex free‐surface flows. However, increasing flow complexity requires appropriate approaches for taking account of turbulent effects, whereas some authors are still working without any turbulence closure in SPH. A review of recently developed turbulence models adapted to the SPH method is presented herein, from the simplistic point of view of a one‐equation model involving mixing length to more sophisticated (and thus realistic) models like explicit algebraic Reynolds stress models (EARSM) or large eddy simulation (LES). Each proposed model is tested and validated on the basis of schematic cases for which laboratory data, theoretical or numerical solutions are available in the general field of turbulent free‐surface incompressible flows (e.g. open‐channel flow and schematic dam break). They give satisfactory results, even though some progress should be made in the future in terms of free‐surface influence and wall conditions. Recommendations are given to SPH users to apply this method to the modelling of complex free‐surface turbulent flows. Copyright © 2006 John Wiley & Sons, Ltd.     

Turbulent flow simulation of liquid jet emanating from pressure-swirl atomizer

The performances of three linear eddy viscosity models (LEVM) and one algebraic Reynolds stress model (ARSM) for the simulation of turbulent flow inside and outside pressure-swirl atomizer are evaluated by comparing the interface position with available experimental data and by comparing the turbulence intensity profiles at the atomizer exit. 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. The turbulence intensity has a significant influence on the global characteristics of the flow field. The turbulence models with better predictions of the turbulence intensity, such as Gatski-Speziale’s ARSM 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, or the discharge coefficient, film thickness and the liquid sheet outer surface position for the atomizer flows. The standard k–ε model predicts stronger turbulence intensity as compared to the other models and therefore yields smaller film thickness and larger liquid sheet outer surface position. In average, the ARSM model gives both quantitatively and qualitatively better results as compared to the standard k–ε model and the low Reynolds number models.