Les turbines aéronautiques doivent satisfaire à des normes d'émissions polluantes toujours en baisse. La qualité du mélange du carburant et de l'air dans la chambre de combustion est responsable de la formation de polluants nocifs pour l'environnement. La simulation aux grandes échelles (LES) permet d'étudier les mécanismes de mélanges turbulents de l'air et du carburant. La prise en compte de l'aspect liquide du carburant injecté devient nécessaire pour prédire correctement l'apparition de vapeur de carburant au sein du foyer. Le but de cette thèse est évaluer la fiabilité des simulations LES Euler-Euler dans une configuration complexe. Les processus d'injection, et d'évaporation du carburant liquide sont analysés et modélisés dans les simulations LES car ils pilotent la formation de vapeur de carburant. Les méthodes numériques pour résoudre les équations continues de la phase dispersée doivent permettre des simulations précises et robustes dans une configuration représentative d'une chambre de combustion. Les simulations présentées dans ces travaux reproduisent l'écoulement diphasique évaporant non-réactif du banc d'essai Mercato. Ce banc est équipé d'un système d'injection d'air vrillé et d'un atomiseur pressurisé-swirlé de kérosène typiques des foyers aéronautiques réels. Dans ces travaux, le modèle pour l'injection de liquide FIM-UR a été développé pour définir les conditions limites conduisant à un spray issu d'un atomiseur préssurisé-swirlé. Le kérosène employé dans les campagnes expérimentales est modélisé dans les simulations par un composé permettant d'obtenir des temps d'évaporation réalistes. Trois stratégies numériques ont été mises en place sur la configuration MERCATO. Les comparaisons des résultats numériques aux mesures expérimentales ont permis d'évaluer la stratégie numérique conduisant à la meilleure précision. L'utilisation du schéma centré TTGC associé à un opérateur de viscosité artificielle localisée par un senseur adapté est optimale lorsque l'équation sur l'énergie décorrélée des gouttes est résolue. Cette stratégie permet de contrôler la localisation et les niveaux de viscosité par rapport à un schéma décentré. Les termes sources liés au mouvement mésoscopique permettent de redistribuer l'énergie dans les zones de compression ou de détente de la phase dispersée, et d'obtenir les bonnes répartitions des fluctuations dans la chambre de combustion. La stratégie retenue est comparée aux statistiques de la dynamique du spray résolu par une approche Lagrangienne employant la même injection monodispersse. Le méthode Euler-Euler conduit à la même précision de la dynamique de la phase dispersée que la méthode Euler-Lagrange. L'accès à l'évolution instationnaire de l'écoulement permet d'identifier les mêmes mécanismes de dispersion et de mélange dans les deux simulations. Des différences sur la répartition de diamètre moyen et de carburant dans la chambre ont été mis en évidence et reliés à la polydispersion locale qui n'est pas résolue dans l'approche Euler-Euler monodisperse et qui apparaît naturellement dans l'approche Euler-Lagrange malgré l'injection monodisperse.

