# Simulation numérique du transfert de chaleur mixte dans une cavité carrée: effet de la présence des chicanes

Mohamed Amine Belmiloud * et Noureddine Sad Chemloul

Laboratoire de la Recherche et des Technologies Industrielles
Département de Génie Mécanique
Université Ibn Khaldoun, B.P. 78, 14000 Tiaret, Algérie.

(reçu le 2 Février 2016 – accepté le 28 Mars 2016)

**Résumé** - Dans cette étude, nous avons simulé numériquement l'influence de la vitesse de rotation $\omega$, du positionnement des chicanes et des cisaillements des parois sur le transfert de chaleur par convection mixte dans une cavité carrée. La paroi verticale gauche est chauffée à une température constante plus élevée que celle de la paroi verticale droite $T_h > T_c$ et les deux parois horizontales sont adiabatiques. La méthode des volumes finis est utilisée pour résoudre les équations gouvernantes dimensionnelles. La vitesse de cisaillement est déterminée pour un nombre de Richardson $Ri = 10$, la vitesse de rotation du corps solide est comprise entre 10 et 40 rd/s. Le nombre de Grashof Gr est fixé à $10^5$ et le nombre de Prandtl Pr est maintenu constant 0.71. Les résultats obtenus dans cette étude montrent que, l'influence de la vitesse de rotation $\omega$ est présente pour les configurations avec chicane, en outre, le nombre de Nusselt moyen $Nu_{moy}$ prend des valeurs maximales pour la configuration sans chicanes.

**Abstract** - In this study, we numerically simulated the influence of the rotational velocity $\omega$, the positioning of baffles and wall shear on heat transfer by mixed convection in a square cavity. The left vertical wall is heated to a constant temperature higher than the right vertical wall $T_h > T_c$ and two horizontal walls are adiabatic. The finite volume method is used to solve the governing equations dimensional. The shear rate is determined for a Richardson number $Ri = 10$, the velocity of rotation of the solid body is between 10 and 40 rd/s. The Grashof number Gr is fixed to $10^5$ and the Prandtl number Pr is kept constant 0.71. The results obtained in this study show that the influence of the velocity of rotation is present for configurations with baffle, in addition, the number of average Nusselt $Nu_{moy}$ takes maximum values for configuration without baffles.

**Mots clés:** Simulation numérique - Convection mixte - Cisaillement des parois - Vitesse de rotation - Chicanes.

## 1. INTRODUCTION

Le problème fondamental de transfert de chaleur par convection mixte dans une cavité carrée avec des parois cisailles a reçu une attention considérable de la part des chercheurs. Ce problème est souvent rencontré dans de nombreuses applications d'ingénierie. Ces types de problèmes se posent également dans les boîtiers électroniques, les fours industriels.

Dans leur recherche, Iwatsu et al., [1] et Iwatsu et al., [2] réalisent respectivement une simulation numérique à deux et à trois dimensions de la convection mixte dans une cavité carrée chauffée dont la paroi supérieure est mobile. La convection mixte dans une cavité avec la variation des différentes conditions aux limites a été étudiée par Hasnaoui et al., [3].

Lee et al., [4] obtiennent des solutions avec la méthode des éléments finis de la convection mixte dans une cavité carrée chauffée par la paroi inférieure. Moraga et al.,

<footnote>* Belmiloud_amine@yahoo.fr</footnote>
[5] effectuent une analyse numérique d'un modèle tridimensionnel de la convection mixte dans une cavité refroidie par air, afin de comparer les variations des différentes propriétés avec les résultats des modèles bidimensionnels.

Sharif [6] a étudié le processus de la convection mixte laminaire dans des cavités rectangulaires bidimensionnelles et pour trois différentes valeurs du nombre de Richardson Ri représentant la dominance de convection forcée, de convection mixte, et de la convection naturelle en utilisant le code commercial Fluent. Les effets de l'inclinaison de la cavité sur les processus résultant de convection sont également étudiés. L'augmentation du nombre de Nusselt avec l'inclinaison de la cavité est moyenne pour la convection forcée et significative pour la convection naturelle.

Oztop *et al.*, [7], ont étudié le transfert thermique par convection mixte dans une cavité à deux parois verticales mobiles. Les résultats obtenus montrent que, le nombre de Richardson influe sur le transfert thermique, où pour des petites valeurs de Ri, le nombre de Nusselt moyen est diminué, par contre, pour Ri augmente le transfert thermique est augmenté. Dans le cas des deux parois verticales montant vers le haut, le transfert est faible par rapport à celui correspondant au cisaillement des parois.

