In this thesis, we construct a new hydrodynamical algorithm able of handling general compressible reactive flow problems, based on a finite-volume method inspired by the original MUSCL scheme of van Leer (1979). The algorithm is of second-order in the smooth part of the flow and avoids dimensional splitting. It uses MPI to achieve parallelism, and includes an astrophysical equation of state and a nuclear reaction network. It proves to be robust to tests cases. In particular it reproduces quite well the reactive and non-reactive results obtained with two different numerical methods (Fryxell & al. 1989, Busegnies & al. 2007). Moreover the time-dependent results are in agreement with the corresponding steady state solution. This gives us confidence in applying it to an astrophysical situation which has never been studied, the propagation of a detonation in conditions relevant to superbursts. The algorithm is described in (Noel & al. 2007). In a firt step we obtain the detonation profiles in pure carbon and in a mixture of carbon and iron. In both cases we underline the large difference between the total reaction length and the length on which some species burn. This difference leads to enormous numerical difficulties because all the length scales cannot be resolved at the same time in a single simulation. We show that the carbon detonation might be studied in a partial resolution approach like the one of Gamezo & al. (1999). In a second step we construct a new reduced nuclear reaction network able to reproduce the energy production due to the photo-disintegrations of heavy elements, like ruthenium, which are thought to occur during superbursts in mixed H/He accreting systems. Using this new nuclear network we simulate detonations in mixture of carbon and ruthenium. An interesting feature is that, in this case, all the reaction lengths can be resolved in the same simulation. This makes the C/Ru detonations easier to study in future multi-dimensional simulations than the pure carbon ones (Noel & al. 2007b). Finally we perform some numerical experiments which show that our algorithm is able to deal with initially inhomogeneous medium, and that the multi-dimensional simulations are attainable even if they are quite computational time consuming. - B. Van Leer, J. Comp. Phys., 21, 101, 1979 - Fryxell, B.A., Muller, E., and Arnett, W.D., Technical report MPA 449, 1989 - Busegnies, Y., Francois, J. and Paulus, G., Shock Waves, 11, 2007 - Gamezo, V.N., Wheeler, J.C., Khokhlov, A.M., and Oran, E.S., ApJ, 512, 827, 1999 - Noël, C., Busegnies, Y., Papalexandris, M.V. & al., A&A, 470, 653, 2007 - Noël, C., Goriely, S., Busegnies, Y. & Papalexandris, M.V., submitted to A&A, 2007b / Un algorithme parallèle basé sur une méthode aux volumes finis inspirée du schéma MUSCL de Van Leer (1979) a été construit. Il a été développé sur base de la méthode de Lappas & al. (1999) qui permet de résoudre simultanément toutes les dimensions spatiales. Cette méthode se base sur la construction de surfaces appropriées dans l'espace-temps, le long desquelles les équations de bilan se découplent en équations plus simples à intégrer. Cet algorithme est actuellement le seul à éviter le "splitting" des dimensions spatiales. Dans les modèles conventionnels (PPM, FCT, etc.), l'intégration spatiale des équations est réalisée de manière unidimensionnelle pour chaque direction. Un réseau de réactions nucléaires ainsi qu'une équation d'état astrophysique ont été inclus dans l'algorithme et celui-ci a ensuite été soumis à une grande variété de cas tests réactifs et non réactifs. Il a été comparé à d'autres codes généralement utilisés en astrophysique (Fryxell & al. 1989, Fryxell & al. 2000, Busegnies & al. 2007) et il reproduit correctement leurs résultats. L'algorithme est décrit dans Noël & al. (2007). Sur base de cet algorithme, les premières simulations de détonation dans des conditions thermodynamiques représentatives des Superbursts ont été réalisées. Différentes compositions du milieu ont été envisagées (carbone pur, mélange de carbone et de fer, mélange de carbone et de cendres du processus rp). Dans la plupart des systèmes où des Superbursts ont été observés, la matière accrétée est un mélange d'hydrogène et d'hélium. Dans ce cas, des phases de combustion précédant le Superburst produisent des nucléides plus lourd que le fer (Schatz & al. 2003). Ces nucléides peuvent être photodésintégrés durant le Superburst. Pour prendre en compte ces réactions endothermiques de photodésintégration, nous avons construit un nouveau réseau réduit de réactions nucléaires qui a été incorporé dans l'algorithme hydrodynamique (Noël & al. 2007b). Ce réseau réduit reproduit globalement l'énergétique d'un réseau complet et a permis de faire la première simulation numérique de détonation dans des conditions caractéristiques de systèmes accréteurs d'un mélange hydrogène-hélium. Finallement quelques simulations multidimensionelles préliminaires ont éte réalisées. - Busegnies, Y., Francois, J. and Paulus, G., Shock Waves, 11, 2007 - Fryxell, B.A., Muller, E., and Arnett, W.D., Technical report MPA 449, 1989 - Fryxell, B.A., Olson, K., Ricker, P. & al., ApJS, 131, 273, 2000 - Lappas, T., Leonard, A. and Dimotakis, P.E., SIAM J. Sci. Comput., 20, 1064, 1999 - Noël, C., Busegnies, Y., Papalexandris, M.V. & al., A&A, 470, 653, 2007 - Noël, C., Goriely, S., Busegnies, Y. & Papalexandris, M.V., submitted to A&A, 2007b - Röpke, F. K. PhD thesis, Technischen Universitat Munchen, 2003 - Schatz, H., Bildsten, L., Cumming, A. and Ouellette, M., Nuclear Physics A, 718, 247, 2003 - Van Leer, B. Comp. Phys., 21, 101, 1979 - Weinberg, N.N. and Bildsten, L., ArXiv e-prints, 0706.3062, 2007

