Pour concevoir des foyers aéronautiques plus fiables et moins polluants, les industriels ont de plus en plus recours à des simulations numériques s'appuyant sur de nombreux modèles physiques. Si l'on s'intéresse en particulier aux problématiques des charges thermiques pariétales et des émissions polluantes, la modélisation des phénomènes couplés de combustion, de formation des suies et de transfert radiatif est nécessaire. Ainsi, cette thèse a pour objectif de développer une méthodologie permettant de simuler ces phénomènes couplés de manière instationnaire, dans un foyer représentatif des systèmes industriels. Un modèle de formation des suies simple et robuste, à caractère empirique, a d'abord été mis au point sur une configuration de flamme 1D laminaire prémélangée. Ce modèle a été choisi car, étant compatible avec des mécanismes réactionnels globaux, il est bien adapté aux simulations instationnaires en géométrie complexe. Dans un deuxième temps, il a été appliqué à la simulation instationnaire de l'écoulement turbulent réactif diphasique dans un foyer doté d'un prototype d'injecteur industriel. Les niveaux de température obtenus ainsi que la topologie du champ de fraction volumique de suies ont été comparés aux résultats expérimentaux. Les écarts constatés ont été analysés et des corrections ont été proposées. Enfin, une stratégie de couplage entre l'approche LES (Large Eddy Simulation) et un outil de simulation des transferts radiatifs basé sur la méthode de Monte Carlo a été mise au point et éprouvée sur le foyer étudié. L'écoulement apparaît peu affecté par le rayonnement, mais en revanche, les transferts radiatifs ont un impact relativement important sur les flux reçus par les parois

