L'objet de cette thèse est l'étude et la mise en œuvre d'une méthode d’éléments finis multi-échelles pour la simulation d'écoulements miscibles en milieux poreux. La définition des fonctions de base multi-échelles suit l'idée introduite par F. Ouaki. La nouveauté de ce travail consiste à combiner cette approche multi-échelle avec des éléments finis de type Galerkine Discontinus (DG) de façon à pouvoir utiliser ces nouveaux éléments sur des maillages non-conformes composés de mailles de formes diverses. Nous rappelons, dans un premier temps, le principe des méthodes DG et montrons comment ces méthodes peuvent être utilisées pour discrétiser une équation de convection-diffusion instationnaire identique à celle rencontrée dans le problème d'écoulement considéré dans ce travail. Après avoir vérifié l'existence et l'unicité d'une solution à ce problème, nous redémontrons la convergence des méthodes DG vers cette solution en établissant une estimation d'erreur a priori. Nous introduisons, ensuite, les éléments finis multi-échelles non conformes et détaillons leur mise en œuvre sur ce problème de convection-diffusion. En supposant les conditions aux limites et les paramètres du problème périodiques, nous montrons une nouvelle estimation d'erreur a priori pour cette méthode. Dans une seconde partie, nous considérons le problème d'écoulement complet où l'équation considérée dans la première partie est résolue de manière couplée avec l'équation de Darcy. Nous introduisons différents cas tests inspirés de modèles d'écoulements rencontrés en géosciences et comparons les solutions obtenues avec les deux méthodes DG, à savoir la méthode classique utilisant un seul maillage et la méthode étudiée ici. Nous proposons de nouvelles conditions aux limites pour la résolution des problèmes de cellule qui permettent, par rapport à des conditions aux limites linéaires plus classiquement utilisées, de mieux reproduire les variations des solutions le long des interfaces du maillage grossier. Les résultats de ces tests montrent que la méthode multi-échelle proposée permet de calculer des solutions proches de celles obtenues avec la méthode DG sur un seul maillage et de réduire, de façon significative, la taille du système linéaire à résoudre à chaque pas de temps. / This work deals with the study and the implementation of a multiscale finite element method for the simulation of miscible flows in porous media. The definition of the multiscale basis functions is based on the idea introduced by F. Ouaki. The novelty of this work lies in the combination of this multiscale approach with Discontinuous Galerkin methods (DG) so that these new finite elements can be used on nonconforming meshes composed of cells with various shapes. We first recall the basics of DG methods and their application to the discretisation of a convection-diffusion equation that arises in the flow problem considered in this work. After establishing the existence and uniqueness of a solution to the continuous problem, we prove again the convergence of DG methods towards this solution by establishing an a priori error estimate. We then introduce the nonconforming multiscale finite element method and explain how it can be implemented for this convection-diffusion problem. Assuming that the boundary conditions and the parameters of the problem are periodic, we prove a new a priori error estimate for this method. In a second part, we consider the whole flow problem where the equation, studied in the first part of that work, is coupled and simultaneously solved with Darcy equation. We introduce various synthetic test cases which are close to flow problems encountered in geosciences and compare the solutions obtained with both DG methods, namely the classical method based on the use of a single mesh and the one studied here. For the resolution of the cell problems, we propose new boundary conditions which, compared to classical linear conditions, allow us to better reproduce the variations of the solutions on the interfaces of the coarse mesh. The results of these tests show that the multiscale method enables us to calculate solutions which are close to the ones obtained withDG methods on a single mesh and also enables us to reduce significantly the size of the linear system that has to be solved at each time step.

Dans cette thèse, nous analysons une méthode de Galerkine discontinue (GD) et deux estimateurs d'erreur a posteriori pour l'équation d'advection-diffusion-réaction linéaire et stationnaire avec diffusion hétérogène. La méthode GD considérée, la méthode SWIP, utilise des moyennes pondérées dont les poids dépendent de la diffusion. L'analyse a priori montre que la convergence est optimale en le pas du maillage et robuste par rapport aux hétérogénéités de la diffusion, ce qui est confirmé par les tests numériques. Les deux estimateurs d'erreur a posteriori sont obtenus par une analyse par résidus et contrôlent la (semi-)norme d'énergie de l'erreur. L'analyse d'efficacité locale montre que presque tous les estimateurs sont indépendants des hétérogénéités. Le deuxième estimateur d'erreur est plus précis que le premier, mais son coût de calcul est légèrement plus élevé. Cet estimateur est basé sur la construction d'un flux H(div)-conforme dans l'espace de Raviart-Thomas-Nédéléc.

