La capacité de simuler la dynamique de la couche limite turbulente représente aujourd'hui un enjeu important pour la prévision de l'aérodynamique instationnaire et de l'aéroacoustique des aéronefs et des véhicules terrestres. Aussi, les travaux présentés dans ce manuscrit proposent une méthode originale de simulation de la dynamique des écoulements turbulents pariétaux, à haut nombre de Reynolds, sur des géométries curvilignes. L'approche ZDES, dans son mode " Wall Modelled Large Eddy Simulation ", est ici retenue. Dans un premier temps, une méthode de réactivation turbulente, par combinaison de l'injection de turbulence synthétique et de l'application de termes de forçage, est développée afin de permettre une résolution WMLES locale des régions d'intérêt au sein d'une simulation majoritairement RANS. Puis l'étude des interactions entre la physique et la méthode de résolution numérique (ZDES), sur le cas d'une couche limite turbulente sans gradient de pression en développement spatial jusqu'à Re13 000, a conduit à une généralisation des conditions d'emploi du mode WMLES de la ZDES. Des post-traitements instationnaires avancés ont permis de démontrer la capacité de la méthode à simuler la dynamique particulière de la zone externe à ces grands nombres de Reynolds. Enfin, les différentes méthodes développées au cours de cette étude ont été appliquées à la simulation d'une manche à air coudée de section rectangulaire. Cette nouvelle méthode a permis de simuler avec succès la dynamique multi-échelles de cet écoulement et des distorsions dynamiques associées, pour un coût environ 50 fois inférieur à celui d'une approche LES classique

Mieux comprendre les spécificités de la dynamique des couches limites à grand nombre de Reynolds malgré les contraintes métrologiques et son coût de simulation numérique est crucial. A titre d'exemple, cette dynamique peut déterminer plus de la moitié de la traînée d'un avion en croisière. Décrire la turbulence pariétale peut guider la résolution numérique d'une partie des fluctuations à un coût maîtrisé par des stratégies WMLES (simulation des grandes échelles avec modèle de paroi). Les présentes analyses physiques de couches limites turbulentes incompressibles à gradient de pression nul et à grand nombre de Reynolds s'appuient sur des simulations numériques avancées. Après validation d'une base de données, le frottement moyen pariétal est décomposé selon l'identité FIK (Fukagata et al. (2002)), dont l'application malgré le développement spatial est discutée. Une analyse spectrale montre que les grandes échelles (\lambda_x > \delta) contribuent à environ la moitié du frottement vers Re_\theta = 10^4. Les limitations de l'identité FIK motivent la dérivation d'une décomposition physique de la génération du frottement dont le comportement asymptotique est alors relié à la production d'énergie cinétique turbulente dans la zone logarithmique. Pour mieux reconstruire les spectres spatiaux, une nouvelle méthode d'estimation de la vitesse de convection turbulente en fonction de la longueur d'onde des fluctuations, adaptée au développement spatial et à des signaux temporels de durée finie, est dérivée, interprétée et évaluée à Re_\theta = 13000. Certaines des conclusions éclairent des modifications d'une stratégie WMLES, le mode III de la méthode ZDES. / Better understanding the specificities of the dynamics of high-Reynolds number boundary layers despite metrological constraints and its numerical simulation cost is crucial. For instance, this dynamics can determine more than half of the drag of a cruising aircraft. Describing wall turbulence can guide the numerical resolution of some of the fluctuations at a limited cost by WMLES strategies (wall-modelled large eddy simulation). The present physical analyses of zero-pressure gradient incompressible turbulent boundary layers at high Reynolds number rely on advanced numerical simulations. After validating a database, mean skin friction is decomposed by means of the FIK identity (Fukagata et al. (2002)), whose application despite the spatial growth is discussed. A spectral analysis shows that the large scales (\lambda_x > \delta) contribute approximately half of the friction near Re_\theta = 10^4. The limitations of the FIK identity motivate the derivation of a physical decomposition of the generation of friction whose asymptotic behaviour is then related to turbulent kinetic energy production in the logarithmic layer. In order to better reconstruct spatial spectra, a new method to estimate the turbulent convection velocity as a function of the wavelength of the fluctuations, adapted to spatial growth and to temporal signals of finite duration, is derived, interpreted, and assessed at Re_\theta = 13000. Some of the conclusions enlighten modifications to a WMLES strategy, mode III of the ZDES method.