Dans un premier tome, l'auteur étudie les impacts géologiques, hydrogéologiques et techniques de l'arrêt de l'exhaure minière dans les dernières concessions exploitées dans le bassin houiller de la Loire. Cette étude permet de définir une stratégie de surveillance de la remonté du niveau d'eau et de déterminer le taux d'agressivité des eaux d'émergences vis-à-vis des bétons. Plusieurs solutions techniques sont aussi envisagées en cas d'apparition d'émergences dans l'ensemble des concessions. Le deuxième tome de cet ouvrage est consacré à la modélisation assistée par ordinateur du milieu. L'auteur a mis au point un logiciel utilisant la méthode des différences finies, pour une modélisation en coupe du milieu qui permet de simuler la remontée du niveau de l'eau. L'outil mis au point peut être transposé et adapté à tout site du même type.

Dans la présente étude, la conductivité thermique et le module d'élasticité de revêtementsd'YPSZ élaborés par projection plasma ont été prédits par modélisations numériques 2D et3D de type différences finies et éléments finis.L'influence de la résolution d'image, de la taille et de la valeur du seuil sur les propriétésprédites du revêtement a été étudiée. En outre, les effets de la méthode numérique et du typede condition aux limites ont été étudiés. En particulier, la quantification de l'effet Knudsen(effet de raréfaction) sur le transfert de chaleur à travers une structure poreuse a été réaliséepar modélisation numérique en combinaison avec l'analyse d'image. Les conductivitéseffectives obtenues par modélisation 3D s'avèrent plus élevées que celles obtenues en 2D, etaussi en meilleur accord avec les résultats mesurés. Une corrélation 2D/3D a été trouvéepour la modélisation de la conductivité thermique : cette corrélation permet de prédire lesvaleurs 3D à partir des valeurs calculées en 2D.

Ce travail de thèse concerne l'étude des écoulements de fluides gazeux compressibles laminaires subsonique-supersonique dans une tuyère de type convergent-divergent. Les écoulements étudiés sont à nombres de Reynolds modérés et s'affranchissent de l'hypothèse de condition adiabatique de paroi couramment utilisée afin de mieux prendre en compte les phénomènes de transfert de chaleur par convection et rayonnement avec le milieu extérieur. Cette étude des phénomènes de transferts conjugués a permis de déterminer le comportement dynamique simultané du fluide et de la paroi de la tuyère. Enfin, compte tenu des niveaux élevés de températures mis en jeu dans ces systèmes, une analyse concernant le comportement thermomécanique de l'ensemble de la structure de paroi avec des matériaux monocouches et multicouches a été réalisé. De nombreuses configurations géométriques, propriétés physiques et conditions aux limites sur le fluide et la paroi ont été analysées. Les résultats présentés montrent, la structure des écoulements à travers les iso-contours de vitesses, des nombres de Mach, des pressions dans le fluide, des températures dans le fluide et dans la paroi ainsi que les déformations et les contraintes de la paroi qui résultent des couplages thermomécaniques. Une analyse des performances de la tuyère, en termes de force de poussée et de coefficient de débit spécifique, est largement discutée dans ce travail.

Les temps de trajet, amplitudes et angles à la source des ondes sismiques sont utilisés dans de nombreuses applications telles que la migration, la tomographie, l'estimation de la sensibilité de détection et la localisation des microséismes. Dans le contexte de la microsismicité, il est nécessaire de calculer en quasi temps réel ces attributs avec précision. Nous avons développé ici un ensemble d'algorithmes rapides et précis en 3D pour des modèles à fort contraste de vitesse.Nous présentons une nouvelle méthode pour calculer les temps de trajet, les amplitudes et les angles à la source des ondes correspondant aux premières arrivées. Plus précisément, nous résolvons l'équation Eikonal, l'équation de transport et l'équation des angles en nous basant sur une approche par différences finies pour des modèles de vitesse en 3D. Nous proposons une nouvelle méthode hybride qui bénéficie des avantages respectifs de plusieurs approches existantes de résolution de l'équation Eikonal. En particulier, les approches classiques proposent généralement de résoudre directement les équations et font l'approximation localement d'une onde plane. Cette approximation n'est pas bien adaptée au voisinage de la source car la courbure du front d'onde est importante. Des erreurs de temps de trajet sont alors générées près de la position de la source, puis propagées à travers tout le modèle de vitesse. Ceci empêche de calculer correctement les amplitudes et les angles à la source puisqu'ils reposent sur les gradients des temps. Nous surmontons cette difficulté en introduisant les opérateurs sphériques ; plus précisément nous reformulons les temps de trajet, amplitudes et angles à la source par la méthode des perturbations.Nous validons nos nouvelles méthodes pour différents modèles à fort contraste de vitesse en 2D et 3D et montrons notre contribution par rapport aux approches existantes. Nos résultats sont similaires à ceux calculés en utilisant la modélisation de la forme d'onde totale alors qu'ils sont bien moins coûteux en temps de calcul. Ces résultats ouvrent donc de nouvelles perspectives pour de nombreuses applications telles que la migration, l'estimation de la sensibilité de détection et l'inversion des mécanismes au foyer.

