Stability analysis for different formulations of the nonlinear term in PN-PN-2 spectral element discretizations of the Navier-Stokes equations

We show that for the P-N - PN-2 spectral element method. in which velocity and pressure are approximated by polynomials of order N and N - 2, respecti vely, numerical instabilities can occur in the spatially discretized Navier -Stokes equations. Both a staggered and nonstaggered arrangement of the N - 2 pressure points are considered. These instabilities can be masked by vis cous damping at low Reynolds numbers. We demonstrate that the instabilities depend on the formulation of the nonlinear term. The numerical discretizat ion is stable for the convective form but unstable for the divergence and t he skew-symmetric form. Further numerical analysis indicates that this inst ability is not caused by nonlinear effects. since it occurs for linearized systems as well. An eigenvalue analysis of the fully discretized shows that an instability is introduced by the formulation of the nonlinear term. We demonstrate that the instability is related to the divergence error of the computed solution at those velocity points at which the continuity equation is not enforced. (C) 2001 Elsevier Science.