<p>Severe accidents in a Light Water Reactor (LWR) have been a subject of the research for the last three decades. The research in this area aims to further understanding of the inherent physical phenomena and reduce the uncertainties surrounding their quantification, with the ultimate goal of developing models that can be applied to safety analysis of nuclear reactors. The research is also focusing on evaluation of the proposed accident management schemes for mitigating the consequences of such accidents.</p><p>During a hypothetical severe accident, whatever the scenario, there is likelihood that the core material will be relocated and accumulated in the lower plenum in the form of a debris bed or a melt pool. Physical phenomena involved in a severe accident progression are complex. The interactions of core debris or melt with the reactor structures depend very much on the debris bed or melt pool thermal hydraulics. That is why predictions of heat transfer during melt pool formation in the reactor lower head are important for the safety assessment.</p><p>The main purpose of the present study is to advance a method for describing turbulent natural convection heat transfer of a melt pool, and to develop a computational platform for cost-effective, sufficiently-accurate numerical simulations and analyses of Core Melt-Structure-Water Interactions in the LWR lower head during a postulated severe core-melting accident.</p><p>Given the insights gained from Computational Fluid Dynamics (CFD) simulations, a physics-based model and computationally-efficient tools are developed for multi-dimensional simulations of transient thermal-hydraulic phenomena in the lower plenum of a Boiling Water Reactor (BWR) during the late phase of an in-vessel core melt progression. A model is developed for the core debris bed heat up and formation of a melt pool in the lower head of the reactor vessel, and implemented in a commercial CFD code. To describe the natural convection heat transfer inside the volumetrically decay-heated melt pool, we advanced the Effective Convectivity Conductivity Model (ECCM), which was previously developed and implemented in the MVITA code. In the present study, natural convection heat transfer is accounted for by only the Effective Convectivity Model (ECM). The heat transport and interactions are represented through an energy-conservation formulation. The ECM then enables simulations of heat transfer of a high Rayleigh melt pool in 3D large dimension geometry.</p><p>In order to describe the phase-change heat transfer associated with core debris, a temperature-based enthalpy formulation is employed in the ECM (the phase-change ECM or so called the PECM). The PECM is capable to represent possible convection heat transfer in a mushy zone. The simple approach of the PECM method allows implementing different models of the fluid velocity in a mushy zone for a non-eutectic mixture. The developed models are validated by a dual approach, i.e., against the existing experimental data and the CFD simulation results.</p><p>The ECM and PECM methods are applied to predict thermal loads to the vessel wall and Control Rod Guide Tubes (CRGTs) during core debris heat up and melting in the BWR lower plenum. Applying the ECM and PECM to simulations of reactor-scale melt pool heat transfer, the results of the ECM and PECM calculations show an apparent effectiveness of the developed methods that enables simulations of long term accident transients. It is also found that during severe accident progression, the cooling by water flowing inside the CRGTs plays a very important role in reducing the thermal load on the reactor vessel wall. The results of the CFD, ECM and PECM simulations suggest a potential of the CRGT cooling as an effective mitigative measure during a severe accident progression.</p>

