Le travail réalisé en collaboration entre l’entreprise CNIM et le laboratoire IUSTI a permis la miseen place d’un outil de simulation numérique afin d’étudier à l’échelle locale différents phénomènescouplés qui se produisent dans un réacteur de gazéification en lit fluidisé dense. L’approche choisie,DEM (Discrete Element Method), est basée sur le suivi de paquets de particules. Les résultats desimulation ont été comparés à des mesures expérimentales réalisées dans les laboratoires IUSTIet LERMAB : mesure de pression dans un lit fluidisé peu profond et caractérisation de la réactionde pyrolyse par suivis temporels de la masse volumique par une méthode innovante et de latempérature de la particule de biomasse ainsi que la composition des gaz produits. À plus grandeéchelle, les prédictions numériques du code ont été comparées à des mesures expérimentales deprofil de vitesse des particules et de mélange et ségrégation issues de la littérature. Les résultatsnumériques du lit fluidisé bisolide sans et avec réactions chimiques, obtenus sont en bon accordavec les mesures expérimentales. Ils ont permis de mettre en évidence le fort couplage entre lesphénomènes hydrodynamiques et thermochimiques.Ces travaux ouvrent de nouvelles perspectives tant sur le plan de l’expérimentation que de la simulation numérique où l’approche DEM retenue a montré un fort potentiel. L’extrapolation dumodèle DEM pour la simulation de la gazéification à des réacteurs industriels reste encore un défiau regard des ressources informatiques. / Energy production from green and renewable resources, such as biomass, are currently experiencinga significant growth. Thermochemical conversion of this biomass by gasification is a process usedfor over a century but still requires significant developments in terms of rentability optimizationand quality improvement of products gases.The work carried out in collaboration between the CNIM company and the IUSTI laboratoryallowed the establishment of a numerical simulation tool to study locally different coupled phenomenaoccurring in a dense fluidized bed gasification reactor. The chosen approach, DEM (DiscreteElement Method), is based on the monitoring of particle packets. The simulation results werecompared to experimental measurements realised in IUSTI and LERMAB laboratories : measuringpressure in a shallow fluidized bed and characterization of pyrolysis reaction with the timetracking of particle density by an innovative method, temperature and the composition of theproducts gases. On a larger scale, the numerical code predictions were compared with velocity,particles mixing and segregation profiles from experimental measurements of the literature. Theobtained numerical results of bi-solid fluidized bed with and without chemical reactions are ingood agreement with the experimental measurements. It helped to highlight the strong couplingbetween hydrodynamic and thermochemical phenomena.This work opens up new perspectives on the experimental plan and numerical simulation whichDEM approach has shown great potential. The extrapolation of the DEM model for the simulationof gasification industrial reactors remains a challenge in terms of computer resources.

Ces travaux portent sur la résolution successive de deux problèmes inverses en transferts thermiques : l'estimation de la densité de flux en surface d'un matériau puis de la conductivité thermique équivalente d'une couche déposée en surface de ce matériau. Le modèle direct est bidimensionnel orthotrope (géométrie réelle d'un matériau composite), instationnaire, non-linéaire et ses équations sont résolues par éléments finis. Les matériaux étudiés sont les composants face au plasma (tuiles composite carbone-carbone) dans le Tokamak JET. La densité de flux recherchée varie avec une dimension spatiale et avec le temps. La conductivité du dépôt de surface varie spatialement et peut également varier au cours du temps pendant l'expérience (toutes les autres propriétés thermophysiques dépendent de la température). Les deux problèmes inverses sont résolus à l'aide de l'algorithme des gradients conjugués associé à la méthode de l'état adjoint pour le calcul exact du gradient. La donnée expérimentale utilisée pour la résolution du premier problème inverse (estimation de flux surfacique) est le thermogramme fourni par un thermocouple enfoui. Le second problème inverse utilise, lui, les variations spatio-temporelles de la température de surface du dépôt inconnu (thermographie infrarouge) pour identifier sa conductivité. Des calculs de confiance associée aux grandeurs identifiées sont réalisés avec la démarche Monte Carlo. Les méthodes mises au point pendant ces travaux aident à comprendre la dynamique de l'interaction plasma-paroi ainsi que la cinétique de formation des dépôts de carbone sur les composants et aideront au design des composants des machines futures (WEST, ITER). / This work deals with the successive resolution of two inverse heat transfer problems: the estimation of surface heat flux on a material and equivalent thermal conductivity of a surface layer on that material. The direct formulation is bidimensional, orthotropic (real geometry of a composite material), unsteady, non-linear and solved by finite elements. The studied materials are plasma facing components (carbon-carbon composite tiles) from Tokamak JET. The searched heat flux density varies with time and one dimension in space. The surface layers conductivity varies spatially and can vary with time during the experiment (the other thermophysical properties are temperature dependent). The two inverse problems are solved by the conjugate gradient method with the adjoint state method for the exact gradient calculation. The experimental data used for the first inverse problem resolution (surface heat flux estimation) is the thermogram provided by an embedded thermocouple. The second inverse problem uses the space and time variations of the surface temperature of the unknown surface layer (infrared thermography) for the conductivity identification. The confidence calculations associated to the estimated values are done by the Monte Carlo approach. The method developed during this thesis helps to the understanding of the plasma-wall interaction dynamic, as well as the kinetic of the surface carbon layer formation on the plasma facing components, and will be helpful to the design of the components of the future machines (WEST, ITER).

