Numerical solution of the incompressible Navier–Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier–Stokes equations is presented. The scheme is constructed on a staggered‐mesh grid system. The convection terms are discretized with a fifth‐order‐accurate upwind compact difference approximation, the viscous terms are discretized with a sixth‐order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth‐order difference approximation on a cell‐centered mesh. Time advancement uses a three‐stage Runge–Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented. Copyright © 1999 John Wiley & Sons, Ltd.     

Implicit weighted essentially non‐oscillatory schemes for the incompressible Navier–Stokes equations

A class of lower–upper/approximate factorization (LUAF) implicit weighted essentially non‐oscillatory (ENO; WENO) schemes for solving the two‐dimensional incompressible Navier–Stokes equations in a generalized co‐ordinate system is presented. The algorithm is based on the artificial compressibility formulation, and symmetric Gauss–Seidel relaxation is used for computing steady state solutions while symmetric successive overrelaxation is used for treating time‐dependent flows. WENO spatial operators are employed for inviscid fluxes and central differencing for viscous fluxes. Internal and external viscous flow test problems are presented to verify the numerical schemes. The use of a WENO spatial operator not only enhances the accuracy of solutions but also improves the convergence rate for the steady state computation as compared with using the ENO counterpart. It is found that the present solutions compare well with exact solutions, experimental data and other numerical results. Copyright © 1999 John Wiley & Sons, Ltd.     

Combined compact difference method for solving the incompressible Navier–Stokes equations

This paper presents a numerical method for solving the two‐dimensional unsteady incompressible Navier–Stokes equations in a vorticity–velocity formulation. The method is applicable for simulating the nonlinear wave interaction in a two‐dimensional boundary layer flow. It is based on combined compact difference schemes of up to 12th order for discretization of the spatial derivatives on equidistant grids and a fourth‐order five‐ to six‐alternating‐stage Runge–Kutta method for temporal integration. The spatial and temporal schemes are optimized together for the first derivative in a downstream direction to achieve a better spectral resolution. In this method, the dispersion and dissipation errors have been minimized to simulate physical waves accurately. At the same time, the schemes can efficiently suppress numerical grid‐mesh oscillations. The results of test calculations on coarse grids are in good agreement with the linear stability theory and comparable with other works. The accuracy and the efficiency of the current code indicate its potential to be extended to three‐dimensional cases in which full boundary layer transition happens. Copyright © 2011 John Wiley & Sons, Ltd.     

Fourth‐order accurate compact‐difference discretization method for Helmholtz and incompressible Navier–Stokes equations

A fourth‐order accurate solution method for the three‐dimensional Helmholtz equations is described that is based on a compact finite‐difference stencil for the Laplace operator. Similar discretization methods for the Poisson equation have been presented by various researchers for Dirichlet boundary conditions. Here, the complicated issue of imposing Neumann boundary conditions is described in detail. The method is then applied to model Helmholtz problems to verify the accuracy of the discretization method. The implementation of the solution method is also described. The Helmholtz solver is used as the basis for a fourth‐order accurate solver for the incompressible Navier–Stokes equations. Numerical results obtained with this Navier–Stokes solver for the temporal evolution of a three‐dimensional instability in a counter‐rotating vortex pair are discussed. The time‐accurate Navier–Stokes simulations show the resolving properties of the developed discretization method and the correct prediction of the initial growth rate of the three‐dimensional instability in the vortex pair. Copyright © 2004 John Wiley & Sons, Ltd.     

A compact finite difference method on staggered grid for Navier–Stokes flows

Compact finite difference methods feature high‐order accuracy with smaller stencils and easier application of boundary conditions, and have been employed as an alternative to spectral methods in direct numerical simulation and large eddy simulation of turbulence. The underpinning idea of the method is to cancel lower‐order errors by treating spatial Taylor expansions implicitly. Recently, some attention has been paid to conservative compact finite volume methods on staggered grid, but there is a concern about the order of accuracy after replacing cell surface integrals by average values calculated at centres of cell surfaces. Here we introduce a high‐order compact finite difference method on staggered grid, without taking integration by parts. The method is implemented and assessed for an incompressible shear‐driven cavity flow at Re = 10³, a temporally periodic flow at Re = 10⁴, and a spatially periodic flow at Re = 10⁴. The results demonstrate the success of the method. Copyright © 2006 John Wiley & Sons, Ltd.     

Comparison of wide and compact fourth‐order formulations of the Navier–Stokes equations

In this study the numerical performances of wide and compact fourth‐order formulation of the steady 2D incompressible Navier–Stokes equations will be investigated and compared with each other. The benchmark driven cavity flow problem will be solved using both wide and compact fourth‐order formulations and the numerical performances of both formulations will be presented and also the advantages and disadvantages of both formulations will be discussed. Copyright © 2008 John Wiley & Sons, Ltd.     

An implicit edge‐based ALE method for the incompressible Navier–Stokes equations

