Mémoire (M.Eng.)--Université du Québec à Chicoutimi, 1999. / Document électronique également accessible en format PDF. CaQCU

Thèse (M.Eng) -- Université du Québec à Chicoutimi, 2004. / Bibliogr.: f. 125-138. Document électronique également accessible en format PDF. CaQCU

Peu de stations au sol mesurent l'éclairement et des méthodes ont été développées pour l'évaluer à partir d'images prises par des satellites géostationnaires ou en orbite polaire. Dans les méthodes actuelles, dites inverses, le signal perçu par le satellite est inversé pour dériver le rayonnement au sol. Il est estimé de manière générale que les limites des méthodes inverses sont certainement atteintes en terme de précision. En conséquence, une nouvelle voie, basée sur la modélisation directe de la propagation du rayonnement solaire dans l'atmosphère, est explorée. Dans cette thèse, nous concevons une nouvelle méthode -dite Heliosat-4- basée sur le code de transfert radiatif libRadtran et tirant profit des produits de meilleure qualité, extraits des nouvelles capacités d'observation de la terre. Dans cet objectif, nous avons - par une analyse de sensibilité, établi les entrées nécessaires pour toute méthode basée sur un code de transfert radiatif, - établi que l'éclairement au sol pour une atmosphère nuageuse peut être considéré comme égal au produit de l'éclairement obtenu par ciel-clair et d'une fonction de l'extinction due aux nuages et de la contribution de la réflectance du sol, - mis en œuvre des approximations de l'équation de transfert radiatif pour des atmosphères claire ou nuageuse, permettant des calculs rapides. Heliosat-4 est affranchi des paramètres empiriques, facile à implémenter et fournit simultanément les composantes directe, diffuse et la distribution spectrale de l'éclairement. Ces caractéristiques, bien que nécessaires pour différentes applications, sont rarement fournies. Les premiers résultats de la validation de cette nouvelle méthode avec les mesures sol sont prometteurs : les performances observées sont proches de celles des méthodes actuelles.

