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.     

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.     

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.     

Simulation of three-dimensional incompressible flows in generalized curvilinear coordinates using a high-order compact finite-difference lattice Boltzmann method

In the present study, a high-order compact finite-difference lattice Boltzmann method is applied for accurately computing 3-D incompressible flows in the generalized curvilinear coordinates to handle practical and realistic geometries with curved boundaries and nonuniform grids. The incompressible form of the 3-D nineteen discrete velocity lattice Boltzmann method is transformed into the generalized curvilinear coordinates. Herein, a fourth-order compact finite-difference scheme and a fourth-order Runge-Kutta scheme are used for the discretization of the spatial derivatives and the temporal term, respectively, in the resulting 3-D nineteen discrete velocity lattice Boltzmann equation to provide an accurate 3-D incompressible flow solver. A high-order spectral-type low-pass compact filtering technique is applied to have a stable solution. All boundary conditions are implemented based on the solution of the governing equations in the 3-D generalized curvilinear coordinates. Numerical solutions of different 3-D benchmark and practical incompressible flow problems are performed to demonstrate the accuracy and performance of the solution methodology presented. Herein, the 2-D cylindrical Couette flow, the decay of a 3-D double shear wave, the cubic lid-driven cavity flow with nonuniform grids, the flow through a square duct with 90° bend and the flow past a sphere at different flow conditions are considered for validating the present computations. Numerical results obtained show the accuracy and robustness of the present solution methodology based on the implementation of the high-order compact finite-difference lattice Boltzman method in the generalized curvilinear coordinates for solving 3-D incompressible flows over practical and realistic geometries.     

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.     

Genuinely multidimensional characteristic‐based scheme for incompressible flows

This paper presents a novel multidimensional characteristic‐based (MCB) upwind method for the solution of incompressible Navier–Stokes equations. As opposed to the conventional characteristic‐based (CB) schemes, it is genuinely multidimensional in that the local characteristic paths, along which information is propagated, are used. For the first time, the multidimensional characteristic structure of incompressible flows modified by artificial compressibility is extracted and used to construct an inherent multidimensional upwind scheme. The new proposed MCB scheme in conjunction with the finite‐volume discretization is employed to model the convective fluxes. Using this formulation, the steady two‐dimensional incompressible flow in a lid‐driven cavity is solved for a wide range of Reynolds numbers. It was found that the new proposed scheme presents more accurate results than the conventional CB scheme in both their first‐ and second‐order counterparts in the case of cavity flow. Also, results obtained with second‐order MCB scheme in some cases are more accurate than the central scheme that in turn provides exact second‐order discretization in this grid. With this inherent upwinding technique for evaluating convective fluxes at cell interfaces, no artificial viscosity is required even at high Reynolds numbers. Another remarkable advantage of MCB scheme lies in its faster convergence rate with respect to the CB scheme that is found to exhibit substantial delays in convergence reported in the literature. The results obtained using new proposed scheme are in good agreement with the standard benchmark solutions in the literature. Copyright © 2007 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.     

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.     

亚、跨、超音速及不可压流动的数值分析方法的研究

为了对亚、跨、超音速及不可压无粘流动进行数值模拟,将LU－SGS方法与预处理方法结合,给出了PLU－SGS方法。方程离散基于有限体积法,采用高阶精度AUSMPW格式。方程求解采用了特征边界条件。通过典型算例的数值试验对比分析,表明PLU－SGS方法可以有效地对亚、跨、超音速及不可压流动进行数值模拟,并具有较高的计算精度和收敛速度。     

Incompressible Roe‐like scheme for turbulent flows

In the current study, numerical investigation of incompressible turbulent flow is presented. By the artificial compressibility method, momentum and continuity equations are coupled. Considering Reynolds averaged Navier–Stokes equations, the Spalart–Allmaras turbulence model, which has accurate results in two‐dimensional problems, is used to calculate Reynolds stresses. For convective fluxes a Roe‐like scheme is proposed for the steady Reynolds averaged Navier–Stokes equations. Also, Jameson averaging method was implemented. In comparison, the proposed characteristics‐based upwind incompressible turbulent Roe‐like scheme, demonstrated very accurate results, high stability, and fast convergence. The fifth‐order Runge–Kutta scheme is used for time discretization. The local time stepping and implicit residual smoothing were applied as the convergence acceleration techniques. Suitable boundary conditions have been implemented considering flow behavior. The problem has been studied at high Reynolds numbers for cross flow around the horizontal circular cylinder and NACA0012 hydrofoil. Results were compared with those of others and a good agreement has been observed. Copyright © 2012 John Wiley & Sons, Ltd.     

