Un micro-drone à voilure tournante est l’appareil aérien optimal pour assister un rover d’exploration à lanavigation sur la planète Mars. Toutefois, les écoulements qu’il rencontre sont compressibles à très faiblenombre de Reynolds, ce qui constitue un domaine de l’aérodynamique inédit et quasiment inexploré à cejour. L’objectif de la thèse est de comprendre, simuler et recréer expérimentalement les phénomènes aérodynamiquesliés au régime inédit des écoulements martiens pour concevoir un système propulsif performant.Après avoir validé les outils de simulation numérique, le comportement instationnaire des écoulements estétudié sur des géométries 2D et 3D. L’écoulement est dominé par la viscosité : les couches limites laminairessont épaisses et le décollement a beaucoup d’influence sur son comportement très instationnaire.Par la suite, plusieurs millions de géométries de profil sont évaluées par un processus d’optimisation basésur un code 2D stationnaire. Les profils optimisés sont fortement cambrés (entre 5.5% et 7%) et de faibleépaisseur relative (e/c ∼ 2%). Le bord d’attaque et le bord de fuite sont très cambrés pour permettrerespectivement l’adaptation à l’écoulement incident et la fixation du point de décollement de la couchelimite. À partir des géométries de profils, l’ensemble du système propulsif est optimisé par intégration despolaires 2D. La théorie des éléments de pale permet de déterminer rapidement les configurations les plusperformantes aérodynamiquement. Et une méthode de sillage libre permet l’optimisation de rotors isoléset de systèmes propulsifs coaxiaux. Les rotors ont des solidités et des vrillages importants, ce qui rappelleles formes d’hélices marines. Les simulations Navier-Stokes 3D mettent en évidence la tridimensionnalitédes écoulements sur la pale, elle est fortement corrélée avec la solidité du rotor et le vrillage de bout depale. La rotation stabilise la couche limite et donne lieu à un décollement stable au bord d’attaque pourcertaines géométries. Le dévrillage en bout de pale permet de stabiliser le tourbillon et de diminuer la perteinduite. Un banc de mesure est placé dans un caisson dépressurisé pour estimer les efforts de poussée et decouple générés par les rotors optimisés en conditions aérodynamiques martiennes. Les essais permettentde valider les tendances d’estimation des codes de simulation ainsi que les processus d’optimisation. Laconfiguration bi-rotors coaxiaux, en comparaison avec une configuration à deux rotors adjacents, permetun gain d’encombrement de moitié pour une perte sur la puissance de seulement 15%. C’est la configurationla plus adaptée pour un micro-drone en atmosphère ténue. Un système propulsif coaxial optimisé (dediamètre 30 cm) permettrait de sustenter un micro-drone d’environ 400 grammes en conditions nominalessur la planète Mars. / A micro-rotorcraft is the most suited aerial vehicle for rover navigation assistance on Mars. The martianatmosphere’s density, being hundred times lower than on Earth, requires the micro-drone to hover at highrotational speed. Hence, flows on the blade are both compressible and at very low Reynolds number (fewthousands). It constitutes a new aerodynamic domain to be explored. The purpose of the dissertation isto understand, simulate and experimentally duplicate the aerodynamic phenomena in a view to design anefficient propulsion system. After a phase of validation of the simulation tools, the flows’ unsteady behavioris studied on 2D and 3D geometries. Wall flow is highly viscous : laminar boundary layers are thick andtheir separation has a huge influence on its unsteadiness. Then, several millions of airfoil geometries areevaluated by an optimization process based on a steady 2D solver. As final result, the optimized airfoildisplays a highly cambered shape (between 5.5% and 7%) with low relative thickness (t/c ∼ 2%). Leadingand trailing edges are strongly cambered, allowing proper incoming flow adaptation and late boundarylayer separation. Based on this airfoil geometry, rotor shapes are optimized by two methods. Blade elementtheory provides a quick investigation of the most aerodynamically efficient configurations. And a free wakesolver is applied for the final design of isolated rotors and coaxial dual-rotors. Optimized shapes exhibitimportant twist and solidity, evoking marine propellers. 3D Navier-Stokes simulations highlight the flow’sthree-dimensional mechanisms on the blade, which are highly correlated to the rotor’s solidity and twist.Rotation stabilizes the blade’s boundary layer and a stable leading edge separation is observed in somecases. Blade tip twist reduction diminishes the tip vortex and the induced loss. A thrust and torquemeasurement setup is placed in a depressurized tank for the evaluation of optimized rotors in martianatmospheric conditions. Performed experiments confirm the solvers’ trend and validate the design process.Therefore, chosen coaxial dual-rotors configuration provides an important size gain for a low correspondingpower loss (∼ 15%) compared to adjacent dual-rotors. Coaxial dual-rotors are the optimal configuration fora micro-rotorcraft in low-pressure atmosphere flight conditions. Such propulsion system (with a diameterof 30 cm) could lift a 400 grams micro-rotorcraft in hover on the planet Mars.

