Flux‐difference splitting‐based upwind compact schemes for the incompressible Navier–Stokes equations

Third‐order and fifth‐order upwind compact finite difference schemes based on flux‐difference splitting are proposed for solving the incompressible Navier–Stokes equations in conjunction with the artificial compressibility (AC) method. Since the governing equations in the AC method are hyperbolic, flux‐difference splitting (FDS) originally developed for the compressible Euler equations can be used. In the present upwind compact schemes, the split derivatives for the convective terms at grid points are linked to the differences of split fluxes between neighboring grid points, and these differences are computed by using FDS. The viscous terms are approximated with a sixth‐order central compact scheme. Comparisons with 2D benchmark solutions demonstrate that the present compact schemes are simple, efficient, and high‐order accurate. Copyright © 2008 John Wiley & Sons, Ltd.     

Analysis of super compact finite difference method and application to simulation of vortex–shock interaction

Turbulence and aeroacoustic noise high‐order accurate schemes are required, and preferred, for solving complex flow fields with multi‐scale structures. In this paper a super compact finite difference method (SCFDM) is presented, the accuracy is analysed and the method is compared with a sixth‐order traditional and compact finite difference approximation. The comparison shows that the sixth‐order accurate super compact method has higher resolving efficiency. The sixth‐order super compact method, with a three‐stage Runge–Kutta method for approximation of the compressible Navier–Stokes equations, is used to solve the complex flow structures induced by vortex–shock interactions. The basic nature of the near‐field sound generated by interaction is studied. Copyright © 2001 John Wiley & Sons, Ltd.     

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.     

Implementation of high‐order compact finite‐difference method to parabolized Navier–Stokes schemes

The numerical solution to the parabolized Navier–Stokes (PNS) and globally iterated PNS (IPNS) equations for accurate computation of hypersonic axisymmetric flowfields is obtained by using the fourth‐order compact finite‐difference method. The PNS and IPNS equations in the general curvilinear coordinates are solved by using the implicit finite‐difference algorithm of Beam and Warming type with a high‐order compact accuracy. A shock‐fitting procedure is utilized in both compact PNS and IPNS schemes to obtain accurate solutions in the vicinity of the shock. The main advantage of the present formulation is that the basic flow variables and their first and second derivatives are simultaneously computed with the fourth‐order accuracy. The computations are carried out for a benchmark case: hypersonic axisymmetric flow over a blunt cone at Mach 8. A sensitivity study is performed for the basic flowfield, including profiles and their derivatives obtained from the fourth‐order compact PNS and IPNS solutions, and the effects of grid size and numerical dissipation term used are discussed. The present results for the flowfield variables and also their derivatives are compared with those of other basic flow models to demonstrate the accuracy and efficiency of the proposed method. The present work represents the first known application of a high‐order compact finite‐difference method to the PNS schemes, which are computationally more efficient than Navier–Stokes solutions. Copyright © 2008 John Wiley & Sons, Ltd.     

Implicit preconditioned high‐order compact scheme for the simulation of the three‐dimensional incompressible Navier–Stokes equations with pseudo‐compressibility method

This paper combines the pseudo‐compressibility procedure, the preconditioning technique for accelerating the time marching for stiff hyperbolic equations, and high‐order accurate central compact scheme to establish the code for efficiently and accurately solving incompressible flows numerically based on the finite difference discretization. The spatial scheme consists of the sixth‐order compact scheme and 10th‐order numerical filter operator for guaranteeing computational stability. The preconditioned pseudo‐compressible Navier–Stokes equations are marched temporally using the implicit lower–upper symmetric Gauss–Seidel time integration method, and the time accuracy is improved by the dual‐time step method for the unsteady problems. The efficiency and reliability of the present procedure are demonstrated by applications to Taylor decaying vortices phenomena, double periodic shear layer rolling‐up problem, laminar flow over a flat plate, low Reynolds number unsteady flow around a circular cylinder at Re = 200, high Reynolds number turbulence flow past the S809 airfoil, and the three‐dimensional flows through two 90°curved ducts of square and circular cross sections, respectively. It is found that the numerical results of the present algorithm are in good agreement with theoretical solutions or experimental data. Copyright © 2011 John Wiley & Sons, Ltd.     

An upwind compact mpact approach with group velocity control for compressible flow fields

