A generalised and low-dissipative multi-directional characteristics-based scheme with inclusion of the local Riemann problem investigating incompressible flows without free-surfaces

摘要

In the present study, we develop a generalised Godunov-type multi-directional characteristics-based (MCB) scheme which is applicable to any hyperbolic system for modelling incompressible flows. We further extend the MCB scheme to include the solution of the local Riemann problem which leads to a hybrid mathematical treatment of the system of equations. We employ the proposed scheme to hyperbolic-type incompressible flow solvers and apply it to the Artificial Compressibility (AC) and Fractional-Step, Artificial Compressibility with Pressure Projection (FSAC-PP) method. In this work, we show that the MCB scheme may improve the accuracy and convergence properties over the classical single-directional characteristics-based (SCB) and non-characteristic treatments. The inclusion of a Riemann solver in conjunction with the MCB scheme is capable of reducing the number of iterations up to a factor of 4.7 times compared to a solution when a Riemann solver is not included. Furthermore, we found that both the AC and FSAC-PP method showed similar levels of accuracy while the FSAC-PP method converged up to 5.8 times faster than the AC method for steady state flows. Independent of the characteristics and Riemann solver-based treatment of all primitive variables, we found that the FSAC-PP method is 7-200 times faster than the AC method per pseudo-time step for unsteady flows. We investigate low- and high-Reynolds number problems for well-established validation benchmark test cases focusing on a flow inside of a lid driven cavity, evolution of the Taylor-Green vortex and forced separated flow over a backward-facing step. In addition to this, comparisons between a central difference scheme with artificial dissipation and a low-dissipative interpolation scheme have been performed. The results show that the latter approach may not provide enough numerical dissipation to develop the flow at high-Reynolds numbers. We found that the inclusion of a Riemann solver is able to overcome this s