Ce travail de thèse a été consacré à la simulation numérique des transferts par convection forcée dans des écoulements turbulents pleinement développés. Deux configurations ont été considérées : la conduite cylindrique et l'espace annulaire. Les champs turbulents sont analysés statistiquement en moyennant dans le temps et dans l'espace les champs instantanés de vitesse et de température afin d'en déduire les profils moyens, les écart-types, les flux de chaleur turbulents ainsi que des statistiques d'ordres plus élevés telles que les coefficients de dissymétrie et d'aplatissement. Une approche plus poussée dans l'analyse des écoulements turbulents est aussi utilisée : les fonctions de densité de probabilité (pdf) et les fonctions de densité de probabilité jointes (jpdf). On s'intéressera à l'influence de la diffusivité thermique du fluide sur les transferts thermo-convectifs dans les deux géométries. On étudiera l'effet de la rotation de la paroi sur le transfert de chaleur dans la conduite cylindrique. Dans la géométrie annulaire, l'influence du rapport des flux de chaleur imposés sur les parois est analysée / This thesis was devoted to the numerical simulation of heat transfer by forced convection in fully developed turbulent flows. Two configurations were considered : the pipe and the annulus flows. The turbulent fields are statistically analyzed by averaging the instantaneous fields of velocity and temperature in order to deduce the mean profiles, root mean square, turbulent heat fluxes as well as higher order statistics such as skewness and flatness coefficients. A more thorough in the analysis of turbulent flows is also used : probability density functions (pdf) and joined probability density functions (jpdf). We look at the influence of the thermal diffusivity of the fluid on the convective heat transfer in both geometry, the influence of the ratio of heat fluxes imposed on the walls is analyzed

Le but de la présente étude est d'élargir la base de connaissances sur l'écoulement et le transfert de chaleur dans les échangeurs de chaleur à plaques. Pour cela, une étude expérimentale et des simulations numériques ont été réalisées. Le montage expérimental a servi à valider le modèle numérique et à explorer le régime de convection mixte. Le modèle numérique représente fidèlement l'ECP utilisé dans ce travail, qui contient deux canaux confinés par des plaques de type chevron. La validation a été faite pour le régime laminaire et pour le régime turbulent, dans lequel un modèle de turbulence a été choisi comme le plus adéquat, suite à une étude comparative entre plusieurs modèles. Le modèle numérique a été utilisé pour compléter les données expérimentales et pour explorer des niveaux intangibles expérimentalement. Ainsi, les champs thermique et d'écoulement dans cet ECP ont été étudiés pour deux combinaisons de fluide : eau - eau et eau - huile de moteur. On a mis en évidence la non-uniformité du champ de températures et le fait que pour ce type particulier d'ECP, les fluides suivent préférentiellement deux passages lisses situés le long des bords des canaux. La dépendance du nombre de Reynolds de ces phénomènes a été mise en évidence. Le modèle numérique a pu être modifié pour évaluer l'effet des canaux lisses périphériques sur le transfert de chaleur et les pertes de pression. Ce modèle consiste de la même géométrie, avec ces canaux périphériques obturés. Cette comparaison montre que l'absence des passages lisses provoque l'augmentation du transfert thermique et aussi, des pertes de pression. Une autre modification portant sur la comparaison des configurations verticale et en diagonale a été faite numériquement. Des différences quantitatives insignifiantes concernant les performances thermiques et hydrauliques ont été remarquées. Une investigation sur l'influence de la convection naturelle sur le transfert de chaleur et les pertes de pression a été faite. Dans ce but, le montage expérimental a été muni d'un pivot avec lequel l'ECP pouvait être inversé. On a conclu que la position de l'ECP influence significativement sur les performances thermiques et sur les pertes de pression dans le régime d'écoulement de faible nombre de Reynolds.

Ce travail de thèse a été consacré à la simulation numérique des transferts par convection forcée dans des écoulements turbulents pleinement développés. Deux configurations ont été considérées : la conduite cylindrique et l'espace annulaire. Les champs turbulents sont analysés statistiquement en moyennant dans le temps et dans l'espace les champs instantanés de vitesse et de température afin d'en déduire les profils moyens, les écart-types, les flux de chaleur turbulents ainsi que des statistiques d'ordres plus élevés telles que les coefficients de dissymétrie et d'aplatissement. Une approche plus poussée dans l'analyse des écoulements turbulents est aussi utilisée : les fonctions de densité de probabilité (pdf) et les fonctions de densité de probabilité jointes (jpdf). On s'intéressera à l'influence de la diffusivité thermique du fluide sur les transferts thermo-convectifs dans les deux géométries. On étudiera l'effet de la rotation de la paroi sur le transfert de chaleur dans la conduite cylindrique. Dans la géométrie annulaire, l'influence du rapport des flux de chaleur imposés sur les parois est analysée

