Lorsqu'un fluide en écoulement ou au repos est chauffé par le bas, l'instabilité thermique apparaît permettant une dissipation du gradient de température au-delà de ce qui est possible par simple conduction ou diffusion. Le fluide est soumis à des forces qui permettent un écoulement convectif lié à la différence de température entre les faces inférieure et supérieure de la couche. Le problème est rendu encore plus complexe par le fait que la distribution de température est elle-même déterminée par l'écoulement convectif. En d'autres termes, la force qui engendre l'écoulement est elle-même modifiée par l'écoulement. Pour un fluide multi composants, et avec prise en compte de l'effet SORET le phénomène devient d'une complexité particulièrement ardue, car la séparation des constituants engendre une compétition entre les forces volumiques d'origine thermique et compositionnelle. Le phénomène admet des développements extrêmement larges et de multiples applications en science et en ingénierie. En effet, l'instabilité thermique peut être étudiée avec tension superficielle variable, sans approximation de BOUSSINESQ, avec génération de chaleur au sein du fluide, avec ou sans confinement, etc. Ici, par souci de simplicité, nous étudierons le cas particulier d'une couche de mélange de fluides binaires entre deux plaques horizontales à températures fixées influencée par un écoulement forcé à faible nombres de PECLET. Nous commencerons par une étude qualitative du phénomène via une analyse de stabilité linéaire marginale et transitoire par éléments finis, laquelle démarche nous a conduit à réduire sous forme diagonale un opérateur linéaire rigoureusement déterminé. Les éléments spectraux de cet opérateur contiennent les informations utiles à déterminer les plus déstabilisants parmi une infinité de modes propres. Ces résultats nous serviront comme guide afin de déterminer une démarche d'expérimentation numérique par usage d'un code 3D en volumes finis hautement précis en espace, et utilisant le schéma d'EULER en temps, ouvrant la porte à des études qualitative spatiale et quantitative ultérieures. <br />L'ensemble de ces opérations nous ont conduit à déterminer le comportement d'un fluide binaire en convection mixte. Des relations littérales déterminant les seuils de transitions entre les différents états du système ont été démontrées. Des diagrammes de stabilité généraux ont été établis ; la mécanique de la convection est mieux éclaircie, l'influence de l'allongement du canal a été étudiée, la propagation dans les deux sens des ondes progressives des rouleaux transversaux ont été mises en évidence, un recueil de paramètres caractéristiques de ces écoulements est donné en termes de nombre d'ondes, pulsations, vitesses de phase et de transferts de chaleur et de masse. Pour finir, nous nous sommes intéressés au cas particulier de la séparation en essayant de mettre au point des idées permettant de séparer des constituants par appui sur les résultats trouvés.

Cette thèse est consacrée à l'étude des phénomènes aéroélastiques dans les tuyères de propulseurs. Elle s'inscrit dans le cadre des activités du groupe de recherche Aérodynamique des Tuyères et Arrière-Corps (ATAC). Un des principaux objectifs de ce travail est de démontrer que les phénomènes aéroélastiques peuvent être nuisibles pour les moteurs fusée (présence de modes instables) au même titre que les charges latérales qu'ils peuvent amplifier. Dans cet objectif, deux modèles permettant de prédire la stabilité aéroélastique de divergents flexibles en présence un choc de décollement sont détaillés. Par ailleurs, un outil numérique, reposant sur un couplage faible des codes fluide et structure dont dispose le LMFN-CORIA, est développé pour étudier ce type de problème. Cet outil performant est ensuite validé sur le cas-test classique du "flutter". Une application et une analyse approfondie des modèles de stabilité aéroélastique sont réalisées. Les modèles sont employés pour étudier la stabilité d'un divergent bidimensionnel plan soumis à un choc de décollement. L'influence de plusieurs paramètres sur la stabilité de la structure est évaluée. Les résultats obtenus sont ensuite confrontés aux résultats de simulations numériques couplées. En dernier lieu, les modèles sont utilisés pour étudier la stabilité de la tuyère LEA souple sur laquelle l'ONERA effectue des mesures expérimentales. Enfin, des travaux complémentaires sur l'amorçage rapide d'une tuyère supersonique sont présentés en annexes.

