A second-order projection method for the incompressible Navier-Stokes equations in arbitrary domains (Q1339528)

We present a projection method for the numerical solution of the incompressible Navier-Stokes equations in an arbitrary domain that is second-order accurate in both space and time. We use the continuity equation to derive a pressure equation to compute the gradient part of the vector field. An integral form of the continuity equation is used to give us a natural way to define the discrete divergence operator for cells near the boundary which ensures the diagonal dominance of the resulting pressure equation. We then use the restarted version of the GMRES method to solve the pressure equation.