# Etude numérique de la convection naturelle au sein d'une cavité carrée avec chauffage et refroidissement variables

Slimane Manar¹, Hicham Rouijaa²*, El Alami Semma¹ et Mustapha El Alami³

¹ Laboratoire de Mécanique, Faculté des Sciences et Techniques
Université Hassan Ier, B.P. 577, Route de Casablanca, Settat, Maroc

² Université Cadi Ayyad, Faculté Poly-disciplinaire de Safi
B.P. 4162, Route Sidi Bouzid, Safi, Maroc

³ Groupe de Thermique, Laboratoire de Physique des Matériaux, Micro Electronique, Automatique et Thermique, (LPMMAT), Département de Physique, Faculté des Sciences Université Hassan II Aïn Chock, B.P. 5366 Maârif, Casablanca, Maroc

(Reçu le 9 Décembre 2013 – accepté le 30 Mai 2014)

**Résumé** - Le travail présenté dans ce papier est une suite d'une série de travaux menés pour l'étude numérique de la convection thermosolutale en phase fluide d'une cavité de croissance cristalline de type Bridgman verticale. Le but est d'étudier numériquement la convection naturelle au sein d'une cavité carrée avec chauffage et refroidissement variables. La cavité contient un fluide binaire caractérisé par un nombre de Prandtl Pr = 0.01 (cristal fluide) Les parois horizontales inférieure et supérieure de la cavité étudiée (Fig. 1) sont portées respectivement aux températures chaude ($T_c$) et froide ($T_f$). Ces températures sont uniformes et leurs amplitudes, variant sinusoidalement dans le temps, oscillent autour d'une valeur moyenne fixe. Les parois verticales sont portées à la même température que celle du bas, sauf une zone considérée adiabatique du côté de la paroi froide. Les équations de Navier-Stokes et l'équation de conservation de l'énergie sont résolues par une méthode de volumes finis de 2ème ordre en temps et en espace. Les résultats, obtenus sont donnés pour $10^3 < Ra < 10^5$, Pr = 0.01, avec les amplitudes des températures chaude et froide ayant été choisies égales à 1. L'accent est mis sur l'effet des fréquences de modulation des températures chaudes et froides sur l'écoulement et le transfert thermique. Ces résultats montrent l'existence d'une fréquence critique de la température chaude minimisant le flux d'écoulement et le transfert thermique. Cette fréquence est indépendante de la fréquence de la température froide.

**Abstract** – The work presented in this paper is a continuation of a series of work for the numerical study of convection thermosolutal phase fluid from a cavity vertical Bridgman crystal growth category. The aim is to study numerically the natural convection in a square cavity with heating and cooling variables. The cavity contains a fluid binary characterized by a number of Prandtl Pr = 0.01 (crystal fluid). The lower and upper horizontal walls of the cavity studied (Fig. 1) are brought respectively to hot ($T_c$) and cold ($T_f$) temperatures. These temperatures are uniform and their amplitudes varying sinusoidally in time, oscillate around a fixed mean value. The vertical walls are brought to the same temperature as that of the bottom, except a given area adiabatic on the side of the wall cold. The Navier-Stokes equation and energy conservation are solved by a finite volume method to the second order in time and space. The results obtained are given for $10^3 < Ra < 10^5$, Pr = 0.01, with the amplitudes of the hot and cold temperatures have been chosen equal to 1. The focus is on the effect of modulation frequencies of hot and cold temperatures on the flow and heat transfer. These results show the existence of a critical frequency of hot temperature minimizing the flow stream and heat transfer. This frequency is independent of the frequency of the cold temperature.

**Mots clés:** Convection naturelle - Simulation numérique - Fréquence de modulation - Croissance cristalline.

<footnote>* hrouijaa@hotmail.com</footnote>
# 1. INTRODUCTION