Al-Amiri *et al.*, [8], ont étudié le transfert thermique par convection mixte dans une cavité à paroi mobile avec une surface de fond ondulée. Les résultats de ce travail ont montré que le nombre de Nusselt moyen augmente avec l'augmentation de l'amplitude de l'ondulation et du nombre de Reynolds Re. Ils ont montré que le transfert de chaleur optimal est atteint pour un petit nombre de Richardson.

Mekroussi *et al.*, [9] ont à étudier numériquement la convection mixte dans une cavité carrée, utilisant une ondulation sur la paroi inférieure de la cavité. Les résultats, lignes de courants et des isothermes, sont donnés pour différents angles (0°, 45°, 90°, 120°, 150°, 180°). Belmiloud *et al.*, [10] ont étudié numériquement la convection mixte couplée à une radiation dans une cavité carrée avec la paroi supérieure mobile. Leurs résultats ont montré que la vitesse de circulation augmente et le nombre de Nusselt convectif diminue, lorsque l'émissivité de la paroi de la cavité augmente.

Le mode de transfert de chaleur dominant est la convection forcée, lorsque les valeurs du nombre de Richardson Ri sont très proches de celles correspondantes à la convection forcée $Ri > 0.1$.

Lorsque le nombre de Richardson Ri augmente, le mode de transfert de chaleur dominant est la convection naturelle. La vitesse du couvercle de la cavité influe sur le transfert thermique ceci est prouvé par la diminution du nombre de Nusselt Nu lorsque le nombre de Richardson Ri augmente.

## 2. FORMULATION DU PROBLEME

La valeur de la vitesse de rotation $\omega$ pour le corps solide est comprise entre 10 et 40 $rd/s$, le nombre de Grashof $Gr = 10^5$ et la vitesse de cisaillement est déterminée à partir du nombre de Richardson $Ri = 10$.

La géométrie considérée est une cavité carrée dont les deux parois horizontales sont adiabatiques, et les deux parois verticales sont différentiellement chauffées $T_c < T_h$ avec un écart de la température de 10°C. La longueur des chicanes est de 0.5H. Le rapport d'aspect de la cavité est de $A = 1$, le rapport entre le corps solide et la cavité est de 0.235 et la température de référence est prise égale à celle du fluide $T_0 = 293.15$ K.
<figure><img src="image_1.png" /><figcaption>Fig. 1: Problème schématisé: -a- Cisaillement horizontal, -b- Cisaillement vertical</figcaption></figure>

Les équations générales de conservation à savoir, l'équation de continuité, les équations de mouvement suivant x et y et l'équation de l'énergie, sont données sous leur forme adimensionnelle comme suit:

$$ \frac{\partial u^+}{\partial x^+} + \frac{\partial v^+}{\partial y^+} = 0 \qquad u \qquad (1) $$

$$ u^+ + \frac{\partial u^+}{\partial x^+} + v^+ \frac{\partial v^+}{\partial y^+} = - \frac{\partial P^+}{\partial x^+} + \frac{1}{\text{Re}} \left( \frac{\partial^2 u^+}{\partial x^{+2}} + \frac{\partial^2 u^+}{\partial y^{+2}} \right) + \text{Ri} T^+ \quad (2) $$

$$ u^+ + \frac{\partial v^+}{\partial x^+} + v^+ \frac{\partial v^+}{\partial y^+} = - \frac{\partial P^+}{\partial y^+} + \frac{1}{\text{Re}} \left( \frac{\partial^2 v^+}{\partial x^{+2}} + \frac{\partial^2 v^+}{\partial y^{+2}} \right) \quad (3) $$

$$ u^+ + \frac{\partial T^+}{\partial x^+} + v^+ \frac{\partial T^+}{\partial y^+} = \frac{1}{\text{Re.Pr}} \left( \frac{\partial^2 T^+}{\partial x^{+2}} + \frac{\partial^2 T^+}{\partial y^{+2}} \right) \quad (4) $$

Dans les équations ci-dessus, les variables sans dimension sont définies par:

