The problem of icing in aircrafts is known as one of the most critical issues for the aviation community as this matter is directly connected to the flight safety. Due to the high cost of flight tests, many experiments and numerical models have been used to model the icing condition and the influence of the flight. In particular, an accurate estimation of the heat transfer is needed for preliminary design of the anti-icing and de-icing systems. Therefore, the goal of this thesis is to propose a correlation for smooth airfoils before icing (anti-icing systems) and rough airfoils after the ice accretion (de-icing systems) to predict the heat transfer Nusselt number and aerodynamics coefficient. The CFD is used as a high-fidelity method to predict the heat transfer on the smooth and rough airfoils, before and after icing conditions. To have a better understanding of CFD and aircraft icing, the literature review of this thesis is dedicated to a discussion about the history of icing studies and the necessity of continuous icing research, in spite of the numerous achievements. Also, the numerical modelling for icing is discussed in this section. The Spalart-Allmaras turbulence model (SA-SMOOTH) for the smooth airfoils and the Spalart-Allmaras turbulence model with the roughness modification (SA-ROUGH) for the rough airfoils are chosen to build the database. Before building the aerodynamics and heat transfer database, the CFD results are validated and verified against prior experimental and numerical results. Also, to ensure the reliability of the results, the quality of the meshes are evaluated with the GCI method. The database is built for the NACA 0009, NACA0012, and NACA0015 airfoils with Reynolds = 0.625 × 10⁶, 1.25 × 10⁶, 2.5 × 10⁶, and 5 × 10⁶ and 0 ͦ ≤ AOA ≤ 13 ͦ for the smooth and 0 ͦ ≤ AOA ≤ 9 ͦ for the rough airfoils. The Nusselt number and the drag coefficient are then plotted as a function of the lift coefficient to derive the correlation. Three drag coefficient correlations for the C_d and C_dmin for the smooth and C_dminR for the rough airfoils are proposed with three Nusselt number correlations including 𝑁𝑢௔௩௚ for the general correlations and Nu_max and Nu_maxR for the smooth and rough airfoils initial values. To verify the proposed heat transfer correlation, the new proposed correlation is compared with a previous work for the smooth surface.

Le problème du givrage sur les avions est connu comme étant l'un des problèmes les plus critiques pour la communauté aéronautique, car il est directement lié à la sécurité des vols. En raison du coût élevé des essais en vol, de nombreuses expériences et des modèles numériques ont été utilisés pour modéliser les conditions de givrage et l'influence sur les conditions de vol. Ainsi, une estimation précise, une estimation précise du transfert de chaleur est nécessaire pour la conception préliminaire des systèmes d'antigivrage et de dégivrage. Par conséquent, l'objectif de cette thèse est de proposer une corrélation pour les profils lisses avant le givrage (systèmes d'antigivrage) et les profils rugueux après l'accrétion de glace (systèmes de dégivrage) afin de prédire le transfert de chaleur, le nombre de Nusselt et le coefficient aérodynamique. La CFD est utilisée comme une méthode haute-fidélité pour prédire le transfert de chaleur sur les ailes lisses et rugueuses, avant et après les conditions de givrage. Pour avoir une meilleure compréhension de la CFD et du givrage des avions, la revue de la littérature de cette thèse est consacrée à une discussion des études sur le givrage et la nécessité de continuer la recherche sur ce sujet, malgré les nombreuses réalisations. La modélisation numérique du givrage est également abordée dans cette section. Le modèle de turbulence Spalart-Allmaras (SA-SMOOTH) pour les profils lisses et le modèle de turbulence SpalartAllmaras avec modification de la rugosité (SA-ROUGH) pour les profils rugueux sont choisis pour construire la base de données. Avant de construire la base de données sur l'aérodynamique et le transfert de chaleur, les résultats de la CFD sont validés et vérifiés par rapport aux résultats expérimentaux et numériques antérieurs. De plus, pour assurer la fiabilité des résultats, la qualité des maillages est évaluée avec la méthode GCI. La base de données est construite pour les profils NACA 0009, NACA0012, et NACA0015 avec Reynolds = 0.625 × 10⁶, 1.25 × 10⁶, 2.5 × 10⁶, 𝑒𝑡 5 × 10⁶ et 0 ͦ ≤ AOA ≤ 13 ͦ pour les profils lisses et 0 ͦ ≤ AOA ≤ 9 ͦ pour les profils rugueux. Le nombre de Nusselt et le coefficient de traînée sont ensuite tracés en fonction du coefficient de portance pour en déduire la corrélation. Trois corrélations de coefficient de traînée pour C_d et C_dmin pour les profils lisses et C_dminR pour les profils rugueux sont proposées avec les trois corrélations de nombre de Nusselt comprenant 𝑁𝑢௔௩௚ pour les corrélations générales et Nu_max et Nu_maxR pour les valeurs initiales des profils lisses et rugueux. Pour vérifier la corrélation de transfert de chaleur proposée, elle est comparée à un travail précédent pour la surface lisse.