A decoupling numerical method for fluid flow

A first-order non-conforming numerical methodology, Separation method, for fluid flow problems with a 3-point exponential interpolation scheme has been developed. The flow problem is decoupled into multiple one-dimensional subproblems and assembled to form the solutions. A fully staggered grid and a conservational domain centred at the node of interest make the decoupling scheme first-order-accurate. The discretization of each one-dimensional subproblem is based on a 3-point interpolation function and a conservational domain centred at the node of interest. The proposed scheme gives a guaranteed first-order accuracy. It is shown that the traditional upwind (or exponentially weighted upstream) scheme is less than first-order-accurate. The pressure is decoupled from the velocity field using the pressure correction method of SIMPLE. Thomas algorithm (tri-diagonal solver) is used to solve the algebraic equations iteratively. The numerical advantage of the proposed scheme is tested for laminar fluid flows in a torus and in a square-driven cavity. The convergence rates are compared with the traditional schemes for the square-driven cavity problem. Good behaviour of the proposed scheme is ascertained.     

A NOTE ON THE APPLICATION OF THE AVERAGE CORRECTION TECHNIQUE

This note develops an average correction technique for accelerating the rate of convergence of the SIMPLE-like algorithm by implementing the average pressure correction method as proposed by Wen and Ingham (Int. j. numer. methods fluids, 17 , 385–400 (1993); 19 , 889–903 (1994)) with an average velocity correction. The technique is illustrated by considering the classical problem of fluid flow over a backward-facing step using (i) no average correction, (ii) an average velocity correction, (iii) an average pressure correction and (iv) both average velocity and pressure corrections. When both average velocity and pressure corrections are employed, it is found that the number of iterations required for convergence is almost independent of the initial guessed values of fluid velocity and pressure and the fastest rate of convergence may be achieved.     

Different modelings of cell‐face velocities and their effects on the pressure–velocity coupling,accuracy and convergence of solution

In this article, a detailed study on the effects of different modelings of cell‐face velocities on pressure–velocity coupling, accuracy and convergence rate of the solution has been conducted. Discussions are focused on the collocated scheme of Schneider and Karimian (Computational Mechanics 1994; 14 : 1–16) in the context of a control‐volume finite‐element Method. In this scheme, variables at the control volume surface are evaluated based on the physics of their governing equations, and the fully coupled system obtained is solved using a direct sparse solver. A special test problem has been defined to check the pressure–velocity coupling for all of the formulations. Other test cases, including Taylor problem, inviscid converging–diverging nozzle and the lid‐driven cavity, have been conducted for different Reynolds numbers, mesh sizes and time steps to investigate the accuracy and the performance of the formulations. Finally, a reliable and efficient scheme for the evaluation of cell‐face velocities is proposed, which can be easily extended to three dimensions. Copyright © 2010 John Wiley & Sons, Ltd.     

An algebraic factorisation scheme for spectral element solution of incompressible flow and scalar transport

A spectral element algorithm for solution of the unsteady incompressible Navier–Stokes and scalar (species/heat) transport equations is developed using the algebraic factorisation scheme. The new algorithm utilises Nth order Gauss–Lobatto–Legendre points for velocity and the scalar, while (N-2)th order Gauss–Legendre points are used for pressure. As a result, the algorithm does not require inter-element continuity for pressure and pressure boundary conditions on solid surfaces. Implementations of the algorithm are performed for conforming and non-conforming grids. The latter is accomplished using both the point-wise matching and integral projection methods, and applied for grids with both polynomial and geometric non-conformities. Code validation cases include the unsteady scalar convection equation, and Kovasznay flow in two- and three-dimensional domains. Using cases with analytical solutions, the algorithm is shown to achieve spectral accuracy in space and second-order accuracy in time. The results for the Boussinesq approximation for buoyancy-driven flows, and the species mixing in a continuous flow micro-mixer are also included as examples of applications that require long-time integration of the scalar transport equations.     

A least-squares finite element method for time-dependent incompressible flows with thermal convection

