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