Aerodynamic Simulations of Helicopter Rotors through the Actuator Line Method

Les hélicoptères se distinguent des autres aéronefs par leur polyvalence inégalée dans leur capacité de vol vertical efficace leur permettant d'accéder facilement à des zones reculées. Bien que les capacités de vol vertical et leurs performances soient d'une importance primordiale, l'hélicoptère d'aujourd'hui doit encore fonctionner efficacement dans tous les régimes de vol, comme le vol vers l'avant, la descente et en interaction étroite avec les obstacles environnants sur des terrains difficiles. Ces défis, en particulier le dernier, s'étendent également à un champ plus large de véhicules à décollage/atterrissage vertical plus généraux tels que les drones sous la forme de micro-véhicules aériens et l'avènement de la mobilité aérienne urbaine. Toutefois, en raison de la grande polyvalence et de l'applicabilité de l'aéronef, il serait impossible pour les concepteurs d'envisager tous les cas de vol et les configurations hors design dans lesquels les opérateurs utilisent l'aéronef. Pour certaines de ces applications spécifiques, il est nécessaire d'évaluer les performances de vol des hélicoptères dans des conditions de vol difficiles. L'objectif principal de cette thèse est d'étudier les calculs aérodynamiques haute-fidélité des rotors d'hélicoptères dans diverses conditions de vol. Des conditions telles que le vol stationnaire, le vol avant, l'effet de sol et l'interaction avec un bâtiment de faible hauteur sont évaluées. Les simulations traditionnelles de haute fidélité basées sur les méthodes Navier-Stokes par Moyennage de Reynolds Instationnaire (U-RANS) peuvent prendre jusqu'à plusieurs semaines sur des centaines de coeurs de CPU pour compléter un calcul pour une configuration, ce qui les rend intensives en temps et en calcul. Pour tenter d'enrayer ce problème, un modèle, toujours basé sur un cadre U-RANS haute fidélité, est exploré. La méthode développée doit simuler et prédire avec précision les principaux paramètres de performance du rotor tels que la poussée, le couple, le facteur de mérite et le sillage général du rotor. Le modèle retenu est la méthode de la ligne actuatrice (ALM). Cette méthode, qui est largement utilisée pour les simulations d'éoliennes, remplace les pales du rotor par des termes sources de quantité de mouvement dans les équations U-RANS. La suppression du maillage des pales réduit considérablement la taille du maillage de calcul, ce qui diminue le coût de calcul. Tout d'abord, afin d'adapter l'ALM aux écoulements de rotor d'hélicoptère, une étude paramétrique est réalisée dans un cadre simplifié. Dans un cadre bidimensionnel, différentes méthodologies d'échantillonnage de la vitesse sont testées avec l'échantillonnage intégral de la vitesse s'avérant supérieur. D'autres sensibilités du modèle liées à la modélisation et à l'écoulement sont évaluées et caractérisées. Ensuite, une simple extension en 3D pour une application d'aile fixe montre l'adéquation de la méthode dans le traitement des écoulements 3D. Ensuite, l'ALM est appliquée à un rotor en vol stationnaire. Un kernel Gaussien tronquénormalisé correctement ajusté donne des résultats supérieurs en termes de capacités de prédiction des coefficients du rotor. L'ALM produit également un sillage tourbillonnaire similaire à celui d'une simulation avec les pales pleinement résolues. Troisièmement, la méthode est appliquée à un cas de vol axial où les capacités de prédiction du modèle ALM sont presque identiques par rapport à la simulation de référence avec pales entièrement résolue. L'ALM est également capable de prédire la diminution linéaire du facteur de mérite pour un rotor en montée avec une grande précision. En descente, l'accord est moins bon en raison de l'absence de régime d'anneau tourbillonnaire sensé être présent tel qu'observé dans les résultats expérimentaux. Enfin, l'ALM est testée dans des conditions d'écoulement de plus en plus complexes, y compris le vol stationnaire en effet de sol, le vol avant en effet de sol, en effet de sol en interaction avec un obstacle en forme de boîte, avec et sans vent incident. Dans tous les cas, l'ALM est capable de prédire les tendances globales de l'écoulement avec une petite erreur de prédiction de la poussée dans les états d'effet de sol les plus extrêmes. La méthode est également capable de prédire la magnitude relative et le changement d'orientation du coefficient de moment. L'ALM est par ailleurs en bon accord avec les simulations de référence avec les pales pleinement résolues. La concordance globale de l'ALM avec les simulations à pales entièrement résolues et les données expérimentales démontre ses capacités à prédire avec précision l'écoulement autour d'un rotor ainsi que ses performances dans des conditions difficiles. La méthode est donc jugée apte à remplacer les simulations à pales entièrement résolues qui sont plus coûteuses en temps de calcul.