Les mesures des niveaux d'expression de tous les gènes d'un génome requièrent une analyse<br />statistique afin d'obtenir des conclusions fiables. Les biologistes ont du mal à faire un choix dans la<br />foule de méthodes existantes. Afin de déterminer quelle méthode est la plus adéquate pour la<br />problématique abordée, des comparaisons de méthodes d'analyse disponibles sont nécessaires.<br />Actuellement les critères de comparaison se révèlent soit lacunaires ou soit non pertinents du point de<br />vue biologique.<br />Nous avons introduit un nouveau critère biologique de comparaison des méthodes d'analyse du<br />transcriptome fondé sur une structure des génomes bactériens : les opérons. Les gènes d'un opéron<br />sont généralement transcrits sur un même ARNm. Si un gène d'un opéron bactérien est identifié, les<br />autres gènes de l'opéron devraient l'être également. Nous avons ainsi comparé des méthodes<br />d'analyse appliquées au transcriptome : l'ACP et l'ACI, respectivement analyses en composantes<br />principales et indépendantes, l'ANOVA, analyse de variance, la régression des moindres carrés<br />partiels PLS et différents t-tests. Chaque méthode aborde le nuage de données d'un point de vue<br />différent ce qui donne des résultats complémentaires. Globalement, l'ACI a fourni les meilleurs<br />résultats tant en sensibilité qu'en terme de précision.<br />Un autre aspect, en plein développement, de l'analyse du transcriptome est la méta-analyse de<br />données d'origines diverses malgré les biais inhérents à cette technologie. Généralement ces métaanalyses<br />visent à préciser les résultats concernant des gènes différentiellement exprimés ou coexprimés.<br />Elles ouvrent également la possibilité d'étudier de nouveaux champs en biologie. Nous<br />avons utilisé des données de transcriptome indépendantes afin d'étudier l'organisation de l'expression<br />des gènes et, ainsi, celle du chromosome bactérien. L'étude du transcriptome de trois bactéries, B.<br />subtilis, E. coli et S. meliloti a révélé des corrélations d'expression à longue distance valables quel que<br />soit le gène étudié. Les structures en opéron se manifestent clairement au travers de cette étude, qui<br />a également permis de préciser que la co-expression de gènes proches s'étend au-delà des opérons<br />dans une région qui se répand jusqu'à une centaine de gènes.<br />Pour conclure, l'analyse du transcriptome n'a pas réellement nécessité la mise au point de méthodes<br />d'analyse statistiques spécifiques. Cependant, elle permet d'aborder de nouveaux horizons dans la<br />biologie, et notamment l'organisation chromosomique du génome bactérien.

Les méthodes actuelles de calcul instationnaire des équations de Navier-Stokes (Simulation Numérique Directe (SND), Simulation des Grandes Echelles (SGE)) demandent encore un temps de calcul excessif pour pouvoir être utilisées à des fins industrielles. D'autre part, les méthodes statistiques (RANS) moins coûteuses, ne permettent pas de rendre compte des phénomènes à large bande fréquentielle, qui sont à l'origine, par exemple, des ondes acoustiques. Le développement de méthodes localement couplées permet de concilier les avantages des deux types d'approches précités (SGE et RANS). Une méthode originale de ce type à été développée sur la base de la Non-Linear Disturbance Equations (NLDE). Elle consiste à résoudre l'intégralité du champ grâce à un calcul statistique. Ensuite, aux endroits critiques, le complément fréquentiel est calculé en SGE, de manière à enrichir le champ. Le code a été implanté sur le NEC SX5 de l'ONERA et testé sur plusieurs cas d'application en incompressible (canal plan bi-périodique) et en compressible (aube de turbomachine, aile hypersustentée). Ses possibilités ont été testées et des conditions limites spécifiques ont été développées. La méthode se comporte bien dans tous ces cas d'application et peut être utilisée de façon locale, permettant ainsi une nette réduction du temps de calcul.

Les écoulements diphasiques turbulents sont présents dans de nombreux systèmes industriels (moteur à piston, turbines à gaz, moteurs fusée...). La compréhension fine de telles configurations s'avèrent de nos jours nécessaire pour limiter notamment les émissions de polluants et de gaz à effet de serre, et la consommation des énergies fossiles. Nous nous intéressons ici à la simulation aux grandes échelles des écoulements diphasiques turbulents, permettant de capturer une large partie du spectre de la turbulence, et ainsi être capable de prédire des phénomènes instables ou transitoires. La phase dispersée est ici modélisée par une approche eulérienne, en raison de ses avantages dans le contexte du calcul haute performance. Le travail de cette thèse a consisté à étendre le formalisme eulérien existant dans le code AVBP à la simulation de sprays polydisperses dans des écoulements turbulents. Pour cela, le Formalisme Eulérien Mésoscopique (FEM) a été couplé à une approche Multi-fluide. Cette nouvelle approche, intitulée Formalisme Eulérien Mésoscopique Multi-fluide (FEMM), a été évaluée sur des cas simples canoniques, permettant de bien caractériser le comportement autant en terme de dynamique turbulente que d'effets polydisperses. Les stratégies numériques disponibles dans le code de calcul AVBP sont aussi analysées, afin d'en cerner les limites pour la simulation eulérienne d'une phase liquide. Ce nouveau formalisme est finalement appliqué à la configuration aéronautique MERCATO, pour laquelle on dispose de résultats numériques obtenus avec d'autres approches (FEM et approche lagrangienne), et de résultats expérimentaux. Un accord satisfaisant avec l'expérience est montré pour toutes les approches, même si le FEM, monodisperse, obtient de moins bon résultats en terme de fluctuations. D'autres résultats expérimentaux s'avèrent nécessaires pour évaluer les approches et déterminer quelle est la plus prédictive pour cette configuration, notamment concernant la fraction massique de kerosene, autant en phase liquide qu'en phase gazeuse.

