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.     

Flow simulation on moving boundary‐fitted grids and application to fluid–structure interaction problems

We present a method for the parallel numerical simulation of transient three‐dimensional fluid–structure interaction problems. Here, we consider the interaction of incompressible flow in the fluid domain and linear elastic deformation in the solid domain. The coupled problem is tackled by an approach based on the classical alternating Schwarz method with non‐overlapping subdomains, the subproblems are solved alternatingly and the coupling conditions are realized via the exchange of boundary conditions. The elasticity problem is solved by a standard linear finite element method. A main issue is that the flow solver has to be able to handle time‐dependent domains. To this end, we present a technique to solve the incompressible Navier–Stokes equation in three‐dimensional domains with moving boundaries. This numerical method is a generalization of a finite volume discretization using curvilinear coordinates to time‐dependent coordinate transformations. It corresponds to a discretization of the arbitrary Lagrangian–Eulerian formulation of the Navier–Stokes equations. Here the grid velocity is treated in such a way that the so‐called Geometric Conservation Law is implicitly satisfied. Altogether, our approach results in a scheme which is an extension of the well‐known MAC‐method to a staggered mesh in moving boundary‐fitted coordinates which uses grid‐dependent velocity components as the primary variables. To validate our method, we present some numerical results which show that second‐order convergence in space is obtained on moving grids. Finally, we give the results of a fully coupled fluid–structure interaction problem. It turns out that already a simple explicit coupling with one iteration of the Schwarz method, i.e. one solution of the fluid problem and one solution of the elasticity problem per time step, yields a convergent, simple, yet efficient overall method for fluid–structure interaction problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Spectral/hp penalty least‐squares finite element formulation for unsteady incompressible flows

In this paper, we present spectral/hp penalty least‐squares finite element formulation for the numerical solution of unsteady incompressible Navier–Stokes equations. Pressure is eliminated from Navier–Stokes equations using penalty method, and finite element model is developed in terms of velocity, vorticity and dilatation. High‐order element expansions are used to construct discrete form. Unlike other penalty finite element formulations, equal‐order Gauss integration is used for both viscous and penalty terms of the coefficient matrix. For time integration, space–time decoupled schemes are implemented. Second‐order accuracy of the time integration scheme is established using the method of manufactured solution. Numerical results are presented for impulsively started lid‐driven cavity flow at Reynolds number of 5000 and transient flow over a backward‐facing step. The effect of penalty parameter on the accuracy is investigated thoroughly in this paper and results are presented for a range of penalty parameter. Present formulation produces very accurate results for even very low penalty parameters (10–50). Copyright © 2008 John Wiley & Sons, Ltd.     

An arbitrary Lagrangian–Eulerian finite element approach to non‐steady state turbulent fluid flow with application to mould filling in casting

This paper presents a two‐dimensional Lagrangian–Eulerian finite element approach of non‐steady state turbulent fluid flows with free surfaces. The proposed model is based on a velocity–pressure finite element Navier–Stokes solver, including an augmented Lagrangian technique and an iterative resolution of Uzawa type. Turbulent effects are taken into account with the k–ε two‐equation statistical model. Mesh updating is carried out through an arbitrary Lagrangian–Eulerian (ALE) method in order to describe properly the free surface evolution. Three comparisons between experimental and numerical results illustrate the efficiency of the method. The first one is turbulent flow in an academic geometry, the second one is a mould filling in effective casting conditions and the third one is a precise confrontation to a water model. Copyright © 2000 John Wiley & Sons, Ltd.     

A high‐order finite difference method for incompressible fluid turbulence simulations

A Hermitian–Fourier numerical method for solving the Navier–Stokes equations with one non‐homogeneous direction had been presented by Schiestel and Viazzo (Internat. J. Comput. Fluids 1995; 24 (6):739). In the present paper, an extension of the method is devised for solving problems with two non‐homogeneous directions. This extension is indeed not trivial since new algorithms will be necessary, in particular for pressure calculation. The method uses Hermitian finite differences in the non‐periodic directions whereas Fourier pseudo‐spectral developments are used in the remaining periodic direction. Pressure–velocity coupling is solved by a simplified Poisson equation for the pressure correction using direct method of solution that preserves Hermitian accuracy for pressure. The turbulent flow after a backward facing step has been used as a test case to show the capabilities of the method. The applications in view are mainly concerning the numerical simulation of turbulent and transitional flows. Copyright © 2003 John Wiley & Sons, Ltd.     

Large eddy simulation of incompressible free round jet with discontinuous Galerkin method

