In the present work, analytical solutions are presented for thermal convection of the linear Phan-Thien-Tanner fluid (LPTT) in slits and tubes of constant wall temperature by taking account of the viscous dissipation term. Unlike the similar previous studies in which the advection term was neglected in the heat transfer equation, it is considered in this investigation. A continuous relation between the Nusselt number and the Brinkman number is obtained. Expressions for the temperature distribution are derived in closed form and in terms of a Frobenius series for the slit and tube flows, respectively. Based on these solutions, the effects of fluid elasticity and Brinkman number on thermal convection of LPTT fluid flows are studied in detail. It is shown that at negative Brinkman numbers (fluid cooling), increasing the Deborah number leads to a decrease in the Nusselt number, but an increase in the centerline temperature. Nonetheless, this trend is opposite for positive Brinkman numbers (fluid heating), i.e., an increase in the Nusselt number and a decrease in the centerline temperature. Also, there is a Brinkman number beyond which the Nusselt number is smaller than zero, meaning that there is weak heat convection in the flow. Also, the results confirm that the extensibility parameter affects the temperature profile in the same way as the Deborah number. Keywords Viscoelastic fluid • Phan-Thien-Tanner model • Forced convection • Constant wall temperature tube and slit • Viscous dissipation List of symbols A Area of cross section, m 2 Br Brinkman number, Br = U 2 ∕k(T w −T m) c p Specific heat, J kg −1 K −1 d h Hydraulic diameter, d h = 2R for tube flow and d h = 4H for slit flow, m Deformation rate tensor, Eq. (8) De Deborah number, defined as De = U∕R for tube case and De = U∕H for slit case F Dimensionless function h Heat transfer coefficient, W m −2 K −1 H Half of the distance between two parallel plates, m J Equals 0 for slit and 1 for tube k Conductivity coefficient, W m −1 K −1 K Equals 1.5 and 2 for slit and tube, respectively Nu Nusselt number, Nu = hd h ∕k p Pressure, Pa P Perimeter of cross section, m R Tube radius, m T Fluid temperature, K u Axial velocity, ms −1 U Mean velocity, ms −1 Velocity vector y Radial (tube) or transverse (slit) direction, m z Axial direction, m Greek symbols Equals 1 for pure entropy elasticity and 0 for pure energy elasticity 1 1 =1.5U N ∕U for slit and 1 =2 U N ∕U for tube Constant, Eq. (12) Extensibility coefficient Constant viscosity coefficient, Pa s Relaxation time, s Density, Kg m −3 ∇ Upper-convected derivative of stress tensor, Eq. (9) Viscous dissipation, Eq. (15)