A new finite volume method for the incompressible Navier–Stokes equations, expressed in arbitrary Lagrangian–Eulerian (ALE) form, is presented. The method uses a staggered storage arrangement for the pressure and velocity variables and adopts an edge‐based data structure and assembly procedure which is valid for arbitrary n‐sided polygonal meshes. Edge formulas are presented for assembling the ALE form of the momentum and pressure equations. An implicit multi‐stage time integrator is constructed that is geometrically conservative to the precision of the arithmetic used in the computation. The method is shown to be second‐order‐accurate in time and space for general time‐dependent polygonal meshes. The method is first evaluated using several well‐known unsteady incompressible Navier–Stokes problems before being applied to a periodically forced aeroelastic problem and a transient free surface problem. Published in 2003 by John Wiley & Sons, Ltd.     

Stretched Cartesian grids for solution of the incompressible Navier–Stokes equations

Two Cartesian grid stretching functions are investigated for solving the unsteady incompressible Navier–Stokes equations using the pressure–velocity formulation. The first function is developed for the Fourier method and is a generalization of earlier work. This function concentrates more points at the centre of the computational box while allowing the box to remain finite. The second stretching function is for the second‐order central finite difference scheme, which uses a staggered grid in the computational domain. This function is derived to allow a direct discretization of the Laplacian operator in the pressure equation while preserving the consistent behaviour exhibited by the uniform grid scheme. Both functions are analysed for their effects on the matrix of the discretized pressure equation. It is shown that while the second function does not spoil the matrix diagonal dominance, the first one can. Limits to stretching of the first method are derived for the cases of mappings in one and two directions. A limit is also derived for the second function in order to prevent a strong distortion of a sine wave. The performances of the two types of stretching are examined in simulations of periodic co‐flowing jets and a time developing boundary layer. Copyright © 2000 John Wiley & Sons, Ltd.     

A Chimera method for the incompressible Navier–Stokes equations

The Chimera method was developed three decades ago as a meshing simplification tool. Different components are meshed independently and then glued together using a domain decomposition technique to couple the equations solved on each component. This coupling is achieved via transmission conditions (in the finite element context) or by imposing the continuity of fluxes (in the finite volume context). Historically, the method has then been used extensively to treat moving objects, as the independent meshes are free to move with respect to the others. At each time step, the main task consists in recomputing the interpolation of the transmission conditions or fluxes. This paper presents a Chimera method applied to the Navier–Stokes equations. After an introduction on the Chimera method, we describe in two different sections the two independent steps of the method: the hole cutting to create the interfaces of the subdomains and the coupling of the subdomains. Then, we present the Navier–Stokes solver considered in this work. Implementation aspects are then detailed in order to apply efficiently the method to this specific parallel Navier–Stokes solver. We conclude with some examples to demonstrate the reliability and application of the proposed method. Copyright © 2014 John Wiley & Sons, Ltd.     

Colocated schemes for the incompressible Navier–Stokes equations on non‐smooth grids for two‐dimensional problems

The accuracy of colocated finite volume schemes for the incompressible Navier–Stokes equations on non‐smooth curvilinear grids is investigated. A frequently used scheme is found to be quite inaccurate on non‐smooth grids. In an attempt to improve the accuracy on such grids, three other schemes are described and tested. Two of these are found to give satisfactory results. Copyright © 2000 John Wiley & Sons, Ltd.     

Schwarz domain decomposition for the incompressible Navier–Stokes equations in general co‐ordinates

This paper describes a domain decomposition method for the incompressible Navier–Stokes equations in general co‐ordinates. Domain decomposition techniques are needed for solving flow problems in complicated geometries while retaining structured grids on each of the subdomains. This is the so‐called block‐structured approach. It enables the use of fast vectorized iterative methods on the subdomains. The Navier–Stokes equations are discretized on a staggered grid using finite volumes. The pressure‐correction technique is used to solve the momentum equations together with incompressibility conditions. Schwarz domain decomposition is used to solve the momentum and pressure equations on the composite domain. Convergence of domain decomposition is accelerated by a GMRES Krylov subspace method. Computations are presented for a variety of flows. Copyright © 2000 John Wiley & Sons, Ltd.     

Velocity–pressure coupling in finite difference formulations for the Navier–Stokes equations

A new numerical procedure for solving the two‐dimensional, steady, incompressible, viscous flow equations on a staggered Cartesian grid is presented in this paper. The proposed methodology is finite difference based, but essentially takes advantage of the best features of two well‐established numerical formulations, the finite difference and finite volume methods. Some weaknesses of the finite difference approach are removed by exploiting the strengths of the finite volume method. In particular, the issue of velocity–pressure coupling is dealt with in the proposed finite difference formulation by developing a pressure correction equation using the SIMPLE approach commonly used in finite volume formulations. However, since this is purely a finite difference formulation, numerical approximation of fluxes is not required. Results presented in this paper are based on first‐ and second‐order upwind schemes for the convective terms. This new formulation is validated against experimental and other numerical data for well‐known benchmark problems, namely developing laminar flow in a straight duct, flow over a backward‐facing step, and lid‐driven cavity flow. Copyright © 2010 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 multigrid method for steady incompressible Navier–Stokes equations based on partial flux splitting

