Dans les écoulements de convection naturelle en canal, de nombreuses structures cohérentes se développent. Leur rôle dans les transferts de masse et de chaleurs reste aujourd’hui mal connu. Des études numériques basées sur des Simulations aux Grandes Echelles d’un canal vertical dont un des côtés est chauffé uniformément à flux constant ont donc été entreprises. Différentes étapes de la transition laminaire-turbulente sont identifiées et deux bandes de fréquences sont observées en proche paroi chauffée. Une méthode, basée sur la Proper Orthogonal Decomposition (POD), a aussi été utilisée et a permis une très nette séparation spatio-temporelle des modes les plus énergétiques de l’écoulement en fonction des fréquences précédemment identifiées. Il devient donc possible d’étudier la contribution de ces modes à l’écoulement. Une analyse POD des champs de vitesses expérimentaux permet l’identification de similarités avec les structures observées numériquement.Les variations du gradient de température extérieur au canal, ont été identifiées comme l’une des sources de ces différences et ont fait l’objet d’études expérimentales et numériques. Il est montré qu’un augmentation du gradient positif de température diminue le débit mais aussi déplace la transition plus bas dans le canal. Les études numériques permettent également d’étudier le cas de gradients de température faibles ou négatifs, gradients qui sont difficile à obtenir expérimentalement. Un modèle théorique de l’influence de la stratification extérieure sur le débit massique induit est aussi développé. L’accord entre les résultats prédits et les résultats expérimentaux et numériques est excellent. Ce dernier point souligne que la stratification thermique est un facteur clé dans ce type d’écoulement et que son influence est précisément quantifiée dans ce travail / In natural convective flows, complex coherent structures develop whose the role in heat and mass transfer are not well understood. A numerical study, based on Large-Eddy-Simulations of a vertical channel with one side uniformly was therefore carried out. Different stages of transitional flow development were identified numerically with two characteristic frequency bands being observed in the flow, near the heated wall. Methods derived from the Proper Orthogonal Decomposition (POD) was also used and allows the most energetic modes to be separated accordingly to two characteristic frequency bands found numerically. As result, the contribution of the two families of modes to the near wall turbulent heat transfer and velocity-temperature correlation has been evaluated. POD was also performed on experimental measurements showing similarities with the numerically observed structures.In this work, variations of the external thermal stratification have been identified as one possible source of these differences and its influence was therefore investigated experimentally and numerically. It is shown that the increase in the positive gradient of the external stratification not only decreases the mass flow rate but also displaces the transition height to a lower location in the channel. Numerical simulations also allow the study of cases of weak and negative thermal stratifications which are difficult to achieve in laboratories. A theoretical model of the influence of the external thermal stratification on the mass flow rate was also developed. There is an excellent agreement between the theoretical predictions and the experimentally and numerically obtained mass flow rates. This clearly highlights that external temperature distributions are key driving factors and their influence is accurately quantified in this work

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.

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.

Une étude expérimentale sur la faisabilité du procédé de dessalement par congélation sur les parois d'un tube cylindrique est présentée. Le procédé se décompose en deux étapes essentielles : la cristallisation suivie par le ressuage. Les cristallisations sont conduites à partir d'une solution stagnante ou agitée par un bullage d'air. Les paramètres opératoires clés sont la rampe de refroidissement et la salinité de l'eau traitée. Leur effet sur la salinité de la glace produite a été quantifié, dans les deux modes de fonctionnement, statique et agité. Le dispositif expérimental permet également de contrôler le gradient thermique à travers la couche de solution pour les cristallisations effectuées en statique. L'absence de gradient thermique et donc de courant de convection mène à des couches de glace très contaminées en sels. Les résultats combinant la cristallisation et le ressuage montrent que dans tous les cas, on peut obtenir de l'eau potable mais avec un rendement plus élevé et une durée réduite si la cristallisation a lieu en mode dynamique. Les conditions opératoires optimales sont proposées. La modélisation des transferts de matière et de chaleur à travers les couches de glace et de solution en régime diffusionnel et convectif a permis d'interpréter l'influence des paramètres opératoires des cristallisations opérées sans agitation. Enfin, l'évaluation économique du procédé complet de dessalement, basées sur des points de fonctionnement expérimentaux, indique que la consommation énergétique d'une petite installation pourrait être très faible en utilisant une machine frigorifique idéale fonctionnant entre l'unité de ressuage et l'unité de congélation.