Les travaux présentés dans ce mémoire s'inscrivent dans la thématique de la transmission d'énergie sans fil, appliquée à l'alimentation à distance de capteurs, de réseaux de capteurs et d'actionneurs à faible consommation. Cette étude porte sur la conception,l'optimisation, la réalisation et la mesure de circuits Rectennas (Rectifying antennas)compacts, à faible coût et à haut rendement de conversion RF-DC.Un outil d'analyse globale, basé sur la méthode des Différences Finies dans le Domaine Temporel (FDTD), a été développé et utilisé pour prédire avec précision la sortie DC des rectennas étudiées. Les résultats numériques obtenus se sont avérés plus précis et plus complets que ceux de simulations à base d'outils commerciaux. La diode Schottky a été rigoureusement modélisée, en tenant compte de ses éléments parasites et de son boîtier SOT23, et introduite dans le calcul itératif FDTD.Trois rectennas innovantes, en technologie micro-ruban, ont été développées,optimisées et caractérisées expérimentalement. Elles fonctionnent à 2.45 GHz et elles ne contiennent ni filtre d'entrée HF ni vias de retour à la masse. Des rendements supérieurs à 80% ont pu être mesurés avec une densité surfacique de puissance de l'ordre de 0.21 mW/cm²(E = 28 V/m). Une tension DC de 3.1 V a été mesurée aux bornes d'une charge optimale de1.05 k_, lorsque le niveau du champ électrique est égal à 34 V/m (0.31 mW/cm²).Des réseaux de rectennas connectées en série et en parallèle ont été développés. Les tensions et les puissances DC ont été doublées et quadruplées à l'aide de deux et de quatre éléments, respectivement.

Ce travail porte sur le problème des partitions minimales, à l'interface entre théorie spectrale et optimisation de forme. Une introduction générale précise le problème et présente des résultats, principalement dûs à B. Helffer, T. Hoffmann-Ostenhof et S. Terracini, qui sont utilisés dans le reste de la thèse.Le premier chapitre est une étude spectrale asymptotique du laplacien de Dirichlet sur une famille de domaines en dimension deux qui tend vers un segment. L'objectif est d'obtenir une localisation des lignes nodales dans la limite des domaines minces. En appliquant les résultats de Helffer, Hoffmann-Ostenhof et Terracini, on montre ainsi que les domaines nodaux des premières fonctions propres forment des partitions minimales.Le deuxième chapitre étudie les valeurs propres de certains opérateurs de Schrödinger sur un domaine plan avec condition au bord de Dirichlet. On considère des opérateurs qui ont un potentiel électrique nul et un potentiel magnétique d'un type particulier, dit d'Aharonov-Bohm, avec des singularités en un nombre fini de points appelés pôles. On démontre que les valeurs propres dépendent continuement des pôles. Dans le cas de pôles distincts et éloignés du bord, on prouve que cette dépendance est analytique lorsque la valeur propre est simple. On exprime de plus une condition suffisante pour que la fonction qui aux pôles associe une valeur propre présente un point critique. On utilise alors la caractérisation magnétique des partitions minimales pour montrer que l'énergie minimale est une valeur critique d'une de ces fonctions.Le troisième chapitre est un article écrit en collaboration avec Virginie Bonnaillie-Noël. Il porte sur une famille d'exemples, les secteurs angulaires de rayon unité et d'ouverture variable, dont on tente de déterminer les partitions minimales. On applique pour cela les théorèmes généraux rappelés dans l'introduction afin de déterminer les partitions nodales qui sont minimales. On s'intéresse plus particulièrement aux partitions minimales en trois domaines. En appliquant les idées du deuxième chapitre, on montre que pour certaines valeur de l'angle, il n'existe aucune partition minimale qui soit symétrique par rapport à la bissectrice du domaine. D'un point de vue quantitatif, on obtient des encadrements précis de l'énergie minimale.Le quatrième chapitre consiste en l'étude des partitions minimales de tores plats dont on fait varier le rapport entre longueur et largeur. On utilise une méthode numérique très différente de celle du troisième chapitre, basée sur un article de B. Bourdin, D. Bucur et É. Oudet. Elle consiste en une relaxation suivie d'une optimisation par un algorithme de gradient projeté. On peut ainsi tester des résultats théoriques antérieurs. Les résultats présentés suggèrent de plus la construction explicite de familles de partitions (en liaison avec des pavages du tore) qui donnent une nouvelle majoration de l'énergie minimale.Un dernier chapitre de perspectives présente plusieurs applications possibles des méthodes décrites dans la thèse.