$$ x^+ = \frac{x}{H}, \quad y^+ = \frac{y}{H}, \quad u^+ = \frac{u}{u_{liq}}, \quad v^+ = \frac{v}{u_{liq}}, \quad P^+ = \frac{P}{\rho \cdot u_{liq}}, \quad T^+ = \frac{T - T_c}{T_h - T_c} \quad (5) $$

Le nombre de Grashof Gr, le nombre de Reynolds Re, le nombre de Richardson Ri et le nombre de Prandtl Pr sont donnés par:

$$ \text{Gr} = \frac{g\beta(T_h - T_c) \cdot H^3}{v^2}, \quad \text{Ri} = \frac{\text{Gr}}{\text{Re}^2}, \quad \text{Re} = \frac{u_{liq} \cdot H}{v}, \quad \text{Pr} = \frac{v}{\alpha} \quad (6) $$

Pour résoudre les équations {Eq. 2} – {Eq. 4}, on considère les conditions aux limites adimensionnelles suivantes:

A la paroi verticale gauche: $T^+ = 1$; $u^+ = 0$; $v^+ = 0$ ou $v^+ = 1$

A la paroi verticale droite: $T^+ = 0$; $u^+ = 0$; $v^+ = 0$ ou $v^+ = -1$

A la paroi horizontale inférieure: $\frac{\partial T^+}{\partial y^+} = 0$; $v^+ = 0$; $u^+ = 0$ ou $u^+ = -1$
Pour déterminer les caractéristiques du transfert de chaleur à l'état stationnaire, il faut prendre en considération la contribution de la convection.

Dans la présente étude, le nombre de Nusselt moyen total est défini comme:

$$Nu = \frac{Q \cdot H}{\lambda \cdot (T_h - T_c)} \qquad (7)$$

La forme adimensionnelle de l'équation {Eq. 7} est donnée par:

$$Nu = \int_{0}^{1} -\frac{1}{T^{+}} \frac{\partial T^{+}}{\partial x^{+}} dy^{+} \qquad (8)$$

### 3. NUMERICAL ANALYSIS

La résolution des équations de conservation de la masse, de mouvement et de l'énergie par l'algorithmes des différences finies appelé 'Simple'. Les détails de cette méthode est décrite par Patankar [11]. Les différentes équations sont discrétisées dans un volume de contrôle. Le schéma différence (PLDS) utilisé pour le calcul des variables scalaires et le schéma 'Quick' pour les variables vectorielles est présenté par Hayase et al., [12]. La valeur de la convergence considérée dans cette étude est:

$$\frac{|\varphi_n - \varphi_{n-1}|_{\max}}{|\varphi_n|} < 10^{-5} \qquad (9)$$

Le **Tableau 1** présente les valeurs du nombre de Nusselt moyen $Nu_{moy}$ déterminées dans le cas des chicanes horizontales et sans effet de cisaillement. On remarque que la différence relative de $Nu_{moy}$ obtenue est très faible. Le maillage utilisé dans tous les calculs ultérieurs est (84 × 84).

<table><caption>**Tableau 1:** Effet du maillage sur les résultats obtenus pour la configuration des chicanes horizontales et sans cisaillement</caption><thead><tr><th></th><th></th><th></th><th></th><th colspan="3">Erreur (abs) %</th></tr></thead><tbody><tr><td>Maillage</td><td>68 × 68</td><td>76 × 76</td><td>84 × 84</td><td>68 × 68</td><td>76 × 76</td><td>84 × 84</td></tr><tr><td>NU<sub>moy</sub></td><td>1.8770</td><td>1.8769</td><td>1.8761</td><td>0.0466</td><td>0.0454</td><td>-</td></tr></tbody></table>

#### 3.1 Validation de code

Pour la validation numérique du code de calcul, nous avons comparé nos résultats avec ceux de Figueredo *et al.*, [13], et Souza [14]. Le tableau présente les valeurs du nombre de Nusselt moyen $Nu_{moy}$ déterminée au niveau de la paroi froide. Nous remarquons que les valeurs du nombre de Nusselt moyen $Nu_{moy}$ obtenues dans la présente étude et celles obtenues par [14] sont très proches avec une erreur maximale de 0.996 %. Par contre la valeur relative maximale obtenue par [13] est 4.502 %.