Severe accidents in a Light Water Reactor (LWR) have been a subject of intense research for the last three decades. The research in this area aims to reach understanding of the inherent physical phenomena and reduce the uncertainties in their quantification, with the ultimate goal of developing models that can be applied to safety analysis of nuclear reactors, and to evaluation of the proposed accident management schemes for mitigating the consequences of severe accidents.  In a hypothetical severe accident there is likelihood that the core materials will be relocated to the lower plenum and form a decay-heated debris bed (debris cake) or a melt pool. Interactions of core debris or melt with the reactor structures depend to a large extent on the debris bed or melt pool thermal hydraulics. In case of inadequate cooling, the excessive heat would drive the structures' overheating and ablation, and hence govern the vessel failure mode and timing. In turn, threats to containment integrity associated with potential ex-vessel steam explosions and ex-vessel debris uncoolability depend on the composition, superheat, and amount of molten corium available for discharge upon the vessel failure. That is why predictions of transient melt pool heat transfer in the reactor lower head, subsequent vessel failure modes and melt characteristics upon the discharge are of paramount importance for plant safety assessment.  The main purpose of the present study is to develop a method for reliable prediction of melt pool thermal hydraulics, namely to establish a computational platform for cost-effective, sufficiently-accurate numerical simulations and analyses of core Melt-Structure-Water Interactions in the LWR lower head during a postulated severe core-melting accident. To achieve the goal, an approach to efficient use of Computational Fluid Dynamics (CFD) has been proposed to guide and support the development of models suitable for accident analysis.   The CFD method, on the one hand, is indispensable for scrutinizing flow physics, on the other hand, the validated CFD method can be used to generate necessary data for validation of the accident analysis models. Given the insights gained from the CFD study, physics-based models and computationally-efficient tools are developed for multi-dimensional simulations of transient thermal-hydraulic phenomena in the lower plenum of a LWR during the late phase of an in-vessel core melt progression. To describe natural convection heat transfer in an internally heated volume, and molten metal layer heated from below and cooled from the top (and side) walls, the Effective Convectivity Models (ECM) are developed and implemented in a commercial CFD code. The ECM uses directional heat transfer characteristic velocities to transport the heat to cooled boundaries. The heat transport and interactions are represented through an energy-conservation formulation. The ECM then enables 3D heat transfer simulations of a homogeneous (and stratified) melt pool formed in the LWR lower head. In order to describe phase-change heat transfer associated with core debris or binary mixture (e.g. in a molten metal layer), a temperature-based enthalpy formulation is employed in the Phase-change ECM (so called the PECM). The PECM is capable to represent natural convection heat transfer in a mushy zone. Simple formulation of the PECM method allows implementing different models of mushy zone heat transfer for non-eutectic mixtures. For a non-eutectic binary mixture, compositional convection associated with concentration gradients can be taken into account. The developed models are validated against both existing experimental data and the CFD-generated data. ECM and PECM simulations show a superior computational efficiency compared to the CFD simulation method. The ECM and PECM methods are applied to predict thermal loads imposed on the vessel wall and Control Rod Guide Tubes (CRGTs) during core debris heatup and melting in a Boiling Water Reactor (BWR) lower plenum. It is found that during the accident progression, the CRGT cooling plays a very important role in reducing the thermal loads on the reactor vessel wall. Results of the ECM and PECM simulations suggest a high potential of the CRGT cooling to be an effective measure for severe accident management in BWRs. / <p>QC 20100812</p>

Cette étude concerne la compréhension des mécanismes de transfert de chaleur et le développement d’un système de stockage pour la valorisation de la chaleur fatale industrielle. L’utilisation de Matériaux à Changement de Phase (MCP) permet d’atteindre une densité énergétique élevée et de restituer la chaleur à température constante. Cependant, leur faible conductivité thermique impose d’améliorer les transferts thermiques, notamment par l’utilisation d’échangeurs à surface augmentée. Le but est de comprendre le comportement de tels échangeurs en régime transitoire au contact de MCP. Une étude expérimentale à basse température, où quatre échangeurs de type tube-calandre ont été testés avec différentes orientations (horizontale/verticale) et injections (haut/bas), a mis en évidence des phénomènes de transfert thermique importants, comme la convection naturelle à la charge et la contraction volumique à la décharge. Ces observations ont été validées par un modèle CFD tridimensionnel. Une méthode de comparaison des performances basée sur un calcul d’énergie par le biais d’un maillage expérimental est proposée et permet de sélectionner un échangeur selon les critères de densités énergétiques, de temps caractéristique et de coût. Trois MCP, envisagés pour l’application, ont alors été testés à température réelle (100-200 °C) au contact d’un échangeur tube inox à ailettes transverses en aluminium pour évaluer leur cyclabilité et comparer leur comportement. Le mélange de sels, H105 (Tfusion = 122 °C), n’est pas retenu pour l’application à cause de sa faible densité énergétique (≈ 56 kWh/m3) et sa plage de fusion trop étalée. L’acide sébacique (Tfusion = 132 °C) a un comportement répétable au cours des cycles et une densité énergétique plus élevée (≈ 66 kWh/m3). L’alcool de sucre, l’érythritol (Tfusion = 118 °C), présente de bonnes thermo-physiques (128 kWh/m3) mais la maîtrise de sa cristallisation est un point clé pour l’utiliser en tant que MCP. / This PhD thesis deals with the understanding of the heat transfer mechanisms and with the development of thermal energy storage system for the industrial waste heat recovery application. The use of Phase Change Materials (PCM) is attractive for its high storage density and its possibility to deliver heat at constant temperature. However, the PCM low thermal conductivity leads to develop heat transfer improvement methods, such as heat exchangers with increased heat transfer surface. The goal is to characterize the behavior of such heat exchangers An experimental study, where four several heat exchangers have been tested with different orientations (horizontal/vertical) and injection types (upward/downward), highlighted the impact of natural convection during the melting process and the volume contraction one during the solidification. These results have been validated through a 3D numerical model. A performance comparison method based on an energy calculation through an experimental mesh is proposed and enables to select a heat exchanger on criteria such as the storage density, the characteristic time and the cost. Three PCM, adapted to our application, have been tested at the intended temperature (100-200 °C) by integrating them into a storage system made of a stainless steel tube with aluminum circular fins. Their ability to resist to repeated cycles has been assessed and their behavior has been compared. The salts mixture, H105 (Tmelting = 122 °C), is not selected for the application because of it low storage density (≈ 56 kWh/m3) and its large melting area. The sebacic acid (Tmelting = 132 °C) has a repeatable behavior with cycles and a higher storage density (≈ 66 kWh/m3) and is appropriate as storage material. The sugar alcohol, erythritol (Tmelting = 118 °C), has good thermo-physical properties (128 kWh/m3) but the crystallization control is a key point to use it as a PCM.