Ce travail concerne la convection naturelle au sein d'un système fluide-poreux en couche horizontale. On présente l'analyse de stabilité linéaire des modèles à un et deux domaines, avec diffusion visqueuse dans le milieu poreux. Nos résultats sont comparés avec ceux du modèle à deux domaines utilisant la formulation de Darcy. Un bon accord est observé entre les résultats des modèles à deux domaines, ce qui indique que le terme de Brinkman joue un rôle secondaire dans la stabilité. On montre que le modèle à un domaine peut conduire à des résultats sensiblement différents lorsque la transition entre fluide et le milieu poreux est décrite par une discontinuité des propriétés. Il faut alors modifier la formulation en effectuant la différentiation au sens des distributions. Ainsi, le modèle à un domaine conduit aux mêmes seuils de stabilité que les formulations à deux domaines. L'influence des paramètres caractéristiques sur la stabilité des systèmes thermique et thermosolutal est discutée.

Ce travail se situé dans la perspective d'un contrôle de la structure de solidification des alliages métalliques sous l'effet de la convection naturelle et forcée afin d'améliorer à terme la maîtrise des microstructures de solidification grâce à un brassage électromagnétique efficace permettant d'avoir une homogénéisation du bain liquide qui par la suite peut améliorer la microstructure finale de l'alliage. La possibilité retenue dans ce travail est de réaliser ce brassage sans contact avec la solution liquide (alliage sous fusion) et sans pollution par d'autres éléments en utilisant un brassage par la force de Lorentz. L'objet de la thèse comporte une étude théorique à la fois expérimentale basée sur une installation expérimentale particulièrement documentée et instrumentée, développée au laboratoire SIMAP/EPM à Grenoble, nommée AFRODITE. Ce dispositif expérimental permet de fournir des données de quantitatives et qualitatives sur le processus de solidification des alliages métalliques. Ces données sont nécessaires à la contribution aux études menées sur la solidification des alliages métallique et enrichir la base des donnée des modèles numériques développés pour prédire les défauts liés à la solidification. L'alliage choisi dans notre travail est l'étain-plomb, vue sa basse température de fusion. Les expériences envisagées visent à étudier l'effet de deux modes de configuration dynamique sur la solidification de l'alliage Sn-Pb: la convection thermosolutale avec la variation de deux paramètres essentiels (la vitesse de refroidissement et la différence de température expérimentale) et la convection forcée par l'utilisation de plusieurs modes de brassage électromagnétique. Cette étude s'intéresse en particulier à la caractérisation des macrostructures et les défauts liés à la macroségrégation. L'originalité de l'étude vise à mesurer in situ les températures instantanées. Ceci nous a permis d'évaluer l'évolution du transfert thermique due à la convection naturelle/forcée, ainsi que leurs influence sur le processus de la solidification sous différents aspects. L'analyse post-mortem de l'alliage métallique, fournit la structure de solidification et la distribution des ségrégations à différentes échelles (mésoscopique et macroscopique).