In this thesis, we construct a new hydrodynamical algorithm able of handling general compressible reactive flow problems, based on a finite-volume method inspired by the original MUSCL scheme of van Leer (1979). The algorithm is of second-order in the smooth part of the flow and avoids dimensional splitting. It uses MPI to achieve parallelism, and includes an astrophysical equation of state and a nuclear reaction network. It proves to be robust to tests cases. In particular it reproduces quite well the reactive and non-reactive results obtained with two different numerical methods (Fryxell & al. 1989, Busegnies & al. 2007). Moreover the time-dependent results are in agreement with the corresponding steady state solution. This gives us confidence in applying it to an astrophysical situation which has never been studied, the propagation of a detonation in conditions relevant to superbursts. The algorithm is described in (Noel & al. 2007).<p><p>In a firt step we obtain the detonation profiles in pure carbon and in a mixture of carbon and iron. In both cases we underline the large difference between the total reaction length and the length on which some species burn. This difference leads to enormous numerical difficulties because all the length scales cannot be resolved at the same time in a single simulation. We show that the carbon detonation might be studied in a partial resolution approach like the one of Gamezo & al. (1999).<p><p>In a second step we construct a new reduced nuclear reaction network able to reproduce the energy production due to the photo-disintegrations of heavy elements, like ruthenium, which are thought to occur during superbursts in mixed H/He accreting systems. Using this new nuclear network we simulate detonations in mixture of carbon and ruthenium. An interesting feature is that, in this case, all the reaction lengths can be resolved in the same simulation. This makes the C/Ru detonations easier to study in future multi-dimensional simulations than the pure carbon ones (Noel & al. 2007b).<p><p>Finally we perform some numerical experiments which show that our algorithm is able to deal with initially inhomogeneous medium, and that the multi-dimensional simulations are attainable even if they are quite computational time consuming.<p><p>- B. Van Leer, J. Comp. Phys. 21, 101, 1979<p>- Fryxell, B.A. Muller, E. and Arnett, W.D. Technical report MPA 449, 1989<p>- Busegnies, Y. Francois, J. and Paulus, G. Shock Waves, 11, 2007<p>- Gamezo, V.N. Wheeler, J.C. Khokhlov, A.M. and Oran, E.S. ApJ, 512, 827, 1999<p>- Noël, C. Busegnies, Y. Papalexandris, M.V. & al. A&A, 470, 653, 2007<p>- Noël, C. Goriely, S. Busegnies, Y. & Papalexandris, M.V. submitted to A&A, 2007b<p><p>/<p><p>Un algorithme parallèle basé sur une méthode aux volumes finis inspirée du schéma MUSCL de Van Leer (1979) a été construit. Il a été développé sur base de la méthode de Lappas & al. (1999) qui permet de résoudre simultanément toutes les dimensions spatiales. Cette méthode se base sur la construction de surfaces appropriées dans l'espace-temps, le long desquelles les équations de bilan se découplent en équations plus simples à intégrer. Cet algorithme est actuellement le seul à éviter le "splitting" des dimensions spatiales. Dans les modèles conventionnels (PPM, FCT, etc.), l'intégration spatiale des équations est réalisée de manière unidimensionnelle pour chaque direction. <p>Un réseau de réactions nucléaires ainsi qu'une équation d'état astrophysique ont été inclus dans l'algorithme et celui-ci a ensuite été soumis à une grande variété de cas tests réactifs et non réactifs. Il a été comparé à d'autres codes généralement utilisés en astrophysique (Fryxell & al. 1989, Fryxell & al. 2000, Busegnies & al. 2007) et il reproduit correctement leurs résultats. L'algorithme est décrit dans Noël & al. (2007).<p><p>Sur base de cet algorithme, les premières simulations de détonation dans des conditions thermodynamiques représentatives des Superbursts ont été réalisées. Différentes compositions du milieu ont été envisagées (carbone pur, mélange de carbone et de fer, mélange de carbone et de cendres du processus rp). Dans la plupart des systèmes où des Superbursts ont été observés, la matière accrétée est un mélange d'hydrogène et d'hélium. Dans ce cas, des phases de combustion précédant le Superburst produisent des nucléides plus lourd que le fer (Schatz & al. 2003). Ces nucléides peuvent être photodésintégrés durant le Superburst. Pour prendre en compte ces réactions endothermiques de photodésintégration, nous avons construit un nouveau réseau réduit de réactions nucléaires qui a été incorporé dans l'algorithme hydrodynamique (Noël & al. 2007b). Ce réseau réduit reproduit globalement l'énergétique d'un réseau complet et a permis de faire la première simulation numérique de détonation dans des conditions caractéristiques de systèmes accréteurs d'un mélange hydrogène-hélium. <p>Finallement quelques simulations multidimensionelles préliminaires ont éte réalisées.<p><p>- Busegnies, Y. Francois, J. and Paulus, G. Shock Waves, 11, 2007<p>- Fryxell, B.A. Muller, E. and Arnett, W.D. Technical report MPA 449, 1989<p>- Fryxell, B.A. Olson, K. Ricker, P. & al. ApJS, 131, 273, 2000<p>- Lappas, T. Leonard, A. and Dimotakis, P.E. SIAM J. Sci. Comput. 20, 1064, 1999<p>- Noël, C. Busegnies, Y. Papalexandris, M.V. & al. A&A, 470, 653, 2007<p>- Noël, C. Goriely, S. Busegnies, Y. & Papalexandris, M.V. submitted to A&A, 2007b<p>- Röpke, F. K. PhD thesis, Technischen Universitat Munchen, 2003<p>- Schatz, H. Bildsten, L. Cumming, A. and Ouellette, M. Nuclear Physics A, 718, 247, 2003<p>- Van Leer, B. Comp. Phys. 21, 101, 1979<p>- Weinberg, N.N. and Bildsten, L. ArXiv e-prints, 0706.3062, 2007 / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

