Calculation of curved open channel flow using physical curvilinear non‐orthogonal co‐ordinates

There are two main difficulties in numerical simulation calculations using FD/FV method for the flows in real rivers. Firstly, the boundaries are very complex and secondly, the generated grid is usually very non‐uniform locally. Some numerical models in this field solve the first difficulty by the use of physical curvilinear orthogonal co‐ordinates. However, it is very difficult to generate an orthogonal grid for real rivers and the orthogonal restriction often forces the grid to be over concentrated where high resolution is not required. Recently, more and more models solve the first difficulty by the use of generalized curvilinear co‐ordinates (ξ,η). The governing equations are expressed in a covariant or contra‐variant form in terms of generalized curvilinearco‐ordinates (ξ,η). However, some studies in real rivers indicate that this kind of method has some undesirable mesh sensitivities. Sharp differences in adjacent mesh size may easily lead to a calculation stability problem oreven a false simulation result. Both approaches used presently have their own disadvantages in solving the two difficulties that exist in real rivers. In this paper, the authors present a method for two‐dimensional shallow water flow calculations to solve both of the main difficulties, by formulating the governing equations in a physical form in terms of physical curvilinear non‐orthogonal co‐ordinates (s,n). Derivation of the governing equations is explained, and two numerical examples are employed to demonstrate that the presented method is applicable to non‐orthogonal and significantly non‐uniform grids. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

A multigrid method with higher‐order discretization schemes

The implementation of the multigrid method into the SIMPLE algorithm presents interesting aspects concerning the mass fluxes conservation on coarser grids, the k–ε turbulence model and the higher‐order discretization schemes. Higher‐order discretization schemes for the convection terms are increasingly used in order to guarantee accuracy in demanding engineering applications. However, when used in single‐grid algorithms, their convergence is considerably slower compared with the first‐order schemes. Unbounded higher‐order schemes offer maximum accuracy, but quite often they do not converge due to their oscillatory behaviour. This paper demonstrates the dual function of the multigrid method: reduction of CPU time and stabilization of the iterating procedure, making it possible to perform computations with the third‐order accurate QUICK scheme in all cases. The method is applied to the calculation of two‐ and three‐dimensional flows with or without turbulence modelling. The results show that the convergence rate of the present algorithm does not deteriorate when QUICK is used and that, if applied on complex engineering cases, large gains in computational time can be achieved. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical prediction of a turbulent curved wake and comparison with experimental data

Numerical studies of the curved wake of a NACA 0012 airfoil of chord length 0.150 m are presented. The airfoil is placed in air at 10 m/s in the straight section of a duct of 0.457 m × 0.457 m cross‐section followed by a 90° bend with a mean radius‐to‐height ratio of 1.17. The trailing edge is located at one chord length upstream of the bend entry plane. The authors' own measurements are used to define the boundary conditions and for comparison with the predicted results. The numerical models are based on the time‐averaged, three‐dimensional conservation equations of fluid flow, incorporating the k–ε, RNG k–ε, realizable k–ε and the Reynolds stress turbulence models. The results show that the models are capable of predicting the effects of curvature on the wake development. However, quantitative differences between prediction and experiment exist. The results obtained using the Reynolds stress model show better agreement with the experimental data, compared with the k–ε based models, but not consistently for all parameters. There are also better predictions by the RNG k–ε and realizable k–ε models compared with the standard k–ε model. The predicted results using the RNG k–ε are closer to experimental data than the realizable k–ε. Copyright © 2005 John Wiley & Sons, Ltd.     

On steady state computation of turbulent flows using k–ε models approximated by the time splitting method

The time splitting method is frequently used in numerical integration of flow equations with source terms since it allows almost independent programming for the source part. In this paper we will consider the question of convergence to steady state of the time splitting method applied to k–ε turbulence models. This analysis is derived from a properly defined scalar study and is carried out with success for the coupled k–ε equations. It is found that the time splitting method does not allow convergence to steady state for any choice of finite values of the time step. Numerical experiments for some typical turbulent compressible flow problems support the fact that the time splitting method is always nonconvergent, while its nonsplitting counterpart is convergent. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical method for calculation of the incompressible flow in general curvilinear co‐ordinates with double staggered grid

A solution methodology has been developed for incompressible flow in general curvilinear co‐ordinates. Two staggered grids are used to discretize the physical domain. The first grid is a MAC quadrilateral mesh with pressure arranged at the centre and the Cartesian velocity components located at the middle of the sides of the mesh. The second grid is so displaced that its corners correspond to the centre of the first grid. In the second grid the pressure is placed at the corner of the first grid. The discretized mass and momentum conservation equations are derived on a control volume. The two pressure grid functions are coupled explicitly through the boundary conditions and implicitly through the velocity of the field. The introduction of these two grid functions avoids an averaging of pressure and velocity components when calculating terms that are generated in general curvilinear co‐ordinates. The SIMPLE calculation procedure is extended to the present curvilinear co‐ordinates with double grids. Application of the methodology is illustrated by calculation of well‐known external and internal problems: viscous flow over a circular cylinder, with Reynolds numbers ranging from 10 to 40, and lid‐driven flow in a cavity with inclined walls are examined. The numerical results are in close agreement with experimental results and other numerical data. Copyright © 2003 John Wiley & Sons, Ltd.     