Dans le cadre de la problématique du refroidissement des déchets radioactifs, ce travail a pour objectifs l'étude de l'écoulement autour d'un cylindre chauffant monté en paroi soumis à un écoulement transverse selon l'approche URANS. Les limitations bien connues de la modélisation au premier ordre conduisent à considérer une approche au deuxième ordre. C'est dans ce sens qu'un modèle d'échange thermique est développé puis validé sur des cas académiques.<br>Dans un premier temps, lorsque le régime de convection mixte est dominant, ces calculs, réalisés par le modèle bas-Reynolds k-w SST, et complétés par un calcul isotherme, mettent en évidence de nombreuses structures tourbillonnaires recensées par la bibliographie. Un calcul avec le modèle haut-Reynolds Rij-epsilon SSG est aussi effectué. Avec le modèle k-w SST, l'échange thermique est convenablement reproduit par rapport à l'expérience VALIDA, conduite par le CEA, alors qu'avec le modèle Rij-epsilon SSG celui-ci est fortement sous-estimé, certainement en raison de l'utilisation de lois de paroi. <br>Dans un second temps, lorsque la convection naturelle est prépondérante, la topologie des écoulements est alors complètement différente et le transfert thermique est bien moins fidèle à l'expérience VALIDA.<br>Le modèle à pondération elliptique EBRSM, issu de l'approche de Durbin, consiste en la prise en compte des effets de paroi, et en particulier le blocage à la paroi. D'après ce formalisme, un modèle algébrique des flux thermiques avec pondération elliptique est développé, l'EBAFM. Avec ce modèle, des tests a priori dans les trois régimes de convection puis les simulations sur les mêmes cas montrent une nette amélioration dans la prédiction des différents écoulements faisant place à une perspective intéressante d'un modèle intermédiaire entre une formulation SGDH et des équations de transport.

Dans le cadre de la problématique du refroidissement des déchets radioactifs, ce travail a pour objectifs l'étude de l'écoulement autour d'un cylindre chauffant monté en paroi soumis à un écoulement transverse selon l'approche URANS. Les limitations bien connues de la modélisation au premier ordre conduisent à considérer une approche au deuxième ordre. C'est dans ce sens qu'un modèle d'échange thermique est développé puis validé sur des cas académiques. <br /><br />Dans un premier temps, lorsque le régime de convection mixte est dominant, ces calculs, réalisés par le modèle bas-Reynolds ßt, et complétés par un calcul isotherme, mettent en évidence de nombreuses structures tourbillonnaires recensées par la bibliographie. Un calcul avec le modèle haut-Reynolds Rij-epsilon SSG est aussi effectué. Avec le modèle k-w SST, l'échange thermique est convenablement reproduit par rapport à l'expérience VALIDA, conduite par le CEA, alors qu'avec le modèle Rij-epsilon SSG celui-ci est fortement sous-estimé, certainement en raison de l'utilisation de lois de paroi. <br /><br />Dans un second temps, lorsque la convection naturelle est prépondérante, la topologie des écoulements est alors complètement différente et le transfert thermique est bien moins fidèle à l'expérience VALIDA.<br /><br />Le modèle à pondération elliptique EB-RSM, issu de l'approche de Durbin, consiste en la prise en compte des effets de paroi, et en particulier le blocage à la paroi. D'après ce formalisme, un modèle algébrique des flux thermiques avec pondération elliptique est développé, l'EB-AFM. Avec ce modèle, des tests a priori dans les trois régimes de convection puis les simulations sur les mêmes cas montrent une nette amélioration dans la prédiction des différents écoulements faisant place à une perspective intéressante d'un modèle intermédiaire entre une formulation SGDH et des équations de transport.

L'objectif de ce travail est de prendre en compte les instationnarités naturelles à grande échelle dans les écoulements décollés et à un coût plus faible que la LES, tout en s'intéressant à la modélisation des effets de paroi par des modèles statistiques au second ordre. S'inspirant des approches de Durbin, le modèle à pondération elliptique EB-RSM reproduit l'effet non-local de blocage, en résolvant une équation différentielle sur le terme de pression. La limite à deux composantes de la turbulence est bien prédite en canal. Ce modèle est appliqué à la marche descendante, dans une approche URANS. Nous avons montré que les erreurs numériques peuvent être suffisantes pour exciter le mode le plus instable de la couche cisaillée, et aboutir à une solution instationnaire. La solution est stationnaire quand on raffine le maillage, rendant l'URANS peu fiable. Récemment, Schiestel \& Dejoan ont proposé le modèle hybride non-zonal PITM. Le coefficient $C_{\e_2}$ de l'équation de la dissipation devient fonction de la coupure dans le spectre, et la valeur $C_{\e_1}=3/2$ est déduite par ces auteurs. Nous avons donné une formulation plus générale où la valeur de $C_{\e_1}$ est quelconque. Pour offrir un formalisme plus cohérent aux modèles hybrides non-zonaux dans les écoulements de paroi, une approche basée sur un filtrage temporel est proposée. Enfin, l'adaptation du modèle EB-RSM dans un cadre hybride a été réalisée. Les résultats en canal sont encourageants : la transition continue d'un modèle RANS en proche paroi à une LES au centre du canal est mise en évidence. Le transfert d'énergie des échelles modélisées vers celles résolues est bien reproduit quand on raffine le maillage.