An inclined jet in crossflow, characteristic of film cooling applications is studied with highly resolved LES. The configuration reproduces an experimental set-up in which three dimensional measurements of mean velocity and mean flow concentration have been performed using MRV/MRC completed with PIV measurements at some specific locations. The simulations are shown to capture with high fidelity both the separated flow inside the film cooling hole and the distribution of scalar concentration downstream of injection. This paper focuses on the gradient-diffusion hypothesis for the scalar fluxes in RANS modeling. In this configuration we show a large misalignment between turbulent scalar fluxes and the scalar concentration gradient in the streamwise direction.