Nous cherchons à mettre en évidence l'influence des hétérogénéités de température, de structure ou de composition des atmosphères sur leurs observations. Dans le années à venir, de plus en plus d'appareils vont permettre l'observation par transmission des atmosphères d'exoplanètes. Toutefois, les outils numériques permettant de contraindre ces dernières reposent sur des modèles simples à une dimension. Ils supposent en effet des atmosphères ne possédant qu'une structure verticale (le climat est le même en tout point de la surface, la composition ou la température n'évolue qu'avec l'altitude). Cette approche a le mérite de permettre des calculs rapides et de contraindre les paramètres globaux de l'atmosphère avec des temps raisonnables. Ceci ne serait pas possible en l'état avec une modélisation en 3 dimensions des atmosphères, même si ce serait beaucoup plus réaliste. Ce que nous cherchons à mettre en évidence, ce sont les limites des techniques actuelles d'inversion et donc, de caractérisation des atmosphères qui seront observées. Pour cela, il fallait mettre au point un logiciel capable de résoudre le transfert radiatif au sein d'une atmosphère en 3 dimensions (et non plus 1 seul). Une fois le logiciel terminé, nous avons éprouvé l'algorithme de traitement du signal TauREx en comparant les résultats qu'il proposait à des simulations atmosphériques parfaitement contrôlées. Nous nous sommes tout principalement arrêté sur les biais découlant d'hétérogénéités de température en simulant des atmosphères avec un fort contraste jour/nuit. Ceci nous a permis de caractériser les biais découlant de ces types d'hétérogénéités, de les quantifier et de mettre l'accent sur un biais jusqu'ici très sous-estimé par la communauté, à savoir celui découlant des hétérogénéités le long de la ligne de visée. Nous avons appuyé nos propos et concentré nos efforts sur l'interprétation de l'inversion d'une simulation complexe de l'atmosphère de GJ 1214 b. La reconstitution de la chaine observationnelle : GCM (LMD), Pytmosph3R (LAB) et TauREx (UCL) ouvre les portes d'un vaste panel d'études envisageables, et notamment tout ce qui va concerner l'identification et la caractérisation des biais systématiques qui incomberont les observations à venir. / Transmission spectroscopy provides us with information on the atmospheric properties at the limb, which is often intuitively assumed to be a narrow annulus aound ther planet. Consequently, the few recent studies on the effect of atmospheric horizontal heterogeneities on transmission spectra have used approaches sensitive to variations along the limb only. Here we demonstrate that the region probed in transmission – the limb – actually extends significantly toward the day and night sides of the planet. Consequently we show that thestrong day-night thermal and compositional gradients expected on synchronous exoplanets create sufficient heterogeneities across the limb to result in important systematic effects on the spectrum and bias its interpretation. To quantify these effects, we developed a 3D radiative transfer model able to generate transmission spectra of atmospheres based on 3D atmospheric structures, whether they come from a Global Climate Model or more parametrized models. We first apply this tool to a simulation of the atmosphere of GJ 1214 b toproduce synethic JWST observations and show that producing a spectrum using only atmospheric columns at the terminator results in errors greater than expected noise. This demonstrates the necessity of a real 3D approach to model data for such precise observatories.Second, we investigate how day-night temperature gradients cause a systematic bias in retrieval analysis performed with 1D forward models. For that purpose we synthesize a large set of forward spectra for prototypical HD209458 b and GJ 1214 b type planets varying the temperatures of the day and night sides as well as the width of the transition region. We then perform typical retrievalanalyses and compare the retrieved parameters to the ground truth of the input model. This study reveals systematic biases on the retrieved temperature (found to be higher than the terminator temperature) and absorber abundances. This is due to the fact that the hotter dayside is more extended vertically and screens the nightside—a result of the nonlinear properties of atmospheric transmission.These biases will be difficult to detect as the 1D profiles used in the retrieval procedure are found to provide an excellent match to theobserved spectra based on standard fitting criteria (chi2, posterior distributions). This fact needs to be kept in mind when interpretingcurrent and future data.