Cette thèse vise à étudier le couplage entre méthodes de Galerkine discontinue (DG) et méthodes de différences finies (DF) en maillages hybrides non structuré / cartésien, en vue d'applications en aéroacoustique numérique. L'idée d'une telle approche consiste à pouvoir tirer profit localement des avantages respectifs de ces méthodes, soit, en d'autres termes, à pouvoir prendre en compte la présence de géométries complexes par une méthode DG en maillage non structuré, et les zones qui en sont suffisamment éloignées par une méthode DF en maillage cartésien, moins coûteuse. Plus précisément, il s'agit de concevoir un algorithme d'hybridation de ces deux types de schémas pour l'approximation des équations d'Euler linéarisées, puis d'évaluer avec attention le comportement numérique des solutions qui en sont issues. De par le fait qu'aucun résultat théorique ne semble actuellement atteignable dans un cas général, cette étude est principalement fondée sur une démarche d'expérimentation numérique. Par ailleurs, l'intérêt d'une telle hybridation est illustré par son application à un calcul de propagation acoustique dans un cas réaliste / This thesis aims at studying coupling techniques between Discontinuous Galerkin (DG) and finite difference (FD) schemes in a non-structured / Cartesian hybrid-mesh context,in the framework of Aeroacoustics computations. The idea behind such an approach is the possibility to locally take advantage of the qualities of each method. In other words, the goal is to be able to deal with complex geometries using a DG scheme on a non-structured mesh in their neighborhood, while solving the rest of the domain using a FD scheme on a cartesian grid, in order to alleviate the needs in computational resources. More precisely, this work aims at designing an hybridization algorithm between these two types of numerical schemes, in the framework of the approximation of the solutions of the Linearized Euler Equations. Then, the numerical behaviour of hybrid solutions is cautiously evaluated. Due to the fact that no theoretical result seems achievable at the present time, this study is mainly based on numerical experiments. What's more, the interest of such an hybridization is illustrated by its application to an acoustic propagation computation in a realistic case

Mes travaux de thèse portent sur la construction de conditions aux limites absorbantes (CLAs) pour des problèmes de propagation d'ondes posés dans des milieux limités par des surfaces régulières. Ces conditions sont nouvelles car elles prennent en compte non seulement les ondes proagatives (comme la plupart des CLAs existantes) mais aussi les ondes évanescentes et rampantes. Elles sont donc plus performantes que les conditions existantes. De plus, elles sont facilement implémentables dans un schéma d'éléments finis de type Galerkine Discontinu (DG) et ne modifie pas la condition de stabilité de Courant-Friedrichs-Lewy (CFL). Ces CLAs ont été implémentées dans un code simulant la propagation des ondes acoustiques ainsi que dans un code simulant la propagation des ondes en régime harmonique. Les comparaisons réalisées entre les nouvelles conditions et celles qui sont les plus utilisées dans la littérature montrent que prendre en compte les ondes évanescentes et les ondes rampantes permet de diminuer les réflexions issues de la frontière artificielle et donc de rapprocher la frontière artificielle du bord de l'obstacle. On limite ainsi les coûts de calcul, ce qui est un des avantages de mes travaux. De plus, compte tenu du fait que les nouvelles CLAs sont écrites pour des frontières quelconques, elles permettent de mieux adapter le domaine de calcul à la forme de l'obstacle et permettent ainsi de diminuer encore plus les coûts de calcul numérique. / In my PhD, I have worked on the construction of absorbing boundary conditions (ABCs) designed for wave propagation problems set in domains bounded by regular surfaces. These conditions are new since they take into account not only propagating waves (as most of the existing ABCs) but also evanescent and creeping waves. Therefore, they outperform the existing ABCs. Moreover, they can be easily implemented in a discontinuous Galerkin finite element scheme and they do not change the Courant-Friedrichs-Lewy stability condition. These ABCs have been implemented in two codes that respectively simulate the propagation of acoustic waves and harmonic waves. The comparisons performed between these ABCs and the ABCs mostly used in the litterature show that when we take into account evanescent and creeping waves, we reduce the reflections coming from the artificial boundary. Therefore, thanks to these new ABCs, the artificial boundary can get closer to the obstacle. Consequently, we reduce the computational costs which is one of the advantages of my work. Moreover, since these new ABCs are written for any kind of boundary, we can adapt the shape of the computational domain and thus we can reduce again the computational costs.