In the paper, discontinuous Galerkin method is applied to simulation of incompressible free round turbulent jet using large eddy simulation with eddy viscosity approach. The solution algorithm is based on the classical projection method, but instead of the solution of the Poisson equation, a parabolic equation is advanced in pseudo‐time, which provides the pressure field ensuring the proper pressure–velocity coupling. For time and pseudo‐time integration, explicit Runge–Kutta method is employed. The computational meshes consist of hexahedral elements with flat faces. Within a given finite element, all flow variables are expressed with modal expansions of the same order (including velocity and pressure). Discretisation of the viscous terms in the Navier–Stokes equations and Laplacian in the Poisson equation is stabilised with mixed finite element approach. The correctness of the solution algorithm is verified in a commonly used test case of laminar flow in 3D lid‐driven cavity. The results of computations of the free jet are compared with experimental and numerical reference data, the latter obtained from the high‐order pseudospectral code. The statistics of centerline flow velocity – mean velocity and its fluctuations – show satisfactory agreement with the reference data. Copyright © 2015 John Wiley & Sons, Ltd.     

Simulating turbulent Dean flow in Cartesian coordinates

A simplified approach to simulate turbulent flows in curved channels is proposed. A set of governing equations of motion in Cartesian coordinates is derived from the full Navier–Stokes equations in cylindrical coordinates. Terms to first order in the dimensionless curvature parameter are retained, whereas higher‐order terms are neglected. The curvature terms are implemented in a conventional Navier–Stokes code using Cartesian coordinates. Direct numerical simulations (DNS) of turbulent flow in weakly curved channels are performed. The pronounced asymmetries in the mean flow and the turbulence statistics observed in earlier DNS studies are faithfully reproduced by the present simplified Navier–Stokes model. It is particularly rewarding that also distinct pairs of counter‐rotating streamwise‐oriented vortices are embedded in the simulated flow field. Copyright © 2008 John Wiley & Sons, Ltd.     

A reduced‐order approach for optimal control of fluids using proper orthogonal decomposition

In this article, a reduced‐order modeling approach, suitable for active control of fluid dynamical systems, based on proper orthogonal decomposition (POD) is presented. The rationale behind the reduced‐order modeling is that numerical simulation of Navier–Stokes equations is still too costly for the purpose of optimization and control of unsteady flows. The possibility of obtaining reduced‐order models that reduce the computational complexity associated with the Navier–Stokes equations is examined while capturing the essential dynamics by using the POD. The POD allows the extraction of a reduced set of basis functions, perhaps just a few, from a computational or experimental database through an eigenvalue analysis. The solution is then obtained as a linear combination of this reduced set of basis functions by means of Galerkin projection. This makes it attractive for optimal control and estimation of systems governed by partial differential equations (PDEs). It is used here in active control of fluid flows governed by the Navier–Stokes equations. In particular, flow over a backward‐facing step is considered. Reduced‐order models/low‐dimensional dynamical models for this system are obtained using POD basis functions (global) from the finite element discretizations of the Navier–Stokes equations. Their effectiveness in flow control applications is shown on a recirculation control problem using blowing on the channel boundary. Implementational issues are discussed and numerical experiments are presented. Copyright © 2000 John Wiley & Sons, Ltd.     

Finite element modified method of characteristics for the Navier–Stokes equations

An algorithm based on the finite element modified method of characteristics (FEMMC) is presented to solve convection–diffusion, Burgers and unsteady incompressible Navier–Stokes equations for laminar flow. Solutions for these progressively more involved problems are presented so as to give numerical evidence for the robustness, good error characteristics and accuracy of our method. To solve the Navier–Stokes equations, an approach that can be conceived as a fractional step method is used. The innovative first stage of our method is a backward search and interpolation at the foot of the characteristics, which we identify as the convective step. In this particular work, this step is followed by a conjugate gradient solution of the remaining Stokes problem. Numerical results are presented for:

• a Convection–diffusion equation. Gaussian hill in a uniform rotating field.
• b Burgers equations with viscosity.
• c Navier–Stokes solution of lid‐driven cavity flow at relatively high Reynolds numbers.
• d Navier–Stokes solution of flow around a circular cylinder at Re=100.

Copyright © 2000 John Wiley & Sons, Ltd.     

A robust numerical method for flow through a pipe driven by an oscillating pressure gradient

The problem of periodic flow of an incompressible fluid through a pipe, which is driven by an oscillating pressure gradient (e.g. a reciprocating piston), is investigated in the case of a large Reynolds number. This process is described by a singularly perturbed parabolic equation with a periodic right‐hand side, where the singular perturbation parameter is the viscosity ν. The periodic solution of this problem is a solution of the Navier–Stokes equations with cylindrical symmetry. We are interested in constructing a parameter‐robust numerical method for this problem, i.e. a numerical method generating numerical approximations that converge uniformly with respect to the parameter ν and require a bounded time, independent of the value of ν, for their computation. Our method comprises a standard monotone discretization of the problem on non‐standard piecewise uniform meshes condensing in a neighbourhood of the boundary layer. The transition point between segments of the mesh with different step sizes is chosen in accordance with the behaviour of the analytic solution in the boundary layer region. In this paper we construct the numerical method and discuss the results of extensive numerical experiments, which show experimentally that the method is parameter‐robust. Copyright © 2006 John Wiley & Sons, Ltd.     