Retracted and replaced: A flow‐condition‐based interpolation finite element procedure for triangular grids

A flow‐condition‐based interpolation finite element scheme is presented for use of triangular grids in the solution of the incompressible Navier–Stokes equations. The method provides spatially isotropic discretizations for low and high Reynolds number flows. Various example solutions are given to illustrate the capabilities of the procedure. This article and been retracted and replaced. See retraction and replacement notice DOI: 10.1002/fld.1247 . Copyright © 2005 John Wiley & Sons, Ltd.     

A fractional two‐step implicit algorithm using curvilinear collected grids

An efficient fractional two‐step implicit algorithm is reported to simulate incompressible fluid flows in a boundary‐fitted curvilinear collocated grid system. Using the finite volume method, the convection terms are discretized by the high‐accuracy Roe's scheme to minimize numerical diffusion. An implicitness coefficient Π is introduced to accelerate the rate of convergence. It is demonstrated that the proposed algorithm links the fractional step method to the pressure correction procedure, and the SIMPLEC method could be considered as a special case of the fractional two‐step implicit algorithm (when Π=1). The proposed algorithm is applicable to unsteady flows and steady flows. Three benchmark two‐dimensional laminar flows are tested to evaluate the performance of the proposed algorithm. Performance is measured by sensitivity analyses of the efficiency, accuracy, grid density, grid skewness and Reynolds number on the solutions. Results show that the model is efficient and robust. Copyright © 2006 John Wiley & Sons, Ltd.     

High‐resolution,monotone solution of the adjoint shallow‐water equations

A monotone, second‐order accurate numerical scheme is presented for solving the differential form of the adjoint shallow‐water equations in generalized two‐dimensional coordinates. Fluctuation‐splitting is utilized to achieve a high‐resolution solution of the equations in primitive form. One‐step and two‐step schemes are presented and shown to achieve solutions of similarly high accuracy in one dimension. However, the two‐step method is shown to yield more accurate solutions to problems in which unsteady wave speeds are present. In two dimensions, the two‐step scheme is tested in the context of two parameter identification problems, and it is shown to accurately transmit the information needed to identify unknown forcing parameters based on measurements of the system response. The first problem involves the identification of an upstream flood hydrograph based on downstream depth measurements. The second problem involves the identification of a long wave state in the far‐field based on near‐field depth measurements. Copyright © 2002 John Wiley & Sons, Ltd.     

An efficient parallel high-order compact scheme for the 3D incompressible Navier–Stokes equations

This article provides a strategy for solving incompressible turbulent flows, which combines compact finite difference schemes and parallel computing. The numerical features of this solver are the semi-implicit time advancement, the staggered arrangement of the variables and the fourth-order compact scheme discretisation. This is the usual way for solving accurately turbulent incompressible flows. We propose a new strategy for solving the Helmholtz/Poisson equations based on a parallel 2d-pencil decomposition of the diagonalisation method. The compact scheme derivatives are computed with the parallel diagonal dominant (PDD) algorithm, which achieves good parallel performances by introducing a bounded numerical error. We provide a new analysis of its effect on the numerical accuracy and conservation features. Several numerical experiments, including two simulations of turbulent flows, demonstrate that the PDD algorithm maintains the accuracy and conservation features, while conserving a good parallel performance, up to 4096 cores.     

Fourth‐order compact formulation of Navier–Stokes equations and driven cavity flow at high Reynolds numbers

A new fourth‐order compact formulation for the steady 2‐D incompressible Navier–Stokes equations is presented. The formulation is in the same form of the Navier–Stokes equations such that any numerical method that solve the Navier–Stokes equations can easily be applied to this fourth‐order compact formulation. In particular, in this work the formulation is solved with an efficient numerical method that requires the solution of tridiagonal systems using a fine grid mesh of 601 × 601. Using this formulation, the steady 2‐D incompressible flow in a driven cavity is solved up to Reynolds number with Re = 20 000 fourth‐order spatial accuracy. Detailed solutions are presented. 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.     

Study of a staggered fourth‐order compact scheme for unsteady incompressible viscous flows