L'objectif industriel de cette thèse est de proposer un modèle de transfert radiatif dans un réacteur de reformage de méthane. Dans ce procédé, des gaz réactifs circulent dans le réacteur tubulaire rempli de pastilles catalytiques.L'empilement de pastilles constitue un milieu poreux où le rapport de la taille caractéristique des pores sur la dimension radiale du réacteur est grand (1=10 à 1=5). De plus, les très forts gradients de porosité dus à l'organisation des pastilles au voisinage des parois ont un impact important sur les transferts thermiques et en particulier les transferts radiatifs.L'objectif scientifique est de développer et valider un modèle de transfert radiatif applicable à des milieux poreux fortement hétérogènes et anisotropes ne suivant pas la loi de Beer. Dans un premier temps, les propriétés radiatives du milieu homogénéisé équivalent au milieu poreux réel sont complètement déterminées par la fonction de distribution cumulée d'extinction Gext, la fonction de phase p et la porosité Π. Ces fonctions, précédemment introduites pour des milieux homogènes éventuellement anisotropes, sont calculées avec une grande précision par une méthode de Monte Carlo. Elles ont été généralisées ici à des milieux hétérogènes. Il a été montré à partir d'un nouveau critère de validité adapté aux milieux hétérogènes que le milieu homogénéisé équivalent ne suit pas la loi de Beer, en particulier au voisinage des parois.De ce fait, l'équation de transfert radiatif généralisée (GRTE) doit prendre en compte l'émission par un milieu non Beerien fortement hétérogène même à la limite optiquement mince : un coefficient d'absorption n'y a pas de sens physique et des corrélations entre émission et transmission apparaissent dues au caractère non Beerien. Le principe de réciprocité et les propriétés des fonctions d'extinction Gext ont permis d'exprimer rigoureusement les termes sources d'émission dans ce type de milieux fortement hétérogènes non Beeriens. Un facteur de corrélation émission-transmission a été introduit. La GRTE, sous forme intégrale, a été résolue par une méthode de transfert de Monte Carlo. Le modèle complet a été appliqué après validation aux réacteurs de reformage de méthane de Air Liquide. / The industrial goal of this work is to propose a radiative transfer model in a tubular reactor of steam methane reforming. During the reforming process, reactive gases are injected in the tubular reactor filled with catalytic pellets. The packed bed of pellets forms a porous medium, and a particular feature of it is that the characteristic pore size is large compared to the reactor inner dimension. In addition, the organization of the pellets in the near-wall region results in important porosity gradients which have a significant effect on the heat transfer, and more specifically on the radiative transfer.The scientific goal is to develop and validate a radiative transfer model applicable to strongly nonhomogeneous, anisotropic and non Beerian porous media.First, the radiative properties of the homogenised phase equivalent to the real porous medium are completely determined by the cumulated distribution function of extinction Gext, the phase function p, and the local porosity Π. These functions, previously introduced for statistically homogeneous and anisotropic porous media, are calculated very accurately by a Monte Carlo method. They have been extended to statistically non-homogeneous porous media. Similarly, the expression of the validity criterion of the Beer law is extended to statistically anisotropic and non-homogeneous porous media : it is proven that for the considered porous media the Beer law is not valid in the homogenised phase, in particular in the vicinity of the walls. As a result, the Generalized Radiative Transfer Equation (GRTE) is needed and the emission source terms must be determined in a strongly nonhomogeneous non Beerian even at the optically thin limit : an absorption coefficient doesn't have any physical meaning and correlations between emission and transmission appear due to the non-Beerian behavior.The reciprocity principle and the properties of the extinction functions Gext allow the emission source terms in this kind of strongly non-homogeneous and non-Beerian media to be accurately determined. A correlation factor emission-transmission has been introduced. The GRTE has been solved by a Monte Carlo method.The complete model is applied, after validation, to the steam methane reformers in use by Air Liquide.

Les phénomènes de diffraction et de diffusion des ondes jouent un rôle important dans l'interprétation de la coda des sismogrammes. Par conséquent, une compréhension approfondie des mécanismes de diffraction et de leurs influences sur la propagation des ondes est une étape fondamentale vers l'identification des propriétés statistiques d'un milieu aléatoire. Cette thèse porte sur la diffraction des ondes élastiques dans des milieux aléatoirement hétérogènes avec un comportement local isotrope. On s'intéresse au régime où: La longueur d'onde est du même ordre de grandeur que la longueur de corrélation, la longueur d'onde est petite comparé à la distance de propagation (haute-fréquence) et l'amplitude des fluctuations est petite. Une approche cinétique basée sur les équations de transfert radiatif des ondes élastiques est adoptée. La première partie de cette thèse décrit une analyse détaillée de l'influence de la structure de corrélation sur les paramètres de diffraction et sur l'établissement d'un régime de diffusion. La seconde partie présente les simulations éléments spectraux à grande échelle des ondes élastiques afin d'observer numériquement l'apparition d'un régime d'équipartition. Des analyses théoriques ainsi que des simulations montrent également une nouvelle approche pour l'identification des propriétés statistiques du milieu. / Scattering and diffusion phenomena play a crucial role in the interpretation of the coda part ofseismograms. Consequently, a profound understanding of scattering mechanisms and their effectson wave propagation is a fundamental step towards the identification of the statistical propertiesof random media. The focus of this work is on the scattering of elastic waves in a randomly heterogeneousmedia with locally isotropic material behavior. The weakly heterogeneous regime isconsidered, in which the wave length is similar to the correlation length, the wave length is smallcompared to the propagation length (high frequency) and the amplitude of the heterogeneities issmall. A kinetic framework based on the transport equations of elastic waves is adopted. Thefirst part of the thesis describes a detailed analysis of the influence of the correlation structure onthe scattering parameters and on the arising of the diffusion regime. The second part presentslarge scale spectral element simulations of elastic waves to observe numerically the onset of theequipartitioning regime. The theoretical analyses and simulations also reveal a novel approach toidentify local properties of the heterogeneous medium.