<table><caption>**Tableau 2:** Nombre de Nusselt moyen $Nu_{moy}$ évalué sur la paroi froide</caption><thead><tr><th></th><th>Présente étude</th><th>Figueredo <em>et al.</em>, [13]</th><th>Erreur (abs) %</th><th>Souza [14]</th><th>Erreur (abs) %</th></tr></thead><tbody><tr><td>34110</td><td>2.941</td><td>2.884</td><td>1.976</td><td>2.912</td><td>0.996</td></tr><tr><td>60000</td><td>3.489</td><td>3.468</td><td>0.605</td><td>3.456</td><td>0.955</td></tr><tr><td>100000</td><td>4.074</td><td>4.160</td><td>2.067</td><td>4.038</td><td>0.891</td></tr><tr><td>136430</td><td>4.475</td><td>4.686</td><td>4.502</td><td>4.440</td><td>0.788</td></tr></tbody></table>
# 4. RESULTATS ET DISCUSSIONS

## 4.1 Lignes des courants et lignes isothermes

La figure 2 présente la variation des contours de vitesse et les lignes isothermes pour les différentes configurations (sans chicane, chicanes horizontales et chicanes verticales) et pour le cisaillement des parois verticales. Nous avons remarqué pour les contours de vitesse, des tourbillons au voisinage du corps solide avec des tailles différentes, cette taille varie en fonction de la position des chicanes. La valeur maximale est obtenue au voisinage du corps solide pour la configuration **-c-** 'chicanes verticales'.

Pour les lignes isothermes, on remarque que, la zone de fluide froide / chaude est plus grande dans la configuration **-a-** 'sans chicane' par rapport aux autres configurations **-b-** et **-c-**. La comparaison entre la configuration **-b-** 'chicanes horizontales' et la configuration **-c-** 'chicanes verticales' montre que la zone intermédiaire entre les deux plaques froide / chaude est plus grande dans le cas de la position des chicanes horizontales.

<figure><img src="image_2.png" /><figcaption>Fig. 2: Contours de vitesse et lignes isothermes en fonction de la position de la chicane pour le cisaillement vertical, et pour un corps solide adiabatique, ($\omega$ = 30 rd/s)</figcaption></figure>
La variation des contours de vitesse et des lignes isothermes pour les deux positions de cisaillement (parois horizontales et parois verticales) et pour la vitesse de rotation ($\omega$=20 rd/s) est représentée en figure 3. On remarque que les tourbillons apparaissent au voisinage du corps solide suivant le mode de cisaillement. La vitesse maximale est obtenue dans le cas où le cisaillement est provoqué par les parois verticales. Les lignes isothermes montrent que le transfert de la chaleur est meilleur pour le mode de cisaillement vertical, ceci est dû à la poussée d'Archimède, où la vitesse de cisaillement et la circulation du fluide ont le même sens.

<figure><img src="image_3.png" /><figcaption>Fig. 3: Variation des contours de vitesse et les lignes isothermes pour les deux positions de cisaillement et pour ($\omega$ = 20 rd/s): -a- horizontal, -b- vertical</figcaption></figure>

## 4.2 Variation du transfert de chaleur

La figure 4 montre la variation du nombre de Nusselt moyen $Nu_{moy}$ en fonction des différentes valeurs de la vitesse de rotation $\omega$ et pour les différentes configurations des cavités (sans chicane, avec chicanes horizontales, avec chicanes verticales). Tous les résultats de cette étude sont obtenus au niveau de la paroi verticale chaude.

Nous remarquons qu'en figure 4-a-, le nombre de Nusselt moyen $Nu_{moy}$ diminue avec l'augmentation de la vitesse de rotation $\omega$, entre autre, la position du cisaillement influe sur le nombre de Nusselt moyen $Nu_{moy}$ avec une valeur maximale déterminée pour le cisaillement vertical.

Les figures 4-b- et 4-c- montrent que le nombre de Nusselt moyen $Nu_{moy}$ augmente avec l'augmentation de la vitesse de rotation $\omega$. Par comparaison entre la configuration -b- et la configuration -c-, la valeur maximale de $Nu_{moy}$ est obtenue
pour la configuration **-b-**, d'autre part, le nombre de Nusselt $Nu_{moy}$ prend une valeur maximale, lorsque la chicane et le mode de cisaillement ont la même position.

La comparaison entre les trois configurations, montre que la configuration sans chicane permet d'obtenir un meilleur transfert de chaleur. Pour des cavités contenant des chicanes, l'augmentation du transfert de chaleur exige que les corps solides soient en rotation.

