An insight into Newton's cooling law using fractional calculus

For small temperature differences between a heated body and its environment, Newton's law of cooling predicts that the instantaneous rate of change of temperature of any heated body with respect to time is proportional to the difference in temperature of the body with the ambient, time being measured in integer units. Our experiments on the cooling of different liquids (water, mustard oil, and mercury) did not fit the theoretical predictions of Newton's law of cooling in this form. The solution was done using both Caputo and Riemann-Liouville type fractional derivatives to check if natural phenomena showed any preference in mathematics. In both cases, we find that cooling of liquids has an identical value of the fractional derivative of time that increases with the viscosity of the liquid. On the other hand, the cooling studies on metal alloys could be fitted exactly by integer order time derivative equations. The proportionality constant between heat flux and temperature difference was examined with respect to variations in the depth of liquid and exposed surface area. A critical combination of these two parameters signals a change in the mode of heat transfer within liquids. The equivalence between the proportionality constants for the Caputo and Riemann-Liouville type derivatives is established.