L’industrie pétrolière s’intéresse désormais à des régions de la terre de plus en plus difficiles d’accès et il est essentiel de proposer des techniques permettant de garantir l’efficaité d’un forage. Parmi ces techniques, la Reverse Time Migration (RTM) est connue pour sa précision. Elle utilise les ondes réfléchies pour reconstruire une carte du sous-sol représentant les interfaces géophysiques. Elle peut être décrite en trois étapes : (i) propager le champs émis par les sources durant la campagne d’acquisition; (ii) pour chaque source, propager le champ enregistré par les récepteurs; (iii) obtenir une image du sous-sol en appliquant une condition d’imagerie à chaque pas de temps et pour chaque source. Cette technique requiert de très grosses capacités de calcul et il est encore difficile d’imager des milieux réalistes 3D, même avec l’aide du calcul haute performance. Nous avons choisi la méthode de Galerkine discontinue pour modéliser la partie propagation car elle permet d’obtenir des solutions précises et est adaptable au calcul parallèle. La quantité d’information à sauvegarder pour faire une corrélation étant importante, on se doit de trouver un algorithme de calcul d’images du sous-sol réduisant ce coût. Nous avons utilisé l’algorithme de Griewank, appelé “Optimal Checkpointing”. Ce problème de coût étant réglé, on se doit de considérer l’efficacité des ondes élastiques incluant des champs multiples pour améliorer la précision de l’image. La condition traditionnelle de J. Claerbout ne prend pas en compte les conversions d’ondes, et n’est alors surtout utile que dans le cas acoustique. De plus, les ondes P et S interagissant entre elles, il est intéressant de trouver une condition d’imagerie utilisant ce fait. Cela a été abordé dans le cadre de la méthode de l’état adjoint dans les travaux de A. Tarantola et J. Tromp et ce travail en propose utilisation dans le cadre de la RTM. Nous proposons une nouvelle condition d’imagerie prenant en compte les paramètres élastiques du milieu considéré et permettant de supprimer les artefacts numériques. Nous illustrons les images sur des cas industriels / Since a large number of sedimentary basins have been explored, oil exploration is now interested in investigating regions of the Earth which are hostile. Among existing methods for seismic imaging, Reverse Time Migration (RTM) is a technique known by industry to be efficient. The RTM uses reflected waves and is able to construct a map of the subsurface which is depicted by the interfaces limiting the geophysical layers. The algorithm of RTM can be described as a three-step procedure: (i) compute the wavefields emitted by the sources used during the seismic acquisition campaign; (ii) for each source, compute the so-called “backpropagated wavefield”, which is the wavefield obtained by using as sources the signals recorded at the receivers during the acquisition campaign and by reversing the time; (iii) get an image of the subsurface by applying an imaging condition combining the propagated and the backpropagated wavefields at each time step of the numerical scheme and for each source. This technique is computationnaly intensive and it is still difficult to image realistic 3D elastic media, even with the help of HPC. We have thus chosen to consider high-order Discontinuous Galerkin Methods which are known to be well-adapted to provide accurate solutions based upon parallel computing. As we need to correlate a lot of wavefields, we need to find an algorithm reducing the CPU time and the storage : this is the Griewank’s algorithm, so-called “Optimal Checkpointing”. The traditional imaging condition, proposed by J. Claerbout, does not take wave conversions into account and since P-wave and S-wave interact with each other, it might be relevant to use an imaging condition including these interactions. In fact, this has been done successfully by A. Tarantola and J. Tromp for seismology applications based upon the inversion of the global Earth. In this work, we propose a new imaging condition using the elastic parameters which attenuates numerical artifacts. We illustrate the properties of the new imaging condition on industrial benchmarks like the Marmousi model. In particular, we compare the new imaging condition with other imaging conditions by using as criteria the quality of the image.

Dans cette thèse nous proposons de développer un système d'imagerie ultrasonore innovante de micro- défauts basé sur l'utilisation conjointe de techniques d'acoustique non linéaire et du concept de "transducteur à cavité chaotique". Ce transducteur correspond à la combinaison d'une céramique piézoélectrique collée sur une cavité de forme chaotique et du principe de retournement temporel. La faisabilité et les performances de ce nouveau système sont explorées par des simulations numériques. Des paramètres optimaux d'utilisation pour une implémentation expérimentale sont proposés. Une grande partie des travaux menés dans le cadre de cette thèse se concentre sur le développement d'outils numériques permettant l'amélioration de telles techniques d'imagerie. Un schéma d'éléments finis de type Galerkin Discontinu (GD) est étendu à l'élastodynamique non linéaire. Un type de zone absorbante parfaitement adaptée, appelée "Nearly Perfectly Matched Layer" (NPML) a aussi été développé. Dans le cas de matériaux orthotropes, comme des problèmes de stabilité apparaissent, un mélange de NPML et de zone atténuante, dont on contrôle la proportion respective, est introduit afin de stabiliser les NPML. Une validation expérimentale du concept de "transducteur à cavité chaotique" pour la focalisation dans un milieu solide, réverbérant ou non, en utilisant une seule source est réalisée. Les méthodes de retournement temporel et de filtre inverse sont présentées et comparées. La démonstration expérimentale qu'un "transducteur à cavité chaotique" peut être utilisé conjointement avec les méthodes d'inversion d'impulsion afin de réaliser une image de non linéarités localisées est présentée

