En mécanique des fluides (CFD), l'adaptation de maillage anisotrope est reconnue pour sa capacité à réduire le ratio entre le nombre de degrés de liberté et la précision du calcul. Cependant, son application dans le cas d'écoulements compressibles avec des chocs pose les problématiques suivantes : (i) les schémas numériques d'ordre élevé de type shock capturing retombent à l'ordre un dans les chocs, (ii) les senseurs utilisés pour l'adaptation prescrivent dans les chocs des tailles qui tendent vers zéro. Il est donc nécessaire de prescrire une taille minimale. On perd alors tout l'intérêt d'une adaptation anisotrope. On apporte une réponse à ces problématiques en considérant une adaptation anisotrope multi-échelles du maillage basée sur le modèle de maillage continu. On alors montre que le processus adaptatif converge dans les chocs si le schéma numérique utilisé est non compressif. La prescription d'une taille minimale n'est plus nécessaire. On retrouve également un ordre deux de convergence dans tout le domaine, même en présence de chocs. Si on se donne des informations supplémentaires (fonctionnelle précise à observer, équation aux dérivées partielles, schéma numérique utilisé pour la résoudre) les méthodes génériques précédentes ne sont plus op- timales dans la distribution des degrés de liberté. On étudie cette problématique dans le cas particulier des équations d'Euler pour des fonctionnelles scalaires. Ce type d'étude est très bien adapté pour le calcul de grandeurs d'intérêt comme la portance ou la traînée en aérodynamique. On propose une estimation d'erreur a priori pour le contrôle de l'erreur d'approximation sur une fonctionnelle. Cette estimation est ensuite minimisée sur l'espace des maillages continus afin de décrire le maillage anisotrope optimal. Enfin, on applique l'adaptation multi-échelles à la prédiction haute-fidélité du bang sonique.

Ce travail concerne la mise en oeuvre d'une méthode numérique d'optimisation de type contrôle optimal appliquée à un problème d'assimilation de données en biogochimie marine. Après avoir présenté le systme d'équations aux dérivés partielles non-linéaires régissant l'évolution en temps et en espace des différentes variables physiques et biologiques, un premier travail, mathématique, a consisté à montrer l'existence, l'unicité et la positivité de la solution du modèle biologique. La seconde partie du travail est numérique. Le modèle est discrétisé par diffrences finies et les codes linéaire tangent et adjoint sont obtenus par différentiation automatique. Ces outils informatiques étant développés, on peut aborder le problème inverse d'assimilation variationnelle de données. Les variables de contrôle sont les paramètres intervenant dans les termes non-linéaires de réactions biologiques. On cherche un jeu de paramètres optimal minimisant une fonction cout. Celle-ci mesure l'écart au sens des moindres carrés entre les observations et les sorties correspondantes du modèle.Une étude de sensibilité préliminaire, utilisant le modèle tangent linéaire, ainsi que des expériences d'identification, utilisant le modèle adjoint, avec données simulées, sont menées. On utilise enfin la méthode pour assimiler des données réelles de la station Dyfamed en Méditerranée Nord-Occidentale.

Ce travail de thèse a pour but l'optimisation de forme en forgeage 3D. Les problèmes à résoudre consistent à chercher la forme optimale du lopin initial ou des outils de préforme afin de minimiser une fonction coût F qui représente une mesure de non-qualité définie par les industriels. Ce sont souvent des problèmes multi optima, et le temps nécessaire pour une évaluation de la fonction coût est très élevé (de l'ordre de la journée). L'objectif de cette thèse est de construire un module d'optimisation automatique qui permet de localiser l'extremum global à un coût raisonnable (moins de 50 calculs de la fonction coût à chaque optimisation). La simulation du procédé est effectuée avec le logiciel éléments finis FORGE3®. Les formes axisymétrique des pièces initiales ou des outils de préforme (dans le cadre du forgeage multi-passes) sont paramétrées en utilisant des polygônes quadratiques ou des courbes Bsplines. Différentes fonctions coûts sont considérées, comme l'énergie totale de forgeage ou la mesure non-qualité de la surface (défaut de repli). Le gradient de ces fonctions coûts est obtenu par la méthode de l'Etat Adjoint combinée avec la méthode de différentiation semi-analytique. Dans ce travail, afin d'aborder une famille de procédés de forgeage plus vaste, ce calcul du gradient (initié dans la thèse de M. Laroussi) a été étendu aux paramètres de forme des outils de préformes dans le cadre du forgeage multi passes.Différents algorithmes d'optimisation ont été étudiés : un algorithme BFGS standard, un algorithme de type asymptotes mobiles, une stratégie d'évolution couplée avec une surface de réponse basée sur le Krigeage et deux nouveaux algorithmes hybrides proposés dans le cadre de ce travail. Cette approche hybride consiste à coupler un algorithme génétique avec une méthode de surface de réponse pour réduire le nombre d'évaluations de la fonction coût. Tous les algorithmes étudiés sont comparés sur deux problèmes caractéristiques de forgeage 3D, respectivement l'optimisation de la géométrie de la préforme et celle des outils de préforme. Les résultats obtenus montrent la faisabilité de l'optimisation de forme en forgeage 3D, c'est-à-dire l'obtention de résultats satisfaisant en moins de 50 simulations 3D et la robustesse des algorithmes à base de méta-modèle.