Nous présentons dans cette étude le développement de la convection à partir de différentes perturbations de l'état conductif d'une couche fluide confinée dans une cavité cylindrique, chauffée par le bas et avec une surface supérieure libre. La discrétisation spatiale du domaine repose sur la méthode des éléments spectraux et les itérations temporelles sont assurées par une méthode splitting.Au déclenchement de la convection, les structures convectives correspondent à des modes de Fourier, et les seuils critiques dépendent du rapport de forme de la cavité, et des nombres de Biotet de Marangoni qui caractérisent la surface libre. Les transitions d'écoulements au-delà du seuil primaire sont caractérisées quantitativement en fonction du nombre de Rayleigh pour différentes valeurs du nombre de Biot et Ma = 0. Les résultats présentés sont obtenus en résolvant l'ensemble des équations non-linéaires de conservation à travers une méthode de continuation. Lorsque la convection se déclenche sous la forme d'un mode axisymétrique m = 0, l'évolution non-linéaire montre la coexistence de différentes structures convectives, des structures axisymétriques avec écoulement montant ou descendant au centre de la cavité et des structures correspondant à des combinaisons de modes qui apparaissent sur des branches secondaires sous-critiques.L'action d'un champ magnétique constant est ensuite étudiée pour des fluides conducteurs dans une même configuration comprenant une surface supérieure libre. Nous montrons l'effet stabilisateur du champ magnétique sur les seuils primaires ainsi que son action sélective sur les différents modes de convection. Nous analysons l'évolution des structures convectives au delà de ces seuils et montrons comment le champ magnétique modifie les transitions entre ces structures.En soumettant le bain fondu à un champ magnétique tournant, le mouvement de rotation du fluide se superpose aux mouvements de convection thermique et on observe une diminution des fluctuations de température et un retard du déclenchement de l'instabilité de Rayleigh-Bénard(lorsque les deux parois haut/bas du bain sont rigides). La rotation influe sur ce déclenchement qui de stationnaire devient oscillatoire, à l'exception du mode m = 0 de Fourier, pour qui la transition reste stationnaire jusqu'à une certaine valeur critique du nombre de Taylor magnétique.La dynamique de l'écoulement axisymétrique de part et d'autre de cette valeur critique sera étudiée en détail. / The growth of thermal convection out of different perturbations of the conductive base state is investigated using a spectral element time-stepping code. The fluid is subject to a vertical heat transfer in a cylindrical cavity with an upper free surface corresponding to the so-called Rayleigh-Bénard-Marangoni situation and the heat exchange through the free surface is evaluated via the Biot number. The results of the stability diagrams show that the evolution of the primary thresholds are largely influenced by the Biot number, the Marangoni number, and the aspect ratio of the cavity. Flow transitions are elucidated in quantitative detail as a function of the Rayleigh number for different Biot numbers in the tension free limit Ma = 0. The results shown are obtained by solving the full nonlinear field equations numerically among a continuation method. When an axisymmetric m = 0 Fourier mode is obtained at onset, the non-linear evolution shows the coexistence of different convective structures, the axisymmetric structures with up-ow or down-ow at the center and mixed-mode structures which appear on secondary subcritical branches. The action of a constant magnetic field is then considered for melts in the same type of configuration with an upper free surface. We show the global stabilizing effect of the magnetic field on the primary bifurcation thresholds and the selective effect on the different instability modes. We also analyze the evolution of the convective structures above the thresholds and show how the magnetic field modifies the transitions between these structures. When applying a magnetic body forcing in the azimuthal direction (RMF), one can damp the unavoidable thermal fluctuations inside the melt and delay the transition to the Rayleigh-Bénard instability (for rigid-rigid circular plates at top and bottom). The rotation effect also changes the transitions from steady to oscillatory, except for the m = 0 Fourier mode where the transitionis first steady until a critical Taylor number and then becomes oscillatory. The dynamics of the transitions to the axisymmetric flow, below and above this value of critical magnetic Taylor number, is particularly interesting and will be described.