La ferrite résiduelle δ est présente dans les microstructures de coulée des aciers inoxydables austénitiques. Elle résulte de la transformation incomplète δ→γ ayant lieu l'étape de solidification. Sa présence peut nuire à la forgeabilité à chaud des aciers inoxydables et peut conduire à la formation de criques de rives et de pailles en J lors du laminage à chaud des brames. Ce travail de thèse a pour but de comprendre les mécanismes de la transformation δ→γ à haute température dans les aciers inoxydables austénitiques via une modélisation expérimentale et numérique. La transformation a été étudié dans un alliage ternaire Fe-Cr-Ni coulé par lingot et de composition proche de celle des alliages industriels. Trois morphologies de ferrite ont été mises en évidence à l'état brut de solidification: lattes au bord du lingot, vermiculaire et lattes au centre. Leur cinétique de dissolution est étudiée à des températures allant de 1140°C à 1340°C et caractérisée en termes de fraction de ferrite et profils de composition du Cr et du Ni. La dissolution de la ferrite vermiculaire comprend trois étapes : une croissance initiale transitoire suivie par deux régimes de dissolution à haute puis à faible taux de transformation. D'un autre côté, il a été possible d'étudier la dissolution de la ferrite dans des microstructures multicouches élaborées par l'empilement de plaques de ferrite et d'austénite du système Fe-Cr-Ni et soudées à l'état solide par Compression Isostatique à Chaud puis réduits en épaisseurs par laminages successifs. L'étude et la caractérisation de la cinétique de dissolution de la ferrite est plus facile dans ces microstructures étant donnée la planéité initiale des interfaces δ/γ. L'analyse des résultats expérimentaux a été menée via le développement d'un modèle numérique, à interface mobile, de la transformation de phases δ→γ pilotée par la diffusion. La diffusion peut être traitée dans les géométries plane, cylindrique et sphérique. En guise de validation, le modèle a été utilisé pour analyser la dissolution de la ferrite dans les microstructures multicouches. Par la suite il a été appliqué au cas de la ferrite vermiculaire en usant d'une approche novatrice où la morphologie des dendrites est approximée par une combinaison de cylindres et de sphères. Malgré la simplicité des hypothèse sous-jacentes, le modèle a permis d'expliquer les mécanismes de croissance initiale et de changement de régime de dissolution. D'autre part, via une étude paramétrique, l'effet des données d'entrée a été étudié et les plus pertinentes d'entre eux en termes de prédiction quantitative ont été mises en avant, en particulier la description thermodynamique du digramme Fe-Cr-Ni, le gradient initial et la distribution des rayons des particules de ferrite. / Residual δ-ferrite is widely encountered in the as-cast microstructure of austenitic stainless steels. It stems from the incomplete high temperature solid-state δ→γ transformation occurring upon the solidification stage. Its presence has a detrimental effect the hot workability of stainless steels, leading to the formation of edge cracks and sliver defects during slabs hot rolling. This PhD work aims at bringing more understanding of the kinetics of high temperature δ→γ transformation in austenitic stainless steels via experimental and numerical modeling. The transformation was studied in a ternary Fe-Cr-Ni ingot-cast alloy with composition close to the industrial alloys. Three ferrite morphologies were identified: lathy at the edge of the ingot, vermicular and lathy at the center. Their dissolution kinetics were established at temperatures ranging from 1140°C to 1340°C and characterized in terms of ferrite fraction and Cr and Ni diffusion. The vermicular ferrite undergoes a transient growth followed by a high then a low rate dissolution regimes. On the other hand, ferrite dissolution was also studied in the multilayered microstructures. such microstructures were elaborated by alternating ferrite and austenite sheets of the Fe-Cr-Ni system, diffusion-bonded by Hot isostatic Pressing and reduced in thickness by successive rollings. Dissolution is easier to handle in such microstructures thanks to the initial planar δ/γ interfaces. Analysis of the experimental results were carried out with a numerical moving-boundary model of diffusion-controlled δ→γ transformation. Diffusion can be treated in the planar, cylindrical and spherical geometries. As a preliminary validation, the model was used to analyze kinetics of ferrite dissolution in the multilayered microstructures. It was then applied to the cast alloy using an original descriptive approach combining spheres and cylinders as equivalent morphology of dendritic ferrite. Although based on simplifying assumptions, the model was able to reproduce experimental results with satisfactory agreement. Mechanisms underlying the initial growth of vermicular ferrite and the transition in dissolution regimes were outlined. The effect of a wide range of input parameters has been considered and relevant parameters for quantitative calculations were brought to light, such as thermodynamical descriptions of the Fe-Cr-Ni system, composition gradients and distribution of ferrite's radii.