Le Cycle Organique de Rankine (abrégé ORC de l’anglais Organic Rankine Cycle) est une technologie permettant la conversion de chaleur basse température en électricité. L’ORC transcritique a été identifié comme une solution prometteuse pour la valorisation de la chaleur fatale. Cependant, peu d’installations expérimentales ont permis de confirmer ces performances. Ce travail de thèse présente le fonctionnement et l’optimisation d’ORC sous-critique et transcritique pour la conversion de chaleur basse température en électricité à partir de différentes sources. Premièrement, les contextes thermodynamique et technologique de l’ORC sont présentés. Des critères de performance énergétiques et exergétiques sont définis et appliqués à une base de données d’installations expérimentales afin d’exposer l’état de l’art actuel des ORC. Deuxièmement, les outils numériques et expérimentaux, spécifiquement développés ou utilisé pour ces travaux, sont présentés. Trois installations expérimentales d’ORC transcritique complet ou incomplet fournissent les données expérimentales. Différents modèles numériques sont utilisés : sous l’environnement Matlab pour la modélisation en permanent, l’analyse des données expérimentales et l’analyse énergétique/exergétique ; L’environnement Modelica/Dymola pour l’analyse des transitoires et de la dynamique du système. Dans un troisième temps, ces différents outils sont utilisés pour étudier quatre différentes problématiques : - Le fonctionnement de la pompe de circulation est étudié, d’un point de vue énergétique et volumétrique. Des modèles semi-empiriques et des corrélations de performance sont présentés. - Les transferts thermiques en supercritique sont examinés, en local et en global. Les coefficients de transfert thermique sont comparés avec différentes corrélations de la littérature. - L’influence de la charge de réfrigérant sur les performances et le comportement de l’ORC est analysée. La charge optimale est estimée pour différentes conditions de fonctionnement et des mécanismes de régulation de la charge sont présentés. - Les performances énergétiques et exergétiques de l’ORC sont comparées avec la base de données. Une analyse exergétique du procédé a permis d’identifier des voies d’amélioration. / The Organic Rankine Cycle (ORC) is a technology used for low-grade thermal energy conversion into electricity. Transcritical ORC has been identified as a solution for efficient waste heat recovery. However, few experimental tests have been conducted to confirm the interest of transcritical ORC and investigate its operational behaviors. The work presented focuses on the operation and the optimization of subcritical and transcritical Organic Rankine Cycles for low-grade heat conversion into electricity from various heat sources (solar, industrial waste heat). First, the thermodynamic framework of ORC technology is presented. Energetic and exergetic performance criteria, appropriate to each type of input source, are introduced and selected. The criteria are later applied to a database of ORC prototypes, in order to objectively analyze the state-of-the-art. In a second step, the experimental and numerical tools, specifically developed or used in the present thesis, are presented. Three subcritical and transcritical ORC test benches (hosted by CEA and AUA) provided experimental data. Numerical models were developed under different environments: Matlab for steady-state modeling, data processing and energy/exergy analysis. The Modelica/Dymola environment for system dynamics and transient operations. Lastly, the different tools are exploited to investigate four different topics: - The ORC pump operation is investigated, both under an energetic and volumetric standpoint, while semi-empirical models and correlations are exposed. - Supercritical heat transfers are explored. Global and local heat transfer coefficients are estimated and analyzed under supercritical conditions, while literature correlations are introduced for comparison. - Working fluid charge influence over the ORC performance and behavior is investigated. Optimal fluid charge is estimated under various operating conditions and mechanisms for charge active regulation are exposed. - ORC system performances and behavior are discussed. Through both an energetic and exergetic standpoint, performances are compared with the state-of-the-art, while optimization opportunities are identified through an exergetic analysis.

