Les modèles emboîtés sont largement utilisés en météorologie et en océanographie. Ils permettent un accroissement local de la résolution, dans les zones où cela semble nécessaire, via l'intégration d'un même modèle sur une hiérarchie de grilles. Dans le cas d'interaction one-way, les conditions aux frontières pour la grille fine proviennent d'une interpolation de la solution obtenue sur la grille à plus faible résolution. Dans le cas d'interaction two-way, une rétroaction de la grille fine vers la grille grossière est ajoutée. Toutefois, le problème de l'assimilation variationnelle de données dans de tels systèmes n'a pas, ou peu, été étudié à ce jour. Ces classes de méthodes, notamment l'algorithme 4D-Var, permettent d'améliorer la solution d'un modèle, jusqu'ici mono-grille, en minimisant une fonctionnelle mesurant l'écart de ce modèle aux observations présentes sur une fenêtre temporelle. Le travail présenté ici vise donc à formuler un algorithme d'assimilation 4D-Var localement multi-grille. Pour le cas général d'une grille haute résolution emboîtée localement dans une autre à plus faible résolution, nous posons les équations du système adjoint dans les deux cas d'interactions one-way et two-way. Nous montrons ainsi que la formulation adjointe fait naturellement apparaître de nouvelles interactions entre les grilles, dans le sens opposé de celles existant dans la formulation directe. De plus, nous proposons différentes variantes à ces algorithmes, réalisant un couplage faible entre les solutions des différents modèles via l'ajout d'un terme de contrôle au niveau des transferts inter-grilles. Nous présentons également l'application d'une méthode multi-grille, le Full Approximation Scheme, à l'assimilation variationnelle de données. Cette approche permet d'obtenir un algorithme d'assimilation multi-grille potentiellement très efficace. Enfin, ces méthodes sont testées sur le cas d'un modèle Saint Venant 2D. Nous constatons une réduction importante des erreurs des solutions multi-grilles, ainsi qu'une accélération de la convergence de ces algorithmes.

Les premiers quatre chapitres traitent de l'interaction fluide-structure stationnaire. On étudie l'interaction évolutive en temps dans les chapitres cinq et six. Les deux derniers chapitres sont consacrés aux écoulements à frontière libre avec tension de surface qui ont certaines similitudes avec les problèmes d'interaction fluide-structure.<br /><br />Le fil directeur de mes travaux est de prendre comme ``contrôle'' une partie des conditions aux limites à l'interface et ``d'observer'' si les conditions de couplage sont vérifiées. En traitant l'observation par la méthode de moindres carrés, on obtient des problèmes de type contrôle optimal. Dans le chapitre 1, on prouve que la fonction coût est semi-continue inférieurement et en conséquence, on peut démontrer l'existence d'un contrôle optimal. On prouve la différentiabilité de la fonction coût, et on donne la forme analytique du gradient dans le chapitre 2. On présente également des résultats numériques. Dans le chapitre 3 on étudie la sensibilité du problème et on donne la forme analytique du gradient sans faire appel à l'état adjoint. Des résultats numériques sont obtenus. Dans le chapitre 4, pour résoudre le problème du fluide, on prescrit la composante normale de la vitesse du fluide et la composante normale des forces de surface. C'est une formulation rarement utilisée pour résoudre les équations de Stokes. On cherche à minimiser la composante tangentielle de la vitesse du fluide à l'interface. On prouve que le problème fluide est bien défini et on présente des résultats numériques. Dans le chapitre 5, on introduit un algorithm où on doit résoudre à chaque pas de temps un problème de minimisation. C'est un algorithme bien adapté notamment quand le fluid est pulsatif. On présente des résultats numériques pour des pas de temps relativement grand.<br /><br />Un résultat de convergence concernant un algorithme pour des maillages dynamiques est présenté dans le chapitre 6. Dans les chapitres 7 et 8, on veut déterminer numériquement l'évolution d'un domaine bidimensionnel avec application au développement cellulaire. L'écoulement du fluide dans le domaine en mouvement dépend de la tension de surface à la frontière libre. Cette tension est proportionnelle à la courbure de la frontière. Les algorithmes employés sont de type ``front-tracking''. Des résultats numériques sont présentés.

