Artificial compressibility method revisited: Asymptotic numerical method for incompressible Navier–Stokes equations

The artificial compressibility method for the incompressible Navier–Stokes equations is revived as a high order accurate numerical method (fourth order in space and second order in time). Similar to the lattice Boltzmann method, the mesh spacing is linked to the Mach number. An accuracy higher than that of the lattice Boltzmann method is achieved by exploiting the asymptotic behavior of the solution of the artificial compressibility equations for small Mach numbers and the simple lattice structure. An easy method for accelerating the decay of acoustic waves, which deteriorate the quality of the numerical solution, and a simple cure for the checkerboard instability are proposed. The high performance of the scheme is demonstrated not only for the periodic boundary condition but also for the Dirichlet-type boundary condition.