L'atmosphère de Saturne subit d'importantes variations saisonnières d'insolation, à cause de son obliquité, de son excentricité et de l'ombre de ses anneaux. Dans la stratosphère (de 20 hPa à 10-4 hPa), les échelles de temps photochimiques et radiatives sont du même ordre de grandeur que la période de révolution de Saturne (29,5 ans). On s'attend donc à mesurer des variations saisonnières et méridiennes significatives de la température et des espèces produites par la photochimie (en particulier C2H6, C2H2 et C3H8) dans cette région. Grâce à sa durée (2004-2017), la mission Cassini est l'occasion inédite de suivre l'évolution saisonnière de l'atmosphère de Saturne.Au cours de ma thèse, j'ai analysé des observations au limbe Cassini/CIRS car elles permettent de sonder à la fois la structure méridienne et verticale de la stratosphère de Saturne. Ainsi, j'ai mesuré les variations saisonnières de la température et des abondances de C2H6, C2H2 et C3H8. J'ai également contribué au développement d'un modèle radiatif-convectif et d'un GCM (Global Climate Model) de l'atmosphère de Saturne. Les prédictions de ces modèles sont comparées avec les températures mesurées avec CIRS, de façon à étudier les processus radiatifs et dynamiques qui contribuent à l'évolution saisonnière. Les simulations numériques réalisées avec ce GCM m'ont également permis d'étudier la propagation des ondes atmosphérique ainsi que les effets de l'ombre des anneaux sur l'atmosphère de Saturne. Par ailleurs, la comparaison entre les distributions de C2H6, C2H2 et C3H8 et des modèles photochimiques (Moses et Greathouse 2005, Hue et al. 2015) donne des indications sur le transport méridien. / Saturn's atmosphere undergoes important seasonal variations of insolation, due to its obliquity, its eccentricity and the shadow of its rings. In the stratosphere (from 20 hPa to 10-4 hPa), radiative and photochemical timescales are in the same order as Saturn's revolution period (29.5 ans). Hence, significative seasonal and meridional variations of temperature and photochemical by-products (especially C2H6, C2H2, and C3H8) are expected. Because of its duration (2004-2017), the Cassini mission is an unprecedented opportunity to monitor the seasonal evolution of Saturn's atmosphere. During my PhD, I analysed Cassini/CIRS limb observations as they probe the meridional and vertical structure of Saturn's stratosphere. Hence, I measured seasonal variations of temperature and abundances of C2H6, C2H2, and C3H8. I also contributed to the development of a radiative-convective model and a GCM (Global Climate Model) of Saturn's atmosphere. The predictions of these models are compared with the temperatures measured from CIRS observations, in order to study the radiative and dynamical processes which contribute to the seasonal evolution. Numerical simulations performed with the GCM also allowed me to study atmospheric waves propagation and the effects of rings shadowing in Saturn's atmosphere. Besides, comparison between C2H6, C2H2, and C3H8 distributions and photochemical models (Moses and Greathouse 2005, Hue et al., 2015) give insights on meridional transport.