Cette thèse porte sur la simulation de la propagation et diffraction d'ondes dans un milieu élastique anisotrope hétérogène fissuré à l'aide de méthodes numériques explicites. L'objectif est de développer une méthode numérique performante capable de prendre en compte les détails géométriques ou singularités de la solution de manière précise. Les deux premières parties sont consacrées à des méthodes de raffinement de maillage spatio-temporel. Adapter le pas de temps localement au pas d'espace permet en même temps de diminuer la dispersion numérique dans la grille grossière et de gagner en temps de calcul. Les méthodes proposées sont conservatives, ce qui garantit la stabilité des schémas numériques. La géométrie des fissures est prise en compte par la méthode des domaines fictifs. La troisième partie présente un nouvel élément fini qui garantit la convergence de cette méthode. La dernière partie décrit le couplage entre les techniques de raffinement et la méthode de domaines fictifs.

Dans ce travail nous avons ouvert des voies nouvelles pour la méthode des équations intégrales de frontière appliquée à la résolution des problèmes tridimensionnels de l'électromagnétisme en basse fréquence. Elles permettent de réaliser un maillage adaptatif associé à des résolutions approchées, locales, du système intégral. Dans le premier chapitre nous discutons brièvement les différentes méthodes numériques appliquées à l'électromagnétisme. Le deuxième chapitre est consacré à la méthode des équations intégrales de frontière. Nous y justifions nos choix particuliers (type et ordre des éléments discrétisant les frontières) et décrivons les algorithmes réalisés, par exemple pour la prise en compte des symétries ou périodicités des systèmes étudiés. Dans le troisième chapitre, nos méthodes de calcul des intégrales élémentaires, combinant l'analytique et le numérique, sont décrites et validées. Dans le chapitre quatre, nous présentons un estimateur d'erreur basé sur l'écart entre la valeur interpolée du potentiel et sa valeur calculée par la méthode intégrale elle-même. Pour finir, nous utilisons dans le chapitre cinq cet estimateur d'erreur pour réaliser automatiquement un. maillage adapté à chaque problème traité, en autorisant aussi un maillage « non-conforme ». Une méthode de résolution locale permet un gain important en temps de calcul pendant cette phase de maillage adaptatif. Le maillage final est séparé en domaines, ce qui permet d'améliorer le résultat par itérations entre ces domaines, et sans recourir à une résolution globale, très coûteuse pour des problèmes de grande taille.

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.

Dans le cadre de la simulation du procédé de consolidation de pièces céramiques ou composites CerMet par frittage, cette thèse se propose de modéliser et de simuler deux aspects particuliers du frittage. Le premier consiste à prendre en compte l'effet d'éléments dopants sur l'évolution microstructurale d'un élément volumique représentatif d'un compact de poudre. Pour cela les chemins de diffusion tels que la diffusion volumique et de surface sont pris en compte. Le second aspect est en relation avec l'élaboration de matériaux composites à dispersoïdes et l'étude de l'évolution de leur microstructure par simulation numérique selon les caractéristiques des inclusions ne participant pas au processus de frittage. Ces simulations reposent sur le développement de la méthode Level-Set qui permet de suivre tous les changements de topologie qui interviennent au cours du frittage à l'échelle des grains.Des résultats concluants ont été obtenus à partir de la simulation de l'effet du dopage, dans le cas particulier de l'alumine dopée magnésie, puis généralisés au cas du codopage. Pour les multimatériaux, la première étape de la modélisation a consisté à considérer des inclusions ayant le même type de loi de comportement que la matrice céramique, à savoir élastique, linéaire et isotrope. Cependant, les propriétés des matériaux seront différentes. Du point de vue numérique, la complexité du problème a consisté à bien gérer les deux phases solides, et en particulier leur interface commune en vue de la résolution de l'équilibre mécanique par éléments finis, en utilisant plusieurs fonctions Level-Set.