La méthode SEA (Statistical Energy Analysis) permet de prévoir les échanges d'énergie entre sous-systèmes d'une structure lorsque celle-ci est soumise à une excitation. La SEA est théoriquement très simple à mettre en place et permet de faire des variations paramétriques sur un modèle de la structure en phase de conception. Cette méthode est utilisée de manière pertinente dans le domaine du bâtiment mais peine à s'implanter dans le secteur des transports où les systèmes étudiés sont beaucoup plus complexes. Plusieurs difficultés apparaissent dans ces cas de figure et empêchent une modélisation correcte de la structure. La première difficulté vient de la nécessité de découper la structure en sous-systèmes respectant les hypothèses de la SEA. Ce découplage, presque trivial dans le bâtiment, devient particulièrement complexe dans le cas d'une caisse de voiture par exemple. L'estimation des puissances injectées par les sources dans la structure est une deuxième difficulté. En effet, les énergies des sous-systèmes sont directement proportionnelles à la puissance injectée. Ainsi, une mauvaise estimation de la puissance entraîne une erreur sur les échanges d'énergies. Le but de ce travail de thèse est double. Dans la première partie, une modèle de puissance injectée dans une plaque rectangulaire soumise à une couche limite turbulente (CLT) établie. Ce modèle simple permet de souligner l'influence des différents paramètres agissant sur la puissance injectée et d'estimer avec le temps de calcul très faible la puissance injectée par bande de fréquence. Une validation expérimentale est aussi exposée. Dans la deuxième partie, une méthode de sous-structuration automatique d'une structure en sous-systèmes SEA est proposée. Elle est basée sur l'analyse et la classification de fonctions de transfert énergétiques simulées par éléments finis. Un algorithme de classification permet de diviser la base de données en N sous-systèmes (N allant de 2 à Nmax). Un indice, validé sur des assemblages simples, indique la sous-structuration optimale pour une application SEA. Des applications industrielles, démontrant l'applicabilité de la méthode à des structures complexes, seront présentées.

