“…For natural convection problem, which has a wide range of applications in many research fields (Wu et al , 2017b, 2016), the density difference in the fluid occurring due to temperature gradient is the driving mechanism of fluid motion[1]. In case that the density variation is small, it can be modeled by using a Boussinesq approximation, which treats the density as a constant but with an added buoyancy force, and most literature studies the constant density natural convection based on the Boussinesq approximation (Boland and Layton, 1990; Du et al , 2015; Feng et al , 2011; Huang et al , 2015, 2013, 2012; Liao, 2012, 2010; Si et al , 2014; Szumbarski et al , 2014; Su et al , 2017a, 2017b, 2014a, 2014b; Sun et al , 2011; Davis, 1983; Wang et al , 2018a, 2018b; Wu et al , 2015b, 2017a, 2016; Zhang et al , 2016, 2018]). However, in most geophysical flows and many other situations, fluid motion is usually driven by large temperature differences, which results in a considerable density change and the Boussinesq approximation is no longer valid.…”