La thèse porte sur la caractérisation des transferts thermiques dans les milieux fluide-particules, notamment en proposant un modèle décrivant le changement de phase au sein de particules sous écoulement fluide. Les transferts thermiques sont modélisés en prenant en compte l'influence de la résistance aux transferts externes (échanges avec le fluide) et internes (conduction à travers la particule et sa paroi, convection naturelle dans la phase liquide de la particule, changement de phase) à la particule. Les échanges externes avec le fluide sont pris en compte à l’aide de corrélations liant un nombre de Nusselt externe aux nombres de Reynolds et de Prandtl. La conduction interne est décrite à l’aide de solutions analytiques. L’effet de la convection naturelle a été étudié sur une particule isolée soumise à un gradient de température sous différents nombres de Rayleigh et de Prandtl permettant son déclenchement. Les résultats obtenus ont permis d’établir une corrélation reliant un nombre de Nusselt interne aux nombres de Prandtl et Rayleigh de la particule. Cette corrélation permet de recalculer l’évolution temporelle de la température moyenne de la particule en prenant en compte l’effet de la convection naturelle. Le changement de phase est décrit grâce à un modèle local basé sur l’approche « Phase Field » moyenné sur l’ensemble de la particule et validé par comparaison avec des résultats numériques et expérimentaux issus de la littérature. Enfin, le modèle complet et l’influence des phénomènes pris en compte sont testés sur un lit fixe de particules à l’échelle moyennée (Discrete Element Method–Computional Fluids Dynamics). La conduction et la convection interne donnent une quantité totale d’énergie stockée relativement similaire dans le lit à nombre de Biot égal, mais dont la distribution est différente. Le changement de phase tend à grandement densifier le stockage. Une augmentation du nombre de Biot tend à augmenter la quantité d’énergie stockée. Enfin, il est montré que les transferts sont dépendants de la distribution de porosité. / This work intends to characterize heat transfer in fluid-particle flows, specifically when phase change occurs inside the particles. The proposed model takes into account the external heat resistance (heat transfer at the particle-fluid interface) and the internal heat resistance (conduction inside and at the wall of the particle, natural convection in the liquid phase of the particle, phase change). External transfer with the surrounding fluid is described by correlations linking an external Nusselt number to Reynolds and Prandtl numbers related to the surrounding fluid. Internal conduction is calculated thanks to analytical solutions. The influence of natural convection was studied on an isolated sphere for different combinations of Rayleigh and Prandtl numbers. A correlation between an internal Nusselt number, and particle Rayleigh and Prandtl numbers was established using these simulations. This correlation allows calculating the transient evolution of the average temperature of the particle when natural convection occurs. Phase change is taken into account by a Phase Field model averaged over the particle and validated by comparison with experimental and numerical studies from the literature. Finally, the whole model and the effects of the different phenomena it describes are tested on a fixed bed of particles at mesoscopic scale using a Discrete Element Method–Computional Fluids Dynamics (DEM-CFD) model. Internal conduction and natural convection gives similar quantities of total energy stored for the same Biot number, however heat transfer distribution is modified. Phase change greatly reduces the volume of storage. Increasing the Biot number leads to a greater amount of energy stored. Finally, heat transfer greatly depends on porosity distribution.