Résumé: Les mécanismes de perfectionnement du transfert thermique des nanofluids sont encore peu clairs. Les études précédentes au sujet des nanofluids ont essayé de résoudre certains des nombreux défis au sujet de la performance thermique et hydrodynamique des nanofluides et de leurs propriétés ; toutefois il reste beaucoup de problèmes non résolus et questions sans réponse certaine. Par conséquent, plus d'études sont nécessaires, qui peuvent être expérimentales, numériques ou théoriques. Dans la présente étude, des nanofluides sont étudiés intensivement en utilisant des approches numériques et analytiques. La partie numérique se compose de trois chapitres et couvre un éventail de problèmes de transfert thermique, incluant; laminaire et turbulente, monophasique et diphasique, aussi bien que, convection mixte et convection forcée. Plusieurs concentrations volumétriques de nanoparticules et nombres de Reynolds sont considérés. Le deuxième chapitre est consacré à la convection laminaire mixte de nanofluide d'Al[indice inférieur 2]O[indice inférieur 3]-eau à l'intérieur d'un tube horizontal. Le flux uniforme de chaleur est appliqué au mur. Deux nombres de Reynolds et trois concentrations volumétriques de nanoparticules sont utilisés, et finalement les résultats numériques thermiques et hydrodynamiques de trois différents modèles diphasiques et du modèle monophasique sont comparés aux données expérimentales. On démontre que les résultats de ces différentes approches sont extrêmement différents. Pour un régime de convection laminaire mixte, les modèles diphasiques sont en meilleur accord avec des données expérimentales. Les résultats de modèles diphasiques sont proches mais loin des résultats du modèle monophasique. Le troisième chapitre évalue la sensibilité de la formulation laminaire sur des combinaisons choisies des expressions pour la conductivité et la viscosité des nanofluids. Deux expressions pour la conductivité et trois pour la viscosité sont choisis, ce qui donne six combinaisons. Ces choix s'avèrent avoir des effets très importants sur les résultats finals. Par conséquent, chaque étude numérique devrait d'abord justifier son choix des corrélations de viscosité-conductivité. En outre, une liste des modèles les plus importants pour la conductivité et la viscosité des nanofluids est recueillie et incluse dans ce chapitre. Le quatrième chapitre évalue les résultats du modèle monophasique et trois différents modèles diphasiques pour la convection forcée turbulente de nanofluide dans un tube horizontal. Le flux uniforme de la chaleur est appliqué au mur. Le modèle turbulent "Realizable k-[epsilon]" est employé, qui est un modèle à deux équations. Deux ensembles de données expérimentales pour différents nanofluides (Al[indice inférieur 2]O[indice inférieur 3]-eau et Cu-eau) sont employés, qui couvrent un éventail des concentrations volumétriques de nanoparticules et de nombres de Reynolds. L'exactitude monophasique des résultats est confirmée avec un choix approprié de combinaisons de conductivité-viscosité. Les résultats des différents modèles diphasiques sont proches; cependant, ils sont très loin des résultats monophasique [i.e. monophasiques] et des données expérimentales. Les modèles diphasiques ne pourraient pas satisfaire les données expérimentales pour le régime convection forcée turbulente de deux nanofluides différents par deux différentes études expérimentales, alors que l'approche monophasique le fait bien. Dans la partie analytique de l'étude, de nouveaux modèles pour la conductivité thermique des nanofluides et le nombre de Nusselt de l'écoulement autour des nanoparticules sont dérivés. Ces modèles tiennent compte de l'effet du mouvement Brownien, de la résistance thermique surfacique, du groupement des particules, de la distribution de taille de ces groupements et de la micro-convection aussi bien que de la concentration des particules, de la dimension particulaire et de la température. Le groupement des particules et la distribution de leur taille sont analysés à l'aide de la théorie fractale. Le modèle proposé pour la conductivité des nanofluides est comparé aux données expérimentales de plusieurs études pour cinq nanofluides différents et différentes concentrations volumétriques de nanoparticules. Ce modèle est également comparé à deux modèles semblables. II montre une très bonne concordance avec l'expérience et une meilleure performance comparé à ces modèles choisis.||Abstract: The mechanisms of nanofluids heat transfer enhancement are still unclear. Previous studies about nanofluids have tried to solve some of many challenges about the thermal and hydrodynamic performance of nanofluids and their properties; however still there are many problems unsolved and questions without a certain answer. Hence, more studies are necessary, which can be experimental, numerical and theoretical. In the present study, nanofluids are investigated intensively using numerical and analytical approaches. The numerical part consists of three chapters and covers a wide range of heat transfer problems, including; laminar and turbulent, single-phase and two-phase as well as mixed convection and forced convection flows. Several particle volume fractions and a large number of Reynolds numbers are considered. Chapter two is dedicated to laminar mixed convection flow of Al 2 O 3 -water nanofluid inside a horizontal tube. Uniform heat flux is applied at the wall. Two Reynolds numbers and three particle volume fractions are used, and finally the thermal and hydrodynamic numerical results from three different two-phase models and the single phase model are compared with experimental data. It is shown that the predictions of these different approaches are extremely different. For a laminar mixed convection flow, two-phase models are in better agreement with a given experimental data. The two-phase models predictions are close but far from single-phase. Chapter three evaluates the sensitivity of the laminar formulation on selected combinations of models for the conductivity and viscosity of nanofluids. Two models for the conductivity and three for the viscosity are chosen, which make six combinations. These choices are found to have very important effects on the final results. Therefore, every numerical study should first justify their choice of viscosity-conductivity correlations. Also, a list of the most important models for the conductivity and viscosity of nanofluids are gathered and included in this chapter. Chapter four evaluates the predictions of single-phase and three different two-phase models for turbulent forced convection inside a horizontal tube. Uniform heat flux is applied at the wall. Realizable k-[varepsilon] turbulent model is used, which is a two-equation model. Two sets of experimental data for different nanofluids (Al 2 O3 -water and Cu-water) are used, which cover a wide range of volume fractions and Reynolds numbers. The single-phase results accuracy is confirmed with an appropriate selection of conductivity-viscosity combination. The results from different two-phase models are found to be very close; however, they were too far from the single-phase predictions and the experimental data. Two-phase models could not satisfy the experimental data for turbulent forced convection flow of two different nanofluids from different experimental studies, while single-phase approach does it well. In the analytical part of the study, new models for the thermal conductivity of nanofluids and the Nusselt number of the flow around the nanoparticles are derived. These models take into account the effect of Brownian motion, interfacial thermal resistance, particles clustering, clusters size distribution and micro-convection as well as particles concentration, particles size and temperature. The clusters size and size distribution are analyzed based on the fractal theory. The proposed model for the conductivity of nanofluids is compared with experimental data from several studies for five different nanofluids and various magnitudes of volume fractions. This model is also compared with two similar models. It shows very good agreement with experiment and better performance compared to those selected models.[symboles non conformes]

