Une galaxie est un système complexe au sens où, autant des phénomènes se produisant à l'échelle du milieu interstellaire, comme des explosions de supernovæ ou l'activité d'un trou noir supermassif, que des interactions entre galaxies au sein de groupes ou d'amas, comme l'effeuillage par effet de marée ou par effet de bélier, influencent et conditionnent l'évolution de la galaxie dans son ensemble. Comme les processus œuvrant dans de tels systèmes font intervenir une gamme d'échelle de temps et de distance considérable, allant de l'étoile individuelle aux amas de galaxies tout entier, leur modélisation constitue un immense défi qui ne peut être relevé ni par une approche purement analytique ni par l'entremise de techniques exclusivement numériques.Cette thèse, à l'interface entre modèles semi-analytique et analyse de simulations numériques, se concentre sur l'étude de l'effeuillage des étoiles des galaxies satellites par effet de marée et sur les rétro-actions induites par les supernovæ.Ce manuscrit présente, d'une part, un modèle d'occupation des halos permettant de contraindre la masse d'étoiles perdue par les galaxies satellites depuis leur entrée dans leur groupe ou leur amas ainsi qu'un modèle d'effeuillage impulsif prédisant la masse stellaire arrachée aux satellites. Ce dernier est confronté, par le truchement du modèle d'occupation des halos, aux observations des fonctions de masses des groupes et des amas.Il expose, d'autre part, l'étude des rétro-actions des supernovæ implémenté dans les simulations numériques du projet NIHAO, conduite en séparant en différentes composantes le gaz des simulations et en comptabilisant les échanges entre ces dernières, laquelle a permis de mettre en évidence trois processus distincts par le biais desquels les supernovæ réduisent ou suppriment la formation stellaire.Enfin, il détaille les améliorations techniques et scientifiques apportées au modèle semi-analytique GalICS. / A galaxy is a complex system since as many phenomena take place at the scale of the interstellar medium, such as supernovae explosions or the activity of supermassive black holes, as interactions between galaxies within groups or clusters, such as tidal or ram pressure stripping, affect and condition the evolution of the galaxy itself as a whole. Because the processes acting in such systems involve a considerable range of times and distances, going from individual stars to entire clusters of galaxies, they modelling constitutes an immense challenge that cannot be met neither by a purely analytical approach nor by solely numerical technics.This thesis, being at the interface between semi-analytical models and the analysis of numerical simulations, focuses on the study of star stripping in satellite galaxies by tidal effects and on the supernovae induced feedback.This manuscript presents, on one hand, an halo occupation model that allows to constrain the stellar mass lost by satellite galaxies since they entered their group or their cluster, and a model of impulsive stripping that predicts the stellar mass ripped out of satellites. The latter is compared, through the halo occupation model, to the observed mass functions of groups and clusters.It exposes, on the other hand, the study of the supernovae feedback implemented in the numerical simulations of the NIHAO project, performed separating the simulated gas into different components and counting the exchanges that take place between them. This allowed for the highlighting of three distinct processes through which supernovae reduce or suppress their star formation.

L'objectif de cette thèse est de comprendre la structure des binaires gamma, binaires à collision de vents composées d'une étoile massive et d'un pulsar jeune. Ces binaires possèdent probablement une structure similaire aux binaires à collision de vents composées de deux étoiles massives, avec des particularités liées à la nature relativiste du vent de pulsar. L'interaction de deux vents supersoniques d'étoiles massives crée une structure choquée qui présente des signatures observationnelles du domaine radio aux rayons X. Plusieurs instabilités ainsi que le mouvement orbital des étoiles influent sur la structure choquée. Afin de comprendre leur impact, j'ai effectué des simulations à haute résolution de binaires à collision de vents à l'aide du code hydrodynamique RAMSES. Ces simulations sont numériquement coûteuses à réaliser, surtout lorsque un des vents domine fortement l'autre. A petite échelle, les simulations soulignent l'importance de l'instabilité de couche mince non-linéaire dans les collisions isothermes alors que l'instabilité de Kelvin-Helmholtz peut fortement modifier la structure choquée dans une collision adiabatique. A plus grande échelle, cette instabilité peut parfois détruire la structure spirale à laquelle on s'attend si la différence de vitesse entre les vents est trop importante. WR 104 est une binaire dont on observe la structure spirale grâce à l'émission de poussières. Les simulations de ce système montrent un bon accord avec la structure observée et indiquent que des processus de refroidissement du gaz sont nécessaires à la formation de poussières. Pour modéliser les vents de pulsar dans les binaires gamma, RAMSES a été étendu à l'hydrodynamique relativiste. J'utilise ce nouveau code pour réaliser des simulations préliminaires de binaires gamma. Elles montrent effectivement une structure similaire aux binaires stellaires, avec de légères corrections relativistes . Ce code est adapté à l'étude de divers systèmes astrophysiques tels que les jets relativistes, les sursauts gamma ou les nébuleuses de pulsar et fera partie de la prochaine version de RAMSES qui sera rendue publique.