L’imagerie moléculaire occupe une place majeure dans le domaine de la recherche préclinique. Parmi les modalités existantes, les techniques optiques fondées sur la détection d’un rayonnement visible ou proche infrarouge sont les plus récentes et sont principalement représentées par les méthodes de tomographie optique de luminescence. Ces méthodes permettent une caractérisation 3D d’un milieu biologique par la reconstruction de cartes de concentration ou la localisation de marqueurs luminescents sensibles à des processus biologiques et chimiques se déroulant à l’échelle de la cellule ou de la molécule. La tomographie de luminescence se fonde sur un modèle de propagation de la lumière dans les tissus, un protocole d’acquisition du signal en surface du milieu et une procédure numérique d’inversion de ces mesures afin de reconstruire les paramètres d’intérêts. Ce travail de thèse s’articule donc autour de ces trois axes et apporte un élément de réponse à chacun des problèmes. L’objectif principal de cette étude est d’introduire et de présenter des outils d’évaluation des performances théoriques d’une méthode de tomographie optique. L’un des aboutissements majeurs est la réalisation de reconstructions tomographiques expérimentales à partir d’images acquises par un imageur optique conçu pour l’imagerie planaire 2D et développé par la société Quidd. Dans un premier temps nous abordons la théorie du transport en milieu diffusant afin de poser les concepts et outils sur lesquels vont s’appuyer l’ensemble des travaux. Nous présentons particulièrement deux modèles de propagation différents ainsi que les méthodes de résolution et les difficultés théoriques qui leur sont liées. Dans une deuxième partie nous introduisons les outils statistiques utilisés pour caractériser les systèmes tomographiques et leur résolution potentielle. Nous définissons une procédure et nous l’appliquons à l’étude de quelques situations simples en tomographie de luminescence. La dernière partie de ce travail présente la mise au point d’une procédure d’inversion. Après avoir présenté le cadre théorique dans lequel cette procédure s’inscrit nous la validons à partir de données numériques avant de l’appliquer avec succès à des mesures expérimentales. / Molecular imaging is a major modality in the field of preclinical research. Among the existing methods, techniques based on optical detection of visible or near infrared radiation are the most recent and are mainly represented by luminescence optical tomography techniques. These methods allow for 3D characterization of a biological medium by reconstructing maps of concentration or localisation of luminescent beacons sensitive to biological and chemical processes at the molecular or cellular scale. Luminescence optical tomography is based on a model of light propagation in tissues, a protocol for acquiring surface signal and a numerical inversion procedure used to reconstruct the parameters of interest. This thesis is structured around these three axes and provides an answer to each problem. The main objective of this study is to introduce and present the tools to evaluate the theoretical performances of optical tomography methods. One of its major outcomes is the realisation of experimental tomographic reconstructions from images acquired by an optical imager designed for 2D planar imaging and developed by the company Quidd. In a first step we develop the theory of transport in scattering medium to establish the concept on which our work will rely. We present two different propagation models as well as resolution methods and theoretical difficulties associated with them. In a second part we introduce the statistical tools used to characterise tomographic systems. We define and apply a procedure to simple situations in luminescence optical tomography. The last part of this work presents the development of an inversion procedure. After introducing the theoretical frameworkwe validate the procedure fromnumerical data before successfully applying it to experimentalmeasurements.