Industriellement, la forte demande en matériaux utilisés dans de nombreux dispositifs en électronique, optique et optoélectronique, a permis l'émergence de nouveaux cristaux adaptés et optimisés vis-à-vis d'une utilisation spécifique. La maîtrise des outils et des conditions de croissance cristalline s'avère primordiale pour l'obtention de cristaux de grande pureté et de haute qualité cristalline répondant aux objectifs recherchés. Le couplage entre la microstructure du cristal et l'écoulement dû à la convection thermosolutale dans la phase liquide au cours de la cristallisation a été démontré dans plusieurs travaux [1-3]. Dans la littérature, l'étude de la convection thermosolutale dans les cavités rectangulaires a été effectuée pour plusieurs configurations géométriques (horizontale, verticale, inclinée, aplatie...). Différents scénarios de conditions aux limites thermiques et solutales (coopérants, opposées ou croisées) ont été investiguées.

Dans cette étude, nous nous sommes limités à l'étude du phénomène au niveau de la zone basse de du procédé Bridgman et à l'influence de la modulation des températures chaude et froide sur l'écoulement et le transfert thermique dans cette zone qui ressemble à une cavité carrée. Ainsi nous réduisons notre investigation à l'étude du phénomène de l'écoulement et le transfert thermique dans cette cavité.

Le problème de la convection naturelle de fluides confinés dans des cavités a fait l'objet de nombreuses études. Une excellente revue des travaux effectués dans le passé dans ce domaine a été présentée par Ostrach [4] et Bejan *et al.*, [5]. Ces travaux indiquent clairement que la plupart des études disponibles dans ce domaine portent sur des écoulements convectifs en régime permanent, ce qui est le cas de la plupart des applications pratiques.

Les premières études portent sur le comportement instationnaire de la convection naturelle concernant le cas de cavités soumises à des conditions thermiques permanentes. Pour cette situation, l'obtention de régime convectif transitoire, observé dans certaines conditions, est due à l'existence d'instabilités dynamiques au voisinage des parois [6-8], ainsi qu'à la formation d'instabilités thermiques [9].

De tels écoulements ont également été observés dans le cas de distributions spéciales de la température sur les parois de la cavité [10, 11] ou pour des cavités de grandes extensions [12]. En ce qui concerne la convection naturelle transitoire induite par l'imposition de conditions thermiques sur la cavité variant dans le temps, il existe relativement peu de travaux disponibles dans la littérature.

Parmi les premières contributions, on peut citer celle de Patterson *et al.*, [11] qui ont étudié le cas d'une couche fluide initialement au repos dans une enceinte dont les parois verticales sont soumises à une variation brusque de la température à un instant donné. L'effet du refroidissement soudain d'une paroi verticale d'une cellule remplie d'un fluide initialement à température uniforme a été considéré par Nicolette *et al.*, [7]. Les résultats numériques sont en bon accord avec les mesures expérimentales obtenues par les mêmes auteurs.

Dans une étude numérique et analytique, Yang *et al.*, [13] ont analysé la réponse de l'écoulement dans une cavité verticale de grande extension dont l'une de ses parois verticales est soumise à une température qui varie périodiquement dans le temps. Pour des fréquences faibles et élevées, ils ont observé que les caractéristiques de l'écoulement s'accordent bien avec celle de la convection naturelle le long d'une plaque plane verticale soumise aux mêmes conditions thermiques.
Kazmierczak et al., [14] ont étudié la convection naturelle transitoire dans une cavité carrée due à une variation sinusoïdale de la température d'une paroi verticale chaude. La paroi opposée (froide) est maintenue à une température constante constamment inférieure à celle de la paroi verticale. Ils ont analysé les effets de la période et l'amplitude de la température oscillante à travers la paroi froide de la cavité. Toutes les solutions transitoires obtenues sont périodiques dans le temps, et montrent que, malgré la dépendance des conditions aux limites vis-à-vis du temps et les variations notables des champs d'écoulement et de température, le transfert de chaleur moyen (dans le temps) évalué sur un cycle est approximativement égal à la valeur obtenue lorsque la paroi chaude est maintenue à une température constante.