Les travaux présentés dans ce manuscrit sont consacrés à l'étudie du un petit instrument sophistiqué de mesure d'épaisseur optique : ODS (Optical Depth Sensor) a été conçu pour étudier l'atmosphère martienne et terrestre. L'une de ses principales missions consiste à fournir des mesures journalières de l'épaisseur optique des aérosols (AOD: Aerosol Optical Depth) et de détecter et caractériser les nuages optiquement minces à l'aube et au crépuscule. Les méthodes d'analyses sont basées sur l'utilisation de tables d'intensité reproduisant les signaux observés ODS en fonction de différents paramètres clés. Ces tables ont été réalisées à l'aide d'un code de transfert radiatif dans une géométrie en plans parallèles pour l'estimation de l'opacité des poussières, tandis que les propriétés des nuages sont issues d'un modèle de type Monte-Carlo en géométrie sphérique. Le premier objectif de ce travail a consisté à développer les algorithmes d'analyses nécessaires à l'étude du signal d'ODS. / The work presented in this manuscript is devoted to the development of reliable retrieval procedures for a lightweight and sophisticated optical depth sensor (ODS) which measures alternatively scattered flux at zenith and the sum of the direct flux and the scattered flux in the blue and red wavelength ranges. The ODS sensor is dedicated to Mars and Earth atmosphere and its principal goals are perform measurements of the daily mean aerosol optical depth (AOD) and retrieve the altitude and optical depth of high altitude clouds at twilight. The retrieval procedure is based on the use of look-up tables of intensities reproducing the signals that should be observed by ODS, as a function of different key parameters. For the estimation of AOD, the look-up tables are obtained by using a plane-parallel radiative transfer code, while for the study of cloud properties a Monte-Carlo radiative transfer code in spherical geometry is used.

