La partie 1 releve de l'Analyse Numerique. Partant de l'interpretation Element Finis Mixtes des schemas volumes finis classiques, l'estimation a posteriori de l'erreur est analysee dans la hierarchie des elements de Raviart-Thomas. Un estimateur calculable est explicite pour ces schemas volumes finis.<br />La partie 2 introduit, d'abord un maillage rectangulaire, puis un maillage structure, une famille de schemas volumes finis de type differences finies. Des essais numeriques sur des problemes modeles montrent que l'ordre prevu par l'analyse peut etre atteint.<br />La partie 3 presente l'application de ces schemas volumes finis a la simulation numerique du comportement d'un bloc de gomme en presence d'une fissure finie. Il s'agit d'un materiau hyperelastique compressible en grandes deformations et differents tenseurs de contraintes, avec tests en quasi-incompressible et des simulations d'endommagement.

Cette thèse propose de nouveaux outils de calcul pour la simulation du moulage de pièces minces par injection de thermoplastiques. Ces travaux s'intègrent au sein d'un code de calcul 3D dédié à la simulation du remplissage de pièces volumiques. La description thermo-mécanique du remplissage du moule est obtenue par une discrétisation éléments finis eulérienne. Nous proposons en première partie une méthode de génération de maillages pour traiter les géométries minces. La solution proposée, issue des méthodes d'optimisations locales, nécessite l'introduction de la notion de métriques. Elle permet de générer et d'adapter des maillages anisotropes bi- ou tridimensionnels non structurés tétraédriques ou triangulaires. En deuxième partie, on a été amené à caractériser l'erreur commise sur le calcul d'un écoulement visqueux incompressible sur de tels maillages. Un nouvel estimateur d'erreur a posteriori pour le problème de Stokes est proposé pour traiter ces maillages. Enfin, l'utilisation de tels maillages peut dégrader, localement, le calcul du remplissage au niveau des surfaces libres. Pour maîtriser cet effet, une méthode d'adaptation de maillage par déformations locales est proposée. S'inspirant des méthodes de régularisation par barycentrage, elle permet d'adapter dynamiquement les mailles au niveau des fronts en mouvement, tout en régularisant le maillage en amont et en aval de ceux-ci. L'ensemble de ces outils permet d'élargir le domaine d'applications du logiciel de simulation existant aux géométries minces, mais aussi d'en améliorer la précision dans le cadre de l'injection multifluides. Des applications, regroupées dans le dernier chapitre, permettent d'illustrer les avancées permises par cette étude.

La théorie de la supraconductivité, modélisée par Ginzburg et Landau, motive les travaux relatifs à l'opérateur de Schrödinger avec champ magnétique. L'objet de cette thèse est d'analyser l'influence de la géométrie du domaine sur l'apparition de la supraconductivité en étendant les résultats existant pour des domaines réguliers à des domaines à coins. L'analyse semi-classique conduit à étudier trois opérateurs modèles : la réalisation de Neumann de l'opérateur de Schrödinger avec champ magnétique constant sur le plan, le demi-plan et les secteurs angulaires. L'étude des deux premiers est bien connue et nous nous concentrons sur le dernier. Après avoir déterminé le bas du spectre essentiel, nous montrons que le bas du spectre est une valeur propre pour un secteur d'angle aigu. Nous explicitons le développement limité de la plus petite valeur propre quand l'angle du secteur tend vers 0 et précisons la localisation de l'état fondamental grâce aux techniques d'Agmon. Nous illustrons et estimons ensuite le comportement des vecteurs et valeurs propres à l'aide d'outils numériques basés sur la méthode des éléments finis. La localisation de l'état fondamental rend le problème discret très mal conditionné mais l'analyse des propriétés de l'opérateur et des défauts des méthodes classiques permet malgré tout de mettre en oeuvre un algorithme robuste et efficace calculant l'état fondamental. Afin d'améliorer les résultats numériques, nous construisons des estimateurs a posteriori pour ce problème aux valeurs propres et utilisons les techniques de raffinement de maillages pour localiser l'état propre dans des domaines généraux et étudier la variation du bas du spectre en fonction de l'angle du secteur.