Un nouveau modèle de turbulence adapté aux écoulements turbulents soumis à laflottabilité a été développé en utilisant la configuration du canal plan vertical différentiellementchauffé comme référence. L’étude des DNS disponibles pour chacun des régimes de convection amontré les défauts des relations constitutives classiques conduisant à la mauvaise représentationdes écoulements de convection naturelle. Ces modèles ne prennent en compte ni le couplage deschamps thermique et dynamique ni l’anisotropie de l’écoulement, tout deux induits par la flottabilité.Une approche algébrique a donc été utilisées. L’hypothèse d’équilibre local a été validéedans une large partie du canal sauf dans la région de paroi et au voisinage d’un gradient devitesse nul, quel que soit le régime de convection. Les modèles homogènes et pariétaux pour lescorrélations de pression ont été étudiés et sélectionnés. Deux modèles EARSM et EAHFMprenant en compte les termes de flottabilité ont été développés. Ces modèles intègrent aussi destraitements spécifiques à la paroi reposant sur la pondération elliptique. Ils sont couplés à unmodèle corrigé pour mieux représenter le pic d’énergie cinétique turbulente prèsde la paroi. Le modèle complet a été confronté aux DNS sur la configuration du canal pourchacun des régimes de convection à travers des tests a priori et des calculs complets montrantdes résultats très encourageants et de meilleures prévisions que les modèles classiques. / A new turbulence model dedicated to buoyant flows is developped using the differentiallyheated vertical plane channel flow configuration as test case. For each convection regime,the examination of available DNS databases pointed out the failure of classical modeling topredict buoyant flows. Neither the coupling between thermics and dynamics nor the anisotropy,both due to the buoyancy, are considered by these models. So, algebraic models are used. Theweak equilibrium assumption is validated in a large part of the channel except in the wall regionand close to zero velocity gradient whatever the convection regime. The wall and homogeneousmodels for the pressure terms are investigated and selected. Then, an EARSM and an EAHFMare developped to include the buoyant terms. These models both include wall treatments. Theyare coupled with a model modified to improve the representation of the turbulentkinetic energy maximum close to the wall. The complete model is finally compared to theDNS on the channel flow configuration for all convection regime thanks to a priori tests andcomplete computations, showing encouraging results and better predictions than classical models

