A Hybrid Reconstructed Discontinuous Galerkin and Continuous Galerkin Finite Element Method for Incompressible Flows on Unstructured Grids

摘要

A hybrid reconstructed discontinuous Galerkin and continuous Galerkin method based on an incremental pressure projection formulation, termed rDG(P_nP_m)+CG(P_n) in this paper, is developed for solving the unsteady incompressible Navier-Stokes equations on unstructured grids. In this method, a reconstructed discontinuous Galerkin method (rDG(P_nP_m)) is used to discretize the velocity and a standard continuous Galerkin method (CG(P_n)) is used to approximate the pressure. The rDG(P_nP_m)-CG(P_n) method is designed to increase the accuracy of the hybrid DG(P_n)-CG(P_n) method and yet still satisfy Ladyzenskaja-Babuska-Brezzi (LBB) condition, thus avoiding the pressure checkerboard instability. An upwind method is used to discretize the nonlinear convective fluxes in the momentum equations in order to suppress spurious oscillations in the velocity field. A number of incompressible flow problems for a variety of flow conditions are computed to numerically assess the spatial order of convergence of the rDG(P_nP_m)+CG(P_m) method. The numerical experiments indicate that both rDG(P_0P_1)+CG(P_1) and rDG(P_1P_2)+CG(P_1) methods can attain the designed second order and third order accuracy in space for the velocity, respectively and the third order rDG(P_1P_2)+CG(P_1) method significantly outperforms its second order rDG(P_1P_1)+CG(P_1) and rDG(P_1P_1)+CG(P_1) counterparts: being able to not only increase the accuracy of the velocity by one order but also improve the accuracy of the pressure.