A quantified parametric study for the heat transfer acting on a hypersonic blunt body with counterflowing jets is presented. Three-dimensional turbulent Navier-Stokes equations are solved to simulate freestream-jet interactive flowfields. The freestream and jet controlling parameters are treated as input sources of variation, and a point-collocation non-intrusive polynomial chaos (NIPC) method is utilized to quantify the variations in the output surface heat flux and total surface heat load acting on the blunt body by identifying the maximum and minimum of surrogate response values predicted by the NIPC. All of the sample cases are confirmed to form steady jet structures. Furthermore, through a global sensitivity analysis, Sobol indices evaluated by the NIPC, are used to rank the contributions of each input parameter to the variation in output quantities of interest. It is found that the designed upstream injection of baseline case effectively reduces the heat transfer to body surface compared with the no-jet case. The variations of input parameters induce remarkable variations of output heat flux and total heat load. The sensitivity analysis indicates that the jet-to-freestream total-pressure ratio is the top contributor to variations in heat flux, followed by the freestream Mach number. The jet total temperature is mainly important on the front part of forebody, while the contributions of jet Mach number and freestream temperature slightly increase downstream. The freestream density has the smallest effect. The sensitivity of total heat load to input parameters coincides with that of heat flux. This parametric study is expected to illustrate the significance of flow-controlling parameters to the heat transfer over blunt body and provide insight for aerothermal management by using counterflowing jets in hypersonic flows. (C) 2017 Elsevier Ltd. All rights reserved.