Mes travaux de thèse portent sur la construction de conditions aux limites absorbantes (CLAs) pour des problèmes de propagation d'ondes posés dans des milieux limités par des surfaces régulières. Ces conditions sont nouvelles car elles prennent en compte non seulement les ondes proagatives (comme la plupart des CLAs existantes) mais aussi les ondes évanescentes et rampantes. Elles sont donc plus performantes que les conditions existantes. De plus, elles sont facilement implémentables dans un schéma d'éléments finis de type Galerkine Discontinu (DG) et ne modifie pas la condition de stabilité de Courant-Friedrichs-Lewy (CFL). Ces CLAs ont été implémentées dans un code simulant la propagation des ondes acoustiques ainsi que dans un code simulant la propagation des ondes en régime harmonique. Les comparaisons réalisées entre les nouvelles conditions et celles qui sont les plus utilisées dans la littérature montrent que prendre en compte les ondes évanescentes et les ondes rampantes permet de diminuer les réflexions issues de la frontière artificielle et donc de rapprocher la frontière artificielle du bord de l'obstacle. On limite ainsi les coûts de calcul, ce qui est un des avantages de mes travaux. De plus, compte tenu du fait que les nouvelles CLAs sont écrites pour des frontières quelconques, elles permettent de mieux adapter le domaine de calcul à la forme de l'obstacle et permettent ainsi de diminuer encore plus les coûts de calcul numérique.

Ce mémoire présente une étude de la modélisation POD-Galerkine réduite dans l'objectif du développement d'algorithmes robustes et efficaces de contrôle actif d'écoulements. Les deux premiers chapitres sont consacrés à la définition et au calcul numérique d'une base de modes POD et à la construction formelle de modèles dynamiques réduits par la méthode de Galerkine à partir des équations de Navier-Stokes. Le cas des écoulements incompressibles est traité en détail. Une étude de la modélisation d'un écoulement tridimensionnel turbulent, pour lequel le problème de la modélisation des petites échelles se pose, est ensuite menée : les interactions entre modes POD et les effets de la réduction de la base modale sont analysés qualitativement et quantitativement, à l'aide notamment d'une paramétrisation visqueuse. Dans le chapitre suivant, des méthodes de calibration, qui<br />reposent sur la résolution d'un problème d'optimisation, sont développées afin de pouvoir calculer automatiquement des modèles<br />réduits fiables pour un coût informatique raisonnable.<br />Enfin, le dernier chapitre est consacré au problème du contrôle d'écoulements par des actionneurs de soufflage ou d'aspiration : les stratégies d'exploitation de la modélisation POD-Galerkine pour le contrôle sont abordées, puis des investigations numériques sont présentées.

Th. doct.--Énerg.--Paris--Éc. centrale, 1996. N°: 1996-45. / 1996 d'après la déclaration de dépôt légal. Bibliogr. p. 215-222. Résumé en français et en anglais.

Dans cette thèse, nous nous intéressons à la modélisation mathématique des hétérogénéités de longueurs caractéristiques beaucoup plus petites que la longueur d'ondes. La thèse consiste en deux parties. La partie théorique est dédiée à l'obtention d'un développement asymptotique raccordé: la solution est décrite à l'aide d'un développement de champ proche au voisinage de l'obstacle et par un développement de champ lointain hors de ce voisinage. Le développement de champ lointain met en jeu des solutions singulières de l'équation des ondes tandis que le champ proche lui est régi par un modèle quasi-statique. Ces deux développements sont alors raccordés dans une zone intermédiaire dite de raccord. Nous obtenons alors des estimations d'erreurs permettant de rendre rigoureux ce développement asymptotique formel. La deuxième partie est numérique. Elle décrit à la fois la méthode de Galerkine discontinue, une méthode de raffinement de maillage espace-temps et propose une discrétisation des modèles asymptotiques obtenues précédemment. Elle est illustrée par un certain nombre de tests numériques. / In this thesis, we focus our attention on the modeling of heterogeneities which are smaller than the wavelength. The document is decomposed into two parts : a theoretical one and a numerical one. In the first part, we derive a matched asymptotic expansion composed of a far-field expansion and a near-field expansion. The terms of the far-field expansion are singular solutions of the wave equation whereas the terms of the near-field expansion satisfy quasistatic problems. These expansions are matched in an intermediate region. We justify mathematically this theory by proving error estimates. In the second part, we describe the Discontinuous Galerkin method, a local time stepping method and the implementation of the matched asymptotic method. Numerical simulations illustrate these results.