L’objectif de ces travaux de thèse est d’améliorer les connaissances sur les réacteurs photochimiques afin de permettre une conception, un dimensionnement et un pilotage plus efficace, compatible avec des contraintes industrielles. Ces travaux se basent sur une méthode d’intensification des procédés et visent la mise en œuvre de réactions de photochimie préparative. La photochimie préparative est une voie de synthèse particulièrement attrayante dans le contexte de la chimie verte. Elle donne accès à des composés hautement fonctionnalisés (molécules à cycle tendu, hétéroatomes etc.) difficilement ou pas accessibles par voie thermique, de manière sélective, en très peu d’étapes réactionnelles et sans ajout des réactifs supplémentaires. Malgré ses potentialités, la transposition industrielle demeure limitée. Les principales contributions de cette thèse sont d’avoir analysé finement les bénéfices de la petite échelle pour la conduite de réactions photochimiques et d’avoir proposé une méthodologie pour transposer des réactions photochimiques dans des équipements industriels. Dans un premier temps, un outil numérique a été développé afin de modéliser le couplage entre les différents phénomènes rencontrés en microphotoréacteur (transfert radiatif, transport des espèces et cinétique photochimique). Une formulation du système d’équations basée sur des nombres sans dimension a été utilisée de façon à accéder à une vision indépendante de l’échelle de l’expérimentation. A cet effet, les nombres sans dimension classiquement rencontrés en génie de la réaction ont dû être adaptés aux spécificités des réactions photochimiques. Quatre nombres sans dimension contrôlant les performances en sortie de microphotoréacteur ont été mis en exergue : les nombres de Damköhler I et II, l’absorbance dans le milieu et le facteur de compétition d’absorption des photons incidents. A partir de ce modèle, une cartographie décrivant les différents régimes de fonctionnement rencontrés en microphotoréacteur a été établie. Une attention particulière a été ensuite portée sur les cas de fonctionnement où la conversion en sortie du réacteur est significativement réduite par la compétition des espèces pour absorber les photons. Cette compétition est liée au nombre de Damköhler II qui permet d’évaluer l’efficacité du mélange diffusif transverse. Ces simulations numériques ont abouti à la construction d’un abaque pour choisir les conditions opératoires à imposer afin de ne pas être limité par le mélange transverse. Dans un second temps, différents outils et dispositifs expérimentaux ont été développés afin de caractériser les réacteurs (notamment par actinométrie) et de suivre en ligne par spectrophotométrie différents systèmes réactionnels. Des microphotoréacteurs de type « capillary tower » et en spirale ont été conçus. Ces dispositifs expérimentaux ont permis d’opérer dans une très large gamme de conditions opératoires (flux de photons, temps de séjour) afin de valider les observations numériques. Pour cela, différents systèmes photochimiques ont été mis en œuvre : une photocycloaddition [2+2] intramoléculaire et deux systèmes photochromiques. Les résultats expérimentaux obtenus, avec ou sans limitation par le transfert diffusif transverse, ont clairement confirmé la pertinence des observations numériques. En outre, la faisabilité d’utiliser un microréacteur comme outils d’acquisition de données cinétiques de réactions photochimiques a été démontrée. Finalement, sur la base des observations expérimentales et numériques, une méthodologie générale est présentée sous forme de logigramme pour déterminer les paramètres de dimensionnement en microphotoréacteur (temps de séjour et densité de flux de photons reçus à la paroi). Des critères ont été proposés pour caractériser les microphotoréacteurs : la productivité, et de manière plus originale, le rendement énergétique global (incluant l’efficacité photonique). / This work aims at improving the knowledge on photochemical reactor engineering in order to propose a methodology to implement photochemical reactions in new continuous intensified technologies. Synthetic organic photochemistry is an extremely powerful method for the conversion of simple substrates into complex products, opening new perspectives. As photochemical substrate activation often occurs without additional reagents, the formation of by-products is also minimized, making photochemistry even more attractive in the modern context of Green Chemistry. The main contributions of this thesis are to finely analyze the benefits of the microreactor technology for performing photochemical reactions, and to propose methodology to transpose photochemical reactions from lab scale to industrial scale. Firstly, a numerical modeling has been proposed to describe the coupling between the different physical phenomena occurring inside a microphotoreactor (radiative transfer, mass and momentum transfers, photochemical kinetic). A formulation of the equation system based on dimensionless numbers has been used to access a generic view, independent of the scale of the experiment. For that, the dimensionless numbers classically encountered in chemical reaction engineering have been adapted to account for photochemical reaction specificity. Four dimensionless numbers controlling the performances at the microreactor’s outlet have been outlined: the Damköhler I and II numbers, the absorbance inside the medium and the competitive absorption factor. From these numbers, a map describing the different zones in which a microphotoreactor can operate has been established. A special attention has been then paid on the cases where the conversion at the microreactor’s outlet is significantly reduced due to the occurrence of a photon competitive absorption between several species. The influence of this competition phenomenon is directly linked to the Damköhler II number which assess for the transverse mixing efficiency. In a second time, various experimental tools and set-up have been developed to characterize photochemical reactors (measurement of the photon flux density received by actinometry) and to implement online analysis by spectrophotometry. A “capillary tower” and a “spiral” microphotoreactors have been developed. Both these microphotoreactors have enabled to operate in a wide range of operating conditions (photon flux, residence time) so as to validate numerical simulations. For that, three photochemical systems have been implemented: an intramolecular [2+2] photocycloaddition and two photochromic systems. The experimental results obtained have confirmed the relevancy of the numerical observations, whether some mass transfer limitations occur or not. Moreover, the feasibility to use a microphotoreactor as a tool for acquiring kinetic data on photochemical reactions has been demonstrated. Finally, based on the numerical and experimental observations, a detailed flow chart has been built to rapidly determine the key parameters for scaling a microphotoreactor (residence time and photon flux density). Some criteria have been proposed to characterize the microphotoreactor: the productivity and, more originally, the global energetic yield (including the photonic efficiency).

