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.

[1]  Marcelo H. Kobayashi,et al.  A fourth-order-accurate finite volume compact method for the incompressible Navier-Stokes solutions , 2001 .

[2]  F. Grinstein,et al.  Recent Progress on MILES for High Reynolds Number Flows , 2002 .

[3]  John A. Ekaterinaris,et al.  A staggered grid, high-order accurate method for the incompressible Navier-Stokes equations , 2006, J. Comput. Phys..

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

[5]  Eric Vedy,et al.  A high‐order finite difference method for incompressible fluid turbulence simulations , 2003 .

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

[7]  R. Moser,et al.  Spectral methods for the Navier-Stokes equations with one infinite and two periodic directions , 1991 .

[8]  C. Pozrikidis A note on the regularization of the discrete Poisson-Neumann problem , 2001 .

[9]  O. Vasilyev,et al.  Fully conservative finite difference scheme in cylindrical coordinates for incompressible flow simulations , 2004 .

[10]  P. Moin,et al.  Fully Conservative Higher Order Finite Difference Schemes for Incompressible Flow , 1998 .

[11]  M. Shoucri,et al.  Solution of the Poisson Equation in an Annulus , 1997 .

[12]  Parviz Moin,et al.  An improvement of fractional step methods for the incompressible Navier-Stokes equations , 1991 .

[13]  P. Moin,et al.  Suitability of upwind-biased finite difference schemes for large-eddy simulation of turbulent flows , 1997 .

[14]  H. Fasel,et al.  A Compact-Difference Scheme for the Navier—Stokes Equations in Vorticity—Velocity Formulation , 2000 .

[15]  P. Moin,et al.  On the Effect of Numerical Errors in Large Eddy Simulations of Turbulent Flows , 1997 .

[16]  Toshio Kobayashi,et al.  Numerical solution of the incompressible Navier-Stokes equations with an upwind compact difference scheme , 1999 .

[17]  Stéphane Abide,et al.  A 2D compact fourth-order projection decomposition method , 2005 .

[18]  M. Piller,et al.  Finite-volume compact schemes on staggered grids , 2004 .

[19]  High‐order upwind compact scheme and simulation of turbulent premixed V‐flame , 2005 .

[20]  Roland Schiestel,et al.  A hermitian-fourier numerical method for solving the incompressible navier-stokes equations , 1995 .

[21]  S. Orszag Accurate solution of the Orr–Sommerfeld stability equation , 1971, Journal of Fluid Mechanics.

[22]  Joel H. Ferziger,et al.  A robust high-order compact method for large eddy simulation , 2003 .

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

[24]  Johan Meyers,et al.  A computational error-assessment of central finite-volume discretizations in large-eddy simulation using a Smagorinsky model , 2007, J. Comput. Phys..

[25]  P. Moin,et al.  The minimal flow unit in near-wall turbulence , 1991, Journal of Fluid Mechanics.

[26]  Miguel R. Visbal,et al.  Turbulence spectra characteristics of high order schemes for direct and large eddy simulation , 2001 .

[27]  Joel H. Ferziger,et al.  Improved Methods for Large Eddy Simulations of Turbulence , 1979 .

[28]  P. Swarztrauber A direct Method for the Discrete Solution of Separable Elliptic Equations , 1974 .

[29]  A. B. Strong,et al.  A Fully Conservative Second-Order Finite Difference Scheme for Incompressible Flow on Nonuniform Grids , 2002 .

[30]  Per Lötstedt,et al.  High order accurate solution of the incompressible Navier-Stokes equations , 2005 .

[31]  Wendy Kress,et al.  A Compact Higher Order Finite Difference Method for the Incompressible Navier–Stokes Equations , 2002, J. Sci. Comput..

[32]  Ulrich Rüde,et al.  Cache-Aware Multigrid Methods for Solving Poisson's Equation in Two Dimensions , 2000, Computing.

[33]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[34]  Hjh Herman Clercx,et al.  Benchmark computations of the normal and oblique collision of a dipole with a no-slip boundary , 2006 .

[35]  A. Veldman,et al.  Symmetry-preserving discretization of turbulent flow , 2003 .

[36]  Towards the ultimate variance-conserving convection scheme , 2004 .

[37]  Franck Nicoud,et al.  Conservative High-Order Finite-Difference Schemes for Low-Mach Number Flows , 2000 .

[38]  P. Moin,et al.  Turbulence statistics in fully developed channel flow at low Reynolds number , 1987, Journal of Fluid Mechanics.

[39]  Mahidhar Tatineni,et al.  Very high-order compact finite difference schemes on non-uniform grids for incompressible Navier-Stokes equations , 2007, J. Comput. Phys..

[40]  F. Nicoud,et al.  Compact finite difference schemes on non-uniform meshes. Application to direct numerical simulations of compressible flows , 1999 .

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

[42]  Koji Fukagata,et al.  Highly energy-conservative finite difference method for the cylindrical coordinate system , 2002 .