Considérant les récents progrès dans le domaine du Calcul Haute Performance, le but ultime des constructeurs aéronautiques tels que Safran Aircraft Engines (SAE) sera de simuler un moteur d'avion complet, à l'échelle 1, utilisant la mécanique des fluides numérique d'ici 2030. Le but de cette thèse de doctorat est donc de donner une contribution scientifique à ce projet. En effet, ce travail est consacré au développement d'une méthode élément finis variationnelle adaptative visant à améliorer la simulation aérothermique du refroidissement des aubes de turbine. Plus précisément, notre objectif est de développer une nouvelle méthode d'adaptation de maillage multi-échelle adaptée à la résolution des transferts thermiques hautement convectifs dans les écoulements turbulents. Pour cela, nous proposons un contrôle hiérarchique des erreurs, basé sur des estimateurs d'erreur sous-échelle de type VMS. La première contribution de ce travail est de proposer une nouvelle méthode d'adaptation de maillage isotrope basée sur ces estimateurs d'erreur sous-échelle. La seconde contribution est de combiner (i) un indicateur d'erreur d'interpolation anisotrope avec (ii) un estimateur d'erreur sous-échelle pour l'adaptation anisotrope de maillage. Les résultats sur des cas analytiques 2D et 3D montrent que la méthode d'adaptation de maillage multi-échelle proposée nous permet d'obtenir des solutions hautement précises utilisant moins d'éléments, en comparaison avec les méthodes d'adaptation de maillage traditionnelles. Enfin, nous proposons dans cette thèse une description des méthodes de calcul parallèle dans Cimlib-CFD. Ensuite, nous présentons les deux systèmes de calcul utilisés pendant le doctorat. L'un d'eux est, en particulier, le super-calculateur GENCI Occigen II qui nous a permit de produire des résultats numériques sur un cas d'aube de turbine complète composé de 39 trous en utilisant des calculs massivement parallèles. / By 2030, considering the progress of HPC, aerospace manufacturers like Safran Aircraft Engines (SAE), hope to be able to simulate a whole aircraft engine, at full scale, using Computational Fluid Dynamic (CFD). The goal of this PhD thesis is to bring a scientific contribution to this research framework. Indeed, the present work is devoted to the development of a variational adaptive finite element method allowing to improve the aerothermal simulations related to the turbine blade cooling. More precisely, our goal is to develop a new multiscale mesh adaptation technique, well suited to the resolution of highly convective heat transfers in turbulent flows. To do so, we propose a hierarchical control of errors based on recently developed subscales VMS error estimators. The first contribution of this work is then to propose a new isotropic mesh adaptation technique based on the previous error estimates. The second contribution is to combine both (i) the coarse scales interpolation error indicator and (ii) the subscales error estimator for anisotropic mesh adaptation. The results on analytic 2D and 3D benchmarks show that the proposed multiscale mesh adaptation technique allows obtaining highly precise solutions with much less elements in comparison with other mesh adaptation techniques. Finally, we propose in this thesis a description of the parallel software capabilities of Cimlib-CFD. Then, we present the two hardware systems used during this PhD thesis. The first one is the lab's cluster allowing the development of numerical methods. The second one however, is the GENCI Occigen II supercomputer which allows producing numerical results using massively parallel computations. In particular, we present a more realistic industrial concerning the cooling of a complete turbine vane composed by 39 holes.