Dans cette thèse, nous nous intéressons à des méthodes numériques pour un modèle d'écoulements incompressibles et miscibles ayant des application dans l'hydrogéologie et l'ingénierie pétrolière. Nous étudions et analysons un schéma numérique combinant une méthode d'éléments finis mixtes (EFM) et une méthode des volumes finis (VF) pour approcher le système couplé entre une équation elliptique (pression-vitesse) et une équation de convection-diffusion-réaction (concentration). Le schéma VF considérée est de type "vertex centred" semi-implicite en temps : explicite pour la convection et implicite pour la diffusion. On utilise un schéma de Godunov pour approcher le terme convectif et une approximation élément fini P1 pour le terme de diffusion. Nous montrons des résultats de stabilité L≂ estimations BV et le principe du maximum discret sous une condition CFL appropriée. Ensuite, nous montrons la convergence de la solution approchée obtenue par le schéma combiné EFM-VF vers la solution du problème couplé. La démonstration de la convergence se fait en plusieurs étapes : premièrement, on déduit la convergence forte de la solution approchée de la concentration dans L2(Q), en utilisant la stabilité L≂, les estimations BV et des arguments de compacité. Dans l'étape suivante, on étudie le schéma découplé EFM, en donnant des résultats de convergence pour la pression et la vitesse. Enfin, le processus de convergence de la solution approchée du schéma combiné EFM-VF vers la solution exacte est obtenu par passage à la limite et par unicité de solution pour le problème continu. Des simulations numériques académiques et réalistes pour des problèmes bidimensionnels confirment la stabilité et l'efficacité du schéma combiné. Enfin, nous étudions des estimateurs d'erreur a posteriori de type résiduel pour une équation de convection-diffusion-réaction discrétisée par un schéma VF "vertex centred" semi-implicite en temps. Nous introduisons deux sortes d'indicateurs. Le premier est local en temps et en espace et constitue un outil efficace pour l'adaptation du maillage à chaque pas de temps. Le second est global en espace mais local en temps et peut être utilisé pour l'adaptation en temps. Nous montrons que l'estimateur est une borne supérieure de l'erreur. Des résultats numériques d'adaptations de maillage sont présentés et montrent l'efficacité de la méthode. La partie logiciels de ce travail porte sur deux volets. Le premier a permis de réaliser un code de calcul 2D, MFlow, écrit en C++, pour la résolution du système des écoulements miscibles considérés dans cette thèse. Le second volet concerne la collaboration avec un groupe de chercheurs pour l'élaboration de la plate-forme Homogenizer++ réalisée dans le cadre du GDR MoMaS (http://momas.univ-lyon1.fr/).

