On the accuracy of stiff-accurate diagonal implicit Runge-Kutta methods for finite volume based Navier-Stokes equations

The paper aims at developing low-storage implicit Runge-Kutta methods which are easy to implement and achieve higher-order of convergence for both the velocity and pressure in the finite volume formulation of the incompressible Navier-Stokes equations on a static collocated grid. To this end, the effect of the momentum interpolation, a procedure required by the finite volume method for collocated grids, on the differential-algebraic nature of the spatially-discretized Navier-Stokes equations should be examined first. A new framework for the momentum interpolation is established, based on which the semi-discrete Navier-Stokes equations can be strictly viewed as a system of differential-algebraic equations of index 2. The accuracy and convergence of the proposed momentum interpolation framework is examined. We then propose a new method of applying implicit Runge-Kutta schemes to the time-marching of the index 2 system of the incompressible Navier-Stokes equations. Compared to the standard method, the proposed one significantly reduces the numerical difficulties in momentum interpolations and delivers higher-order pressures without requiring additional computational effort. Applying stiff-accurate diagonal implicit Runge-Kutta (DIRK) schemes with the proposed method allows the schemes to attain the classical order of convergence for both the velocity and pressure. We also develop two families of low-storage stiff-accurate DIRK schemes to reduce the storage required by their implementations. Examining the two dimensional Taylor-Green vortex as an example, the spatial and temporal accuracy of the proposed methods in simulating incompressible flow is demonstrated.

[1]  P. Houwen.  Explicit Runge-Kutta formulas with increased stability boundaries , 1972 .

[2]  B. P. Leonard,et al.  A stable and accurate convective modelling procedure based on quadratic upstream interpolation , 1990 .

[3]  S. G. Rubin,et al.  A diagonally dominant second-order accurate implicit scheme , 1974 .

[4]  A. Murua,et al.  Partitioned half-explicit Runge-Kutta methods for differential-algebraic systems of index 2 , 1997, Computing.

[5]  C. Rhie,et al.  Numerical Study of the Turbulent Flow Past an Airfoil with Trailing Edge Separation , 1983 .

[6]  Santiago Badia,et al.  Segregated Runge–Kutta methods for the incompressible Navier–Stokes equations , 2016 .

[7]  Shashi M. Aithal,et al.  Influence of momentum interpolation methods on the accuracy and convergence of pressure-velocity coupling algorithms in OpenFOAM® , 2017, J. Comput. Appl. Math..

[8]  Barry Koren,et al.  Accuracy analysis of explicit Runge-Kutta methods applied to the incompressible Navier-Stokes equations , 2012, J. Comput. Phys..

[9]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[10]  James M. Ortega,et al.  Iterative solution of nonlinear equations in several variables , 2014, Computer science and applied mathematics.

[11]  M. Darwish,et al.  Erratum to: The Finite Volume Method in Computational Fluid Dynamics , 2016 .

[12]  Laurent O. Jay,et al.  Convergence of a class of runge-kutta methods for differential-algebraic systems of index 2 , 1993 .

[13]  Ernst Hairer,et al.  The numerical solution of differential-algebraic systems by Runge-Kutta methods , 1989 .

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

[15]  R. Lewis,et al.  Low-storage, Explicit Runge-Kutta Schemes for the Compressible Navier-Stokes Equations , 2000 .

[16]  J. Butcher.  Implicit Runge-Kutta processes , 1964 .

[17]  Sijun Zhang,et al.  Generalized formulations for the Rhie-Chow interpolation , 2014, J. Comput. Phys..

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

[19]  Seok-Ki Choi,et al.  Note on the use of momentum interpolation method for unsteady flows , 1999 .

[20]  Christopher A. Kennedy,et al.  Diagonally implicit Runge–Kutta methods for stiff ODEs , 2019 .

[21]  L. Kantorovich,et al.  Functional analysis in normed spaces , 1952 .

[22]  Christophe Duwig,et al.  On the implementation of low-dissipative Runge–Kutta projection methods for time dependent flows using OpenFOAM® , 2014 .

[23]  J. Fromm.  A method for reducing dispersion in convective difference schemes , 1968 .

[24]  H. Jasak,et al.  Consistent second-order time-accurate non-iterative PISO-algorithm , 2018 .

[25]  Nikolay Nikitin,et al.  Third‐order‐accurate semi‐implicit Runge–Kutta scheme for incompressible Navier–Stokes equations , 2006 .

[26]  Bo Yu,et al.  DISCUSSION ON MOMENTUM INTERPOLATION METHOD FOR COLLOCATED GRIDS OF INCOMPRESSIBLE FLOW , 2002 .

[27]  P. Moin,et al.  Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations , 1984 .

[28]  E. Hairer,et al.  Half-explicit Runge-Kutta methods for differential-algebraic systems of index 2 , 1993 .

[29]  Valerio D'Alessandro,et al.  On the development of OpenFOAM solvers based on explicit and implicit high-order Runge-Kutta schemes for incompressible flows with heat transfer , 2018, Comput. Phys. Commun..

[30]  A. Pascau,et al.  Cell face velocity alternatives in a structured colocated grid for the unsteady Navier–Stokes equations , 2011 .

[31]  Daniele Cavaglieri,et al.  Low-storage implicit/explicit Runge-Kutta schemes for the simulation of stiff high-dimensional ODE systems , 2015, J. Comput. Phys..

[32]  X. Systemsofinde,et al.  Convergence of a Class of Runge-kutta Methods for Differential-algebraic Systems of Index 2 , 2005 .