Depth‐averaged two‐dimensional curvilinear explicit finite analytic model for open‐channel flows

A depth‐averaged two‐dimensional model has been developed in the curvilinear co‐ordinate system for free‐surface flow problems. The non‐linear convective terms of the momentum equations are discretized based on the explicit–finite–analytic method with second‐order accuracy in space and first‐order accuracy in time. The other terms of the momentum equations, as well as the mass conservation equation, are discretized by the finite difference method. The discretized governing equations are solved in turn, and iteration in each time step is adopted to guarantee the numerical convergence. The new model has been applied to various flow situations, even for the cases with the presence of sub‐critical and supercritical flows simultaneously or sequentially. Comparisons between the numerical results and the experimental data show that the proposed model is robust with satisfactory accuracy. Copyright © 2000 John Wiley & Sons, Ltd.     

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.     

A fully coupled Navier–Stokes solver for calculation of turbulent incompressible free surface flow past a ship hull

This paper deals with the calculation of free surface flow of viscous incompressible fluid around the hull of a boat moving with rectilinear motion. An original method used to avoid a large part of the theoretical problems connected with free surface boundary conditions in three‐dimensional Navier–Stokes–Reynolds equations is proposed here. The linearised system of convective equations for velocities, pressure and free surface elevation unknowns is discretised by finite differences and two methods to solve the fully coupled resulting matrix are presented here. The non‐linear convergence of fully coupled algorithm is compared with the velocity–pressure weakly coupled algorithm SIMPLER. Turbulence is taken into account through Reynolds decomposition and k–ε or k–ω model to close the equations. These two models are implemented without wall function and numerical calculations are performed up to the viscous sub‐layer. Numerical results and comparisons with experiments are presented on the Series 60 CB=0.60 ship model for a Reynolds number Rn=4.5×10⁶ and a Froude number Fn=0.316. Copyright © 1999 John Wiley & Sons, Ltd.     

Verification testing in computational fluid dynamics: an example using Reynolds‐averaged Navier–Stokes methods for two‐dimensional flow in the near wake of a circular cylinder

Verification testing was performed for various Reynolds‐averaged Navier–Stokes methods for uniform flow past a circular cylinder at Re= 5232. The standard and renormalized group (RNG) versions of the k–ε method were examined, along with the Boussinesq, Speziale and Launder constitutive relationships. Wind tunnel experiments for flow past a circular cylinder were also performed to obtain a comparative data set. Preliminary studies demonstrate poor convergence for the Speziale relationship. Verification testing with the standard and RNG k–ε models suggests that the simulations exhibit global monotonic convergence for the Boussinesq models. However, the global order of accuracy of the methods was much lower than the expected order of accuracy of 2. For this reason, pointwise convergence ratios and orders of accuracy were computed to show that not all sampling locations had converged (standard k–ε model: 19% failed to converge; RNG k–ε model: 14% failed to converge). When the non‐convergent points were removed from consideration, the average orders of accuracy are closer to the expected value (standard k–ε model: 1.41; RNG k–ε model: 1.27). Poor iterative and global grid convergence was found for the RNG k–ε/Launder model. The standard and RNG k–ε models with the Boussinesq relationship were compared with experimental data and yielded results significantly different from the experiments. Copyright © 2003 John Wiley & Sons, Ltd.     

Incompressible flows with combustion simulated by a preconditioning method using multigrid acceleration and MUSCL reconstruction

A numerical algorithm for the steady state solution of three‐dimensional incompressible flows is presented. A preconditioned time marching scheme is applied to the conservative form of the governing equations. The preconditioning matrix multiplies the time derivatives of the system and circumvents the eigenvalue‐caused stiffness at low speed. The formulation is suitable for constant density flows and for flows where the density depends on non‐passive scalars, such as in low‐speed combustion applications. The k–ε model accounts for turbulent transport effects. A cell‐centred finite volume formulation with a Runge–Kutta time stepping scheme for the primitive variables is used. Second‐order spatial accuracy is achieved by developing for the preconditioned system an approximate Riemann solver with MUSCL reconstruction. A multi‐grid technique coupled with local time stepping and implicit residual smoothing is used to accelerate the convergence to the steady state solution. The convergence behaviour and the validation of the predicted solutions are examined for laminar and turbulent constant density flows and for a turbulent non‐premixed flame simulated by a presumed probability density function (PDF) model. Copyright © 2001 John Wiley & Sons, Ltd.     