Pour concevoir des foyers aéronautiques plus fiables et moins polluants, les industriels ont de plus en plus recours à des simulations numériques s’appuyant sur de nombreux modèles physiques. Si l’on s’intéresse en particulier aux problématiques des charges thermiques pariétales et des émissions polluantes, la modélisation des phénomènes couplés de combustion, de formation des suies et de transfert radiatif est nécessaire. Ainsi, cette thèse a pour objectif de développer une méthodologie permettant de simuler ces phénomènes couplés de manière instationnaire, dans un foyer représentatif des systèmes industriels. Un modèle de formation des suies simple et robuste, à caractère empirique, a d’abord été mis au point sur une configuration de flamme 1D laminaire prémélangée. Ce modèle a été choisi car, étant compatible avec des mécanismes réactionnels globaux, il est bien adapté aux simulations instationnaires en géométrie complexe. Dans un deuxième temps, il a été appliqué à la simulation instationnaire de l’écoulement turbulent réactif diphasique dans un foyer doté d’un prototype d’injecteur industriel. Les niveaux de température obtenus ainsi que la topologie du champ de fraction volumique de suies ont été comparés aux résultats expérimentaux. Les écarts constatés ont été analysés et des corrections ont été proposées. Enfin, une stratégie de couplage entre l’approche LES (Large Eddy Simulation) et un outil de simulation des transferts radiatifs basé sur la méthode de Monte Carlo a été mise au point et éprouvée sur le foyer étudié. L’écoulement apparaît peu affecté par le rayonnement, mais en revanche, les transferts radiatifs ont un impact relativement important sur les flux reçus par les parois / Numerical simulations, involving numerous physical models, are more and more employed to design more reliable and less pollutant industrial combustors. In particular, focusing on wall thermal load and pollutant emission issues, coupled phenomena such as combustion, soot formation and radiative transfer have to be modelled. Thus, the aim of this PhD thesis is to develop a methodology to simulate these coupled phenomena in an unsteady way, in an industrial-like combustion chamber. A simple and robust empirical soot formation model has been developed and applied in a first step to a 1D laminar premixed flame configuration. This model was chosen because it is well adapted to unsteady simulations in complex geometries, as it is compatible with global reaction mechanisms. In a second step it was applied to the unsteady simulation of the two-phase turbulent reactive flow in a combustor equipped with an industrial injector prototype. Predicted temperature levels and topology of the soot volume fraction field were compared to experimental results. Some discrepancies were observed: they were analysed and corrections were proposed. Finally, a coupling strategy between the LES (Large Eddy Simulation) approach and a radiative transfer simulation tool based on the Monte Carlo method was developed and tested on the same combustor. It appears that radiative transfer does not greatly modify the flow, but has a relatively important effect on wall heat fluxes.

