Cette thèse présente une étude théorique d'un modèle d'atmosphère stellaire, modélisée comme une couche plan-parallèle irradiée sur une face, avec des électrons non thermalisés a priori. Les électrons sont caractérisés par leur fonction de distribution des vitesses (fdv), que l'on cherche à calculer en même temps que les autres grandeurs de l'atmosphère. Notre principal objectif est de comprendre le mécanisme de thermalisation des électrons, qui tend à rapprocher leur fdv de la fonction de Maxwell-Boltzmann lorsque les collisions élastiques dominent les interactions inélastiques des électrons avec le milieu ambiant, une hypothèse universellement admise en théorie des atmosphères stellaires. Les processus inélastiques (collisionnels ou radiatifs) perturbent cet équilibre, et la fdv des électrons peut s'écarter considérablement de l'équilibre maxwellien aux hautes énergies. De tels écarts modifient fortement les populations atomiques et le champ radiatif. Les calculs numériques consistent en la comparaison de trois modèles d'atmosphères: en équilibre thermodynamique local (ETL), hors ETL avec électrons thermalisés, et hors ETL avec électrons non thermalisés a priori. Nous avons résolu ce problème dans un plasma d'hydrogène pur en prenant en compte les principaux types d'interaction présents dans les atmosphères stellaires. L'équation cinétique des électrons a été résolue en calculant son terme de collision élastique à l'aide d'un modèle BGK longuement justifié dans la thèse. Notre principale contribution se situe au niveau du transfert de rayonnement. Nous avons utilisé, et surtout développé, les codes de l'équipe "Transfert" de l'Observatoire de Lyon. Les calculs montrent que la fdv des électrons s'écarte considérablement d'une maxwellienne dans la région hors ETL de l'atmosphère stellaire. Pour conclure, nous envisageons quelques extensions possibles de ce travail et certaines applications astrophysiques.

Cette thèse porte sur l'étude des atmosphères de Titan et de la Terre primitive avec des modèles de circulation générale (GCM). Tout d'abord, j'ai analysé la structure thermique et la dynamique de la basse troposphère de Titan. Cette étude a abouti à une caractérisation complète de la couche limite et a révélé l'existence d'une circulation de couche limite, qui impacte tous les aspects de la météorologie titanienne (régimes de vents, ondes, formation des dunes et des nuages, échanges de moment cinétique et superrotation). A partir de cette analyse, j'ai proposé une nouvelle hypothèse pour expliquer l'orientation vers l'est des dunes de Titan grâce à un couplage entre les orages tropicaux et la superrotation. Ceci a été validé par des simulations méso-échelles et a permis de proposer un schéma global expliquant la formation des dunes et leurs différentes caractéristiques. J'ai ensuite participé au développement d'un GCM générique, conçu pour étudier tout type d'atmosphère. Je l'ai appliqué aux paléoclimats de Titan pour simuler une période où l'atmosphère a pu être dépourvue de méthane. Dans ce cas, le climat devait être différent d'aujourd'hui avec potentiellement des conséquences géologiques fondamentales notamment pour l'érosion et l'âge de la surface. Finalement, j'ai appliqué ce GCM générique au cas de la Terre primitive. J'ai montré que, malgré un soleil moins lumineux qu'aujourd'hui et des quantités de gaz à effet de serre contraintes par les archives minéralogiques, le climat de la Terre Archéenne a pu être tempérée. En particulier, grâce à une rétro-action nuageuse, la Terre aurait pu éviter une glaciation globale et rester propice au développement de la vie. / This thesis focuses on the study of the atmospheres of Titan and the early Earth with Global Climate Models (GCM). First, I analysed the thermal structure and the dynamics of Titan's lower troposphere. This analysis allowed a full caracterization of the planetary boundary layer and revealed the existence of a boundary layer circulation which impacts every aspect of Titan's weather (wind patterns, atmospheric waves, dune and cloud formation, exchange of momentum with the surface, and development of the superrotation). Thanks to this study, I proposed a new hypothesis to explain the eastward orientation of Titan's dunes that implies a coupling between tropical storms and the superrotation. This has been validated with mesoscale simulations and provided a general framework to explain Titan's dune formation and features. Then, I participated to the development of a generic GCM, designed to study any kind of atmosphere. I applied it to Titan's paleoclimates, when the atmosphere was depleted of methane. In such a case, the climate should have been different from today, with potentially fundamental geological consequences, in particular for the erosion and the age of the surface. Finally, I applied this GCM to the case of the early Earth using greenhouse gas abundances constrained by mineralogical data. I showed that despite a weaker solar insolation, the Archean Earth's climate may have been temperate. In particular, the Earth may have avoided a full glaciation and remained suitable for the development of life thanks to cloud feedback, even assuming a amount of CO2 just a little larger than today.

