We propose a simple and efficient finite-difference method for the incompressible MHD equation. The numerical method combines the advantage of the MAC scheme for the Navier–Stokes equation and Yee's scheme for the Maxwell equation. In particular, the semi-discrete version of our scheme introduces no numerical dissipation and preserves the energy identity exactly.

A very simple and efficient finite element method is introduced for two and three dimensional viscous incompressible flows using the vorticity formulation. This method relies on recasting the traditional finite element method in the spirit of the high order accurate finite difference methods introduced by the authors in another work. Optimal accuracy of arbitrary order can be achieved using standard finite element or spectral elements. The method is convectively stable and is particularly suited for moderate to high Reynolds number flows.

In this paper, we provide stability and convergence analysis for a class of finite difference schemes for unsteady incompressible Navier-Stokes equations in vorticity-stream function formulation. The no-slip boundary condition for the velocity is converted into local vorticity boundary conditions. Thom's formula, Wilkes' formula, or other local formulas in the earlier literature can be used in the second order method; while high order formulas, such as Briley's formula, can be used in the fourth order compact difference scheme proposed by E and Liu. The stability analysis of these long-stencil formulas cannot be directly derived from straightforward manipulations since more than one interior point is involved in the formula. The main idea of the stability analysis is to control local terms by global quantities via discrete elliptic regularity for stream function. We choose to analyze the second order scheme with Wilkes' formula in detail. In this case, we can avoid the complicated technique necessitated by the Strang-type high order expansions. As a consequence, our analysis results in almost optimal regularity assumption for the exact solution. The above methodology is very general. We also give a detailed analysis for the fourth order scheme using a 1-D Stokes model.

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.