An algorithm for rotating axisymmetric flows: model,validation and industrial applications

The study of axisymmetric flows is of interest not only from an academic point of view, due to the existence of exact solutions of Navier–Stokes equations, but also from an industrial point of view, since these kind of flows are frequently found in several applications. In the present work the development and implementation of a finite element algorithm to solve Navier–Stokes equations with axisymmetric geometry and boundary conditions is presented. Such algorithm allows the simulation of flows with tangential velocity, including free surface flows, for both laminar and turbulent conditions. Pseudo‐concentration technique is used to model the free surface (or the interface between two fluids) and the k–ε model is employed to take into account turbulent effects. The finite element model is validated by comparisons with analytical solutions of Navier–Stokes equations and experimental measurements. Two different industrial applications are presented. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical investigation of fully developed turbulent fluid flow and heat transfer in a square duct

Three-dimensional fully developed turbulent fluid flow and heat transfer in a square duct are numerically investigated with the author's anisotropic low-Reynolds-number k-ε turbulence model. Special attenton has been given to the regions close to the wall and the corner, which are known to influence the characteristics of secondary flow a great deal. Hence, instead of the common wall function approach, the no-slip boundary condition at the wall is directly used. Velocity and temperature profiles are predicted for fully developed turbulent flows with constant wall temperature. The predicted variations of both local wall shear stress and local wall heat flux are shown to be in close agreement with available experimental data. The present paper also presents the budget of turbulent kinetic energy equation and the systematic evaluation for existing wall function forms. The commonly adopted wall function forms that are valid for two-dimensional flows are found to be inadequate for three-dimensional turbulent flows in a square duct.     

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.     

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.     

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.     

An implicit mixed finite-volume–finite-element method for solving 3D turbulent compressible flows

The development of new aeronautic projects require accurate and efficient simulations of compressible flows in complex geometries. It is well known that most flows of interest are at least locally turbulent and that the modelling of this turbulence is critical for the reliability of the computations. A turbulence closure model which is both cheap and reasonably accurate is an essential part of a compressible code. An implicit algorithm to solve the 2D and 3D compressible Navier–Stokes equations on unstructured triangular/tetrahedral grids has been extended to turbulent flows. This numerical scheme is based on second-order finite element–finite volume discretization: the diffusive and source terms of the Navier–Stokes equations are computed using a finite element method, while the other terms are computed with a finite volume method. Finite volume cells are built around each node by means of the medians. The convective fluxes are evaluated with the approximate Riemann solver of Roe coupled with the van Albada limiter. The standard k–ϵ model has been introduced to take into account turbulence. Implicit integration schemes with efficient numerical methods (CGS, GMRES and various preconditioning techniques) have also been implemented. Our interest is to present the whole method and to demonstrate its limitations on some well-known test cases in three-dimensional geometries. © 1997 John Wiley & Sons, Ltd.     

Parallel large eddy simulation of turbulent flow around MIRA model using linear equal‐order finite element method

A parallel large eddy simulation code that adopts domain decomposition method has been developed for large‐scale computation of turbulent flows around an arbitrarily shaped body. For the temporal integration of the unsteady incompressible Navier–Stokes equation, fractional 4‐step splitting algorithm is adopted, and for the modelling of small eddies in turbulent flows, the Smagorinsky model is used. For the parallelization of the code, METIS and Message Passing Interface Libraries are used, respectively, to partition the computational domain and to communicate data between processors. To validate the parallel architecture and to estimate its performance, a three‐dimensional laminar driven cavity flow inside a cubical enclosure has been solved. To validate the turbulence calculation, the turbulent channel flows at Re_τ = 180 and 1050 are simulated and compared with previous results. Then, a backward facing step flow is solved and compared with a DNS result for overall code validation. Finally, the turbulent flow around MIRA model at Re = 2.6 × 10⁶ is simulated by using approximately 6.7 million nodes. Scalability curve obtained from this simulation shows that scalable results are obtained. The calculated drag coefficient agrees better with the experimental result than those previously obtained by using two‐equation turbulence models. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical modelling of 3D stratified free surface flows: a case study of sediment dumping