Ce travail de thèse, réalisé dans le cadre d'une convention CIFRE entre TURBOMECA et le CERFACS et en partenariat avec l'IVK, se place dans un contexte d'amélioration des performances des turbines axiales équipant les turboréacteurs d'hélicoptère. Un des points critiques du dimensionnement de tels moteurs est la maitrise de la durée de vie des pales de la turbine haute pression qui font face à de très hautes températures provenant de la chambre de combustion. Les prédictions numériques de l'environnement aérothermique des pales (écoulements dans la veine et système de refroidissement) sont réalisées aujourd'hui dans le milieu industriel à l'aide de la modélisation Reynolds Averaged Navier-Stokes (RANS). Grâce à des capacités de calculs grandissantes, l'approche Simulation aux Grandes Echelles (SGE) offre désormais un nouveau potentiel de prédictions d'écoulements. Les travaux de cette thèse s'intéressent ainsi à la capacité de la SGE à prédire l'écoulement du circuit de refroidissement interne d'une pale de turbine. Pour simplifier l'analyse de ce problème ou plusieurs phénomènes physiques sont en jeu, une progression en trois parties est proposée. La première s'intéresse à l'étude aérothermique de géométries simplifiées de canaux de refroidissement (coude à 180° et canal avec promoteurs de turbulence) en configuration statique. Aux régimes d'écoulement considérés, une approche résolue en paroi avec maillage non-structuré hybride est proposée et validée en vue d'une application industrielle facilitée. La seconde partie étend l'analyse de l'écoulement à un cas de canal avec promoteurs de turbulence en rotation utilisant une méthode de résolution numérique dans un repère absolu. Les investigations des résultats de la SGE fournissent des prédictions moyennes et instationnaires en bon accord avec les expériences disponibles et les travaux précédents aussi bien pour la dynamique de l'écoulement que les transferts de chaleur. Enfin, une troisième partie présente une application de la méthode sur un cas de pale réelle avec couplage thermique entre le circuit de refroidissement et le solide de la pale. Cette dernière partie classée confidentielle n'est pas présente dans le manuscrit disponible publiquement. Les résultats de l'approche résolue en paroi et de la rotation dans le repère absolu comparés aux résultats RANS disponibles pour le cas applicatif montrent d'importante différences locales et ainsi le potentiel de la méthode proposée.

La réduction de la consommation de carburant et des émissions de dioxyde de carbone par les transports, notamment automobiles, conduit au développement de groupes motopropulseurs (GMP) plus économes. Parmi les technologies existantes, l'hybridation électrique des moteurs à combustion interne s'est fortement développée ces dernières années. La minimisation de la consommation d'un GMP de type hybride électrique requiert notamment la connaissance de l'impact des différents transferts thermiques y ayant lieu, notamment après un démarrage à température ambiante. Dans ce contexte, la thèse présentée dans ce mémoire a pour objet l'étude par simulation des transferts thermiques au sein des différents organes du GMP (moteur à combustion interne, moteur électrique, batteries). Elle repose tout d'abord sur un examen bibliographique détaillé des différentes architectures hybrides électriques, des transferts thermiques et les modèles associés, et des différentes stratégies de fonctionnement. Le développement d'un modèle zéro-dimensionnel d'un GMP hybride prenant en compte les différents transferts thermiques est ensuite présenté. Les choix des différents modèles de transferts thermiques sont présentés et justifiés, en termes de précision et de temps de calcul. Pas la suite, l'impact des transferts thermiques est étudié grâce au simulateur, notamment en termes de consommation instantanée et cumulée sur différents cycles de conduite. L'impact des transferts thermiques pour différents niveaux d'hybridation est également quantifié. Finalement, plusieurs études paramétriques et prospectives sont réalisées à la fin de ce mémoire, permettant d'envisager des moyens de réduction de la consommation de carburant.