Ma thèse a été consacrée à la modélisation du rayonnement thermique, abordée d'un point de vue électromagnétique. Plus précisément, elle s'applique à l'étude de l'excitation thermique d'ondes de surface (plasmon-polaritons ou phonon-polaritons) susceptibles d'exister pour certains matériaux métalliques ou diélectriques. L'originalité de ce travail a consisté à évaluer la contribution de ces ondes de surface à l'émission thermique de matériaux micro- ou nanostructurés (réseaux) ; aux transferts d'énergie par rayonnement dans des systèmes de taille micro- ou nanométrique ; aux forces de Casimir dans la limite électrostatique. Au cours de cette étude, nous avons détaillé le formalisme qui nous a permis d'aborder le rayonnement thermique d'un point de vue électromagnétique. Nous avons pu donner une expression exacte et explicite de la densité d'énergie électromagnétique (due au rayonnement thermique) au-dessus d'une interface plane séparant un milieu quelconque du vide, en explicitant notamment la contribution des ondes évanescentes. Nous avons tenté de définir une quantité qui serait l'analogue, pour les ondes évanescentes, de l'émissivité monochromatique directionnelle. L'expression de la densité d'énergie obtenue montre qu'en champ proche, il est théoriquement possible de déterminer localement la constante diélectrique d'un substrat en mesurant le champ d'émission thermique qu'il "rayonne" en champ proche. Par ailleurs, nous montrons qu'un microscope optique en champ proche, détectant le champ d'émission thermique d'un substrat est l'analogue optique d'un STM (Scanning Tunneling Microscope) pour les électrons. Nous pourrions ainsi mesurer la densité d'états locale du champ électromagnétique. Utilisant le phénomène d'onde de surface dans l'infrarouge pour le SiC, nous avons été capable de dimensionner une source thermique présentant une certaine cohérence spatiale : dans une direction fixée, elle émet préférentiellement à une longueur d'onde et pour une longueur d'onde fixée, elle émet à l'intérieur d'un lobe très étroit angulairement. Les mesures expérimentales confirment avec un excellent accord cette prédiction théorique. Ainsi, la source que nous avons réalisée avec ce réseau de SiC est le premier exemple expérimental de source thermique patiellement cohérente spatialement. Par ailleurs, cette source a un spectre d'émission qui dépend de l'angle d'émission, c'est l'"effet Wolf". Les développements actuels au laboratoire concernent la mise en œuvre expérimentale d'une expérience de mesure d'émissivité infrarouge avec un spectromètre à transformée de Fourier. Nous abordons ensuite le problème du transfert radiatif entre deux milieux semi-infinis séparés par une faible épaisseur de vide, de 10 nm à plusieurs dizaines de microns et montrons que l'approche radiométrique du rayonnement thermique n'est pas valable lorsque les échelles caractéristiques deviennent du même ordre de grandeur que la longueur d'onde du rayonnement. Nous avons montré alors que les ondes de surface (dont nous ne pouvons rendre compte qu'avec une approche électromagnétique), et plus particulièrement les phonon-polaritons de surface dans le cas du SiC, jouent un rôle fondamental dans le transfert d'énergie électromagnétique. Premièrement, celui-ci est (\it quasi) monochromatique et présente un pic très prononcé à la fréquence de résonance du polariton de surface qui donne la contribution majeure au transfert. Une expression asymptotique en champ proche du coefficient de transfert radiatif en fonction de la fréquence est également donnée. Deuxièment, le transfert est amplifié de plusieurs ordres de grandeur à courte distance, l'amplification variant comme l'inverse du carré de la distance séparant les deux milieux. Pour le cas de deux milieux de SiC séparés d'une distance de 10 nm, la contribution du rayonnement à l'échange d'énergie est du même ordre de grandeur que celle de la conduction (pour un gaz dans les conditions normales). Il est donc essentiel de ne pas négliger cette contribution du rayonnement dans ce type de systèmes, comme cela est fait couramment. Une application envisagée de ce transfert radiatif concerne l'effet thermo-photovoltaïque amplifié. Dans une configuration où nous avons considéré une particule sub-longueur d'onde (approximation dipolaire) située au-dessus d'une interface plane, nous avons mis en évidence le même type d'effets, lorsque la particule et la surface supportent des ondes de surface. La particule absorbe ainsi de manière importante autour de deux fréquences bien précises, correspondant aux fréquences de résonance des ondes de surface. La puissance absorbée par la particule augmente de plusieurs ordres de grandeur en champ proche, variant comme l'inverse de la puissance troisième de la distance entre la particule et le substrat. De la même manière, nous avons réalisé une cartographie de la distribution spatiale de la puissance absorbée par unité de volume à l'intérieur du substrat lorsque celui-ci est éclairé par le champ d'émission thermique de la particule. Nous avons montré que dans le cas d'une particule de SiC au-dessus d'un échantillon de SiC, la densité de puissance pouvait atteindre 100 MW/m3. De plus, la puissance déposée est confinée sur une zone dont la taille typique est de l'ordre de la distance entre la particule et le substrat. Une application de ce travail se situe dans le stockage haute-densité sur des matériaux magnétiques ou à transition de phase. Dans le même thème, nous nous sommes intéressés au cas de nanoparticules métalliques enfouies dans une matrice diélectrique. Le but est d'expliquer certains résultats expérimentaux obtenus lorsque ces nanoparticules sont illuminées par des impulsions laser femtosecondes. Le problème est alors d'étudier la dynamique de relaxation électronique en régime basse fluence (variation de la température électronique faible) et haute fluence (variation importante) et notamment l'influence de la taille de la particule. Lors de cette 'étude, nous avons montré que l'introduction phénoménologique d'un mécanisme supplémentaire d'échange d'énergie entre les électrons et la matrice diélectrique permet de rendre compte de comportements observés expérimentalement. Le premier concerne la dépendance du temps de relaxation électronique en fonction de la taille de la particule lorsque nous sommes dans un régime où la température électronique n'est pas trop élevée (< 400 K). Le temps de relaxation diminue alors avec le rayon de la particule, montrant que ce mécanisme supplémentaire est essentiellement surfacique. Nous avons alors appliqué notre modèle dans un régime où les températures électroniques sont plus élevées (> 400 K), correspondant à des expériences où les nanoparticules sont illuminées par des impulsions de forte fluence. Dans ce régime, la prise en compte du mécanisme de surface conduit à deux conclusions : la première est que l'extrapolation des valeurs calculées pour de fortes fluences au régime des faibles fluences montre que le temps de relaxation est relativement indépendant de la taille de la particule ; la seconde est que la variation du temps de relaxation est fortement non linéaire en fonction de la fluence incidente. Enfin, nous avons étudié l'"effet Casimir" en choisissant une méthode de calcul qui nous permet une interprétation physique aisée des résultats. Nous avons alors été capable de tracer le "spectre" de la force de Casimir entre deux milieux semi-infinis de SiC à température nulle. Nous avons montré qu'en champ proche, la contribution majeure (pics dans le spectre) était donnée par les polaritons de surface, notamment celui qui existe dans l'ultraviolet pour le SiC. De plus, la "relation de dispersion" de la force de Casimir montre que cette contribution des ondes de surface couplées peut être elle-même divisée en deux : une contribution due à ce que nous avons appelé des modes "liants", qui sont attractifs et une contribution due à des modes "anti-liants" qui sont répulsifs. Nous avons identifié, pour les ondes évanescentes, dans la formule analytique de la force les termes correspondant à chacun de ces modes.