In this paper we construct an upwind compact finite difference scheme with group velocity control for better simulation of compressible flow fields. Compared with traditional difference schemes, compact schemes have higher accuracy for the same stencil width. By means of the characteristic analysis of the operators, the group velocity of wave packets will be controlled to suppress the non‐physical oscillations in numerical solutions. In numerical simulation of the 3D compressible flow fields the third‐order accurate upwind compact operator is used to approximate the derivatives in the convection terms of the compressible N–S equations, the traditional finite difference scheme is used to approximate the viscous terms. Numerical solutions indicate that the method is satisfactory. Copyright © 2004 John Wiley & Sons, Ltd.     

Implementing a high‐order accurate implicit operator scheme for solving steady incompressible viscous flows using artificial compressibility method

This paper uses a fourth‐order compact finite‐difference scheme for solving steady incompressible flows. The high‐order compact method applied is an alternating direction implicit operator scheme, which has been used by Ekaterinaris for computing two‐dimensional compressible flows. Herein, this numerical scheme is efficiently implemented to solve the incompressible Navier–Stokes equations in the primitive variables formulation using the artificial compressibility method. For space discretizing the convective fluxes, fourth‐order centered spatial accuracy of the implicit operators is efficiently obtained by performing compact space differentiation in which the method uses block‐tridiagonal matrix inversions. To stabilize the numerical solution, numerical dissipation terms and/or filters are used. In this study, the high‐order compact implicit operator scheme is also extended for computing three‐dimensional incompressible flows. The accuracy and efficiency of this high‐order compact method are demonstrated for different incompressible flow problems. A sensitivity study is also conducted to evaluate the effects of grid resolution and pseudocompressibility parameter on accuracy and convergence rate of the solution. The effects of filtering and numerical dissipation on the solution are also investigated. Test cases considered herein for validating the results are incompressible flows in a 2‐D backward facing step, a 2‐D cavity and a 3‐D cavity at different flow conditions. Results obtained for these cases are in good agreement with the available numerical and experimental results. The study shows that the scheme is robust, efficient and accurate for solving incompressible flow problems. Copyright © 2010 John Wiley & Sons, Ltd.     

Application of a preconditioned high‐order accurate artificial compressibility‐based incompressible flow solver in wide range of Reynolds numbers

In the present study, the preconditioned incompressible Navier‐Stokes equations with the artificial compressibility method formulated in the generalized curvilinear coordinates are numerically solved by using a high‐order compact finite‐difference scheme for accurately and efficiently computing the incompressible flows in a wide range of Reynolds numbers. A fourth‐order compact finite‐difference scheme is utilized to accurately discretize the spatial derivative terms of the governing equations, and the time integration is carried out based on the dual time‐stepping method. The capability of the proposed solution methodology for the computations of the steady and unsteady incompressible viscous flows from very low to high Reynolds numbers is investigated through the simulation of different 2‐dimensional benchmark problems, and the results obtained are compared with the existing analytical, numerical, and experimental data. A sensitivity analysis is also performed to evaluate the effects of the size of the computational domain and other numerical parameters on the accuracy and performance of the solution algorithm. The present solution procedure is also extended to 3 dimensions and applied for computing the incompressible flow over a sphere. Indications are that the application of the preconditioning in the solution algorithm together with the high‐order discretization method in the generalized curvilinear coordinates provides an accurate and robust solution method for simulating the incompressible flows over practical geometries in a wide range of Reynolds numbers including the creeping flows.     

High‐order numerical simulations of flow‐induced noise

In this paper, the flow/acoustics splitting method for predicting flow‐generated noise is further developed by introducing high‐order finite difference schemes. The splitting method consists of dividing the acoustic problem into a viscous incompressible flow part and an inviscid acoustic part. The incompressible flow equations are solved by a second‐order finite volume code EllipSys2D/3D. The acoustic field is obtained by solving a set of acoustic perturbation equations forced by flow quantities. The incompressible pressure and velocity form the input to the acoustic equations. The present work is an extension of our acoustics solver, with the introduction of high‐order schemes for spatial discretization and a Runge–Kutta scheme for time integration. To achieve low dissipation and dispersion errors, either Dispersion‐Relation‐Preserving (DRP) schemes or optimized compact finite difference schemes are used for the spatial discretizations. Applications and validations of the new acoustics solver are presented for benchmark aeroacoustic problems and for flow over an NACA 0012 airfoil. Copyright © 2010 John Wiley & Sons, Ltd.     