Les infrastructures de transport jouent un rôle majeur dans le développement socio-économique des régions nordiques. La construction de remblais combinés aux changements climatiques engendrent des impacts négatifs sur le pergélisol sous-jacent, causant la dégradation des infrastructures. Des techniques de mitigation ont été proposées et testées pour limiter la dégradation du pergélisol. Toutefois, il y a peu d'informations disponibles sur les procédures de conception ou sur des lignes directrices. Le but de cette recherche est de développer des outils d'ingénierie améliorés pour les techniques de stabilisation convectives, se concentrant sur le remblai à convection d'air (ACE) et le drain thermique. Plus spécifiquement, l'approche du bilan thermique est proposée pour déterminer la condition thermique de remblais conventionnels et pour permettre la sélection de la technique de mitigation la plus appropriée pour extraire la chaleur en excès transmise au sol, si le système est actuellement instable, ou pour donner un facteur de sécurité considérant une instabilité future estimée. Quatre modèles thermiques ont été développés et calibrés à partir de données mesurées sur des sites expérimentaux. Un abaque de bilan thermique pour les remblais conventionnels ainsi que plusieurs abaques de capacité d'extraction de chaleur pour l'ACE et le drain thermique ont été développés et validés utilisant, respectivement, des simulations numériques et les données de sites d'essais du Yukon et du Nord du Québec, Canada. Ces abaques permettent aux concepteurs et ingénieurs de concevoir des remblais à convection d'air et des drains thermiques optimisés pour limiter, ou même éviter, le dégel du pergélisol à des sites spécifiques. / Transportation infrastructure plays a vital role in the social and economic development of northern regions. The construction of embankments and climate change can lead to negative impacts on the underlying permafrost, causing degradation of the infrastructure. Mitigation techniques have been proposed and tested to limit permafrost degradation. However, there is limited information on the design procedures or guidelines. The purpose of this research is to develop improved engineering tools for convective stabilization techniques, focussing on air convection embankment (ACE) and heat drain. More specifically, the heat balance approach is proposed to determine the thermal condition of conventional embankments and to allow the selection of suitable mitigation techniques to extract the amount of extra heat flowing into the foundation, if the system is currently unstable, or give a safety factor considering estimated future instability. Four thermal models have been built and calibrated with field data from experimental sites. One heat balance chart for conventional embankments and, several heat extraction capacity charts for ACE and heat drain have been developed using numerical simulations and have been validated using data from test sites in Yukon and Northern Quebec, Canada. These charts allow designers and engineers to design optimized ACE and heat drain to limit or even avoid the thawing of permafrost at specific sites.

Une étude expérimentale sur la faisabilité du procédé de dessalement par congélation sur les parois d'un tube cylindrique est présentée. Le procédé se décompose en deux étapes essentielles : la cristallisation suivie par le ressuage. Les cristallisations sont conduites à partir d'une solution stagnante ou agitée par un bullage d'air. Les paramètres opératoires clés sont la rampe de refroidissement et la salinité de l’eau traitée. Leur effet sur la salinité de la glace produite a été quantifié, dans les deux modes de fonctionnement, statique et agité. Le dispositif expérimental permet également de contrôler le gradient thermique à travers la couche de solution pour les cristallisations effectuées en statique. L'absence de gradient thermique et donc de courant de convection mène à des couches de glace très contaminées en sels. Les résultats combinant la cristallisation et le ressuage montrent que dans tous les cas, on peut obtenir de l'eau potable mais avec un rendement plus élevé et une durée réduite si la cristallisation a lieu en mode dynamique. Les conditions opératoires optimales sont proposées. La modélisation des transferts de matière et de chaleur à travers les couches de glace et de solution en régime diffusionnel et convectif a permis d’interpréter l’influence des paramètres opératoires des cristallisations opérées sans agitation. Enfin, l’évaluation économique du procédé complet de dessalement, basées sur des points de fonctionnement expérimentaux, indique que la consommation énergétique d’une petite installation pourrait être très faible en utilisant une machine frigorifique idéale fonctionnant entre l’unité de ressuage et l’unité de congélation. / An experimental study on the feasibility of sea water desalination by indirect freezing is presented. The whole process of desalination involves a freezing step, followed by a purification of the ice layer by sweating. Crystallization is led from quiescent solution or from agitated solution by air bubbling. The key operating parameters are the cooling rate, and solution salinity. Their effect on ice purity has been quantified in the static and agitated modes. The experimental setup enables as well the control of temperature gradient through the solution when the freezing step is conducted in the static mode. In the absence of temperature gradient and consequently the absence of convection currents, the ice layers formed are very contaminated in salts. The results combining crystallization and sweating show that in all cases, we can obtain drinking water, but with high yield and reduced process time when the crystallization is led in agitated mode. The optimal operating conditions are presented. The modelling of mass and heat transfers in the two phases in the diffusional and convective regimes has enabled the interpretation of the influence of operating parameters of crystallizations operated without agitation. At last, the economic evaluation of the whole desalination process, based on experimental operating points, indicates that energetic consumption of a small scale unit maybe low when using an ideal refrigerating machine which works between the sweating and the crystallization unities.