File:C:\Users\kop's landry\Desktop\SOUTENANCE EES\gear 3.

EES 06/06/2024 15:03:31 Page 1

EES Ver. 9.478: #2363: The Sabotage Rebellion hackers, [Link]

"! A- hypothèses de calcul"

"!#1 les calculs s'effectueront en regime permanant"

"!#2 Le transfert de chaleur se fait de manière unidimensionnelle"

"! #3 Pas de génération interne de chaleur"

"! #4 Coefficients de convection constants"

"! #5 On néglige les résistances thermiques de contact"

"! #6 Pas de changement de phase"

"! #7 Résistance par conduction sur la couche négligeable "

"! #8 Propriétés constantes des materiaux constants"

"! #9 Nous limiterons les pertes de chaleur le long de la conduite de NF DTU 45.2 de 30% (du flux total)
de NF DTU 45.2"

"! #10 Le calcul de l'évaluation de l'epaisseur s'effectuera pour une conduite de longueur 1 m"

"! #11 la nature de l'isoalnt thermique choisi pour notre dimensionnement sera le STYROTEC cpa "

"! #12 pour l'acheminement de l'eau glacée dans les terminaux le type de conduite utilisé sera du
Kryoclim"

"! B- données du calcul"

D_kryo_out = 0,160 "

Diamètre externe de la conduite (m)"
D_kryo_in=0,137 "
Diamètre interne de la conduite (m)"
e_tuyau= D_kryo_out - D_kryo_in "
Épaisseur de la paroi de la conduite (m)"
X_moyenne_conduite = (D_kryo_out-D_kryo_in) {Épaisseur moyenne de la conduite}
L=1 "
Longueur de la conduite (m)"
S_lat_moyenne_kryo = pi * ((D_kryo_out + D_kryo_in) / 2) * L " Surface
latérale moyenne (m^2)"
S_lat_interne_kryo = pi * D_kryo_in * L " Surface
latérale interne (m^2)"
S_lat_externe_kryo = pi * D_kryo_out * L " Surface
latérale externe (m^2)"

"! valeur du flux thermique"

PHI_froid = 194000 " Flux
File:C:\Users\kop's landry\Desktop\SOUTENANCE EES\gear [Link] 06/06/2024 15:03:31 Page 2
EES Ver. 9.478: #2363: The Sabotage Rebellion hackers, [Link]

froid à travers la conduite (W)"

PHI_perdu = alpha* PHI_froid
alpha=0,3 "
pourcentage de réduction des pertes thermiques de NF DTU 45.2 de 30% en moyenne"
PHI_surf_rayonnement = 179 " Densité
du flux solaire incident par rayonnement (W/m^2)"

"! caracteristiques des materiaux"

"! Nature de la conduite : kryoclim"
lambda_kryo = 0,17 "
Conductivité thermique de la conduite (W/m·K)"
"nature de l'isolant : STYROTEC cpa"
lambda_isolant = 0,031 "
Conductivité thermique de l'isolant (W/m·K)"
"! caracteriqtiques externes "
"! Nature du fluide : air"
P_atm = 101,3 " Pression
atmosphérique (kPa)"
T_air = 33 "
Température ambiante de l'air (°C)"
T_air_k = T_air+273 "
Température ambiante de l'air (K)"
sigma = 5,67e-8 "
Constante de Stefan-Boltzmann (W/m^2·K^4)"
epsilon_aluminium = 0,9 "
Émissivité thermique de la couche externe de l'isolant"

"! C- caracteristiques du fluide dans la conduite"

"! Nature du fluide : eau"

Q_dot_veau = 0,0094 " Débit
volumique du fluide frigoporteur (m^3 /s)"
T_eau_in = 7 "
Température d'entrée de l'eau (°C)"
T_eau_out = 12 "
Température de sortie de l'eau (°C)"
T_moyenne_eau = (T_eau_out + T_eau_in) / 2 "
Température moyenne de l'eau (°C)"
x_eau = 0 " Titre en
eau dans la conduite"