Un bassin sédimentaire est un milieu poreux de grande dimension (plusieurs dizaines de kilomètres de long et de large pour une profondeur d'environ cinq kilomètres) qui évolue au cours du temps par les effets de compaction et de sédimentation. Au cours de cette évolution, des hydrocarbures vont se former et s'écouler dans le bassin. On établit alors un modèle permettant de simuler cette évolution de bassin ainsi que la création, la migration et le piégeage des hydrocarbures dans des roches appelées roches magasins. Ces phénomènes se déroulant sur des centaines de millions d'années, on s'est attaché à étudier principalement une discrétisation temporelle de ces équations. On a ainsi mis en avant un raffinement local du pas de temps dont le principe est de recalculer la solution sur une zone jugée "mauvaise". A l'extérieure de cette zone, la solution est admissible. La difficulté vient de la détermination de la zone qui doit - être suffisamment "grande" pour avoir une bonne qualité de la solution, mais suffisamment "petite" pour obtenir un gain calcul. Les estimateurs a posteriori permettent de contourner cette difficulté. On a donc entrepris une étude théorique de ces estimateurs a posteriori dans le cas des équations linéaires elliptique et parabolique. Des simulations numériques montrent l'efficacité de ces estimateurs dans des cas académiques.

Le but de cette thèse est d'étudier et de calculer les pertes par courants de Foucault dans les transformateurs de distribution, et principalement dans leurs cuves, par la méthode des éléments finis en trois dimensions. La modélisation de ces appareils s'est heurtée à des difficultés liées : - aux faibles épaisseurs des tôles : ainsi, leur maillage avec des éléments finis volumiques conduit à un grand nombre d'éléments. - aux effets de peau qui peuvent être prononcés. Les difficultés de maillage sont encore accrues, à cause des variations rapides des champs dans la peau. - aux non linéarités de la caractéristique magnétique des matériaux. Compte tenu de ces difficultés, nous avons dû développer des méthodes numériques adaptées, qui sont : - différents types d'éléments spéciaux, permettant de diminuer les problèmes de maillage. En magnétostatique, ils permettent de décrire des entrefers minces et des régions ferromagnétiques minces. En magnétodynamique, ils permettent la modélisation de régions minces conductrices présentant une épaisseur de peau même faible par rapport à l'épaisseur de la région, de conducteurs massifs dans lesquels l'effet de peau est prononcé, décrits avec la notion d'"impédance de surface", de régions minces faiblement conductrices et de régions conductrices de type filaire. - des méthodes prenant en compte la saturation des matériaux magnétiques en régime harmonique (sinusoïdal). Une première méthode, pour des éléments volumiques, s'appuie sur une équivalence énergétique et permet de calculer uniquement des grandeurs globales. Une autre méthode, utilise la formule de l'impédance de surface calculée pour une caractéristique du matériau magnétique supposée en échelon. Les éléments spéciaux et les méthodes de prise en compte des non linéarités ont été validés par comparaison à des solutions analytiques, à des modélisations bidimensionnelles ou axisymétriques. La simulation d'un transformateur réel a été traitée pour déterminer les pertes par courants de Foucault dans la cuve.

Cette thèse porte sur le développement de schémas semi-implicites à pas fractionnaires pour les équations de Navier-Stokes compressibles ; ces schémas entrent dans la classe des méthodes de correction de pression.La discrétisation spatiale choisie est de type "à mailles décalées :éléments finis mixtes non conformes (éléments finis de Crouzeix-Raviart ou Rannacher-Turek) ou schéma MAC classique.Une discrétisation en volumes finis décentrée amont du bilan de masse garantit la positivité de la masse volumique.La positivité de l'énergie interne est obtenue en discrétisant le bilan d'énergie interne continu, par une méthode de volumes finis décentrée amont, enfin, et en couplant ce bilan d'énergie interne discret à l'étape de correction de pression.On effectue une discrétisation particulière en volumes finis sur un maillage dual du terme de convection de vitesse dans le bilan de quantité de mouvement et une étape de renormalisation de la pression; ceci permet de garantir le contrôle au cours du temps de l'intégrale de l'énergie totale sur le domaine.L'ensemble de ces estimations a priori implique en outre, par un argument de degré topologique, l'existence d'une solution discrète. L'application de ce schéma aux équations d'Euler pose une difficulté supplémentaire.En effet, l'obtention de vitesses de choc correctes nécessite que le schéma soit consistant avec l'équation de bilan d'énergie totale, propriété que nous obtenons comme suit. Tout d'abord, nous établissons un bilan discret (local) d'énergie cinétique.Ce dernier comporte des termes sources, que nous compensons ensuite dans le bilan d'énergie interne. Les équations d'énergie cinétique et interne sont associées au maillage dual et primal respectivement, et ne peuvent donc être additionnées pour obtenir un bilan d'énergie totale ; cette dernière équation est toutefois retrouvée, sous sa forme continue, à convergence : si nous supposons qu'une suite de solutions discrètes converge lorsque le pas de temps et d'espace tendent vers 0,, nous montrons en effet, en 1D au moins, que la limite en satisfait une forme faible.Ces résultats théoriques sont confortés par des tests numériques.Des résultats similaires sont obtenus pour les équations de Navier-Stokes barotropes.

Un certain nombre d'applications électrostatiques, notamment les dépoussiéreurs électrostatiques, mettent en oeuvre le phénomène de décharge couronne dans les gaz à partir d'électrodes ayant une forte courbure. Ces électrodes injectent des charges dans le gaz et la charge d'espace qui en résulte modifie la distribution du champ électrique. Le but de notre travail a été de développer un code numérique pour déterminer les distributions de champ électrique et de densité de charge dans la configuration lame-plan, la création de charge se faisant par effet couronne sur le bord de la lame. <br />Nous utilisons la méthode des éléments finis MEF pour résoudre l'équation de Poisson et la méthode des caractéristiques MC pour résoudre l'équation de conservation de la charge. Les deux équations couplées sont résolues par approximations successives en redéfinissant le maillage structuré à chaque itération pour éliminer la diffusion numérique. La redéfinition du maillage structuré est la partie la plus importante du travail. L'algorithme converge bien pour différentes distributions de la charge définies sur la lame injectrice. Les solutions numériques obtenues pour une loi d'injection imposée sur la lame plate donnent des résultats en très bon accord avec les mesures de densité de courant sur l'électrode plane d'un dispositif lame-plan.

La modélisation numérique des décharges plasma joue un rôle important dans la compréhension des mécanismes physiques ou chimiques ayant lieu dans les dispositifs assistés par plasma. Une grande partie de ces mécanismes est déjà prise en compte dans les codes actuels. En revanche, beaucoup d'entre eux ne permettent pas de travailler avec des géométries complexes. Cette limitation provient essentiellement de l'utilisation de maillages structurés, cartésiens. Ceux-ci ne sont pas bien adaptés aux géométries courbes. Les calculs en maillages structurés deviennent rapidement compliqués et spécifiques à une géométrie donnée. Notre travail concerne la modélisation de décharge pour un réacteur de traitement à la pression atmosphérique développé par Dow Corning. Sa configuration complexe ainsi que ses grandes dimensions nous ont incités à faire un nouveau code fonctionnant en maillages non structurés. Celui-ci doit être capable de s'adapter à la présence d'une pointe, d'arrondis et de multiples diélectriques mais aussi permettre le passage rapide à de nouvelles géométries. De plus ses grandes dimensions nécessitent l'utilisation de maillages raffinés uniquement aux endroits nécessaires (pointe, surfaces des diélectriques...). Le modèle mathématique utilisé est basé sur l'équation de Poisson couplée aux équations de transport de type dérive-diffusion. Plusieurs discrétisations numériques ont été testées dans des configurations physiques différentes. Nous présentons et validons les méthodes numériques choisies. Les résultats obtenus pour le réacteur Dow Corning sont alors exposés et commentés.

Dans ce rapport nous décrivons les nouvelles possibilités offertes par les algorithmes d'optimisation génétiques dans le domaine de l'électrotechnique. Après avoir analysé les différentes méthodes d'optimisation existantes, nous mettons en évidence leur points forts et leurs points faibles en les comparant sur différents cas tests. Les conclusions et les constations issues de ces confrontations nous ont guidé pour développer un algorithme d'optimisation perfonnant c'est-à-dire à la fois capable de localiser l'optimum global et peu coûteux en nombre d'évaluations de la fonction à optimiser. L'introduction d'infonnations supplémentaires concernant la "nature" des paramètres du problème traité s'est avérée fondamentale pour les algorithmes génétiques et ce point a été abordé car nous n'avons pas voulu limiter nos optimisations à un domaine bien particulier de l'électrotechnique. Les algorithmes d'optimisation mis au point ont été validés sur trois applications distinctes: - Optimisation de la forme d'un refroidisseur pour composant de puissance. - Optimisation de maillages tridimensionnels pour des logiciels éléments fInis . - Optimisation de la forme d'un dispositif électromagnétique composé de bobines supraconductrices.

Les transferts de chaleur par convection forcée turbulente dans une conduite plane partiellement remplie par un milieu poreux sont étudiés numériquement. L'étude concerne l'analyse de l'encrassement dans un canal plan représentatif d'un échangeur de chaleur à plaques. Un fluide, ayant un fort pouvoir encrassant, est considéré en régime turbulent. L'objectif de cette étude est de proposer une technique qui repose sur l'utilisation de matériaux fibreux comme capteur de particules pouvant réduire les méfaits de l'encrassement. Cela consiste à essayer de réduire la résistance d'encrassement en agissant sur les propriétés thermiques du dépôt. L'étude de la cinétique de l'encrassement permet de déterminer la loi de variation de l'épaisseur du dépôt au cours du temps. Cette équation est couplée aux équations de conservation. Un modèle de conductivité thermique effective (fluide, dépôt, fibres poreuses) a été choisi et le phénomène de colmatage de la matrice poreuse est considéré. L'apport du milieu poreux sur les performances de l'échangeur est analysé