In this paper we study the performances of generalized minimal residual method (GMRES) preconditioned with geometric multigrid (GMG), applied to steady and unsteady buoyancy driven flow problems, discretized with the finite element method. For the unsteady case, the second order Crank–Nicolson method is used for the temporal discretization. At each geometric multigrid level, we use Richardson iterative solvers preconditioned with different combinations of physics-based and domain decomposition preconditioners. Three different preconditioners are considered: incomplete LU decomposition (ILU), overlapping Vanka-type domain decomposition for additive Schwarz method (ASM), and field split (FS) physics-based decomposition. We also analyze the effect on the smoother of how the variables are ordered, and in particular whether the leading variable is the velocity or the temperature, resulting in six classes of preconditioners: ILU_VT, ILU_TV, ASM_VT, ASM_TV, FS_VT and FS_TV. The eigenvalue analysis for the six preconditioners is conducted to study the rate of GMRES convergence under several Prandtl numbers. The numerical performances of nested combinations of the above preconditioners are compared. Numerical results show that the pair of FS_VT and FS_TV preconditioners works better than the other two pairs, and that the FS_TV preconditioner always performs the best in terms of the computational time for all the steady and unsteady cases.
- Buoyancy driven flow problems
- Domain decomposition and physics-based decomposition preconditioners
- Eigenvalue analysis