Finite element implementation of two‐equation and algebraic stress turbulence models for steady incompressible flows

The main purpose of this paper is to describe a finite element formulation for solving the equations for k and ε of the classical k–ε turbulence model, or any other two‐equation model. The finite element discretization is based on the SUPG method together with a discontinuity capturing technique to deal with sharp internal and boundary layers. The iterative strategy consists of several nested loops, the outermost being the linearization of the Navier–Stokes equations. The basic k–ε model is used for the implementation of an algebraic stress model that is able to account for the effects of rotation. Some numerical examples are presented in order to show the performance of the proposed scheme for simulating directly steady flows, without the need of reaching the steady state through a transient evolution. Copyright © 1999 John Wiley & Sons, Ltd.     

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.     

Experimental validation of two depth‐averaged turbulence models

A finite volume turbulence model for the resolution of the two‐dimensional shallow water equations with turbulent term is presented. After making a finite volume discretization of the depth‐averaged k–ε equations in conservative form, the q–r equations, that give stability to the process, are obtained. Wall and inlet boundary conditions for the turbulent equations and wall conditions for the hydrodynamic equations are discussed. A comparison between the k–ε and q–r models and some experimental results is made. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical study of turbulent flow over two‐dimensional surface‐mounted ribs in a channel

This study accurately predicts the cases of turbulent flow around a surface‐mounted two‐dimensional rib with varying lengths. The numerical method employs a differencing scheme for integrating the elliptic Reynolds‐averaged Navier–Stokes equations and the continuity equation. A two‐equation k–ε turbulence model is employed to simulate the turbulent transport quantities and close the solving problem. The near‐wall regions of the separated sides of the rib are resolved by a near‐wall model of a two‐layer approach instead of the wall function approximation. Computations for flow over a surface‐mounted rectangular rib are conducted for the variations in the rib lengths. Results indicate that upstream of the obstacle, the length of the recirculating region remains unchanged with varying rib lengths; while the downstream length of the recirculating region is a strong function of rib length and changes nearly linearly for the varying lengths of B/H=0.1 to B/H=4.0. Reattachment on top of the rib, owing to its increasing length, affects the downstream boundary layer development. Copyright © 1999 John Wiley & Sons, Ltd.     

On coupling the Reynolds‐averaged Navier–Stokes equations with two‐equation turbulence model equations

Two methods for coupling the Reynolds‐averaged Navier–Stokes equations with the q–ω turbulence model equations on structured grid systems have been studied; namely a loosely coupled method and a strongly coupled method. The loosely coupled method first solves the Navier–Stokes equations with the turbulent viscosity fixed. In a subsequent step, the turbulence model equations are solved with all flow quantities fixed. On the other hand, the strongly coupled method solves the Reynolds‐averaged Navier–Stokes equations and the turbulence model equations simultaneously. In this paper, numerical stabilities of both methods in conjunction with the approximated factorization‐alternative direction implicit method are analysed. The effect of the turbulent kinetic energy terms in the governing equations on the convergence characteristics is also studied. The performance of the two methods is compared for several two‐ and three‐dimensional problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical simulation of the vertical structure of discontinuous flows

A numerical method to solve the Reynolds‐averaged Navier–Stokes equations with the presence of discontinuities is outlined and discussed. The pressure is decomposed into the sum of a hydrostatic component and a hydrodynamic component. The numerical technique is based upon the classical staggered grids and semi‐implicit finite difference methods applied for quasi‐ and non‐hydrostatic flows. The advection terms in the momentum equations are approximated in order to conserve mass and momentum following the principles recently developed for the numerical simulation of shallow water flows with large gradients. Conservation of these properties is the most important aspect to represent near local discontinuities in the solution, following from sharp bottom gradients or hydraulic jumps. The model is applied to reproduce the flow over a step where a hydraulic jump forms downstream. The hydrostatic pressure assumption fails to represent this type of flow mainly because of the pressure deviation from the hydrostatic values downstream the step. Fairly accurate results are obtained from the numerical model compared with experimental data. Deviation from the data is found to be inherent to the standard k–ε model implemented. Copyright © 2001 John Wiley & Sons, Ltd.     

An application of Roe's flux-difference splitting for k-ϵ turbulence model

In this paper Roe's flux-difference splitting is applied for the solution of Reynolds-averaged Navier-Stokes equations. Turbulence is modelled using a low-Reynolds number form of the k-? tubulence model. The coupling between the turbulence kinetic energy equation and the inviscid part of the flow equations is taken into account. The equations are solved with a diagonally dominant alternating direction implicit (DDADI) factorized implicit time integration method. A multigrid algorithm is used to accelerate the convergence. To improve the stability some modifications are needed in comparison with the application of an algebraic turbulence model. The developed method is applied to three different test cases. These cases show the efficiency of the algorithm, but the results are only marginally better than those obtained with algebraic models.