"! D- Determination des caracteristiques internes"

"! 1 #masse volumique"

rho_eau = Density(Water; T = T_moyenne_eau; x = x_eau) " Masse

volumique de l'eau (kg/m^3)"

mu_eau = Viscosity(Water; T =T_moyenne_eau; x = x_eau) "

Viscosité dynamique de l'eau (Pa·s) mu"
File:C:\Users\kop's landry\Desktop\SOUTENANCE EES\gear [Link] 06/06/2024 15:03:31 Page 3
EES Ver. 9.478: #2363: The Sabotage Rebellion hackers, [Link]

nu_eau = mu_eau / rho_eau "Viscosité

cinématique de l'eau (m^2/s)"

mu_eau2=10^(-4)*((17,9-0,07377*T_moyenne_eau+0,0003354*(T_moyenne_eau^2))/(1+0,03032*
T_moyenne_eau+8,765*10^(-5)*T_moyenne_eau^2))
" Viscosité dynamique de l'eau (Pa·s) "
lambda_eau = Conductivity(Water; T = T_moyenne_eau; x = x_eau) "
Conductivité thermique de l'eau (W/m·K)"

u_eau = (4 * Q_dot_veau) / (pi * D_kryo_in^2) " Vitesse

de passage de l'eau dans la conduite (m/s)"

"! #2 regime d'eau"

Re_eau = (rho_eau * u_eau * D_kryo_in) / nu_eau "! Nombre de Reynolds de l'eau"

Pr_eau = Prandtl(Water; T = T_moyenne_eau; x = x_eau) "!

Nombre de Prandtl de l'eau"

pH=((13,06+1,387*T_moyenne_eau-0,0037*(T_moyenne_eau^2))/(1+0,12407*T_moyenne_eau+0,005297*
T_moyenne_eau^2))

Nussel_eau = 0,023 * Re_eau^0,8 * Pr_eau^o "!

Nombre de Nusselt de l'eau Corrélation de Dittus-Boelter"
o=0,4 " O=0,4
car T_p>T_moyenne_eau"

"! #3 Coefficient d'echange convectif"

h_conv_eau = (Nussel_eau * lambda_eau) /D_kryo_in "!

