Etude du modèle mécaniste dit à quatre (4) flux proposée par Yeoh et al (2008) pour la prédiction du transfert de chaleur dans la région d'ébullition sous saturée.

4.4. Organigramme de programme :

Entré de données géométriques du
canal et les conditions initiales

Lecture des données

Calcul des propriétés thermo-physiques du
fluide à l'état de saturation pour la pression
P.

j = 1

l'état de saturation : P, _TSAT, VF, VG, CpF, CpG, HF, HG, KL, KG,

Impression des propriétés thermo-physiques du fluide à

ìL, ìG

A

G (J) = 100 + (j-1) * 100

Calcul des propriétés thermo-hydrauliques à la
vitesse G (J) : ReL, PrL, ReyL et hFC.

1

N = 1

B

N = N + 1

ÖW (N) =100 + (N-1)×100

K = K+1

C

K = 1

Z (K) = Z₀ + (K -- 1) × n

n : Le pas (n = 0.01) [m]

Calcul de _ÖONB (J, K)

0_w(N) > ⁰0NB(J,K)

Oui Non

Existence du régime diphasique

Transfert en convection forcée
simple phase liquide

Calcul les cotes des differences regions :

z_nb(J, N), zFDB(J, N) et z(J, N)

3

2

Calcul des propriétés thermo-physiques du
l'eau et sa vapeur à la cote Z(K) du canal pour
la pression P.
HFZ, HGZ, TFZ, VFZ, VGZ, SIGMAZ, visfz,

Calcul des propriétés thermo-physiques du
l'eau à la cote Z(K) du canal pour la
pression P.
HFZ, TFZ, VGZ, SIGMAZ, visfz ,....

4

5

Region I
(monophasique)

Même procedure
que simple

h

Calcul HFC (J, N, K)
(Par la correlation de
Bittus Boeter)

Z(K) > Z_nb(J, N )

Oui Non

Region II (diphasique)

Tw(N ) *

~~, ~~, Rf , f, ...

Tw(J, N ,K) *

(I) * (I) * (I)^* + (I) * (DW EI) *

fc e tc tcEs1

(DW (J, N ) * = (DW (N )

Oui

Non

6

Tw(N ) = Tw(J, N ) *

5 6

Impression des résultats:

z_nbE(J, N), zFDB (J, N) et z_sc(J, N)

TW (J, N, K), h_Yeoh (J, N, K).

Impression des résultats:

z_nb(J, N), zFDB(J, N) et z_sc(J, N)

TW (J, N, K), Hfc (J, N, K).

Impression des résultats:
TW (J, N, K), HFC (J, N, K).

C

Non

Z() Zsc. (J, N)

Oui

B

N = 6

Non

A

Oui

Non

J = 12

Oui

Fin

Figure 4-4 : Organigramme du programme