The Euler implicit/explicit scheme for the Boussinesq equations

In this article, we consider the stability and convergence of the first-order implicit/explicit scheme for the Boussinesq equations. The finite element spatial discretization is based on a MINI element for the velocity and pressure, which satisfies the discrete inf-sup condition, and a linear polynomial for the temperature. The temporal terms are treated by the Euler implicit/explicit scheme, which is implicit for the linear terms and explicit for the nonlinear terms. The advantage of using the implicit/explicit scheme is that a linear system with constant coefficient matrix is obtained, which can save a lot of computational cost. The main novelties of this work are the stability of numerical solutions under the conditions k1Δt≤1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$k_{1}\Delta t\leq1$\end{document} and k2Δt≤1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$k_{2}\Delta t\leq1$\end{document} with two positive constants k1,k2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$k_{1},k_{2}$\end{document} and the optimal error estimates of numerical solutions in different norms. Finally, some numerical results are provided to verify the performances of the Euler implicit/explicit scheme.