An implicit high order finite volume scheme for the solution of 3D Navier–Stokes equations with new discretization of diffusive terms

The computation of three-dimensional viscous flows on complex geometries requiring distorted meshes is of great interest. This paper presents a finite volume solver using a quadratic reconstruction of the unknowns for the advective fluxes computation, and a conservative and consistent discretization of the diffusive terms, based on an extended version of the Coirier's diamond path. A fully implicit time integration procedure is employed with a preconditioned matrix-free GMRES solver.