L'objectif de ce travail est d’évaluer le positionnement et le potentiel des récepteurs à lit fluidisé à changement de section par rapport aux autres méthodes de chauffage de gaz à haute température par voie solaire. A cette fin, une connaissance approfondie des phénomènes thermiques et hydrodynamiques du récepteur est nécessaire. Pour acquérir cette connaissance, nous avons modélisé les transferts thermiques dans le récepteur en portant une attention particulière sur les transferts radiatifs en prenant en compte les diffusions multiples de la lumière dans le milieu particulaire, les effets de parois sur les transferts radiatifs et la directionnalité du rayonnement solaire concentré. La détermination précise de la distribution de particules dans le ciel du lit fluidisé s'est avérée un paramètre critique pour le calcul des transferts thermiques. Ces modèles, plus tard affinés par une confrontation avec des références expérimentales, nous ont permis d'explorer l'effet de la géométrie sur les transferts thermiques dans le récepteur. Il ont permis entre autres de mettre en évidence l'intérêt d'utiliser une colonne de fluidisation à changement de section et l'importance de l'optimisation du couple concentrateur solaire / récepteur afin d'éviter d'éventuelles surchauffes au niveau des parois du récepteur. De même, il semble que l'homogénéisation de la température dans le lit fluidisé contenu dans le récepteur soit favorable à son rendement. / The aim of this work is to evaluate the position and the potential of solar fluidized bed receivers compared to other methods for the solar heating of gases at high temperature. To this end, a thorough knowledge of the heat transfer and hydrodynamic of the receiver is necessary. To acquire this knowledge, we modeled the heat transfer in the receiver with a focus on the radiative transfer by taking into account the multiple scattering of light in the particle medium, the effect of walls on radiative heat transfer and the directionality of the concentrated solar radiation. The accurate determination of the distribution of particles within the fluidized bed has been a critical parameter for the calculation of heat transfer. With these models, later refined by a confrontation with experimental references, we have studied the effect of geometry on heat transfer in the receiver. This study highlighted the necessity to use a switching section fluidization column and the importance to optimize the pair : solar concentrator / receiver to avoid any overheating at the walls of the receiver. Moreover, it appears that the homogenization of the temperature in the fluidized bed of the receiver increase its performance.