En géosciences, les applications associées au calage de modèles d'écoulement nécessitent d'appeler plusieurs fois un simulateur au cours d'un processus d'optimisation. Or, une seule simulation peut durer plusieurs heures et l'exécution d'une boucle complète de calage peut s'étendre sur plusieurs jours. Diminuer le temps de calcul global à l'aide des techniques de bases réduites (RB) constitue l’objectif de la thèse.Il s'agit plus précisément dans ce travail d'appliquer ces techniques aux écoulements incompressibles diphasiques eau-huile en milieu poreux. Ce modèle, bien que simplifié par rapport aux modèles utilisés dans l'industrie pétrolière, constitue déjà un défi du point de vue de la pertinence de la méthode RB du fait du couplage entre les différentes équations, de la forte hétérogénéité des données physiques, ainsi que du choix des schémas numériques de référence.Nous présentons d'abord le modèle considéré, le schéma volumes finis (VF) retenu pour l'approximation numérique, ainsi que différentes paramétrisations pertinentes en simulation de réservoir. Ensuite, après un bref rappel de la méthode RB, nous mettons en oeuvre la réduction du problème en pression à un instant donné en suivant deux démarches distinctes. La première consiste à interpréter la discrétisation VF comme une approximation de Ritz-Galerkine, ce qui permet de se ramener au cadre standard de la méthode RB mais n'est possible que sous certaines hypothèses restrictives. La seconde démarche lève ces restrictions en construisant le modèle réduit directement au niveau discret.Enfin, nous testons deux stratégies de réduction pour la collection en temps de pressions paramétrées par les variations de la saturation. La première considère le temps juste comme un paramètre supplémentaire. La seconde tente de mieux capturer la causalité temporelle en introduisant les trajectoires en temps paramétrées. / In geosciences, applications involving model calibration require a simulator to be called several times with an optimization process. However, a single simulation can take several hours and a complete calibration loop can extend over serval days. The objective of this thesis is to reduce the overall simulation time using reduced basis (RB) techniques.More specifically, this work is devoted to applying such techniques to incompressible two-phase water-oil flows in porous media. Despite its relative simplicity in comparison to other models used in the petroleum industry, this model is already a challenge from the standpoint of reduced order modeling. This is due to the coupling between its equations, the highly heterogeneous physical data, as well as the choice of reference numerical schemes.We first present the two-phase flow model, along with the finite volume (FV) scheme used for the discretization and relevant parameterizations in reservoir simulation. Then, after having recalled the RB method, we perform a reduction of the pressure equation at a fixed time step by two different approaches. In the first approach, we interpret the FV discretization as a Ritz-Galerkine approximation, which takes us back to the standard RB framework but which is possible only under severe assumptions. The second approach frees us of these restrictions by building the RB method directly at the discrete level.Finally, we deploy two strategies for reducing the collection in time of pressuresparameterized by the variations of the saturation. The first one simply considers time as an additional parameter. The second one attempts to better capture temporalcausality by introducing parameterized time-trajectories.

Le but de cette étude est d'améliorer les performances, en termes d'espace mémoire et de temps de calcul, des simulations actuelles de l'Interaction mécanique Pastille-Gaine (IPG), phénomène complexe pouvant avoir lieu lors de fortes montées en puissance dans les réacteurs à eau sous pression. Parmi les méthodes de raffinement de maillage, méthodes permettant de simuler efficacement des singularités locales, une approche multi-grille locale a été choisie car elle présente l'intérêt de pouvoir utiliser le solveur en boîte noire tout en ayant un faible nombre de degrés de liberté à traiter par niveau. La méthode Local Defect Correction (LDC), adaptée à une discrétisation de type éléments finis, a tout d'abord été analysée et vérifiée en élasticité linéaire, sur des configurations issues de l'IPG, car son utilisation en mécanique des solides est peu répandue. Différentes stratégies concernant la mise en oeuvre pratique de l'algorithme multi-niveaux ont également été comparées. La combinaison de la méthode LDC et de l'estimateur d'erreur a posteriori de Zienkiewicz-Zhu, permettant d'automatiser la détection des zones à raffiner, a ensuite été testée. Les performances obtenues sur des cas bidimensionnels et tridimensionnels sont très satisfaisantes, l'algorithme proposé se montrant plus performant que des méthodes de raffinement h-adaptatives. Enfin, l'algorithme a été étendu à des problèmes mécaniques non linéaires. Les questions d'un raffinement espace/temps mais aussi de la transmission des conditions initiales lors du remaillage ont entre autres été abordées. Les premiers résultats obtenus sont encourageants et démontrent l'intérêt de la méthode LDC pour des calculs d'IPG. / The aim of this study is to improve the performances, in terms of memory space and computational time, of the current modelling of the Pellet-Cladding mechanical Interaction (PCI),complex phenomenon which may occurs during high power rises in pressurised water reactors. Among the mesh refinement methods - methods dedicated to efficiently treat local singularities - a local multi-grid approach was selected because it enables the use of a black-box solver while dealing few degrees of freedom at each level. The Local Defect Correction (LDC) method, well suited to a finite element discretisation, was first analysed and checked in linear elasticity, on configurations resulting from the PCI, since its use in solid mechanics is little widespread. Various strategies concerning the implementation of the multilevel algorithm were also compared. Coupling the LDC method with the Zienkiewicz-Zhu a posteriori error estimator in orderto automatically detect the zones to be refined, was then tested. Performances obtained on two-dimensional and three-dimensional cases are very satisfactory, since the algorithm proposed is more efficient than h-adaptive refinement methods. Lastly, the LDC algorithm was extended to nonlinear mechanics. Space/time refinement as well as transmission of the initial conditions during the remeshing step were looked at. The first results obtained are encouraging and show the interest of using the LDC method for PCI modelling.