La caractérisation des circulations thermo-convectives dans les milieux naturels est particulièrement étudiée pour ses applications industrielles, notamment en géothermie. Cette caractérisation correspond à deux enjeux liés à (i) la détection des remontées hydrothermales et (ii) la quantification de l'énergie évacuée par le système actif. La méthodologie proposée par cette thèse comporte ainsi deux volets. Le premier volet vise la localisation des systèmes actifs par caractérisation géophysique de sub-surface. Les techniques de prospection mises en œuvre sont la tomographie de résistivité électrique (TRE), la mesure in situ de température et la mesure de polarisation spontanée (PS). L'objet étudié dans cette thèse est la remontée hydrothermale située au niveau de la limite supérieure de l'effondrement de Rina Grande, au sommet du volcan Stromboli. L'approche géophysique de sub-surface comprend une cartographie spatiale (à l’échelle métrique) en TRE, PS et température et un suivi temporel des signaux PS et thermique. La cartographie spatiale nous renseigne sur l'extension spatiale du système hydrothermal (TRE) et le caractère actif des circulations (PS et température). Le modèle de résistivité électrique a été amélioré en tenant compte de la topographie, issue d'un Modèle Numérique de Terrain généré par photogrammétrie. Le traitement des données de surveillance montre qu'une relation linéaire lie le signal thermique et le signal PS, une relation qui ne peut être purement thermo-électrique. Le jeu de données temporel met l'accent sur la dynamique du système convectif et notamment la probable dépendance de la vigueur du système convectif avec les variations saisonnières de température atmosphérique. Le second volet vise à quantifier l'énergie évacuée par des systèmes thermo-convectifs à l'aide d'une approche numérique. L'approche géophysique met en évidence le fait que les zones endommagées telles que des limites structurales favorisent la circulation des fluides hydrothermaux. Ceci s'explique par le fait que fractures, failles ouvertes, zones d'altération etc. changent à petite échelle la perméabilité du milieu. L'approche numérique est donc vouée à la simulation 3D de l'écoulement et du transfert de chaleur dans des milieux poreux fracturés. Le travail s'est basé sur un code d'écoulement en milieux poreux fracturé pré-existant. Ce code a été adapté afin de résoudre le problème thermique. L'étude théorique vise à quantifier l'influence des paramètres du réseau de fractures sur l'énergie dissipée. L'influence de la densité de fractures, de la transmissivité de fractures et l'anisotropie du réseau de fractures sont évaluées. L'étude compare également les flux de chaleur dégagés par des modèles où les fractures sont insérées de manière discrète avec des modèles homogènes et de même propriétés macroscopiques (approche effective). Les résultats montrent que la validité de l'approche effective est fortement dépendante des paramètres du réseau de fracture. / The characterization of thermo-convective circulations in natural environments is particularly studied especially for geothermal applications. This characterization corresponds to two issues related to (i) the detection of the hydrothermal flows(ii) the quantification of the energy discharged by the active system. The methodology proposed by this thesis thus comprises two parts. The first component aims at locating active systems with geophysical methods. The prospecting techniques used are electrical resistivity tomography (ERT), in situ temperature measurement and self potential (SP) measurements. The object studied in this thesis is the hydrothermal system located at the upper limit of the Rina Grande sector collapse, at the top of the Stromboli volcano. The sub-surface geophysical approach includes spatial (metric) mapping in ERT, SP and temperature, and monitoring of SP and thermal signals. Spatial mapping informs us about the spatial extension of the hydrothermal system (ERT) and the active character of the circulations (SP and temperature). The model of electrical resistivity has been improved taking into account the topography, resulting from a Digital Elevation Model generated by photogrammetry. The treatment of the monitoring data shows that a linear relationship links the thermal and the PS signals, a relationship that can not be entirely due to thermoelectric effect. The temporal dataset focuses on the dynamics of the convective system and in particular the probable dependence of the vigor of the convective system on seasonal variations in atmospheric temperature.The second part aims at quantifying the energy released by thermo-convective systemsusing a numerical approach. The geophysical approach highlights the fact that damaged areas such as structural boundaries favor the circulation of hydrothermal fluids. This is due to the fact that fractures, open faults, alteration zones, etc. change the permeability of the medium on a small scale. The numerical approach is therefore dedicated to the 3D simulation of the flow and heat transfer in fractured porous media. The work is based on a pre-existing code in fracturing porous media, able to solve the flow problem. This code has been adapted to solve the thermal problem also. The theoretical study aims to quantify the influence of fracture network parameters on the released energy. The influence of fracture density, fracture transmissivity and fracture network anisotropy are evaluated. The study also compares the heat fluxes generated by models where the fractures are discretely inserted with homogeneous models with the same macroscopic properties (effective approach). The results show that the validity of the effective approach is highly dependent on the parameters of the fracture network.

Le présent travail apporte une contribution à la compréhension des mécanismes des transferts combinés de chaleur et de masse en convection naturelle lors de l’évaporation d’un film liquide d’eau d’épaisseur négligeable dans un canal vertical ondulé. L’écoulement est laminaire et bidimensionnel. Les équations régissant le phénomène sont résolues à l’aide d’une méthode aux volumes finis et le traitement du couplage vitesse-pression est réalisé par la méthode de projection. Les influences de la densité de flux de chaleur, de la température ainsi que l’humidité de l’air à l’entrée et la forme de la paroi du canal sur les transferts sont étudiées. Les résultats sont présentés sous la forme de ligne de courant, d’isothermes et d’iso-concentrations.Les simulations numériques effectuées ont permis l’étude détaillée de la structure de l’écoulement ainsi que des champs thermiques et massiques. Nous représentons également, les nombres de Nusselt et de Sherwood. / The present work is a contribution to the understanding of the mechanisms of combined heat and mass transfers in natural convection during the evaporation of a liquid film with negligible thickness in a wavy vertical channel. The flow is laminar and two-dimensional. The equations governing the phenomenon are resolved using the finite volumes method and the treatment of the coupling between velocity and pressure is carried out by the projection method. The influences of the heat flux density, the temperature and the humidity of the inlet air and the shape of the channel wall on the transfers are studied. The results are presented in the form of cstreamlines, isotherms and iso-concentrations.The numerical simulations carried out have allowed the detailed study of the flow structure as well as the thermal and mass fields. We also represent the Nusselt and Sherwood numbers.

