Ce travail de thèse concerne la prédiction des structures en béton soumises à des chargements extrêmes. Il s'intéresse plus particulièrement au comportement du béton sous fort confinement où la contrainte peut atteindre des niveaux de l'ordre du giga Pascal. La modélisation de ce comportement doit être capable de reproduire la compaction irréversible. Pour ce faire, deux lois de comportement élasto-plastique - endommageable ont été développées et implantées dans un code aux éléments discrets. Les paramètres utilisés dans ces lois sont calibrés par les simulations des essais de traction/compression uniaxial, des essais hydrostatiques et triaxiaux. Une fois les paramètres calibrés, la loi montrant le meilleur agrément avec l'expérience a été choisie pour la prédiction de la réponse du béton sous différents niveaux de confinement. Les résultats du modèle sont analysés non seulement à l'échelle macroscopique mais également à l'échelle de l'élément discret. La nécessité de prendre en compte une loi d'interaction de type élasto-plastique-endommageable est aussi montrée. La deuxième partie du travail de thèse développe une méthode de couplage entre le modèle éléments discrets et un modèle d'écoulement compressible en tenant compte des mécanismes physiques fondamentaux d'interaction entre l'écoulement interne et les particules solides d'un matériau poreux. Le problème d'écoulement est résolu par une méthode en volumes finis, où le volume est discrétisé en tétraèdres issus d'une triangulation régulière de Delaunay. Notre modèle est une adaptation aux fluides compressibles d'un modèle développé initialement pour les écoulements incompressibles. Ce couplage a été utilisé pour simuler le comportement triaxial des bétons humides et saturés sous différents niveaux de confinement. Les résultats nous montrent une bonne reproduction du comportement non-drainé du béton saturé sous faible confinement. Pour fort confinement, les simulations ne se rapprochent des résultats expérimentaux qu'au prix d'une compressibilité du fluide plus faible que celle de l'eau. Par ailleurs, la contrainte effective était une variable pertinente pour décrire le comportement du béton humide par un état limite intrinsèque indépendant du degré de saturation.

L'ozone (O3) possède une anomalie isotopique en oxygène qui est unique et caractéristique, offrant ainsi un précieux traceur des processus oxydatifs à l'œuvre dans l'atmosphère moderne mais aussi ceux ayant eu lieu dans le passé. Cette signature isotopique, dénotée Δ17O, se propage au sein du cycle atmosphérique de l'azote réactif (NOx = NO + NO2) et est préservée lors du dépôt du nitrate (NO3-) présent dans l'aérosol, par exemple. L'anomalie isotopique en oxygène portée par le nitrate, Δ17O(NO3-), représente ainsi un traceur de l'importance relative de l'ozone ainsi que celle d'autres oxydants dans le cycle des NOx. Ces dernières années, de nombreux travaux de recherche ont été dédiés à l'interprétation des mesures de Δ17O(NO3-). Pourtant, les processus atmosphériques responsables du transfert de l'anomalie isotopique de l'ozone vers le nitrate ainsi que leur influence globale sur la composition isotopique du nitrate à différentes échelles spatiales et temporelles sont encore mal compris. De plus, la magnitude absolue ainsi que la variabilité spatio-temporelle de Δ17O(O3) sont peu contraintes, car il est difficile d'extraire de l'ozone de l'air ambiant. Cet obstacle technique contrecarre l'interprétation des mesures de Δ17O depuis plus d'une décennie. Les questions scientifiques posées au cours de ce travail de thèse ont été choisies dans le but de combler ces lacunes. Le principal outil d'analyse utilisé dans ce travail est la " méthode bactérienne " associée à la spectrométrie de masse en flux continu (CF-IRMS), une combinaison de techniques qui permet l'analyse de la composition isotopique totale du nitrate (c'est-à-dire, la mesure de δ15N, δ18O et Δ17O). Cette méthode a été employée pour l'analyse isotopique d'échantillons de nitrate obtenus pour deux cas d'études : (i) une étude des variations spatiales de la composition isotopique du nitrate atmosphérique sur la côte californienne à l'échelle journalière; et (ii) une étude du transfert du nitrate et de sa composition isotopique à l'interface entre l'air et la neige à l'échelle saisonnière sur le plateau Antarctique. En outre, cette méthode a été adaptée à la caractérisation isotopique de l'ozone via la conversion chimique de ses atomes d'oxygène terminaux en nitrate. Au cours de cette thèse, un important jeu de données rassemblant de nombreuses mesures troposphériques de Δ17O(O3) a été obtenu, incluant une année entière de mesures à Grenoble, France (45 °N) ainsi qu'un transect latitudinal de collecte dans la couche limite au-dessus de l'océan Atlantique, entre 50 °S to 50 °N. Ces observations ont permis de doubler le nombre de mesures troposphériques de Δ17O(O3) existantes avant cette thèse et d'accroître de manière conséquente notre représentation globale de cette variable isotopique essentielle. Enfin, les deux cas étudiés et présentés dans ce document révèlent des aspects nouveaux et inattendus de la dynamique isotopique du nitrate atmosphérique, avec d'importantes conséquences potentielles pour la modélisation de la qualité de l'air et l'interprétation de l'information isotopique contenue dans les carottes de glace prélevées aux pôles.