L’ébullition est un mode de transfert de chaleur intervenant dans de nombreux systèmes thermiques ou énergétiques de par son efficacité. Dans certains, elle se produit à basse pression. La pression statique induite par la colonne de liquide au-dessus de la surface de formation des bulles n’est alors pas négligeable devant la pression de saturation à la surface libre. Dès lors, la pression et le sous-refroidissement induit ne peuvent plus être considérés homogènes autour des bulles, d’où des inhomogénéités des propriétés thermophysiques dans le fluide. Les influences relatives des forces s’exerçant sur une bulle pendant sa croissance sont modifiées par rapport aux pressions plus élevées : une dynamique de bulles différente apparaît. Ces conditions particulières affectent également les transferts thermiques. L’influence de la pression sur l’ébullition en vase a été étudiée expérimentalement dans le régime de bulles isolées en site unique. L’eau a d’abord été testée sur une large gamme de pressions subatmosphériques. Quatre comportements de dynamique de bulles ont été identifiés d’après la visualisation par caméra rapide. Plusieurs paramètres de la dynamique ont été quantifiés grâce à un traitement d’images adapté appliqué aux vidéos. Pour généraliser le concept d’ébullition à « basse pression » et mieux en appréhender les phénomènes fondamentaux, de nouveaux essais ont été réalisés avec un second fluide, le cyclohexane, choisi pour sa similitude thermodynamique avec l’eau bouillant en deçà de la pression atmosphérique. La comparaison des comportements des deux fluides a permis d’identifier certains paramètres responsables des spécificités du phénomène. En outre, de nouvelles fonctionnalités sont apportées au dispositif expérimental pour – notamment – effectuer la mesure rapide du flux transféré sous la bulle pendant sa croissance, synchroniser ces mesures thermiques avec l’acquisition d’images et étudier des surfaces d’ébullition structurées. Les résultats obtenus sont encourageants pour l’analyse des comportements spécifiques de l’ébullition à basse pression et ses applications. / Boiling is an efficient heat transfer mode used in numerous thermal or energy systems. In some systems boiling takes place at low pressure. The static head of the liquid column over the wall where bubbles nucleate is then not negligible against the saturation pressure at the free surface level. The pressure and the induced subcooling degree therefore cannot be considered as homogeneous around growing bubbles, resulting in non-homogeneous thermophysical properties in the fluid. The relative influence of the forces acting on a growing bubble differs from higher pressure conditions, yielding specific bubble dynamics features. Heat transfer is consequently also affected. The effect of the pressure on pool boiling was experimentally investigated during the isolated bubbles regime taking place from a single activated nucleation site. Experiments were first conducted with water for a wide range of subatmospheric pressures. Four distinct bubble dynamics behaviors were identified through high-speed camera visualizations. An adapted image processing of the recordings enabled the measurement of several bubble dynamics characteristics. In order to generalize the concept of pool boiling at "low pressure" and to get a better understanding of the related fundamental phenomena, new experiments were performed with a second fluid, cyclohexane, chosen from original thermodynamic similarity with water boiling at pressures lower than atmospheric. The comparison of fluids’ behaviors made possible the identification of parameters governing the specific phenomena occurring during boiling at low pressure. Besides, the experimental facility was improved to provide new functionalities. The high-speed measurement of the heat flux transferred under the growing bubble, its synchronization with the high-speed videos images and the study of boiling on enhanced surfaces are in particular made possible. Results are encouraging for a better understanding of the specific behaviors of low pressure boiling and for its future implementation in practical applications.