La Société le Nickel ’SLN’ est l’un des plus grands exploitants de nickel en Nouvelle Calédonie. Le Stockage des déchets de l’exploitation minière sur le site d’exploitation se fait sous forme de tas manufacturés appelés versants. Pour assurer la stabilité de ces versants, il est nécessaire d’assurer un stockage hors d’eau de ces stériles miniers. Pour se faire, une mise en place des ouvrages en enrochement enterrés se fait au fur et à mesure de l’avancement des travaux d’installation des versants. Le rôle de ces ouvrages enterrés est de maintenir l’écoulement d’eau à travers le milieu poreux. Le premier objectif de la thèse est donc de mieux comprendre les interactions hydrodynamiques dans les milieux très perméables dans le cas d’écoulements pleinement turbulents. La première partie de ce travail concerne une étude expérimentale menée à l’Institut National Agronomique de Tunisie (INAT) mais aussi à l’Institut de mécanique de fluide de Toulouse (IMFT) dont le but est d’étudier à la fois l’évolution de la ligne d’eau dans les milieux macroporeux et la variation de la perméabilité de ces derniers en fonction du nombre de Reynolds. Des simulations avec des modèles de perte de charge les plus utilisés dans la littérature ont été effectuées, pour tester leur capacité à reproduire la ligne d’eau dans les conditions de nos expériences. Les deux relations de Forchheimer (1901) et d’Ergun (1952) donnent des résultats proches de nos expériences avec des valeurs de coefficients adaptés pour notre cas, alors que la relation de Barree et Conway (2004) ne semble pas reproduire les résultats attendus. L’eau qui circule dans les ouvrages enterrés est une eau chargée d’où l’intérêt d’étudier le transport sédimentaire dans ces milieux macroporeux. Pour se faire, une étude expérimentale a été réalisée sur le transport et le dépôt des particules dans des milieux synthétiques constitués de particules de diamètre moyen de 15 cm

L’objet de cette thèse est de modéliser et de simuler des écoulements turbulents diphasiques incompressibles à phases non miscibles. La modélisation et la simulation de ce type d’écoulements sont traitées dans le cadre des méthodes de Simulation des Grandes Echelles (SGE) ou Large Eddy Simulation (LES) en anglais qui consistent à calculer directement les plus grandes structures de l’écoulement et à modéliser les plus petites. Ces méthodes adaptées aux écoulements turbulents monophasiques sont étendues au cadre des écoulements turbulents diphasiques. Pour cela, elles sont couplées avec une méthode eulérienne de type ’ Volume Of Fluid’ (VOF) spécifique au caractère diphasique de l’écoulement. La pertinence du couplage entre les modélisations SGE et VOF est testée sur la configuration industrielle proposée par le CEA-CESTA: l’impact d’un jet rond turbulent sur une surface libre eau/air dans une cavité. Des mesures expérimentales de vitesse (Particle Image Velocimetry PIV) réalisées au CEA-CESTA sont disponibles pour valider les résultats numériques issus des simulations. / The scope of this dissertation is to model and simulate non-miscible two-phase turbulent and incompressible flows. The modeling and the simulation of this kind of flows are carried out in the framework of the Large Eddy Simulation (LES) which consists in calculating directly the largest structures of the flow and in modeling the finest ones. These numerical methods, applied usually to the simulation of single-phase turbulent flows, are extended to the simulation of two-phase turbulent flows in this work. Hence, the LES methods are coupled with an Eulerian ’Volume of Fluid’ (VOF) approach which is particularly adapted to interfacial flows. The relevance of this numerical coupling bewtween LES and VOF methods is validated in the following industrial configuration of the CEA-CESTA: the impact of a turbulent round jet on a free water/air surface in a cavity. Some experimental velocity measurements (Particle Image Velocimetry PIV), carried out at the CEA-CESTA, are available to validate the numerical results.