Accurate iteration-free mixed-stabilised formulation for laminar incompressible Navier–Stokes: Applications to fluid–structure interaction

Stabilised mixed velocity–pressure formulations are one of the widely-used finite element schemes for computing the numerical solutions of laminar incompressible Navier–Stokes. In these formulations, the Newton–Raphson scheme is employed to solve the nonlinearity in the convection term. One fundamental issue with this approach is the computational cost incurred in the Newton–Raphson iterations at every load/time step. In this paper, we present an iteration-free mixed finite element formulation for incompressible Navier–Stokes that preserves second-order temporal accuracy of the generalised-alpha and related schemes for both velocity and pressure fields. First, we demonstrate the second-order temporal accuracy using numerical convergence studies for an example with a manufactured solution. Later, we assess the accuracy and the computational benefits of the proposed scheme by studying the benchmark example of flow past a fixed circular cylinder. Towards showcasing the applicability of the proposed technique in a wider context, the inf–sup stable P2–P1 pair for the formulation without stabilisation is also considered. Finally, the resulting benefits of using the proposed scheme for fluid–structure interaction problems are illustrated using two benchmark examples in fluid-flexible structure interaction.     

Numerical simulation of incompressible two‐phase flows with a Boussinesq–Scriven interface stress tensor

We consider the numerical simulation of a three‐dimensional two‐phase incompressible flow with a viscous interface. The simulation is based on a sharp interface Navier–Stokes model and the Boussinesq–Scriven constitutive law for the interface viscous stress tensor. In the recent paper [Soft Matter 7, 7797–7804, 2011], a model problem with a spherical droplet in a Stokes Poiseuille flow with a Boussinesq–Scriven law for the surface viscosity has been analyzed. In that paper, relations for the droplet migration velocity are derived. We relate the results obtained with our numerical solver for the two‐phase Navier–Stokes model to these theoretical relations. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical simulation of turbulent flows in trapezoidal meandering compound open channels

Three‐dimensional (3D) numerical study is presented to investigate the turbulent flow in meandering compound open channels with trapezoidal cross‐sections. The flow simulation is carried out by solving the 3D Reynolds‐averaged continuity and Navier–Stokes equations with Reynolds stress equation model (RSM) for steady‐state flow. Finite volume method (FVM) is applied to numerically solve the governing equations of fluid flow. The velocity magnitude, tangential velocity, transverse velocity and Reynolds stresses are calculated for various flow conditions. Good agreement between the simulated and available laboratory measurements was obtained, indicating that the RSM can accurately predict the complicated flow phenomenon. Comparison of the calculated secondary currents of four cases (one being inbank flow and other three being overbank flow) with different water depths reveals that (i) the inbank flow exhibits different flow behaviors from that of the overbank flow does and (ii) the water depth has significant effects on the magnitude and direction of secondary currents. Copyright © 2010 John Wiley & Sons, Ltd.     

Finite element solution of three‐dimensional turbulent flows applied to mold‐filling problems

This paper presents a finite element solution algorithm for three‐dimensional isothermal turbulent flows for mold‐filling applications. The problems of interest present unusual challenges for both the physical modelling and the solution algorithm. High‐Reynolds number transient turbulent flows with free surfaces have to be computed on complex three‐dimensional geometries. In this work, a segregated algorithm is used to solve the Navier–Stokes, turbulence and front‐tracking equations. The streamline–upwind/Petrov–Galerkin method is used to obtain stable solutions to convection‐dominated problems. Turbulence is modelled using either a one‐equation turbulence model or the κ–ε two‐equation model with wall functions. Turbulence equations are solved for the natural logarithm of the turbulence variables. The change of dependent variables allows for a robust solution algorithm and good predictions even on coarse meshes. This is very important in the case of large three‐dimensional applications for which highly refined meshes result in untreatable large numbers of elements. The position of the flow front in the mold cavity is computed using a level set approach. Finally, equations are integrated in time using an implicit Euler scheme. The methodology presents the robustness and cost effectiveness needed to tackle complex industrial applications. Copyright © 2000 John Wiley & Sons, Ltd.     

A priori tests on numerical errors in large eddy simulation using finite differences and explicit filtering