La transduction électrostatique est utilisée depuis plusieurs décennies dans les fréquences du domaine audible, principalement sous la forme de microphones membranaires. La transposition du même principe de transduction, mais dans un domaine de fréquence au-dessus du MHz, et par l’utilisation de dispositifs micro-usinés, c'est-à-dire produits à l’aide de technologies de photolithographie, a été proposée à partir de la fin des années 1990. Ces transducteurs, désignés sous l’acronyme cMUT (capacitive micromachine ultrasonic transducers), se composent d’un assemblage de transducteurs élémentaires, chacun possédant une partie mobile conventionnellement appelée diaphragme ou membrane, actionnée par la pression électrostatique. Cette thèse s’inscrit dans le développement de transducteurs de ce type destinés au domaine de l’imagerie médicale ultrasonore. Ce secteur d’application utilise actuellement des transducteurs basés sur des céramiques (ou, dans certains cas précis, des polymères) piézoélectriques. Le cMUT est intéressant dans certains sous-domaines d’application des ultrasons médicaux en raison de sa bonne adaptation à une production en grande série, de son intégration plus facile avec des éléments électroniques, de son faible échauffement et de l’absence de matériaux toxiques dans son processus de fabrication. La partie théorique de cette thèse repose sur une approche de modélisation par différences finies. Un modèle basé sur la théorie des plaques minces est développé pour prendre en compte la mécanique du transducteur élémentaire cMUT (c'est-à-dire d’un seul diaphragme). Ce modèle est ensuite complété par l’intégration de l’effet d’un chargement acoustique par un fluide. De façon à modéliser un transducteur entier, il est nécessaire de prendre en compte le couplage acoustique existant entre les différentes membranes. Pour rendre cela possible, un circuit équivalent, permettant de réduire chaque membrane à un système à un seul degré de liberté, est développé. Il est validé en le comparant au modèle de différences finies dans des cas où celui-ci peut être utilisé. Les travaux expérimentaux ont fait appel principalement aux deux techniques de caractérisation suivantes : les mesures d’impédance électrique, et les mesures de déplacement effectuées par interférométrie laser. Ces mesures ont été utilisées dans une double optique. D’une part, dans un objectif de caractérisation, ils ont permis de vérifier la fonctionnalité des dispositifs fabriqués et d’évaluer leurs performances. D’autre part, en comparant différentes configurations entre elles, ils ont rendu possible une validation expérimentale du modèle qui a été mis au point. / Electro-acoustic transduction based on electrostatic force has a long history in the range of audible frequencies, mainly as membrane-based microphones. Starting in the late 1990’s, it has been proposed to use the same principle in the multi-MHz frequency domain, thanks to micro-machined devices – meaning they’re produced through lithography technology. Such transducers, known as cMUT for capacitive micromachine ultrasonic transducers, are made of an assembly of elementary vertically mobile cells, usually designated as membranes, driven by electrostatic force. This PhD work is part of the development of this kind of transducers designed for medical imaging applications. This area currently uses transducers based on piezoelectric ceramics (or piezoelectric polymers for some peculiar cases). CMUT is an interesting alternative for some subdomains of medical ultrasound applications, due to its volume production ability, its easier integration with electronic elements, its low heat dissipation and the absence of toxic materials.