Cette thèse contribue à la simulation numérique des écoulements d'un fluide compressible modélisé par les équations de Euler et Navier-Stokes: étude d'un schéma d'ordre élévé basé sur une matrice de masse, modélisation des écoulement turbulents compressibles à très haut Reynolds, développement des estimateurs et correcteurs d'erreurs a posteriori et a priori, et adaptation de maillage anisotrope pour les fonctionnelles d'observation. Concernant la prédiction des écoulements turbulents, on s'est intéressé aux modèles hybrides de type RANS/LES comportant les nouveautés suivantes: traitement des tourbillons de grande échelle utilisant la formulation VMS (Variational Multi-Scale) et du RANS employé sur la paroi sur une distance imposée via une zone de protection conçue afin d'éviter le phénomène assez commun apelé "grid induced model depletion". Le niveau de viscosité du modèle VMS-LES est de plus controllé par un procédé de double filtre dynamique. La seconde partie concerne l'adaptation de maillage anisotrope pour mieux observer une fontionnelle d'observation. Les estimations a priori sont réalises pour le modèle des équations d'Euler et Navier-Stokes en instationnaire en 2D et 3D. A partir de ces estimations on sait définir les maillages optimaux au cours du calcul instationnaire, en fonction de l'état et de l'état adjoint. Le système d'optimalité est discrétisé et résolu à l'aide d'une méthode de point fixe instationnaire global, comportant une stratégie de stockage/recalcul pour le couplage état/ état adjoint. Des applications à la propagation d'ondes de choc et d'ondes acoustiques sont présentées.

Le but de cette étude est d'améliorer les performances, en termes d'espace mémoire et de temps de calcul, des simulations actuelles de l'Interaction mécanique Pastille-Gaine (IPG), phénomène complexe pouvant avoir lieu lors de fortes montées en puissance dans les réacteurs à eau sous pression. Parmi les méthodes de raffinement de maillage, méthodes permettant de simuler efficacement des singularités locales, une approche multi-grille locale a été choisie car elle présente l'intérêt de pouvoir utiliser le solveur en boîte noire tout en ayant un faible nombre de degrés de liberté à traiter par niveau. La méthode Local Defect Correction (LDC), adaptée à une discrétisation de type éléments finis, a tout d'abord été analysée et vérifiée en élasticité linéaire, sur des configurations issues de l'IPG, car son utilisation en mécanique des solides est peu répandue. Différentes stratégies concernant la mise en oeuvre pratique de l'algorithme multi-niveaux ont également été comparées. La combinaison de la méthode LDC et de l'estimateur d'erreur a posteriori de Zienkiewicz-Zhu, permettant d'automatiser la détection des zones à raffiner, a ensuite été testée. Les performances obtenues sur des cas bidimensionnels et tridimensionnels sont très satisfaisantes, l'algorithme proposé se montrant plus performant que des méthodes de raffinement h-adaptatives. Enfin, l'algorithme a été étendu à des problèmes mécaniques non linéaires. Les questions d'un raffinement espace/temps mais aussi de la transmission des conditions initiales lors du remaillage ont entre autres été abordées. Les premiers résultats obtenus sont encourageants et démontrent l'intérêt de la méthode LDC pour des calculs d'IPG.