摘要

This work introduces a novel pressure-correction (PC) projection finite element method for solving the incompressible natural convection problem with variable density. The main novel ideas of this work are to overcome the stability issue due to the nonlinear inertial term and the hyperbolic term for conventional finite element methods. We choose the conforming finite element pair (P-2, P-2, P-1, P-2) to approximate the density, the velocity, the pressure, and the temperature field, respectively. Moreover, the stability analysis of the first- order PC projection finite element method is given. Finally, ample numerical results are presented to demonstrate the efficacy and accuracy of the proposed method.