Laminar free convection flow in warm bathing water had just been investigated numerically for a range of Reynolds number 0 ≤ Re ≤ 110 keeping Fr = 2.5 and Pr = 9.5 fixed throughout the study. The result showed that cabbeling changes slightly with Re owing to the increasing Re, even though, the velocity in such free convection flows are usually small. The result showed that volume of hot water at the upper section at some point in time depleted completely and in turn induces the entire ambient water temperature which later became the same temperature everywhere without any external influence. These results require just a little mixing for temperature between 0◦C and 10◦C to attain Tm. The time taken for dense fluid to sink to the bed for smaller Re is slightly longer as compared to the time it takes as Re increases. Temperature profiles were also analysed at some points (X, 69; 30; 10) below the contact layer and plotted against the x-coordinate. In a similar manner, profiles of x-component and y-component velocities were also determined at (X, 69) for the various Re cases as considered and plotted against the x-coordinate. Fluctuations in the temperature profiles describes the convection precess as both hot and cold fluid mixes in all direction and continue to deplete further. Meanwhile, fluctuations in the curves for the y-component velocity profile indicates that even as the dense fluid continue to descend, fluid that is still positively buoyant moves upwards. The maximum time taken for descending dense fluid to reach domain floor and the time taken to attain that depth were also considered. From the empirically determined data set, we could identify a single regime of Re-dependence and shown by the straight line in Figs. 7, which represent best fit power law obtained by linear regression of logRe on logτn (see equation ).