Cette thèse est consacrée à la simulation d'écoulements turbulents, incompressibles ou à faible nombre de Mach pour des applications touchant à la sûreté nucléaire. En particulier, nous nous concentrons sur le développement et l'analyse mathématique de schémas numériques performants pour la méthode dite de Simulation des Grandes Echelles. Ces schémas sont basés sur des méthodes à pas fractionnaires de type correction de pression et des éléments finis non conformes de bas degré. Deux arguments semblent essentiels à la construction de tels schémas: le contrôle de l'énergie cinétique et la précision pour des écoulements à convection dominante. Concernant la discrétisation en temps, nous proposons un schéma de type Crank-Nicolson et nous montrons qu'il satisfait un contrôle de l'énergie cinétique. Ce schéma présente de plus l'avantage d'être peu dissipatif numériquement (résidu d'ordre deux en temps). Concernant le défaut de précision de la discrétisation par l'élément fini de Rannacher-Turek, nous envisageons deux approches. La première consiste à construire un schéma pénalisé contraignant les degrés de liberté tangents aux faces des cellules à s'écrire comme combinaison linéaire des degrés de liberté normaux alentour. La deuxième approche repose sur l'enrichissement de l'espace discret d'approximation pour la pression. Enfin, différents tests numériques sont présentés en dimensions deux et trois et pour des maillages généraux, afin d'illustrer les capacités des schémas étudiés et de confronter les résultats théoriques et expérimentaux.

Une méthode de frontières immergées est développée pour la simulation d'écoulements compressibles et validée au travers de cas-tests spécifiques (réflexion d'ondes acoustiques et quantification de la conservation de la masse dans des canaux inclinés). La simulation aux grandes échelles (LES) d'une cavité transsonique est ensuite présentée. Le bouclage aéro-acoustique, très sensible aux conditions aux limites, est reproduit avec précision par la LES dans le cas où les parois sont immergées dans un maillage structurée. La comparaison des stratégies de modélisation de sous-maille pour cet écoulement transsonique et l'adaptation des filtres en présence de frontières immergées sont également discutées. Le rôle, souvent sous-estimé, du schéma de viscosité artificiel, est quantifié.Dans la dernière partie du manuscrit, des études sont réalisées pour aider au dimensionnement d'un nouveau concept de chambre de combustion où la flamme est stabilisée par la recirculation de gaz brûlés dans une cavité (chambre TVC pour Trapped Vortex Combustor). La modélisation de la combustion turbulente est basée sur une chimie tabulée, couplée à une fonction densité de probabilité présumée (PCM-FPI). L'étude de la dynamique de la flamme est réalisée pour diverses conditions de fonctionnement (débit de l'écoulement principal et présence ou non d'un swirl). Les spécificités de mise en œuvre de la simulation d'un écoulement de ce type sont discutées et un soin particulier est apporté au traitement de la condition de sortie, qui constitue un point sensible de la chaîne de modélisation. Les phénomènes d'instabilités et de retour de la flamme sont mis en évidence ainsi que les modifications à apporter au dispositif afin de minimiser ces effets. L'existence d'un cycle limite acoustique est souligné et une formule permettant d'anticiper le niveau des fluctuations de pression est proposée et validée. Une correction au modèle PCM-FPI est présentée afin de préserver la vitesse de flamme et d'assurer une reproduction plus précise de la dynamique de flamme.