Des questions fondamentales concernant les propriétés de diffusion des systèmes biologiques dans des conditions isothermes et non-isothermes restent en suspens en raison de l’absence de techniques expérimentales capables de visualiser et de mesurer les phénomènes de diffusion avec une très bonne précision. Il existe en conséquence un besoin de développer de nouvelles techniques expérimentales permettant d’approfondir notre compréhension des phénomènes de diffusion. La convection naturelle en cavité tridimensionnelle inclinée est elle-aussi très peu étudiée. Cette inclinaison de la cavité peut correspondre à un léger défaut expérimental ou être imposée volontairement. Dans cette thèse, nous étudions les phénomènes de transport de chaleur et de masse en cavité parallélépipédique, nous intéressant particulièrement à la thermodiffusion en situation sans convection et à la convection naturelle en fluide pur (sans thermodiffusion). La diffusion de masse est étudiée à l’aide d’une nouvelle technique optique, tandis que la convection naturelle est tout d’abord étudiée en détails avec une méthode numérique sophistiquée, puis visualisée expérimentalement à l’aide du même système optique que pour les mesures de diffusion. Nous présentons l’interféromètre optique de haute précision développé pour les mesures de diffusion. Cet interféromètre comprend un interféromètre polarisé de Mach–Zehnder, un polariseur tournant, une caméra CCD et un algorithme de traitement d’images original. Nous proposons aussi une méthode pour déterminer le coefficient de diffusion isotherme en fonction de la concentration. Cette méthode, basée sur une analyse inverse couplée à un calcul numérique, permet de déterminer les coefficients de diffusion à partir des profils de concentration transitoires obtenus par le système optique. Mentionnons de plus que c’est la première fois que la thermodiffusion est visualisée dans des solutions aqueuses de protéines. La méthode optique proposée présente trois avantages principaux par rapport aux autres méthodes similaires : (i) un volume d’échantillon réduit, (ii) un temps de mesure court, (iii) une stabilité hydrodynamique améliorée. Toutes ces méthodes ont été validées par des mesures sur des systèmes de référence. La technique optique est d’abord utilisée pour étudier la diffusion isotherme dans des solutions de protéines : (a) dans des solutions binaires diluées, (b) dans des solutions binaires sur un large domaine de concentration, (c) dans des solutions ternaires diluées. Les résultats montrent que (a) le coefficient de diffusion isotherme dans les systèmes dilués décroit avec la masse moléculaire, comme prédit grossièrement par l’équation de Stokes-Einstein ; (b) la protéine BSA a un comportement diffusif de type sphère dure et la protéine lysozyme de type sphère molle ; (c) l’effet de diffusion croisée est négligeable dans les systèmes ternaires dilués. La technique optique est aussi utilisée (d) dans des solutions binaires diluées non-isothermes, révélant que les molécules d’aprotinin (6.5 kDa) et de lysozyme (14.3 kDa) sont, respectivement, thermophiliques et thermo-phobiques, quand elles sont en solutions aqueuses à température ambiante. Enfin, la technique optique est utilisée pour l’étude de la convection de Rayleigh-Bénard en cavité cubique horizontale. Puisque la convection peut aussi être étudiée de façon réaliste en utilisant les équations de Navier-Stokes, une analyse numérique de bifurcation est proposée, permettant une étude approfondie de la convection naturelle dans des cavités tridimensionnelles parallélépipédiques. Pour cela, une méthode de continuation a été développée à partir d’un code aux éléments finis spectraux. La méthode numérique proposée est particulièrement bien adaptée aux études de convection correspondant à des diagrammes de bifurcation complexes. [...] / Fundamental questions concerning the mass diffusion properties of biological systems under isothermal and non-isothermal conditions still remain due to the lack of experimental techniques capable of visualizing and measuring mass diffusion phenomena with a high accuracy. As a consequence, there is a need to develop new experimental techniques that can deepen our understanding of mass diffusion. Moreover, steady natural convection in a tilted three-dimensional rectangular enclosure has not yet been studied. This tilt can be a slight defect of the experimental device or can be imposed on purpose. In this dissertation, heat and mass transfer phenomena in parallelepiped enclosures are studied focusing on convectionless thermodiffusion and on natural convection of pure fluids (without thermodiffusion). Mass diffusion is studied with a novel optical technique, while steady natural convection is first studied in detail with an improved numerical analysis and then with the same optical technique initially developed for diffusion measurements. A construction of a precise optical interferometer to visualize and measure mass diffusion is described. The interferometer comprises a polarizing Mach–Zehnder interferometer, a rotating polariser, a CCD camera, and an original image-processing algorithm. A method to determine the isothermal diffusion coefficient as a function of concentration is proposed. This method uses an inverse analysis coupled with a numerical calculation in order to determine the diffusion coefficients from the transient concentration profiles measured with the optical system. Furthermore, thermodiffusion of protein molecules is visualized for the first time. The proposed method has three main advantages in comparison to similar methods: (i) reduced volume sample, (ii) short measurement time, and (iii) increased hydrodynamic stability of the system. These methods are validated by determining the thermophysical properties of benchmark solutions. The optical technique is first applied to study isothermal diffusion of protein solutions in: (a) dilute binary solutions, (b) binary solutions with a wide concentration range, and (c) dilute ternary solutions. The results show that (a) the isothermal diffusion coefficient in dilute systems decreases with molecular mass, as roughly predicted by the Stokes-Einstein equation; (b) BSA protein has a hard-sphere-like diffusion behaviour and lysozyme protein a soft sphere characteristic; and (c) the cross-term effect between the diffusion species in a dilute ternary system is negligible. The optical technique is then applied to (d) non-isothermal dilute binary solutions, revealing that that the aprotinin (6.5 kDa) and lysozyme (14.3 kDa) molecules are thermophilic and thermophobic, respectively, when using water as solvent at room temperature. Finally, the optical technique is applied to study Rayleigh-Bénard convection in a horizontal cubical cavity. Since natural convection can be studied in more depth by solving the Navier-Stokes equations, a bifurcation analysis is proposed to conduct a thorough study of natural convection in three-dimensional parallelepiped cavities. Here, a continuation method is developed from a three-dimensional spectral finite element code. The proposed numerical method is particularly well suited for the studies involving complex bifurcation diagrams of three-dimensional convection in rectangular parallelepiped cavities. This continuation method allows the calculation of solution branches, the stability analysis of the solutions along these branches, the detection and precise direct calculation of the bifurcation points, and the jump to newly detected stable or unstable branches, all this being managed by a simple continuation algorithm. This can be used to calculate the bifurcation diagrams describing the convection in tilted cavities. [...]