Coefficient de convection de l'eau (W/m^2·K)"
{PHI_surf_froid = PHI_froid / S_lat_interne_kryo " ! Flux
surfacique froid (W/m^2)"}

"! #4 Temperatures sur la conduite"

T_p_kryo_in = T_moyenne_eau + (PHI_perdu / (h_conv_eau * S_lat_interne_kryo)) "

Température interne de la paroi (°C)"
T_p_kryo_out = T_p_kryo_in + (PHI_perdu * e_tuyau / (lambda_kryo * S_lat_externe_kryo)) "
Température externe de la paroi (°C)"

"!# 5 Flux de chaleur par conduction a travers la parois de la conduite"

PHI_conduction_kryo = (lambda_kryo * S_lat_moyenne_kryo * (T_p_kryo_out - T_p_kryo_in)) / e_tuyau

"! E- Determination des facteurs mis en jeu lors de l'

apport de chaleur externe"

PHI_surf_rayonnement= epsilon_aluminium * sigma * ((T_parois_alu_rayonnement)^4 - (T_air_k)^4) " La

densité flux par rayonnement (W/m^2)"
File:C:\Users\kop's landry\Desktop\SOUTENANCE EES\gear [Link] 06/06/2024 15:03:31 Page 4
EES Ver. 9.478: #2363: The Sabotage Rebellion hackers, [Link]

"! #1 Temperature sur la paroi en contact avec l'exterieur"

T_parois_alu_rayonnement = T_parois_alu + 273 " Température externe de la paroi (k)"

T_moyenne_1 = (T_air + T_parois_alu) / 2
T_moyenne_2 = T_moyenne_1 + 273

"! #2 le coeficient de transfert radiatif"

h_ray_air = 4 * epsilon_aluminium * sigma * (T_parois_alu_rayonnement)^3

g = 9,81 " La
connstante gravitationnelle"

D_ext =0,1608

rho_air = Density(Air; T = T_moyenne_1; P = P_atm) " Masse

volumique de l'air (kg/m^3)"

beta = 1 / T_moyenne_2 "

Coefficient de dilatation thermique (K^-1) "

DELTAT = T_parois_alu - T_air

mu_air = 10^(-5) * (0,0046 * T_moyenne_1 + 1,7176) "

Viscosité dynamique de l'air (kg/m^3)"

lambda_air=Conductivity(Air;T=T_moyenne_1) "
Conductivité thermique de l'air (W/m·K)"

Pr_air = Prandtl(Air_ha; T= T_moyenne_2; P = P_atm) " Nombre

de Prandtl de l'air"
Gr_air = (rho_air^2 * g * beta * DELTAT * D_ext^3) / mu_air^3 " Nombre
de Grashof de l'air"

Ra = Gr_air * Pr_air " Nombre

de Rayleigh"
C = 0,1
n =0,33
nussel_air=C* Ra^n "!
Nombre de Nusselt de l'air"

"! #3 Coefficient d'echage convectif"

h_conv_air = (Nussel_air * lambda_air) / D_ext "
Coefficient de convection de l'air (W/m^2·K)"

T_isolant_out = T_parois_alu

phi_surf_convectif = h_conv_air * (T_parois_alu - T_air)

"! le coefficient d’échange superficiel"

h_ray_con = h_ray_air + h_conv_air " le

coefficient d’échange superficiel externe global (W/m^2·K)"
File:C:\Users\kop's landry\Desktop\SOUTENANCE EES\gear [Link] 06/06/2024 15:03:31 Page 5
EES Ver. 9.478: #2363: The Sabotage Rebellion hackers, [Link]

"! G- Dimensions de l'isolation"

"! #1 equations generales de transfert thermique"

"PHI_perdu = (pi*(T_moyenne_eau - T_air)) / (R_convective_eau+R_conduction_conduite +R_c
onduction_isolant + R_convective_externe )"

"où"

R_convective_externe =1 /(D_ext *h_ray_con) "! La résistance superficielle au contact de l’air"

R_convective_eau =1/(D_kryo_in*h_conv_eau) "! La résistance superficielle au contact du fluide"
R_conduction_isolant = ln((D_ext) / (D_kryo_out)) / (2*lambda_isolant) "! La
résistance de la couche isolante"
R_conduction_conduite=e_tuyau/(2*lambda_kryo) "! La résistance de la paroi métalique de la conduite"

"! #2 valeur de l'epaisseur de l'isolation"

X_isolant= (pi*(T_moyenne_eau - T_air)) / (PHI_perdu*(R_convective_eau+ R_conduction_conduite+

R_conduction_isolant+ R_convective_externe))+D_kryo_out/2

"! 3 la valeur du diametre "

diamètre_isolant =D_kryo_out+2*X_isolant
"! 4 la valeur de la surface radiale de l'isolant "

S_isolant=(pi*((diamètre_isolant)^2-(D_kryo_out)^2))/4

{phi_perdu = ((pi)*(T_moyenne_eau - T_air))/((R_convective_eau)+(R_conduction_conduite)+(ln(D_ext/

D_isolant_in)/(2*lambda_isolant))+(R_convective_externe))}

A- hypothèses de calcul

#1 les calculs s'effectueront en regime permanant

#2 Le transfert de chaleur se fait de manière unidimensionnelle

#3 Pas de génération interne de chaleur

#4 Coefficients de convection constants

#5 On néglige les résistances thermiques de contact

#6 Pas de changement de phase

#7 Résistance par conduction sur la couche négligeable

#8 Propriétés constantes des materiaux constants