Les écoulements vrillés, qui sont largement utilisés dans les turbines à gaz, sont connus pour être sujet à des bifurcations entre différentes topologies (grandes reconfigurations de l'écoulement) qui peuvent affecter les performances et la sécurité du moteur. Ce travail se concentre sur l'étude de ces bifurcations en utilisant la Simulation aux Grandes Echelles (SGE). Cette étude montre qu'un petit changement dans les conditions dynamique du fluide, induite par les différents modèles de sous-maille utilisés, peut provoquer une transition entre deux régimes d'écoulement distincts lorsque l'écoulement tourbillonnaire est proche des conditions critiques de transition. La sensibilité de la SGE aux modèles de sous-maille est également identifiée comme le résultat d'un manque de résolution à certains endroits critiques, un problème qui est analysé en utilisant une méthode d'adaptation de maillage. L'adaptation de maillage est testée sur des cas académiques et industriels. Ici, par ajustement de la résolution du maillage sur la base des caractéristiques de l'écoulement étudié (raffinement et grossissement de la grille en maintenant constant le coût numérique), des améliorations substantielles peuvent être obtenues, en terme de prédictions de la SGE. Ce travail peut être considéré comme une des premières étapes vers la mise en place d'une procédure standard (reproductible et indépendante de l'utilisateur) de maillage pour la SGE.

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 progrès en imagerie 3D ont contribué à l'augmentation de leurs utilisations dans de nombreuses applications. L'évolution des débits des lignes internet ont rendu les transmissions d'objets 3D de plus en plus importantes. Les travaux sur la protection des objets 3D sont donc un sujet d'actualité. Cette thèse s'inscrit dans le cadre d'une collaboration entre la société STRATEGIE et le LIRMM (Laboratoire d'informatique, de robotique et de microélectronique de Montpellier). Nous abordons deux sujets manipulant des objets 3D. Le premier sujet abordé est la digitalisation de forme intérieure de chaussure. L'algorithme de digitalisation repose sur une méthode de \textit{shape from silhouettes}. L'acquisition des silhouettes se fait à l'aide d'un appareil photo numérique et d'un plateau rotatif sur lequel repose l'objet à digitaliser. Le second sujet abordé est le tatouage d'objet 3D. Nous proposons, dans cette thèse, deux méthodes d'insertion de données cachées ne s'appuyant pas sur la modification des sommets des objets 3D pour réaliser l'insertion. Ces méthodes permettent d'insérer des méta-données dans des objets 3D comme le nom de l'auteur, la date de création de l'objet, un logo 2D ou 3D, la couleur de l'objet ou encore la texture de l'objet lui même. L'idée principale des deux méthodes est de trouver et de synchroniser des zones particulières pouvant être utilisées pour insérer le message. L'insertion de données s'appuie sur la modification de la connexité des arêtes dans des zones sélectionnées. Ces modifications ont pour conséquence de modifier la structure des triangles construits dans ces zones. Ces méthodes de dissimulation présentent l'avantage de ne pas modifier la position initiale des points du modèle 3D. Cet invariant sur la position des points nous permet de rendre l'insertion robuste aux transformations affines de type rotation, translation ou changement d'échelle.

Le traitement la forme d'un produit et la formalisation de l'information impliquée dans le processus de liaison entre les activités de Conception et d'Analyse Mécanique sont essentiels. Dans ce travail, un cadre générique la définition d'un modèle intermédiaire est proposé. Le modèle de forme de référence est polyédrique mais également enrichi avec une description de plus haut niveau. La représentation mixte de forme qui en découle permet de transférer quelques hypothèses et objectifs du problème considéré en termes de contraintes qui doivent satisfaites pendant la génération du modèle intermédiaire, où, plus tard dans le processus, pendant la génération modèle de simulation où effectuer l'analyse mécanique. De plus, la relation qu'il y a entre l'application d'un ( mécanique et la suppression d'un sous-domaine de forme est mise en évidence. En particulier, la formulation critère mécanique a posteriori supportant les transformations de forme est fournie.