The time-dependent Navier–Stokes equations and the energy balance equation for an incompressible, constant property fluid in the Boussinesq approximation are solved by a least-squares finite element method based on a velocity–pressure–vorticity–temperature–heat-flux ( u –P–ω–T– q ) formulation discretized by backward finite differencing in time. The discretization scheme leads to the minimization of the residual in the l²-norm for each time step. Isoparametric bilinear quadrilateral elements and reduced integration are employed. Three examples, thermally driven cavity flow at Rayleigh numbers up to 10⁶, lid-driven cavity flow at Reynolds numbers up to 10⁴ and flow over a square obstacle at Reynolds number 200, are presented to validate the method.     

Three‐dimensional transient Navier–Stokes solvers in cylindrical coordinate system based on a spectral collocation method using explicit treatment of the pressure

A spectral collocation method is developed for solving the three‐dimensional transient Navier–Stokes equations in cylindrical coordinate system. The Chebyshev–Fourier spectral collocation method is used for spatial approximation. A second‐order semi‐implicit scheme with explicit treatment of the pressure and implicit treatment of the viscous term is used for the time discretization. The pressure Poisson equation enforces the incompressibility constraint for the velocity field, and the pressure is solved through the pressure Poisson equation with a Neumann boundary condition. We demonstrate by numerical results that this scheme is stable under the standard Courant–Friedrichs–Lewy (CFL) condition, and is second‐order accurate in time for the velocity, pressure, and divergence. Further, we develop three accurate, stable, and efficient solvers based on this algorithm by selecting different collocation points in r‐, ? ‐, and z‐directions. Additionally, we compare two sets of collocation points used to avoid the axis, and the numerical results indicate that using the Chebyshev Gauss–Radau points in radial direction to avoid the axis is more practical for solving our problem, and its main advantage is to save the CPU time compared with using the Chebyshev Gauss–Lobatto points in radial direction to avoid the axis. Copyright © 2010 John Wiley & Sons, Ltd.     

Preconditioned conjugate gradients for solving the transient Boussinesq equations in three-dimensional geometries

In this paper we present a new version of the ‘modified finite element method’ (MFEM) presented by Gresho, Chan, Lee and Upson.¹ The main modification of the original algorithm is the introduction of a cost-effective and memory-saving iterative solver for the discretized Poisson equation for the pressure. The vectorization of the preconditioner has been especially considered. For low Prandtl number problems we also split the advection-diffusion operator of the energy equation into explicit and implicit parts. In that sense the present approach is related to the recent implicitization of the diffusive terms introduced by Gresho and Chan² and by Gresho.³ The algorithm is applied to the study of buoyancy-driven flow oscillations occuring in a horizontal crucible of molten metal under the action of a horizontal temperature gradient.     

Pressure boundary condition for the time‐dependent incompressible Navier–Stokes equations

In Gresho and Sani (Int. J. Numer. Methods Fluids 1987; 7 :1111–1145; Incompressible Flow and the Finite Element Method, vol. 2. Wiley: New York, 2000) was proposed an important hypothesis regarding the pressure Poisson equation (PPE) for incompressible flow: Stated there but not proven was a so‐called equivalence theorem (assertion) that stated/asserted that if the Navier–Stokes momentum equation is solved simultaneously with the PPE whose boundary condition (BC) is the Neumann condition obtained by applying the normal component of the momentum equation on the boundary on which the normal component of velocity is specified as a Dirichlet BC, the solution (u, p) would be exactly the same as if the ‘primitive’ equations, in which the PPE plus Neumann BC is replaced by the usual divergence‐free constraint (? · u = 0), were solved instead. This issue is explored in sufficient detail in this paper so as to actually prove the theorem for at least some situations. Additionally, like the original/primitive equations that require no BC for the pressure, the new results establish the same thing when the PPE approach is employed. Copyright © 2005 John Wiley & Sons, Ltd.     

A collocated grid,projection method for time‐accurate calculation of low‐Mach number variable density flows in general curvilinear coordinates