The objective of this paper is the development and assessment of a fourth‐order compact scheme for unsteady incompressible viscous flows. A brief review of the main developments of compact and high‐order schemes for incompressible flows is given. A numerical method is then presented for the simulation of unsteady incompressible flows based on fourth‐order compact discretization with physical boundary conditions implemented directly into the scheme. The equations are discretized on a staggered Cartesian non‐uniform grid and preserve a form of kinetic energy in the inviscid limit when a skew‐symmetric form of the convective terms is used. The accuracy and efficiency of the method are demonstrated in several inviscid and viscous flow problems. Results obtained with different combinations of second‐ and fourth‐order spatial discretizations and together with either the skew‐symmetric or divergence form of the convective term are compared. The performance of these schemes is further demonstrated by two challenging flow problems, linear instability in plane channel flow and a two‐dimensional dipole–wall interaction. Results show that the compact scheme is efficient and that the divergence and skew‐symmetric forms of the convective terms produce very similar results. In some but not all cases, a gain in accuracy and computational time is obtained with a high‐order discretization of only the convective and diffusive terms. Finally, the benefits of compact schemes with respect to second‐order schemes is discussed in the case of the fully developed turbulent channel flow. Copyright © 2008 John Wiley & Sons, Ltd.     

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.     

A COMPUTATIONAL METHOD FOR THE NAVIER-STOKES EQUATIONS AT ALL SPEEDS

ＩｎｔｒｏｄｕｃｔｉｏｎＦｏｒｃｏｍｐｕｔａｔｉｏｎｏｆｔｈｅｖｉｓｃｏｕｓｃｏｍｐｒｅｓｓｉｂｌｅｆｌｏｗｓｔｈｅＮａｖｉｅｒ＿Ｓｔｏｋｅｓｅｑｕａｔｉｏｎｓａｒｅｇｅｎｅｒａｌｌｙｒｅｐｒｅｓｅｎｔｅｄｉｎｔｈｅｃｏｎｓｅｒｖａｔｉｏｎｌａｗｆｏｒｍａｓａｈｙｐｅｒｂｏｌｉｃｓｙｓｔｅｍ .Ｌａｃｋｉｎｇｔｈｅｍａｔｈｅｍａｔｉｃａｌｔｏｏｌｔｏａｎａｌｙｚｅｔｈｉｓｎｏｎｌｉｎｅａｒｓｙｓｔｅｍ ,ｔｈｅｎｕｍｅｒｉｃａｌｍｅｔｈｏｄｓｕｓｅｄｉｎｓｏｌｖｉｎｇｔｈｅｎｏｎｌｉｎｅａｒｈｙｐｅ…     

A hybrid pressure‐based solver for nonideal single‐phase fluid flows at all speeds

This paper describes the implementation of a numerical solver that is capable of simulating compressible flows of nonideal single‐phase fluids. The proposed method can be applied to arbitrary equations of state and is suitable for all Mach numbers. The pressure‐based solver uses the operator‐splitting technique and is based on the PISO/SIMPLE algorithm: the density, velocity, and temperature fields are predicted by solving the linearized versions of the balance equations using the convective fluxes from the previous iteration or time step. The overall mass continuity is ensured by solving the pressure equation derived from the continuity equation, the momentum equation, and the equation of state. Nonphysical oscillations of the numerical solution near discontinuities are damped using the Kurganov‐Tadmor/Kurganov‐Noelle‐Petrova (KT/KNP) scheme for convective fluxes. The solver was validated using different test cases, where analytical and/or numerical solutions are present or can be derived: (1) A convergent‐divergent nozzle with three different operating conditions; (2) the Riemann problem for the Peng‐Robinson equation of state; (3) the Riemann problem for the covolume equation of state; (4) the development of a laminar velocity profile in a circular pipe (also known as Poiseuille flow); (5) a laminar flow over a circular cylinder; (6) a subsonic flow over a backward‐facing step at low Reynolds numbers; (7) a transonic flow over the RAE 2822 airfoil; and (8) a supersonic flow around a blunt cylinder‐flare model. The spatial approximation order of the scheme is second order. The mesh convergence of the numerical solution was achieved for all cases. The accuracy order for highly compressible flows with discontinuities is close to first order and, for incompressible viscous flows, it is close to second order. The proposed solver is named rhoPimpleCentralFoam and is implemented in the open‐source CFD library OpenFOAM^®. For high speed flows, it shows a similar behavior as the KT/KNP schemes (implemented as rhoCentralFoam‐solver, Int. J. Numer. Meth. Fluids 2010), and for flows with small Mach numbers, it behaves like solvers that are based on the PISO/SIMPLE algorithm.