La génération actuelle des sondeurs infrarouges avancés constitue l’une des sources les plus importantes d’observation dans les systèmes d’assimilation de données dans les modèles de la Prévision Numérique du Temps (PNT). Cependant la richesse d’informations fournies par ce type de capteur avec son grand nombre de canaux et sa couverture globale est loin d’être complètement exploitée. La présence de nuages dans le champ de vision de l’instrument, qui affecte la majorité des observations, est l’une des raisons pour lesquelles les centres de PNT rejettent une grande quantité des observations des sondeurs. Les centres de PNT ont cependant commencé à assimiler au-dessus des océans les radiances affectées par les nuages en utilisant des canaux dont les effets radiatifs nuageux sont modélisés par un modèle de nuage simple. Certains de ces algorithmes de détection sont évalués dans ce manuscrit, et leurs limitations sont explicitées. Afin d’accroître la quantité de données assimilées, il est nécessaire de mieux représenter les nuages et leurs effets radiatifs. Depuis quelques années, des études ont été menées pour mieux représenter leurs effets dans les modèles de transfert radiatif ([Faijan et al., 2012] ; [Martinet et al., 2013]) ; et utiliser dans l’assimilation de nouveaux canaux infrarouges informatifs sur les hydrométéores nuageux. ([Martinet et al., 2014]). Cette thèse se concentre principalement sur ces méthodes de détection de scènes homogènes en consacrant sa majeur partie à l’établissement, l’évaluation et l’amélioration d’algorithme de détection de scènes homogènes en se basant sur la colocalistion d’observation avec d’autres sondeurs. Ces études sont rendus possibles par la prise en compte des champs d’hydrométéores fournis par les schémas convectifs du modèle ARPEGE en entrée du modèle de transfert radiatif nuageux RTTOV-CLD. Une partie validation des simulations est opérée dans cette thèse, en comparant l’apport les forces et faiblesses du schéma convectif en opérationnel ainsi que PCMT. Par la suite, différents tests, ou critères, de détection sont proposés, et en réalisant des expériences d’assimilation et en évaluant l’impact de ces ces critères de sélection proposés sur la qualité des prévisions à longues échéances, un des tests parmi ceux proposés se démarque des autres en conservant une quantité importante d’observation ciel clair et démontre des impacts neutres à légèrement positifs sur les prévisions. Les nouvelles méthodes de sélection de scènes homogènes proposées dans cette études permettent d’envisager une amélioration significative du contrôle de qualité des observation IASI en ciel clair. Cela ouvre ainsi donc la voie à une utilisation, plus maîtrisée, des scènes nuageuses. Nous expliquons dans ce manuscrit pourquoi il serait imprudent de précéder à des assimilation de radiances infrarouge contaminées par la présence de nuages. Pour contourner cette difficulté, une technique d’assimilation en deux étapes déjà utilisé pour l’assimilation des réflectivité radar ([Wattrelot et al., 2014]) dans AROME est évaluée. Cette méthode basée sur l’inversion bayésienne a récemment été adaptée pour les observations microondes satellitaire ([Duruisseau et al., 2018]). Dans cette étude, nous explorons cette technique pour les observations IASI. Plusieurs tests de sensibilité sont effectués sur les différents paramètres de l’algorithme, avec pour objectif de préparer de futurs travaux sur l’assimilation all-sky infrarouges, explicités dans les perspectives de ce manuscrit. / The current generation of advanced infrared sounders is one of the most important sources of observations in data assimilation systems in numerical weather prediction (NWP) models. However, the total amount of information provided by this type of sensor, with its large number of channels and its global coverage, is far from being fully exploited. The presence of clouds in the instrument’s field of view, which affects the majority of observations, is one of the reasons why NWP centers reject a large amount of observations from sounders. NWP centers, however, have begun to assimilate cloud-affected radiances over the oceans using channels whose cloudy radiative effects are modeled by a simple cloud model. Some of these detection algorithms are evaluated in this manuscript, and their limitations are clarified. In order to increase the amount of assimilated data, it is necessary to better represent clouds and their radiative effects in the models. For several years, studies have been conducted to better represent their effects in radiative transfer models ([Faijan et al., 2012] ; [Martinet et al., 2013]) ; and to use new informative infrared channels of cloudy hydrometeors in the assimilation. [Martinet et al., 2014]. This thesis focuses on several approaches for the assimilation of cloudy radiances. In the first part, the characterization of the cloud parameters currently used for the assimilation of cloudy radiances was evaluated in the global and regional scale models. In addition, as part of the "all-sky" assimilation, which considers both clear and cloudy radiances, the evaluation and improvement of homogeneous scene detection algorithms based on the colocation of observations with other imagers was studied. These studies are made possible by taking into account the hydrometeorological fields provided by the convective schemes of the ARPEGE model as the input of the RTTOV-CLD cloud radiative transfer model. Part of this thesis concerns the validation of simulations, by comparing the contribution of the new convective PCMT scheme to the one used in operational applications. Subsequently, different criteria for selecting homogeneous scenes are proposed. By conducting assimilation experiments and evaluating the impact of these proposed selection criteria on the quality of long-term forecasts, one of the proposed tests stands out from the others by keeping a significant amount of clear sky observations and demonstrating neutral to slightly positive impacts on the forecasts. These new methods for selecting homogeneous scenes proposed in this study allows the consideration of improving the quality control of IASI observations in clear sky. To address the issue of all-sky radiance data assimilation, the two-step assimilation technique, already used for radar reflectivity assimilation in AROME ([Wattrelot et al., 2014]), was evaluated for IASI radiances in the ARPEGE model in a case study. This method based on Bayesian inversion has recently been adapted for satellite microwave observations ([Duruisseau et al., 2018]). Several sensitivity tests were carried out on the different parameters of the algorithm, with the objective of preparing for future work on infrared all-sky assimilation, as explained in the perspectives of this manuscript.

