In transonic high-pressure turbine stages, oblique shocks originated from vane trailing edges impact the rear suction side of each adjacent vane. High-pressure vanes are usually cooled to tolerate the combustor exit temperature levels, which would reduce dramatically the residual life of a solid vane. Then, it is highly probable that shock impingement will occur in proximity of one of the coolant rows. It has already been observed that the presence of an adverse pressure gradient generates non-negligible effects on heat load due to the increase in boundary layer thickness and turbulence level, with a detrimental impact on the local adiabatic effectiveness values. Furthermore, the generation of a tornado-like vortex has been recently observed that could further decrease the efficacy of the cooling system by moving cold flow far from the vane wall. It must be also underlined that manufacturing deviations and in-service degradation are responsible for the stochastic variation of geometrical parameters. This latter phenomenon greatly alters the unsteady location of the shock impingement and the time-dependent thermal load on the vane. Present work starts from what is shown in literature and provides a highly-detailed description of the aero-thermal field that occurs on a model that represents the flow conditions occurring on the rear suction side of a cooled vane. The numerical model is initially validated against the experimental data obtained by the University of Karlsruhe during TATEF2 EU project, and then an uncertainty quantification methodology based on the probabilistic collocation method and on Padè's polynomials is used to consider the probability distribution of the geometrical parameters. The choice of aleatory unknowns allows to consider the mutual effects between shock-waves, trailing edge thickness and hole diameter. Turbulence is modelled by using the Reynolds Stress Model already implemented in ANSYS ® Fluent ® . Special attention is paid to the description of the flow field in the shock/boundary layer interaction region, where the presence of a secondary effects will completely change the local adiabatic effectiveness values.