Similarity solution and Runge Kutta method to a thermal boundary layer model at the entrance region of a circular tube

Purpose The purpose of this paper is to re-examine the assumptions implicit in Leveque’s approximation, and the variation of the temperature and the thickness of the boundary layer were illustrated using the developed solution. The analytical solutions are then checked against numerical solution programming by FORTRAN code obtained via using Runge–Kutta fourth-order (RK4) method. Finally, other important thermal results obtained from this analysis, such as approximate Nusselt number in the thermal entrance region, was discussed in detail. After that, the analytical results of the present paper are validated with certain previous investigations which were found in the specialized literature. Design/methodology/approach By defining a similarity variable, the governing equations are reduced to a dimensionless equation with an analytic solution in the entrance region. This paper gives justification for the similarity variable via scaling analysis, details the process of converting to a similarity form and presents a similarity solution. The calculation methodology for numerical resolution is based on the RK4 technique. Findings The profiles of the solutions are provided from which the authors infer that the numerical and exact solutions agreed very well. Another result that the authors obtained from this paper is the number of Nusselt in the thermal entrance region for which a parametric study was carried out and discussed well for the impact of scientific contribution. Originality/value The novelty of this paper is the application of the RK4 with a step size control, as a sequential numerical method of a ODEs system compared with the exact similarity solution of the thermal boundary layer problem.