Lage et al., [15] ont étudié numériquement et théoriquement la convection naturelle dans une cavité carrée refroidie isothermiquement sur un coté vertical alors que la paroi opposée est chauffée par un flux thermique pulsatoire. Ils ont montré que, pour une fréquence qui impose le flux de chaleur pulsé, la convection naturelle induite dans la cavité entre en résonance avec la température aux limites. Cette résonance est caractérisée par des fluctuations maximales du taux de transfert de chaleur à travers un plan vertical passant par le centre géométrique de la cellule. Ils ont observé également, en se basant sur les résultats numériques, que le phénomène de résonance existe dans la marge [0.01 – 7] du nombre de Prandtl et qu'il s'amplifie à mesure que le nombre de Rayleigh augmente.

En ce qui concerne le cas d'une cellule chauffée par le bas d'une façon transitoire, il semble que la première investigation faite sur cette configuration soit due à Mantle-Miller et al., [16]. Dans cette étude expérimentale, les auteurs ont examiné l'écoulement et le transfert de chaleur en convection naturelle dans le cas où la température de la paroi chaude varie périodiquement dans le temps autour d'une valeur moyenne, la paroi supérieure froide étant maintenue à une température constante et les parois verticales adiabatiques.

L'amplitude des oscillations de la température imposée varie entre 3 et 0% de la différence des températures moyennes sur les parois thermiquement actives. Pour les variations faibles de la température excitatrice, il est montré que le transfert de chaleur moyen sur un cycle est comparable au transfert de chaleur permanent calculé à partir de la température moyenne de la paroi chaude. Cependant, pour les grandes modulations de la température, on peut observer les différences de l'ordre de 12 % entre ces deux valeurs.

Plus récemment, Lakhal et al., [17, 18] ont étudié numériquement l'effet d'un chauffage variable sur la convection naturelle transitoire conduite dans une cavité carrée chauffée partiellement par le bas. L'effet de plusieurs paramètres dominants, à savoir: la position et la longueur de l'élément chauffant, l'amplitude et la période de la température excitatrice et le nombre de Rayleigh sur le transfert de chaleur et sur l'écoulement du fluide est analysé. Les conclusions principales concernent la période et l'amplitude de la quantité de chaleur qui quitte le système à travers sa paroi froide.

Contrairement aux résultats obtenus par Kazmierczak et al., [14], il semble que les solutions trouvées n'oscillent avec la même période que celle de la température excitatrice que pour les nombre de Rayleigh modérés. De plus pour une combinaison donnée des différents paramètres géométriques et thermodynamiques, une double périodicité de l'écoulement du fluide et du transfert de chaleur a été obtenue.
## 2. CONFIGURATION GEOMETRIQUE ET METHODE NUMERIQUE

La configuration étudiée dans ce travail est schématisée par la figure.1. Les parois horizontales inférieure et supérieure de la cavité sont portées respectivement aux températures chaude et froide. Les parois verticales sont portées à la même température que celle du bas, sauf une zone considérée adiabatique du côté de la paroi froide. Cette zone de hauteur égale à $L_{ref}/4$ est conçue pour relaxer le déplacement de l'interface solide/liquide dans les configurations de croissance cristalline à tirage vertical. L'étude suivante portera sur un fluide, de nombre de Prandlt $Pr = 0.01$ et la valeur de Rayleigh prise entre $10^3$ et $10^5$.

<figure><img src="image_1.png" /><figcaption>Fig. 1: a- Configuration géométrique et conditions aux limites
b- Profil des températures froide et chaude appliquées aux parois</figcaption></figure>

### 2.1 Mise en équations

L'écoulement résultant est supposé bidimensionnel (troisième dimension de la cavité suffisamment grande) et laminaire, et les propriétés du fluide sont considérées constantes. On néglige l'échauffement du à la viscosité dans l'équation d'énergie car les vitesses en convection naturelle sont modérées. La variation de la masse volumique suivant la température est calculée en utilisant l'approximation de Boussinesq.

