Artificial Compressibility Method for the Incompressible Navier–Stokes Equations With Variable Density

We introduce a novel artificial compressibility technique to approximate the incompressible Navier–Stokes equations with variable fluid properties such as density and dynamical viscosity. The proposed method uses the couple pressure and momentum, equal to the density times the velocity, as primary unknowns. It also involves an adequate treatment of the diffusive operator such that treating the nonlinear convective term explicitly leads to a method with time‐independent stiffness matrices that is suitable for pseudo‐spectral methods. The stability and temporal convergence of a semi‐implicit version of the method are established under the hypothesis that the density is approximated with a method that conserves the minimum‐maximum principle. Numerical illustrations confirm that both the semi‐implicit and explicit methods are stable and converge with order one under the classic CFL condition. Moreover, the proposed method is shown to perform better than a momentum‐based pressure projection method, previously introduced by one of the authors, on setups involving gravitational waves and immiscible multi‐fluids in a cylinder.