#9 Nous limiterons les pertes de chaleur le long de la conduite de NF DTU 45.2 de 30% (du flux total) de NF DTU 45.2

#10 Le calcul de l'évaluation de l'epaisseur s'effectuera pour une conduite de longueur 1 m

#11 la nature de l'isoalnt thermique choisi pour notre dimensionnement sera le STYROTEC cpa

#12 pour l'acheminement de l'eau glacée dans les terminaux le type de conduite utilisé sera du Kryoclim
 File:C:\Users\kop's landry\Desktop\SOUTENANCE EES\gear [Link] 06/06/2024 15:03:31 Page 6
EES Ver. 9.478: #2363: The Sabotage Rebellion hackers, [Link]

B- données du calcul

D kryo;out = 0,16 Diamètre externe de la conduite (m)

D kryo;in = 0,137 Diamètre interne de la conduite (m)

e tuyau = D kryo;out – D kryo;in Épaisseur de la paroi de la conduite (m)

X moyenne;conduite = D kryo;out – D kryo;in

L = 1 Longueur de la conduite (m)

D kryo;out + D kryo;in
S lat;moyenne;kryo = p · · L Surface latérale moyenne (m2)
2

S lat;interne;kryo = p · D kryo;in · L Surface latérale interne (m2)

S lat;externe;kryo = p · D kryo;out · L Surface latérale externe (m2)

valeur du flux thermique

F froid = 194000 Flux froid à travers la conduite (W)

F perdu = a · F froid

a = 0,3 pourcentage de réduction des pertes thermiques de NF DTU 45.2 de 30% en moyenne

F surf;rayonnement = 179 Densité du flux solaire incident par rayonnement (W/m2)

caracteristiques des materiaux

Nature de la conduite : kryoclim

l kryo = 0,17 Conductivité thermique de la conduite (W/m·K)

nature de l'isolant : STYROTEC cpa

l isolant = 0,031 Conductivité thermique de l'isolant (W/m·K)

caracteriqtiques externes

Nature du fluide : air

P atm = 101,3 Pression atmosphérique (kPa)

T air = 33 Température ambiante de l'air (°C)

T air;k = T air + 273 Température ambiante de l'air (K)

s = 5,67 x 10 Constante de Stefan-Boltzmann (W/m2·K4)

–8

e aluminium = 0,9 Émissivité thermique de la couche externe de l'isolant

C- caracteristiques du fluide dans la conduite

Nature du fluide : eau

Q veau = 0,0094 Débit volumique du fluide frigoporteur (m3 /s)

File:C:\Users\kop's landry\Desktop\SOUTENANCE EES\gear [Link] 06/06/2024 15:03:31 Page 7
EES Ver. 9.478: #2363: The Sabotage Rebellion hackers, [Link]

T eau;in = 7 Température d'entrée de l'eau (°C)

T eau;out = 12 Température de sortie de l'eau (°C)

T eau;out + T eau;in
T moyenne;eau = Température moyenne de l'eau (°C)
2

x eau = 0 Titre en eau dans la conduite

D- Determination des caracteristiques internes

1 #masse volumique

r eau = r water ; T = T moyenne;eau ; x = x eau Masse volumique de l'eau (kg/m3)

m eau = Visc water ; T = T moyenne;eau ; x = x eau Viscosité dynamique de l'eau (Pa·s) mu

m eau
n eau = Viscosité cinématique de l'eau (m2/s)
r eau

17,9 – 0,07377 · T moyenne;eau + 0,0003354 · T moyenne;eau

2
m eau2 = 10 ·
– 4
Viscosité dynamique de l'eau (Pa·s)
1 + 0,03032 · T moyenne;eau + 8,765 · 10 · T moyenne;eau
– 5 2

l eau = k water ; T = T moyenne;eau ; x = x eau Conductivité thermique de l'eau (W/m·K)

4 · Q veau
u eau = Vitesse de passage de l'eau dans la conduite (m/s)
p · D kryo;in
2