Comme le groupe de symétries de Lie des équations aux dérivées partielles représentent les propriétés physiques intrinsèques contenues dans les équations, il offre un outil efficace pour étudier et modéliser les phénomènes physiques. Ainsi, dans cette thèse, on se propose d'appliquer la théorie du groupe de symétries de Lie à la modélisation des écoulements anisothermes.On calcule alors des lois de paroi, et, plus généralement des lois d'échelle, pour la vitesse et la température dans le cas d'un écoulement parallèle. En fait, ces lois d'échelle se révèlent être simplement des solutions auto-similaires des équations de Navier-Stokes moyennées par rapport aux symétries des équations.Ensuite, par l'approche de la théorie des groupes de Lie, on construit une classe de modèles de sous-maille qui sont invariants par les symétries des équations de Navier-Stokes anisothermes.Ces modèles ont l'avantage de respecter les propriétés physiques des équations qui sont contenues dans les symétries. De plus, par cette approche, le modèle de flux de chaleur apparaît naturellement,sans qu'on ait besoin de faire appel à la notion de nombre de Prandtl de sous-maille,ce qui augmente la portée de ces modèles par rapport à la plupart des modèles existants. Par ailleurs, le comportement proche de la paroi de certains des modèles proposés est étudié. Enfin,des tests numériques en convection naturelle et en convection mixtes sont effectués.

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.

Afin de pouvoir simuler la formation des principaux polluants au sein de foyers aéronautiques réalistes, un modèle de réduction de la chimie détaillée (FPI), basé sur la construction de tables à partir de calculs de flammes de prémélange laminaires élémentaires, est adapté et couplé au code d'aérothermochimie CEDRE de l'ONERA. Après une brève validation de ce modèle via la simulation de flammes laminaires canoniques, les interactions chimie/turbulence sont modélisées sous l'hypothèse des flammelettes, en approchant les PDF des paramètres d'entrée des tables par des fonctions beta. Cette approche complète est appliquée à la simulation numérique de l'écoulement au sein d'une configuration plus appliquée : la chambre PRECCINSTA. Ce cas bien connu a permis notamment l'évaluation des capacités du modèle dans un contexte plus industriel par comparaison des résultats de calcul aux données expérimentales disponibles. Il a en particulier permis de tester l'approche FPI étendue à la modélisation de la combustion partiellement prémélangée. Par ailleurs, l'utilisation d'un modèle de chimie réduite s'avère particulièrement appropriée pour prédire l'émission de substances polluantes, par exemple CO. Cependant, lorsque l'on considère la formation de NO, FPI ne peut pas être utilisé directement du fait de la lente dynamique chimique de cette espèce.Pour pallier à cette limitation, deux approches permettant de modéliser la production de NO au sein d'écoulements complexes sont proposées, fondées sur l'utilisation des tables chimiques FPI. Les capacités de ces modèles sont finalement analysées à l'aide de calculs effectués sur la configuration PRECCINSTA.

Cette thèse propose une modélisation aussi générique que possible de la combustion dans les moteurs automobiles dans un cadre de simulation aux grandes échelles. Une première étude aborde la fermeture du terme de transport non résolu pour la flamme de prémélange. Par la suite, un couplage entre les modèles ecfm-les et pcm-fpi est proposé et validé pour intégrer les effets de chimie complexe à la simulation de la flamme de prémélange. Ce travail est étendu par l'adjonction de modélisations spécifiques à l'allumage par bougie et de l'auto-inflammation, toujours en intégrant les effets de chimie détaillée. Des calculs d'application aux combustions anormales dans les moteurs à allumage commandé concluent ce travail.