Lors de cette thèse, je me suis attaché à améliorer notre connaissance des planètes géantes, depuis notre voisine Jupiter jusqu’aux exoplanètes lointaines : les Jupiter chauds. Grâce aux nouvelles observations gravitationnelles extrêmement fines du satellite Juno, entré en orbite autour de Jupiter en juillet 2016, il est possible d’améliorer significativement les modèles de structure interne de la planète. Cependant, cela ne peut se faire qu’à condition d’avoir une méthode suffisamment précise pour exploiter au maximum les données. J’ai donc étudié la méthode des sphéroides de Maclaurin concentriques et ses limitations. A l’aide des connaissances contemporaines sur les équations d’état, les propriétés diffusives et les transition ou séparation de phase entre l’Hydrogène et l’Hélium, il m’a alors été possible de produire de nouveaux modèles de Jupiter. Arriver à combiner les observations gravitationnelles de Juno et les abondances d’éléments observées par Galiléo n’a pu se faire qu’en décomposant Jupiter en au moins 4 zones, de l’enveloppe externe jusqu’au coeur compact. J’ai montré que la taille de ce coeur compact était dégénérée avec la variation d’entropie à l’intérieur de la planète.La structure interne des Jupiter chauds quant à elle est très dépendante de leur dynamique atmosphérique, qui entraîne une inflation de leur rayon. J’ai étudié les atmosphères de ces planètes à l’aide du modèle de circulation globale de l’Université d’Exeter et d’un code linéaire que j’ai développé, appelé ECLIPS3D. La caractéristique la plus importante de la circulation atmosphérique est la présence d’un jet superrotatif, étendu en latitude.J’ai donc étudié la création de ce jet à l’aide d’arguments théoriques pour s’assurer de sa pertinence physique. L’étude de la solution linéaire dépendante du temps, associée à des arguments numériques sur la convergence de quantité de mouvement par les vents verticaux m’ont permis d’établir une compréhension globale, cohérente de l’accélération de la superrotation dans l’atmosphère de ces planètes.Avec ce travail, j’ai amélioré ma compréhension théorique des planètes géantes et développé des codes qui peuvent être utilisés pour améliorer nos connaissances sur la structure interne et la dynamique atmosphérique des planètes géantes, que ce soit Jupiter, Saturne ou les Jupiter chauds. / Through this thesis, I have been motivated by the will to improve our knowledge of giant planets, from our neigh- bouring Jupiter to the far away worlds across the galaxy: hot Jupiters.With the latest, extremely precise observations of the satellite Juno, new models of the interior of Jupiter can be derived. A precise enough method is required to take full advantage of these outstanding data, and I therefore studied the concentric Maclaurin spheroid method and its limitations.With contemporary understanding on the equations of state, diffusive properties and phase transition/separation of hydrogen and helium, I could then focus on producing new interior models of Jupiter. Combining the gravitational observations of Juno with the elemental observations of Galileo has proven to be a complicated task, which required to decompose the planet into at least four regions from the outer envelope to the inner, compact core. I have shown that the size of the compact core is degenerated with the entropy variation within the planet.Concerning hot Jupiters, I have reminded of the need to understand their atmospheric dynamics to constrain their interior structure, as the wind circulation can lead to an inflation of their radius. Studying numerically their at- mospheric dynamics was performed with the University of Exeter’s global circulation model as well as with the development of a linear solver that I called ECLIPS3D. An important, robust feature is the presence of a broad equatorial superrotation in the atmosphere of these planets.Finally, I have explored the spin up of this superrotation on theoretical grounds, to assess its physical relevance. I have calculated the linear time dependent solution to show the importance of differential drag and radiative damp- ing, and have used numerical simulations to highlight the importance of vertical momentum acceleration. Globally, a coherent picture of the initial spin up of superrotation was obtained.Through this work, I have improved my theoretical understanding of giant planets and developed various codes that can be used to study and improve our knowledge of the interior structure and atmospheric dynamics of giant planets, from Jupiter and Saturn to hot Jupiters.