#2 regime d'eau

r eau · u eau · D kryo;in

Re eau = Nombre de Reynolds de l'eau
n eau

Pr eau = Pr water ; T = T moyenne;eau ; x = x eau Nombre de Prandtl de l'eau

13,06 + 1,387 · T moyenne;eau – 0,0037 · T moyenne;eau

2
pH =
1 + 0,12407 · T moyenne;eau + 0,005297 · T moyenne;eau
2

Nussel eau = 0,023 · Re eau · Pr eau Nombre de Nusselt de l'eau Corrélation de Dittus-Boelter
0,8 o

o = 0,4 O=0,4 car Tp>T,moyenne,eau

#3 Coefficient d'echange convectif

Nussel eau · l eau

h conv;eau = Coefficient de convection de l'eau (W/m2·K)
D kryo;in

#4 Temperatures sur la conduite

F perdu
T p;kryo;in = T moyenne;eau + Température interne de la paroi (°C)
h conv;eau · S lat;interne;kryo

e tuyau
T p;kryo;out = T p;kryo;in + F perdu · Température externe de la paroi (°C)
l kryo · S lat;externe;kryo

# 5 Flux de chaleur par conduction a travers la parois de la conduite

File:C:\Users\kop's landry\Desktop\SOUTENANCE EES\gear [Link] 06/06/2024 15:03:31 Page 8
EES Ver. 9.478: #2363: The Sabotage Rebellion hackers, [Link]

l kryo · S lat;moyenne;kryo · T p;kryo;out – T p;kryo;in

F conduction;kryo =
e tuyau

E- Determination des facteurs mis en jeu lors de l'apport de chaleur externe

F surf;rayonnement = e aluminium · s · T parois;alu;rayonnement – T air;k La densité flux par rayonnement (W/m2)

4 4

#1 Temperature sur la paroi en contact avec l'exterieur

T parois;alu;rayonnement = T parois;alu + 273 Température externe de la paroi (k)

T air + T parois;alu
T moyenne;1 =
2

T moyenne;2 = T moyenne;1 + 273

#2 le coeficient de transfert radiatif

h ray;air = 4 · e aluminium · s · T parois;alu;rayonnement

3

g = 9,81 La connstante gravitationnelle

D ext = 0,1608

r air = r Air ; T = T moyenne;1 ; P = P atm Masse volumique de l'air (kg/m3)

1
b = Coefficient de dilatation thermique (K-1)
T moyenne;2

DT = T parois;alu – T air

m air = 10 · 0,0046 · T moyenne;1 + 1,7176 Viscosité dynamique de l'air (kg/m3)

– 5

l air = k Air ; T = T moyenne;1 Conductivité thermique de l'air (W/m·K)

Pr air = Pr Air ha ; T = T moyenne;2 ; P = P atm Nombre de Prandtl de l'air

r air · g · b · DT · D ext
2 3
Gr air = Nombre de Grashof de l'air
m air
3

Ra = Gr air · Pr air Nombre de Rayleigh

C = 0,1

n = 0,33

nussel air = C · Ra Nombre de Nusselt de l'air

n

#3 Coefficient d'echage convectif

nussel air · l air

h conv;air = Coefficient de convection de l'air (W/m2·K)
D ext

T isolant;out = T parois;alu

f surf;convectif = h conv;air · T parois;alu – T air

le coefficient d échange superficiel

h ray;con = h ray;air + h conv;air le coefficient d échange superficiel externe global (W/m2·K)

 File:C:\Users\kop's landry\Desktop\SOUTENANCE EES\gear [Link] 06/06/2024 15:03:31 Page 9
EES Ver. 9.478: #2363: The Sabotage Rebellion hackers, [Link]

G- Dimensions de l'isolation

#1 equations generales de transfert thermique

PHIperdu = (pi*(Tmoyenne,eau - Tair)) / (Rconvective,eau+R,conduction,conduite +Rconduction,isolant + Rconvective,externe )