Dans le contexte énergétique actuel, le développement de l'énergie solaire passe par son développement à grande échelle en milieu urbain. Les solutions actives et intégrées telles que les double-façades photovoltaïques permettent d'une part d'envisager une valorisation de la chaleur et d'autre part de valoriser l'enveloppe des bâtiments en tant que surface de captation. De plus, cette configuration limite la perte de rendement et l'accélération du vieillissement liées à la surchauffe des cellules. En effet, les composants photovoltaïques sont séparés du bâtiment par une lame d'air et la convection naturelle développée dans cet espace favorise l'évacuation de la chaleur absorbée par les panneaux. Cependant, le comportement global d'une double-façade photovoltaïque est très complexe et encore mal compris. Cette thèse se focalise sur l'expérimentation du phénomène de la convection naturelle. Un banc d'essais, développé dans le cadre de la thèse de Daverat, modélise la double-façade par un canal vertical avec chauffage pariétal. Le fluide d'étude est de l'eau afin de s'affranchir du rayonnement entre parois et d'étudier le couplage conduction-convection. Le banc d'essais est constitué d'un canal vertical de 65cm de haut, chauffé à une densité de flux uniforme, plongé dans une cellule de 1,5m de haut remplie d'eau. L'écoulement est observé par des mesures de densité de flux de chaleur, de température pariétale et de vitesse, ces dernières étant réalisées par vélocimétrie par images de particules et par vélocimétrie laser Doppler. Deux configurations de chauffage sont étudiées : symétrique et asymétrique. Pour l'étude de la configuration avec un chauffage symétrique, les données expérimentales issues du travail de thèse de Daverat sont analysées par une approche zonale. Un découpage du canal est proposé et une analyse des ordres de grandeur permet de comprendre les observations expérimentales. Ainsi, cette analyse se focalise sur le transfert thermique aux parois, l'évolution de la pression dans le canal et sur les fluctuations, en particulier, sur l'évaluation du terme croisé du tenseur de Reynolds. La configuration d'un chauffage asymétrique est étudiée expérimentalement pour la phase de démarrage et le régime stationnaire. Pour la phase de démarrage, la caractérisation des premières heures à partir de la mise en chauffage de la paroi est réalisée en termes de température pariétale, de vitesse et de comportement bi-/tri-dimensionnel de l'écoulement. Un écoulement de retour de grande échelle est mis en évidence. L'impact de la stratification thermique externe du canal est également étudié. Pour le régime stationnaire, les analyses sont réalisées sur des données pour lesquels le régime est considéré établi. Les comportements thermique et dynamique sont étudiés et mettent en évidence deux états distincts. Des écoulements de retour sont également observés. Enfin, les influences du nombre de Prandtl et de l'écartement du canal sont discutées. / Under the current energy context, the development of solar energy goes through its large scale development, especially in urban areas. Active and integrated solutions, such as photovoltaic double-skins, allow both the heat recovery for building needs and the exploitation of building envelope as collecting surface. Furthermore, this configuration limits the efficiency loss and the acceleration of ageing process due to overheating of the solar cells. Indeed, in this configuration, the photovoltaic panels are separated from the building by an air channel and the thermally driven flow that develops in the channel helps to evacuate the heat absorbed by the panels. However, the global behaviour of the photovoltaic double-skin is complex and not yet fully understood. This study is part of a scientific project on multi-scale modelling and experimentation of BIPV components. It focuses on the experimentation of the natural convection phenomenon in double-façades. An experimental apparatus developed during C. Daverat's thesis models the double-façade by a vertical channel with wall heating. The working fluid is water so as to avoid radiative heat transfer between walls. The experimental apparatus consists of a 65 cm high vertical channel, with isoflux heating, immersed in a 1.5 m high water tank. The fluid flow is instrumented with measurements of heat flux, wall temperature and velocity. The velocity measurements are made by Particle Image Velocimetry and Laser Doppler Velocimetry. Two heating configurations are studied: symmetrical (same heat flux is imposed on both walls) and asymmetrical (only one wall is heated, the other one being adiabatic). For the study with symmetrical heating configuration, experimental data from Daverat's thesis are analysed by a zonal approach. The channel is split into several zones and a scaling analysis is conducted in order to explain experimental observations. Thus, the study focuses on heat transfer, pressure evolution in the channel and fluctuations, in particular, the evaluation of the maximum Reynolds stress. The asymmetrical configuration is studied for the transient regime and the steady regime. For the transient regime, early stage (first hours) is characterised in terms of wall temperature, velocity and two- or three-dimensional flow behaviour. A large-scale reversal flow is observed. The impact of external thermal stratification is also studied. For the steady regime, thermal and dynamic behaviours are studied and two different states are distinguished and characterised. Reversal flows are also observed. Finally, the influence of the Prandtl number and the channel width is discussed.