Les protubérances solaires sont de fascinantes structures magnétiques à grande échelle dans l'atmosphère solaire. Elles sont l'objet de recherches depuis des dizaines d'années mais des questions comme leur formation, leur stabilité et leur éruption ne sont toujours pas bien comprises. De grands progrès ont été obtenus dans ce domaine et celui des éruptions en particulier avec une combinaison d'observations synoptiques et continues depuis l'espace (voir SoHO/EIT, STEREO/SECCHI/EUVI, et SDO/AIA) et en spectro-imagerie multilongueurs d'onde. Depuis le lancement du satellite IRIS en 2013, beaucoup de résultats d'observations et de modélisations ont été obtenus grâce à la très haute résolution spectrale et spatiale d'IRIS tant en spectroscopie qu'en imagerie. Dans cette thèse, nous nous focalisons sur les signatures observationnelles des processus mis en avant pour expliquer les éruptions de protubérances. Nous essayons aussi de déterminer les variations de conditions physiques d'une protubérance éruptive et d'estimer à la fois les masses de matière qui quittent le Soleil et celles qui retournent au Soleil pendant une éruption. En ce qui concerne les vitesses, nous parvenons à construire le vecteur vitesse en chaque pixel de la protubérance observée, grâce à la combinaison d'une méthode de "flot optique" appliquée aux images AIA 304 et IRIS Mg II h&k ce qui donne accès aux vitesses dans le plan du ciel, d'une part, et grâce à la technique Doppler appliquée aux profils IRIS Mg II h&k ce qui permet de calculer les vitesses le long de la ligne de visée, d'autre part. En ce qui concerne les densités et températures, nous comparons les intensités absolues observées avec les valeurs déduites de calculs de transfert radiatif Hors-Equilibre Thermodynamique Local, ce qui nous permet de construire des cartes de ces quantités. Les densités électroniques ainsi déduites sont dans la gamme 1.3E9 - 7.0E10 cm⁻³ et les densités totales d'hydrogène dans la gamme 1.5E9 – 3.6E11 cm⁻³ dans les diverses régions de la protubérance. La température moyenne est d'environ 1.1E4 K une valeur plus élevée que dans les protubérances quiescentes. Le degré d'ionisation varie entre 0.1 et 10. Les flux de masse à l'intérieur de la protubérance et leurs variations temporelles ont pu ainsi être calculés. La masse totale ainsi déduite est dans la gamme 1.3E14 - 3.2E14 g. La perte totale de masse de la protubérance vers la surface solaire évaluée sur la durée d'observation avec IRIS est d'un ordre de grandeur plus faible que la masse totale de la protubérance. Nous explorons aussi les corrélations entre indices spectraux observables dans les raies h et k de Mg II et des quantités physiques comme la densité et la mesure d'émission (ME). Nous avons choisi de calculer des modèles uni-dimensionnels (1D) isothermes et isobares en utilisant le code PRODOP_Mg NLTE disponible à MEDOC (IAS) et en procédant au calcul exact du rayonnement incident. Nous en déduisons des corrélations entre les intensités émergentes dans les raies h et k d'une part, et les densités et les ME d'autre part. Moyennant quelques hypothèses sur la température, nous établissons une relation entre les intensités k (et h) et la ME, une relation qui devrait être utile pour déterminer les densités d'hydrogène et d'électrons aussi bien que l'épaisseur effective d'une protubérance observée. Ainsi donc, l'évolution des propriétés morphologiques et thermodynamiques d'une protubérance éruptive ont été étudiées de dans cette thèse. Ces travaux conduisent à une meilleure compréhension de quelques aspects des protubérances (éruptives), tels la distribution et l'évolution des densités, de la température, des vitesses et du degré d'ionisation. Nous avons ainsi fourni des contraintes utiles à la modélisation des protubérances. Dans notre conclusion, nous résumons nos résultats et proposons quelques suggestions pour de futures analyses et observations et pour des capacités instrumentales optimisées. / Solar prominences are fascinating, large-scale magnetic structures in the solar atmosphere. They have been investigated for many decades, but the issues of their formation, stability, and eruption are still not well understood. Much progress has been made in our knowledge of prominences and their eruptions with both synoptic measurements from space (with SoHO/EIT, STEREO/SECCHI/EUVI, and SDO/AIA) and multiwavelength spectro-imaging. Since the launch of IRIS in 2013, a lot of results have been obtained in both observational and modeling domains with IRIS high spectral and spatial resolution imaging and spectroscopy. In this thesis, we focus on the observational signatures of the processes which have been put forward for explaining eruptive prominences. We also try to figure out the variations of physical conditions of the eruptive prominence and estimate the masses leaving the Sun vs. the masses returning to the Sun during the eruption. As far as velocities are concerned, we derive a full velocity vector for each pixel of the observed prominence by combining an optical flow method on the AIA 304Å and IRIS Mg II h&k images in order to derive the plane-of-sky velocities and a Doppler technique on the IRIS Mg II h&k profiles to compute the line-of-sight velocities. As far as densities and temperatures are concerned, we compare the absolute observed intensities with values derived from Non-Local Thermodynamic Equilibrium (NLTE) radiative transfer computations to build maps of these quantities. The derived electron densities range from 1.3E9 to 7.0E10 cm⁻³ and the derived total Hydrogen densities range from 1.5E9 to 3.6E11 cm⁻³ in different regions of the prominence. The mean temperature is around 1.1E4 K which is higher than in quiescent prominences. The ionization degree is in the range of 0.1 to 10. The mass flows in the prominence and their variations with time are consequently computed. The total mass is 1.3E14 to 3.2E14 g. The total mass drainage from the prominence to the solar surface during the observation of IRIS is about one order of magnitude smaller than the total mass of prominence. We also explore the correlations between the observable spectral features in h and k lines of Mg II to physical quantities such as the density and the Emission Measure (EM). We choose to compute one-dimensional (1D) isothermal and isobaric models using the PRODOP_Mg NLTE code available at MEDOC (IAS) with the exact computation of the incident radiation. Then we derive correlations between the k and h emergent intensities on one hand and the densities and EM on the other hand. With some assumptions on the temperature, we obtain a unique relation between the k (and h) intensities and the EM that should be useful for deriving either the hydrogen and electron densities or the effective thickness of an observed prominence. Thus, the evolution of the morphology and thermodynamic properties of an erupting prominence have been studied in the thesis. These investigations lead to our understanding in some aspects of prominences, e.g., the distribution and evolution of densities, temperatures, velocities and ionization degree. These could be useful constraints for theoretical prominence models. In the conclusion, we summarize our results and provide some suggestions for future analysis, observations and ideal observing capabilities.