La comprehension des processus conduisant à la formation des étoiles est l'un des enjeux majeurs de l'astrophysique contemporaine. Au sein des nuages conduisant à la formation d'étoiles, les conditions de température, pression, etc... sont telles qu'il est impossible de les reproduire par l'expérience. C'est pourquoi la simulation numérique reste le seul moyen d'étudier les phénomènes physiques intervenant dans le processus de formation des étoiles et ainsi de vérifier la théorie. Ma thèse est axée autour des méthodes numériques utilisées dans le contexte de la formation d'étoiles, phénomène multi-échelles et hautement non-linéaire, nécessitant l'utilisation d'outils bien adaptés. Dans cette thèse autour de l'étude des premières phases de l'effondrement de coeurs denses préstellaires, mon travail s'est divisé en 4 parties liées. Dans une première étude, j'ai utilisé un code lagrangien 1D à symétrie sphérique (Audit et al. 2002) pour comparer plusieurs modèles traitant plus ou moins précisément le transfert radiatif et l'interaction matière-rayonnement. Cette comparaison est basée sur des calculs simples d'effondrement gravitationnel conduisant à la formation du premier coeur de Larson. J'ai aussi tiré bénéfice de ce premier travail pour étudier les propriétés du choc d'accrétion sur le premier coeur de Larson. Nous avons développé un modèle semi-analytique permettant de reproduire les propriétés de saut au choc en partant d'hypothèses bien connues. Ayant validé les méthodes utilisées précédemment, nous avons retenu l'approche de diffusion à flux limité que j'ai ensuite intégrée avec les équations de l'hydrodynamique radiative dans le code AMR RAMSES (Teyssier 2002). Après validation des schémas implémentés, nous avons utilisé RAMSES pour réaliser des effondrements multidimensionnels avec champ magnétique et transfert radiatif. Nous avons ainsi réalisé les premières simulations combinant les effets du champ magnétique et du transfert radiatif aux petites échelles avec une grande précision. Nos résultats montrent que le transfert radiatif à un impact significatif sur la fragmentation au cours de l'effondrement des coeurs denses préstellaires. Enfin, j'ai réalisé une comparaison du code RAMSES (approche eulérienne) et du code SPH DRAGON (Goodwin 2004, approche lagrangienne). Nous avons étudié l'impact de la résolution numérique sur la conservation du moment angulaire et la fragmentation. Nous avons montré qu'en utilisant des critères de résolution forts et bien supérieurs aux critères usuels de la littérature, les deux outils convergent et semblent donc bien adaptés à la formation d'étoiles.

Dans un premier temps, nous nous intéressons à la modélisation de la couleur de la peau. Cette dernière est une structure multicouche, chacune contenant des centres diffuseurs de différentes tailles. L'interaction lumière/matière dans un tel milieu est modélisée par l'équation de transfert radiatif, résolue par la méthode de la fonction auxiliaire. Les propriétés optiques des centres diffuseurs sont soient déterminées par la théorie de Mie, soient trouvées dans la littérature. Un modèle optique est développé et validé par des spectres de réflexion diffuse mesurés sur des peaux réelles. L'influence des paramètres physiologiques est ensuite étudiée. Enfin, le problème est inversé et ces paramètres physiologiques sont déterminés à partir d'un spectre mesuré. Dans un second temps, nous nous intéressons aux carnations dans les peintures de chevalet occidentales. Une étude bibliographique et expérimentale sur de vraies œuvres révèle des similitudes entre les carnations et la peau réelle. Il en ressort aussi que les pigments ont peu changé au cours des siècles, contrairement aux techniques picturales (caractérisées principalement par le liant). L'influence des liants sur l'aspect visuel des peintures, c'est-à-dire la brillance (diffusion de surface) et la couleur (diffusion de volume), est étudiée expérimentalement sur des échantillons, réalisés avec 5 liants et 4 pigments. On discrimine trois types de liants: les liants aqueux, la tempera à l'œuf et l'huile de carthame. La brillance vient principalement du taux d'évaporation du solvant des liants et la couleur provient principalement des indices de réfraction des liants. Ces conclusions peuvent être légèrement modifiées selon les pigments.

