A symmetry-preserving second-order time-accurate PISO-based method

A new conservative symmetry-preserving second-order time-accurate PISO-based pressure-velocity coupling for solving the incompressible Navier-Stokes equations on unstructured collocated grids is presented in this paper. This new method for implicit time stepping is an extension of the conservative symmetry-preserving incremental-pressure projection method for explicit time stepping and unstructured collocated meshes of Trias et al. (2014). In order to assess and compare both methods, we have implemented them within one unified solver in the open source code OpenFOAM. We combine both methods with a Butcher tableau for a family of explicit and diagonally implicit Runge-Kutta temporal schemes. We assess the energy conservation properties of the implemented discretisation methods and the temporal consistency of the selected Runge-Kutta schemes using Taylor-Green vortex and lid-driven cavity flow test cases. Although both implemented methods are based on a symmetry-preserving discretisation, we show that both methods still produce a small amount of numerical dissipation when the total pressure is directly solved from a Poisson equation. This numerical dissipation is mainly caused by the corresponding pressure error which is of $O(\Delta t \Delta h^2)$. When an incremental-pressure approach is used, where a pressure correction is solved from a Poisson equation, the pressure error reduces to $O(\Delta t^2 \Delta h^2)$, yielding better conservation properties: both methods are then effectively fully-conservative. Furthermore, we conclude that all selected explicit and implicit higher order temporal schemes suffer from a reduction of the temporal order to approximately one when the pressure Poisson equation is based on the total pressure due to the presence of a pressure error of $O(\Delta t \Delta h^2)$.

[1]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[2]  Valerio D’Alessandro,et al.  Analysis of the numerical dissipation rate of different Runge-Kutta and velocity interpolation methods in an unstructured collocated finite volume method in OpenFOAM® , 2020, Comput. Phys. Commun..

[3]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[4]  Weeratunge Malalasekera,et al.  An introduction to computational fluid dynamics - the finite volume method , 2007 .

[5]  Shashank,et al.  A co-located incompressible Navier-Stokes solver with exact mass, momentum and kinetic energy conservation in the inviscid limit , 2010, J. Comput. Phys..

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

[7]  Jason E. Hicken,et al.  A shift transformation for fully conservative methods: turbulence simulation on complex, unstructured grids , 2005 .

[8]  A. Chorin Numerical solution of the Navier-Stokes equations , 1968 .

[9]  J. Butcher On Runge-Kutta processes of high order , 1964, Journal of the Australian Mathematical Society.

[10]  V. S. Vaidhyanathan,et al.  Transport phenomena , 2005, Experientia.

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

[12]  F. Krogh,et al.  Solving Ordinary Differential Equations , 2019, Programming for Computations - Python.

[13]  J. Kan A second-order accurate pressure correction scheme for viscous incompressible flow , 1986 .

[14]  A. Kouzmin,et al.  IT Development: , 1966, Current History.

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

[16]  John Kim,et al.  DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOWS UP TO RE=590 , 1999 .

[17]  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..

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

[19]  Barry Koren,et al.  Quasi-DNS capabilities of OpenFOAM for different mesh types , 2014 .

[20]  E. Hairer,et al.  Solving Ordinary Differential Equations II , 2010 .

[21]  Arthur Veldman,et al.  Direct Numerical Simulation of Driven Cavity Flows , 1993 .

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

[23]  F. Archambeau,et al.  Code Saturne: A Finite Volume Code for the computation of turbulent incompressible flows - Industrial Applications , 2004 .

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

[25]  H. Bijl,et al.  Analysis and application of high order implicit Runge–Kutta schemes to collocated finite volume discretization of the incompressible Navier–Stokes equations , 2015 .

[26]  Chun Shen,et al.  Implementation of density-based solver for all speeds in the framework of OpenFOAM , 2014, Comput. Phys. Commun..

[27]  A. Gosman,et al.  Solution of the implicitly discretised reacting flow equations by operator-splitting , 1986 .

[28]  Steven J. Ruuth,et al.  Implicit-explicit Runge-Kutta methods for time-dependent partial differential equations , 1997 .

[29]  Sergio Pirozzoli,et al.  A low-dissipative solver for turbulent compressible flows on unstructured meshes, with OpenFOAM implementation , 2017 .

[30]  Giacomo Castiglioni,et al.  A numerical dissipation rate and viscosity in flow simulations with realistic geometry using low-order compressible Navier–Stokes solvers , 2015 .

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

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

[33]  A. W. Vreman,et al.  Comparison of direct numerical simulation databases of turbulent channel flow at Re τ = 180 , 2014 .

[34]  Thomas S. Lund,et al.  Kinetic energy conservation issues associated with the collocated mesh scheme for incompressible flow , 2006, J. Comput. Phys..

[35]  F. Xavier Trias,et al.  Symmetry-preserving discretization of Navier-Stokes equations on collocated unstructured grids , 2014, J. Comput. Phys..

[36]  Thilo Knacke,et al.  Improving LES with OpenFOAM by minimising numerical dissipation and use of explicit algebraic SGS stress model , 2019 .

[37]  Bernard J. Geurts,et al.  A quantification method for numerical dissipation in quasi-DNS and under-resolved DNS, and effects of numerical dissipation in quasi-DNS and under-resolved DNS of turbulent channel flows , 2017, J. Comput. Phys..

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

[39]  S. Castaño,et al.  Assessment of Solution Algorithms for LES of Turbulent Flows Using OpenFOAM , 2019, Fluids.