Unconditionally stable, linearised IMEX schemes for incompressible flows with variable density

For the incompressible Navier--Stokes system with variable density and viscosity, we propose and analyse an IMEX framework treating the convective and diffusive terms semi-implicitly. This extends to variable density and second order in time some methods previously analysed for variable viscosity and constant density. We present three new schemes, both monolithic and fractional-step. All of them share the methodological novelty that the viscous term is treated in an implicit-explicit (IMEX) fashion, which allows decoupling the velocity components. Unconditional temporal stability is proved for all three variants. Furthermore, the system to solve at each time step is linear, thus avoiding the costly solution of nonlinear problems even if the viscosity follows a non-Newtonian rheological law. Our presentation is restricted to the semi-discrete case, only considering the time discretisation. In this way, the results herein can be applied to any spatial discretisation. We validate our theory through numerical experiments considering finite element methods in space. The tests range from simple manufactured solutions to complex two-phase viscoplastic flows.
Comment: 23 pages