Les traînées de condensation, ces panaches blancs souvent observés dans le sillage des avions, contribuent à la perturbation du bilan radiatif terrestre. Composées de cristaux de glace, elles peuvent en effet persister et entraîner une augmentation de la nébulosité. Le but de cette thèse est de mieux comprendre et décrire le processus de formation des cristaux de glace dans le champ proche d’un avion. L’écoulement de sillage est constitué d’une phase gazeuse (air et vapeur) et d’une phase solide (particules) en interaction. La simulation numérique de cet écoulement repose sur une méthode mixte eulérienne (gaz) et lagrangienne (particules). La simulation des grandes échelles est utilisée pour résoudre les équations de Navier-Stokes relatives au fluide turbulent, instationnaire et compressible. La croissance de chaque particule par condensation de la vapeur d’eau est calculée grâce au modèle microphysique de Fukuta-Walter. Deux régimes sont considérés : le régime jet et le régime d’interaction entre le jet et le tourbillon de sillage. Dans une première partie sont présentées la modélisation et les méthodes numériques employées. La deuxième partie consiste à relater les résultats obtenus lors des simulations. Le régime jet est d’abord traité afin d’étudier le rôle des paramètres caractérisant l’écoulement sur la taille et la densité numérique des cristaux de glace. Un champ tourbillonnaire est ensuite superposé à l’écoulement de jet pour examiner l’influence du tourbillon de sillage sur les processus microphysiques. Les simulations ont démontré l’importance du choix des paramètres décrivant l’écoulement pour caractériser avec précision les traînées de condensation.

Le contexte général de cette étude concerne l’impact des traînées de condensation, ces fameux panaches blancs fréquemment observés dans le sillage des avions, sur l’atmosphère. D’un point de vue aérodynamique, la formation des traînées de condensation se caractérise par l’interaction entre un jet turbulent et un tourbillon de sillage. L’objectif de cette thèse est de contribuer à une meilleure compréhension des phénomènes thermiques et dynamiques mis en jeu dans cet écoulement. Ce travail repose sur la résolution numérique des équations de Navier-Stokes, et de conservation de l’énergie, formulées pour le cas d’un écoulement tridimensionnel, instationnaire et compressible. Deux approches sont considérées: la simulation numérique directe et la simulation des grandes échelles. Une simulation temporelle de la transition à la turbulence d’un jet rond non-isotherme est effectuée sans tenir compte du champ tourbillonnaire. A l’issue de cette simulation, un modèle de tourbillon de sillage est superposé à l’écoulement de jet. La première partie de cette thèse présente les deux approches, les différents modèles de sous-maille choisis pour les simulations des grandes échelles, ainsi que les méthodes numériques employées. La deuxième partie est consacrée à la simulation de l’écoulement de jet, et ici, l’objectif est de déterminer le modèle de sous-maille approprié à cette configuration d’écoulement. La troisième partie est dédiée à la simulation de l’interaction entre le jet et le tourbillon de sillage. Les résultats sont comparés à ceux issue d’une campagne d’essais. Les simulations ont démontré le développement de grosses structures de la turbulence autour du cœur tourbillonnaire, dans lesquelles se concentre le champ de température

La prise en compte des traînées de condensations (contrails) dans les modèles de grande échelle requiert une paramétrisation. Cependant, les données actuelles sont insuffisantes pour cette paramétrisation. La simulation numérique détaillée d'un contrail complète les données et elle contribue à mieux définir ces paramétrisations. Les études des aéronefs ont mis en évidence le rôle de la turbulence atmosphérique dans la dynamique des tourbillons de sillage. Pourtant, la majorité des simulations numériques de contrails réalisées à ce jour utilisent une représentation simplifiée et paramétrée de cette turbulence, plutôt que de la résoudre explicitement. Le travail réalisé ici consiste à mettre en place une simulation de contrail avec une résolution tridimensionnelle de la turbulence atmosphérique. Dans un premier temps, une méthode de forçage stochastique a été mise en place pour engendrer des écoulements turbulents dans un milieu stratifié, et on montre la capacité de cette méthode à reproduire la turbulence de l'atmosphère. Dans un second temps, une simulation de contrail a été mise en place en utilisant la turbulence atmosphérique engendrée précédemment. Les effets du niveaux de turbulence, de la température et de la saturation en vapeur d'eau sur le contrail sont étudiés pour les quatre premières minutes : la turbulence contrôle la destruction des tourbillons de sillage, alors que la saturation et la température contrôlent la persistance et le taux de formation de glace. Les propriétés microphysiques et optiques obtenues sont en accord avec les observations et des lois sont proposées pour adapter ces propriétés aux paramétrisations des modèles de grande échelle.

