Robust Globally Divergence-Free Weak Galerkin Finite Element Methods for Natural Convection Problems. (2019). Communications in Computational Physics, 26(4), 1039-1070. https://doi.org/10.4208/cicp.OA-2018-0107