Les discrétisations adaptatives sont importantes dans les problèmes de fluxcompressible/incompressible puisqu'il est souvent nécessaire de résoudre desdétails sur plusieurs niveaux, en permettant de modéliser de grandes régionsd'espace en utilisant un nombre réduit de degrés de liberté (et en réduisant letemps de calcul).Il existe une grande variété de méthodes de discrétisation adaptative, maisles grilles cartésiennes sont les plus efficaces, grâce à leurs stencilsnumériques simples et précis et à leurs performances parallèles supérieures.Et telles performance et simplicité sont généralement obtenues en appliquant unschéma de différences finies pour la résolution des problèmes, mais cetteapproche de discrétisation ne présente pas, au contraire, un chemin faciled'adaptation.Dans un schéma de volumes finis, en revanche, nous pouvons incorporer différentstypes de maillages, plus appropriées aux raffinements adaptatifs, en augmentantla complexité sur les stencils et en obtenant une plus grande flexibilité.L'opérateur de Laplace est un élément essentiel des équations de Navier-Stokes,un modèle qui gouverne les écoulements de fluides, mais il se produit égalementdans des équations différentielles qui décrivent de nombreux autres phénomènesphysiques, tels que les potentiels électriques et gravitationnels. Il s'agitdonc d'un opérateur différentiel très important, et toutes les études qui ontété effectuées sur celui-ci, prouvent sa pertinence.Dans ce travail seront présentés des approches de différences finies et devolumes finis 2D pour résoudre l'opérateur laplacien, en appliquant des patchsde grilles superposées où un niveau plus fin est nécessaire, en laissant desmaillages plus grossiers dans le reste du domaine de calcul.Ces grilles superposées auront des formes quadrilatérales génériques.Plus précisément, les sujets abordés seront les suivants:1) introduction à la méthode des différences finies, méthode des volumes finis,partitionnement des domaines, approximation de la solution;2) récapitulatif des différents types de maillages pour représenter de façondiscrète la géométrie impliquée dans un problème, avec un focussur la structure de données octree, présentant PABLO et PABLitO. Le premier estune bibliothèque externe utilisée pour gérer la création de chaque grille,l'équilibrage de charge et les communications internes, tandis que la secondeest l'API Python de cette bibliothèque, écrite ad hoc pour le projet en cours;3) la présentation de l'algorithme utilisé pour communiquer les données entreles maillages (en ignorant chacune l'existence de l'autre) en utilisant lesintercommunicateurs MPI et la clarification de l'approche monolithique appliquéeà la construction finale de la matrice pour résoudre le système, en tenantcompte des blocs diagonaux, de restriction et de prolongement;4) la présentation de certains résultats; conclusions, références.Il est important de souligner que tout est fait sous Python comme framework deprogrammation, en utilisant Cython pour l'écriture de PABLitO, MPI4Py pour lescommunications entre grilles, PETSc4py pour les parties assemblage et résolutiondu système d'inconnues, NumPy pour les objets à mémoire continue.Le choix de ce langage de programmation a été fait car Python, facile àapprendre et à comprendre, est aujourd'hui un concurrent significatif pourl'informatique numérique et l'écosystème HPC, grâce à son style épuré, sespackages, ses compilateurs et pourquoi pas ses versions optimisées pour desarchitectures spécifiques. / Adaptive discretizations are important in compressible/incompressible flow problems since it is often necessary to resolve details on multiple levels,allowing large regions of space to be modeled using a reduced number of degrees of freedom (reducing the computational time).There are a wide variety of methods for adaptively discretizing space, but Cartesian grids have often outperformed them even at high resolutions due totheir simple and accurate numerical stencils and their superior parallel performances.Such performance and simplicity are in general obtained applying afinite-difference scheme for the resolution of the problems involved, but this discretization approach does not present, by contrast, an easy adapting path.In a finite-volume scheme, instead, we can incorporate different types of grids,more suitable for adaptive refinements, increasing the complexity on thestencils and getting a greater flexibility.The Laplace operator is an essential building block of the Navier-Stokes equations, a model that governs fluid flows, but it occurs also in differential equations that describe many other physical phenomena, such as electric and gravitational potentials, and quantum mechanics. So, it is a very importantdifferential operator, and all the studies carried out on it, prove itsrelevance.In this work will be presented 2D finite-difference and finite-volume approaches to solve the Laplacian operator, applying patches of overlapping grids where amore fined level is needed, leaving coarser meshes in the rest of the computational domain.These overlapping grids will have generic quadrilateral shapes.Specifically, the topics covered will be:1) introduction to the finite difference method, finite volume method, domainpartitioning, solution approximation;2) overview of different types of meshes to represent in a discrete way thegeometry involved in a problem, with a focuson the octree data structure, presenting PABLO and PABLitO. The first one is anexternal library used to manage each single grid’s creation, load balancing and internal communications, while the second one is the Python API ofthat library written ad hoc for the current project;3) presentation of the algorithm used to communicate data between meshes (beingall of them unaware of each other’s existence) using MPI inter-communicators and clarification of the monolithic approach applied building the finalmatrix for the system to solve, taking into account diagonal, restriction and prolongation blocks;4) presentation of some results; conclusions, references.It is important to underline that everything is done under Python as programmingframework, using Cython for the writing of PABLitO, MPI4Py for the communications between grids, PETSc4py for the assembling and resolution partsof the system of unknowns, NumPy for contiguous memory buffer objects.The choice of this programming language has been made because Python, easy to learn and understand, is today a significant contender for the numerical computing and HPC ecosystem, thanks to its clean style, its packages, its compilers and, why not, its specific architecture optimized versions.

