“…The numerical approach used here follows the recent work of the authors, 62 where steady‐state natural convection around embedded heated bodies in a square enclosure was investigated using the streamfunction–velocity formulation on compact nonuniform Cartesian grids (Figure 2). Equation (1) is discretized using the scheme of Kumar and Kalita, 63 which is given as The expressions for , and are detailed in Kumar and Kalita 63 .…”