Les changements du champ magnétique terrestre sur une grande gamme d'échelles spatiales et temporelles reflètent les processus variés de la éeodynamo. Je propose un modèle simplifié de la dynamique rapide du noyau, adapté à l'étude des variations du champ magnétique de l'année au siècle ; la variation séculaire. L'hypothèse quasi-géostrophique du modèle est basée sur la prépondérance des forces de rotation par rapport aux forces magnétiques à ces échelles de temps. La partie axisymétrique correspond au formalisme d'ondes de torsion d'Alfvén. La dynamique se place dans le plan équatorial. A la frontière noyau-manteau, l'écoulement interagit avec le champ magnétique radial via la composante radiale de l'équation d'induction. Cette partie du modèle connecte la dynamique et les observations. L'assimilation variationnelle de données permet d'interpréter la variation séculaire en terme de dynamique. Une fonction objectif est minimisée en calculant sa sensibilité par rapport aux variables de contrôle via l'intégration du modèle adjoint. J'illustre cette inversion par des expériences jumelles pour un écoulement stationnaire dans le noyau et pour des ondes de torsion. On accède ainsi à des variables d'état qui ne sont pas directement observées. En utilisant comme observations des écoulements reconstruits à la surface du noyau, cette méthode permet de déduire que la tension magnétique dans le noyau, force de rappel des ondes de torsion, correspond à un champ magnétique fort, au minimum 3-4 mT. De telles ondes de torsion rapides sont cohérentes avec un signal à 6 ans dans les données de variations de la longueur du jour.

<p>In recent years, there has been a great interest towards the higher order numerical methods for singularly perturbed problems. As compared to their lower order counterparts, they provide better accuracy with fewer mesh points. Construction and/or implementation of direct higher order methods is usually very complicated. Thus a natural choice is to use some convergence acceleration techniques, e.g., Richardson extrapolation, defect correction, etc. In this thesis, we will consider various classes of problems described by singularly perturbed ordinary and partial differential equations. For these problems, we design some novel numerical methods and attempt to increase their accuracy as well as the order of convergence. We also do the same for existing numerical methods in some instances. We &macr / nd that, even though the Richardson extrapolation technique always improves the accuracy, it does not perform equally well when applied to different methods for certain classes of problems. Moreover, while in some cases it improves the order of convergence, in other cases it does not. These issues are discussed in this thesis for linear and nonlinear singularly perturbed ODEs as well as PDEs. Extrapolation techniques are analyzed thoroughly in all the cases, whereas the limitations of the defect correction approach for certain problems is indicated at the end of the thesis</p>

No description available.

Conservation laws play an important role in science. The aim of this thesis is to provide an overview and develop new methods for constructing conservation laws using Lie group theory. The derivation of conservation laws for invariant variational problems is based on Noether’s theorem. It is shown that the use of Lie-Bäcklund transformation groups allows one to reduce the number of basic conserved quantities for differential equations obtained by Noether’s theorem and construct a basis of conservation laws. Several examples on constructing a basis for some well-known equations are provided. Moreover, this approach allows one to obtain new conservation laws even for equations without Lagrangians. A formal Lagrangian can be introduced and used for computing nonlocal conservation laws. For self-adjoint or quasi-self-adjoint equations nonlocal conservation laws can be transformed into local conservation laws. One of the fields of applications of this approach is electromagnetic theory, namely, nonlocal conservation laws are obtained for the generalized Maxwell-Dirac equations. The theory is also applied to the nonlinear magma equation and its nonlocal conservation laws are computed.

No description available.

No description available.

No description available.