Corrigé de l'examen nal du module CFD et Logiciels-2021
Exercice n°1 :
1. La CFD est une science qui à pour objet l'étude de la dynamique des uides par le biais de l'outil numérique.(0.5
pts)
2. a) ϕ est une propriété physique du uide qui dépend du temps et de l'espace et qui véri e les lois de conser-
vation de la physique. ϕ peut représenter la vitesse ou la température par exemple. (0.5 pts)
b) Γ est le coe cient de diusion associé à la propriété physique ϕ (0.5 pts)
c) A est le terme instationnaire, B est le terme convectif, C est le terme de diusion, D est le terme source.
(1.0 pt)
d) L'équation de conservation de la quantité de mouvement (0.5 pts)
Exercice n°2 :
1. Divisez le domaine en 12 volumes de contrôle uniformes avec ∆x = ∆y = 0, 1m en les numérotant de gauche à
droite et de bas en haut. (1.0 pt)
Figure 1 Discrétisation du Domaine géométrique
2. L'intégration de l'équation diérentielle sur un volume de contrôle quelconque donne :
ZZZ ZZZ ZZZ
∂ ∂T ∂ ∂T
(k )dV + (k )dV = QdV (1)
V C ∂x ∂x V C ∂y ∂y VC
L'élément dV = dxdydz , l'intégrale précédente s'écrit donc :
Z ε Z n Z e Z ε Z w Z n Z e Z n Z ε
∂T ∂T
dz dy ∂(k )+ dz dx ∂(k )=Q dx dy dz (2)
0 s w ∂x 0 e s ∂y w s 0
soit : Z e Z n
∂T ∂T
ε∆y ∂(k ) + ε∆x ∂(k ) = Qε∆y∆x (3)
w ∂x s ∂y
Après intégration l'équation précédente s'écrit comme suit :
dT dT dT dT
ε∆y[k ]e − ε∆y[k ]w + ε∆x[k ]n − ε∆x[k ]s = Qε∆y∆x (4)
| {zdx } | {zdx } | {zdx } | {zdx }
1 2 3 4
Appliquons maintenant l'équation 4 pour tous les n÷uds du domaine géométrique :
∂T ∂T
n÷ud 1 en remarquant que −k = q1 au niveau de la frontière Ouest et −k = 0 au niveau de la
∂x ∂y
frontière Sud nous obtenons :
ε k∆y k∆x
∆x [T2 − T1 ] + ε∆yq1 + ε ∆y [T5 − T1 ] + 0 = Qε∆x∆y
Sachant que ∆x = ∆y l'équation précédente s'écrit :
εk[T2 − T1 ] + εk[T5 − T1 ] = Qε∆x∆y − ε∆yq1
Application numérique : εk = 10−2 .390 = 3.9 , Qε∆x∆y = 5.103 .10−2 .10−1 .10−1 = 0.5 et ε∆yq1 =
10−2 .10−1 .200 = 0.2 Nous obtenons donc pour le n÷ud 1
−7.8T1 + 3.9T2 + 3.9T5 = 0.3(1.0 pt) (5)
1
n÷ud 2
3.9[T3 − T2 ] − 3.9[T2 − T1 ] + 3.9[T6 − T2 ] = 0.5 soit :
3.9T1 − 11.7T2 + 3.9T3 + 3.9T6 = 0.5(1.0 pt) (6)
n÷ud 3
3.9[T4 − T3 ] − 3.9[T3 − T2 ] + 3.9[T7 − T3 ] − 2(3.9)TB = 0.5 soit
3.9T2 − 19.5T3 + 3.9T4 + 3.9T7 = −389.5(1.0 pt) (7)
n÷ud 4
3.9T3 − 15.6T4 + 3.9T8 = 0.5(1.0 pt) (8)
n÷ud 5
3.9T1 − 11.7T5 + 3.9T6 + 3.9T9 = 0.3(1.0 pt) (9)
n÷ud 6
3.9T2 + 3.9T5 − 15.6T6 + 3.9T7 + 3.9T10 = 0.5(1.0 pt) (10)
n÷ud 7
3.9T3 + 3.9T6 − 15.6T7 + 3.9T8 + 3.9T11 = 0.5(1.0 pt) (11)
n÷ud 8
3.9T4 + 3.9T7 − 11.7T8 + 3.9T12 = 0.6(1.0 pt) (12)
n÷ud 9
3.9T5 − 15.6T9 + 3.9T10 = −1559.7(1.0 pt) (13)
n÷ud 10
3.9T6 + 3.9T9 − 19.5T10 + 3.9T11 = −1559.5(1.0 pt) (14)
n÷ud 11
3.9T7 + 3.9T10 + 3.9T12 − 11.7T11 = 0.5(1.0 pt) (15)
n÷ud 12
3.9T8 + 3.9T11 − 7.8T12 = 0.6(1.0 pt) (16)
3. Assemblage des équations de tout les n÷uds sous forme matricielle (4.0 pts) :
−7.8 3.9 0 0 3.9 0 0 0 0 0 0 0 T1 0.3
3.9 −11.7 3.9 0 0 3.9 0 0 0 0 0 0 T2 0.5
0
3.9 −19.5 3.9 0 0 3.9 0 0 0 0 0 T3 −389.5
0
0 3.9 −15.6 0 0 0 3.9 0 0 0 0 T4 0.5
3.9
0 0 0 −11.7 3.9 0 0 3.9 0 0 0 T5 0.3
0
3.9 0 0 3.9 −15.6 3.9 0 0 3.9 0 0 T6 = 0.5
0
0 3.9 0 0 3.9 −15.6 3.9 0 0 3.9 0 T7 0.5
0
0 0 3.9 0 0 3.9 −11.7 0 0 0 3.9
T8 0.6
0
0 0 0 3.9 0 0 0 −15.6 3.9 0 0 T9 −1559.7
0
0 0 0 0 3.9 0 0 3.9 −19.5 3.9 0 T10 −1559.5
0 0 0 0 0 0 3.9 0 0 3.9 −11.7 3.9 T11 0.5
0 0 0 0 0 0 0 3.9 0 0 3.9 −7.8 T12 0.6