@Article{CiCP-16-3, author = {}, title = {Stability of Projection Methods for Incompressible Flows Using High Order Pressure-Velocity Pairs of Same Degree: Continuous and Discontinuous Galerkin Formulations}, journal = {Communications in Computational Physics}, year = {2014}, volume = {16}, number = {3}, pages = {817--840}, abstract = {
This paper presents limits for stability of projection type schemes when
using high order pressure-velocity pairs of same degree. Two high order $h/p$ variational methods encompassing continuous and discontinuous Galerkin formulations
are used to explain previously observed lower limits on the time step for projection
type schemes to be stable [18], when $h$- or $p$-refinement strategies are considered. In
addition, the analysis included in this work shows that these stability limits depend not only on the time step but on the product of the latter and the kinematic viscosity, which is of particular importance in the study of high Reynolds number flows.
We show that high order methods prove advantageous in stabilising the simulations
when small time steps and low kinematic viscosities are used.
Drawing upon this analysis, we demonstrate how the effects of this instability can
be reduced in the discontinuous scheme by introducing a stabilisation term into the
global system. Finally, we show that these lower limits are compatible with Courant-Friedrichs-Lewy (CFL) type restrictions, given that a sufficiently high polynomial order or a mall enough mesh spacing is selected.