Cette thèse se penche sur la propagation d’ondes au sein du système coupléTerre-océan-atmosphère. La compréhension de ces phénomènes a une importance majeure pour l’étude de perturbations sismiques et d’explosions atmosphériques notamment dans le cadre de missions spatiales planétaires. Les formes d’ondes issues du couplage fluide-solide permettent d’obtenir de précieuses informations sur la source du signal ou les propriétés des milieux de propagation. On développe donc deux outils numériques d’ordre élevé pour la propagation d’ondes acoustiques et de gravité. L'u en différences finies et se concentre sur le milieu atmosphérique et la propagation d’ondes linéaires dans un milieu stratifié visqueux et avec du vent. Cette méthode linéaire est validée par des solutions quasi-analytiques reposant sur les équations de dispersion dans une atmosphère stratifiée. Elle est aussi appliquée à deux cas d’études : la propagation d’ondes liée à l’impact d’une météorite à la surface de Mars (mission NASA INSIGHT), et la propagation d’ondes atmosphériques liées au tsunami de Sumatra en 2004. La seconde méthode résout la propagation non-linéaire d’ondes gravito-acoustiques dans une atmosphère couplée, avec topographie, à la propagation d’ondes élastiques dans un solide visco-élastique. Cette méthode repose sur sur le couplage d’une formulation en éléments finis discontinus, pour résoudre les équations de Navier-Stokes la partie fluide, par éléments finis continus pour résoudre les équations de l’élastodynamique dans la partie solide. Elle a été validée grâce à des solutions analytiques ainsi que par des comparaisons avec les résultats de la méthode par différences finies. / This thesis deals with the wave propagation problem within the Earth-ocean-atmosphere coupled system. A good understanding of the these phenomena has a major importance for seismic and atmospheric explosion studies, especially for planetary missions. Atmospheric wave-forms generated by explosions or surface oscillations can bring valuable information about the source mechanism or the properties of the various propagation media. We develop two new numerical full-wave high-order modeling tools to model the propagation of acoustic and gravity waves in realistic atmospheres. The first one relies on a high-order staggered finite difference method and focus only on the atmosphere. It enables the simultaneous propagation of linear acoustic and gravity waves in stratified viscous and windy atmosphere. This method is validated against quasi-analytical solutions based on the dispersion equations for a stratified atmosphere. It has also been employed to investigate two cases : the atmospheric propagation generated by a meteor impact on Mars for the INSIGHT NASA mission and for the study of tsunami-induced acoutic and gravity waves following the 2004 Sumatra tsunami. The second numerical method resolves the non-linear acoustic and gravity wave propagation in a realistic atmosphere coupled, with topography, to the elastic wave propagation in a visco-elastic solid. This numerical tool relies on a discontinuous Galerkin method to solve the full Navier-Stokes equations in the fluid domain and a continuous Galerkin method to solve the elastodynamics equations in the solid domain. It is validated against analytical solutions and numerical results provided by the finite-difference method.