Test de validation

 

Test de validation

 

 

Après avoir écrit et inséré toutes les subroutines, on s'est intéressé à la validation du code qui est une étape indispensable avant toute utilisation. Notre test de validation consiste à comparer nos résultats avec les résultats obtenus par De Vahl Davis et Le Quéré .

La configuration du test est la suivante:

Nous considérons une cavité carrée( H= b=0.05 (m)  )dans laquelle on impose un gradient de température horizontale et telle que ρ0=1.2928 (kg/m3 )     μ=1.719.10-5 (kg/m.s)   k=1.83.10-5(m2/s) . ∆=Th-Tc=278.4 -273 (K). Cela donne Pr=0.72 et Ra=105 .

Le tableau suivant récapitule les résultats de De Vahl et le Quéré.

 

 

 

  103 104 105
Umax(1/2,z) 3.634 16.18 34.8
z 0.81 0.82 0.85
vmax(x,1/2) 3.7 19.5 68.6
x 0.18 0.12 0.066

Umax(1/2,z)= correspont à la vitesse horizontale maximale sur une droite verticale à mi cavité.

Vmax(x,1/2)= correspont à la vitesse verticale maximale sur une droite horizontale à mi cavité.

Ps:Les vitesse donnés sont adimensionnées : division sur k/H .

Résultats obtenus ( Ra= 105):

Le graphe suivant présente les résultats obtenus avec notre code:

 

 

 

 

 

 

 

 

 

 

 

Fig 1 : Profile de U à mi cavité (droite verticale)                                                                                                                         Fig 2: Profil de V à mi cavité(droite horizontale)

 

  Résultats du code Résultats De Vahl / Le Queré
Umax 34.97 34.8
z 0.86 0.85
Vmax 68.30 68.6
x 0.07 0.066

Donc : on peut affirmer que le code donne de bon résultats avec une erreur < 6 % .