High‐order compact finite difference schemes for the vorticity–divergence representation of the spherical shallow water equations

This work is devoted to the application of the super compact finite difference method (SCFDM) and the combined compact finite difference method (CCFDM) for spatial differencing of the spherical shallow water equations in terms of vorticity, divergence, and height. The fourth‐order compact, the sixth‐order and eighth‐order SCFDM, and the sixth‐order and eighth‐order CCFDM schemes are used for the spatial differencing. To advance the solution in time, a semi‐implicit Runge–Kutta method is used. In addition, to control the nonlinear instability, an eighth‐order compact spatial filter is employed. For the numerical solution of the elliptic equations in the problem, a direct hybrid method, which consists of a high‐order compact scheme for spatial differencing in the latitude coordinate and a fast Fourier transform in longitude coordinate, is utilized. The accuracy and convergence rate for all methods are verified against exact analytical solutions. Qualitative and quantitative assessments of the results for an unstable barotropic mid‐latitude zonal jet employed as an initial condition are addressed. It is revealed that the sixth‐order and eighth‐order CCFDMs and SCFDMs lead to a remarkable improvement of the solution over the fourth‐order compact method. Copyright © 2015 John Wiley & Sons, Ltd.     

A high‐order compact finite‐difference lattice Boltzmann method for simulation of steady and unsteady incompressible flows

A high‐order compact finite‐difference lattice Boltzmann method (CFDLBM) is proposed and applied to accurately compute steady and unsteady incompressible flows. Herein, the spatial derivatives in the lattice Boltzmann equation are discretized by using the fourth‐order compact FD scheme, and the temporal term is discretized with the fourth‐order Runge–Kutta scheme to provide an accurate and efficient incompressible flow solver. A high‐order spectral‐type low‐pass compact filter is used to stabilize the numerical solution. An iterative initialization procedure is presented and applied to generate consistent initial conditions for the simulation of unsteady flows. A sensitivity study is also conducted to evaluate the effects of grid size, filtering, and procedure of boundary conditions implementation on accuracy and convergence rate of the solution. The accuracy and efficiency of the proposed solution procedure based on the CFDLBM method are also examined by comparison with the classical LBM for different flow conditions. Two test cases considered herein for validating the results of the incompressible steady flows are a two‐dimensional (2‐D) backward‐facing step and a 2‐D cavity at different Reynolds numbers. Results of these steady solutions computed by the CFDLBM are thoroughly compared with those of a compact FD Navier–Stokes flow solver. Three other test cases, namely, a 2‐D Couette flow, the Taylor's vortex problem, and the doubly periodic shear layers, are simulated to investigate the accuracy of the proposed scheme in solving unsteady incompressible flows. Results obtained for these test cases are in good agreement with the analytical solutions and also with the available numerical and experimental results. The study shows that the present solution methodology is robust, efficient, and accurate for solving steady and unsteady incompressible flow problems even at high Reynolds numbers. Copyright © 2014 John Wiley & Sons, Ltd.     

Verification of a mixed high‐order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high‐order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier–Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity–velocity formulation for a two‐dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high‐order compact and non‐compact finite‐differences from fourth‐order to sixth‐order of accuracy. The other scheme used spectral methods instead of finite‐difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed‐order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright © 2009 John Wiley & Sons, Ltd.     

A fourth‐order compact finite difference scheme for the steady stream function–vorticity formulation of the Navier–Stokes/Boussinesq equations

A fourth‐order compact finite difference scheme on the nine‐point 2D stencil is formulated for solving the steady‐state Navier–Stokes/Boussinesq equations for two‐dimensional, incompressible fluid flow and heat transfer using the stream function–vorticity formulation. The main feature of the new fourth‐order compact scheme is that it allows point‐successive overrelaxation (SOR) or point‐successive underrelaxation iteration for all Rayleigh numbers Ra of physical interest and all Prandtl numbers Pr attempted. Numerical solutions are obtained for the model problem of natural convection in a square cavity with benchmark solutions and compared with some of the accurate results available in the literature. Copyright © 2003 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.     

A higher‐order Boussinesq‐type model with moving bottom boundary: applications to submarine landslide tsunami waves