Cette thèse présente les résultats de deux études : la première concerne la convection naturelle turbulente dans une cavité rectangulaire chauffée uniformément par le bas et remplie d’un nanofluide et la seconde concerne l’investigation du transfert de chaleur conjugué dans un échangeur de chaleur à tubes ailetés.L’enceinte de la première étude a un faible rapport d’aspect. Ses parois gauche, droite et supérieure sont maintenues à une température relativement basse. Le fluide de travail est un nanofluide constitué d’eau et de nanoparticules, soient d’alumine (Al2O3), ou de cuivre (Cu) ou d’oxyde de cuivre (CuO). L’influence des paramètres tels que le nombre de Rayleigh (basé sur la hauteur H de la cavité et la densité de flux de chaleur), le type de nanofluide et la fraction volumique des nanoparticules sur la performance de refroidissement est présentée. Les équations de Navier-Stokes et les équations de conservation de la masse et de l'énergie sont résolues pour une géométrie bidimensionnelle par la méthode numérique des volumes finis. L'algorithme SIMPLE est utilisé pour le couplage pression-vitesse. La discrétisation des termes convectifs est faite avec le schéma QUICK. Le modèle de turbulence k-epsilon standard est utilisé. Le maillage du domaine simulé est généré par le code Gambit. Les résultats montrent que pour toutes les valeurs de Ra, le nombre de Nusselt moyen augmente d’une façon linéaire et monotone avec l’augmentation de la concentration des nanoparticules dans le fluide de base. Le flux de chaleur moyen prend des valeurs qui diminuent en fonction de l’ordre suivant : Cu, CuO et Al2O3.La deuxième étude est une investigation numérique des caractéristiques dynamique et thermique d'un échangeur de chaleur. Les calculs supposent un transfert de chaleur et un écoulement en régime permanent. Le nombre de Nusselt et le coefficient de frottement qui caractérisent l'échangeur de chaleur sont déterminés pour différentes valeurs du nombre de Reynolds. L’équation de conservation de l'énergie dans le fluide et l’équation de conduction de la chaleur dans le solide en trois dimensions ont été résolues avec les équations de la conservation de la masse et de la quantité de mouvement afin de déterminer ces caractéristiques. Les deux régimes d'écoulement laminaire et turbulent sont considérés. L'effet de la modélisation de la turbulence a été étudié en utilisant trois modèles différents (Spalart-Allmaras modèle de turbulence à une équation, le modèle k-epsilon ; standard et le modèle RSM). La validation du modèle a été effectuée en comparant les facteurs de frottement, f, et le facteur, j, de Colburn avec les données expérimentales trouvées dans la littérature. Les résultats tracés ont montré un bon accord qualitatif entre les résultats numériques et les données expérimentales. Les résultats montrent également que le plus simple des trois modèles de turbulence testés (à savoir. Spalart-Allmaras) donne les valeurs les plus proches des données expérimentales. / This thesis presents the results of two studies: the first concerns natural turbulent convection in a rectangular cavity heated from the bottom wall and filled with a nanofluid and the second relates to the investigation of conjugate heat transfer in a fin-and-tube heat exchanger.The cavity of the first study is tall and has a heat source embedded on its bottom wall, while its left, right and top walls are maintained at a relatively low temperature. The working fluid is a water based nanofluid having three nanoparticle types: alumina, copper and copper oxide. The influence of pertinent parameters such as the Rayleigh number, the type of nanofluid and solid volume fraction of nanoparticles on the cooling performance is studied. Steady forms of twodimensional Reynolds-Averaged-Navier-Stokes equations and conservation equations of mass and energy, coupled with the Boussinesq approximation, are solved by the volume control based on the discretisation method employing the SIMPLE algorithm for pressure-velocity coupling. Turbulence is modeled by using the standard k-epsilon model. The Rayleigh number, Ra, is varied from 2.49xE09 to 2.49xE11. The volume fractions of nanoparticles where varied too. Stream lines, isotherms, velocity profiles and temperature profiles are presented for various combinations of Ra, the type of nanofluid and solid volume fraction of nanoparticles. The results are reported in the form of average Nusselt number on the heated wall. It is shown that for all values of Ra, the average heat transfer rate from the heat source increases almost linearly and monotonically as the solid volume fraction increases. Finally the average heat transfer rate takes on values that decrease according to the ordering Cu, CuO and Al2O3.In the second study We determined the heat transfer and friction characteristics of a realistic fin-and-tube heat exchanger. The computations assume steady-state heat transfer and fluid flow. Nusselt number and friction factor characteristics of the heat exchanger are presented for various values of Reynolds numbers. The energy conservation and the heat conduction equations in 3 dimensions have been solved in the fluid and the solid respectivelyalong with the mass and momentum conservation equations in order to determine these characteristics. Both laminar and turbulent flow regimes are considered. The effect of turbulence modeling was investigated using three different models (the one equation Spalart-Allmaras turbulence model, the standard k-epsilon; model and the RSM model). The computations allowed the determination of the dynamic and thermal fields. Model validation was carried out by comparing the calculated friction factor f and Colburn j-factor to experimental results found in the literature. The plotted results showed a qualitatively good agreement between numerical results and experimental data. The results obtained also showed that the simplest of the three turbulence models tested(i.e. Spalart-Allmaras) gives the closest values to the experimental data.

