A parametric reduced order model based on proper orthogonal decomposition with Galerkin projection has been developed and applied for the modeling of heat transport in T-junction pipes which are widely found in nuclear power reactor cooling systems. Thermal mixing of different temperature coolants in T-junction pipes leads to temperature fluctuations and this could potentially cause thermal fatigue in the pipe walls. The novelty of this paper is the development of a parametric ROM considering the three dimensional, incompressible, unsteady Navier-Stokes equations coupled with the heat transport equation in a finite volume regime. Two different parametric cases are presented in this paper: parametrization of the inlet temperatures and parame...