A time‐accurate algorithm is proposed for low‐Mach number, variable density flows on curvilinear grids. Spatial discretization is performed on collocated grid that offers computational simplicity in curvilinear coordinates. The flux interpolation technique is used to avoid the pressure odd–even decoupling of the collocated grid arrangement. To increase the stability of the method, a two‐step predictor–corrector time integration scheme is employed. At each step, the projection method is used to calculate the hydrodynamic pressure and to satisfy the continuity equation. The robustness and accuracy of the method is illustrated with a series of numerical experiments including thermally driven cavity, polar cavity, three‐dimensional cavity, and direct numerical simulation of non‐isothermal turbulent channel flow. Copyright © 2012 John Wiley & Sons, Ltd.     

A high‐order mass‐lumping procedure for B‐spline collocation method with application to incompressible flow simulations

This paper presents new developments of the staggered spline collocation method for cost‐effective solution to the incompressible Navier–Stokes equations. Maximal decoupling of the velocity and the pressure is obtained by using the fractional step method of Gresho and Chan, allowing the solution to sparse elliptic problems only. In order to preserve the high‐accuracy of the B‐spline method, this fractional step scheme is used in association with a sparse approximation to the inverse of the consistent mass matrix. Such an approximation is constructed from local spline interpolation method, and represents a high‐order generalization of the mass‐lumping technique of the finite‐element method. A numerical investigation of the accuracy and the computational efficiency of the resulting semi‐consistent spline collocation schemes is presented. These schemes generate a stable and accurate unsteady Navier–Stokes solver, as assessed by benchmark computations. Copyright © 2003 John Wiley & Sons, Ltd.     

Pseudostress–velocity formulation for incompressible Navier–Stokes equations

This paper presents a numerical algorithm using the pseudostress–velocity formulation to solve incompressible Newtonian flows. The pseudostress–velocity formulation is a variation of the stress–velocity formulation, which does not require symmetric tensor spaces in the finite element discretization. Hence its discretization is greatly simplified. The discrete system is further decoupled into an H ( div ) problem for the pseudostress and a post‐process resolving the velocity. This can be done conveniently by using the penalty method for steady‐state flows or by using the time discretization for nonsteady‐state flows. We apply this formulation to the 2D lid‐driven cavity problem and study its grid convergence rate. Also, computational results of the time‐dependent‐driven cavity problem and the flow past rectangular problem are reported. Copyright © 2009 John Wiley & Sons, Ltd.     

A consistent splitting scheme for unsteady incompressible viscous flows I. Dirichlet boundary condition and applications

A well‐recognized approach for handling the incompressibility constraint by operating directly on the discretized Navier–Stokes equations is used to obtain the decoupling of the pressure from the velocity field. By following the current developments by Guermond and Shen, the possibilities of obtaining accurate pressure and reducing boundary‐layer effect for the pressure are analysed. The present study mainly reports the numerical solutions of an unsteady Navier–Stokes problem based on the so‐called consistent splitting scheme (J. Comput. Phys. 2003; 192 :262–276). At the same time the Dirichlet boundary value conditions are considered. The accuracy of the method is carefully examined against the exact solution for an unsteady flow physics problem in a simply connected domain. The effectiveness is illustrated viz. several computations of 2D double lid‐driven cavity problems. Copyright © 2005 John Wiley & Sons, Ltd.     

A stable unstructured finite volume method for parallel large-scale viscoelastic fluid flow calculations

