In this work, we deal with high-order solver for incompressible flow based on velocity correction scheme with discontinuous Galerkin discretized velocity and standard continuous approximated pressure. Recently, small time step instabilities have been reported for pure discontinuous Galerkin method, in which both velocity and pressure are discretized by discontinuous Galerkin. It is interesting to examine these instabilities in the context of mixed discontinuous Galerkin–continuous Galerkin method. By means of numerical investigation, we find that the discontinuous Galerkin–continuous Galerkin method shows great stability at the same configuration. The consistent velocity divergence discretization scheme helps to achieve more accurate results at small time step size. Since the equal order discontinuous Galerkin–continuous Galerkin method does not satisfy inf-sup stability requirement, the instability for high Reynolds number flow is investigated. We numerically demonstrate that fine mesh resolution and high polynomial order are required to obtain a robust system. With these conclusions, discontinuous Galerkin–continuous Galerkin method is able to achieve high-order spatial convergence rate and accurately simulate high Reynolds flow. The solver is tested through a series of classical benchmark problems, and efficiency improvement is proved against pure discontinuous Galerkin scheme.
A high-order semi-explicit discontinuous Galerkin solver for 3D incompressible flow with application to DNS and LES of turbulent channel flow