Les équations gouvernant ce problème sont ceux de Navier-Stokes associées à l'équation de l'énergie en régime non stationnaire.

Les nombres de Prandtl $Pr$ et de Rayleigh $Ra$ sont définis par:

$$Pr = v / \alpha \qquad Ra = ( g \times \beta \times \Delta T \times H'3 ) / ( v \times \alpha)$$

-**Equation de la conservation de la masse**

$$\frac{\partial U}{\partial x} + \frac{\partial V}{\partial y} = 0 \qquad (1)$$

-**Equation de la quantité de mouvement**

$$\frac{\partial U}{\partial t} + U \times \frac{\partial U}{\partial x} + V \times \frac{\partial U}{\partial y} = \frac{\partial P}{\partial x} + Pr \times \left( \frac{\partial^2 U}{\partial x^2} + \frac{\partial^2 U}{\partial y^2} \right) \qquad (2)$$

$$\frac{\partial V}{\partial t} + U \times \frac{\partial V}{\partial x} + V \times \frac{\partial V}{\partial y} = - \frac{\partial P}{\partial x} + Pr \times \left( \frac{\partial^2 U}{\partial x^2} + \frac{\partial^2 U}{\partial y^2} \right) \qquad (3)$$

-**Equation d'énergie**
$$ \frac{\partial T}{\partial t} + U \times \frac{\partial T}{\partial x} + V \times \frac{\partial T}{\partial y} = - \frac{\partial^2 T}{\partial x^2} + \frac{\partial^2 T}{\partial y^2} \quad (4) $$

## 2.2 Méthode numérique

Les équations (1) à (4) ont été adimensionnées en utilisant comme échelle caractéristique la hauteur H' de la cellule, la différence de température ($\Delta T = T'_c - T'_f$) entre les parois thermiques actives, la vitesse ($\alpha / H'$) et le temps ($H'_\alpha / \alpha$), les équations primées désignent les grandeurs dimensionnelles.

Les conditions aux limites hydrodynamiques sont caractérisées par le non glissement des particules fluides sur les parois solides et l'imperméabilité de celle-ci, ce qui donne $u = v = 0$ sur toutes les parois.

Les conditions aux limites relatives à ce problème sont données par:

$T = T_c \text{ pour- } y = 0; x = 0, 1 \text{ et } y \ [0, 0.75]$

$T = T_c \text{ pour- } y = 1$

$\partial T / \partial x = 0 \text{ pour- } x = 0, 1 \text{ et } y \ [0, 0.75]$

$\partial T / \partial n = 0$

Pour les portions adiabatiques,

$u = v = 0 \text{ pour- } x = (0, 1) \text{ et } y \ [0,1]$

$u = v = 0 \text{ pour- } y = (0, 1) \text{ et } x \ [0,1]$

Les températures $T_f$ et $T_c$ sont définies par:

$$ T_f = A_f \times \sin(2\pi \times f_f \times t + \varphi_f), \text{ et } T_c = A_c \times \sin(2\pi \times f_c \times t + \varphi_c) $$

Avec, $A(f,c)$, $f(f,c)$ et $\varphi(f,c)$ sont successivement les amplitudes, les fréquences et les déphasages des températures froide et chaude.

La résolution numérique du système d'équations est effectuée par une méthode de volumes finis utilisant un schéma d'Euler implicite d'ordre 2 pour la discrétisation temporelle. Certaines variables peuvent présenter de forts gradients, ce qui nécessite l'utilisation d'un schéma de discrétisation spatiale qui soit à la fois stable et non diffusif.

Pour cette raison, nous avons utilisé un schéma de type *Ultimate2D* qui associe un schéma convectif décentré quadratique avec une stratégie de limitation de flux. La résolution du couplage vitesse - pression est assurée par l'algorithme *Simplec*. Le choix d'un pas de temps $\Delta t = 10^{-3}$ et un maillage de 64×64 dans le plan (x, y) a permis d'avoir une description précise des phénomènes thermo-convectifs au sein de la cavité.

A chaque pas de temps, les champs dynamique et thermique sont calculés itérativement jusqu'à ce que le critère de convergence suivant soit satisfait: valeur maximale relative entre deux itérations successives inférieure à 10.

# 3. RESULTATS ET DISCUSSIONS

Les résultats sont présentés en termes de structure d'écoulement et transfert de chaleur (Nombre de Nusselt).
## 3.1 Résultats relatifs au mode stationnaire

Avant de présenter les résultats correspondant au régime oscillatoire le champ du transfert thermique et de l'écoulement en régime stationnaire sont décrites brièvement.

Les lignes de champs et les courbes relatives au transfert thermique et la fonction courant maximale sont présentés dans la figure 2 pour différentes valeurs de Ra ; 10³, 10⁴ et 10⁵.

On note que lorsque les conditions aux limites sont introduites les solutions obtenues atteignent des régimes d'écoulement établi après une phase de transition dont la durée dépend du nombre de Rayleigh. Ils passent par un état stationnaire intermédiaire.

L'écoulement reste bicellulaire (Fig. 2), faiblement convectif pour les valeurs modérées de Rayleigh (Ra =10³) avec l'augmentation de celui-ci. L'écoulement devient notable avec les grandes valeurs de Ψmax liée à la transition de l'écoulement bicellulaire au monocellulaire dominé par une cellule tournant dans le sens trigonométrique.

Pour tous les nombres de Rayleigh comprise entre 10⁴ et 10⁵, l'écoulement stationnaire est formé d'une cellule négative, tournant dans le sens trigonométrique. Au delà de cette limite supérieure de Ra, une cellule de recirculation positive, représentant l'écoulement secondaire et tournant dans le sens horaire, apparaît sur l'un des coins inférieure ou supérieure de la cavité.

Cependant, son intensité reste très faible devant celle de la cellule principale dans la marge de Ra confédérée (Ra < 10⁶).

## 3.2 Résultats relatifs au mode oscillatoire

C'est dans le but de concentrer l'étude sur l'influence des fréquences de modulation des températures chaudes et froides et du nombre de Rayleigh sur l'écoulement et le transfert thermique que les valeurs de A<sub>f</sub>, A<sub>c</sub> et Pr sont considérées constantes et égales respectivement à 1, 1 et 0.01.

Les résultats seront présentés en termes de Ψ<sub>max</sub>(t) et Nu(t) désignant respectivement, la fonction de courant maximale et le nombre de Nusselt.

Dans le but d'étudier l'effet des oscillations des températures chaude et froide appliquées successivement aux parois inférieure et supérieure sur le transfert thermique et l'écoulement, les équations suivantes sont utilisées:

$$G(\Phi) = \Phi_{\text{moy}} = (\Phi_{\text{min}} + \Phi_{\text{max}})/2$$

$$A(\Phi) = \Phi_{\text{amp}} = (\Phi_{\text{max}} - \Phi_{\text{min}})/2$$

Où Φ représente une grandeur physique quelconque.

Les équations de Navier-Stocks couplées à l'équation de l'énergie ont été résolues numériquement en utilisant la technique du volume fini.

Dans cette section, les valeurs de A<sub>f</sub> et A<sub>c</sub> sont maintenus, le nombre de Rayleigh variant de 10³ à 10⁵, on considère l'effet de différentes valeurs des fréquences des températures froide et chaude sur l'écoulement et le transfert thermique.

Quand nous observons les courbes relatives à la fonction courant maximal en fonction du temps représentés dans la figure 4 pour différentes valeurs de la fréquence de modulation des températures froides et chaudes, nous constatons que toutes les solutions ont la même forme indépendamment de ces fréquences.
Ils passent tous par un régime transitoire caractérisé par un état intermédiaire stationnaire avant de passer au régime établi. La durée de ce régime dépend des valeurs des fréquences de modulation des températures $T_c$ et $T_f$. Concernant le temps adimensionnel nécessaire au développement du régime permanent, il est montré que ce dernier dépend de $R_a$ et des fréquences de modulation des températures de chauffe et de refroidissement.

Les courbes indiquent également que toutes les solutions sont périodiques et oscillent autour de leurs valeurs moyennes égale à celle du régime stationnaire. La fréquence de ces oscillations est égale au minimum des deux fréquences $f_c$ et $f_f$, tandis que le maximum assurent la modulation en amplitude de la solution.

Nous présentons dans la figure 3, l'écoulement sous forme de lignes de courant pour le cas correspondant à la fréquence de modulation de la température froide $f_c=10$ et ce pour différentes valeurs de la fréquence de modulation de la température chaude $f_c$.

Nous constatons que pour $R_a=10^3$, les lignes de courants sont sous forme de deux cellules symétriques tournantes dans deux sens contraires.

Pour $R_a=10^4$, l'écoulement est dominé par une seul cellule dite cellule principale tournant dans le sens trigonométrique pour les valeurs basse de la fréquence $f_c$, au fur et à mesure que cette dernière augmente, nous assistons à l'apparition d'une cellule secondaire dans le coin inférieur gauche de la cavité. Cette cellule secondaire prenne de l'ampleur et augmente au détriment de la cellule principale qui disparait pour laisser toute la place à celle-ci pour les valeurs élevées de la fréquence $f_c$.

En ce qui concerne le cas du $R_a=10^5$, la structure du flux est caractérisée par la prédominance d'une seule cellule principale avec l'apparition des petites cellules dans les quatre coins de la cavité. Ces cellules disparaissent complètement pour les valeurs très élevées de la fréquence $f_c > 60$.

On présente dans la figure (3), les variations de $Nu(t)$ et $\Psi_{max}(t)$ pour $R_a = 10^4$.

Sur les courbes, on note, en régime oscillatoire établi, les solutions correspondantes à $\Psi_{max}(t)$ et $Nu(t)$ oscillent de manière sinusoïdal dans le temps avec une valeur moyenne et une amplitude qui varie en fonction des fréquences des températures $f_c$ et $f_f$.

Nous constatons que pour les valeurs basses et élevées de $f_c$ et quelque soit la valeur de $f_f$, la valeur moyenne de $\Psi_{max}$ (Fig. 4a) est égale à sa valeur en mode stationnaire.

L'amplitude des oscillations de la fonction courant maximale varie en fonction de $f_c$ et $f_f$. Pour les valeurs fixes de la fréquence $f_f$, cette amplitude diminue avec la fréquence $f_c$.

Concernant le nombre de Nusselt, l'amplitude des oscillations diminue avec la fréquence $f_c$ pour $f_c < 20$. Au-delà de cette valeur, l'amplitude reste quasiment constante et indépendante de $f_c$, tandis qu'elle augmente avec la fréquence $f_f$, quelque soit la valeur de $f_c$.
$Ra = 10^3$

$Ra = 10^4$

$Ra = 10^5$

<figure><img src="image_2.png" /><figcaption>(a)</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>(b)</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>Fig. 2: a- Forme des cellules en régime stationnaire en fonction de Rayleigh pour $f_f = f_c = 0$. b- $\Psi_{max}(t)$ et $Nu(t)$ en régime stationnaire pour $f_f = f_c = 0$ et différentes valeurs de $R_a$</figcaption></figure>

$f_c = 1$

<figure><img src="image_5.png" /></figure>

$f_c = 5$
<figure><img src="image_6.png" /></figure>

$f_c = 10$

<figure><img src="image_7.png" /></figure>

$f_c = 50$

<figure><img src="image_8.png" /><figcaption>Fig. 3: Lignes de courant pour Ra = 10³, 10⁴ et 10⁵
pour f<sub>f</sub> = 10 et différentes valeurs de f<sub>c</sub></figcaption></figure>

$f_f = 0$

<figure><img src="image_9.png" /></figure>

$f_f = 10$
<figure><img src="image_10.png" /></figure>

$f_f = 50$

<figure><img src="image_11.png" /></figure>

$f_f = 100$

<caption>Fig. 4: Variations de $\Psi_{max}(t)$ et $Nu(t)$ dans le temps pour $Ra = 10^4$ et pour $f_f = 10$ et différentes valeurs de $f_f$ et $f_c$</caption>

## 4. CONCLUSION

L'étude numérique de la convection naturelle au sein d'une cavité carrée avec chauffage et refroidissent variables a permis de formuler les conclusions suivantes :

1- Toutes les solutions obtenues sont périodiques dans le temps, leurs fréquence est égale au minimum des fréquences des températures $T_c$ et $T_f$, $min(f_c, f_f)$.

2- Les solutions passent par un régime transitoire, caractérisé par un état intermédiaire stationnaire avant d'atteindre le régime établi.

3- La durée du régime transitoire dépend de la valeur du nombre de Rayleigh et des fréquences de modulation des températures de chauffe et de refroidissement.

4- Généralement les valeurs moyennes du transfert thermique et de la fonction courant maximale sont différentes de celles obtenues en régime stationnaire.

## NOMENCLATURE

<table><tr><td>g</td><td>Accélération gravitationnelle</td><td>Pr</td><td>Nombre de Prandtl</td></tr><tr><td>t</td><td>Temps adimensionné, (t' / (H<sup>2</sup> / α))</td><td>Ψ</td><td>Fonction de courant adimensionné, (Ψ' / α)</td></tr><tr><td>P</td><td>Pression</td><td>τ</td><td>Période de la température</td></tr></table>
<table><tr><td>t</td><td>Temps adimensionnel</td><td>Nu</td><td>Nombre de Nusselt</td></tr><tr><td>T</td><td>Température adimensionnelle</td><td>Q</td><td>Flux thermique, W/m²</td></tr><tr><td>RaT</td><td>Nombre de Rayleigh thermique</td><td>h</td><td>Coefficient de convection local, W/(m².K)</td></tr><tr><td>Ra</td><td>Nombre de Rayleigh, (gβ ΔT' H³ / 9α)</td><td>U</td><td>Vitesse adimensionnelle suivant la direction x</td></tr><tr><td>V</td><td>Vitesse adimensionnelle suivant la direction y</td><td>Ω</td><td>Vorticoté adimensionnée (W H² / α)</td></tr><tr><td>x, y</td><td>Coordonnées cartésiennes adimensionnelles</td><td>μ</td><td>Viscosité dynamique, (Kg/ms)</td></tr><tr><td>H'</td><td>Hauteur de la cavité</td><td>ω'</td><td>Viscosité (s⁻¹)</td></tr><tr><td>L'</td><td>Largeur de la cellule</td><td>θ</td><td>Viscosité cinématique, (m²/s)</td></tr><tr><td>QF</td><td>Densité de flux thermique dimensionnée à travers paroi froide</td><td>C</td><td>Fraction massique</td></tr><tr><td></td><td></td><td>D</td><td>Coefficient de diffusion rapide</td></tr><tr><td>T</td><td>Température adimensionnée (T' - T'f) / ΔT'</td><td>ΔT'</td><td>Différence de température T'c - T'f</td></tr><tr><td>a, b</td><td>Diffusivité thermique Coefficient d'expansion thermique ou solutal</td><td>c, f</td><td>Chaud, Froid</td></tr><tr><td></td><td></td><td>T</td><td>Thermique</td></tr></table>

REFERENCES

[1] G. Hensen, A. Hellawell, S.Z. Lu and R.S. Steube, ‘Some Consequences of Thermosolutal Convection: the Grain Structure of Castings’, Metallurgical and Mate

rials Transactions A., 27A, N°3, pp. 569 – 581, 1996.

[2] C.J. Chang and R.A. Brown, ‘Radial Segregation Induced by Natural Convection and Melt-Solid Interface Shape in Vertical Bridgman Growth’, Journal Crystal Growth, Vol. 61, N°1, pp. 343 – 364, 1983.

[3] A. Bejan, ‘Convection Heat Transfer’, 477 p., Wiley, New York, 1984.

[4] S. Ostrach, ‘Natural Convection in Enclosures’, Advances in Heat Transfer, Vol. 8, pp. 161 – 227, 1972.

[5] N.N. Lin and A. Bejan, ‘Natural Convection in a Partially Divided Enclosure’, International Journal of Heat and Mass Transfer, Vol. 26, N°12, pp. 1867 – 1878, 1983.

[6] S. Kimura and A. Bejan, ‘The Boundary Layer Natural Convection Regime in a Rectangular Cavity with Uniform Heat Flow from Side’, Journal of Heat Transfer, Vol. 106, N°1, pp. 96 – 103, 1984.

[7] V.F. Nicolette, K.T. Yand and J.R. Lloyd, ‘Transient Cooling by Natural Convection in a Two-Dimensional Square Enclosure’, International Journal of Heat and Mass Transfer, Vol. 28, N°9, pp. 1721 – 1732, 1985.

[8] D.G Briggs and D.N. Jones, ‘Two-Dimensional Periodic Natural Convection in a Rectangular Enclosure of Aspect Ratio One’, Journal of Heat Transfer, Vol. 107, N°4, pp. 850 – 854, 1985.

[9] S. Kimura, K. Kamomia and T. Sato, ‘Oscillatory Convection in Rectangular Cavity’, Bulletin of the Japon Society of Mechanical Engineers, Vol. 106, pp. 104 – 108, 1984.
[10] D. Poulikakos, 'Natural Convection in a Confined Fluid Space Driven by a Single Vertical Wall with Warm and Cold Region', Journal of Heat Transfer, Vol. 107, pp. 867–876, 1985.

[11] J.C. Patterson, 'On the Existence of an Oscillatory Approach to Steady Natural Convection in Cavities', Journal of Heat Transfer, Vol. 106, pp. 104–108, 1990.

[12] J.C. Patterson and J. Imberger, 'Unsteady Natural Convection in a Rectangular Cavity', International Journal of Fluid Mechanics, Vol. 100, N°1, pp. 65–86, 1990.

[13] H.Q. Yang, K.T. Yang and Q. Xia, 'Periodic Laminar Convection in a Vertical Cavity', International Journal of Heat and Mass Transfer, Vol. 32, N°11, pp. 2199–2207, 1989.

[14] M. Kazmierczak and Z. Chinoda, 'Buyabcy-Driven Flow in an Enclosure with Time Periodic Boundary Conditions', International Journal of Heat and Mass Transfer, Vol. 32, N°6, pp. 1507–1519, 1992.

[15] J.I. Lage and A. Bejan, 'The Resonance of Natural Convection in a Horizontal Enclosure Heated Periodically from the Side', International Journal of Heat and Mass Transfer, Vol 36, N°6, pp. 2027–2038, 1992.

[16] W.J. Mantle-Miller, M. Kazmierczak and B. Hiawy, 'Natural Convection in a Horizontal Enclosure with Periodically Changing Bottom Wall Temperature', ASME 28th National Heat Transfer Conference, HTD-Vol, 198, Natural Convection in Enclosure, pp. 49 - 56, 1992.

[17] E.K. Lakhal, M. Hasnaoui, P. Vasseur et E. Bilgen, 'Convection Naturelle dans une Cellule Carrée Chauffée Périodiquement par le Bas: Etude Numérique', Revue Générale de Thermique, Fr. pp. 392–393, 1994, pp. 480–485, 1994.

[18] E.K. Lakhal, M. Hasnaoui, P. Vasseur and E. Bilgen, 'Natural Convection in a Square Enclosure Heated Periodically From Part of Bottom Wall', Numerical Heat Transfer, Part A, Vol. 27, N°3, pp. 319–333, 1995.

[19] H. Rouijaa, M. El Alami, E.A. Semma and M. Najam, 'Natural Convection in an Inclined T Shaped Cavity', FDMP Journal, Vol. 235, N°1, pp. 1 - 14, 2010.