A two‐dimensional depth‐integrated numerical model is developed using a fourth‐order Boussinesq approximation for an arbitrary time‐variable bottom boundary and is applied for submarine‐landslide‐generated waves. The mathematical formulation of model is an extension of (4,4) Padé approximant for moving bottom boundary. The mathematical formulations are derived based on a higher‐order perturbation analysis using the expanded form of velocity components. A sixth‐order multi‐step finite difference method is applied for spatial discretization and a sixth‐order Runge–Kutta method is applied for temporal discretization of the higher‐order depth‐integrated governing equations and boundary conditions. The present model is validated using available three‐dimensional experimental data and a good agreement is obtained. Moreover, the present higher‐order model is compared with fully potential three‐dimensional models as well as Boussinesq‐type multi‐layer models in several cases and the differences are discussed. The high accuracy of the present numerical model in considering the nonlinearity effects and frequency dispersion of waves is proven particularly for waves generated in intermediate and deeper water area. Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

A transformation‐free HOC scheme for incompressible viscous flows on nonuniform polar grids

We recently proposed a transformation‐free higher‐order compact (HOC) scheme for two‐dimensional (2‐D) steady convection–diffusion equations on nonuniform Cartesian grids (Int. J. Numer. Meth. Fluids 2004; 44 :33–53). As the scheme was equipped to handle only constant coefficients for the second‐order derivatives, it could not be extended directly to curvilinear coordinates, where they invariably occur as variables. In this paper, we extend the scheme to cylindrical polar coordinates for the 2‐D convection–diffusion equations and more specifically to the 2‐D incompressible viscous flows governed by the Navier–Stokes (N–S) equations. We first apply the formulation to a problem having analytical solution and demonstrate its fourth‐order spatial accuracy. We then apply it to the flow past an impulsively started circular cylinder problem and finally to the driven polar cavity problem. We present our numerical results and compare them with established numerical and analytical and experimental results whenever available. This new approach is seen to produce excellent comparison in all the cases. Copyright © 2009 John Wiley & Sons, Ltd.     

Comments on Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD

The explicit compact difference scheme,proposed in Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD by Lin et al.,published in Applied Mathematics and Mechanics (English Edition),2007,28(7),943-953,has the same performance as the conventional finite difference schemes.It is just another expression of the conventional finite difference schemes. The proposed expression does not have the advantages of a compact difference scheme. Nonetheless,we can more easily obtain and implement compared with the conventional expression in which the coefficients can only be obtained by solving equations,especially for higher accurate schemes.     

Hybrid finite‐volume finite‐difference scheme for the solution of Boussinesq equations

A hybrid scheme composed of finite‐volume and finite‐difference methods is introduced for the solution of the Boussinesq equations. While the finite‐volume method with a Riemann solver is applied to the conservative part of the equations, the higher‐order Boussinesq terms are discretized using the finite‐difference scheme. Fourth‐order accuracy in space for the finite‐volume solution is achieved using the MUSCL‐TVD scheme. Within this, four limiters have been tested, of which van‐Leer limiter is found to be the most suitable. The Adams–Basforth third‐order predictor and Adams–Moulton fourth‐order corrector methods are used to obtain fourth‐order accuracy in time. A recently introduced surface gradient technique is employed for the treatment of the bottom slope. A new model ‘HYWAVE’, based on this hybrid solution, has been applied to a number of wave propagation examples, most of which are taken from previous studies. Examples include sinusoidal waves and bi‐chromatic wave propagation in deep water, sinusoidal wave propagation in shallow water and sinusoidal wave propagation from deep to shallow water demonstrating the linear shoaling properties of the model. Finally, sinusoidal wave propagation over a bar is simulated. The results are in good agreement with the theoretical expectations and published experimental results. Copyright © 2005 John Wiley & Sons, Ltd.     

High‐order accurate numerical solutions of incompressible flows with the artificial compressibility method

A high‐order accurate, finite‐difference method for the numerical solution of incompressible flows is presented. This method is based on the artificial compressibility formulation of the incompressible Navier–Stokes equations. Fourth‐ or sixth‐order accurate discretizations of the metric terms and the convective fluxes are obtained using compact, centred schemes. The viscous terms are also discretized using fourth‐order accurate, centred finite differences. Implicit time marching is performed for both steady‐state and time‐accurate numerical solutions. High‐order, spectral‐type, low‐pass, compact filters are used to regularize the numerical solution and remove spurious modes arising from unresolved scales, non‐linearities, and inaccuracies in the application of boundary conditions. The accuracy and efficiency of the proposed method is demonstrated for test problems. Copyright © 2004 John Wiley & Sons, Ltd.