La gestion thermique des batteries Li-ion pour le véhicule électrique est essentielle, pour assurer une autonomie et une durée de vie optimales de ces batteries. Habituellement, des circuits d'air ou de liquide de refroidissement sont utilisés comme systèmes de gestion thermique. Cependant, ces systèmes sont coûteux en termes d'investissement et d'exploitation et doivent être dimensionnés sur la puissance maximale à extraire. L'utilisation de matériaux à changement de phase (MCP) pour l’absorption sous forme de chaleur latente de la chaleur à dissiper peut représenter une alternative moins coûteuse et plus facile à utiliser. En effet, les MCP peuvent stocker passivement la chaleur excédentaire produite et être utilisés en tant que systèmes passifs. Cependant, les MCP présentent de nombreux inconvénients comme la difficulté de décharger l’énergie thermique stockée, ce qui limite l’aptitude du système au cyclage, ou encore leur conductivité thermique peu élevée qui limite les capacités d’échange. Pour résoudre le problème de la régénération des MCP, un système actif supplémentaire peut être ajouté, dimensionné sur une puissance modérée; l'ensemble devient alors un système semi-passif. Dans cette étude, un système de gestion thermique composé d'un MCP et d’air en convection forcée est évalué. Ce système permet de coupler les avantages de ces deux techniques. Une modélisation du système est développée pour une cellule de batterie. Une comparaison avec de l’air uniquement, en convection forcée, montre l'utilité du MCP. Pour augmenter la capacité d’échange du MCP, un matériau à haute conductivité thermique peut être ajouté au MCP, ce qui permet d’obtenir un composite ayant une conductivité thermique plus élevée. Des composites basés sur les MCP étudiés et des nanostructures de carbone sont élaborés, leur conductivité thermique est mesurée. Ensuite, un système expérimental simulant la dissipation d’une cellule de batterie est construit et utilisé pour évaluer le MCP seul, le MCP inclus dans une mousse métallique et le meilleur composite obtenu. Enfin, pour se rapprocher des conditions réelles, un modèle représentant un stack entier de batterie est développé, des simulations sont produites et les résultats obtenus sont commentés. / Li-ion battery thermal management is essential for electric vehicles (EVs), to ensure an optimal autonomy and lifespan of those batteries. Usually, air or coolant circuits are employed as thermal management systems. However, those systems are expensive in terms of investment and operating costs and must be dimensioned on the maximal power to be extracted. The use of phase change materials (PCMs) as latent heat storage medium allowing the absorption of the heat to be dissipated as latent heat may represent an alternative cheaper and easier to operate. In fact, PCMs can passively store the excess heat produced by a device and be used as passive systems. However, PCMs have several drawbacks like the difficulty to discharge the stored thermal load which limits the system’s cyclability or their low thermal conductivity which limits their heat transfer capacity. To solve the problem of the PCM regeneration, an additional active system can be added, dimensioned on a moderate power; the whole becomes a semi-passive system. In this study, a thermal management system composed of a PCM and forced air convection is evaluated. This system permits to combine the respective advantages of the two techniques. A model of the system is developed for one battery cell. A comparison with forced air convection only points out the usefulness of the PCM. To overcome the PCM low thermal conductivity, a highly conductive material can be added to the PCM permitting to obtain a composite with a higher thermal conductivity. Composites based on the PCMs studied and carbon nanostructures are elaborated, and their thermal conductivity is measured. Then, an experimental system permitting to simulate the dissipation of a battery cell is build and used to evaluate the PCM alone, the PCM embedded in metal foam and the better obtained composite. Finally, to be closer to the real conditions, a model representing an entire battery stack is developed, simulations are produced and the obtained results are discussed.