Dans un contexte d'accroissement des marges de sûreté des réacteurs de générations II et III vis-à-vis des scénarios accidentels, des efforts importants de recherche sont consacrés à l'amélioration de la microstructure du combustible MOX à l'issue de son procédé de fabrication. Les deux caractéristiques microstructurales recherchées sont l'accroissement de l'homogénéité de répartition du plutonium et l'augmentation de la taille de grain. Dans cette optique, une solution envisagée est l'ajout lors du procédé de fabrication et sans modification de celui-ci, de sesquioxyde de chrome Cr2O3. Une précédente thèse sur le sujet a permis de proposer un modèle d'homogénéisation de la répartition du plutonium suite à l'ajout de Cr2O3. L'auteur a souligné l'importance de la formation du précipité PuCrO3 aux joints de grains lorsque la solubilité du chrome dans la matrice (U,Pu)O2 est atteinte. Cependant, les mécanismes d'action du chrome n'ont été étudiés que pour une atmosphère de frittage unique. Plusieurs points restent également à approfondir, notamment la solubilité du chrome et les conditions optimales de formation du précipité PuCrO3. Dans un premier temps, une étude de la spéciation du chrome solubilisé et précipité dans l'oxyde mixte (U,Pu)O2 a été réalisée. Les techniques ayant permis d'analyser directement le chrome sont la microsonde électronique et la spectroscopie d'absorption des rayons X. Il a été montré que le degré d'oxydation et l'environnement du chrome solubilisé sont indépendants de la pression partielle d'oxygène imposée lors du frittage et de la teneur en plutonium de l'oxyde mixte. La nature des précipités et la solubilité du chrome dépendent, quant à eux, de la variable thermodynamique et de la teneur en Pu. Sur la base de ces résultats, un modèle de solubilité du chrome dans l'oxyde mixte (U,Pu)O2-x a été construit. Ce modèle a été réalisé en fonction de la teneur en plutonium y de la solution solide (U1-yPuy)O2-x (y = 0,11 ; 0,275 et 1) et sur la gamme de potentiel d'oxygène d'intérêt pour le frittage du combustible (-445 kJ.mol -1< µO2 < -360 kJ.mol -1). Outre l'optimisation du dopage, ce modèle permet de définir les conditions optimales de formation du précipité PuCrO3 en fonction de la teneur en plutonium et de l'atmosphère de frittage. Dans un second temps, nous avons regardé si les conditions d'obtention du précipité PuCrO3 correspondaient à un accroissement de l'homogénéité de répartition du plutonium et une taille de grains maximale. Pour ce faire, des échantillons fabriqués avec ou sans présence de chrome et frittés sous différentes atmosphères ont été étudiés. Il a été mis en évidence que la cinétique d'interdiffusion U-Pu est complètement modifiée en présence de chrome. De plus, suite à l'ajout de chrome, les conditions permettant d'accroître la cinétique d'interdiffusion U-Pu ne sont pas forcément associées à une taille de grain maximale.A partir de ces résultats, des préconisations pour la mise en œuvre industrielle sont proposées. Elles concernent le choix de l'atmosphère de frittage et la teneur en chrome nécessaire à l'optimisation de la microstructure.