A new stable unstructured finite volume method is presented for parallel large-scale simulation of viscoelastic fluid flows. The numerical method is based on the side-centered finite volume method where the velocity vector components are defined at the mid-point of each cell face, while the pressure term and the extra stress tensor are defined at element centroids. The present arrangement of the primitive variables leads to a stable numerical scheme and it does not require any ad-hoc modifications in order to enhance the pressure–velocity–stress coupling. The log-conformation representation proposed in [R. Fattal, R. Kupferman, Constitutive laws for the matrix–logarithm of the conformation tensor, J. Non-Newtonian Fluid Mech. 123 (2004) 281–285] has been implemented in order improve the limiting Weissenberg numbers in the proposed finite volume method. The time stepping algorithm used decouples the calculation of the polymeric stress by solution of a hyperbolic constitutive equation from the evolution of the velocity and pressure fields by solution of a generalized Stokes problem. The resulting algebraic linear systems are solved using the FGMRES(m) Krylov iterative method with the restricted additive Schwarz preconditioner for the extra stress tensor and the geometric non-nested multilevel preconditioner for the Stokes system. The implementation of the preconditioned iterative solvers is based on the PETSc library for improving the efficiency of the parallel code. The present numerical algorithm is validated for the Kovasznay flow, the flow of an Oldroyd-B fluid past a confined circular cylinder in a channel and the three-dimensional flow of an Oldroyd-B fluid around a rigid sphere falling in a cylindrical tube. Parallel large-scale calculations are presented up to 523,094 quadrilateral elements in two-dimension and 1,190,376 hexahedral elements in three-dimension.     

A high-order parallel Eulerian-Lagrangian algorithm for advection-diffusion problems on unstructured meshes

In this paper, we present a high-order discontinuous Galerkin Eulerian-Lagrangian method for the solution of advection-diffusion problems on staggered unstructured meshes in two and three space dimensions. The particle trajectories are tracked backward in time by means of a high-order representation of the velocity field and a linear mapping from the physical to a reference system, hence obtaining a very simple and efficient strategy that permits to follow the Lagrangian trajectories throughout the computational domain. The use of an Eulerian-Lagrangian discretization increases the overall computational efficiency of the scheme because it is the only explicit method for the discretization of convective terms that admits large time steps without imposing a Courant-Friedrichs-Lewy–type stability condition. This property is fully exploited in this work by relying on a semi-implicit discretization of the incompressible Navier-Stokes equations, in which the pressure is discretized implicitly; thus, the sound speed does not play any role in the restriction of the maximum admissible time step. The resulting mild Courant-Friedrichs-Lewy stability condition, which is based only on the fluid velocity, is here overcome by the adoption of the Eulerian-Lagrangian method for the advection terms and an implicit scheme for the diffusive part of the governing equations. As a consequence, the novel algorithm is able to run simulation with a time step that is defined by the user, depending on the desired efficiency and time scale of the physical phenomena under consideration. Finally, a complete Message Passing Interface parallelization of the code is presented, showing that our approach can reach up to 96% of scaling efficiency.     

Numerical solution of three‐dimensional velocity–vorticity Navier–Stokes equations by finite difference method

This paper describes the finite difference numerical procedure for solving velocity–vorticity form of the Navier–Stokes equations in three dimensions. The velocity Poisson equations are made parabolic using the false‐transient technique and are solved along with the vorticity transport equations. The parabolic velocity Poisson equations are advanced in time using the alternating direction implicit (ADI) procedure and are solved along with the continuity equation for velocities, thus ensuring a divergence‐free velocity field. The vorticity transport equations in conservative form are solved using the second‐order accurate Adams–Bashforth central difference scheme in order to assure divergence‐free vorticity field in three dimensions. The velocity and vorticity Cartesian components are discretized using a central difference scheme on a staggered grid for accuracy reasons. The application of the ADI procedure for the parabolic velocity Poisson equations along with the continuity equation results in diagonally dominant tri‐diagonal matrix equations. Thus the explicit method for the vorticity equations and the tri‐diagonal matrix algorithm for the Poisson equations combine to give a simplified numerical scheme for solving three‐dimensional problems, which otherwise requires enormous computational effort. For three‐dimensional‐driven cavity flow predictions, the present method is found to be efficient and accurate for the Reynolds number range 100?Re?2000. Copyright © 2004 John Wiley & Sons, Ltd.     

The (9,5) HOC formulation for the transient Navier–Stokes equations in primitive variable