L'objectif de cette thèse est de comprendre la structure des binaires gamma, binaires à collision de vents composées d'une étoile massive et d'un pulsar jeune. Ces binaires possèdent probablement une structure similaire aux binaires à collision de vents composées de deux étoiles massives, avec des particularités liées à la nature relativiste du vent de pulsar. L'interaction de deux vents supersoniques d'étoiles massives crée une structure choquée qui présente des signatures observationnelles du domaine radio aux rayons X. Plusieurs instabilités ainsi que le mouvement orbital des étoiles influent sur la structure choquée. Afin de comprendre leur impact, j'ai effectué des simulations à haute résolution de binaires à collision de vents à l'aide du code hydrodynamique RAMSES. Ces simulations sont numériquement coûteuses à réaliser, surtout lorsque un des vents domine fortement l'autre. A petite échelle, les simulations soulignent l'importance de l'instabilité de couche mince non-linéaire dans les collisions isothermes alors que l'instabilité de Kelvin-Helmholtz peut fortement modifier la structure choquée dans une collision adiabatique. A plus grande échelle, cette instabilité peut parfois détruire la structure spirale à laquelle on s'attend si la différence de vitesse entre les vents est trop importante. WR 104 est une binaire dont on observe la structure spirale grâce à l'émission de poussières. Les simulations de ce système montrent un bon accord avec la structure observée et indiquent que des processus de refroidissement du gaz sont nécessaires à la formation de poussières. Pour modéliser les vents de pulsar dans les binaires gamma, RAMSES a été étendu à l'hydrodynamique relativiste. J'utilise ce nouveau code pour réaliser des simulations préliminaires de binaires gamma. Elles montrent effectivement une structure similaire aux binaires stellaires, avec de légères corrections relativistes . Ce code est adapté à l'étude de divers systèmes astrophysiques tels que les jets relativistes, les sursauts gamma ou les nébuleuses de pulsar et fera partie de la prochaine version de RAMSES qui sera rendue publique. / The aim of this thesis is to understand the structure of colliding wind binaries composed of a massive star and a young pulsar, called gamma-ray binaries. They are expected to display a similar structure to colliding wind binaries composed of massive stars, with some particularities due to the relativistic nature of the pulsar wind. The interaction of the supersonic winds from massive stars creates a shocked structure with observational signatures from the radio domain to the X-rays. The structure is affected by various instabilities and by the orbital motion of the stars. To understand their impact, I carried out high resolution simulations of colliding wind binaries with the hydrodynamical code RAMSES. They are computationally demanding, especially when one of the winds strongly dominates the other one. Small scale simulations highlight the importance of the Non-linear Thin Shell Instability in isothermal collisions while the Kelvin-Helmholtz instability may strongly impact the dynamics of adiabatic collisions. I found that, at larger scales, this instability can destroy the expected large scale spiral structure when there is an important velocity gradient between the winds. WR 104 is a system that displays a spiral structure with important dust emission. The simulation of this system shows a good agreement with the observed structure and indicates cooling processes are necessary to enable dust formation. To model the pulsar wind in gamma-ray binaries, an extension of RAMSES has been developed, that incorporates relativistic hydrodynamics. I used this new relativistic code to perform preliminary simulations of gamma-ray binaries. They display a similar structure to colliding wind binaries with small relativistic corrections. We expect to use this code to perform large scale simulations of gamma-ray binaries. It will be part of the next public release of RAMSES and is suited for the study of many astrophysical problems such as relativistic jets, pulsar wind nebulae or gamma-ray bursts.

