A conservative, second order, unconditionally stable artificial compression method

This report presents a new artificial compression method for incompressible, viscous flows. The method has second order consistency error and is unconditionally, long time, energy stable for the velocity and, weighted by the timestep, for the pressure. It uncouples the pressure and velocity and requires no artificial pressure boundary conditions. When the viscosity. = 0 the method also exactly conserves a system energy. The method is based on a Crank-Nicolson Leapfrog time discretization of the slightly compressible model (1-epsilon(1) grad div)u(t) + u .del u + 1/2(div u)u - v Delta + del p = f and epsilon(2pt) + div u = 0. This report presents the method, gives a stability analysis, presents numerical tests and gives a preliminary analysis with tests of the non-physical acoustic waves generated. Consideration of the physical fidelity of the artificial compression method leads to a related method. (C) 2017 Elsevier B.V. All rights reserved.