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