Flux splitting is applied to the convective part of the steady Navier–Stokes equations for incompressible flow. Partial upwind differences are introduced in the split first-order part, while central differences are used in the second-order part. The discrete set of equations obtained is positive, so that it can be solved by collective variants of relaxation methods. The partial upwinding is optimized in the same way as for a scalar convection–diffusion equation, but involving several Peclet numbers. It is shown that with the optimum partial upwinding accurate results can be obtained. A full multigrid method in W-cycle form, using red–black successive under-relaxation, injection and bilinear interpolation, is described. The efficiency of this method is demonstrated.     

Solution of the incompressible Navier–Stokes equations by the method of lines

The numerical method of lines (NUMOL) is a numerical technique used to solve efficiently partial differential equations. In this paper, the NUMOL is applied to the solution of the two‐dimensional unsteady Navier–Stokes equations for incompressible laminar flows in Cartesian coordinates. The Navier–Stokes equations are first discretized (in space) on a staggered grid as in the Marker and Cell scheme. The discretized Navier–Stokes equations form an index 2 system of differential algebraic equations, which are afterwards reduced to a system of ordinary differential equations (ODEs), using the discretized form of the continuity equation. The pressure field is computed solving a discrete pressure Poisson equation. Finally, the resulting ODEs are solved using the backward differentiation formulas. The proposed method is illustrated with Dirichlet boundary conditions through applications to the driven cavity flow and to the backward facing step flow. Copyright © 2015 John Wiley & Sons, Ltd.     

Multigrid solution of the incompressible Navier–Stokes equations for three‐dimensional recirculating flow: coupled and decoupled smoothers compared

A comparison of multigrid methods for solving the incompressible Navier–Stokes equations in three dimensions is presented. The continuous equations are discretised on staggered grids using a second‐order monotonic scheme for the convective terms and implemented in defect correction form. The convergence characteristics of a decoupled method (SIMPLE) are compared with those of the cellwise coupled method (SCGS). The convergence rates obtained for computations of the three‐dimensional lid‐driven cavity problem are found to be very similar to those obtained for computations of the corresponding two‐dimensional problem with comparable grid density. Although the convergence rate of SCGS is thus superior to that of SIMPLE, the decoupled method is found to be more efficient computationally and requires less computing time for a given level of convergence. The linewise implementation of the coupled method (CLGS) is also investigated and shown to be more efficient than SCGS, although the convergence rate and computing time required per cycle are both found to depend on the direction of sweep. The optimal implementation of CLGS is found to be only marginally more effective than SIMPLE, but a change to the structure of the data storage would increase the advantage. Copyright © 1999 John Wiley & Sons, Ltd.     

A high‐order scheme for the incompressible Navier–Stokes equations with open boundary condition

We present a numerical scheme to solve the incompressible Navier–Stokes equations with open boundary condition. After replacing the incompressibility constraint by the pressure Poisson equation, the key is how to give an appropriate boundary condition for the pressure Poisson equation. We propose a new boundary condition for the pressure on the open boundary. Some numerical experiments are presented to verify the accuracy and stability of scheme. Copyright © 2013 John Wiley & Sons, Ltd.     

Newton linearization of the incompressible Navier–Stokes equations

The present study aims to accelerate the convergence to incompressible Navier–Stokes solution. For the sake of computational efficiency, Newton linearization of equations is invoked on non‐staggered grids to shorten the sequence to the final solution of the non‐linear differential system of equations. For the sake of accuracy, the resulting convection–diffusion–reaction finite‐difference equation is solved line‐by‐line using the proposed nodally exact one‐dimensional scheme. The matrix size is reduced and, at the same time, the CPU time is considerably saved due to the decrease of stencil points. The effectiveness of the implemented Newton linearization is demonstrated through computational exercises. Copyright © 2004 John Wiley & Sons, Ltd.     

Steady‐state solution of incompressible Navier–Stokes equations using discrete least‐squares meshless method

A fractional step method for the solution of the steady state incompressible Navier–Stokes equations is proposed in this paper in conjunction with a meshless method, named discrete least‐squares meshless (DLSM). The proposed fractional step method is a first‐order accurate scheme, named semi‐incremental fractional step method, which is a general form of the previous first‐order fractional step methods, i.e. non‐incremental and incremental schemes. One of the most important advantages of the proposed scheme is its capability to use large time step sizes for the solution of incompressible Navier–Stokes equations. DLSM method uses moving least‐squares shape functions for function approximation and discrete least‐squares technique for discretization of the governing differential equations and their boundary conditions. As there is no need for a background mesh, the DLSM method can be called a truly meshless method and enjoys symmetric and positive‐definite properties. Several numerical examples are used to demonstrate the ability and the efficiency of the proposed scheme and the discrete least‐squares meshless method. The results are shown to compare favorably with those of the previously published works. Copyright © 2010 John Wiley & Sons, Ltd.