<figure><img src="image_4.png" /><figcaption>Fig. 4: Variation du nombre de Nusselt moyen $Nu_{moy}$ en fonction de $\omega$ pour les configurations: **-a-** sans chicane, **-b-** chicanes horizontales, **-c-** chicanes verticales</figcaption></figure>

La figure 5 montre la variation du nombre de Nusselt $Nu$, le long de $y^+$ pour les différentes vitesses $\omega$, les résultats sont obtenus au niveau de la paroi chaude. Nous remarquons qu'en figure 5-a-, les courbes du nombre de Nusselt $Nu$ sont presque confondues.

En figure 5-b-, les valeurs du nombre de Nusselt $Nu$ sont égales pour les différentes vitesses $\omega$, lorsque $-0.5 \le y^+ \le -0.35$. Pour $-0.35 \le y^+ \le 0.1$, le nombre de Nusselt $Nu$ augmente avec l'augmentation de $\omega$, au-delà les courbes presque confondues.

La figure 5-c-, montre que le nombre de Nusselt $Nu$ augmente avec l'augmentation de $\omega$ pour $-0.5 \le y^+ \le 0.25$, au-delà les trois courbes de $\omega = 20, 30, 40$ rd/s sont confondues. On constate que l'influence de la vitesse de rotation est significative pour le cas des chicanes verticales.
<figure><img src="image_5.png" /><figcaption>Fig. 5: Variation du nombre de Nusselt le long de $y^+$ pour les différentes vitesses de rotation $\omega$ et pour cisaillement horizontal
-a- sans chicane, -b- chicanes horizontales, -c- chicane verticales</figcaption></figure>

La variation du nombre de Nusselt moyen $Nu_{moy}$ en fonction des différentes valeurs de $\omega$ pour cavité sans chicane et pour les différentes températures du corps solide est représentée en figure 6.

<figure><img src="image_6.png" /><figcaption>Fig. 6: Variation du nombre de Nusselt moyen en fonction de $\omega$, pour les différentes températures du corps solide</figcaption></figure>
Par comparaison à un corps solide adiabatique, nous remarquons que la plus grande valeur du nombre de Nusselt moyen $Nu_{moy}$ est obtenue pour une température du corps $T_s = 288$ K. Le nombre de Nusselt moyen $Nu_{moy}$ diminue lorsque la température du corps solide augmente. On peut donc déduire que le meilleur transfert thermique est obtenu lorsque la température du corps est petite.

## 5. CONCLUSION

Le but de cette étude est d'analyse l'effet du cisaillement crée par le mouvement des parois, de la présence des chicanes et leurs positions, et l'influence de la vitesse de rotation du corps solide sur la convection mixte dans une cavité carrée. D'après les résultats obtenus, on conclut:

- Sans la présence des chicanes, le nombre de Nusselt moyen $Nu_{moy}$ diminue avec l'augmentation de la vitesse de rotation $\omega$, et la valeur maximale de $Nu_{moy}$ est obtenue pour un cisaillement vertical.

- En présence des chicanes, le nombre de Nusselt moyen $Nu_{moy}$ augmente avec l'augmentation de la vitesse de rotation $\omega$, et la valeur maximale de Nusselt $Nu_{moy}$ est obtenue dans le cas des chicanes horizontales.

- La comparaison entre les trois configurations montre que le bon transfert de chaleur est obtenu pour le cas sans chicanes. L'influence de la vitesse de rotation est significative dans le cas des chicanes verticales. La plus grande valeur du nombre de Nusselt moyen $Nu_{moy}$ est obtenue lorsque la température du corps est égale à 288 K.

## NOMENCLATURE

<table><tr><td>A, rapport d'aspect de la cavité, L/H</td><td>g, Accélération de gravité, m/s</td></tr><tr><td>β, Coef. dilatation thermique (fluide)</td><td>α, Diffusivité thermique (fluide), m²/s</td></tr><tr><td>Gr, Nombre de Grashof</td><td>ν, viscosité cinématique, m²/s</td></tr><tr><td>Nu, Nombre de Nusselt</td><td>λ, Conductivité thermique, W/mK</td></tr><tr><td>H, Hauteur de la cavité, m</td><td>ω, Vitesse de rotation, rd/s</td></tr><tr><td>Pr, Nombre de Prandtl</td><td>Ri, Nombre de Richardson</td></tr><tr><td>T, Température du fluide, K</td><td>T₀, Température de référence, K</td></tr><tr><td>u₁ᵢₗ, Vitesse de la paroi mobile, m/s</td><td>+, Valeurs adimensionnelle</td></tr><tr><td>x, y, Composantes de vitesse, m/s</td><td>x, y, Composantes cartésiennes, m/s</td></tr><tr><td>h, Chaud - c, Froid</td><td>h, moyen</td></tr></table>

