Annexe 1
La variation de la viscosité des gaz avec la température.
La viscosité des corps purs
Annexe 2
Abaques de calcul de facteur de compressibilité
L’équation non linéaire suivante doit être résolue pour obtenir la masse volumique
réduite Y ; pour cela, on utilise généralement la méthode de NEWTON-RAPHSON.
2 𝐘+𝐘 2 +𝐘 3 −𝐘 4
𝐅(𝐘) = −0.06125𝐏𝐏𝐑 ∗ 𝝈 ∗ 𝐞[−1.2(1−𝝈) ] + (1−𝐘)3
– (14.76 ∗ 𝝈 − 9.76 ∗ 𝝈2 +
3 )𝐘 (2.18+2.82∗𝝈)
4.58 ∗ 𝝈3 )𝐘 2 + (90.7 ∗ 𝝈 − 242.2 ∗ 𝝈2 + 42.4 ∗ 𝝈 =0
Diagramme de moody pour le calcul du coefficient des pertes de charge
Méthode de Méthode de Méthode de Méthode de
Dranchuk- Dranchuk- Hankinson- Hankinson-
Abou Kassem Purvis- Thomas- Thomas-
Robinson Philips Philips
0.4≤ 𝑷𝑷𝑹 ≤5 5≤ 𝑷𝑷𝑹 ≤15
𝑨1 0.3265 0.31506237 0.001290236 0.0014507882
𝑨2 -1.0700 -1.0467099 0.38193005 0.37922269
𝑨3 -0.5339 -0.57832729 0.022199287 0.024181399
𝑨4 0.01569 0.53530771 0.12215481 0.11812287
𝑨5 -0.05165 -0.61232032 -0.015674794 0.037905663
𝑨6 0.5475 -0.10488813 0.027271364 0.19845016
𝑨7 -0.7361 0.68157001 0.023834219 0.048911693
𝑨8 0.1844 0.68446549 0.43617780 0.0631425417
𝑨9 0.1056 ****** ****** ******
𝑨10 0.6134 ****** ****** ******
𝑨11 0.7210 ****** ****** ******
Les valeurs numériques des coefficients des différentes méthodes des corrélations
empiriques.
Propriétés des gaz d’hydrocarbure
Le calcul thermodynamique
Les constantes pour le calcul de coefficient de compressibilité :
Xi Tc Pc ai bi xi*ai^0,5 xi*bi
C1 0,8488 343,04 666 7072,4521 0,4788436 71,38227 0,406442
4 9 8 43 52
C2 0,0737 549,75 707 21660,835 0,7228842 10,84688 0,053276
6 8 64 57
C3 0,0168 665,67 617 40044,397 1,0029883 3,361864 0,016850
6 6 8 18 2
nC4 0,0035 765,32 548,8 63807,085 1,2964309 0,884102 0,004537
4 7 7 26 51
iC4 0,0022 734,67 527,9 59889,237 1,2937751 0,538390 0,002846
4 6 11 31
nC5 0,0008 845,37 488,1 91997,908 1,6101123 0,242649 0,001288
8 26 09
iC5 0,0008 828,79 490,4 87143,117 1,5711341 0,236160 0,001256
2 9 11 91
nC6 0,001 913,37 436,9 124712,28 1,9435009 0,353146 0,001943
4 8 9 27 5
N2 0,0482 226,98 493 3402,4499 0,4280148 2,811531 0,020630
1 5 21 32
CO2 0,0042 547,54 1071 14155,454 0,4752762 0,499702 0,001996
2 7 12 16
am 8309,545099
bm 0,511068078
A 0,163385247
B 0,055924897
La solution de l’équation de Redlich Kong
Calcul de la chaleur spécifique
M(kg/kmol) Cp(J/kg,K) Mi*XI cp*XI
16,43 2205,44 13,945784 1871,97747
30,07 1707,75 2,216159 125,861175
44,097 1626,45 0,7408296 27,32436
58,124 1618,99 0,203434 5,666465
58,124 1626,45 0,1278728 3,57819
72,151 1625,77 0,0577208 1,300616
72,151 1609,74 0,0577208 1,287792
86,178 1617,4 0,086178 1,6174
28,013 1039,7 1,3502266 50,11354
44,01 833,55 0,184842 3,50091
𝑪𝒑 = 2100 𝑴𝒎𝒐𝒚 = 19
Calcul de la viscosité
𝝁 𝝁 *yi*Mi^0,5 yi*Mi^0,5
0,0000134 4,6103E-05 3,44052052
0,0000089 3,5969E-06 0,4041422
0,0000075 8,3671E-07 0,11156136
0,0000073 1,9479E-07 0,02668368
0,0000071 1,1909E-07 0,0167726
0,0000066 4,4849E-08 0,00679534
0,0000066 4,4849E-08 0,00679534
0,0000063 5,8484E-08 0,00928321
0,0000176 4,4899E-06 0,25510963
0,0000148 4,1237E-07 0,02786281
𝝁 1.29. 105
Calcul du nombre de Reynolds et de coefficient des pertes de charge
Re 13788778,27
lamda 0,032803323
Annexe 3
Les Organigrammes
Notation des symboles des organigrammes
𝑆𝑦𝑚𝑏𝑜𝑙𝑒 Définition
𝑀 Débit massique
𝑑 Diamètre de la conduite
𝑡𝑠 Température du sol
𝑇0 Température initiale
𝑃0 Pression initiale
𝐿 Longueur de la conduite
𝐶𝑝 Chaleur spécifique
𝛾 Exposant polytropique
𝑘 Coefficient du transfert de chaleur globale
𝜆 Coefficient de pertes de charge
𝑍𝑚0 facteur de compressibilité moyen initial
𝑝𝑖 Pression initiale
𝑛 Nombre de points de discrétisation
𝜇 Viscosité dynamique
𝑟𝑔 Constante de gaz
𝑎1 Coefficient de Chokove
𝑑𝑥 Pas d’espace
𝑠 Section de la conduite
𝑖 Point de discrétisation
𝑇(𝑖) Température à chaque pas d’espace
𝑡𝑚 Température moyenne
𝑃(𝑖) Pression à chaque pas d’espace
𝑝𝑚 Pression moyenne
𝜌(𝑖) Masse volumique à chaque pas d’espace
𝑍𝑚 Facteur de compressibilité moyen
𝑢(𝑖) Vitesse d’écoulement de fluide à chaque pas d’espace
𝑎(𝑖) Vitesse de propagation d’onde à chaque pas d’espace
𝑑𝑗 Coefficient de joule Thomson
𝑡𝑚𝑎𝑥 Température maximale
𝑡𝑠𝑟 Température du sol réduite
𝑑𝑡 Pas du temps
𝑑𝑥 Pas d’espace
𝑎𝑚 Vitesse de propagation d’onde moyenne
Les équations utilisées pour le calcul des constantes thermodynamiques :
4,58. 1011 . 𝑃𝑚
𝑐𝑝 = (48,13 + ) . 𝑇𝑚0.665
𝑇𝑚5
𝜇 = (0.316. 𝑇𝑚 + 0.0175. 𝑃𝑚 + 1.63). 10−6
5650
𝐷𝑗 = √224 − 𝑃𝑚
𝑇𝑚2.13
∆1.5 𝑃𝑚
𝑍 = 1 − 2. 106 .
𝑇𝑚1.5