Cette thèse porte sur l'étude des atmosphères de Titan et de la Terre primitive avec des modèles de circulation générale (GCM). Tout d'abord, j'ai analysé la structure thermique et la dynamique de la basse troposphère de Titan. Cette étude a abouti à une caractérisation complète de la couche limite et a révélé l'existence d'une circulation de couche limite, qui impacte tous les aspects de la météorologie titanienne (régimes de vents, ondes, formation des dunes et des nuages, échanges de moment cinétique et superrotation). A partir de cette analyse, j'ai proposé une nouvelle hypothèse pour expliquer l'orientation vers l'est des dunes de Titan grâce à un couplage entre les orages tropicaux et la superrotation. Ceci a été validé par des simulations méso-échelles et a permis de proposer un schéma global expliquant la formation des dunes et leurs différentes caractéristiques. J'ai ensuite participé au développement d'un GCM générique, conçu pour étudier tout type d'atmosphère. Je l'ai appliqué aux paléoclimats de Titan pour simuler une période où l'atmosphère a pu être dépourvue de méthane. Dans ce cas, le climat devait être différent d'aujourd'hui avec potentiellement des conséquences géologiques fondamentales notamment pour l'érosion et l'âge de la surface. Finalement, j'ai appliqué ce GCM générique au cas de la Terre primitive. J'ai montré que, malgré un soleil moins lumineux qu'aujourd'hui et des quantités de gaz à effet de serre contraintes par les archives minéralogiques, le climat de la Terre Archéenne a pu être tempérée. En particulier, grâce à une rétro-action nuageuse, la Terre aurait pu éviter une glaciation globale et rester propice au développement de la vie.

Les atmosphères de Jupiter chaudes reçoivent des quantitées extrêmes de radiation de leur étoile hôte, ce qui nous permet d’étudier les processus chimiques et dynamiques d’atmosphères dans un régime d’énergie qui n’est pas accessible par l’étude du système solaire. La spectroscopie de transit et d’éclipse secondaire permet la caractérisation de ces atmosphères. Toutefois, ces méthodes ne considèrent pas la variation spatiale des abondances moléculaires, du profil de température et des nuages dans l’atmosphère. À l’aide de la méthode de cartographie par éclipse secondaire à plusieurs longueurs d’ondes, il est possible de mesurer la distribution spatiale des paramètres atmosphériques. À ce jour, cette méthode n’a été appliquée qu’à la Jupiter chaude HD 189733 b à partir d’observations photométriques du télescope spatial Spitzer à 8 μm. Avec le lancement prochain du JWST, des observations à haute précision spectrophotométrique sont attendues, ce qui renforce le besoin pour une méthode permettant la simulation et l’analyse statistique d’observations d’éclipse secondaire à plusieurs longueurs d’ondes. Nous présentons une nouvelle méthode afin de simuler des courbes de lumière à plusieurs longueurs d’ondes, ainsi qu’une méthode statistique pour déterminer la distribution atmosphérique à partir d’observations. La simulation des courbes de lumière à plusieurs longueurs d’ondes est discrétisée pour permettre le calcul rapide des modèles nécessaire à l’analyse statistique d’observations. Nous appliquons cette méthode à des observations simulées de JWST NIRSpec pour la planète WASP-33 b à diverses distributions spatiales de coefficient de dissociation de l’eau. Nos résultats indiquent que 5 observations d’éclipse secondaire sont suffisantes pour mesurer la distribution spatiale du coefficient de dissociation de l’eau en longitude et latitude. L’application de la méthode de cartographie par éclipse secondaire à plusieurs longueurs d’ondes va permettre d’appronfondir d’avantage notre connaissance des processus dynamiques des atmosphères de Jupiter chaudes. / Hot Jupiter atmospheres, receiving extreme amounts of irradiation from their host star, enable us to study the atmospheric chemistry and dynamics of atmospheres at regimes that would otherwise be inaccessible from solar system observations. Transmission and emission spectroscopy of these atmospheres enables their characterization. However, these methods do not consider the spatial distribution of the molecular abundances, temperature profile, and clouds. Using multi-wavelength secondary eclipse mapping, it is possible to measure the spatial distribution of the atmospheric parameters. Secondary eclipse mapping has only been applied for photometric Spitzer/IRAC 8 μm observations of HD 189733 b so far. With the imminent launch of JWST, high-precision spectrophotometric observations of hot Jupiters are expected, reinforcing the need for a framework that can simulate and retrieve secondary eclipse mapping observations. We present a novel method to simulate secondary eclipse light curves for non-uniform atmospheres at multiple wavelengths, as well as a statistical method to retrieve the distribution from observations. The simulation of the multi-wavelength light curves is discretized to enable the fast computation necessary for statistical analysis of observations. We apply these methods to simulated JWST NIRSpec observations of WASP-33 b for various distributions of the water dissociation coefficient and conclude that 5 secondary eclipse observations are sufficient to constrain significant spatial variations in longitude and latitude. The application of the multi-wavelength secondary eclipse mapping to JWST will further deepen our understanding of the atmospheric dynamics processes to which hot Jupiter atmospheres are subjected.