où

1
R convective;externe = La résistance superficielle au contact de l air
D ext · h ray;con

1
R convective;eau = La résistance superficielle au contact du fluide
D kryo;in · h conv;eau

D ext
ln
D kryo;out La résistance de la couche isolante
R conduction;isolant =
2 · l isolant

e tuyau
R conduction;conduite = La résistance de la paroi métalique de la conduite
2 · l kryo

#2 valeur de l'epaisseur de l'isolation

p · T moyenne;eau – T air D kryo;out

X isolant = +
F perdu · R convective;eau + R conduction;conduite + R conduction;isolant + R convective;externe 2

3 la valeur du diametre

diamètre isolant = D kryo;out + 2 · X isolant

4 la valeur de la surface radiale de l'isolant

p · diamètre isolant – D kryo;out

2 2
S isolant =
4

SOLUTION
Unit Settings: SI C kPa kJ mass deg
a = 0,3 b = 0,003131 [C ]
-1

C = 0,1 DT = 26,86 [C]

diamètreisolant = 0,3073 [m] Dext = 0,1608 [m]
Dkryo,in = 0,137 [m] Dkryo,out = 0,16 [m]
ealuminium = 0,9 etuyau = 0,023
g = 9,81 [N/kg] Grair = 5,807E+11
hconv,air = 114,3 [W/m2*K] hconv,eau = 426220 [W/m2*K]
hray,air = 7,528 [W/m2*K] hray,con = 121,8 [W/m2*K]
L = 1 [m] lair = 0,02709 [W/m*K]
leau = 0,5664 [W/m*K] lisolant = 0,031 [W/m*K]
lkryo = 0,17 [W/m*K] mair = 0,00001931 [m /s]
2

meau = 0,001326 [kg/m-s] meau2 = 0,001329 [kg/m*s]

n = 0,33 nusselair = 678,5
Nusseleau = 103088 neau = 0,000001326 [m /s]
2

o = 0,4 pH = 9,75
Fconduction,kryo = 54017 [W/m] Ffroid = 194000 [W/m]
Fperdu = 58200 [W/m] fsurf,convectif = 3070 [W/m ]
2

Fsurf,rayonnement = 179 [W] Prair = 0,7027

Preau = 9,804 Patm = 101,3
File:C:\Users\kop's landry\Desktop\SOUTENANCE EES\gear [Link] 06/06/2024 15:03:32 Page 10
EES Ver. 9.478: #2363: The Sabotage Rebellion hackers, [Link]

Qveau = 0,0094 [m3/s] Ra = 4,081E+11

Reeau = 6,587E+07 rair = 1,104 [kg/m ]
3

reau = 999,7 [kg/m ] Rconduction,conduite = 0,06765 [m2*K/W]

3

Rconduction,isolant = 0,08044 [m2*K/W] Rconvective,eau = 0,00001713 [m2*K/W]

Rconvective,externe = 0,05104 [m2*K/W] s = 5,670E-08
Sisolant = 0,05404 [m2] Slat,externe,kryo = 0,5027 [m2]
Slat,interne,kryo = 0,4304 [m2] Slat,moyenne,kryo = 0,4665 [m2]
Tair = 33 [C] Tair,k = 306 [K]
Teau,in = 7 [C] Teau,out = 12 [C]
Tisolant,out = 59,86 [C] Tmoyenne,1 = 46,43 [C]
Tmoyenne,2 = 319,4 [C] Tmoyenne,eau = 9,5 [C]
Tparois,alu = 59,86 [C] Tparois,alu,rayonnement = 332,9 [K]
Tp,kryo,in = 9,817 [C] Tp,kryo,out = 15675 [C]
ueau = 0,6377 [m/s] xeau = 0
Xisolant = 0,07363 [m] Xmoyenne,conduite = 0,023 [m]

No unit problems were detected.