A preconditioning for the spectral solution of incompressible variable-density flows

L. Reynier (1) , Bastien Di Pierro (1) , Frédéric Alizard (1) , Anne Cadiou (1) , Lionel Le Penven (1) , Marc Buffat (1)

1 LMFA - Laboratoire de Mecanique des Fluides et d'Acoustique

In the present study, the efficiency of preconditioners for solving linear systems associated with the discretized variable-density incompressible Navier-Stokes equations with semiimplicit second-order accuracy in time and spectral accuracy in space is investigated. The method, in which the inverse operator for the constant-density flow system acts as preconditioner, is implemented for three iterative solvers: the General Minimal Residual, the Conjugate Gradient and the Richardson Minimal Residual. We discuss the method, first, in the context of the one-dimensional flow case where a top-hat like profile for the density is used. Numerical evidence shows that the convergence is significantly improved due to the notable decrease in the condition number of the operators. Most importantly, we then validate the robustness and convergence properties of the method on two more realistic problems: the two-dimensional Rayleigh-Taylor instability problem and the three-dimensional variable-density swirling jet.

Keywords: Variable-density flows; incompressible Navier-Stokes equations; direct numerical simulation; preconditioning; elliptic solver

Domains: Fluid Dynamics [physics.flu-dyn]

L. Reynier, Bastien Di Pierro, Frédéric Alizard, Anne Cadiou, Lionel Le Penven, et al.. A preconditioning for the spectral solution of incompressible variable-density flows. Computers and Fluids, In press, 266, pp.106024. ⟨10.1016/j.compfluid.2023.106024⟩. ⟨hal-04183956v2⟩