Le but de cette thèse est de raffiner et de mieux comprendre l'utilisation de la méthode spectroscopique, qui compare des spectres visibles de naines blanches à atmosphère riche en hydrogène (DA) à des spectres synthétiques pour en déterminer les paramètres atmosphériques (température effective et gravité de surface). Notre approche repose principalement sur le développement de modèles de spectres améliorés, qui proviennent eux-mêmes de modèles d'atmosphère de naines blanches de type DA. Nous présentons une nouvelle grille de spectres synthétiques de DA avec la première implémentation cohérente de la théorie du gaz non-idéal de Hummer & Mihalas et de la théorie unifiée de l'élargissement Stark de Vidal, Cooper & Smith. Cela permet un traitement adéquat du chevauchement des raies de la série de Balmer, sans la nécessité d'un paramètre libre. Nous montrons que ces spectres améliorés prédisent des gravités de surface qui sont plus stables en fonction de la température effective. Nous étudions ensuite le problème de longue date des gravités élevées pour les DA froides. L'hypothèse de Bergeron et al., selon laquelle les atmosphères sont contaminées par de l'hélium, est confrontée aux observations. À l'aide de spectres haute résolution récoltés au télescope Keck à Hawaii, nous trouvons des limites supérieures sur la quantité d'hélium dans les atmosphères de près de 10 fois moindres que celles requises par le scénario de Bergeron et al. La grille de spectres conçue dans ces travaux est ensuite appliquée à une nouvelle analyse spectroscopique de l'échantillon de DA du SDSS. Notre approche minutieuse permet de définir un échantillon plus propre et d'identifier un nombre important de naines blanches binaires. Nous déterminons qu'une coupure à un rapport signal-sur-bruit S/N > 15 optimise la grandeur et la qualité de l'échantillon pour calculer la masse moyenne, pour laquelle nous trouvons une valeur de 0.613 masse solaire. Finalement, huit nouveaux modèles 3D de naines blanches utilisant un traitement d'hydrodynamique radiative de la convection sont présentés. Nous avons également calculé des modèles avec la même physique, mais avec une traitement standard 1D de la convection avec la théorie de la longueur de mélange. Un analyse différentielle entre ces deux séries de modèles montre que les modèles 3D prédisent des gravités considérablement plus basses. Nous concluons que le problème des gravités élevées dans les naines blanches DA froides est fort probablement causé par une faiblesse dans la théorie de la longueur de mélange. / The goal of this thesis is to refine and to understand better the spectroscopic method, which compares optical spectra of hydrogen-atmosphere white dwarfs (DA) with synthetic spectra to determine the atmospheric parameters (effective temperature and surface gravity). Our approach rests mainly on the development of improved model spectra, which come themselves from DA model atmospheres. We present a new grid of DA synthetic spectra with the first consistent implementation of the non-ideal gas theory of Hummer & Mihalas and the unified theory of Stark broadening from Vidal, Cooper & Smith. This allows for an adequate treatment of the quenching effects in Balmer lines, without the need of a free parameter. We show that these improved spectra predict surface gravities that are much more stable as a function of the effective temperature. We then study the long-standing problem that surface gravities in cool DA stars are significantly higher than those found in hotter DA white dwarfs. The hypothesis of Bergeron et al., according to which the atmospheres are contaminated by small amounts of helium, is constrained with observations. Using high-resolution spectra collected at the Keck in Hawaii, we find superior limits on the helium abundances in the atmospheres that are nearly 10 times lower than those required to sustain the Bergeron et al. scenario. The grid of spectra calculated in this work is then applied to a new spectroscopic analysis of the DA in the SDSS. Our careful analysis allows us to define a cleaner sample and to identify a large number of double degenerates. We find that a cutoff at a signal-to-noise ratio S/N > 15 optimizes the size and quality of the sample for computing the average mass, for which we find a value of 0.613 solar mass. Finally, eight new 3D white dwarf models with a radiative-hydrodynamics treatment of the convection are presented. We also calculated models with the same physics, except for a treatment of the convection with the standard mixing-length theory. A differential analysis between these two sets of models shows that the 3D models predict considerably lower surface gravities. We conclude that the high-log g problem in cool DA white dwarfs is caused by a weakness in the mixing-length theory.