Ce travail est consacré à l'étude de la diffusion de la lumière dans les tissus biologiques, en vue du développement d'une méthode de biopsie optique pour la détection précoce du cancer des organes creux. D'un point de vue théorique nous avons développé un modèle permettant de rendre compte de la diffusion de la lumière dans un système multistratifié, hétérogène, diffusant, dont la géométrie est proche de celle des tissus biologiques. Ce modèle est basé sur la résolution de l'Equation de Transfert Radiatif Vectorielle (ETRV) par la méthode des Ordonnées Discrètes. Il permet de calculer la répartition angulaire et les variations spectrales de la lumière diffusée suivant différents états de polarisation. Nos simulations du processus de cancérisation à l'intérieur d'un tissu épithélial ont montré que l'on peut espérer distinguer un tissu sain d'un tissu cancéreux à partir de l'analyse de l'intensité et de l'état de polarisation de la lumière diffusée en fonction de sa direction de diffusion. L'élément morphologique discriminant étant le doublement de la taille du noyau des cellules, principal diffuseur de la lumière, dans le tissu dysplasique. D'un point de vue expérimental nous avons réalisé des mesures de diffusion de lumière sur des échantillons école (lait et intralipides) à l'aide d'un gonio-spectro-photomètre construit au laboratoire. Elles ont permis la validation de notre modèle. Dans la continuité de ces expériences nous avons étudié des tissus reconstitués dont les cellules présentaient des tailles de noyaux différentes et des peaux de souris dont l'une présentait vraisemblablement un début de tumeur. Les résultats de ces mesures sont en accord quantitatif (pour les tissus de culture) et qualitatif (pour les essais préliminaires sur la peau de souris) avec les prédictions du modèle. Ils sont encourageants pour poursuivre notre étude de la diffusion élastique de la lumière comme outil de diagnostic.

Ce travail porte sur l'étude du plasma source d'un réacteur de type "hélicon" pouvant être utilisé pour la gravure des semiconducteurs. Nous utilisons un mélange Ar-He-F2 à des pressions variant de 4 à 20 mTorr avec des conditions de mélanges variables. L'interférométrie microondes (68 GHz) a été utilisée pour mesurer la densité électronique dont les valeurs vont de 0.4 à 2 10^12 cm-3 selon les conditions de pression et de puissance RF. En se basant sur la mesure de l'intensité des raies spectrales de l'He, nous avons établi un modèle collisionnel radiatif de la décharge. Il nous a permis de caractériser la fonction de distribution en énergie des électrons (fdee). Un modèle à deux températures a été choisi pour décrire cette fdee. La variation de ces températures en fonction des conditions de décharge est présentée. La température du gaz, trouvée de 850 K, a été mesurée par absorption laser qui a permis de mesurer la densité des niveaux 4s d'Ar et de valider la technique d'auto-absorption. La cinétique des raies fluor a été étudiée et a permis d'obtenir pour la première fois les coefficients d'excitation électroniques du fluor. Ces données ainsi que la mise en évidence des voies de population des états excités d'Ar et d'He montrent les limites de l'applicabilité de l'actinométrie. Une approche théorique est présentée pour montrer l'importance des différents processus responsables de l'entretien de la décharge et du dépôt d'énergie.