A three‐dimensional numerical model has been developed to simulate stratified flows with free surfaces. The model is based on the Reynolds‐averaged Navier–Stokes (RANS) equations with variable fluid density. The equations are solved in a transformed σ‐coordinate system with the use of operator‐splitting method (Int. J. Numer. Meth. Fluids 2002; 38 :1045–1068). The numerical model is validated against the one‐dimensional diffusion problem and the two‐dimensional density‐gradient flow. Excellent agreements are obtained between numerical results and analytical solutions. The model is then used to study transport phenomena of dumped sediments into a water body, which has been modelled as a strongly stratified flow. For the two‐dimensional problem, the numerical results compare well with experimental data in terms of mean particle falling velocity and spreading rate of the sediment cloud for both coarse and medium‐size sediments. The model is also employed to study the dumping of sediments in a three‐dimensional environment with the presence of free surface. It is found that during the descending process an annulus‐like cloud is formed for fine sediments whereas a plate‐like cloud for medium‐size sediments. The model is proven to be a good tool to simulate strongly stratified free surface flows. Copyright © 2005 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.     

A Two-Time-Scale Model for Turbulent Mixing Flows Induced by Rayleigh–Taylor and Richtmyer–Meshkov Instabilities

A two-time-scale closure model for compressible flows previously developed is extended to turbulent Rayleigh-Taylor and Richtmyer-Meshkov driven flows where mixing coexists with mean pressure gradients. Two model coefficients are calibrated with the help of Canuto-Goldman's model. For several Rayleigh-Taylor configurations, it is shown that the characteristic lengths scale as t ² while the kinetic energies and spectral transfers behave as t ² and t, respectively. The computed phenomenological coefficients of Youngs' scaling law are compared with experimental data ones. Comparisons with Youngs' three-dimensional numerical simulation (The Physics of Fluids A 3 (1991) 1312) are also performed. Finally three shock tube experiments, where the Richtmyer-Meshkov instability initiates the mixing, are simulated. The mixing thickness evolution is well reproduced while the turbulence levels seem to be overestimated with such first order models. The capability of the two-time-sale model to recover available data for different turbulent flows allows us to conclude to a more universal behavior in comparison with single-time-scale models. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

Shock-unsteadiness model applied to oblique shock wave/turbulent boundary-layer interaction

Reynolds-averaged Navier–Stokes prediction of shock wave/turbulent boundary layer interactions can yield significant error in terms of the size of the separation bubble. In many applications, this can alter the shock structure and the resulting surface properties. Shock-unsteadiness modification of Sinha et al. (Physics of Fluids, Vol.15, No.8, 2003) has shown potential in improving separation bubble prediction in compression corner flows. In this article, the modification is applied to oblique shock wave interacting with a turbulent boundary layer. The challenges involved in the implementation of the shock-unsteadiness correction in the presence of multiple shock waves and expansion fans are addressed in detail. The results show that a robust implementation of the model yields appreciable improvement over standard k–ω turbulence model predictions.     

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 modeling of cohesive sediments dynamics in estuaries: Part I—Description of the model and simulations in the Po River Estuary

The present work contributes to the numerical modeling of complex turbulent multiphasic fluid flows occurring in estuarine channels. This research finds its motivation in the increasing need for efficient management of estuaries by taking into account the complex turbulent stratified flows encountered in estuaries and costal zones. A time‐dependent, 3D finite element model of suspended sediment transport taking into account the effects of cohesiveness between sediments is presented. The model estuary is the forced time‐dependent winds, time elevation at open boundaries and river discharge. To cope with the stiffness problems a decoupling method is employed to solve the shallow‐water equations of mass conservation, momentum and suspended sediment transport with the conventional hydrostatic pressure. The decoupling method partitions a time step into three subcycles according to the physical phenomena. In the first sub‐cycle the pure hydrodynamics including the k–ε turbulence model is solved, followed by the advection–diffusion equations for pollutants (salinity, temperature, suspended sediment concentration, (SSC)), and finally the bed evolution is solved. The model uses a mass‐preserving method based on the so‐called Raviart–Thomas finite element on the unstructured mesh in the horizontal plane, while the multi‐layers system is adopted in vertical with the conventional conforming finite element method, with the advantage that the lowermost and uppermost layers of variable height allow a faithful representation of the time‐varying bed and free surface, respectively. The model has been applied to investigate the SSC and seabed evolution in Po River Estuary (PRE) in Italy. The computed results mimic the field data well. Copyright © 2007 John Wiley & Sons, Ltd.     