L'eau est un constituant essentiel de l'atmosphère de supergéantes rouges (RSG), mais dont l'influence reste mal comprise. Le spectre observé de l'eau de ces étoiles ne peut être reproduit que par l'ajout d'une coquille de gaz moléculaire, les MOLsphères. Cependant, l'hypothèse des MOLsphères reste fragile et sujette à caution. Dans le but de mieux interpréter les spectres observés, la synthèse de spectres hors équilibre thermodynamique local est une approche potentiellement importante. Les effets hors ETL étant potentiellement fort, ils pourraient être en mesure de lever les problèmes de l'interprétation des raies de l'eau sans ajout de MOLsphère et impliquer un rôle important dans la dynamique de l'atmosphère. Pour cela, nous avons développé une méthode originale en mesure de résoudre l'équation de transfert pour les nombreuses transitions radiatives de l'eau sans approximationETL. Cette méthode a été mise en oeuvre via le développement d'un code de transfert radiatif parallèle. Les premiers résultats montrent que les effets hors ETL dans l'atmosphère des RSG, et leur impact sur le spectre comme sur certaines observables utilisées pour sonder ces étoiles, sont importants / Water is an important constituent of the atmosphere of red supergiant stars (RSG), which influence remains however poorly understood. The water spectrum of these stars can apparently be only reproduced through the addition of a detached shell of cool molecular gas, the so-called MOLspheres. However, this hypothesis is still cautious. In order to better interpret observed spectra, non local thermodynamic equilibrium (NLTE) spectrum synthesis may be potentially important. NLTE effects being potentially important, they may alleviate the problems in the interpretation of water spectra, and affect the atmosphere dynamics. We thus developed an original method to solve the radiative transfer equation, adapted to the numerous water transitions and without the LTE approximation. This method has been implemented in an original parallel code. Preliminary results show that NLTE effects in RSG atmospheres and their impact on observables such as the emergent spectrum are very important.

Le but de cette thèse est de raffiner et de mieux comprendre l'utilisation de la méthode spectroscopique, qui compare des spectres visibles de naines blanches à atmosphère riche en hydrogène (DA) à des spectres synthétiques pour en déterminer les paramètres atmosphériques (température effective et gravité de surface). Notre approche repose principalement sur le développement de modèles de spectres améliorés, qui proviennent eux-mêmes de modèles d'atmosphère de naines blanches de type DA. Nous présentons une nouvelle grille de spectres synthétiques de DA avec la première implémentation cohérente de la théorie du gaz non-idéal de Hummer & Mihalas et de la théorie unifiée de l'élargissement Stark de Vidal, Cooper & Smith. Cela permet un traitement adéquat du chevauchement des raies de la série de Balmer, sans la nécessité d'un paramètre libre. Nous montrons que ces spectres améliorés prédisent des gravités de surface qui sont plus stables en fonction de la température effective. Nous étudions ensuite le problème de longue date des gravités élevées pour les DA froides. L'hypothèse de Bergeron et al., selon laquelle les atmosphères sont contaminées par de l'hélium, est confrontée aux observations. À l'aide de spectres haute résolution récoltés au télescope Keck à Hawaii, nous trouvons des limites supérieures sur la quantité d'hélium dans les atmosphères de près de 10 fois moindres que celles requises par le scénario de Bergeron et al. La grille de spectres conçue dans ces travaux est ensuite appliquée à une nouvelle analyse spectroscopique de l'échantillon de DA du SDSS. Notre approche minutieuse permet de définir un échantillon plus propre et d'identifier un nombre important de naines blanches binaires. Nous déterminons qu'une coupure à un rapport signal-sur-bruit S/N > 15 optimise la grandeur et la qualité de l'échantillon pour calculer la masse moyenne, pour laquelle nous trouvons une valeur de 0.613 masse solaire. Finalement, huit nouveaux modèles 3D de naines blanches utilisant un traitement d'hydrodynamique radiative de la convection sont présentés. Nous avons également calculé des modèles avec la même physique, mais avec une traitement standard 1D de la convection avec la théorie de la longueur de mélange. Un analyse différentielle entre ces deux séries de modèles montre que les modèles 3D prédisent des gravités considérablement plus basses. Nous concluons que le problème des gravités élevées dans les naines blanches DA froides est fort probablement causé par une faiblesse dans la théorie de la longueur de mélange. / The goal of this thesis is to refine and to understand better the spectroscopic method, which compares optical spectra of hydrogen-atmosphere white dwarfs (DA) with synthetic spectra to determine the atmospheric parameters (effective temperature and surface gravity). Our approach rests mainly on the development of improved model spectra, which come themselves from DA model atmospheres. We present a new grid of DA synthetic spectra with the first consistent implementation of the non-ideal gas theory of Hummer & Mihalas and the unified theory of Stark broadening from Vidal, Cooper & Smith. This allows for an adequate treatment of the quenching effects in Balmer lines, without the need of a free parameter. We show that these improved spectra predict surface gravities that are much more stable as a function of the effective temperature. We then study the long-standing problem that surface gravities in cool DA stars are significantly higher than those found in hotter DA white dwarfs. The hypothesis of Bergeron et al., according to which the atmospheres are contaminated by small amounts of helium, is constrained with observations. Using high-resolution spectra collected at the Keck in Hawaii, we find superior limits on the helium abundances in the atmospheres that are nearly 10 times lower than those required to sustain the Bergeron et al. scenario. The grid of spectra calculated in this work is then applied to a new spectroscopic analysis of the DA in the SDSS. Our careful analysis allows us to define a cleaner sample and to identify a large number of double degenerates. We find that a cutoff at a signal-to-noise ratio S/N > 15 optimizes the size and quality of the sample for computing the average mass, for which we find a value of 0.613 solar mass. Finally, eight new 3D white dwarf models with a radiative-hydrodynamics treatment of the convection are presented. We also calculated models with the same physics, except for a treatment of the convection with the standard mixing-length theory. A differential analysis between these two sets of models shows that the 3D models predict considerably lower surface gravities. We conclude that the high-log g problem in cool DA white dwarfs is caused by a weakness in the mixing-length theory.

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.

