“…Only in few recent numerical works, the two-dimensional and three-dimensional character of the solar pond have been considered while studying the same problem of heat diffusion [12,13]. In particular, the thermal effect Indices a ambient condition g ground r reference condition (293.15 K) s surface due to the underground was also taken into account [13]. With regard to the pond stability that constitutes, as stated before, a major factor in its operation and performance, it has been studied by several researchers who employed, in most cases, the well-known linear perturbations theory (see in particular [14][15][16][17][18][19]).…”