When low‐order finite‐difference methods are applied in large eddy simulation (LES), the magnitude of the numerical error may be larger than that of the subgrid‐scale (SGS) term. In this paper, the effect of explicit filtering on the numerical error related to the spatial discretization of the convection term and the exact SGS term is studied a priori in the turbulent fully developed channel flow. As the filter width is increased the grid resolution is kept constant. Also filtering in the inhomogeneous wall‐normal direction is discussed. The main conclusions are related to two approaches to explicit filtering. In the traditional approach, the whole velocity field is filtered explicitly while in the alternative approach, only the non‐linear convection term of the Navier–Stokes equations is filtered explicitly. Based on the results presented in the paper it seems that the first approach leads to an unphysical situation. However, the later approach works in the desired way, and the numerical error becomes clearly smaller than the SGS term. The main difference between the two approaches seems to be the interpretation of the resolved non‐linear term in the filtered Navier–Stokes equations. Copyright © 2005 John Wiley & Sons, Ltd.     

Calculation of three‐dimensional turbulent flow with a finite volume multigrid method

A numerical method for the efficient calculation of three‐dimensional incompressible turbulent flow in curvilinear co‐ordinates is presented. The mathematical model consists of the Reynolds averaged Navier–Stokes equations and the k–ε turbulence model. The numerical method is based on the SIMPLE pressure‐correction algorithm with finite volume discretization in curvilinear co‐ordinates. To accelerate the convergence of the solution method a full approximation scheme‐full multigrid (FAS‐FMG) method is utilized. The solution of the k–ε transport equations is embedded in the multigrid iteration. The improved convergence characteristic of the multigrid method is demonstrated by means of several calculations of three‐dimensional flow cases. Copyright © 1999 John Wiley & Sons, Ltd.     

Numerical solution of single and multi‐phase internal transonic flow problems

This paper presents numerical methods for solving turbulent and two‐phase transonic flow problems. The Navier–Stokes equations are solved using cell‐vertex Lax–Wendroff method with artificial dissipation or cell‐centred upwind method with Roe's Riemann solver and linear reconstruction. Due to a big difference of time scales in two‐phase flow of condensing steam a fractional step method is used. Test cases including 2D condensing flow in a nozzle and one‐phase transonic flow in a turbine cascade with transition to turbulence are presented. Copyright © 2005 John Wiley & Sons, Ltd.     

A finite element model for free surface flows on fixed meshes

In this paper, we present a finite element model for free surface flows on fixed meshes. The main novelty of the approach, compared with typical fixed mesh finite element models for such flows, is that we take advantage of the particularities of free surface flow, instead of considering it a particular case of two‐phase flow. The fact that a given free surface implies a known boundary condition on the interface, allows us to solve the Navier–Stokes equations on the fluid domain uncoupled from the solution on the rest of the finite element mesh. This, together with the use of enhanced integration allows us to model low Froude number flows accurately, something that is not possible with typical two‐phase flow models applied to free surface flow. Copyright © 2007 John Wiley & Sons, Ltd.     

Effect of duct velocity profile and buoyancy‐induced flow on efficiency of transient hydrodynamic removal of a contaminant from a cavity

This paper describes a preliminary numerical analysis of the effect of duct velocity profile and buoyancy‐induced flow generated by the heat source on hydrodynamic removal of contaminants contained in cavities. The process of fluid renewal in a cavity is modelled via a numerical solution of the Navier–Stokes equations coupled with the energy equation for transient flows. The foulant has the same density as the fluid in the duct and the duct velocity profile is considered to be Poiseuille flow and Couette flow, respectively. The results show that the change in Grashof number and duct flow velocity profile causes a dramatic difference in the observed flow patterns and cleaning efficiency. From a cleaning perspective, the results suggest that Couette flow at higher value of Grashof number becomes more effective in further purging of contaminated fluid from a cavity. Copyright © 2004 John Wiley & Sons, Ltd.     

A projection method for the spectral solution of non‐homogeneous and incompressible Navier–Stokes equations

This paper is devoted to the development of a parallel, spectral and second‐order time‐accurate method for solving the incompressible and variable density Navier–Stokes equations. The method is well suited for finite thickness density layers and is very efficient, especially for three‐dimensional computations. It is based on an exact projection technique. To enforce incompressibility, for a non‐homogeneous fluid, the pressure is computed using an iterative algorithm. A complete study of the convergence properties of this algorithm is done for different density variations. Numerical simulations showing, qualitatively, the capabilities of the developed Navier–Stokes solver for many realistic problems are presented. The numerical procedure is also validated quantitatively by reproducing growth rates from the linear instability theory in a three‐dimensional direct numerical simulation of an unstable, non‐homogeneous, flow configuration. It is also shown that, even in a turbulent flow, the spectral accuracy is recovered. Copyright © 2012 John Wiley & Sons, Ltd.