Molten salts have attracted a spate of quantity of interests in energy and chemistry fields due to their large volumetric capacity, and thermodynamic stability at elevated temperature. However, there are significant uncertainties in the physical properties of molten salts, such as density, specific heat capacity, dynamic viscosity, and thermal conductivity, due to inexperienced measurement techniques and technology. These uncertainties could have certain impact on the evaluation of molten salts heat transfer characteristics. Therefore, an uncertainty analysis is performed regarding FLiNaK thermophysical properties on convective heat transfer characteristics in circular tube. Then, a polynomial chaos expansion (PCE) method is performed to study the uncertainty affecting the heat transfer characteristic. In the process of uncertainty analysis, tensor product quadrature nodes are used to calculate representative sample points and reduce computational costs. Three primary physical properties of molten salt are taken into account as the input parameters. The Sobol composition method is also used to analyze the contributions of each parameter to heat transfer. The results of the uncertainty analysis suggest that the uncertainties in dynamic viscosity, thermal conductivity, and specific heat capacity have a significant impact on the heat transfer of molten salt, contributing to 80 %, 14 %, and 6 % of the total variability, respectively. It also suggests that the polynomial chaos expansion methodology is both novel and reliable when applied to uncertainty analysis of molten salt.