La variabilité intrasaisonnière de la mousson d'Afrique de l'Ouest se caractérise par une alternance de phases sèches et humides, dont les impacts pe uvent être dramatiques sur les populations locales. Cette variabilité met en jeu un grand nombre d'échelles spatiales et temporelles, rendant difficile sa compréhension, sa modélisation et sa prévision. Cette thèse propose quelques éclairages sur ces différentes thématiques. La dépression thermique saharienne est un acteur majeur de la mousson africaine. La caractérisation de sa variabilité intrasaisonnière a permis de mettre en évidence, à l'échelle de 15 jours, l'existence d'interactions entre les latitudes moyennes et l'Afrique de l'Ouest. Lors de son passage au-dessus de l'Atlantique et la Méditerranée, un train d'ondes de Rossby module les ventilations de la dépression thermique, et donc sa structure. Les anomalies de circulation, de température et d'humidité, ainsi induites sur le Sahel, pourraient alors expliquer une partie des fluctuations intrasaisonnières de la convection, notamment celles qui naissent sur l'est du Sahel, et qui se propagent ensuite vers l'ouest. L'état moyen et la variabilité intrasaisonnière de la mousson africaine restent un défi pour les modèles de climat, même pour la dernière génération, qui a participé à l'exercice d'intercomparaison CMIP3. La variabilité à haute fréquence de la convection est un élément particulièrement difficile à modéliser. Toutefois, la meilleure prise en compte de facteurs inhibant le développement de la convection pourrait être une étape importante pour améliorer la modélisation de la mousson et la prévision de ses fluctuations intrasaisonnières

La conception des ouvrages de génie civil dans les massifs rocheux fracturés nécessite une connaissance de leur comportement lors des travaux. Les ma ssifs fracturés sont généralement très complexes et représentés par des matrices et des discontinuités rocheuses. Cette thèse a pour objectif d'étudier le comportement mécanique du massif du futur tunnel de St Béat. La thèse est divisée en trois parties. Nous étudions tout d'abord les phénomènes prépondérants et les facteurs agissant sur le comportement mécanique des massifs rocheux. Dans la seconde partie, nous effectuons séparément des études expérimentales sur le comportement mécanique de la matrice et des discontinuités rocheuses. Enfin, dans la dernière partie, nous présentons une modélisation en 2D du comportement du massif fracturé lors du creusement du tunnel. L'étude expérimentale s'appuie sur des marbres rencontrés sur le site du tunnel. Le comportement de la matrice rocheuse est étudié à partir des essais de compression triaxiale sur des échantillons de roche saine. L'évolution de l'endommagement est caractérisée à l'aide des mesures de vitesses de propagation des ondes élastiques. Une enveloppe de rupture selon le critère linéaire de Mohr-Coulomb est proposée pour la matrice rocheuse. Le comportement au cisaillement des discontinuités naturelles est ensuite étudié sous différentes conditions de chargement (contrainte normale constante et rigidité normale imposée). La mesure de topographie des surfaces des discontinuités avant et après chaque essai mécanique permet de déterminer les paramètres statistiques de la rugosité. L'influence de la contrainte normale, la rigidité normale imposée, la rugosité initiale et la vitesse de cisaillement sur le comportement des discontinuités est mise en évidence. Des lois de comportement pertin entes sont proposées pour chaque type d'essais.Les caractéristiques mécaniques obtenues sont alors introduites dans le code de calcul UDEC afin de modéliser le comportement du massif en présence du tunnel, sous différentes conditions aux limites. La réponse en terme de déformations et de contraintes induites autour du tunnel est analysée