Cette thèse traite de la modélisation des marées fluides des planètes telluriques du système solaire et des systèmes exoplanétaires.En premier lieu, nous examinons la réponse de marée des couches atmosphériques, soumises au potentiel de marée gravifique et au forçage thermique de l’étoile hôte du système. Nous proposons un nouveau modèle global prenant en compte les processus dissipatifs avec un refroidissement newtonien, modèle à partir duquel nous traitons la dynamique des ondes de marées engendrées par ces forçages, et quantifions leur dissipation, le nombre de Love et le couple de marée exercé sur la couche atmosphérique en fonction de la fréquence de forçage. Ceci nous permet d'étudier l'ensemble des configurations possibles depuis les planètes au voisinage de la synchronisation telles que Vénus jusqu'aux rotateurs rapides tels que la Terre.En second lieu, nous développons une approche similaire pour les océans de planètes de type terrestre, où la friction visqueuse effective de la topographie est prise en compte, à partir de laquelle nous quantifions la réponse de marée d’un océan global potentiellement profond et sa dépendance à la fréquence d’excitation. Dans ce cadre, et ce grâce à des modèles locaux, nous caractérisons de manière détaillée les propriétés des spectres en fréquence de la dissipation engendrée par les ondes de marées au sein des couches fluides planétaires (et stellaires) en fonction des paramètres structurels et dynamiques de ces dernières (rotation, stratification, viscosité et diffusivité thermique). / This thesis deals with the modeling of fluid tides in terrestrial planets of the Solar system and exoplanetary systems.First, we examine the tidal response of atmospheric layers, submitted to the tidal gravitational potential and the thermal forcing of the host star of the system. We propose a new global model taking into account dissipative processes with a Newtonian cooling, model that we use to treat the dynamics of tidal waves generated by these forcings, and to quantify their dissipation, the Love number and the tidal torque exerted on the atmospheric layer as a function of the forcing frequency. This allows us to study possible configurations from planets close to synchronization such as Venus to rapid rotators such as the Earth.Second, we develop a similar approach for the oceans of terrestrial planets where the action of topography is taken into account thanks to an effective viscous friction. From this modeling, we quantify the tidal response of a potentially deep global ocean and its dependence of the tidal frequency. In this framework, and by using local models, we characterize in detail the properties of the frequency spectra of dissipation generated by tidal waves within fluid planetary (and stellar) layers as functions of the structural and dynamical parameters of these latters (rotation, stratification, viscosity and thermal diffusivity).