A fourth-order accurate compact scheme for the solution of steady Navier–Stokes equations on non-uniform grids

Abstract This paper deals with the formulation of a higher-order compact (HOC) scheme on non-uniform grids in complex geometries to simulate two-dimensional (2D) steady incompressible viscous flows governed by the Navier–Stokes (N–S) equations. The proposed scheme which is spatially fourth-order accurate is then tested on three nonlinear problems, namely (i) a problem governed by N–S equations with a constructed analytical solution, (ii) lid-driven cavity flow problem, and (iii) constricted channel flow problem. In the process, we have also expanded the scope of fourth-order 9-point compact schemes to geometries beyond rectangular. It is seen to efficiently capture steady-state solutions of the N–S equations with Dirichlet as well as Neumann boundary conditions. In addition to this, it captures viscous flows involving free and wall bounded shear layers which invariably contain spatial scale variations. Our results are in excellent agreement with analytical and numerical results whenever available and they clearly demonstrate the superior scale resolution of the proposed scheme.

[1]  M. D. Deshpande,et al.  FLUID MECHANICS IN THE DRIVEN CAVITY , 2000 .

[2]  C. Kelley.  Iterative Methods for Linear and Nonlinear Equations , 1987 .

[3]  Murli M. Gupta,et al.  A single cell high order scheme for the convection‐diffusion equation with variable coefficients , 1984 .

[4]  P. R. Garabedian,et al.  Estimation of the relaxation factor for small mesh size , 1956 .

[5]  Y. Adam,et al.  Highly accurate compact implicit methods and boundary conditions , 1977 .

[6]  J. C. Kalita,et al.  Fully compact higher-order computation of steady-state natural convection in a square cavity. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[7]  Roland Hunt,et al.  Fourth‒order method for solving the Navier–Stokes equations in a constricting channel , 1997 .

[8]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[9]  J. C. Kalita,et al.  A transformation‐free HOC scheme for steady convection–diffusion on non‐uniform grids , 2004 .

[10]  Ayodeji O. Demuren,et al.  Higher-Order Compact Schemes for Numerical Simulation of Incompressible Flows , 1998 .

[11]  G. Carey,et al.  High‐order compact scheme for the steady stream‐function vorticity equations , 1995 .

[12]  William D. Henshaw,et al.  A Fourth-Order Accurate Method for the Incompressible Navier-Stokes Equations on Overlapping Grids , 1994 .

[13]  Tapan K. Sengupta,et al.  Analysis of central and upwind compact schemes , 2003 .

[14]  Z. Kopal,et al.  Numerical analysis , 1955 .

[15]  C. Bruneau,et al.  The 2D lid-driven cavity problem revisited , 2006 .

[16]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[17]  Jie Wang,et al.  High order compact computation and nonuniform grids for streamfunction vorticity equations , 2006, Appl. Math. Comput..

[18]  Paulo F. A. Mancera,et al.  A study of a numerical solution of the steady two dimensions Navier-Stokes equations in a constricted channel problem by a compact fourth order method , 2003, Appl. Math. Comput..

[19]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[20]  O. Botella,et al.  BENCHMARK SPECTRAL RESULTS ON THE LID-DRIVEN CAVITY FLOW , 1998 .

[21]  Robert L. Street,et al.  The Lid-Driven Cavity Flow: A Synthesis of Qualitative and Quantitative Observations , 1984 .

[22]  B. Fornberg,et al.  A compact fourth‐order finite difference scheme for the steady incompressible Navier‐Stokes equations , 1995 .

[23]  Luigi Vigevano,et al.  Accurate ω-ψ spectral solution of the singular driven cavity problem , 2002 .

[24]  A. I. van de Vooren,et al.  On the 9-point difference formula for Laplace's equation , 1967 .

[25]  Scott E. Sherer,et al.  High-order compact finite-difference methods on general overset grids , 2005 .

[26]  John C. Strikwerda.  HIGH‐ORDER‐ACCURATE SCHEMES FOR INCOMPRESSIBLE VISCOUS FLOW , 1997 .

[27]  Jun Zhang,et al.  High accuracy iterative solution of convection diffusion equation with boundary layers on nonuniform grids , 2001 .

[28]  Tao Tang,et al.  A Compact Fourth-Order Finite Difference Scheme for Unsteady Viscous Incompressible Flows , 2001, J. Sci. Comput..

[29]  H. K. Moffatt.  Viscous and resistive eddies near a sharp corner , 1964, Journal of Fluid Mechanics.

[30]  Murli M. Gupta.  High accuracy solutions of incompressible Navier-Stokes equations , 1991 .

[31]  W. Spotz.  Formulation and experiments with high‐order compact schemes for nonuniform grids , 1998 .

[32]  H. B. Keller,et al.  Driven cavity flows by efficient numerical techniques , 1983 .

[33]  Miguel R. Visbal,et al.  On the use of higher-order finite-difference schemes on curvilinear and deforming meshes , 2002 .

[34]  Charles-Henri Bruneau,et al.  An efficient scheme for solving steady incompressible Navier-Stokes equations , 1990 .

[35]  P. Roache.  Fundamentals of computational fluid dynamics , 1998 .

[36]  E. Barragy,et al.  STREAM FUNCTION-VORTICITY DRIVEN CAVITY SOLUTION USING p FINITE ELEMENTS , 1997 .

[37]  Shlomo Ta'asan,et al.  Finite difference schemes for long-time integration , 1994 .

[38]  Jiten C. Kalita,et al.  A class of higher order compact schemes for the unsteady two‐dimensional convection–diffusion equation with variable convection coefficients , 2002 .

[39]  Ayodeji O. Demuren,et al.  HIGHER-ORDER COMPACT SCHEMES FOR NUMERICAL SIMULATION OF INCOMPRESSIBLE FLOWS, PART I: THEORETICAL DEVELOPMENT , 2001 .

[40]  R. Pletcher,et al.  Computational Fluid Mechanics and Heat Transfer. By D. A ANDERSON, J. C. TANNEHILL and R. H. PLETCHER. Hemisphere, 1984. 599 pp. $39.95. , 1986, Journal of Fluid Mechanics.

[41]  Jiten C. Kalita,et al.  Erratum: Fully compact higher-order computation of steady-state natural convection in a square cavity [Phys. Rev. E 64, 066703 (2001)] , 2002 .

[42]  R. Hirsh,et al.  Higher order accurate difference solutions of fluid mechanics problems by a compact differencing technique , 1975 .

[43]  Robert J. MacKinnon,et al.  Differential‐equation‐based representation of truncation errors for accurate numerical simulation , 1991 .