La raie hyperfine à 21-cm de l'hydrogène atomique émise dans l'univers à grand redshift (6 < z < 20) permettra de sonder une nouvellère de la cosmologie, l'époque de la réionisation (EoR). Elle fournira des informations plus détaillées, moins ambigües et plus complètes en 3D que les autres observations possibles de l'EoR (comme les raies d'absorption des quasars, ou le CMB). Son intensité dépend de nombreux processus physiques. Les radiotélescopes de nouvelle génération, SKA et ses précurseurs, entreront en service dans les prochaines années pour observer ce signal. La simulation numérique du 21-cm est importante pour optimiser le design des instruments et interpréter les observations futures. <br />Dans ce travail, nous avons développé un module de transfert radiatif pour le continuum ionisant dans le code LICORICE dans le but de simuler l'époque de la réionisation, pour laquelle le transfert radiatif est un outil essentiel. Nous utilisons un algorithme de ray-tracing de type Monte Carlo sur une grille adaptative. Plusieurs tests sont réalisés pour valider le code dans deux cas, celui d'un champ de densité statique et dans le cas de l'hydrodynamique radiative. <br />Nous produisons ensuite le signal à 21-cm émis pendant l'EoR, qui nous fournit une sonde directe de la réionisation et contient de nombreuses informations sur les sources d'ionisation et de chauffage. Habituellement, on suppose que le flux Ly-α, qui influence la puissance du signal à 21-cm, est fort et homogéne pendant toute la réionisation, mais cette hypothèse n'est pas valable durant la première phase de l'EoR. Nous évaluons exactement le flux Ly-α local en utilisant le module de transfert radiatif de la raie Lyman α et montrons que ce flux modifie le spectre de puissance du signal à 21-cm, particuliérement pendant la première phase de la réionisation. Nous trouvons également que le signal est en absorption forte ce qui ne peut pas être simulé avec les hypothèses habituelles.<br /> Même une faible quantité de rayons X peut affecter l'état physique du gaz neutre dans le milieu inter galactique, puisque ceux-ci ont un libre parcours moyen très long. Le chauffage par les rayons X est un facteur important pour évaluer la puissance du signal à 21-cm dans les régions neutres, aussi avons nous inclus le ray-tracing correspondant dans LICORICE. Nous constatons que les rayons X ont besoin de temps pour élever la température du gaz au-dessus de la température du CMB. Ainsi, avec un niveau de rayons X raisonnable, nous continuons à observer le signal en absorption pendant la premi`ere phase de la réionisation.<br />De plus, nous étudions la formation des galaxies dans un contexte cosmologique. Le but est d'explorer la physique des galaxies primordiales en considérant la rétroaction radiative sur l'hydrodynamique. Nos résultats préliminaires montrent que le calcul de l'ionisation hors équilibre modifie la température du gaz diffus, et son histoire d'accrétion.

La vapeur d'eau atmosphérique est une ressource en eau sous exploitée. Les matériaux à forte émissivité infrarouge (rayonnement thermique) sont utilisés en climatisation passive de jour. De nuit, ils se refroidissent et sont une technologie bon marché, robuste et efficace pour condenser cette ressource renouvelable avec des rendements mesurés supérieurs à 0,6 L m-2 par nuit. L'enjeu de ces travaux de recherches était d'avancer la compréhension de ces systèmes, de mettre au point les outils et de développer les techniques permettant la réalisation de condenseurs de grandes dimensions.<br /> Des matériaux (films polymères et peintures) originaux et innovants ont été développés pour ces applications : ils conjuguent une grande émissivité IR et un état de surface optimisé pour un coût inférieur à 1 € m-². Des simulations numériques CFD (Phoenics) ont été programmées dans le but d'optimiser la mise en œuvre des condenseurs de rosée sur des systèmes réels complexes et de grande superficie. De nouvelles méthodes d'estimation (mesure de volume ou estimation indirecte) de la ressource potentielle en eau de rosée sont discutées en fonction de l'influence des paramètres météorologiques pour les îles de Biševo (Croatie), Tahiti et Tikehau (P. Française) et les sites de Panandhro et Yellagiri Hills (Inde). <br />Enfin, deux systèmes réels novateurs ont été réalisés : une toiture pérenne de démonstration de <br />15,1 m² installée à Biševo (rendements moyens de 0,181 mm pour un maximum de 0,526 mm). Le second est la première usine à rosée, réalisée avec l'IIMA, dans l'état du Gujarat (Inde). La surface actuelle est de 850 m² et sera étendue à 15 000 m² pour produire jusqu'à 5 m3 d'eau pure par nuit.

Présentés dans un cadre de sécurité pyrotechnique, ces travaux ont pour principal objectif de simuler les effets thermiques de la combustion de blocs de propergol solide aluminisé, dispersés après l'éclatement d'un propulseur.<br />Un état de l'art des mécanismes mis en jeu aux niveau micro- et macroscopiques a permis de mieux appréhender les transferts de chaleur et de masse se produisant lors de ce type de combustion. des essais à grande échelle ont été réalisés afin d'obtenir de nouvelles données quantitatives, nécessaires au développement du modèle théorique. Une comparaison des résultats pour deux compositions de propergol, l'une hautement chargée en aluminium (20% en masse), l'autre plus faiblement (4%), a montré le rôle prépondérant que les gouttes d'aluminium puis les particules d'alumine jouent dans les transferts radiatifs vers l'environnement extérieur.<br />Cette tendance est également obtenue par un calcul numérique, basé sur le code de simulation d'incendie FDS (Fire Dynamics Simulator), qui permet de calculer la dispersion des gouttes d'aluminium par une approche lagrangienne. La combustion des gouttes A1/A12O3 est prédite par un modèle de type fraction de mélange, couplé à une loi d'évaporation, et leur rayonnement est calculé à partir de la théorie de Mie. Une étude de sensibilité a en outre permis de déduire un jeu de paramètres optimal pour lequel les résultats numériques sont en bon accord avec les mesures expérimentales. Dès lors, un scénario accidentel complet peut être simulé dans un bâtiment d'assemblage de gros propulseurs.