Dans cette thèse, on présente une étude numérique de la convection forcée turbulente dans un dissipateur thermique muni d'une chicane transversale dans le by-pass. Le premier modèle est composé d’ailettes planes et le second consiste à ajouter des ailettes broches entre les ailettes planes. Les équations gouvernantes basées sur le modèle de turbulences k- SSt sont discrétisées et résolues par la méthode des volumes finis et l'algorithme SIMPLE. Les résultats dynamiques sont présentés en en termes de champs de vitesse, des profils de vitesse axiales dans des sections choisies ainsi que la perte de charge. L'étude thermique est présentée en terme de champs de température et de distribution du nombre de Nusselt. Un rapport entre les performances thermique et dynamique est présenté pour évaluer les différents dissipateurs thermiques. / In this thesis, we present an numerical study of turbulent forced convection in a heat sink provided with a transverse baffle in the bypass. The first model is composed of plates fins and the second consists of adding pin fins between the plates fins. The governing equations, based on the k- SSt turbulence model, are disscredized and solved by the finite volume method and the SIMPLE algorithm. Dynamic results are presented in terms of velocity fields, profiles of the axial velocities in selected sections and pressure drop. The thermal study is presented in terms of temperature fields and the distribution of Nusselt number. A ratio between the thermal and dynamic performances is presented to evaluate the different heat sinks.

Cette thèse présente une étude numérique d’un séchoir solaire indirect à convection naturelle destiné à sécher les produits agricoles (les figues). La première partie analyse un séchoir solaire contenant une chambre de séchage couplée à un absorbeur inversé et une cheminée solaire. Les simulations ont été faites afin de déterminer les champs dynamique et thermique sous l’influence de la variation de la configuration de la cheminée solaire et la taille de l’ouverture d’admission. Les équations de conservations basées sur le modèle de turbulence k-ε standard sont résolues par la méthode des volumes finis à l’aide du code commercial ANSYS-Fluent. La deuxième partie analyse l’effet d’ajouter un stockage thermique sous forme d’un lit en gravier dans la chambre de séchage. Le lit est modélisé comme un milieu poreux. En plus, cette étude propose l’utilisation d’une deuxième entrée d'air dans la chambre de séchage afin d’assurer une distribution thermique plus homogène au niveau des claies et de garantir un séchage plus uniforme. / This thesis presents a numerical study of an indirect natural convection solar dryer for drying agricultural products (Figs). The first part analyzes a solar dryer containing a drying chamber coupled to a reversed absorber and a solar chimney. Simulations were made to determine the dynamic and thermal fields under the influence of the variation of the solar chimney configuration and the size of the inlet opening. The governing equations based on the standard k-ε turbulence model are solved by the finite volume method using the ANSYS-Fluent commercial code. The second part analyzes the effect of adding a thermal storage in the form of a gravel bed in the drying chamber. The bed is modeled as a porous medium. Furthermore, this study proposes the use of a second air inlet in the drying chamber in order to ensure a more homogeneous thermal distribution at the level of the racks and to guarantee a more uniform drying.