La réduction des émissions de gaz à effet de serre (GES) passe par la réduction des consommations d’énergie. Le stockage de la chaleur dans les parois des bâtiments permet de réduire la consommation d'énergie. Parmi les techniques de stockage, le stockage latent a la capacité de stocker une quantité d’énergie par unité de volume plus importante qu’un système sensible. Le projet INERTRANS a proposé le développement d’une façade associant une isolation translucide et le stockage latent avec un matériau à changement de phase (MCP). La fusion du MCP s’accompagne de la convection naturelle et l’absorption ou transmission du rayonnement. Le couplage de l’ensemble de ces phénomènes n’a pas été étudié dans la littérature. Dans cette thèse un modèle numérique 2D pour l’étude du changement de phase a été développé. Ce modèle utilise la méthode de Boltzmann sur réseau (LB) à temps de relaxation multiple (MRT), pour la résolution du champ de vitesse et la méthode des différences finies, pour la résolution du champ de températures. Le changement de phase a été traité par la formulation enthalpique. L’originalité est l’application de ce modèle au problème de changement de phase avec convection naturelle, d’une part, et au changement de phase avec convection naturelle et rayonnement, d’autre part. Pour vérifier notre modèle sans rayonnement, un cas de référence de la littérature a été simulé. Il s’agit de la fusion des deux MCP, l’étain et l’octadécane, à faible et fort nombre de Prandtl, respectivement. La simulation de l’étain a confirmé un écoulement multicellulaire. La simulation de l’octadécane a montré une forte influence de la convection avec un front de fusion qui se déforme sur toute la cavité. Le nombre de Nusselt pour l’octadécane avec convection est plus de trois fois le Nusselt sans convection. La simulation de l’acide gras de la brique INERTRANS a montré que la convection ne doit pas être négligée, car le flux prédit avec convection peut être jusqu’à trois fois plus grand que le flux prédit sans convection. La fraction fondue est près du double qu’en conduction seule. La méthode LB appliquée aux transferts radiatifs a été étudiée. Il se trouve, qu’à l’état actuel cette méthode n’est pas compétitive par rapport à une méthode classique des ordonnées discrètes (MOD). Enfin, nous avons couplé la MOD pour le calcul du flux radiatif avec la méthode LB pour le calcul du champ de vitesses et des différences finies pour l’équation de l’énergie. Le rayonnement grande longueur d’onde n’a pas d’influence notable sur les transferts thermiques. Le rayonnement courte longueur d’onde augmente les transferts thermiques, pourtant, cet effet n’est pas aussi important que l’augmentation due à la convection pour le matériau choisi. Puisqu’aucune solution de référence n’existe dans la bibliographie, nos résultats peuvent désormais servir d’éléments de comparaison pour de futurs travaux. Une validation expérimentale constituerait une perspective nécessaire. / Reduction of greenhouse gas emissions requires reduction of energy consumption. Energy storage on building walls allows reduction in energy consumption. Among storage techniques, latent heat storage offers higher energy storage density than sensible heat storage. INERTRANS project has proposed the development of an innovative facade, coupling transparent insulation and energy storage with a fatty acid phase change materials (PCM). Melting of PCM comprises different phenomena, namely, natural convection in the liquid phase and radiation absorption or transmission. The coupling of all this phenomena is not still studied in scientific literature. In this thesis, a 2D numerical model for studying phase change has been developed. This model uses the lattice Boltzmann method (LBM) with multiple relaxation time (MRT) to resolve velocity field, and finite differences for the temperature field. Phase change is treated with the enthalpy formulation. The original contribution is application of this hybrid approach to the phase change with natural convection, on the one hand, and to the phase change with natural convection and radiation, on the other hand. To verify the model without radiation, a test case taken from literature has been simulated. It concerns the melting of two PCM with a low and high Prandtl number, the tin and octadecane, respectively. Tin melting simulation confirms multiple cells flow, starting with four rolls which merges in three then two rolls. Octadecane simulation shows high convection effect, with a melting front deforming all along the cavity height. Nusselt number plot for octadecane melting with convection is more than three times with conduction only. INERTRANS’ fatty acid simulation shows that convection shall not be neglected, because predicted heat flux with convection may be up to three times that predicted with conduction only. Melted fraction is almost twice than with conduction only. The lattice Boltzmann method applied to radiative heat transfer has also been explored. It turns out that in its current state, this method is not competitive compared to a conventional discrete ordinates method (DOM). Finally, we coupled the DOM for radiation heat flux, with the LBM for velocity field calculation and finite differences for the energy equation to solve the coupling between phase change, convection and radiation. Long wavelength radiation has no noticeable effect on heat transfer. Short wavelength radiation increases heat transfer, however, this increase is not as important as that produced by convection for this kind of material. Since no reference solution exists in the literature, our results can now serve as a basis for future work. An experimental validation would be a necessary perspective.