De nombreuses applications en océanographie et en météorologie côtière et opérationnelle nécessitent la mise en place de modèles locaux haute-résolution pour lesquels les interactions océan-atmosphère doivent être correctement représentées. Dans ce cas, la principale difficulté est de raccorder de manière consistante la solution des deux modèles à l'interface air-mer. Le travail présenté ici vise à adapter les méthodes de décomposition de domaine de type algorithmes de Schwarz pour résoudre ce problème d'une manière mathématiquement satisfaisante. Une difficulté importante provient notamment du caractère fondamentalement turbulent des couches limites près de l'interface air-mer. Pour aborder cet aspect, nous présentons dans un premier temps une synthèse des divers schémas de paramétrisation utiles. Puis nous étudions une formulation idéalisée du problème, sous la forme d'un couplage simplifié de deux équations de diffusion modélisant le mélange turbulent dans les couches limites, les coefficients de diffusion turbulente étant fournis par des paramétrisations usuelles. Afin d'assurer une convergence rapide de la méthode, nous recherchons des conditions de transmission optimisées. Nous déterminons celles-ci de manière analytique dans le cas de coefficients de diffusion constants et discontinus à l'interface, puis nous étendons ces résultats dans le cas de coefficients variables à l'aide d'une approche nouvelle. Dans une dernière partie, nous montrons comment les méthodes usuelles de couplage océan-atmosphère peuvent être décrites dans le formalisme des méthodes de Schwarz. Puis nous proposons une première application réaliste (formation et propagation d'un cyclone tropical) et présentons des résultats numériques préliminaires, obtenus avec une méthode non optimisée.

Un enjeu majeur des modèles océaniques est de représenter fidèlement les circulations méso- et subméso-échelles afin de simuler leur importante contribution dans la circulation générale et dans le budget énergétique de l'océan. La poursuite de cet objectif se traduit par une augmentation de la résolution spatiale et temporelle à la fois des modèles et des réseaux d'observation de l'océan. Cependant, à ces petites échelles, la dynamique de l'écoulement revêt un caractère fortement turbulent ou non-linéaire. Dans ce contexte, les méthodes actuelles d'assimilation de données (AD), variationnelles en particulier, sont généralement moins performantes que dans un contexte (quasi-) linéaire. L'objectif de cette thèse est d'explorer sous divers aspects le comportement des méthodes variationnelles d'AD dans un modèle d'océan non-linéaire. Pour ce faire, nous avons réalisé une série d'expériences dites jumelles en assimilant des données altimétriques simulées suivant les caractéristiques des satellites altimétriques Jason-1 et SARAL/AltiKA . À l'aide de ces expériences, nous analysons sous différents angles les problématiques posées par les non-linéarités à l'AD. Enfin, nous ouvrons plusieurs pistes d'amélioration de l'efficacité du système d'AD dans ce contexte. Ce travail est basé sur le logiciel de modélisation océanique NEMO, incluant la configuration de bassin océanique turbulent idéalisé SEABASS, à différentes résolutions spatiales. Dans la continuité de la plateforme de recherche en AD avec NEMO, NEMO-ASSIM, nous avons utilisé et contribué au développement de cet ensemble d'outil, comprenant, entre autre, opérateur d'observation, modèles linéaire tangent et adjoint de NEMO, permettant de mener à bien notre étude. Le système d'AD variationnelle utilisé est le logiciel NEMOVAR. Les résultats présentés tentent de lier les échelles caractéristiques des structures d'erreurs d'analyse et l'activité aux petites échelles. Pour ce faire, nous avons utilisé une large gamme de diagnostics, e.g. erreur quadratique moyenne spatiale et temporelle, caractéristiques des fonctions coûts, caractérisation de l'hypothèse linéaire tangente, PSD des champs d'erreurs d'analyse. Nos expériences montrent que le 4DVAR incrémental contrôle efficacement la trajectoire analysée au 1/4° pour de longues fenêtres d'AD (2 mois). Lorsque la résolution augmente, la convergence de l'algorithme apparaît plus lente voire inexistante sous certaines conditions. Cependant, l'algorithme permet encore de réduire convenablement l'erreur d'analyse. Enfin, l'algorithme 3DFGAT se révèle beaucoup moins performant, quelle que soit la résolution. De plus, nous montrons également l'importance de l'adéquation entre la circulation simulée et l'échantillonnage altimétrique, en terme d'échelles spatiales représentées, pour obtenir de meilleures performances. Enfin, nous avons exploré la stratégie de minimisation dite progressive, permettant d'accélérer la convergence du 4DVAR à haute résolution.

Le carbone-suie, ou “black carbon” (BC), contribue au réchauffement climatique avec un forçage positif de l’ordre de +1.1 W m-2 dont l’incertitude reste haute (de l’ordre de 90%). Ce forçage s’effectue à travers l’interaction aérosol-radiation et l’interaction aérosols-nuage. Ces deux mécanismes sont affectés par le degré de mélange des particules du BC avec divers matériaux non-réfractaires et non-absorbants. Cependant, les estimations du forçage radiatif considèrent rarement les effets du mélange interne. Par ailleurs le rôle du BC comme noyau glaçogène qui influence l’interaction aérosol-nuage est largement inconnu. L’objectif de cette thèse est de mieux comprendre les mécanismes par lesquels le degré de mélange interne du BC influence la variabilité des propriétés optiques du BC et les propriétés d’activation des noyaux glaçogènes contenant du BC.Dans le premier chapitre de cette thèse, nous avons exploré la variabilité spatiale et saisonnière du coefficient d’absorption massique -mass absorption cross-section (MAC)- dans l’atmosphère en Europe. Les valeurs de MAC sont déterminées à partir de concentrations de carbone élémentaire et de coefficients d’absorption observée à différentes stations d’observation européenne du réseau ACTRIS (Aerosol, Cloud and Trace gases Research InfraStructure). Les résultats montrent une faible variabilité spatiale du MAC avec une moyenne de 10 ± 2.5 m2 g-1 à 637 nm de longueur d’onde qui peut être considérée comme représentative du BC en Europe. Le cycle saisonnier du MAC est probablement lié à la composition chimique de l’aérosol et son état de mélange, qui provoque une augmentation du MAC.Dans le second chapitre on s’est intéressé au lien entre l’absorption du BC et son état de mélange après transport sur longue-distance. Ce travail se base sur des mesures effectuées dans le cadre du projet CLIMSLIP (CLimate IMpact of Short-Lived Pollutants and methane in the Arctic). Une campagne de mesure a été conduite sur la station de recherche Zeppelin au Svalbard, Norvège en Avril 2012. Les données acquises avec un Single Particle Soot Photometer (SP2) révélaient que le BC est généralement présent en mélange interne dont l’épaisseur moyenne de la couche superficielle de matériel non-absorbant est de 47 nm pour des particules de BC de diamètre compris entre 170 et 280 nm. Ce mélange interne conduit à une augmentation d’absorption de 46%. Elle entraîne cependant une diminution relativement faible de l’albédo de simple diffusion, de l’ordre de 1%.Enfin, la capacité du BC à agir comme noyaux glaçogène pour la formation de cristaux de glace a été étudiée sur le site de haute altitude du Jungfraujoch (Suisse) dans le cadre du “cloud and aerosol characterization experiment” (CLACE) en 2013. Les différents éléments du nuage étaient séparé à partir d’une prise d’entrée type ice-CVI connectée au SP2. Ce dispositif permet de sélectionner uniquement les cristaux de glace et quantifier la fraction de BC activée. Une réduction de la présence de BC dans les résidus de glace a été observée. Des mesures de l’épaisseur de la couche de mélange interne des particules contentant du BC ont montré que les résidus de cristaux de glace présentaient des enrobages bien plus épais comparée à l’aérosol total.Les résultats obtenus au cours de ce travail ont permis de mieux comprendre l’impact du degré de mélange interne sur les propriétés optiques du BC et sur son rôle dans la formation de cristaux de glace. Les propriétés optiques du BC évoluent en fonction de la saison, tandis que la formation d’une couche superficielle amplifie sa capacité d’absorption du rayonnement solaire. De plus, cette étude souligne l’importance du vieillissement atmosphérique du BC sur sa capacité à servir de noyau de nucléation de la glace. Enfin, il fournit une avancée au sujet des propriétés sensibles mesurées dans l’atmosphère avec des techniques innovantes qui permettront la simulation plus précise du forçage radiatif. / Black carbon (BC) induces a warming effect (RFBC = +1.1 W m-2 ± 90%) through two main pathways: aerosol-radiation interaction (RFari) and aerosol-cloud interaction (RFaci). Both BC-radiation and BC-cloud interaction are affected by the mixing of black carbon with other non-refractory and non-absorbing matter present in the atmosphere. An estimation of the global radiative forcing of BC rarely accounts for internal mixing of BC while the net global cloud radiative forcing is sensitive to assumptions in the initiation of cloud glaciation, which is mostly unknown for black carbon particles. Within this thesis we investigated the variability of the light absorbing properties of black carbon, the mixing of black carbon, and the impact on light absorption and ice activation.In the first part of this thesis we investigated the spatial and seasonal variability of the mass absorption cross section (MAC) over Europe. MAC values were determined from ambient observations of elemental carbon mass concentrations (mEC) and absorption coefficients (σap). The data had been acquired during several years at different background ACTRIS supersites spread over Europe. Site specific MAC values were found to be spatially homogeneous, suggesting that the overall MAC average 9.5 ± 1.9 m2 g-1 at a wavelength of 637 nm might be representative of BC at European background locations. The MAC values showed a distinct seasonal cycle at every station. This seasonality might be related to chemical composition and aging. We observed that the MAC value has a linear and positive proportionality with the non-absorbing matter mass fraction.The second part of the work focuses on the coating acquisition of BC and the induced absorption enhancement after long-range transport. Within the CLIMSLIP (climate impact of short-lived pollutants and methane in the Arctic) project field experiments were conducted at the Zeppelin research site in Svalbard, Norway, during the Arctic spring. SP2 data were used to characterize the BC size distribution and mixing. BC containing particles having a core diameter between 170 and 280 nm were found to have a median coating thickness of 47 nm. The relationship between coating thickness and BC absorption was simulated. The observed coating thickness enhanced the mass absorption cross section by 46%, which led to a decrease of less than 1% in the single scattering albedo.In the final part of this work, the role of black carbon as ice nuclei in mixed phase clouds was investigated at the high elevation measuring site Jungfraujoch (Switzerland) during the cloud and aerosol characterization experiment (CLACE) held in 2013. The ice-CVI inlet and a single particle soot photometer were used to select and quantify the ice activated BC particles. According to the observations, BC containing particles were depleted in the ice residuals. The activation efficiency showed a size dependency, with larger BC containing particles being activated more efficiently compared to smaller ones. Activated BC cores having a diameter between 170 and 240 nm showed a larger coating thickness (median = 53 nm) compared to the total aerosol (median = 16 nm).The results obtained in this thesis shed new light on the effect of the mixing state on the optical properties and cloud activation of black carbon particles. Absorbing properties of BC showed a distinct seasonal pattern, while aging was found to consistently increase its absorption behavior. However, black carbon was found not to act as ice nuclei in low tropospheric mixed-phase clouds, where the coating thickness might play a role in the activation efficiency. This work provides freshly determined physical properties derived from ambient observations that will improve the accuracy of future aerosol and cloud radiative forcing estimations.

No description available.

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.