Folowing [1]}, the fluid velocity and pressure fields in Equation (REF ) are discretized by
the Brezzi-Douglas-Marini (BDM) element [2]}
and the piecewise constants element, respectively, to ensure local mass conservation.
The temperature field in Equation (REF ) is discretized by the enriched Galerkin (EG) method [3]}, [4]}.