Les thermoplastiques et les composites thermoplastiques sont généralement mis en œuvre par thermoformage et sont alors le plus souvent préchauffés en utilisant un chauffage IR. L’avantage du chauffage radiatif est qu'il permet de chauffer les polymères à cœur grâce au caractère semi-transparent des polymères. Néanmoins, dans le cas des polymères semi-cristallins, le chauffage radiatif est affecté par la structure cristalline et cette thèse a donc eu pour objectif d’améliorer la compréhension de l'interaction entre la structure cristalline et les propriétés optiques dans le but de proposer un modèle prédictif de chauffage de thermoplastiques semi-cristallins.Cette étude répond à une problématique industrielle relative au contrôle de la température des thermoplastiques semi-cristallins dans les procédés recourant au chauffage radiatif. L’optimisation de ces procédés requiert un code de calcul suffisamment robuste pour permettre une bonne prédiction du champ de température tout en conservant des temps de calcul acceptables. Une approche combinée expérimentale et numérique a ainsi été proposée dans le but de modéliser la capacité d’absorption du rayonnement thermique de milieux polymères semi-cristallins et le transfert de chaleur par rayonnement avec changement des phases de cristaux/amorphe. Ces travaux se concentrent sur le PEHD, qui présente un intérêt particulier pour l’entreprise Procter&Gamble.Dans cette thèse, après avoir établi une revue bibliographique mettant en avant les couplages existants entre les phénomènes de diffusion optique, la microstructure des polymères semi-cristallins et la température, une caractérisation et une analyse poussées des propriétés radiatives de deux polyéthylènes sont proposées. Les analyses morphologiques et optiques ont été réalisées à température ambiante et dans des conditions de chauffage afin d’identifier les formations cristallines à l’origine de la diffusion optique dans des polymères semi-cristallins et l’évolution de ce couplage au cours du chauffage. A travers ce travail de recherche, un coefficient d’extinction spectral thermo-dépendant a été proposé afin de décrire le caractère optiquement hétérogène du milieu semi-cristallin par un milieu homogène équivalent. Sur la base de la caractérisation de la capacité d'absorption du rayonnement thermique, un modèle thermique conducto-radiatif thermo-dépendant a été développé. Afin d’évaluer la précision de la modélisation, une méthodologie expérimentale spécifique a été proposée pour mesurer la température de surface par thermographie IR dans le cas du PEHD semi-transparent. L’étape finale a consisté à confronter les résultats issus des simulations numériques basées sur cette modélisation à plusieurs campagnes de mesures expérimentales. Les résultats de ces travaux démontrent la forte influence de la structure morphologique des polymères semi-cristallins sur les transferts de chaleur radiatifs. / Thermoplastics and thermoplastic composites are promising candidates for manufacturing highly cost- effective and environmental-friendly components in terms of rapid forming and recyclability. Thermoforming is extensively used for the processing of thermoplastics where IR heating is widely applied. The major advantage of radiative heating is that the significant portion of radiation penetrates into the semi-transparent polymer media.This thesis focuses on understanding of IR heating of semi-crystalline thermoplastics which aims to analyze the driven mechanisms for radiation transport in optically heterogeneous unfilled semi-crystalline polymer media. Considering the relatively narrow thermoforming window of semi-crystalline thermoplastics, accurate temperature control and close monitoring of temperature field is crucially important for successful forming process. It is thus required to build a numerical model robust enough to allow a good prediction of the temperature field while maintaining acceptable calculation times. In this research work, a combined experimental-numerical approach has been proposed which enables both to characterize the radiation absorption capacity of semi-crystalline polymer media and, to model the radiation heat transfer considering the crystalline/amorphous phases change under heating. This research focuses on a particular polymer - highly crystalline HDPE- which is supported by Procter & Gamble.In this thesis, the literature was reviewed at first for highlighting the existing coupled relation between the optical properties and the crystalline structure of semi-crystalline polymers. The role of crystalline morphology on the optical properties and optical scattering of two type of polyethylene, namely HDPE and LLDPE, were addressed. More specifically, the morphological and optical analyses were performed at room temperature and under heating to determine: which crystalline formations are responsible for optical scattering in semi-crystalline polymer media and, how does their coupled relationship evolve under heating conditions? Hence, one of the key contributions of this research is on establishing a temperature-dependent spectral extinction coefficient of HDPE allowing to describe temperature- dependent radiation absorption capacity of its semi-crystalline medium and, to model radiative transfer considering an equivalent homogeneous medium. Based on the characterization of radiation absorption capacity of semi-crystalline media, a temperature-dependent conduction-radiation model was developed. In order to assess the modeling accuracy, an experimental methodology was proposed for non-invasive surface temperature measurements via IR thermography on semi-transparent polymer media. The final step was to compare the results of numerical simulations with the several IR heating experiments to prove the strong influence of the crystalline morphology on heat transfer.