## REFERENCES

[1] R. Iwatsu, J.M. Hyun and K. Kuwahara, 'Mixed Convection in a Driven Cavity with a Stable Vertical Temperature Gradient', International Journal of Heat and Mass Transfer, Vol. 36, N°6, pp. 1601-1608, 1995.

[2] R. Iwatsu and J.M. Hyun, 'Three-Dimensional Driven Cavity Flows with a Vertical Temperature Gradient', International Journal of Heat and Mass Transfer, Vol. 38, N°18, pp. 3319-3328, 1995.
[3] M. Hasnaoui, E. Bilgen and P. Vasseur, ‘Natural Convection Above an Array of Open Cavities Heated from Below’, AIAA Journal Thermophysics and Heat Transfer, Vol. 6, N°2, pp. 255 - 264, 1990.

[4] S.C. Lee and C.K. Chen, ‘Finite Element Solutions of Laminar and Turbulent Mixed Convection in a Driven Cavity’, International Journal for Numerical Methods in Fluids, Vol. 23, pp. 47–64, 1996.

[5] N. Moraga and S. López, ‘Numerical Simulation of Three Dimensional Mixed Convection in an Air-Cooled Cavity’, Numerical Heat Transfer, Part A, Vol. 45, N°8, pp. 811–824, 2004.

[6] M.A.R. Sharif, ‘Laminar Mixed Convection in Shallow Inclined Driven Cavities with Hot Moving Lid on Top and Cooled From Bottom’, Applied Thermal Engineering, Vol. 7, N°5-6, pp. 1036-1042, 2007.

[7] H.F. Oztop and I. Dagtekin, ‘Mixed Convection in Two-Sided Lid-Driven Differentially Heated Square Cavity’, International Journal of Heat and Mass Transfer, Vol. 47, N°4, pp. 1761–1769, 2004.

[8] S. Mekroussi, D. Nehari, M. Bouzit and N. Sad Chemloul, ‘Analysis of Mixed Convection in an Inclined Lid-Driven Cavity with a Wavy Wall’, Journal of Mechanical Science and Technology, Vol. 27, N°7, pp. 2181-2190, 2013.

[9] A. Al-Amiri, K. hanafer, J. Bull and I. Pop, ‘Effect of Sinusoidal Wavy Bottom Surface on Mixed Convection Heat Transfer in a Lid-Driven Cavity', International Journal of Heat and Mass Transfer, Vol. 50, N°9-10, pp. 1771 - 1780, 2007.

[10] M.A. Belmiloud and N. Sad Chemloul, ‘Numerical Study of Mixed Convection Coupled to Radiation in a Square Cavity with a Lid-Driven’, World Academy of Science, Engineering and Technology, International Journal of Mechanical, Aerospace, Industrial, Mechatronic and Manufacturing Engineering, Vol. 9, N°10, pp. 1758 – 1764, 2015

[11] S.V. Patankar, ‘Numerical Heat Transfer and Fluid Flow’, Hemisphere / McGraw-Hill, Washington D.C, 1980.

[12] T. Hayase, J. C. Humphrey and R. Greif, ‘A Consistently Formulated Quick Scheme for Fast and Stable Convergence Using Finite-Volume Iterative Calculation Procedures’, Journal of Computational Physics Archive, Vol. 98, N°1, pp. 108 – 118, 1992.

[13] J.R. Figueredo, M.M. Ganzarolli and P.I.F. Almeida, ‘Convecção Natural em Cavidades Retangulares – Solução Numérica’, II Congresso Latino-Americano de Transferência de Calor e Matéria, São Paulo, pp. 62 - 73, 5-10 October 1986.

[14] J.J. Souza, ‘Simulação Numérica da Transferência de Calor por Convecção Forçada, Natural e Mista numa Cavidade Retangular', M.S. Thesis, Federal University of Itajubá, Itajubá, Brazil, 2006.