L'utilisation des meubles frigorifiques de vente fermés a augmenté régulièrement dans les supermarchés en raison des économies d'énergie réalisées par rapport aux meubles ouverts. D’où la nécessité d'élargir les connaissances scientifiques à ce type de meuble. Pour l’instant la plupart des études de la littérature ont porté sur l'amélioration de l'efficacité énergétique après l'installation de portes, tandis que les études sur les mécanismes de transfert de chaleur dans les meubles frigorifiques fermés sont rares. De plus, il convient d’étudier également les circulations d’air car ils influencent les échanges de chaleur entre l’air et les produits, donc leur température. Cette thèse a pour objectif de mieux comprendre les phénomènes d’écoulement d’air et de transfert de chaleur dans les meubles frigorifiques fermés par la mise en œuvre d'approches expérimentales et numériques.Des investigations expérimentales ont été menées sur un meuble frigorifique fermé placé dans une cellule d’essai à température contrôlée. Quatre-vingts thermocouples calibrés répartis dans le meuble ont permis d’observer les évolutions spatio-temporelles des températures de l’air et du produit dans différentes conditions de fonctionnement. Nous avons fait varier la température de l'air ambiant (15, 19, 24 et 29 ° C), le volume occupé par le produit (meuble vide ou plus ou moins chargé de blocs de méthylcellulose) et la fréquence d'ouverture des portes (0 : fermé en permanence, 10, 20, 40, 60 Ouvertures Par Heure - OPH). Des expériences ont également été réalisées dans une configuration ouverte du meuble frigorifique (portes complètement retirées) afin de déterminer les avantages des portes sur les performances thermiques.Les résultats ont montré que lorsque les portes étaient fermées en permanence, les champs de température étaient similaires quels que soient la température ambiante et le pourcentage de volume occupé : la température la plus élevée se situe à l’avant de l'étagère supérieure alors que la température la plus basse est observée à l’arrière de l'étagère inférieure. Cette répartition change lorsque les portes sont périodiquement ou définitivement ouvertes : la position où la température est la plus élevée migre vers l’avant de l'étagère du milieu. Même pour une fréquence très élevée d’ouverture des portes (60 OPH), les températures de l’air et des produits dans le meuble fermé sont restées plus basses d’au moins 1,0 ° C par rapport au meuble sans portes.Des mesures de vitesse d'air par anémométrie à fil chaud au niveau du rideau d’air ont permis d'observer son allure générale tandis que les mesures effectuées dans le conduit arrière ont permis de quantifier la répartition du flux d'air à travers le plaque arrière perforée. L’utilisation d’une technique de vélocimétrie par image de particules (PIV) a permis de caractériser le rideau d’air avec une résolution et une précision spatiale supérieures. Les résultats ont notamment montré l’existence d’une zone de recirculation de l'air dans la partie supérieure du meuble où l’on observe également une infiltration d'air ambiant chaud par les fentes présentes autour des portes. Ceci induit une augmentation de la température du rideau d'air. Des simulations numériques bidimensionnelles d’écoulement (CFD) avec le modèle de turbulence k-ε ont permis de reproduire les principaux phénomènes d'écoulement observés par PIV et de voir leur influence sur la distribution de température dans le meuble.Enfin, un modèle simplifié des transferts de chaleur a été développé par une approche zonale en régime permanent et en régime transitoire ce qui permet respectivement de prévoir les températures moyennes de l’air et du produit et les fluctuations de température en fonction des cycles marche / arrêt du compresseur. Le modèle transitoire a été résolu avec une approche spectrale. L'influence de différents paramètres sur les niveaux de température et l'amortissement des fluctuations peut ainsi être identifiée. / The use of closed refrigerated display cabinets in supermarkets has been increased steadily because of the potential energy savings compared to open ones. This growing trend has contributed to the necessity to expand research in the field of retail refrigeration. Most studies in literature, however, focused on the improvements of energy efficiency after door installation while studies on the mechanism of heat transfer and airflow within closed display cabinets are still limited. In fact, the airflow pattern influences the heat exchange between air and products, thus, product temperature. This PhD thesis aims to gain an insight into the mechanism of airflow and heat transfer in closed refrigerated display cabinets by the implementation of experimental and numerical approaches.Experimental investigations were conducted in a closed refrigerated display cabinet (an integral type with a single band air curtain and two double-glazing doors) located in a controlled-temperature test room. Air/product temperatures and air velocity are the main parameters taken into investigations. Eighty calibrated thermocouples distributed throughout the cabinet made it possible to observe the spatial and temporal evolutions of the air and product temperatures under different operating conditions. These conditions were ambient air temperature (15, 19, 24 and 29 °C), product-occupied volume (unloaded, half-loaded and full-loaded with test packages made of methylcellulose), door-opening frequency (0 - permanently closed, 10, 20, 40, 60 Openings Per Hour - OPH) and opening duration (15s and 30s). An automatic door opening system was developed and allowed to apply the opening regime as prescribed in the standard test (EN ISO 23953-2, 2015). The experiment was also conducted in an open configuration of the cabinet (doors were completely removed) to determine the benefits of the doors on the temperature performance. The results showed that when the doors were permanently closed, the temperature distribution in the cabinet was similar whatever the ambient temperatures and occupied-volume percentages – the highest temperature position at the front-top shelf and the lowest temperature position at the rear-bottom shelf. The temperature distribution changed when the doors were periodically or permanently open – the front of the middle shelf became the highest temperature position while the lowest temperature position remained at the rear-bottom shelf. However, the air and product temperatures in the cabinet with doors remained lower despite a very high door-opening frequency (i.e. 60 OPH, product temperatures at least 1.0 °C lower), compared to the case without doors.Air velocity measurement using a hot-wire anemometer at the front of the cabinet from the discharge to the return air grilles allowed to observe the shape of the air curtain, while the measurement in the rear duct allowed to quantify the air flow distribution over the perforated back panel of different shelves. The use of a Particle Image Velocimetry (PIV) technique allowed the characterization of the air curtain with higher spatial resolution and accuracy. The result showed a zone of air recirculation at the upper part of the cabinet where warm ambient air infiltration through the door gaps was also observed, leading to an increase in the air curtain temperature. A 2D-CFD k-ε turbulence model was developed to reproduce the main flow phenomena observed by PIV so that its influence on the internal temperature distribution can be examined.Finally, a simplified heat transfer model was developed based on a zonal approach in both static and dynamic regimes which permits, respectively, the predictions of time-averaged air and product temperatures and temperature fluctuations according to the on/off cycle of the compressor regulation. The dynamic model was solved with a spectral approach, thus the influence of different parameters on the damping of the temperature fluctuations can be identified.