Les recherches de SNe ont commencé il y a plus de 100 ans. Depuis, il a été possible de collecter beaucoup de données d'observations astronomiques. Les astronomes ont développé une classification détaillée et ont abouti un relatif consensus sur la nature physique de ces événements très différents. Néanmoins, beaucoup de questions restent sans réponse. En résumé, les supernovæ de type II (riche en hydrogène) résultent de l'éjection l'enveloppe des supergéantes rouges (SGR). Les principales sources de connaissance sur ces objets sont l'évolution de leur luminosité en fonction du temps (`courbes de lumière') et leurs spectres observés à différentes époques. La méthode la plus utilisée pour extraire les informations des données d'observation est la modélisation des courbes de lumières et des spectres des supernovæ. Dans le Chapitre 1 (Introduction), nous présentons successivement l’évolution stellaire, la physique des explosions et l’évolution des éjectas. Nous décrivons aussi les différents types de supernova ; l’état actuel des connaissances sur les CCSNe ainsi que ces limitations. Nous discutons de la théorie de l'évolution stellaire. Nous décrivons notre approche numérique au Chapitre 2 (Supernova modelling). Elle consiste en trois étapes principales : la modélisation de l'évolution stellaire, l'explosion de l'étoile SGR résultante, et la modélisation de l'évolution des éjectas. Nous présentons la structure modélisée des étoiles SGR ; ces modèles et techniques de calcul sont similaires aux modèles utilisés dans les chapitres suivants. Nous discutons notre méthode d'explosion d'un progéniteur quand son noyau dégénéré commence à s'effondrer. Dans le Chapitre 3 (Observational properties), nous discutons les propriétés observées en photométrie et spectroscopie des CCSNe. Nous extrayons les propriétés statistiques de l'échantillon existant. En utilisant la technique présentée, nous avons effectué une étude détaillée de SN 2008bk, une supernova bien observée (Chapitre 4). Nous pouvons contraindre les propriétés du progéniteur et des paramètres d'explosion. Notre modélisation nous permet de comparer non seulement les propriétés de base telles que la luminosité, mais aussi à analyser en détail les caractéristiques spectrales, telles que la présence de certaines raies spectrales ainsi que leur morphologie. Nous montrons qu'une étoile de 12M⊙ sur la séquence principale est un bon candidat au progéniteur de SN 2008bk. Aussi, nous discutons de la forme asymétrique de la raie Hα et concluons qu'elle provient du chevauchement avec la raie forte du Ba II 6596.9 Å. SN 2008bk, avec quelques dizaines d'autres objets, forme une sous-classe importante de CCSNe — supernovae à faible luminosité. Nous avons consacré une attention particulière à cette classe d'objets, dont l'uniformité et les données d'observation nous permet de tirer des conclusions importantes. Au Chapitre 5, nous étudions l'échantillon de 17 SNe de faible luminosité et analysons la possibilité que ces événements résultent d'explosions de progéniteurs de petite ou de grande masse. Il n'y a pas d'accord solide dans la communauté astronomique sur les progéniteurs possibles des explosions de SNe à faible luminosité. Notre analyse montre que les étoiles massives de masse inférieure (~12 M⊙) sont de bons candidats pour les progéniteurs de cette sous-classe de SNe. De plus, nos simulations d'étoiles de masse élevée (25 et 27 M⊙) montrent qu'une explosion ayant une luminosité aussi faible aurait des propriétés d’observation remarquables qui ne sont pas présentes dans les données. Dans le Chapitre 6, nous étendons notre étude sur toutes les CCSNe, en utilisant des modèles plus énergétiques que dans les Chapitres 4 et 5. Nous fournissons des preuves que ce qui différencie la diversité de SNe II est l'énergie d'explosion plutôt que la masse des éjectas (ou plus précisément la masse de l'enveloppe riche en H de progéniteur). / Dedicated SN searches started over 100 years ago. Over that time, astronomers have collected large sets of observational data. They have developed detailed classification and achieved general agreement on the nature of these events. Nevertheless, a lot of questions remain unanswered. In short, most Type II SNe (hydrogen-rich SNe) are terminal explosions of red supergiant (RSG) stars. The main source of knowledge about these objects are the way their luminosity changes with time (`light curves') and how their radiation is distributed in wavelength. One of the widely used methods to extract the information from the observational data is computer modelling. The largest part of our work lays in the numerical simulations. In Chapter 1 (Introduction), we present succinctly the necessary theory which includes stellar evolution, explosion physics and ejecta evolution. We discuss different types of SNe; the modern knowledge on CCSNe and its problems. We discuss stellar evolution theory. We describe the nucleosynthesis that takes place in the cores of massive stars and gives rise to their final chemical stratification. We describe our numerical approach in Chapter 2 (Supernova modelling). It includes three major steps: stellar evolution modelling, explosion of the resulting RSG star, and ejecta evolution modelling. We present modelled structure of RSG stars; these models and computational techniques are similar to models used in subsequent chapters. We then discuss our numerical methods of exploding a SN once its degenerate core starts collapsing. We discuss explosive nucleosynthesis and its impact on the progenitor composition, production of unstable isotopes and the basic physics of radioactive decay. In Chapter 3 (Observational properties), we discuss the photometric and spectral observational properties of core-collapse SNe. We extract statistical properties of the existing sample. Using the presented technique, we performed a detailed study of the well observed object SN 2008bk (Chapter 4). We are able to constrain its progenitor and explosion properties. Our modelling allows us to compare not only the basic properties such as luminosity, but also to analyze in detail the spectral features, such as line identification and morphology. We show that a 12 M⊙ star on the main sequence is a good candidate for the progenitor of SN 2008bk. Also we discuss the asymmetric shape of the Hα line and conclude that it stems from the overlap with the strong Ba II 6596.9 Å line. SN 2008bk, together with about 20 objects, form a subclass of low-luminosity CCSNe Type II. We devoted a particular attention to this class of objects, whose uniformity and observational data allows us to draw important conclusions. In Chapter 5 (Low-luminosity Type II-P SNe), we study the sample of 17 low-luminosity SNe and analyze the possibility that these events are the result of explosions of low- and high-mass progenitors. There is no solid agreement in the astronomical community on the possible progenitors of the low-luminosity explosions of Type II SNe. Our analysis shows that lower-mass massive stars (~12 M⊙) are good candidates for the progenitors of this subclass of SNe. Moreover, our simulations of high-mass stars (25 and 27 M⊙) show that such low brightness of the explosion of such a massive object would have notable observational properties which are not present in the data. In Chapter 6 (Kinetic energy variation), we extend our study further on the whole class of hydrogen-rich core-collapse SN, using more energetic models than in Chapters 4 and 5. We provide evidences that what differentiates the diversity of SNe II is the explosion energy rather than the ejecta mass (or H-rich envelope mass of the progenitor).

La résolution numérique directe des problèmes de transport de photons en interaction avec un milieu matériel est très coûteuse en mémoire et temps CPU. Pour pallier ce problème, une méthode consiste à construire des modèles réduits dont la résolution est moins coûteuse. La littérature abonde de ce genre de modèles : modèles probabilistes (Monte-Carlo), modèles aux moments (M₁, PN), modèles aux ordonnées discrètes (SN), modèles de diffusion... Dans cette thèse, nous nous intéressons aux modèles PN dans lesquels l'opérateur de transport est approché par projections sur une base tronquée d'harmoniques sphériques. Ces modèles ont l'avantage d'être arbitrairement précis sur la dimension angulaire et ne présentent pas les défauts connus des autres méthodes (bruit stochastique, "effets de raies") pouvant briser les éventuelles symétries du problème. Ce dernier point est capital pour la simulation d'expériences de fusion par confinement inertiel (FCI) où la symétrie sphérique joue un rôle important dans la précision des résultats. Nous étudions donc dans cette thèse la structure mathématique des modèles PN ainsi que leur discrétisation dans le cas d'une géométrie 1D sphérique.Nous commençons par le cas du transport linéaire dans le vide. Même dans ce cas simple, les équations du modèle PN contiennent des termes sources d'origine géométrique dont la discrétisation s'avère délicate. Jusqu'à présent, les différents schémas utilisés étaient insatisfaisants pour les raisons suivantes : (1) mauvais comportement au voisinage de r = 0 (phénomène de "flux-dip"), (2) non préservation des équilibres stationnaires, (3) pas de preuve formelle de stabilité. À la lumière de récents travaux, nous proposons une nouvelle discrétisation qui capture exactement les états d'équilibres. Nous démontrons en particulier la stabilité en norme L² du schéma. Nous étendons par la suite ce schéma au cas du transport de photons dans un milieu matériel figé et nous nous intéressons au comportement du schéma en limite diffusion (propriété "asymptotic-preserving").Dans un second temps, nous nous intéressons au couplage entre rayonnement et hydrodynamique. Devant l'absence de consensus sur les modèles "transport" d'hydrodynamique radiative issus de la littérature, nous établissons une étude comparative de ceux-ci basée sur leurs propriétés mathématiques. Nous nous intéressons particulièrement aux propriétés suivantes : (1) conservation de l'énergie et de l'impulsion, (2) précision des effets comobiles, (3) existence d'une entropie mathématiques compatible et (4) restitution de la limite diffusion. Notre étude se réduit aux modèles dits "mixed-frame" et une attention particulière est toujours portée sur l'approximation "PN" de l'opérateur de transport. Nous identifions des défauts (conservation ou entropie) sur des modèles existants et proposons une correction entropique conduisant à un modèle PN satisfaisant toutes les propriétés mathématiques listées ci-dessus. / Computational costs for direct numerical simulations of photon transport problemsare very high in terms of CPU time and memory. One way to tackle this issue is todevelop reduced models that a cheaper to solve numerically. There exists number of these models : moments models, discrete ordinates models (SN), diffusion-like models... In this thesis, we focus on PN models in which the transport operator is approached by mean of a truncated development on the spherical harmonics basis. These models are arbitrary accurate in the angular dimension and are rotationnaly invariants (in multiple space dimensions). The latter point is fundamental when one wants to simulate inertial confinment fusion (ICF) experiments where the spherical symmetry plays an important part in the accuracy of the numerical solutions. We study the mathematical structure of the PN models and construct a new numerical method in the special case of a one dimensionnal space dimension with spherical symmetry photon transport problems. We first focus on a linear transport problem in the vacuum. Even in this simple case, it appears in the PN equations geometrical source terms that are stiff in the neighborhood of r = 0 and thus hard to discretise. Existing numerical methods are not satisfactory for multiple reasons : (1) unaccuracy in the neighborhood of r = 0 ("flux-dip"), (2) do not capture steady states (well-balanced scheme), (3) no stability proof. Following recent works, we develop a new well-balanced scheme for which we show the L² stability. We then extend the scheme for photon transport problems within a no moving media, the linear Boltzmann equation, and interest ourselves on its behavior in the diffusion limit (asymptotic-preserving property). In a second part, we consider radiation hydrodynamics problems. Since modelisation of these problems is still under discussion in the litterature, we compare a set of existing models by mean of mathematical analysis and establish a hierarchy. For each model, we focus on the following mathematical properties : (1) energy and impulsion conservation, (2) accuracy of the comobile effects, (3) existence of a mathematical entropy and (4) behavior in the diffusion limit. Our study reduces to « laboratory frame » models and we are still interested in the PN approximation of the transport operator. We identify defects in entropy structure of existing models and propose an entroy correction which leads to PN-based radiation hydrodynamics models which satisfy all the properties listed above.