Cette étude a pour but de prédire les trajectoire de jets issus des inverseurs de poussée. Ces dispositifs équipent les avions et génèrent une contre-poussée qui permet de freiner l'appareil lorsque celui-ci est au sol en phase d'atterrissage. L'écoulement sort des inverseurs de poussée sous forme de jet. Un tel écoulement est complexe et s'étend sur une grande distance. Nous avons alors opté pour la méthode particulaire de type « Vortex Blob ». Le domaine d'étude est infini, la résolution est instationnaire et la discrétisation limitée aux zones rotationnelles. Celles-ci se retrouvent sous forme de particules porteuses d'une information tourbillonnaire que l'on suit dans leur mouvement de manière lagrangienne. Dans un premier temps, nous passons en revue tous les types d'écoulement de jet et nous ressortons les principales caractéristiques, notamment pour l'écoulement de jet à buse débouchant dans un écoulement transverse où une paire de tourbillons contra-rotatifs apparaît et persiste loin en aval. L'activité tourbillonnaire est intense et l'évolution du jet varie selon le rapport d'injection du jet (vitesse du jet rapporté à la vitesse de l'écoulement transverse). Nous abordons ensuite l'aspect numérique de l'étude présentant la méthode particulaire. Nous utilisons une méthode particulaire tridimensionnelle à laquelle nous avons adapté un algorithme rapide de type « Tree code », un remmaillage global de l'écoulement ainsi qu'un schéma de diffusion de type « Particle Strenght Exchange » qui peut prendre en compte un modèle de Simulation des Grandes Échelles (SGE). Le code est de plus entièrement parallélisé grâce aux librairies MPI. Puis, nous traitons de manière académique les simulations de jets transverses. Nous observons et étudions l'évolution des structures tourbillonnaires existantes dans de tels écoulements en fonction du rapport d'injection du jet, d'une variation temporelle de la vitesse initiale du jet, etc... Dans un dernier temps, nous présentons une méthode hybride de couplage entre des simulations eulériennes et nos simulation lagrangiennes pour traiter des cas industriels. Après une validation de ce méthode de couplage avec décomposition de domaine, nous présentons des simulation complètes d'un avion en phase d'atterrissage. Enfin, nous mettons en évidence un cas de ré ingestion déjà observé expérimentalement en soufflerie.

LA TRANSITION A LA TURBULENCE D'UN SILLAGE DERRIERE UN OBSTACLE (BI OU TRIDIMENSIONNEL) EST ETUDIEE PAR SIMULATION NUMERIQUE DES EQUATIONS DE NAVIER-STOKES INSTATIONNAIRES. LE CODE DE CALCUL QUI A ETE CONCU DANS CE BUT UTILISE UNE METHODE MIXTE PSEUDO-SPECTRALE/DIFFERENCES FINIES D'ORDRE ELEVE. DANS UN PREMIER TEMPS, ON ETUDIE DE FACON DETAILLEE LE DEVELOPPEMENT INITIAL LINEAIRE DES INSTABILITES INFLEXIONNELLES ET LE ROLE JOUE PAR UNE STRATIFICATION STABLE EN DENSITE. LES CARACTERISTIQUES DU CHAMP D'ONDES DE GRAVITE INTERNES, A L'EXTERIEUR DU SILLAGE, SONT EXAMINEES A LA LUMIERE DES RESULTATS DE LA THEORIE LINEAIRE DANS L'APPROXIMATION WKB. ON A PU, DE LA SORTE, IDENTIFIER DIFFERENTS REGIMES HYDRODYNAMIQUES EN FONCTION DU NOMBRE DE FROUDE. LE DEVELOPPEMENT ULTERIEUR NON LINEAIRE DES INSTABILITES ET LA FORMATION DE STRUCTURES TOURBILLONNAIRES DISTINCTES SONT ANALYSES TANT D'UN POINT DE VUE TOPOLOGIQUE QUE QUANTITATIF. L'ACCENT EST MIS SUR LES MECANISMES D'INSTABILITE SECONDAIRE QUI ABOUTISSENT A L'EMERGENCE DE STRUCTURES LONGITUDINALES DANS UNE ALLEE TOURBILLONNAIRE DE KARMAN. CETTE ETUDE EST REALISEE EN TRANSITION NATURELLE OU FORCEE, AINSI QU'A GRAND NOMBRE DE REYNOLDS A L'AIDE D'UNE MODELISATION ORIGINALE DES ECHELLES SOUS-MAILLE. EN PRESENCE D'UNE STRATIFICATION, ON S'EST INTERESSE AU CHAMP D'ONDE DE GRAVITE EMIS PAR UNE ALLEE TOURBILLONNAIRE DANS LE SILLAGE D'UNE SPHERE. LA SIGNATURE DU CHAMP D'ONDE AINSI QUE LES POSSIBILITES DE PROPAGATION A L'EXTERIEUR DU SILLAGE ONT ETE CARACTERISEES POUR DIFFERENTS NOMBRES DE FROUDE