L’objet de ce travail est de proposer un modèle adapté pour la simulation du remplissage de moules qui réponde au meilleur compromis entre temps de calcul et précision des résultats. La difficulté est double. Il faut prendre en compte le phénomène de remplissage qui est un problème complexe à frontières libres et les spécificités liées au Verre : viscosité fortement thermodépendante et température de fusion élevée qui nécessite de prendre en compte le rayonnement. Le Chapitre I est consacrée à la partie écoulement du Verre liquide. La bibliothèque numérique Aquilon/Thétis, adaptée pour traiter ce type de problèmes et les couplages thermique air/verre/parois, a été utilisée (Méthode V.O.F pour le suivi de l’interface, méthodes de type Lagrangien augmenté/Projection vectorielle pour le couplage Vitesse-Pression). Pour l’aspect radiatif, différentes approches sont proposées : conductivité radiative équivalente (Chapitre II), méthode explicite directe pour la validation (Chapitre III) et méthode d’harmoniques sphériques ou méthode PN (Chapitre IV). Dans le Chapitre V, la méthode PN retenue est validée dans des cas simples et est appliquée ensuite à des cas avec couplage convectif en géométries complexes et obstacles semi-transparents (1D, 2D et 3D, 2D axi-symétrique et milieu non gris). Une version P1 modifiée est présentée. Les résultats sont assez proches de ceux donnés par la méthode P3 avec des temps de calcul modestes. L’intérêt de ce modèle est qu’il est facilement intégrable dans des codes numériques existants : une seule équation différentielle du second ordre stationnaire à résoudre en 3D / The aim of this study is to propose an adapted model for the simulation of mould filling that must be a compromise solution between computational time and results accuracy. The double difficulty is to take into account the filling phenomenon that is a complex problem due to the presence of free boundaries and to the Glass specificities: viscosity that is highly thermal dependant and high melting temperature that requires taking into account radiation effects. Chapter I is devoted to the melting Glass flow. The numerical libraries Aquilon/Thétis, adapted for solving such type of problems and the thermal coupling between Air/Glass/Walls, has been used. (V.O.F method for front tracking, Augmented Lagrangian/Vector Projection methods for solving Pressure/Velocity coupling). For radiative aspect, different approaches are proposed: equivalent radiative conductivity (Chapter II), direct explicit method for validation (Chapter III) and spherical harmonics method or PN method (Chapter IV). In the Chapter V, the selected PN method is validated through simple cases and is then applied in other cases with convective coupling in complex geometries including semi-transparent inclusions (1D, 2D and 3D, 2D axi-symmetric and non grey medium). A P1 modified version is presented. The results are close to those given by P3 method but with reduced computational time. The main interest of this model is that it can be easily implemented in existing numerical codes: a single stationary second order partial differential equation to solve in 3D