Alors que les oscillations solaires n'ont pas fini de nous révéler tous leurs secrets, des oscillations acoustiques analogues sont détectées dans un nombre croissant d'étoiles. Comme sur le Soleil, ces oscillations (dite de type solaire) sont amorties par des mécanismes complexes et encore mal connus et excitées par la turbulence dans l'enveloppe convective supérieure des étoiles. Grâce à la qualité photométrique exceptionnelle des missions spatiales CoRoT (CNES) et Kepler(NASA) ainsi qu'à la continuité long terme des observations qu'elles fournissent, on mesure maintenant précisément fréquences, amplitudes et durées de vie de ces oscillations dans une variété d'étoiles dotées de caractéristiques diverses concernant leur stade évolutif, paramètres fondamentaux, composition chimique, champ magnétique, rotation ... etc. Plus que ne le fait la mesure de leurs fréquences, la mesure des amplitudes et durées vies des modes de type solaire nous fournit des contraintes sur les propriétés statiques et dynamiques de la convection, sur la physique des modes et enfin sur la stratification en surface des étoiles. Le jeux conséquent d'étoiles pulsantes détectées par CoRoT et Kepler nous révéle aussi que les amplitudes et durées de ces oscillations varient d'une étoile à l'autre selon des lois d'échelles caractéristiques qui dépendent d'un nombre restreint de paramètres stellaires (masse, luminosité, température effective ... etc). Ce mémoire de thèse résume les travaux que j'ai menés dans ce contexte depuis plus de dix ans en collaboration avec mes collègues et avec les étudiants que j'ai encadré. Ces travaux ont cherché à comprendre et mieux modéliser les amplitudes des oscillations excitées par la convection turbulente, notamment les lois d'échelles observées. Ce faisant, ils ont visé à établir des diagnostics sur les propriétés statiques et dynamiques des régions convectives, avec pour objectif à plus long terme d'améliorer la modélisation des processus de transport convectif dans les intérieurs stellaires.

Les simulations hydrodynamiques pour la physique de haute densité d'énergie ainsi que pour la fusion par confinement inertiel exigent une description détaillée de flux d'énergie. Le mécanisme principal est le transport électronique, qui peut être un phénoméne non local qui doit être décrit avec des modèles de Fokker-Planck, stationnaires et simplifiés dans les codes hydrodynamiques à grande échelle. Mon travail thèse est consacré au développement d'un nouveau modèle de transport non local basé sur l'utilisation d'une méthode de fermeture entropique pour la résolution des premiers moments de l'équation de Fokker-Planck agrémentée d'un opérateur de collision dédié. Une telle fermeture permet une bonne résolution des fortes anisotropies de la fonction de distribution électronique dans les régimes où le développement d'instabilités électrostatiques à petite échelle le requiert. Ce modèle aux moments (M1) est comparé avec succès au modèle de Schurtz, Nicolaï et Busquet (SNB), référent dans le domaine du transport électronique non local. Ce modèle, basé sur l'hypothèse d'une faible anisotropie de la fonction de distribution sous-jacente induisant une relation de fermeture polynomiale (P1), utilise un opérateur de collision simplifié dont nous avons proposé une amélioration. Après avoir considéré plusieurs configurations typiques de transport de chaleur, nous avons montré que le modèle M1 ultidimensionnel peut prendre naturellement en compte des effets d'un plasmas magnétisés sur le transport électronique. De plus, ce modèle permet de calculer des fonctions de distribution utiles aux études cinétiques comme la stabilité du plasma dans la zone de transport. Nous confirmons avec notre modèle que le transport d'énergie électronique peut fortement modifier l'amortissement des ondes de Langmuir et des ondes acoustiques ; contrairement aux modèles non locaux simplifiés, M1 décrit les modifications de la fonction de distribution et l'amortissement des ondes du plasma. La structure du modèle permet également de prendre en compte naturellement des champs magnétiques autogénérés, qui jouent un rôle crucial dans des simulations multidimensionnelles. Ces champs magnétiques pourraient également être étudiés pour concentrer l'énergie dans les schémas d'ignition. Enfin, nous montrons que le modèle M1 reproduit les résultats de la théorie locale élaborée par Braginskii pour tous les niveau de magnétisation et propose de nouveaux résultats pour le régime non local. Ce travail constitue une première validation de l'utilisation des fermetures entropiques, dans les régimes de faibles anisotropies, qui va s'ajouter aux tests dans les régimes fortement anisotropes. / Hydrodynamic simulations in high-energy-density physics and inertial con nement fusion require a detailed description of energy uxes. The leading mechanism is the electron transport, which can be a nonlocal phenomenon that needs to be described with quasistationary and simplified Fokker-Planck models in large scale hydrodynamic codes. My thesis is dedicated to the development of a new nonlocal transport model based on a fast-moving-particles collision operator and on a first moment Fokker-Planck equation, simplified with an entropic closure relation. Such a closure enables a better description of the electron distribution function in the limit of high anisotropies, where small scale electrostatic instabilities could be excited. This new model, so called M1, is successfully compared with the well known nonlocal electron transport model proposed by Schurtz, Nicolaï and Busquet, using different collision operators, and with the reduced Fokker-Planck model, based on a small-anisotropies polynomial closure relation (P1). Several typical configurations of heat transport are considered. We show that the M1 entropic model may operate in two and three dimensions and is able to account for electron transport modifications in external magnetic fields. Moreover, our model enables to compute realistic electron distribution functions, which can be used for kinetic studies, as for the plasma stability in the transport zone. It is demonstrated that the electron energy transport may strongly modify damping of Langmuir and ion acoustic waves, while the simplified nonlocal transport models are not able to describe accurately the modifications of the distribution function and plasma wave damping. The structure of the M1 model allows to naturally take into account self-generated magnetic fields, which play a crucial role in multidimensional simulations. Moreover, magnetic fields could also be used for the focusing of energetic particles in alternative ignition schemes. The M1 model reproduces the results of the local transport theory in plasma, developed by Braginskii, in a broad range of degrees of magnetization and predicts new results in the nonlocal regime. This work constitutes a first validation of the entropic closure assumption in the weakly-anisotropic regime. It can be added to the existing tests, in the strongly-anisotropic regimes.