Finite element simulation of turbulent Couette–Poiseuille flows using a low Reynolds number k–ϵ model

Developing Couette–Poiseuille flows at Re=5000 are studied using a low Reynolds number k–ϵ two‐equation model and a finite element formulation. Mesh‐independent solutions are obtained using a standard Galerkin formulation and a Galerkin/least‐squares stabilized method. The predictions for the velocity and turbulent kinetic energy are compared with available experimental results and to the DNS data. Second moment closure's solutions are also compared with those of the k–ϵ model. The deficiency of eddy viscosity models to predict dissymmetric low Reynolds number channel flows has been demonstrated. Copyright © 1999 John Wiley & Sons, Ltd.     

A numerical investigation of 2d,steady free surface flows

Systematic tests have been performed to study the behaviour of a numerical method developed to calculate 2D, steady free surface flows. The Reynolds equations are solved in the physical space by employing a non–orthogonal staggered grid, while the k-ε model is adopted to approximate the Reynolds stresses. The free surface is calculated following an iterative procedure and various parameters that affect convergence and accuracy of the numerical solution have been examined. Calculated results are compared with measured data for two cases, i.e. the wave generation above a bottom topography at various Froude numbers and the free surface formation above a submerged hydrofoil. © 1997 John Wiley & Sons, Ltd.     

Stereoscopic PIV: validation and application to an isotropic turbulent flow

 A new stereoscopic PIV system to measure the three velocity components is developed and applied to grid turbulence flows. This system uses two CCD cameras coupled with an accurate cross-correlation calculation method. An experimental test (based upon three-dimensional displacements) has been carried out to demonstrate the capability of this process to locate the maximum of correlation, and to detect accurately the 3D displacements. Experiments in a well-established turbulent flow have validated the method for quantitative measurements and a comparison with LDV results showed a good agreement in terms of mean and fluctuating velocities. Combined PIV and stereoscopic PIV measurements on a turbulent flow revealed the need to the stereoscopic systems to measure accurate 2D velocity fields. It has been shown that an error of up to 10% in the velocity fluctuation measured by conventional PIV could be attained due to 3D effects in highly turbulent cases. Finally, the digital cross-correlation technique adapted to the determination of small displacements seems to be the most suitable technique for stereoscopic PIV. Received: 22 July 1997/Accepted: 27 January 1998     

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

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

An unsteady incompressible Navier–Stokes solver for large eddy simulation of turbulent flows

An unsteady incompressible Navier–Stokes solver that uses a dual time stepping method combined with spatially high‐order‐accurate finite differences, is developed for large eddy simulation (LES) of turbulent flows. The present solver uses a primitive variable formulation that is based on the artificial compressibility method and various convergence–acceleration techniques are incorporated to efficiently simulate unsteady flows. A localized dynamic subgrid model, which is formulated using the subgrid kinetic energy, is employed for subgrid turbulence modeling. To evaluate the accuracy and the efficiency of the new solver, a posteriori tests for various turbulent flows are carried out and the resulting turbulence statistics are compared with existing experimental and direct numerical simulation (DNS) data. Copyright © 1999 John Wiley & Sons, Ltd.     

The calculation of incompressible separated turbulent boundary layers

The algebraic turbulent model of Baldwin and Lomax was incorporated into the incompressible full Navier–Stokes code FIDAP. This model was extensively tested in the past in finite difference codes. We believe that the incorporation of the model also into the finite element code has resulted in a practical method to compute a variety of separated turbulent 2D flows. Firstly, we use the model to compute the attached flow about an aerofoil. Next, the application of the model to separated flows is presented by computing the flows at high angles of attack up to maximum lift. It is shown that the model is capable of predicting separation, steady stall and C_Lmax. As a difficult test of the model we compute the laminar separation bubble development directly using the full Navier–Stokes finite element code. As far as we know, this approach has not yet been reported. The importance of using an appropriate upwinding is discussed. When possible, comparison of computed results with experiments is presented and the agreement is good.