Cette thèse est consacrée à l'étude des oscillations périodiques détectées lors des sursauts X des étoiles à neutrons, dans des binaires X de faible masse. Ces oscillations offrent un moyen de sonder l'intérieur de ces objets, en mesurant notamment leur masse et leur rayon, pour ainsi contraindre l'équation d'état de la matière dense. J'ai développé des méthodes de détection et d'analyse de ces signaux, de leurs propriétés temporelles et leur dépendance en énergie. J'ai analysé les oscillations détectées dans tous les sursauts X de type 1 (ainsi qu'un super-sursaut) de 3 étoiles à neutrons observées avec l'instrument Rossi X-ray Timing Explorer/Proportional Counter Array . Sur les courbes de lumière des sursauts, j'ai sélectionné les segments donnant la meilleure signification statistique, pour construire un catalogue de profils moyens d'oscillations. La forme des profils varie grandement d'un sursaut à un autre, pour une même source. Un grand nombre de paramètres peuvent affecter les oscillations. J'ai élaboré un modèle de tache chaude à la surface de l'étoile en rotation rapide pour caractériser l'émission du sursaut X, dans un espace-temps relativiste. En utilisant les chaînes de Markov Monte Carlo pour explorer efficacement un espace des paramètres conséquent, les ajustements sur un échantillon de sursauts ont démontré l'applicabilité du modèle. Par contre, les contraintes obtenues sur la masse et le rayon de l'étoile sont limitées par la qualité des données de l'instrument utilisé. Enfin, des simulations révèlent que des mesures précises sur les paramètres sont possibles en augmentant la surface collectrice des détecteurs, comme le proposent les observatoires X du futur.

Ce manuscrit rassemble une série de travaux observationnels et phénoménologiques relatifs à des objets compacts variables du centre de notre Galaxie, à savoir le trou noir supermassif central, Sagittarius A*, et des étoiles à neutrons hébergées par des sursauteurs X. La première partie traite de la source Sgr A*, sujette à des éruptions quotidiennes, dont les mécanismes déclencheurs et les processus de rayonnement sont encore inconnus. Cette activité éruptive a été sondée par l'intermédiaire de plusieurs vastes campagnes d'observations multi-longueurs d'onde (en rayons gamma, rayons X, infrarouge et submillimétrique) étalées entre 2007 et 2009. Des données recueillies simultanément par les instruments XMM-Newton/EPIC, INTEGRAL/ISGRI+JEM-X, Fermi/LAT, VLT/NACO+VISIR et APEX/LABOCA, lors de plusieurs nouvelles éruptions majeures, ont ainsi permis de caractériser en détail le comportement spectro-temporelle de ces dernières et de contraindre les modèles d'émissions non-thermiques du milieu radiatif (synchrotron, Compton inverse, plasmoïde en expansion). Dans un second temps, une vingtaine de sursauts X de type I en provenance de deux binaires X de faible masse, transitoires, du noyau Galactique, GRS 1741.9-2853 et AX J1745.6-2901, ont été examinés à travers les données de différents satellites X de basse énergie (2-30 keV). Ces observations ont été discutées dans le cadre théorique, relativement bien établi, d'explosions thermonucléaires d'un mélange hydrogène-hélium, amassé à la surface d'étoiles à neutrons accrétantes.