We recently proposed an improved (9,5) higher order compact (HOC) scheme for the unsteady two‐dimensional (2‐D) convection–diffusion equations. Because of using only five points at the current time level in the discretization procedure, the scheme was seen to be computationally more efficient than its predecessors. It was also seen to capture very accurately the solution of the unsteady 2‐D Navier–Stokes (N–S) equations for incompressible viscous flows in the stream function–vorticity (ψ – ω) formulation. In this paper, we extend the scope of the scheme for solving the unsteady incompressible N–S equations based on primitive variable formulation on a collocated grid. The parabolic momentum equations are solved for the velocity field by a time‐marching strategy and the pressure is obtained by discretizing the elliptic pressure Poisson equation by the steady‐state form of the (9,5) scheme with the Neumann boundary conditions. In particular, for pressure, we adopt a strategy on the collocated grid in conjunction with ideas borrowed from the staggered grid approach in finite volume. We first apply this extension to a problem having analytical solution and then to the famous lid‐driven square cavity problem. We also apply our formulation to the backward‐facing step problem to see how the method performs for external flow problems. The results are presented and are compared with established numerical results. This new approach is seen to produce excellent comparison in all the cases. Copyright © 2007 John Wiley & Sons, Ltd.     

AN APPROXIMATE PROJECTION SCHEME FOR INCOMPRESSIBLE FLOW USING SPECTRAL ELEMENTS

An approximate projection scheme based on the pressure correction method is proposed to solve the Navier–Stokes equations for incompressible flow. The algorithm is applied to the continuous equations; however, there are no problems concerning the choice of boundary conditions of the pressure step. The resulting velocity and pressure are consistent with the original system. For the spatial discretization a high-order spectral element method is chosen. The high-order accuracy allows the use of a diagonal mass matrix, resulting in a very efficient algorithm. The properties of the scheme are extensively tested by means of an analytical test example. The scheme is further validated by simulating the laminar flow over a backward-facing step.     

Finite difference scheme for the solution of fluid flow problems on non‐staggered grids

A new finite difference methodology is developed for the solution of computational fluid dynamics problems that do not require the use of staggered grid systems. Previous successful and robust non‐staggered methods, which used primitive variables and mass conservation in order to solve the pressure field, either interpolate cell‐face velocities or interpolate the pressure gradients in a special way, usually with an upwind‐bias to avoid the problem of odd–even coupling between the velocity and pressure fields. The new methodology presented does not detail a ‘special interpolation procedure for a primitive variable’, however, it manages to avoid the problem of odd–even coupling. The odd–even coupling is avoided by applying fourth‐order dissipation to the pressure field. It is shown that this approach can be regarded as a modified Rhie and Chow scheme. The method is implemented using a SIMPLE‐type algorithm and is applied to two test problems: laminar flow over a backward‐facing step and laminar flow in a square cavity with a driven lid. Good agreement is obtained between the numerical solutions and the corresponding benchmark solutions. The pressure dissipation term was found to successfully suppress wiggles in the pressure field. Copyright © 2000 John Wiley & Sons, Ltd.     

A stable second‐order mass‐weighted upwind scheme for unstructured meshes

In this paper, an original second‐order upwind scheme for convection terms is described and implemented in the context of a Control‐Volume Finite‐Element Method (CVFEM). The proposed scheme is a second‐order extension of the first‐order MAss‐Weighted upwind (MAW) scheme proposed by Saabas and Baliga (Numer. Heat Transfer 1994; 26B :381–407). The proposed second‐order scheme inherits the well‐known stability characteristics of the MAW scheme, but exhibits less artificial viscosity and ensures much higher accuracy. Consequently, and in contrast with nearly all second‐order upwind schemes available in the literature, the proposed second‐order MAW scheme does not need limiters. Some test cases including two pure convection problems, the driven cavity and steady and unsteady flows over a circular cylinder, have been undertaken successfully to validate the new scheme. The verification tests show that the proposed scheme exhibits a low level of artificial viscosity in the pure convection problems; exhibits second‐order accuracy for the driven cavity; gives accurate reattachment lengths for low‐Reynolds steady flow over a circular cylinder; and gives constant‐amplitude vortex shedding for the case of high‐Reynolds unsteady flow over a circular cylinder. Copyright © 2005 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.