Influences of heat release by the hydrogen combustion in supersonic turbulent boundary layers are numerically studied using Reynolds-averaged Navier-Stokes equations. The adopted Reynolds-averaged Navier-Stokes methodology is first validated by comparing the numerical results with the existing experimental data. Studies on the effects of the flame perpendicular position inside the boundary layer reveal that, while the flame is restricted around the edge of the boundary layer, the heat release may slightly reduce rather than increase the wall heat flux because of the suppression effect on the turbulent energy transport due to heat release. However, as the flame moves toward the wall, the skin-friction reduction effect would not be obviously strengthened, but the wall heat flux could be dramatically enhanced by the increase of near-wall chemical reactions. At a given hydrogen mass flow rate, the injection scheme with a higher injection height and a lower injection velocity could be helpful to achieve a larger skinfriction reduction while maintaining a lower wall heat flux. Finally, analysis of the heat-release effects on the velocity law of the wall shows that van Driest's velocity law largely deviates from the computational results, whereas White's velocity law remains close to the numerical results within a region of approximately y < 300. Nomenclature B = law-of-the-wall constant (5.0) c f = local skin friction coefficient c p = specific heat at constant pressure, J∕kg · K D = mass diffusivity, m 2 ∕s E = total energy, J∕kg H = total enthalpy, J∕kg h s = static enthalpy for species s, J∕kg h = hydrogen injection height at the entrance, m k = turbulent kinetic energy, m 2 ∕s 2 M = Mach number Pr = Prandtl number Pr t = turbulent Prandtl number p = pressure, MPa q w = local heat flux, MW∕m 2 r= Pr 1∕3 , recovery factor Sc t = turbulent Schmidt number T = temperature, K u = velocity, m∕s u τ = friction velocity, m∕s Y hw = hydrogen mass fraction at the wall Y s = mass fraction of species s δ = boundary-layer thickness, m κ = von Kármán constant (0.41) λ = molecular thermal conductivity, W∕m · K μ = molecular viscosity, N · s∕m 2 μ t = turbulent eddy viscosity, N · s∕m 2 ρ = density, kg∕m 3 τ = shear stress, N∕m 2 ω = specific turbulent dissipation, 1∕s _ ω = mass production per unit volume per unit time, kg∕m 3 · s Subscripts e = edge of boundary layer j = jet s = species s w = wall ∞ = freestream Superscripts t = turbulent = Reynolds averagẽ = Favrè average