In E & Liu (SIAM J Numer. Anal., 1995), we studied convergence and the structure of the error for several projection methods when the spatial variable was kept continuous (we call this the semi-discrete case). In this paper, we address similar questions for the fully discrete case when the spatial variables are discretized using a staggered grid. We prove that the numerical solution in velocity has full accuracy up to the boundary, despite the fact that there are numerical boundary layers present in the semi-discrete solutions.

In this paper we discuss the derivation and use of local pressure boundary conditions for finite difference schemes for the unsteady incompressible Navier-Stokes equations in the velocity-pressure formulation. Their use is especially well suited for the computation of moderate to large Reynolds number flows. We explore the similarities between the implementation and use of local pressure boundary conditions and local vorticity boundary conditions in the design of numerical schemes for incompressible flow in 2D. In their respective formulations, when these local numerical boundary conditions are coupled with a fully explicit convectively stable time stepping procedure, the resulting methods are simple to implement and highly efficient. Unlike the vorticity formulation, the use of the local pressure boundary condition approach is readily applicable to 3D flows. The simplicity of the local pressure boundary condition approach and its easy adaptation to more general flow settings make the resulting scheme an attractive alternative to the more popular methods for solving the Navier-Stokes equations in the velocity-pressure formulation. We present numerical results of a second-order finite difference scheme on a nonstaggered grid using local pressure boundary conditions. Stability and accuracy of the scheme applied to Stokes flow is demonstrated using normal mode analysis. Also described is the extension of the method to variable density flows.

We study a fourth order finite difference method for the unsteady incompressible Navier-Stokes equations in vorticity formulation. The scheme is essentially compact and can be implemented very efficiently. Either Briley’s formula, or a new higher order formula, which will be derived in this paper, can be chosen as the vorticity boundary condition. By formal Taylor expansion, the new formula for the vorticity on the boundary gives 4th order accuracy; while Briley’s formula provides only 3rd order accuracy. However, the use of either formula results in a stable method and achieves full 4th order accuracy. The convergence analysis of the scheme with our new formula will be given in this paper, while that with Briley’s formula has been established in earlier literature. The consistency analysis is easier than that of Briley’s formula, no Strang type analysis is needed. In the stability analysis part, we adopt the technique of controlling some local terms by the diffusion term via discrete elliptic regularity. Physical no-slip boundary conditions are used throughout.

A class of upwind flux splitting methods in the Euler equations of compressible flow is considered in this paper. Using the property that Euler flux F(U) is a homogeneous function of degree one in U, we reformulate the splitting fluxes with F^{+}=A^{+} U, F^{-}=A^{-} U, and the corresponding matrices are either symmetric or symmetrizable and keep only non-negative and non-positive eigenvalues. That leads to the conclusion that the first order schemes are positive in the sense of Lax-Liu [18], which implies that it is L^2-stable in some suitable sense. Moreover, the second order scheme is a stable perturbation of the first order scheme, so that the positivity of the second order schemes is also established, under a CFL-like condition. In addition, these splitting methods preserve the positivity of density and energy.

A new class of implicit high-order non-oscillatory time integration schemes is introduced in a method-of-lines framework. These schemes can be used in conjunction with an appropriate spatial discretization scheme for the numerical solution of time dependent conservation equations. The main concept behind these schemes is that the order of accuracy in time is dropped locally in regions where the time evolution of the solution is not smooth. By doing this, an attempt is made at locally satisfying monotonicity conditions, while maintaining a high order of accuracy in most of the solution domain. When a linear high order time integration scheme is used along with a high order spatial discretization, enforcement of monotonicity imposes severe time-step restrictions. We propose to apply limiters to these time-integration schemes, thus making them non-linear. When these new schemes are used with high order spatial discretizations, solutions remain non-oscillatory for much larger time-steps as compared to linear time integration schemes. Numerical results obtained on scalar conservation equations and systems of conservation equations are highly promising.