The non-orthogonal boundary-fitted coordinate transformation method is applied to the solution of steady three-dimensional momentum and energy equations in laminar flow to obtain temperature field and Nusselt numbers in the thermal entry region of straight ducts of different cross sectional geometries. The conservation equations originally written in Cartesian coordinates are parabolized in the...