Des percées récentes ont été faites dans l’étude des trous noirs supermassifs (SMBH), grâce en grande partie à l’équipe du télescope de l’horizon des évènements (EHT). Cependant, déterminer la masse de ces entités colossales à des décalages vers le rouge élevés reste un défi de taille pour les astronomes. Il existe diverses méthodes directes et indirectes pour mesurer la masse de SMBHs. La méthode directe la plus précise consiste à résoudre la cinématique du gaz moléculaire, un traceur froid, dans la sphère d’influence (SOI) du SMBH. La SOI est définie comme la région où le potentiel gravitationnel du SMBH domine sur celui de la galaxie hôte. Par contre, puisque la masse d’un SMBH est négligeable face à la masse d’une galaxie, la SOI est, d’un point de vue astronomique, très petite, typiquement de quelques dizaines de parsecs. Par conséquent, il faut une très haute résolution spatiale pour étudier la SOI d’un SMBH et pouvoir adéquatement mesurer sa masse. C’est cette nécessité d’une haute résolution spatiale qui limite la mesure de masse de SMBHs à de plus grandes distances. Pour briser cette barrière, il nous faut donc trouver une manière d’améliorer la résolution spatiale d’objets observés à un plus au décalage vers le rouge. Le phénomène des lentilles gravitationnelles fortes survient lorsqu’une source lumineuse en arrière-plan se trouve alignée avec un objet massif en avant-plan, le long de la ligne de visée d’un observateur. Cette disposition a pour conséquence de distordre l’image observée de la source en arrière-plan. Puisque cette distorsion est inconnue et non-linéaire, l’analyse de la source devient nettement plus complexe. Cependant, ce phénomène a également pour effet d’étirer, d’agrandir et d’amplifier l’image de la source, permettant ainsi de reconstituer la source avec une résolution spatiale considérablement améliorée, compte tenu de sa distance initiale par rapport à l’observateur. L’objectif de ce projet consiste à développer une chaîne de simulations visant à étudier la faisabilité de la mesure de la masse d’un trou noir supermassif (SMBH) par cinéma- tique du gaz moléculaire à un décalage vers le rouge plus élevé, en utilisant l’apprentissage automatique pour tirer parti du grossissement généré par la distorsion d’une forte lentille gravitationnelle. Pour ce faire, nous générons de manière réaliste des observations du gaz moléculaire obtenues par le Grand Réseau d’Antennes Millimétrique/Submillimétrique de l’Atacama (ALMA). Ces données sont produites à partir de la suite de simulations hydrody- namiques Rétroaction dans des Environnements Réalistes (FIRE). Dans chaque simulation, l’effet cinématique du SMBH est intégré, en supposant le gaz moléculaire virialisé. Ensuite, le flux d’émission du gaz moléculaire est calculé en fonction de sa vitesse, température, densité, fraction de H2, décalage vers le rouge et taille dans le ciel. Le cube ALMA est généré en tenant compte de la résolution spatiale et spectrale, qui dépendent du nombre d’antennes, de leur configuration et du temps d’exposition. Finalement, l’effet de la forte lentille gravi- tationnelle est introduit par la rétro-propagation du faisceau lumineux en fonction du profil de masse de l’ellipsoïde isotherme singulière (SIE). L’exploitation de ces données ALMA simulées est testée dans le cadre d’un problème de régression directe. Nous entraînons un réseau de neurones à convolution (CNN) à apprendre à prédire la masse d’un SMBH à partir des données simulées, sans prendre en compte l’effet de la lentille. Le réseau prédit la masse du SMBH ainsi que son incertitude, en supposant une distribution a posteriori gaussienne. Les résultats sont convaincants : plus la masse du SMBH est grande, plus la prédiction du réseau est précise et exacte. Tout comme avec les méthodes conventionnelles, le réseau est uniquement capable de prédire la masse du SMBH tant que la résolution spatiale des données permet de résoudre la SOI. De plus, les cartes de saillance du réseau confirment que celui-ci utilise l’information contenue dans la SOI pour prédire la masse du SMBH. Dans les travaux à venir, l’effet des lentilles gravitationnelles fortes sera introduit dans les données pour évaluer s’il devient possible de mesurer la masse de ces mêmes SMBHs, mais à un décalage vers le rouge plus élevé. / Recent breakthroughs have been made in the study of supermassive black holes (SMBHs), thanks largely to the Event Horizon Telescope (EHT) team. However, determining the mass of these colossal entities at high redshifts remains a major challenge for astronomers. There are various direct and indirect methods for measuring the mass of SMBHs. The most accurate direct method involves resolving the kinematics of the molecular gas, a cold tracer, in the SMBH’s sphere of influence (SOI). The SOI is defined as the region where the gravitational potential of the SMBH dominates that of the host galaxy. However, since the mass of a SMBH is negligible compared to the mass of a galaxy, the SOI is, from an astronomical point of view, very small, typically a few tens of parsecs. As a result, very high spatial resolution is required to study the SOI of a SMBH and adequately measure its mass. It is this need for high spatial resolution that limits mass measurements of SMBHs at larger distances. To break this barrier, we need to find a way to improve the spatial resolution of objects observed at higher redshifts. The phenomenon of strong gravitational lensing occurs when a light source in the back- ground is aligned with a massive object in the foreground, along an observer’s line of sight. This arrangement distorts the observed image of the background source. Since this distor- tion is unknown and non-linear, analysis of the source becomes considerably more complex. However, this phenomenon also has the effect of stretching, enlarging and amplifying the image of the source, enabling the source to be reconstructed with considerably improved spatial resolution, given its initial distance from the observer. The aim of this project is to develop a chain of simulations to study the feasibility of measuring the mass of a supermassive black hole (SMBH) by kinematics of molecular gas at higher redshift, using machine learning to take advantage of the magnification generated by the distortion of a strong gravitational lens. To this end, we realistically generate observations of molecular gas obtained by the Atacama Large Millimeter/Submillimeter Antenna Array (ALMA). These data are generated from the Feedback in Realistic Environments (FIRE) suite of hydrodynamic simulations. In each simulation, the kinematic effect of the SMBH is integrated, assuming virialized molecular gas. Next, the emission flux of the molecular gas is calculated as a function of its velocity, temperature, density, H2 fraction, redshift and sky size. The ALMA cube is generated taking into account spatial and spectral resolution, which depend on the number of antennas, their configuration and exposure time. Finally, the effect of strong gravitational lensing is introduced by back-propagating the light beam according to the mass profile of the singular isothermal ellipsoid (SIE). The exploitation of these simulated ALMA data is tested in a direct regression problem. We train a convolution neural network (CNN) to learn to predict the mass of an SMBH from the simulated data, without taking into account the effect of the lens. The network predicts the mass of the SMBH as well as its uncertainty, assuming a Gaussian a posteriori distribution. The results are convincing: the greater the mass of the SMBH, the more precise and accurate the network’s prediction. As with conventional methods, the network is only able to predict the mass of the SMBH as long as the spatial resolution of the data allows the SOI to be resolved. Furthermore, the network’s saliency maps confirm that it uses the information contained in the SOI to predict the mass of the SMBH. In future work, the effect of strong gravitational lensing will be introduced into the data to assess whether it becomes possible to measure the mass of these same SMBHs, but at a higher redshift.