Notre investigation porte sur la simulation numérique des échanges thermo-convectifs dans un canal vertical ouvert à flux imposé. Cette étude rentre dans le cadre des recherches sur le rafraîchissement passif des composants PV intégrés au bâtiment. À cet effet, un code numérique en Différences Finies est utilisé pour résoudre les équations de Navier-Stokes et simuler la convection naturelle dans un canal. Ce problème reste difficile à résoudre parce que l'écriture des conditions aux limites d'entrée et de sortie reste un problème ouvert. Notre travail consiste d'abord en étude des différentes conditions aux limites pour le benchmark numérique AMETH. Les travaux réalisés ont permis de faire un premier choix sur les conditions aux limites. L'étude s'oriente ensuite sur la qualification et la quantification numériques et expérimentales pour deux fluides : l'air (convection-rayonnement) et l'eau (convection pure). Les résultats numériques/expérimentaux ont été comparés et les discordances ont été analysées. Plusieurs aspects phénoménologiques (rayonnement entre surfaces, variation des propriétés thermo-physiques, variation du nombre de Prandtl) ont été abordés afin de caractériser leurs influences respectives sur l'écoulement et le transfert thermique. Enfin, dans le but d'apporter des éléments de réponses sur les conditions aux limites dynamiques, nous avons simulé la convection naturelle d'un canal dans une cavité et tenté une modélisation. / The present investigation deals with natural convection flow in a vertical open-ended channel with wall constant heat flux. This study falls under the framework of research on passive cooling of building integrated PV components. For this purpose, a numerical code developed with Finite Differences scheme is used to solve Navier-Stokes equations and simulate the natural convection in a channel. This problem is difficult to solve because the writing of inlet/outlet boundary conditions remains an open problem. First, our work consists of studying different boundary conditions for the the numerical benchmark AMETH. The work carried out has enabled a first choice of boundary conditions. The study then focuses on numerical and experimental quantification and qualification for two fluids : air ( convection - radiation) and water ( pure convection) . Experimental and numerical results were compared and discrepancies were analyzed. Several phenomenological aspects ( surface radiation, thermophysical properties variation, Prandtl number variation ) were discussed in order to characterize their influence on flow and heat transfer. Finally, in order to provide some answers on dynamical boundary conditions, we simulated natural convection of a channel inside a cavity and tried a modeling.