bn3050

Réf.
: BN3050 V2
Thermohydraulique des
Date de publication :
10 mars 2018 réacteurs à eau sous
pression
Cet article est issu de : Énergies | Génie nucléaire
par Christophe HERER
Mots-clés Résumé La chaleur produite dans le cœur d’un réacteur nucléaire est transférée au fluide
sûreté | écoulement du circuit primaire. Dans les générateurs de vapeur, le circuit primaire échange cette
diphasique |
thermohydraulique | taux de chaleur en vaporisant l'eau du circuit secondaire. Il est nécessaire d’évacuer une certaine
vide quantité d’énergie à l’extérieur. Ces exemples relèvent tous de phénomènes analysés par
la thermohydraulique, à savoir l’hydraulique avec transfert de chaleur et possible
changement de phase (condensation ou ébullition). Cet article traite des domaines
d’application de la thermohydraulique. Les notions nécessaires pour le dimensionnement
et pour les analyses de sûreté des réacteurs à eau sous pression sont explicitées.
Keywords Abstract Heat generated in the core of a nuclear power plant is transferred to the water.
safety | two-phase flows | In the steam generators, this heat is used to vaporize the water in the secondary circuit.
thermal hydraulics | void
fraction Heat resulting from the second law of thermodynamics, or from resulting power, has to be
dumped outside. These examples are typical of thermalhydraulics applications, hydraulics
with heat transfer and possible phase change (boiling and condensation). This article
presents fundamentals of thermalhydraulics. Basics for design and for safety analyses of
pressurized water reactors are presented.
Pour toute question :

Service Relation clientèle
Techniques de l’Ingénieur
Immeuble Pleyad 1 Document téléchargé le : 12/05/2020
39, boulevard Ornano
93288 Saint-Denis Cedex Pour le compte : 7200104319 - altran technologies france // matthieu MARTINELLI // 84.14.209.116
Par mail :
[email protected]
Par téléphone :
00 33 (0)1 53 35 20 20 © Techniques de l'Ingénieur | tous droits réservés
Ce document a ete delivre pour le compte de 7200104319 - altran technologies france // matthieu MARTINELLI // 84.14.209.116
Thermohydraulique des réacteurs

à eau sous pression
Ingénieur chargé d’expertises
Institut de Radioprotection et de Sûreté Nucléaire, Fontenay-aux-Roses, France
Avec la participation de Tony Glantz (IRSN) pour la rédaction des sous-sections 3.4 et 5.2.
1. Notions élémentaires .............................................................................. BN 3 050v2 - 2

1.1 Diamètre hydraulique .............................................................................. — 2
1.2 Valeurs moyennes et maximales des paramètres ................................ — 2
Parution : mars 2018 - Ce document a ete delivre pour le compte de 7200104319 - altran technologies france // matthieu MARTINELLI // 84.14.209.116
2. Caractérisation de l’état du fluide .......................................................... — 3

2.1 État énergétique ....................................................................................... — 3
2.2 Grandeurs diphasiques ........................................................................... — 5
2.3 Régimes d’écoulements .......................................................................... — 7
3. Modélisations........................................................................................... — 9
3.1 Démarche proposée et hypothèses associées ...................................... — 9
3.2 Modèle homogène à l’équilibre (HEM) .................................................. — 9
3.3 Modèle à dérive de flux ........................................................................... — 9
3.4 Modèles à deux fluides............................................................................ — 11
4. Corrélations spécifiques.......................................................................... — 11
4.1 Grandeurs thermohydrauliques et incertitudes .................................... — 12
4.2 Évaluation du taux de vide...................................................................... — 12
4.3 Évaluation des pertes de charge............................................................. — 14
4.4 Évaluation du débit critique .................................................................... — 23
4.5 Transfert thermique ................................................................................. — 26
5. Sûreté des REP......................................................................................... — 33
5.1 Crise d’ébullition ...................................................................................... — 33
5.2 Accident de perte de réfrigérant du circuit primaire (APRP) ................ — 38
6. Conclusion ................................................................................................ — 42
7. Glossaire ................................................................................................... — 42
8. Sigles, notations et symboles ................................................................ — 42
Pour en savoir plus .................................................................................. Doc. BN 3 050v2
L a thermohydraulique est l’étude du transfert de la chaleur couplé à

l’hydraulique des écoulements, avec la prise en compte d’une masse volu-
mique variable résultant par exemple de l’échauffement ou du refroidissement
du fluide. Cependant, l’aspect le plus complexe demeure le changement de
phase et le traitement d’écoulements diphasiques eau-vapeur. Dans un réac-
teur à eau pressurisée, l’eau du circuit primaire a un double rôle de caloporteur
et de modérateur. Outre les enjeux de la maîtrise des transferts d’énergie et
donc des températures dans les systèmes et composants d’un réacteur à eau
sous pression, sa température ou la quantité de vapeur présente a un impact
primordial sur la réactivité et donc la neutronique du cœur.
La thermohydraulique s’intéresse à l’état du fluide, à son mode de circulation
et aux échanges thermiques dans le réacteur, dans les situations normales
Copyright © – Techniques de l’Ingénieur – Tous droits réservés

BN 3 050v2 – 1
tiwekacontentpdf_bn3050 v2 Ce document a ete delivre pour le compte de 7200104319 - altran technologies france // matthieu MARTINELLI // 84.14.209.116
THERMOHYDRAULIQUE DES RÉACTEURS À EAU SOUS PRESSION ____________________________________________________________________________
d’exploitation, d’incidents et d’accidents y compris les plus graves afin que les
critères de sûreté se trouvent respectés en toutes circonstances.
Les méthodes et modèles présentés dans cet article constituent les principes de
base pour permettre de mieux appréhender les problèmes liés à la mécanique des
fluides avec échange de chaleur ou en écoulement à deux phases (liquide et
vapeur). On adoptera une approche empirique, les modèles théoriques et équa-
tions sont du domaine des ouvrages spécialisés. Ces modèles empiriques sont
principalement issus de l’hydraulique en charge et de la thermique dans les
conduites de section circulaire. Cet exposé n’a pas vocation à couvrir l’ensemble
des connaissances, notamment les phénomènes liés à l’analyse des accidents
graves ne sont pas abordés. On trouvera dans cet article une introduction aux
phénomènes physiques, à leur interprétation et à leur modélisation.
1. Notions élémentaires
1.1 Diamètre hydraulique
Une des spécificités de l’industrie nucléaire et du réacteur à eau

Sous-canal
sous pression (REP) réside, en particulier, dans des géométries typique
complexes dans le cœur, la cuve et l’ensemble des composants
du réacteur où l’eau s’écoule dans un milieu confiné mais ouvert
dans lequel les échanges transverses sont importants.
Sous-canal
Une approche simplifiée communément utilisée dans l’hydrau- tybe guide
lique (par rapport à la mécanique des fluides) considère des écou-
lements unidirectionnels dans des conduites circulaires. Afin
d’utiliser pour toute géométrie d’écoulement les résultats de
l’hydraulique établis en conduite circulaire, on utilise la formula-
tion de Chézy qui détermine un diamètre hydraulique Dh équiva-
lent prenant en compte la surface de passage du fluide A et le
périmètre mouillé . Chézy propose la relation simple :
(1) pas
Si le concept du diamètre hydraulique est acceptable, son appli-

cation, et notamment avec la relation (1), introduit des erreurs et Figure 1 – Surface élémentaire dite « sous-canal »
des incertitudes supplémentaires dans les modèles utilisés. Ces
erreurs et incertitudes peuvent être réduites par la détermination Exemple 1
d’une relation pour le diamètre hydraulique plus spécifique et Soient D = 9,5 mm et Dtg = 12,45 mm avec un pas de 12,6 mm.
adaptée à une géométrie donnée.
Surface de passage
Dans le cœur d’un REP, en fonctionnement nominal, la vitesse
– sous canal « typique » : A = 87,88 mm2
du fluide est à direction privilégiée verticale ascendante le long
– sous canal « tube-guide » : A = 75,16 mm2
des assemblages de combustible. Le combustible est constitué de
barreaux cylindriques à section circulaire (de 9,5 mm de diamètre Périmètre mouillé
le plus souvent et exclusivement sur le parc français). Usuelle- – sous canal « typique » :
ment, on définit une surface caractéristique, dénommée « sous- – sous canal « tube-guide » :
canal », qui soit suffisamment petite pour représenter la physique
mais compatible avec les capacités de modélisation et de simula- Diamètre hydraulique
tion numérique. Une telle surface, (illustrée à la figure 1 dans le – sous canal « typique » : Dh = 11,78 mm
cas d’assemblages de combustible à pas carré caractéristique des – sous canal « tube-guide » : Dh = 9,35 mm
réacteurs à eau sous pression à l’exception des réacteurs de Pour un assemblage entier (carré de 214 mm de côté avec
conception soviétique et des réacteurs à eau lourde) représente la 264 crayons, 24 tubes-guide et un tube d’instrumentation de même
surface comprise entre 4 crayons (sous-canal dit « typique ») ou diamètre), la surface de passage est de 0,02404 m2, le périmètre
3 crayons et le tube-guide (sous-canal dit « tube-guide »). mouillé vaut 8,86 mètres ; ce qui donne un diamètre hydraulique
La section de passage A et le périmètre mouillé dépendent équivalent de 10,86 mm.
du type de sous-canal considéré :
– pour un sous-canal « typique » : 1.2 Valeurs moyennes et maximales
des paramètres
– pour un sous-canal « tube-guide » : Les caractéristiques de fonctionnement relatives au circuit pri-
maire définissent les paramètres moyens. À partir des valeurs
moyennes, des valeurs maximales enveloppes sont estimées à
l’aide de facteurs multiplicatifs établis en fonction du type de réac-

BN 3 050v2 – 2
_____________________________________________________________________________ THERMOHYDRAULIQUE DES RÉACTEURS À EAU SOUS PRESSION
teur. Les valeurs moyennes (cf. exemple 2) ne prennent pas en La densité de puissance, ou flux de chaleur moyen φ, est la
compte les hétérogénéités dans le cœur. D’une part la puissance puissance du crayon rapportée à l’unité de surface (produit de la
n’est pas distribuée également dans le cœur, d’autre part l’inser- hauteur active et du périmètre du crayon) soit :
tion des barres de commande modifie ces distributions. Radiale-
ment, on définit le facteur d’élévation d’enthalpie , rapport de
l’élévation d’enthalpie du canal considéré à l’élévation d’enthalpie
moyenne. Axialement, la puissance est distribuée suivant des
formes variables en fonction de l’épuisement du combustible et
de la position des grappes de régulation. Une forme usuellement 2. Caractérisation de l’état
adoptée, permettant d’englober de nombreuses situations, suit
une loi en cosinus. Cependant, pour certains transitoires de fonc- du fluide
tionnement, d’autres distributions axiales de puissance, non
symétriques (avec un axial offset non nul) sont considérées.
Les ordres de grandeurs des débits et vitesses dans les diffé- 2.1 État énergétique
rents composants d’un réacteur sont facilement évalués à l’aide
des caractéristiques de fonctionnement. Le débit volumique 2.1.1 Enthalpie du fluide
dans la cuve du réacteur est celui en provenance de toutes les
boucles froides (sur le parc en France trois boucles pour les L’enthalpie d’une substance pure est fonction de deux variables
réacteurs 900 MWe, 4 pour les réacteurs de puissance supé- thermodynamiques. Par exemple, considérant la température T et
rieure ou égale à 1 300 MWe). Cependant une petite partie de ce la pression P, on écrit :
débit (environ 5 %) contourne le cœur soit directement par le
haut vers le dôme, soit dans le cloisonnement ou encore dans
les tubes-guide des assemblages combustible. Du fait de la loi

de la conservation de la masse, c’est le débit massique qui est
conservé constant dans tout le circuit. Les valeurs moyennes de avec le coefficient de dilatation volu-
débit massique qui passe dans un assemblage ou dans un tube
de GV permettent de calculer les vitesses massiques mique isobare exprimé en K–1.
moyennes ; rapport du débit massique à la section de passage.
À pression constante, le second terme s’annule et l’augmenta-
La vitesse massique est une grandeur plus adaptée à la ther-
tion de la température T1 à T2 d’une unité de masse d’eau se tra-
mohydraulique que la vitesse moyenne, usuellement utilisée en duit par une augmentation de l’enthalpie h suivant :
hydraulique, qui, dans le cas d’un écoulement à masse volu-
mique variable, n’est pas constante tout au long de l’écoule- (2)
ment. La vitesse massique, notée G, est le produit de la vitesse
U du fluide par sa masse volumique ρ et permet de vérifier, en avec cp (P, T) la capacité thermique massique à pression
écoulement monodimensionnel permanent, l’équation de constante, exprimée en J/(kg.K),
conservation de la masse (ρ U = constante).
h exprimée en J/kg.
On évalue aussi la puissance moyenne, la puissance linéique
Si la plage de variation de température est limitée, en prenant la
moyenne et le flux surfacique moyen d’un crayon d’assemblage
valeur moyenne de la capacité massique, cette expression se sim-
combustible à partir du nombre total d’assemblages dans le cœur
plifie selon :
et du nombre de crayons combustible dans un assemblage.
(3)
Exemple 2
On utilise les valeurs du palier 1 450 MWe [1] et les grandeurs ther- Cette simplification n’est valable que si l’évolution de température
modynamiques du tableau 1 à 15,5 MPa, 292,2 °C et 329,6 °C. ne passe pas par la température de saturation, en raison des
grandes variations des valeurs de la capacité thermique entre la
Le débit nominal par pompe dans une boucle est de 24 500 m3/h phase liquide et vapeur. Si le fluide atteint la température de satura-
dans un réacteur à 4 boucles. Le débit dans la cuve est de 98 000 m3/h tion, il peut devenir diphasique. Les grandeurs thermodynamiques
(20 199 kg/s) soit 93,6 kg/s dans chacun des 205 assemblages consti- de chacune des deux phases ne sont alors plus fonction que d’une
tuant le cœur (avec une température d’entrée de 292,2 °C et 95 % du seule variable ; la présence d’une phase supplémentaire ayant fait
débit passant effectivement dans les assemblages). Dans un des perdre un degré de liberté à l’eau qui est alors à saturation.
5 599 tubes d’un générateur de vapeur (GV), le débit massique est de
0,902 kg/s. À ce débit correspond une « vitesse massique » moyenne Lorsque h(P, T) atteint la valeur de l’enthalpie à saturation
(rapport du débit à la surface de passage) de 4 035 kg/m2.s, dans un ou ( ) , il y a potentiellement passage
tube de GV de 16,87 mm de diamètre intérieur. La vitesse massique d’une fraction de masse de l’eau de l’état liquide à l’état vapeur.
moyenne dans un assemblage est du même ordre de grandeur Ce passage requiert une quantité d’énergie supplémentaire,
(3 894 kg/m2.s) avec une surface de passage dans un assemblage de
0,24 m2 (cf. exemple 1). Dans un GV, entre l’entrée (329,6 °C) et la sor- l’enthalpie (chaleur latente) de vaporisation notée (en J/kg)
tie (292,2 °C) d’un tube, la vitesse moyenne passe de 6,18 à 5,44 m/s. avec :
La vitesse moyenne de l’eau dans un assemblage est de 5,25 m/s en
entrée et de 5,97 m/s en sortie. (4)
La puissance moyenne d’un crayon combustible est de 78 529 W Notons aussi qu’il y a relation biunivoque entre psat et Tsat. Des
sachant que 4 250 MW – thermiques – sont générés par les tables [2] fournissent les grandeurs thermodynamiques de l’eau
205 assemblages qui contiennent chacun 264 crayons. tant à l’état monophasique (liquide ou vapeur) qu’à saturation.
On définit la puissance moyenne linéique comme la puissance Cette configuration cesse d’exister au-delà du point critique
moyenne par unité de longueur. La hauteur active étant de 4 267 mm, (22,06 MPa ou 374 °C).
on obtient : Dans le domaine de fonctionnement du réacteur, les principales
grandeurs thermodynamiques ont les valeurs indiquées dans les
tableaux 1 et 2.
Copyright © – Techniques de l’Ingénieur – Tous droits réservés.

BN 3 050v2 – 3
Tableau 1 – Grandeurs thermodynamiques de l’eau liquide
Coefficient Viscosité
Masse volumique
Pression [MPa] Température [°C] Enthalpie [kJ/kg] de dilatation [10–3 K–1] dynamique
[kg/m3]
relation (126) [μPa.s]
15,5 292,2 1 296 2,585 742 91,6
15,5 329,6 1 514 4,774 653 76,3
7,3 229,5 989 1,622 832 117
Tableau 2 – Grandeurs thermodynamiques de l’eau à saturation
Température Masse Masse Tension de

Pression Enthalpie liquide Enthalpie vapeur
de saturation volumique volumique vapeur σ
[MPa] [kJ/kg] [kJ/kg]
[°C] liquide [kg/m3] vapeur [kg/m3] [N/m]
2,0 212,4 850 10,0 908,6 2 798 0,03483
7,3 288,7 734 38,3 1 283 2 769 0,01697
10,0 311,0 688 55,5 1 408 2 725 0,01186
15,5 344,8 594 102 1 630 2 596 0,004669
Les tables d’eau [2] sont maintenant disponibles sous de nom-

breuses formes numériques, soit avec des applications Internet ou
smartphone, des logiciels ou encore des bibliothèques de pro- Q
grammes.
Hentréee Hsortie
2.1.2 Bilan enthalpique • •
m m
Le premier principe de la thermodynamique est établi à partir
de la conservation de l’énergie (cf. [BE8005] ou [24]). Dans un sys-
tème ouvert représenté schématiquement à la figure 2, on peut
écrire, en supposant les variations d’énergies cinétique et poten-
tielle négligeables devant la chaleur et le travail échangé par le Figure 2 – Bilan enthalpique avec transfert de chaleur et sans travail
échangé
système avec l’extérieur
À noter que s’il n’y a pas de changement de phase et pour des

Dans le cas où il n’y a pas d’échange de travail, par exemple écarts de température limités, le bilan enthalpique peut se simpli-
entre l’entrée et la sortie du fluide primaire dans le cœur ou le fier comme suit : .
générateur de vapeur, alors la relation précédente se simplifie
selon : Exemple 3 : débit primaire
En prenant comme hypothèse, à une pression de 15,5 MPa, une
température d’entrée du cœur de 295 °C et une température de sor-
soit par unité de temps : tie de 329 °C,
La quantité de chaleur par unité de temps correspond à la puis-

sance échangée entre le système et l’extérieur, ce raisonnement Si la puissance thermique nominale est de 4 300 MW (thermiques)
étant aussi valable si cet échange se fait avec un échange de tra- alors le débit nominal sera :
vail.
Pour un système échangeant une puissance Π avec l’extérieur,
le bilan enthalpique relie la variation d’enthalpie à la puissance
échangée selon : À l’inverse, si les températures sont mesurées en entrée et sortie du
GV, en l’absence de perte thermique dans le GV, la connaissance de
la puissance thermique du GV permet de déduire le débit primaire
(5)
passant dans le GV.

BN 3 050v2 – 4
Exemple 4 : besoin en eau d’une centrale porelles et spatiales (voir par exemple [3] [4]) liées aux grandeurs
Dans un cycle thermodynamique, une quantité importante d’énergie diphasiques locales définies ci-après pour des écoulements per-
est perdue à la source froide ; dans le cas d’un circuit ouvert qui évacue manents à une direction. Nous décrivons ici des expressions sim-
cette énergie directement à la source froide (océan par exemple), on plifiées des différentes variables utilisées pour caractériser les
peut calculer le débit nécessaire au condenseur suivant la relation sui- écoulements diphasiques. Il s’agit de déterminer la proportion
vante intégrant le rendement η du cycle thermodynamique. relative des phases, la vitesse, la pression et l’état énergétique de
chaque phase. Au contraire de présentations plus théoriques où
Puissance électrique et puissance thermique sont reliés selon : l’on considère des grandeurs tridimensionnelles avec des phases
indifférenciées, les grandeurs monodimensionnelles introduites
par la suite caractérisent préférentiellement la vapeur (ou le gaz),
La puissance perdue à la source froide Πp est la puissance totale repérée par l’indice « v », par rapport au liquide (indice « l »), sans
(thermique) moins la puissance électrique : opérateurs de moyenne temporelle ou spatiale, approche moins
rigoureuse mais qui a le mérite de la simplicité.
En utilisant la relation (3) pour évaluer la variation d’enthalpie au tra- 2.2.1 Proportion des phases : taux de vide (taux
vers de la différence de température de l’eau de la source froide de présence de la phase vapeur) et qualité
entre l’entrée et la sortie du condenseur ;
Le titre thermodynamique à l’équilibre est issu d’un bilan éner-
gétique. Cependant sa connaissance ne donne qu’une évaluation
grossière de la présence effective de vapeur dans le liquide. Deux
Pour une augmentation de la température issue de la source froide grandeurs supplémentaires permettent de définir un peu plus pré-
de 10 °C, le cp peut être pris constant égal à 4,18 kJ/kg K. Pour une cisément la présence des phases :
centrale produisant une puissance de 1 400 MW avec un rendement

– la qualité statique (quelquefois en thermodynamique appelé
de 33 %, le débit provenant de la source froide est alors :
titre massique en vapeur) xst, rapport de la masse en phase vapeur
à la masse totale de fluide ;
– le taux de présence de la phase vapeur α, souvent dénommé
soit environ 68 m3/s. taux de vide, rapport entre le volume occupé par la phase vapeur
au volume total occupé par le fluide :
2.1.3 Titre thermodynamique à l’équilibre
(7)
L’enthalpie est une grandeur thermodynamique extensive, elle
est proportionnelle à la masse, à la différence de la pression ou de (8)
la température, grandeurs intensives. Lorsque l’enthalpie d’une
unité de masse de liquide est inférieure à l’enthalpie à saturation
(toujours donnée dans les tables d’eau par unité de masse), le Dans les écoulements diphasiques en conduite à une direc-
liquide est dit sous-saturé. L’augmentation de la température du tion, le taux de présence de la phase vapeur, dénommé par la
liquide ou la diminution de sa pression le porte à saturation. Tout suite taux de vide α, peut être défini comme un rapport de sur-
apport d’énergie complémentaire entraîne alors une vaporisation : faces, rapport de la section courante occupée par la phase
le fluide est saturé diphasique (figure 3). vapeur à la surface totale de la section.
Ces étapes sont résumées ci-dessous, pour une unité de masse
du fluide :
On obtient alors les relations suivantes (les surfaces sont notées
– fluide à l’état monophasique liquide sous-saturé : A pour area afin de ne pas provoquer de confusion avec S pour
– fluide saturé diphasique : et slip introduit ultérieurement)
– fluide à l’état de vapeur surchauffée : (9)

À saturation, on a alors :
La masse volumique moyenne est alors définie comme étant :
ou encore : (10)
(6) 2.2.1.1 Débits et vitesses

Cependant si la vapeur occupe une proportion α du volume
total, cela est insuffisant pour connaître le débit de vapeur, lié à la
χ est appelé le titre thermodynamique à l’équilibre, avec vitesse de l’écoulement. Le taux de débit volumique de vapeur β,
l’hypothèse que les phases liquide et vapeur sont toutes deux à rapport du débit volumique de vapeur au débit total, est une
la température de saturation. caractéristique supplémentaire de l’écoulement :
C’est en fait une hypothèse assez forte et non vérifiée dans les (11)
transitoires. Le titre thermodynamique à l’équilibre permet d’éva- On définit le facteur d’écoulement C comme :
luer le potentiel énergétique et le taux de passage du fluide de
l’état liquide à l’état vapeur. (12)
La vitesse moyenne de l’écoulement est le rapport du débit

2.2 Grandeurs diphasiques volumique à la surface de passage. Dans un écoulement dipha-
sique à une direction, la vitesse moyenne de chaque phase corres-
La description complète de l’écoulement requiert des informa- pond donc au rapport du débit de la phase considérée à la surface
tions tridimensionnelles dépendantes du temps et de l’espace. Les de passage occupée par la phase. Le taux de vide α étant par
écoulements diphasiques se caractérisent par des moyennes tem- extension le rapport de la surface occupée par la phase vapeur à

BN 3 050v2 – 5
χ>1
Monophasique vapeur
sat
h ≥ hl + hlv
Écoulement
surchauffé
Z [m]
Témpérature de saturation
Écoulement
saturé
Ébullition franche
sat 0≤χ≤1
h ≥ hl + χhlv
Témpérature à la paroi
Témpérature moyenne du fluide
Ébullition nucléée
χ<0
sat
h < hl
T [°C]
Monophasique liquide Écoulement

sous saturé
Figure 3 – Écoulement dans un tube chauffé
la surface totale, on en déduit alors pour chaque phase, les Et l’on peut former aussi les relations relatives à chaque phase
vitesses moyennes respectives de la vapeur et du liquide : selon :
(13) (19)
(14) (20)
La vitesse « superficielle » (appelée aussi flux volumique) repré- Un écoulement à masse volumique variable (et donc particuliè-
sente le rapport du débit volumique à la surface totale, grandeurs rement un écoulement diphasique) est donc susceptible d’avoir de
« virtuelles » qui permettent de s’affranchir de la connaissance du larges variations de débit volumique et en conséquence de vitesse
taux de vide selon : de chaque phase.
On peut introduire la notion de qualité pour évaluer les débits
(15) massiques vapeur et liquide selon :
(16) (21)
Le débit volumique étant rapporté à la surface totale, la relation (22)

suivante est toujours valable :
Cette notion peut être reliée à la définition de la qualité statique
(ou titre massique en vapeur) xst mais prend en compte l’aspect
(17) dynamique de l’écoulement. En écoulement permanent et lorsque
les vitesses du liquide et de la vapeur sont identiques, les valeurs
Le produit du débit volumique et de la masse volumique forme le de la qualité « x » et du titre massique en vapeur « xst » peuvent
débit massique qui, conformément à la loi de conservation de la être considérées comme identiques. Il convient de souligner la dif-
masse, est constant tout au long de l’écoulement : férence de définition entre termes débitants (liés aux débits mas-
siques) et statiques (liés aux masses). Cependant, afin de ne pas
(18) alourdir l’exposé et les relations proposées par la suite, nous ne

BN 3 050v2 – 6
singularisons pas les grandeurs par des indices relatifs au carac- des phases liquide et vapeur. De plus si l’on suppose ces phases,
tère statique ou débitant de leur définition. toutes les deux, à la température de saturation, les enthalpies
L’enthalpie du mélange s’écrit alors : vapeur et liquide de la relation (23) sont alors les enthalpies à
saturation. Cependant, si le titre χ est défini tant en écoulement
(23) sous saturé (χ < 0) ou surchauffé (χ > 1), les valeurs respectives de
la qualité sont imposées en écoulement sous-saturé ou surchauffé
La « vitesse massique moyenne » de deux phases G, introduite à (respectivement x = 0 et x = 1). Une qualité négative n’a pas de
la section 1, comme le rapport du débit massique à la surface de sens physique (ce qui indiquerait une masse négative) mais un
passage s’écrit aussi : titre négatif indique le degré de sous-saturation. Plus le titre sera
négatif, plus le liquide sera éloigné de la saturation. Dans le
même ordre d’idée, une qualité supérieure à 1 n’a pas de sens
(24) (l’intégralité du fluide ayant été vaporisée) tandis qu’un titre supé-
rieur à 1 indique le degré de surchauffe de la vapeur.
Les vitesses massiques de chaque phase suivent les expressions :
La relation (28) permet de relier le taux de vide à la qualité
selon :
(30)
(25) v
l
(26)
ou
v
2.2.2 Glissement entre phases
l (31)
L’écoulement chauffé s’accélère car sa masse volumique dimi-
( v )
nue (typiquement au cours de la traversée du cœur en raison de l
l’apport d’énergie au fluide). À débit massique constant, la masse
volumique diminue donc la vitesse augmente pour respecter la loi Ces égalités impliquent la connaissance du glissement, qui est
de conservation de la masse. En cas d’ébullition, la phase vapeur, abordé dans la section suivante.
moins dense, acquiert plus de vitesse que la phase liquide. La
poussée d’Archimède (relative à la masse volumique très diffé- Exemple 5
rente entre les phases liquide et vapeur), le gradient de pression Pour une qualité de 0,1 (10 % du débit masse en phase vapeur), le
et la plus forte résistance à la paroi de la phase liquide participent taux de vide vaut entre 0,3 (pour un glissement de 1,5 et une pres-
également à l’accélération de la phase vapeur. L’écoulement pré- sion de 15,5 MPa) et 0,7 (pour un glissement de 1,0 et une pression
sente alors une phase vapeur ayant une vitesse différente (le plus de 7,0 MPa). La vapeur occupe donc entre 30 à 70 % du volume
souvent supérieure) de celle de la phase liquide. total. Plus la pression est basse (et donc plus le rapport des masses
volumiques liquide et vapeur est élevé) plus le volume occupé par la
vapeur est important. La variation du glissement (ici 1,0 ou 1,5) inter-
Le glissement noté S (conformément au terme anglais slip vient dans des proportions moindres. Les évolutions sont données à
ratio) est défini comme étant le rapport de la vitesse de la la figure 4.
phase vapeur à celle de la phase liquide :
2.3 Régimes d’écoulements

(27)
Outre le volume occupé par la vapeur, il est nécessaire de
avec S ≥ 1 dans une très large majorité des cas et pour un taux de connaître sa répartition dans l’écoulement, d’autant plus dans une
vide significativement non nul. démarche qui consiste à moyenner les grandeurs.
Les relations entre vitesses (13) et (14), débits massiques (19) et En circulation verticale ascendante, on distingue principalement
(20), qualité (21) et (22) permettent d’établir la relation : les écoulements (figure 5) :
– à bulles (bubble) ;
(28) – à bouchons ou poches (slug) ;
– agité (churn) ;
En remarquant qu’avec les relations (11) et (12), on obtient l’égalité : – annulaire (annular).
Pourront s’ajouter des formes additionnelles (annulaire dis-
persé, brouillard) ou spécifiques aux écoulements horizontaux
(ondulé, stratifié…), mais les conditions pour lesquelles l’écoule-
alors ment suit tel ou tel régime dépendent de la géométrie et des para-
mètres de l’écoulement (notamment les débits massiques et
volumiques de chaque phase et la pression). De nombreuses
(29) cartes ont été établies (voir par exemple [5] ou [BE8236]), à partir
d’expériences sur des écoulements adiabatiques diphasiques eau-
Au fur et à mesure que le fluide s’approche du point critique air principalement en conduite circulaire de petit diamètre, mais
(22 MPa ou 374 °C), les caractéristiques de ses phases liquide et aussi quelques expériences eau-vapeur. On trouve aussi des cor-
vapeur deviennent de plus en plus proches. S tend vers 1 quand rélations empiriques qui délimitent les frontières de différents
le fluide s’approche de son point critique. régimes d’écoulement. Elles n’ont qu’une valeur indicative car
correspondant rarement aux conditions exactes recherchées. En
2.2.3 Relation entre qualité et titre écoulement diabatique (avec chauffage), les frontières entres les
régimes d’écoulement sont encore plus imprécises. L’apparition
Des relations (6) et (23), on peut relier le titre thermodynamique de bulles de vapeur en paroi chauffante et leur migration vers le
à l’équilibre χ et la qualité x dans le cas de présence simultanée centre de l’écoulement font que l’écoulement n’est pas développé.

BN 3 050v2 – 7
plus probable, puis à vérifier la cohérence des résultats obtenus et

leur conservatisme par comparaison avec le test des autres
1,00 régimes possibles. À titre indicatif, pour un taux de vide inférieur
à 0,20, on peut raisonnablement considérer un écoulement à
0,90
bulles. Un écoulement annulaire est probable pour des valeurs du
0,80 taux de vide supérieures à 0,80.
L’analyse devient encore plus complexe dans le cas où les deux
0,70
phases s’écoulent dans des directions opposées (contre-courant) :
Taux de vide [–]
0,60 – dans des écoulements verticaux, lorsque la phase vapeur

s’élève tandis que la phase liquide tente de descendre. L’engorge-
0,50 ment, ou remontée du liquide (en anglais flooding) se produit
0,40
lorsque la vitesse superficielle du gaz est trop importante et le
liquide ne peut plus s’écouler à contre-courant. La situation
0,30 inverse, retombée du liquide ou en anglais flow reversal, étant
aussi possible ;
0,20 – dans conduites horizontales ou inclinées, des relations sont
0,10
proposées [6] [7] pour déterminer les conditions d’établissement
Pressions 15,5 et 7 MPa de tels régimes d’écoulement.
0,00 Pour des écoulements verticaux dans des conduites circulaires de
0,00 0,10 0,20 0,30 0,40 0,50 0,60 0,70 0,80 0,90 1,00
diamètre compris dans les limites de 2 à 40 fois l’échelle de Laplace
Qualité massique [–]
7 MPa S = 1 15,5 MPa S = 1 soit , la corrélation de Wallis éta-

7 MPa S = 1,5 15,5 MPa S = 1
blie en 1961, relativement simple à mettre en œuvre, relie les
vitesses superficielles adimensionnelles des deux phases selon :
Figure 4 – Évolution du taux de vide en fonction de la qualité pour
deux valeurs de glissement et deux valeurs de pression
avec m = 1 et C = 1 lorsque l’effet de la gravité est prépondérant

par rapport à l’effet de la viscosité, et
Les grandeurs C et m de la relation de Wallis ci-dessus peuvent

évoluer selon les auteurs et sont généralement comprises entre
0,7 et 1.
Exemple 6
Pour un écoulement vertical eau-vapeur en conduite de 6 cm de
diamètre, à 2 MPa, le débit massique ascendant de vapeur est de
a écoulement à bulles b écoulement à bouchons 0,120 kg/s, on cherche la valeur limite du débit liquide en dessous de
laquelle le liquide ne peut descendre dans la conduite.
Pour C = m = 1, la valeur limite est égale à :
Les propriétés physiques de l’eau à saturation sont les suivantes

(tableau 2) :
Les grandeurs thermohydrauliques suivantes sont calculées :

c écoulement agité d écoulement annulaire
La relation avec C = 1 et m = 1 conduit à :

Figure 5 – Régimes d’écoulements liquide-gaz
Comme indiqué au paragraphe 2.1.3, l’équilibre thermodyna-

mique est rarement vérifié. En outre les informations sur les incer- À noter que représente l’échelle de Laplace
titudes sont manquantes.
Cependant, l’identification du régime d’écoulement est rare- et .
ment indissociable du choix des modèles à utiliser et de la fiabilité
des évaluations qui en résultent. La démarche consiste donc à Le diamètre D = 0,06 m est bien compris entre 2 et 40 fois
essayer en premier lieu le modèle correspondant au régime le l’échelle de Laplace.

BN 3 050v2 – 8
3. Modélisations Quelles que soient les corrélations, lois de fermeture ou

modèles utilisés, il convient d’avoir à l’esprit les hypothèses utili-
sées afin d’appréhender les limitations qui en découlent. Lois de
fermeture et corrélations ne sont applicables que sur des
3.1 Démarche proposée et hypothèses domaines d’application limités, le plus souvent associés à des
associées régimes d’écoulements spécifiques. En outre, il convient d’estimer
les incertitudes associées à l’application des corrélations et lois de
Formellement, les écoulements de fluides multiphasiques fermeture. Des corrélations plus complexes peuvent permettre de
suivent les relations de la mécanique des fluides. Pour mémoire, réduire les incertitudes au risque de n’être applicables que sur un
les relations de conservation de la masse, de la quantité de mou- domaine restreint, entraînant alors des erreurs importantes en cas
vement et de l’énergie (équations de Navier Stokes) peuvent d’application en dehors du domaine de validité.
s’écrire pour un fluide newtonien selon les relations simplifiées
(32) données à titre d’illustration car comportant des hypothèses
complémentaires. On trouvera dans [3] les équations complètes. 3.2 Modèle homogène à l’équilibre (HEM)
En écoulement à direction privilégiée (vertical ascendant), les
deux phases étant en équilibre thermodynamique et mécanique,
le modèle homogène à l’équilibre (HEM) se révèle satisfaisant et
(32) particulièrement utile. Ce modèle correspond dans la littérature
anglo-saxonne au sigle HEM (Homogeneous Equilibrium Model).
Les hypothèses principales considèrent :
– des vitesses identiques pour les deux phases;
– un profil de vitesse uniforme sur la section, seule alors la
Au-delà de la complexité inhérente aux équations aux dérivées

partielles et au couplage de la pression et des composantes de la vitesse moyenne est considérée ;
vitesse, le traitement des écoulements turbulents fait apparaître, – les deux phases sont à la même température correspondant à
du fait de l’utilisation de moyennes temporelles, des termes sup- la température de saturation imposée par la pression locale.
plémentaires (tenseur de Reynolds) pour lesquels des relations de Les vitesses des phases étant identiques, le glissement (27) est
fermeture sont indispensables [BE8157]. Dans le cas de coexis- égal à 1 et la relation (30) se simplifie pour donner une relation
tence de plusieurs phases, les interactions entre phases requièrent directe entre la qualité (par hypothèse strictement identique au
des relations de fermeture complémentaires, notamment relatives titre thermodynamique à l’équilibre χ) et le taux de vide :
au traitement des interfaces. Le modèle le plus adapté à la des-
cription des écoulements est le modèle à deux fluides (§ 3.4) mais
qui demande des outils et moyens de calcul importants. Lorsque (33)
l’on ne recherche pas l’évaluation locale des paramètres et qu’une
estimation en entrée/sortie est suffisante, des corrélations empi- Avec cette expression, la relation (10) se simplifie et s’écrit :
riques, associées la plupart du temps aux modèles présentés dans
la suite, ont été établies (voir § 4). Ces corrélations permettent
(34)
d’évaluer les ordres de grandeur sans imposer la mise en œuvre
de moyens lourds.
On obtient ainsi l’égalité :
Outre le recours aux moyennes temporelles des paramètres,
une source importante d’imprécisions est relative aux moyennes
spatiales. Pour le taux de vide, par exemple dans une conduite cir-
culaire, les valeurs maximales peuvent se situer soit « près » de la
paroi, soit au centre, selon le régime d’écoulement. Pour ce qui Le modèle HEM permet alors de considérer tout écoulement
concerne la vitesse ou la température, on obtient des variations diphasique comme un écoulement simple phase ayant des pro-
importantes entre la valeur « à la paroi » (nulle pour la vitesse), la priétés physiques particulières. La masse volumique étant définie
valeur maximale et la valeur moyenne. Dans le cas d’écoulement (34), diverses expressions de la viscosité et de la conductivité
diphasique, certains modèles font l’hypothèse d’un profil uni- thermique ont été proposées (ces expressions ne sont pas spéci-
forme (vitesse moyenne sur toute la section) soit avec la vitesse fiques au modèle HEM) et sont rassemblées dans [8] et [28]. Par
de chaque phase identique (modèle HEM § 3.2) soit avec des analogie, on peut considérer les relations suivantes qui ont le
vitesses différentes. Soit on postule un profil radial de vitesse, mérite de converger vers les propriétés physiques des phases
identique ou non pour les deux phases (modèle à dérive de flux, seules liquide ou vapeur lorsque le taux de vide tend respective-
§ 3.3). ment vers zéro ou un :
Des configurations ou situations peuvent s’appuyer sur des
hypothèses simplificatrices telles que : (35)
– en écoulement stationnaire, les dérivées partielles par rapport
au temps sont nulles lorsque l’écoulement est établi et stable ;
– les écoulements à direction privilégiée peuvent soit être consi- (36)
dérés comme monodimensionnels, soit conduire à des équations
simplifiées en ne considérant que les termes prédominants ; Il est alors possible d’utiliser toutes les relations et corrélations
– dans un système à plusieurs composants type eau-air, les établies en écoulement simple phase et notamment en conduite
transferts de masse entre les phases sont nuls. Dans un système à circulaire (voir [BE8161]) sans recourir à des lois de fermeture spé-
un composant eau-vapeur en écoulement adiabatique et en cifiques.
l’absence de forts gradients de pression, les transferts de masse
peuvent être négligés ; titre, taux de vide et masse volumique
peuvent alors considérés comme constants ; 3.3 Modèle à dérive de flux
– l’équilibre thermodynamique entraîne l’égalité entre la qualité
et le titre thermodynamique et l’égalité entre les températures des L’approche de Zuber et Finlay [9] prend en compte d’une part
deux phases à la température de saturation. des vitesses locales des phases différentes au travers de leur dif-

BN 3 050v2 – 9
férence (et non plus de leur rapport comme dans l’expression posées afin de pouvoir évaluer les ordres de grandeurs, il est cepen-
(27)) et considère d’autre part un profil de vitesse non uniforme. dant indispensable de vérifier les domaines d’applications de ces
Sur une section, les valeurs moyennes des paramètres, dans relations ainsi que les incertitudes qui sont associées à la détermina-
l’expression (37) f à titre de terme générique, sont notées : tion des valeurs. Le lecteur est invité pour cela à consulter les réfé-
rences citées.
(37) Ishii [14] propose pour le paramètre d’écoulement C0 :
La vitesse de dérive locale modélise la différence entre la valeur (45)

locale de la vitesse de la phase vapeur Uv, relation (13), et de la
vitesse superficielle (ou flux volumique) du mélange diphasique j, avec pour la dérive de flux
relation (17) :
(38) (46)
On en déduit la dérive de flux vapeur, c’est-à-dire le terme dfv

(drift flux) qui indique le « surplus » de vapeur (volumique) qui
passe par rapport au mélange.
(47)
(39)
De la même façon, la dérive de flux liquide dfl peut être définie
selon (48)
(40) Maier et Coddington [15] comparent des modèles adaptés aux

écoulements en grappe simulant un assemblage de REP et pro-
posent des relations simples à partir de régressions linéaires sur
et on obtient alors dfv = –dfl.
les essais :
On définit la vitesse de dérive effective comme l’intégrale de la
dérive de flux normalisé sur le taux de vide, la vitesse de dérive (49)
prend en compte les différences de vitesses entre les phases :
(41)
(50)
On définit le paramètre C0 à partir des paramètres locaux de
l’écoulement intégrés sur la section de passage A, C0 prend en
compte le profil de vitesse et de taux de vide dans la section : Enfin, la corrélation de Bestion [16] est assez simple à mettre en
œuvre et donne des résultats très satisfaisants, particulièrement
(42) dans les conditions des réacteurs à eau sous pression, et notam-
ment dans le cœur.
On obtient alors une relation entre le taux de vide, la vitesse de
dérive et le paramètre C0 selon : (51)
(43)
Exemple 7 : ordre de grandeur de C0 et ud
Si la vitesse de dérive ud est nulle, alors en utilisant l’expression Pour l’eau et aux pressions de 20 et 100 bar, les valeurs des
(17), on remarque que le paramètre C0 est l’inverse du paramètre masses volumiques liquide et vapeur et de la tension superficielle
d’écoulement C défini en (12). sont données dans le tableau 2.
Avec les relations (25) et (26), les opérateurs de moyenne ayant été Certaines relations comportent de plus les valeurs du diamètre ou
supprimés pour faciliter la lisibilité de la relation, on peut écrire : diamètre hydraulique équivalent, la vitesse massique ainsi que la
valeur du taux de vide α (dans ce cas, il est nécessaire d’itérer jusqu’à
(44) convergence). Nous prendrons :
• D (ou Dh) = 0,005 et 0,5 m ;
• G = 20 et 2 000 kg/m2.s.
Il reste à déterminer les paramètres C0 et ud pour que cette expres-
sion soit entièrement explicitée. À partir des résultats d’essais, de Les valeurs du paramètre d’écoulement (tableau 3) sont proches de 1
nombreuses relations ont été proposées, certaines en fonction du et les valeurs de la vitesse de dérive sont données dans le tableau 4.
régime d’écoulement, comportant aussi la valeur du taux de vide (ce
qui demande des itérations et la convergence vers une valeur stable).
La principale limitation réside dans le caractère monodimensionnel
(en conduite ou tube par exemple) des conditions dans lesquelles ces Tableau 3 – Paramètres d’écoulement C0
relations ont été établies, on se réfère alors au modèle à dérive de (sans dimension)
flux à une dimension. Une comparaison exhaustive de huit modèles
peut être trouvée dans [11]. Le modèle de Chexal-Lellouche [12] est Pression Ishii Maier-
considéré comme le plus précis tout en ayant le plus large domaine Bestion
[MPa] (α > 0,30) Coddington
d’application ; ce modèle a ensuite été amélioré et est présenté dans
[13]. Les relations étant relativement complexes, le lecteur est invité à 2 1,18 1,01 1
se reporter aux références citées, en accès libre sur Internet. Nous
présentons quelques relations simples relatives au modèle à dérive 10 1,14 1,03 1
de flux à une dimension pour calculer C0 et ud. Ces relations sont pro-

BN 3 050v2 – 10
rentes. C’est notamment le cas pour l’air avec un liquide ou l’eau

Tableau 4 – Vitesse de dérive (en m/s) avec sa vapeur (hors voisinage du point critique).
Une des limitations liée à l’utilisation de ce modèle réside dans
P [MPa] ; l’indispensable opérateur de moyenne (temporelle et spatiale)
« Ishii » « Ishii » Maier-
D [m] ; Bestion associé à chaque paramètre. Ce sont donc les paramètres moyens
(bouchons) (agité) Coddington
G [kg/m2 s] de l’écoulement qui sont calculés. Si pour des modélisations à
une dimension, une relativement bonne connaissance des phéno-
2 ; 0,005 ; 20 0,077 0,200 0,594 0,381 mènes, obtenue à l’aide de nombreux supports expérimentaux,
permet d’atteindre un niveau satisfaisant de modélisation pour
2 ; 0,5 ; 20 0,771 0,200 0,594 3,808 élaborer les relations de fermeture, des efforts doivent encore être
fournis pour certains transitoires ou dès que les phénomènes
2 ; 0,005 ; 2 000 0,077 0,200 2,330 0,381 deviennent bi- ou tridimensionnels. D’autre part, il est apparu que
le modèle à deux fluides et six équations pouvait manquer de pré-
2 ; 0,5 ; 2 000 0,771 0,200 2,330 3,808 cision lors de la simulation de certains types d’écoulements,
comme par exemple dans le cas d’écoulement mixte annulaire/
10 ; 0,005 ; 20 0,074 0,158 0,138 0,141 gouttes dans lequel le liquide s’écoulant en film en paroi se com-
porte différemment des gouttes se trouvant dans le cœur de
10 ; 0,5 ; 20 0,743 0,158 0,138 1,407 l’écoulement vapeur. Afin de bien décrire ce type d’écoulement,
une solution consiste à diviser la phase liquide en deux « sous-
10 ; 0,005 ; 0,074 0,158 0,606 0,141 phases », chacune possédant son propre taux de vide, son enthal-
2 000 pie et sa vitesse. Cela revient à augmenter le nombre d’équations
du modèle à deux fluides : on parle alors d’un modèle à deux
fluides et neuf équations ou trois champs (bilans de masse,
10 ; 0,5 ; 2 000 0,743 0,158 0,606 1,407

d’énergie et de quantité de mouvement soit trois équations pour
décrire un champ liquide monophasique, trois équations pour un
champ vapeur et trois équations pour un champ gouttes).
On peut donc résumer l’évolution et le futur des outils de simu-
lation appliqués aux écoulements en thermohydraulique :
3.4 Modèles à deux fluides 1/ Simulation monodimensionnelle à trois équations avec équilibre
thermique et mécanique (HEM) ;
Dans le modèle à deux fluides, chaque phase est considérée 2/ Simulation monodimensionnelle à quatre ou cinq équations avec
séparément et décrite avec les équations de conservation de la description de certains déséquilibres thermique ou mécanique. Par
masse, de la quantité de mouvement (vectorielle) et de l’énergie exemple pour la conservation de la masse et ou l’énergie, une équa-
avec des termes d’échange entre les deux phases, soit six équa- tion pour chaque phase avec échanges entre phase, et pour la quan-
tions couplées. On trouvera dans [4] ou [10] les équations consti- tité de mouvement une équation pour le mélange diphasique ;
tutives de tels modèles. Une grande partie des modèles postulent 3/ Simulation monodimensionnelle à six équations avec descrip-
l’équilibre des pressions mais certains auteurs proposent aussi tion des déséquilibres mécanique et thermique ;
d’introduire des pressions distinctes pour chaque phase. Les six 4/ Prise en compte des gaz incondensables tels que l’air ou l’hydro-
équations permettent de déterminer six paramètres inconnus, à gène (ajout d’équations complémentaires) ;
savoir la pression, les vitesses et températures (ou enthalpies) de 5/ Simulation bi- ou tridimensionnelle, avec simplification de cer-
chaque phase et le taux de vide. Cependant, sept relations de fer- tains termes, pour des composants du circuit primaire (downco-
meture relatives au frottement visqueux et au transfert de chaleur mer de la cuve, cœur) ;
dans chaque phase et au transfert de masse, quantité de mouve- 6/ Simulation complète tridimensionnelle à neuf équations ou plus
ment et énergie entre phases sont nécessaires. (phase liquide, phase vapeur, phase dispersée).
Des hypothèses complémentaires peuvent réduire le nombre de Les avancées dans les outils de simulation ne sont pas forcé-
relations de fermeture à modéliser. Dans le modèle complet à ment synonymes de gain en précision ou en robustesse (stabilité
deux fluides et six équations : numérique) mais ont l’avantage d’être plus générales. Il faut de
– l’équilibre thermodynamique relatif à chaque phase est pos- plus évoquer les développements associés à la prise en compte
tulé, par contre il peut y avoir un déséquilibre entre les deux des incertitudes pour associer à la détermination des valeurs des
phases, et donc des températures de chaque phase différentes, paramètres d’intérêt un intervalle de confiance.
pouvant aussi être distinctes de la température de saturation cor-
respondant à la pression locale. Une hypothèse simplificatrice pos-
tule alors que l’une des phases (ou les deux) est à la température
de saturation ; 4. Corrélations spécifiques
– les vitesses de chaque phase sont indépendantes. Si l’on
impose des vitesses égales ou reliées par le glissement ou par un
La prise en compte des variations spatiales et temporelles des
modèle de dérive, le nombre d’inconnues et de relations de ferme-
paramètres d’intérêts ne requiert pas obligatoirement le recours
ture se réduit. À noter que le modèle HEM revient à imposer tem-
aux modèles locaux et ces variations peuvent être évaluées de
pératures et vitesses identiques, soit trois équations et trois
façon macroscopique au travers de corrélations empiriques qui
inconnues (pression, vitesse, taux de vide) et deux relations de fer-
permettent de relier entre elles les grandeurs moyennes à l’équi-
meture (pour la viscosité et le transfert thermique).
libre (titre, qualité, température, enthalpie…). C’est l’objet des
La prise en compte des écoulements hors équilibre et des inte- relations proposées dans cette section qui permettent d’estimer
ractions entre phases (par exemple l’écoulement de la vapeur sur- des ordres de grandeur. On ne cherche pas ici à faire un état de
chauffée le long des crayons combustibles dans le cas d’une l’art des modèles et corrélations pouvant être utilisés mais à pro-
grosse brèche en phase de renoyage cf. 5.2.2) est une avancée poser des relations simples à mettre en œuvre. Ces relations sont
significative de la qualité des outils de calcul par rapport aux cependant sujettes à des fortes incertitudes (de l’ordre de 20 %
modèles simplifiés, HEM notamment. Le modèle à deux fluides mais pouvant atteindre 50 %, voire 100 %). Il est donc recom-
est bien approprié quand les deux fluides sont nettement différen- mandé de bien vérifier le domaine d’application (de validité) de
ciés, notamment lorsque leurs densités respectives sont très diffé- ces relations et les incertitudes associées.

BN 3 050v2 – 11
4.1 Grandeurs thermohydrauliques Pour évaluer le taux de vide correspondant à un titre donné,
pour une pression imposée, on calcule le glissement, les tables
et incertitudes d’eau permettent la connaissance des masses volumiques liquide
et vapeur à saturation, et la relation (30) relie taux de vide et titre
La détermination des grandeurs thermohydrauliques ne se ou qualité (dont les valeurs sont supposées identiques).
limite pas aux valeurs les plus probables mais impose d’évaluer
les incertitudes ou les plages de valeurs possibles. Les sources Il n’est pas indispensable de passer par l’évaluation du glisse-
d’incertitudes proviennent entre autres des hypothèses simplifica- ment pour relier le taux de vide à la qualité. Il est possible par
trices adoptées, des modèles sélectionnés et des incertitudes ou exemple de relier le taux de vide et la qualité en utilisant le facteur
erreurs sur les variables d’entrée. Il convient alors de soigneuse- d’écoulement C et en appliquant la relation :
ment vérifier si les hypothèses choisies sont pertinentes, de sélec-
tionner les modèles adaptés au problème et d’identifier les plages
de validité de ces modèles et les incertitudes associées. Ces infor- (54)
mations ne sont pas toujours disponibles ou suffisamment docu-
mentées. Quelquefois difficilement justifiables, dans la mesure du C est le facteur d’écoulement pour lequel Bankoff propose l’équa-
possible, des recommandations ou ordres de grandeur seront tion suivante :
associés aux corrélations présentées.
Le facteur d’écoulement est semblable au paramètre d’écoule-

4.2 Évaluation du taux de vide ment C0 du modèle à dérive de flux (cf. § 3.3). L’évaluation du
taux de vide peut s’appuyer sur le modèle de taux de vide à partir
La relation (30) permet la détermination directe du taux de
de la relation (44) et d’expressions de la vitesse de dérive et du
vide à partir de la qualité, la connaissance du glissement étant
paramètre d’écoulement adaptées au type d’écoulement attendu.
supposée. Si l’on se place dans les conditions d’équilibre ther-

modynamique, alors la qualité est identique au titre thermody- La référence [17] présente et compare 68 relations permettant
namique à l’équilibre défini dans la relation (6). À l’équilibre de calculer le taux de vide à partir de la qualité. Selon le degré de
thermodynamique, qualité et titre thermodynamiques sont simplicité recherché, les paramètres d’écoulement (pression,
égaux pour des valeurs comprises entre 0 et 1 puisque le titre χ vitesse, configuration), et la précision souhaitée, de nombreuses
indique la fraction d’enthalpie de changement d’état cédée au relations peuvent être choisies. Il est recommandé d’une part de
fluide, donc la fraction de masse d’eau liquide passée à l’état vérifier la formulation des relations et d’autre part d’utiliser plu-
vapeur. Une qualité négative n’a pas de sens physique (ce qui sieurs relations afin d’établir une zone d’incertitude des valeurs
indiquerait une masse négative) mais un titre négatif indique le obtenues ; l’incertitude obtenue n’a pas de justification formelle
degré de sous-saturation. Plus le titre sera négatif, plus le mais donne un ordre de grandeur des variations possibles.
liquide sera éloigné de la saturation. Dans le même ordre
d’idée, une qualité supérieure à 1 n’a pas de sens (le fluide Exemple 8
ayant été intégralement vaporisé) tandis qu’un titre supérieur à On présente aux figures 6, 7 et 8 les variations du taux de vide
1 indique le degré de surchauffe de la vapeur. avec la qualité massique pour divers modèles avec dérive de flux (Ishii
On se placera donc dans la plupart des cas dans l’hypothèse à bulles et à bouchons, Marcian et Coddington et Bestion), le modèle
d’équilibre thermodynamique (les températures des deux phases HEM et les modèles avec glissement de Zivi, Thom et Chilsom.
étant toutes les deux égales à la température de saturation et on À 2 MPa, le taux de vide évolue très rapidement (figure 6a) et vaut,
utilisera la relation entre taux de vide et titre, à saturation). Les quels que soient les modèles, plus de 0,70 pour une qualité massique
aspects hors saturation seront abordés ensuite. de 0,1. Applicable pour les écoulements annulaires, le modèle de Zivi
s’écarte sensiblement des autres pour un taux de vide inférieur à 0,60
Rappelons que l’on propose d’évaluer la valeur moyenne (1D
et n’est pas recommandé dans cette zone. Concernant les autres
sur la section de passage) du taux de vide, et que notamment
modèles, un zoom jusqu’à 1 % en qualité (figure 6b) permet de diffé-
dans le cas d’un écoulement refroidi ou chauffé, le taux de vide
rencier les différents modèles, à qualité identique, les variations de
peut varier sensiblement, en particulier au voisinage de la paroi
taux de vide sont de 0,10 à 0,20 en valeur absolue.
échangeant avec l’extérieur.
La vitesse massique influence significativement les valeurs de taux de
vide dans les modèles de dérive (dans une moindre mesure le modèle
4.2.1 Écoulement saturé de Ishii à « bulles ») ; on remarque que les modèles proposés ne tendent
pas vers un taux de vide de 1,00 pour une qualité massique de 1,00.
4.2.1.1 Écoulement à l’équilibre thermodynamique Pour le modèle de Bestion, seul le domaine aux faibles vitesses mas-
siques, à basse pression, et pour un grand diamètre est concerné.
Dans le cas d’équilibre thermodynamique, le taux de vide peut
être calculé directement à partir du titre thermodynamique (par Pour les modèles de dérive de flux, l’influence de la vitesse mas-
exemple évalué à l’aide d’un bilan enthalpique). Il faut cependant sique se limite aux relativement basses vitesses massiques ; le taux
évaluer le glissement. de vide tend vers une asymptote autour de 1 000 kg/m2.s (figure 7).
Zivi, Thom et Chilsom proposent des expressions très simples, Les modèles de Bestion et d’Ishii à « bouchons » sont aussi fonc-
la première adaptée aux écoulements annulaires sans entraîne- tion du diamètre et son influence est limitée aux très faibles valeurs
ment (Zivi), les deux autres aux écoulements eau/vapeur : de qualité massique (figure 8).
4.2.1.2 Écoulement hors équilibre thermodynamique

(52) Titre, qualité sont des notions particulièrement sujettes à confu-
sion notamment parce qu’elles peuvent être interprétées comme
des paramètres statiques (considérant le volume) ou dynamiques
(considérant l’écoulement), relatives à des masses, des vitesses
massiques ou des enthalpies, indépendamment des opérateurs de
(53) moyenne que l’on peut appliquer. L’évaluation de la quantité en
masse de vapeur (qualité massique) est l’élément primordial du
calcul de taux de vide, et quelles que soient les relations propo-

BN 3 050v2 – 12
1,00 1,00
Pression 2 MPa ; diamètre 0,1 m ; vitesse massique 2 000 kg/m2s
0,90 0,90
X X X X X X X X X X X X X X V I
X X X X X X X X
X X X X X X
X X X X
X X
X X
0,80 X
X
X
X X
0,80
X
X
X X X
X X X
X
X X X
X X X
X X
0,70 0,70
X
X X
X X
X X
X
X
X X
X X
X
X
X
X
X
X
X
Taux de vide [–]
Taux de vide [–]

X
X
X
0,60 0,60 X
X
X
X
X
X
X
X
X
X
X
0,50 0,50 X
X
X
X
X
X
X
X
X
0,40 0,40
X
X
X
X
X
X
X
X
X
0,30 0,30
X
X
X
X
X
X
X
X
0,20 0,20
X
X
X
X
Pression 2 MPa ; diamètre 0,1 m ; vitesse massique 2 000 kg/m2s
X
X
X
X
0,10 0,10
X
X
X
X
X
0,00 0,00
0,00 0,10 0,20 0,30 0,40 0,50 0,60 0,70 0,80 0,90 1,00 0,00 0,01 0,02 0,03 0,04 0,05 0,06 0,07 0,08 0,09 0,10
Qualité massique (–) Qualité massique (–)
Bestion Bestion
Maier Coddington Maier Coddington

Ishii-bulles Ishii-bulles
Ishii-bouchons
Ishii-bouchons
HEM
HEM
Zivi
Thom Thom
Chilsom Chilsom
a sur tout le domaine de qualité massique b à faible qualité massique
Figure 6 – Évolution du taux de vide en fonction de la qualité pour différents modèles (HEM, avec glissement et à dérive) à basse pression et
vitesse massique moyenne
0,35 0,35
0,30 0,30
0,25 0,25
Taux de vide [–]
Taux de vide [–]
0,20 0,20
0,15 0,15
0,10 0,10
Pression 2 MPa ; qualité massique 0,005 ; diamètre 0,1 m Pression 2 MPa ; qualité massique 0,005 ; vitesse massique 2 000 kg/m2s
0,05 0,05
0,00 0,00
0 500 1 000 1 500 2 000 2 500 3 000 3 500 4 000 4 500 5 000 0,000 0,100 0,200 0,300 0,400 0,500 0,600 0,700 0,800 0,900 1,000
Vitesse massique (kg/m2 s) Diamètre [m]
Bestion Bestion
Maier Coddington Maier Coddington
Ishii-bulles Ishii-bulles
Ishii-bouchons Ishii-bouchons
HEM HEM
Zivi Zivi
Thom Thom
Chilsom Chilsom
Figure 7 – Évolution du taux de vide en fonction de la vitesse mas- Figure 8 – Évolution du taux de vide en fonction du diamètre de la
sique pour différents modèles (HEM, avec glissement et à dérive) à conduite pour différents modèles (HEM, avec glissement et à
basse pression et à très faible qualité massique dérive) à basse pression et à vitesse massique moyenne

BN 3 050v2 – 13
sées au paragraphe 4.2.1.1. L’approche consiste à évaluer la quan- de Nu = 455, les bulles se détachent de la paroi. Pour des
tité en masse de vapeur avec la variable χ de la relation (6), à nombres de Péclet élevés, s’ajoutent les conditions d’écoulement,
l’aide de l’enthalpie du mélange que l’on peut estimer avec un l’ensemble est représenté par le nombre de Stanton
bilan enthalpique (cf. § 2.1.2) et les enthalpies à saturation déter-
minées par la pression. Cependant la vapeur peut être à une tem- ; pour St = 1/154, les bulles se détachent.
pérature plus élevée que la température de saturation et dans ce
cas son enthalpie est supérieure à l’enthalpie vapeur à saturation Pour calculer le taux de vide, les relations établies précédem-
donnée par les tables d’eau. Ce peut être le cas aussi pour la ment restent valables, la difficulté provient de l’évaluation de la
phase liquide, et les températures des deux phases peuvent aussi valeur de la qualité. Un bilan enthalpique conduit à une valeur
être différentes de la température de saturation. On peut donc moyenne dans la section du titre thermodynamique inférieure à
généraliser la relation (6) et calculer un titre hors équilibre selon la zéro. Partant de l’hypothèse qu’au point de début d’ébullition zONB
relation ci-dessous en considérant non pas les enthalpies à satura- la qualité massique est tout juste égale à zéro et devient positive,
tion mais les enthalpies réelles des phases liquide et vapeur, on peut estimer que :
(55)
Les enthalpies réelles sont des inconnues supplémentaires qui soit

peuvent être calculées par exemple grâce au modèle à deux
fluides (cf. § 3.4), ce qui impose une mise en œuvre de moyens de
calcul bien plus importants que les relations proposées au para-
graphe 4.2.1.1. Le calcul du titre hors équilibre (et donc du taux de Tout au long de l’écoulement, l’écart entre qualité en masse de
vide en situation hors équilibre) ne peut donc être estimé à l’aide vapeur et titre thermodynamique à l’équilibre se réduit jusqu’à
de relations simples. devenir nul à une distance L estimée (par exemple en fin de
d’écoulement). L’évolution de la qualité peut alors s’écrire selon :
4.2.2 Écoulement non saturé
On observe à la figure 3 la zone d’ébullition nucléée, depuis (58)
l’écoulement simple phase liquide jusqu’à l’écoulement en ébulli-
tion franche. Dans cette zone, des bulles de vapeur sont créées On calcule alors le taux de vide à l’aide de la qualité massique et
bien que l’enthalpie moyenne de l’écoulement ne soit pas à satu- des relations du paragraphe 4.2.1.1.
ration (et que la température moyenne du fluide soit en dessous
de la température de saturation). Dans une première zone, les Exemple 9
bulles sont créées à la paroi (qui du fait du chauffage est à une Le tableau 5 rassemble les données d’entrée et les valeurs calcu-
température supérieure à la température moyenne du fluide), un lées relatives aux conditions utilisées à la figure 9.
peu plus en aval, ces bulles commencent à se détacher de la paroi La valeur de 300 °C correspond à une enthalpie de 1 342 kJ/kg soit
mais coalescent au sein de l’écoulement. On peut donc identifier un titre thermodynamique d’environ –0,05, δONB = +0,05.
deux points remarquables :
– le point de début d’ébullition en paroi (ONB pour « onset of La prise en compte du déséquilibre thermique (températures de
nucleate boiling » en anglais) ; phases différentes de la température de saturation) et du déséqui-
– le point d’apparition significative de vapeur (OSV pour « onset libre mécanique (vitesse des phases différentes) relève plus d’une
of significant void » en anglais), assimilable au point de détache- utilisation dans les logiciels de calcul que de corrélations qui
ment des bulles. pourraient être facilement applicables. Il est alors illusoire d’indi-
Les évolutions le long d’un tube chauffé du taux de vide, du quer des incertitudes à associer aux relations présentées.
titre thermodynamique et de la qualité peuvent être schématisées Il faut noter aussi qu’en cas d’ébullition sous-saturée, une partie
(figure 9). La puissance de chauffage est uniforme tout au long du du flux de chaleur sert à la vaporisation.
tube et a été choisie telle que le titre évolue de –0,20 à 1,00 sur la
hauteur (de l’entrée à 1,2 m). La qualité massique a été estimée, à
titre illustratif, positive dès le point ONB et rejoint le titre thermo- 4.3 Évaluation des pertes de charge
dynamique en fin de longueur. Les évolutions du taux de vide à
partir du titre thermodynamique à l’équilibre et de la qualité esti- Pour vaincre la résistance à l’écoulement d’un fluide, il est indis-
mée peuvent être comparées. pensable de lui fournir une certaine quantité d’énergie pour le
Saha et Zuber [32] ont établi un critère assez simple qui définit
les conditions pour lesquelles les bulles se détachent et de ce fait
amènent à un taux de vide non nul (OSV). Lorsque la température Tableau 5 – Données d’entrées et valeurs calculées
moyenne de l’eau atteint une température limite , le taux de relatives aux conditions utilisées à la figure 9
vide atteint une valeur significativement non nulle :
μ
G φ λ cp
P [MPa] D [m] [10–6 Pa
(56) [kg/m2 s] [kW/m2] [W/m/K] [kJ/kg/K]
s]
10 2 000 0,01 825 86,5 0,548 5,682
(57)
Les propriétés physiques de l’eau sont calculées à la température
Cette limite est conditionnée par la valeur du nombre de Péclet limite d’ONB (300 °C).
(relation (115), Pe, produit des nombres de Reynolds et de Prandtl,
voir § 4.5.1). Aux faibles nombres de Péclet (G faible), cette tem-
Re Pr Pe
pérature limite est gouvernée par les phénomènes de transfert de
chaleur au travers du nombre de Nusselt ; à partir 231 000 0,896 207 300 311 °C 33 °C 11 °C 300 °C

BN 3 050v2 – 14
fluide semble s’écouler par « lames » sans échanger de manière

significative avec les « lames » voisines, d’où la terminologie
1,00
d’écoulement laminaire ;
10 MPa 2 000 kg/m2s – un second type d’écoulement apparaît dès lors que la vitesse
moyenne du fluide augmente, ou dans des conduites de large dia-
0,80
mètre, pour des fluides moins visqueux. Cet écoulement est très
Taux de vide et qualité massique
perturbé, et diffuse très vite, ce sont des écoulements dits turbu-

lents. Entre le laminaire et le turbulent, il existe une zone de transi-
0,60 tion instable.
Outre la mise en évidence de ces deux types d’écoulement,
Reynolds a proposé un critère simple et efficace qui permet de
0,40
déterminer à l’avance le type d’écoulement. En deçà d’une cer-
taine valeur du nombre sans dimension de Reynolds – rela-
tion (113) –, l’écoulement est laminaire, il est turbulent au-delà. La
0,20 valeur précise dépend de la configuration d’écoulement (externe,
interne) et de paramètres secondaires tels que l’état de surface et
les conditions expérimentales (notamment la présence de micro-
0,00 vibrations). Usuellement en écoulement en conduite circulaire
0,000 0,200 0,400 0,600 0,800 1,000 1,200 industrielle, la valeur communément adoptée est 2000 mais on
peut rencontrer des écoulements qui restent laminaires jusqu’à un
–0,20 nombre de Reynolds de 6000 environ, voire 10000.
Hauteur [m] Il n’en demeure pas moins que les deux types d’écoulements
Taux de vide calculé à partir de la qualité massique répondent aux mêmes équations (Navier-Stokes), ce qui peut
estimée apparaître paradoxal. L’explication provient de l’application d’opé-
rateurs de moyenne temporelle sur les paramètres de l’écoule-
Titre thermodynamique à l’équilibre
ment. En effet, en un point fixe et donné de l’écoulement, la
Qualité massique estimée vitesse et la pression ne sont pas strictement constantes et fluc-
Taux de vide calculé à partir du titre tuent autour d’une valeur moyenne. Soit ces fluctuations ne sont
thermodynamique à l’équilibre pas significatives, et dans ce cas l’écoulement est laminaire, soit
les fluctuations ne peuvent être ignorées et du fait du caractère
non linéaire des équations de Navier-Stokes, des termes supplé-
Figure 9 – Évolution du taux de vide en écoulement chauffé sous
saturé mentaires apparaissent qu’il faut modéliser à l’aide de modèles
dits « de fermeture » car ne relevant pas des équations standards
(voir [BE8157] et [A1870]). On aura donc dans le cas d’écoule-
mettre en mouvement. Cette contribution est clairement mise en ments laminaires, la possibilité de résoudre exactement les équa-
évidence dans le bilan de conservation de quantité de mouvement tions de Navier-Stokes et éventuellement d’établir des solutions
des équations de Navier-Stokes (voir par exemple [BE8157]). Son analytiques. Au contraire, pour les écoulements turbulents, plu-
évaluation requiert la résolution d’un système d’équations aux sieurs options s’offrent dans le choix des modèles de fermeture
dérivées partielles. Il est cependant possible de résoudre un cer- (et donc de leur adéquation à l’écoulement considéré). Bien
tain nombres de problèmes plus directement. L’approche propo- entendu, les écoulements multiphasiques turbulents ouvrent la
sée ici est de présenter des modèles simples, empiriques ou voie à des problèmes de complexité fortement accrue. Nous
physiques, pour l’évaluation des pertes d’énergie dues à la visco- allons proposer ici des modèles empiriques adaptés aux écoule-
sité du fluide. ments en conduite, ou des écoulements internes assimilables à
des écoulements en conduite circulaire par l’utilisation du dia-
Osborne Reynolds en 1883 en observant le mouvement des mètre équivalent introduit au paragraphe 1.1. On distinguera les
fluides, a mis en évidence deux types d’écoulements très différents cas laminaire et turbulent dans les écoulements simple phase (cf.
illustrés en figure 10. Dans cette expérience, un colorant est injecté § 4.3.1) et diphasique (cf. § 4.3.2). Une comparaison des modèles
en amont et permet une visualisation directe de sa diffusion en aval. est établie dans [18] et permet d’estimer les incertitudes associées
De ces observations, on définit la classification suivante : à chaque modèle selon les écoulements attendus.
– un type d’écoulement, peu agité, obtenu pour des petits dia- Outre les pertes de charge dues à la viscosité, on observe des
mètres de conduite, des basses vitesses et hautes viscosités. Le pertes de charges dues aux singularités présentes dans l’écoule-
ment. Ces singularités sont les coudes, les vannes, les obstacles
divers et dissipent de manière irréversible l’énergie de l’écoule-
ment par des tourbillons localisés au voisinage de la singularité
(cf. § 4.3.3).
La pression interne dans l’écoulement varie aussi avec la hau-
teur (énergie potentielle) et la vitesse moyenne de l’écoulement
(énergie cinétique) comme établi par l’équation de Bernoulli que
l’on peut introduire sous sa forme généralisée selon :
(59)
La pression totale dans l’écoulement varie donc, d’une part de

façon irréversible suite au frottement et aux tourbillons créés par
les singularités, ce sont des pertes de charge qui dans la relation
sont globalisées dans le terme hμ et, d’autre part, la pression
totale dans l’écoulement varie de façon réversible suite aux varia-
tions de hauteur ou de vitesse (pertes de pression par gravité (cf.
Figure 10 – Écoulement laminaire (en haut) et turbulent (en bas) § 4.3.4) ou par accélération (cf. § 4.3.5).

BN 3 050v2 – 15
Notons que la masse volumique intervient aussi dans le calcul Le coefficient de frottement est alors déterminé analytiquement et
des pertes de charge. Dans le cas d’un écoulement à masse volu- s’écrit :
mique variable, cet élément doit être pris en compte (cf. § 4.3.6).
(66)
4.3.1 Pertes de charge par frottement en
écoulement simple phase Dans le cas d’une conduite quelconque que l’on peut assimiler à
une conduite circulaire à l’aide du diamètre hydraulique (cf. § 1.1),
Les bases théoriques sont présentées dans [BE8161] et [A1870]. on peut alors écrire :
Le modèle de Darcy-Weisbach est communément utilisé et relie la
perte de charge à la vitesse moyenne du fluide (ou son débit volu- (67)
mique) d’après :
avec Ud la vitesse débitante (en m/s).
(60)
4.3.1.2 Coefficient de frottement en écoulement turbulent
La perte de charge par frottement, normalisée par la pression
Lorsque l’écoulement devient turbulent (Re > 8 000), le coeffi-
dynamique ρU2/2, est proportionnelle à la longueur L sur laquelle
cient de perte de charge est déterminé par des relations empi-
la perte de charge par frottement est calculée. Cette longueur est
riques qui prennent en compte l’état de surface de la conduite
normalisée par le diamètre hydraulique. Le facteur de proportion-
avec une rugosité adimensionnelle ε/D formée par le rapport de la
nalité f est le facteur de Darcy (sans dimension) dont les
variation moyenne de l’épaisseur de la conduite au diamètre total
méthodes d’évaluation sont explicités ci-dessous.
de celle-ci. Si les variations sont nulles, la conduite est considérée
Dans le cas d’une conduite circulaire, la section de passage est comme une conduite lisse. Divers auteurs ont proposé des rela-
πD2/4 et l’on peut alors relier la perte de charge par frottement au
tions empiriques parmi lesquelles nous citons les plus couram-

débit volumique selon : ment utilisées.
Pour une conduite lisse et un nombre de Reynolds pas trop
(61) élevé, Blasius propose :
Dans un écoulement à masse volumique variable, on utilise préfé- (68)

rentiellement la vitesse massique ou le débit massique, la relation
(60) s’écrit alors : Pour une conduite lisse et un nombre de Reynolds élevé, on peut
utiliser la relation de Mc Adams :
(62) (69)
Enfin pour des conduites rugueuses, la relation de Colebrook est
Et dans le cas d’une conduite circulaire en prenant le débit mas- très souvent proposée bien que complexe à utiliser du fait d’itéra-
sique comme variable : tions nécessaires jusqu’à obtention de l’égalité des deux termes
de la relation (70) :
(63)
(70)
La formulation de Darcy-Weibach est semblable à une formulation
proposée par Fanning ou Moody et qui relie la perte de charge au Les méthodes de calcul des pertes de charge reposent sur de
coefficient de trainée (noté usuellement Cx ou CD), coefficient de nombreuses mesures expérimentales et une longue pratique.
trainée qui est formé par le rapport du cisaillement en paroi à la Cependant, une incertitude d’au moins 10 % est à considérer sans
pression dynamique 2τp/ρv2 oublier que la détermination de la rugosité, par exemple, n’est pas
immédiate et peut varier avec le temps compte tenu de l’encrasse-
(64) ment de la conduite. De plus dans des géométries différentes de
la conduite cylindrique, l’évaluation du diamètre hydraulique équi-
Dans de nombreux ouvrages cependant le coefficient de trainée valent telle que proposée par Chézy (relation (1)) est assez som-
CD est noté « f » ce qui induit une confusion avec la formulation maire. Cependant l’incertitude sur les pertes de charge ne se
de Darcy. traduit pas directement sur le débit, comme explicité au para-
graphe 4.3.7.
Quelle que soit l’approche suivie (Darcy ou Fanning), il reste à
déterminer le facteur « f ». Les valeurs proposées sont relatives à
l’équation de Darcy-Weisbach. 4.3.2 Pertes de charge par frottement
en écoulement diphasique
4.3.1.1 Coefficient de frottement en écoulement laminaire
En écoulement diphasique, les pertes de charges sont plus éle-
Dans un écoulement laminaire en conduite circulaire de rayon vées qu’en écoulement simple phase, le frottement interfacial
R, compte tenu des symétries et des conditions aux limites impo- s’ajoutant au frottement pariétal. Le taux de présence de chaque
sant une vitesse nulle à la paroi, le profil de l’écoulement est para- phase ainsi que les configurations d’écoulement sont alors des
bolique et unique fonction de la distance « r » entre le centre de la grandeurs d’influence. On introduit comme paramètre significatif
conduite (où la vitesse est maximale) et la paroi (où r = R) : la qualité massique « x » et les modèles présentés seront adaptés
à des configurations d’écoulement spécifiques. Dans un écoule-
ment diphasique, l’écoulement laminaire peut exister mais est
(65)
rarement rencontré dans les applications nucléaires. Si tel était le
cas, les méthodes présentées pour les écoulements turbulents
avec Q débit volumique (en m3/s). s’appliquent, en modifiant ce qu’il convient de modifier. Le
Cette relation est une relation analytique qui découle des équa- modèle HEM (§ 3.2) permet une évaluation simple et rapide des
tions de Navier-Stokes (elle est donc valable tout le temps tant pertes de charge. Adapté aux écoulements à bulles, il ne peut être
que l’écoulement est laminaire). étendu aux autres écoulements sans introduire des erreurs signifi-

BN 3 050v2 – 16
catives. Pour les écoulements annulaires, on donnera la préfé- Il suffit de calculer la perte de charge en écoulement simple
rence aux modèles à fluides séparés (§ 4.3.2.2). Enfin des facteurs phase (liquide ou vapeur), à la vitesse massique correspondante,
multiplicatifs empiriques ont été établis et permettent une évalua- selon les méthodes présentées au paragraphe 4.3.1 et d’appliquer
tion rapide, quoique entachée de fortes incertitudes. Ces trois le facteur multiplicatif. Le facteur multiplicatif est déterminé empi-
approchent sont détaillées dans la suite. On se place dans un riquement selon les formules générales suivantes :
écoulement à qualité constante ; si la qualité massique varie dans
l’écoulement, on se reporte au cas d’un écoulement avec masse (76)
volumique variable traité au paragraphe 4.3.6.
(77)
4.3.2.1 Utilisation du modèle HEM
Le modèle HEM, particulièrement adapté aux écoulements à C est une constante, X est la variable de la relation (73). Cette
bulles (un taux de vide α < 0,25-0,30 est un ordre de grandeur variable peut s’exprimer en fonction du titre massique et des rap-
acceptable), permet de considérer un écoulement diphasique ports des coefficients de frottement et des masses volumiques
comme un écoulement simple phase d’un fluide équivalent ayant liquide et vapeur à saturation selon :
des propriétés physiques spécifiques.
La masse volumique de ce fluide équivalent est calculée à l’aide de
la relation (34) établie analytiquement. La viscosité peut être calculée (78)
à l’aide de la relation (35) mais d’autres relations sont proposées par
divers auteurs sans qu’il soit possible d’en privilégier une.
Le calcul des pertes de charge par frottement se déroule Par exemple, pour une conduite lisse en écoulement des phases
comme dans le cas simple phase : « séparées » liquide et vapeur turbulent à fort nombre de Rey-
– calcul du nombre de Reynolds considérant la vitesse massique nolds, on peut utiliser la relation de Mac Adams (69). Alors la
et la viscosité du fluide équivalent : perte de charge de la phase « séparée » liquide s’écrit (relation
(62)) :
(79)
– selon la valeur du nombre de Reynolds et dans le cas d’un
écoulement turbulent, de la rugosité de la conduite, calcul du coef- et pour la phase « séparée » vapeur :
ficient de frottement f selon les relations (67) à (70) ;
– application de la relation de Darcy-Weisbach (62) avec la
vitesse massique G. (80)
4.3.2.2 Modèle à fluides séparés La variable X, ou paramètre de Martinelli, racine de la variable

Ce modèle est présenté dans [BE8236], on en donne ci-dessous exprimée dans la relation (73) qui forme le rapport de la perte de
les principales caractéristiques. L’hypothèse principale de cette charge liquide à celle vapeur, s’écrit
méthode, proposée par Lockhart et Martinelli, est celle d’un écou-
lement séparé des phases liquide et vapeur. La perte de charge de (81)
l’écoulement diphasique est calculée en fonction de la perte de
charge d’une phase isolée (liquide ou vapeur) ayant le même
débit que dans l’écoulement diphasique. On établit donc les rap- Pour des écoulements turbulents en conduite lisse de nombre
ports de la perte de charge par unité de longueur (dp/dz) en dipha- de Reynolds modéré, on utilise la relation de Blasius (68) et la
sique à celle en simple phase vapeur ou liquide selon les variables variable X s’écrit :
suivantes ainsi définies :
(82)
(71)
La variable X étant ainsi déterminée (ou plus généralement à
l’aide du rapport des pertes de charge de chaque phase s’écoulant
séparément), une valeur de la constante C, évaluée par les expé-
riences, a été proposée. En écoulement turbulent pour les deux
(72)
phases, on indice la variable X selon Xtt et on a :
(83)
Suivant le type d’écoulement des phases liquide et vapeur (lami-
(73) naire ou turbulent), quatre combinaisons de la variable X sont
possibles, et la constante C a les valeurs différentes données dans
le tableau 6.
Et la perte de charge en écoulement diphasique est reliée à celle
de l’écoulement en phase liquide selon : On peut aussi évaluer le taux de vide à partir des grandeurs C et
X selon :
(74)
(84)
ou alors en considérant que la phase vapeur s’écoule seule :
4.3.2.3 Facteurs multiplicatifs
(75) Pour évaluer les pertes de charge d’un écoulement diphasique,

il est aussi possible de calculer la perte de charge dans un écoule-

BN 3 050v2 – 17
Supposant l’égalité des coefficients de frottement en simple phase

Tableau 6 – Valeurs du coefficient C du modèle fl et en diphasique f2φ, on a :
de Lockhart and Martinelli
(88)
Écoulement Écoulement
C
de la phase liquide de la phase vapeur
On peut aussi supposer des modèles identiques du coefficient de
Turbulent Turbulent 20 frottement en simple phase et diphasique en fonction du nombre
de Reynolds. Dans ce cas, prenant par exemple un modèle utili-
Laminaire Turbulent 12 sant l’expression de Blasius (68), avec pour l’expression de la vis-
cosité dynamique diphasique la relation (35), alors :
Turbulent Laminaire 10
Laminaire Laminaire 5 (89)
Enfin, Martinelli-Nelson et Thom [23] proposent des valeurs

ment liquide seul (indice « lo ») pour « liquid only » et d’appliquer tabulées (tableau 7) en fonction de la qualité pour différentes
un facteur multiplicatif. Pertes de charges diphasique et simple pressions. Les valeurs proposées par Thom reproduites dans le
phase liquide sont reliés selon l’expression : tableau 7 sont celles révisées sur la base de nombreux essais réa-
lisés à Cambridge et sont plus élevées que les valeurs d’origine,
(85) notamment à basse pression. Cependant, les valeurs obtenues
avec Martinelli-Nelson sont toujours plus élevées. Les autres
À partir de la relation de Darcy-Weisbach (60), le facteur peut

modèles (Thom ou relations (88) ou (89)) sont plus proches des
s’écrire :
multiples données expérimentales réalisées ultérieurement.
Les incertitudes liées à l’évaluation des pertes de charge en
(86) écoulement diphasique sont considérables, les pertes de charge
étant particulièrement sensibles à la proportion de vapeur, à la
Si l’on prend comme expression de la masse volumique dipha- configuration d’écoulement (interfaces) et aux hypothèses asso-
sique celle établie avec l’hypothèse HEM (34), alors ciées choisies par l’utilisateur. Un ordre de grandeur de l’incerti-
tude peut être estimé en considérant successivement plusieurs
(87) hypothèses et modèles et faisant varier les paramètres d’entrée
afin d’obtenir une plage de variation.
Tableau 7 – Valeurs du coefficient multiplicatif de perte de charge

pour les modèles de Martinelli-Nelson et de Thom (d’après [5])
Modèle de Martinelli Nelson (1948) Modèle Thom (1964)
Pression …...... [MPa]

4,1 8,6 14,5 20,7 4,1 8,6 14,5 20,7
Qualité …...... [%]
1 1,76 1,48 1,18 1,05 1,46 1,10
5 4,96 3,00 1,69 1,17 2,86 1,62 1,21 1,02
10 8,20 4,40 2,31 1,30 4,78 2,39 1,48 1,08
20 14,68 6,85 3,04 1,51 8,42 3,77 2,02 1,24
30 20,72 9,20 3,76 1,68 12,1 5,17 2,57 1,40
40 26,24 11,40 4,46 1,83 15,8 6,59 3,12 1,57
50 31,32 13,45 5,15 1,97 19,5 8,03 3,69 1,73
60 35,88 15,25 5,81 2,10 23,2 9,49 4,27 1,88
70 39,96 16,75 6,43 2,23 26,9 10,19 4,86 2,03
80 43,46 17,85 6,99 2,35 30,7 12,4 5,45 2,18
90 42,86 17,65 7,20 2,38 34,5 13,8 6,05 2,33
100 27,00 11,80 5,46 2,15 38,3 15,3 6,66 2,48

BN 3 050v2 – 18
4.3.3 Pertes de charge singulières pression selon la relation (59). On a alors, en négligeant la perte
de charge singulière et pour un écoulement à vitesse constante :
Les pertes de charges dues aux singularités présentes dans
l’écoulement sont calculées de façon similaire aux pertes de (91)
charges par frottement. L’expression générale suit la formulation
de la relation (62) selon : Pour un écoulement à masse volumique constante, le calcul est
immédiat. Tout comme les pertes de pression par accélération, la
(90) variation de pression est réversible selon que l’écoulement est
ascendant ou descendant.
k est le coefficient de perte de charge déterminé expérimentalement
ou estimé par exemple à l’aide de principes proposés par Idel’chik
dans [19] et [20]. La vitesse moyenne, et donc le débit, peuvent aussi
4.3.5 Pertes de pression par accélération
être déterminés à l’aide de la mesure de la perte de charge singulière, Toute variation de la surface de passage induit une variation de
après calibration. Le diaphragme est un appareil de mesure couram- la vitesse moyenne et de la vitesse massique du fluide pour res-
ment utilisé. Le tableau 8 et [A1870] proposent des ordres de gran- pecter la loi de conservation de la masse. De la même manière,
deurs de singularités usuelles. À noter que dans le cas d’une
toute variation de la masse volumique entraîne une variation de la
réduction ou d’un élargissement brusque, comme la section de pas-
vitesse moyenne, mais la vitesse massique reste constante. Dans
sage varie, la vitesse du fluide varie en proportion du rapport des sur-
un écoulement en conduite, considérons deux points (1) et (2) où
faces de passage et il convient de prendre en compte la variation de
l’énergie cinétique (cf. § 4.3.5). Le tableau 8 indique uniquement la les sections de passage valent respectivement A1 et A2 et les
perte de charge due à la singularité (et donc irréversible) et relative à masses volumiques du fluide ρ1 et ρ2, alors :
la vitesse dans la section la plus petite.
Selon le type de singularité, les incertitudes sur les valeurs des (92)
coefficients de perte de charge sont variables (de quelques pour-
cents à plusieurs dizaines de pourcent). Une source supplémen- Outre la conservation de la masse, la conservation de l’énergie
taire et significative d’incertitude provient de la succession de implique une variation de pression consécutive à la variation de
plusieurs singularités, chaque singularité en amont perturbant vitesse et donc d’énergie cinétique (relation (59)). On a alors, en
l’écoulement et donc la perte de charge relative à la singularité en négligeant la perte de charge singulière et pour un écoulement
aval. Plus généralement, les coefficients proposés sont proposées horizontal :
pour des écoulements établis, ce qui n’est pas toujours validé
dans des circuits hydrauliques industriels. (93)
Si la masse volumique ne varie pas, alors la relation (92) devient :
4.3.4 Pertes de pression par gravité
Lorsque l’écoulement n’est pas horizontal, l’énergie potentielle
varie et la conservation de l’énergie implique une variation de (94)
Connaissant le rapport des surfaces, la mesure de la différence de

Tableau 8 – Ordre de grandeur des coefficients pression au droit des sections (1) et (2) permet d’évaluer la vitesse
de perte de charge pour des singularités usuelles moyenne et donc le débit ; c’est le principe du venturi.
De petites variations de la masse volumique (écoulement
Type de singularité Coefficient k simple phase) peuvent être prises en compte en les linéarisant. Le
cas avec masse volumique variable avec variation de la qualité est
Coude 90o 0,3 – 1,5 détaillé au paragraphe 4.3.6.
Coude 45o 0,2 – 0,4 À noter que l’on parle ici de perte de pression ou perte de
charge réversible car au contraire des pertes de charges par frot-
Coude 180o 0,2 – 1,5 tement, ces variations de pressions peuvent conduire à une baisse
ou une augmentation de la pression locale.
Vanne ouverte 100 % 0,05 – 10
4.3.6 Cas de la masse volumique variable
Vanne ouverte 20 % 0,20 – 200
Un écoulement diphasique chauffé ou refroidi voit sa masse
Embranchement en « T » (90°) dans le sens de 0,2 – 0,9 volumique varier. Le fluide peut potentiellement changer de phase
l’écoulement et avec l’évolution du taux de vide, la masse volumique varie
selon la relation (10). Dans ce cas, les relations établies entre qua-
Embranchement en « T » (90°) hors sens de 1 – 2 lité et taux de vide (cf. 4.2.1) permettent de calculer la masse volu-
l’écoulement mique et donc les pertes de charge.
Élargissement brusque de 100 % 0,50 – 0,60 En partant de la conservation de la quantité de mouvement (32),
en écoulement permanent unidimensionnel, on obtient :
Élargissement brusque 20 % 0,10 – 0,15
Rétrécissement brusque de 100 % 0,15 – 0,20

La variation de pression selon z peut s’interpréter comme la
Rétrécissement brusque de 20 % 0,03 – 0,06
somme de trois composantes, l’une irréversible et due à la visco-
sité (perte de charge par frottement), une autre provenant de la
Rétrécissement ou élargissement inférieur à 5 % Négligeable
gravité et la dernière dérivée de l’accélération.

BN 3 050v2 – 19
4.3.6.1 Perte de charge par frottement Des abaques comportant les valeurs de r3 et les valeurs numé-
Pour un écoulement simple phase avec masse volumique riques sont données dans [23]. Une partie des valeurs numériques
variable, on peut utiliser les relations établies au paragraphe 4.3.1 est reproduite dans le tableau 9 avec les valeurs correspondantes
en prenant sur une longueur adaptée, la masse volumique évaluées à l’aide de la relation (95). Les écarts ne dépassent pas
moyenne. Lorsque l’écoulement est saturé, l’évolution linéaire de 10 % et sont dans la majorité des cas à +/–2 %.
la masse volumique n’est plus vérifiée (cf. relation (10)) et on a les
conditions aux limites suivantes : 4.3.6.2 Variation de pression due à l’accélération
La variation de pression due à l’accélération a ici des origines
différentes de celle évaluée au paragraphe 4.3.5. L’accélération du
fluide est causée non par une variation de la section de passage
mais par la variation de la masse volumique, et par conséquent,
de la vitesse, la vitesse massique restant constante. En écoule-
ment unidimensionnel permanent, la contribution élémentaire à la
Avec une évolution linéaire de la qualité et la relation (34) valable variation de pression peut s’écrire en prenant la valeur moyenne
pour l’hypothèse HEM – ce qui revient à considérer une évolution de la vitesse massique tant en écoulement simple phase que
linéaire du volume spécifique –, on peut simplement prendre la valeur diphasique :
moyenne de la masse volumique calculée selon :
(96)
Une autre option considère les phases liquide et vapeur séparé-

Pour évaluer le frottement avec une variation linéaire de la qualité ment :
et l’hypothèse HEM, on obtient la perte de charge intégrée sur la hau-

teur L pour laquelle la qualité passe de zéro à xs en sortie selon :
(95) (97)
Thom [23] propose un facteur global r3 suivant :

L’utilisation de la relation (96) permet une intégration immédiate
en prenant les volumes spécifiques (inverse de la masse volu-
mique). Dans la partie simple phase, entre les sections 1 et 2, la
Tableau 9 – Valeurs du coefficient multiplicatif intégré de perte de charge par frottement

avec masse volumique variable pour le modèle de Thom (1964) et utilisant la relation (95)
Coefficients de Thom (1964) Relation (95)

1,7 4,1 8,6 14,5 20.7 1,7 4,1 8,6 14,5 20,7
Qualité [%] (qualité en sortie, en
entrée la qualité est exactement zéro)
1 1,49 1,11 1,03 1,00 1,49 1,19 1,07 1,03 1,01
5 3,71 2,09 1,31 1,10 3,45 1,93 1,36 1,14 1,04
10 6,30 3,11 1,71 1,21 1,06 5,89 2,86 1,72 1,28 1,07
20 11,40 5,08 2,47 1,46 1,12 10,79 4,72 2,44 1,57 1,15
30 16,20 7,00 3,20 1,72 1,18 15,68 6,59 3,15 1,85 1,22
40 21,00 8,80 3,89 2,01 1,26 20,57 8,45 3,87 2,14 1,29
50 25,90 10,60 4,55 2,32 1,33 25,46 10,31 4,59 2,42 1,36
60 30,50 12,40 5,25 2,62 1,41 30,36 12,17 5,31 2,71 1,44
70 35,20 14,20 6,00 2,93 1,50 35,25 14,03 6,02 2,99 1,51
80 40,10 16,00 6,75 3,23 1,58 40,14 15,90 6,74 3,28 1,58
90 45,00 17,80 7,50 3,53 1,66 45,04 17,76 7,46 3,56 1,66
100 49,93 19,65 8,17 3,83 1,74 49,93 19,62 8,18 3,85 1,73

BN 3 050v2 – 20
différence de pression due à la variation de la masse volumique évaluant le taux de vide à l’aide du modèle HEM, on obtient alors
(terme d’accélération) donne : les mêmes valeurs que celles calculées avec la relation (99). Le
tableau 10 compare les deux approches et permet d’estimer
l’impact du glissement dans l’évaluation de la variation de pres-
(98)
sion due à l’accélération.
Dès l’ébullition en masse, dans la section 1, la qualité x devient
4.3.6.3 Variation de pression due à la gravité
non nulle et la masse volumique est celle du liquide à saturation.
Avec la relation (34) du modèle HEM et considérant dans la sec- Pour ce qui concerne la prise en compte de la gravité,
tion 2 la qualité x2, on obtient : l’approche est similaire au cas à masse volumique constante
(§ 4.3.4). Tant que l’écoulement reste simple phase, on peut
prendre la masse volumique moyenne entre les deux hauteurs et
(99) le calcul de la variation de pression est direct. Pour un écoulement
diphasique avec variation de la qualité, l’approche HEM est suffi-
Le cas particulier de l’évaporation d’un fluide sur une hauteur L sante et permet d’intégrer directement grâce au volume spéci-
depuis l’état de liquide saturé vers l’état diphasique de qualité x2 fique et l’on obtient :
et de taux de vide α2 est présenté dans [5] pour le cas d’une évo-
lution linéaire de la qualité mais sans considérer l’hypothèse HEM
pour le calcul du taux de vide. La perte de charge par accélération
s’écrit :
(100)
À noter que Thom [23] propose un facteur r4 compris entre 0,1

Thom [23] propose un facteur global r2 suivant :
et 1 :
(101)
Les valeurs tabulées présentées en [23] et reprises en [5] sont cal- Le tableau 11 compare les deux approches et permet d’estimer
culées à l’aide de la relation (100), le taux de vide étant calculé à l’impact du glissement imposé à 1 (HEM) dans l’évaluation de la
l’aide de la relation (31) avec le glissement évalué selon (53). En variation de pression due à la gravité.
Tableau 10 – Valeurs du coefficient multiplicatif intégré de perte de charge par accélération

avec masse volumique variable pour les modèles avec glissement et avec l’hypothèse HEM
Modèles avec glissement Modèle HEM sans glissement

Relation (100) ou coefficients de Thom (1964) Relation (99)

1,7 4,1 8,6 14,5 20,7 1,7 4,1 8,6 14,5 20,7
Qualité [%] (qualité en sortie, en entrée
la qualité est exactement zéro)
1 0,41 0,20 0,10 0,04 0,01 0,98 0,37 0,14 0,06 0,01
5 2,17 1,04 0,49 0,22 0,07 4,89 1,86 0,72 0,28 0,07
10 4,62 2,17 1,00 0,44 0,13 9,79 3,72 1,44 0,57 0,15
20 10,39 4,68 2,10 0,91 0,27 19,57 7,45 2,87 1,14 0,29
30 17,30 7,54 3,29 1,41 0,41 29,36 11,17 4,31 1,71 0,44
40 25,37 10,75 4,58 1,94 0,55 39,14 14,90 5,74 2,28 0,58
50 34,58 14,30 5,97 2,49 0,70 48,93 18,62 7,18 2,85 0,73
60 44,93 18,21 7,45 3,07 0,85 58,71 22,34 8,61 3,42 0,88
70 56,44 22,46 9,03 3,68 1,00 68,50 26,07 10,05 3,99 1,02
80 69,09 27,06 10,79 4,31 1,16 78,28 29,79 11,48 4,55 1,17
90 82,90 32,01 12,48 5,07 1,32 88,07 33,52 12,92 5,12 1,31
100 98,10 37,30 14,34 5,66 1,48 97,86 37,24 14,35 5,69 1,46

BN 3 050v2 – 21
Tableau 11 – Valeurs du coefficient multiplicatif intégré de perte de charge par gravité avec masse volu-
mique variable pour le modèle de Thom (1964) et utilisant la relation intégrée à partir du modèle HEM(95)
Coefficients de Thom (1964) Relation HEM

1,7 4,1 8,6 14,5 20,7 1,7 4,1 8,6 14,5 20,7
Qualité [%] (qualité en sortie, en entrée
la qualité est exactement zéro)
1 0,84 0,92 0,96 0,99 0,99 0,70 0,85 0,93 0,97 0,99
5 0,55 0,70 0,82 0,91 0,97 0,36 0,56 0,75 0,88 0.97
10 0,40 0,55 0,71 0,84 0,94 0,24 0,42 0,62 0,79 0,93
20 0,27 0,40 0,56 0,72 0,89 0,15 0,29 0,47 0,67 0,88
30 0,21 0,32 0,47 0,64 0,84 0,12 0,22 0,39 0,58 0,83
40 0,17 0,27 0,41 0,58 0,80 0,09 0,19 0,33 0,52 0,79
50 0,14 0,23 0,36 0,53 0,77 0,08 0,16 0,29 0,47 0,75
60 0,12 0,20 0,32 0,48 0,73 0,07 0,14 0,26 0,43 0,72
70 0,11 0,18 0,29 0,44 0,70 0,06 0,13 0,24 0,40 0,69
80 0,10 0,16 0,27 0,41 0,68 0,06 0,12 0,22 0,38 0,66
90 0,09 0,15 0,25 0,39 0,65 0,05 0,11 0,20 0,35 0,64
100 0,08 0,14 0,23 0,36 0,63 0,05 0,10 0,19 0,33 0,62
4.3.7 Évaluation du débit

Le débit de l’écoulement est directement lié à la différence de 1,00
pression entre l’amont et l’aval. La différence de pression peut
0,90
être d’origine gravitaire (amont surélevé par rapport à l’aval) ou
apportée par une pompe. Les pertes de charge dans l’écoule- 0,80
ment vont influencer la valeur du débit de l’écoulement. Intuiti-
vement, d’importantes pertes de charge réduiront 0,70
Pression [MPa]
considérablement le débit. Les pertes de charges sont directe-

0,60
ment proportionnelles au carré de la vitesse moyenne (et donc
du débit) selon les relations générales (60) et (90). Si l’on consi- 0,50
dère pour les pertes de charges par frottement la dépendance
en Reynolds (et donc en vitesse moyenne), la variation globale 0,40
de la perte de charge en fonction du débit ou de la vitesse 0,30
moyenne suit une loi en puissance 1,80 environ. D’un autre Débit
côté, soit la charge est constante (différence de hauteur entre 0,20 maximal
l’amont et l’aval, différence de température pour le thermosi- estimé
0,10 Débit minimal estimé (débit
phon – cf. § 4.5.7.1), soit elle est apportée par la pompe. Il faut
alors connaître la caractéristique de la pompe, c’est-à-dire la (débit thermohydraulique) mécanique)
0,00
courbe de variation du débit en fonction de la charge ou hauteur 0 0,1 0,2 0,3 0,4 0,5 0,6Débit estimé 0,9 1
manométrique apportée par la pompe. aux mieux
Débit volumique [m3/s]
On obtient alors l’évolution, présentée figure 11, de la variation
de la hauteur manométrique apportée par la pompe en fonction Circuit
Pompe
du débit délivré. Sur la même figure 11, on peut tracer l’évolu-
tion de la perte de charge en fonction du débit dans le circuit. Le
débit qui s’établira dans le circuit correspondra alors à l’inter- Figure 11 – Caractéristique de pompe et du circuit
section des deux courbes. À titre illustratif, les incertitudes sur
les pertes de charges ont été intégrées dans cette même
figure 11. la puissance apportée au circuit par la pompe est égale au produit
La puissance réelle de la pompe est liée au rendement hydrau- du débit volumique et de la charge :
lique et mécanique, ces informations étant fournies par le fabri-
cant de la pompe. On peut considérer en première approche que

BN 3 050v2 – 22
4.3.8 Risque de cavitation

Tableau 12 – Vitesse du son dans l’eau
Si la pression du fluide est inférieure à sa pression de satura- (liquide ou vapeur)
tion, il y a risque de cavitation, le fluide passant à l’état vapeur et
pouvant endommager les pompes, les vannes ou tout autre com- Pression Température Vitesse
posant dans le circuit. Ce risque est évalué en comparant le NPSH État
[MPa] [°C] du son [m/s]
disponible, acronyme anglais pour Net Positive Suction Head, à la
pression nécessaire à l’entrée de la pompe (« NPSH required »). 0,1 50 Liquide 1 545
Cette dernière valeur est déterminée par le fabricant de la pompe.
Le NPSH disponible (available ; NPSHA) est déterminé par l’éva- 1 100 Liquide 1 547
luation des pertes de charge du circuit, des pertes de pression
(gravité notamment) et de la pression de saturation. Le NPSHA 2 200 Liquide 1 333
doit toujours être supérieur au NPSHR pour assurer un fonction-
nement de la pompe sans cavitation. 5 200 Liquide 1 344
5 250 Liquide 1 154

4.4 Évaluation du débit critique
5 275 Vapeur 513
En cas d’accident de perte de réfrigérant primaire, une brèche
apparaît sur une conduite à haute pression (autour de 155 bar et 5 300 Vapeur 539
entre 300 et 330 °C pour le circuit primaire lorsque le réacteur est
en puissance). Au début de l’accident, dans la phase de décom- 15,5 300 Liquide 970
pression (cf. § 5.2.2.1), l’eau est soumise à une différence de pres-
sion entre l’amont et l’aval, qui est autour de la pression

atmosphérique.
rapport à sa vitesse dans la conduite et à la brèche. Dans ce cas,
À partir de la relation de Bernoulli, on peut établir la relation
si les conditions amont (dans la cuve ou le réservoir par exemple)
entre la différence de pression et la vitesse moyenne débitante
sont indicées « 0 », on a :
suivant :
(102)
avec C une constante de débit inférieure à un, liée à la perte de La conservation de l’énergie (Bernoulli) et l’expression de l’enthal-
charge. Cette relation n’est valable qu’en cas d’écoulement consi- pie d’un gaz parfait en fonction de la chaleur spécifique à pression
déré comme incompressible, si la vitesse du fluide est faible constante donnent :
devant la vitesse du son.
La vitesse du son « c » se calcule selon la relation (103) à partir
de la dérivée partielle de la pression par rapport à la masse volu-
mique à entropie constante
Le cp s’exprime selon :
(103)
Le tableau 12 regroupe quelques valeurs de la vitesse du son avec Reau la constante universelle des gaz parfaits rapportée à
dans l’eau, liquide ou vapeur. Ces valeurs sont calculées à partir l’eau (0,4615 kJ/kg K).
de la dérivée partielle de la pression par rapport à la masse volu-
mique à entropie constante (relation (103)). On a donc, d’une part venant de la conservation de l’énergie :
Lors de la détente du fluide, le débit à la brèche est limité par

l’atteinte de la vitesse sonique, rarement en phase liquide du fait
de faibles vitesses débitante et d’une vitesse du son élevée. La
chute de pression provoque aussi le passage potentiel de d’une et d’autre part, lorsque la vitesse atteint la vitesse du son (relation
partie du fluide en vapeur. En prenant pour calculer la vitesse (103)) :
débitante, un coefficient de débit égal à 1 et la pression atmosphé-
rique en aval, à 5 MPa pour un fluide à 300 °C (vapeur), l’expres-
sion (102) conduit à une vitesse de 673 m/s supérieure à la vitesse
sonique (539 m/s, cf. tableau 12). Dans ce cas, le débit est maxi-
mum lors de l’atteinte de la vitesse du son, reste contant et ne soit
dépend plus des conditions aval. C’est le débit critique. Lorsque le
fluide est diphasique, le calcul de la vitesse sonique n’a plus beau-
coup de sens. En pratique, on cherche directement l’expression du
débit critique.
La vitesse massique critique correspondante est le produit de la

4.4.1 Débit critique simple phase vapeur vitesse et de la masse volumique au niveau de l’atteinte de la
L’expression du débit critique en phase vapeur seule peut se vitesse sonique :
déterminer avec l’hypothèse d’une détente adiabatique de gaz
parfait. Les conditions de référence sont considérées en amont,
avec l’hypothèse que le mouvement du fluide est négligeable par

BN 3 050v2 – 23
soit encore on obtient :
(107)
On calculera la masse volumique à la brèche avec l’hypothèse
HEM selon la relation (34) et l’enthalpie à la brèche selon (23). La
proportion de vapeur est déduite de l’hypothèse de la détente
isentropique.
et finalement
(104) Les entropies à saturation sont considérées à la pression au

niveau de la brèche. Cette pression n’est pas connue, il faut donc
itérer sur des valeurs estimées de cette pression jusqu’à obtenir le
débit maximal.
4.4.2 Débit critique en diphasique Enfin, on peut utiliser le modèle de débit à l’équilibre (ERM pour
Très rapidement, l’eau liquide passe tout ou partie en vapeur à « equilibrium rate model ») de Fauske selon :
la brèche. Par nature on est en situation de déséquilibre et la pro-
portion de vapeur est difficile à estimer. On cherche alors à éva- (108)
luer directement le débit à la brèche sans déterminer la vitesse du c p, l
son. Selon les hypothèses simplificatrices relatives aux déséqui-
libres thermique (températures des phases vapeur et liquide) et/ avec hlv chaleur latente (J/kg) à P0,
ou dynamique (vitesses), diverses corrélations plus ou moins υlv différence des volumes spécifiques (m3/kg) vapeur et
complexes sont proposées. liquide à saturation à P0,
Lorsque en amont (dans le réservoir), le fluide est saturé dipha- Tsat température de saturation (K) à P0,
sique, une première approche, ou « equilibrium frozen flow
(EFF) », consiste à considérer que la proportion de vapeur à la cp, l chaleur spécifique (J/kg/K) à saturation à P0.
brèche est identique à celle dans le réservoir. On démontre [25] Cette relation peut estimer le débit critique lorsque le fluide est
une relation proche de la relation (104) mais valable pour un écou- saturé en amont. La relation suivante permet de tenir compte de
lement diphasique selon : la sous-saturation :
(109)
0 0
(105)
avec Psat la pression de saturation (Pa) correspond à la
température en amont T0,
Cette relation surestime largement le débit critique lorsque le titre
massique est faible (on s’éloigne du gaz parfait), la référence [25] ρl masse volumique (kg/m3) en amont à la pression ρ0 et
cite une relation proposée par divers auteurs (Starkman, Henry, la température T0.
Fauske) : Le tableau 13 propose les relations à utiliser selon les condi-
tions du fluide dans le réservoir.
(106) La corrélation de Gros d’Aillon utilisée dans le code CATHARE
[29] est valable pour toute condition du fluide dans le réservoir.
Elle comprend de plus un terme lié au déséquilibre dans l’écoule-
avec v volume spécifique pondéré du mélange tel que : ment et à la brèche lié au rapport L/D de la distance L parcourue
par le fluide dans la conduite entre le « réservoir » au diamètre D
de la brèche. Il est possible de considérer la perte de charge sin-
gulière à la brèche (K). La corrélation s’écrit :
Le fluide étant saturé diphasique, une seconde approche utilise le (110)

modèle à l’équilibre homogène (HEM) (voir § 3.2). En utilisant la
conservation d’énergie (Bernoulli) :
En première approximation le rapport L/D et K peuvent être pris
nuls.
Pour des conditions en amont sous-saturées :
Tableau 13 – Relations à utiliser selon les conditions du fluide dans le réservoir
Conditions du Saturé liquide Saturé titre très Titre faible

Sous-saturé Titre fort Saturé vapeur
fluide en amont titre nul faible à moyen
Relations (107) (109) (110) (107) (108) (110) (106) (107) (110) (105) (106) (110) (105) (106) (110) (104) (105)
valables (106) (110)
Commentaire (107) sous évalue le débit et demande

des itérations

BN 3 050v2 – 24
• ; • ΔP = 0,1 P0 (pression en Pa).

• (conditions amont) ; Exemple 10
• (les pressions doivent être exprimées Pour des conditions amont sous-saturées, les diverses relations
en Pa). proposées donnent les valeurs du débit critique données dans le
tableau 14.
La relation (110) se simplifie alors selon :
Le débit critique HEM résulte d’itérations sur la pression estimée à
la brèche, le tableau 14 regroupe les valeurs à convergence.
À saturation, on obtient les valeurs données en fonction du titre
thermodynamique à l’équilibre regroupées dans le tableau 15.
Pour des conditions en amont saturées :
Pour de la vapeur saturée ou surchauffée, les valeurs de débit cri-
• en première approximation, on utilise les relations HEM pour
tique calculé selon les corrélations applicables sont rassemblées dans
le calcul de la masse volumique et on utilise le titre thermo-
le tableau 16.
dynamique à l’équilibre χ ;
• avec et α calculé avec la Les valeurs proposées dans cet exemple montrent une grande
variabilité du débit critique. Une source significative d’incertitude
relation (33) ; supplémentaire réside dans les grandeurs estimées en amont
• ρm calculé avec la relation (34) (inverse du volume spéci- (pression, titre), d’autant qu’en cas de brèche sur le circuit, ces
fique) ; valeurs évoluent rapidement. De par son application possible pour
toute condition au réservoir, la corrélation de Gros d’Aillon donne
des ordres de grandeur satisfaisants.
Tableau 14 – Débits critiques d’après diverses relations pour des conditions amont sous-saturées
Conditions amont Pression à la brèche

(pression 5 MPa (HEM) Débit critique [kg/m2.s]
263,94 °C) [MPa]
P0 [MPa] T0 [°C] (107) HEM (108) (109) ERM (110) GDA
5 240 3,32 53 410 58 125 56 929
5 260 4,63 25 145 34 147 34 908
5 263 4,21 22 090 28 284 29 734
5 263,90 4,16 21 658 26 239 27 968
5 Tsat 4,16 21 639 26 140 27 883
Tableau 15 – Débits critiques d’après diverses relations à saturation
Conditions
Pression à la
amont (pression
brèche (HEM) Débit critique [kg/m2.s]
5 MPa
[MPa]
263,94 °C)
χ (105) EFF (106) SHF (107) HEM (110) GDA
0,00 4,15 21 639 27 881
0,02 3,9 58 311 28 571 19 257 24 693
0,10 3,5 26 078 21 900 15 329 18 093
0,50 3,0 11 662 11 425 9 394 9 936
0,90 2,9 8 693 8 676 7 405 7 624
0,99 2,9 8 288 8 290 7 110 7 295

BN 3 050v2 – 25
Tableau 16 – Débits critiques d’après diverses relations pour de la vapeur saturée ou surchauffée
Pression à la brèche
Conditions amont
(HEM) Débit critique [kg/m2.s]
(pression 5 MPa)
[MPa]
χ (104) Gaz parfait (105) EFF (107) HEM (110) GDA
1,00 2,92 8 250 8 250 7 079 7 260
1,1 7 866 6 933
4.5 Transfert thermique L est une échelle de longueur caractéristique, en écoulement en

conduite circulaire on prend le diamètre de la conduite. Dans un
Le transfert de la chaleur s’opère selon trois types de proces- écoulement interne de conduite de toute géométrie, le diamètre
sus, la conduction (au sein d’un même matériau), la convection hydraulique (relation (1)) ou le diamètre chauffant (pour lequel le
(impliquant un fluide en mouvement) et le rayonnement. Le périmètre mouillé est remplacé par la longueur sur laquelle les
paramètre dominant pour la conduction est la conductivité du échanges thermiques peuvent se produire) peuvent être utilisés.
Cependant on introduit de nouvelles sources possibles d’incerti-
matériau et pour le rayonnement la température du corps. Le
tudes et d’erreurs liés aux relations utilisées pour déterminer le
rayonnement n’est pas un acteur significatif dans les problèmes
diamètre équivalent. Dans les relations proposées dans la suite
de thermohydraulique et n’est pas abordé dans cet article. La

les coefficients empiriques peuvent alors varier si la géométrie
convection relève plus particulièrement du domaine de la ther- prise en compte n’est plus une conduite circulaire.
mohydraulique puisqu’intervient la vitesse du fluide dans l’éva-
luation du transfert thermique. Intuitivement, dans la Le nombre de Nusselt indique le rapport de la quantité de chaleur
convection, plus la vitesse du fluide est importante, meilleur échangée par convection à celle échangée par conduction. L’analyse
sera le transfert thermique. Les bases théoriques peuvent être dimensionnelle permet d’identifier les nombres adimensionnels qui
consultées dans [BE8205]. On indique ici les principales sont en relation avec le nombre de Nusselt dans les problèmes de
méthodes applicables aux écoulements en conduite. On consi- transfert thermique. Dans cette analyse, on distingue deux types de
convection selon l’origine de la mise en mouvement du fluide. Soit
dérera tout d’abord les écoulements simple phase puis ensuite
la vitesse du fluide est imposée (par exemple par une pompe) et la
les écoulements diphasiques. Les relations proposées seront
convection est forcée, soit le fluide est mis en mouvement naturelle-
toutes relatives à celle d’un écoulement établi, des corrélations
ment compte tenu des différences de températures au sein du fluide
particulières (cf. [BE8206]) pouvant être utilisées pour prendre qui engendrent les forces de flottabilité relatives aux différences de
en compte les effets en entrée d’une conduite. masse volumique au sein du fluide. Ce dernier aspect (convection
Dans les situations de transfert thermique, le type de conditions naturelle) est détaillé au paragraphe 4.5.7.
aux limites dicte les phénomènes et paramètres. Si l’on impose Pour la convection forcée, les nombres adimensionnels identi-
les températures, un bon transfert thermique conduit à un fiés sont les nombres de Reynolds, de Prandtl et de Peclet.
échange (flux de chaleur) important. Si l’on impose le flux de cha-
leur (par une puissance générée par exemple), alors un bon trans- Le nombre de Reynolds :
fert thermique conduit à de faibles différences de températures.
Par températures, on entend la température « moyenne » du (113)
fluide et la température en paroi. En effet on se place dans une
approximation à une dimension et l’objectif est de calculer la tem- est le rapport des forces d’inertie ρ U2/2 à celles dues à la visco-
pérature de paroi. Typiquement, une situation à températures sité μU/D.
imposées correspond au générateur de vapeur ou au condenseur ; Le nombre de Prandtl :
une situation à flux de chaleur imposé correspond au cœur du
réacteur. Cependant, quel que soit le type de conditions aux
limites, le transfert thermique par convection répond à la loi de (114)
Newton qui s’écrit :
donne le rapport de la diffusivité de la quantité de mouvement ν à
(111) la diffusivité thermique λ/ρcp.
Le nombre de Peclet est le produit du nombre de Reynolds et
avec φ (W/m2) le flux de chaleur surfacique entre la paroi solide et
de Prandtl :
le fluide et hc (W/m2 °C) le coefficient de transfert thermique par
convection. On présente dans la suite diverses méthodes pour
(115)
évaluer le coefficient de transfert thermique selon le type d’écou-
lement.
Ce nombre joue le même rôle dans l’équation de conservation de
l’énergie que le nombre de Reynolds dans la conservation la
4.5.1 Écoulement simple phase quantité de mouvement. Le nombre de Péclet est le rapport des
transferts de chaleur par convection et par conduction.
Les méthodes les plus efficaces pour évaluer le transfert ther-
mique par convection s’appuient sur l’analyse dimensionnelle. Le En écoulement laminaire (cf. 4.3) dans un écoulement déve-
loppé dans l’hydraulique et la thermique, des solutions analy-
nombre de Nusselt fait intervenir directement le coefficient de
tiques sont possibles et l’on trouve en écoulement dans une
transfert selon :
conduite circulaire :
– à températures imposées, Nu = 3,66 ;
(112)
– à flux de chaleur imposé, Nu = 48/11 = 4,36.

BN 3 050v2 – 26
En écoulement turbulent, pour des écoulements simple phase la température de film, moyenne de la température moyenne et
en conduite circulaire de diamètre D, la relation de Dittus-Boelter de la température en paroi qu’il faut évaluer. Avec des conditions
(telle qu’introduite par McAdams [38]) est la plus souvent utilisée : limites imposées sur le flux de chaleur, des itérations sont néces-
saires jusqu’à convergence vers la valeur finale de la température
de paroi.
(116)
Cependant, dès que la température en paroi atteint la tempéra-
ture de saturation, l’ébullition locale est proche et les relations du
Les caractéristiques physiques sont prises à la température type de celle de Dittus-Boelter sont alors à utiliser avec précaution
moyenne du fluide qui peut être estimée par un bilan thermique. et on préférera des formulations prenant en compte l’ébullition
La vitesse de référence est la vitesse débitante (vitesse moyenne), locale (§ 4.5.2).
il est préférable d’utiliser la vitesse massique G qui est constante
dans la conduite. Cette équation n’est valable que si :
– le nombre de Prandtl est limité (entre 0,7 et 100) ; 4.5.2 Écoulement avec ébullition locale
– le fluide est chauffé (si ce n’est pas le cas – fluide refroidi –, la
puissance du Pr est de 0,3 au lieu de 0,4) ; Lorsque la température de paroi approche la température de
– l’écoulement est en conduite (ou assimilée) lisse (si ce n’est saturation (déterminée par la connaissance de la pression locale),
pas le cas, le facteur de proportionnalité 0,023 peut aller vers des un changement de phase peut se produire. Dans le cas d’un fluide
valeurs autour de 0,055) ; chauffé par transfert de chaleur à la paroi, une légère surchauffe
– l’écoulement est établi, permanent (en entrée d’une conduite, en paroi ΔTp est nécessaire avant ébullition. Pour évaluer la sur-
des relations spécifiques doivent être utilisées voir par exemple chauffe en paroi, les relations les plus utilisées sont celles de :
[21] et [BE8206]) ; • Jens et Lottes :
– il n’y a pas de fortes variations de températures dans le fluide. Le
cas échéant, on peut utiliser le facteur correctif de Seider-Tate et alors (121)
(117) • Thom :
μp est la viscosité du fluide considérée à la température de paroi, 22,65 (122)

qui est l’inconnue. Il faut donc itérer sur des estimations de cette
température jusqu’à convergence. avec T température (en °C),
La relation de Dittus-Boelter peut aussi s’écrire en utilisant le
p pression (en MPa),
nombre de Stanton :
φ flux à la paroi (en MW/m2).
(118) Il n’y a pas de dépendance à la vitesse du fluide comme en
convection forcée car c’est le transfert thermique par évaporation
Le nombre de Stanton exprime le rapport entre le transfert de cha- du liquide qui prédomine, ce qui explique la dépendance en pres-
leur par convection au transfert de chaleur évalué selon un bilan sion relative à la chaleur latente décroissante avec la pression.
thermique (relation (5)). La relation de Dittus-Boelter s’écrit alors : L’équation de Thom est préférable à celle de Jens et Lottes bien
que proche sur le domaine de validité en pression (3,4 à 13,8 MPa)
(119) car elle semble valable jusqu’à 0,4 MPa. L’incertitude de la rela-
tion de Thom est estimée à 15 %.
On remarque que, pour un nombre de Prandtl égal à l’unité, cette
équation est similaire à la relation de MacAdams (69) qui lie le coeffi- La connaissance du flux φ permet de déduire ΔT et ainsi h sui-
cient de frottement au nombre de Reynolds et un facteur 8 relie ce vant la loi de Newton (relation (111)).
coefficient au nombre de Stanton. Ceci indique ainsi des processus
analogues entre le transport de la quantité de mouvement et celui de
la chaleur (analogie de Reynolds). On démontre d’ailleurs exactement
4.5.3 Écoulement saturé
la relation entre le nombre de Stanton et le coefficient de trainée : Lorsque l’écoulement est saturé, jusqu’à des valeurs de qualité
pas trop élevées, les corrélations de Jens et Lottes (121) et Thom
(120) (122) peuvent encore être utilisées. Une autre méthode consiste à
considérer le coefficient de transfert de chaleur comme la somme
L’incertitude sur le nombre de Nusselt évalué à l’aide de la rela- d’un terme dû à l’ébullition nucléée et d’un autre terme dû à la
tion de Dittus-Boelter peut atteindre 25 %, d’autres corrélations convection. Le premier terme est fonction du nombre d’ébullition,
permettent de réduire cette incertitude à 10 % par exemple la cor- rapport du flux au produit de la vitesse massique et de la chaleur
rélation de Gnielinski : latente φ/Ghlv, le second terme est fonction du paramètre de Mar-
tinelli Xtt (81) ou (82), reprenant le principe des fluides séparés
Nu présenté au paragraphe 4.3.2.2. Le coefficient de transfert suit une
relation du type :
f est le coefficient de frottement de Darcy qui peut être estimé à

l’aide d’une relation complémentaire (123)
Les coefficients ont des valeurs ajustées différentes selon divers

auteurs et sont indiquées dans le tableau 17. Le coefficient de
transfert relatif au liquide seul est calculé avec la relation clas-
ou encore avec les relations de Blasius (68) ou McAdams (69). sique de Dittus-Boelter (116), les propriétés du liquide considérées
à saturation.
D’après [40], les propriétés physiques de l’eau sont considérées
à la température moyenne mais d’autres sources, avec des formu- Suivant le même principe de décomposition, la corrélation de
lations quelquefois un peu différentes, préconisent de considérer Chen est préférable car plus proche des mesures expérimentales.

BN 3 050v2 – 27
4.5.4 Transfert de chaleur post-assèchement

Tableau 17 – Valeurs des coefficients du transfert
de chaleur en écoulement saturé utilisant Des situations sont rencontrées dans lesquelles un film de
la relation (123) vapeur isole la paroi chauffante du liquide, en cas de crise d’ébul-
b lition (cf. § 5.1). Dans ce cas le transfert de chaleur se dégrade for-
a2 tement, la vapeur d’eau ayant des propriétés isolantes par rapport
Auteurs a1 ébullition puissance
convection au liquide. La température en paroi est trop importante pour
de Xtt
qu’elle puisse être ‹ remouillée › par le liquide encore présent dans
Dengler et Valeur de a1 nulle ; la 3,5 –0,5 l’écoulement. L’écoulement est dit ‹ annulaire inversé ›, annulaire
Addoms contribution due à car les deux phases sont séparées, inversé car c’est la vapeur qui
l’ébullition est négligée est en paroi et le liquide au centre de l’écoulement. Une autre
Bennett et al. étant considérée 2,9 –2/3 situation de transfert thermique fortement dégradé apparaît
comme non lorsque le film liquide en paroi est totalement asséché ou
significative devant la ‹ entraîné › (arraché) par la vapeur et le liquide se retrouve en
convection gouttes dispersées dans la vapeur. On peut alors estimer le coeffi-
cient de transfert de chaleur soit en supposant l’écoulement tota-
Schrock et 170 1,11 –2/3
7 391 lement en vapeur qui devient surchauffée (et dans ce cas la
Grossman.
0,023 relation de Dittus-Boelter (116) s’applique), soit l’on suppose que
des petites gouttelettes de liquide subsistent dans l’écoulement et
Collier et 6 700 2,34 –2/3 peuvent participer au transfert de chaleur, enfin le transfert de
Pulling chaleur par rayonnement peut être significatif. Il apparaît cepen-
dant que l’écoulement est sujet à des forts déséquilibres ther-
miques, la vapeur étant surchauffée et l’évaporation des
Les modèles pour les contributions dues à l’ébullition et la gouttelettes d’eau subsistant dans l’écoulement retardé. Du
convection sont cependant plus complexes. liquide subsiste (la proportion en masse de vapeur n’est alors pas
La part d’ébullition nucléée se base sur l’équation de Forster- de 100 %) tandis que le titre thermodynamique à l’équilibre peut
Zuber : être largement supérieur à l’unité. On peut estimer la proportion
réelle en masse vapeur xa à partir du titre thermodynamique à
l’équilibre χ – relation (6) – selon :
avec ΔT la surchauffe en paroi soit Tp –

Tsat (à itérer si la température Une relation empirique [39] permet de calculer le rapport de la
en paroi est inconnue),
qualité réelle au titre à l’équilibre selon :
différence, en Pa, entre les
pressions de saturation
correspondant, respectivement
à la température de paroi et à la
température moyenne du fluide
(qui est à saturation). Les enthalpies à saturation sont calculées à partir de la pression
Les unités à utiliser sont les unités SI, notamment chaleur locale, Pc est la pression critique (22,1 MPa pour l’eau), la gran-
latente en J/kg, pressions en Pa et chaleur spécifique en J/kg/K. deur Td exprime le degré normalisé de surchauffe de la vapeur
La part convective se base sur l’équation de Dittus-Boelter (116) :
Pour des valeurs de titre thermodynamique à l’équilibre ne dépas-

sant pas 0,9, l’estimation du coefficient de transfert de chaleur
peut se faire indépendamment de l’estimation de la surchauffe de
Le coefficient total est la combinaison des deux contributions
selon la vapeur. La corrélation de Groeneveld [26] est applicable en tube
et en espace annulaire, avec des jeux de coefficients différents
rassemblés dans la tableau 18 :
avec (124)
Le nombre de Reynolds ‹ homogène › est calculé en prenant en

compte le taux de vide homogène αh – relation (33) – et le titre
massique à l’équilibre χ selon :
(la définition de Xtt est donnée en (81)),
Le nombre de Prandtl est calculé en prenant les propriétés du

et fluide en phase vapeur à la paroi :

BN 3 050v2 – 28
tion), les propriétés en phase vapeur y compris la température

Tableau 18 – Valeurs des coefficients sont prises à la température de saturation – imposée par la pres-
du transfert de chaleur en écoulement saturé sion ambiante (la chaleur latente est évaluée à cette pression). Les
utilisant la relation (124) nombres de Jakob Ja et de Rayleigh Ra sont aussi rappelés en fin
de l’article.
Erreur Le débit-masse de condensat par unité de longueur de tube Γ
a b c d affichée (en kg/s par m) peut être calculé à partir du bilan thermique en
([26]) considérant la chaleur latente « corrigée » h’lv. On obtient alors la
relation suivante :
Tube 0,00109 0,989 1,41 – 1,15 11,5 %
. ’ ρ λ3 (ρ – ρ ) 1/4
Espace 0,0520 0,688 1,26 – 1,06 6,9 % m πD ghlυ ι ι ι υ
Γ = = 0,73
annulaire L hlυ [1 + 0,68Ja] μιD(Tυ – Tp)
Ensemble 0,00327 0,901 1,32 – 1,50 12,4 % Dans l’enceinte, il est nécessaire d’évaluer le taux de condensa-
tion qui dicte l’évolution de la pression. Pour des structures verti-
cales ou inclinées, les relations doivent être modifiées. La
Le terme Y est un facteur de correction pour prendre en compte la condensation peut aussi s’effectuer par formation de gouttes, pro-
présence simultanée de liquide et vapeur. Ce terme s’écrit en cessus plus efficace que la condensation en film. En outre, la pré-
fonction du titre à l’équilibre χ : sence de gaz incondensables (air, hydrogène) perturbe
sensiblement la condensation de la vapeur. On se reportera aux
ouvrages généraux de référence [5], [6] ou [7] (et [BE8238]) qui
décrivent l’ensemble des divers modes de condensation.
La référence [26] précise le domaine de validité et l’erreur relative 4.5.6 Profil de température dans l’assemblage
entre les essais et la relation proposée, cette même référence pré- combustible
sente aussi de nombreuses autres corrélations.
Le lecteur trouvera en référence [22] un grand nombre d’infor-
Pour des titres (à l’équilibre) approchant ou dépassant l’unité,
mations sur le combustible nucléaire. Le combustible des REP est
diverses approches, dont des tables fonction de la pression,
constitué de pastilles d’oxyde d’uranium (UO2) et, éventuellement
vitesse massique, titre à l’équilibre et surchauffe, sont disponibles
dans le cas d’oxyde mixte (MOX) de quelques pourcents d’oxyde
[27]. Ces relations doivent être appliquées avec précaution aux
de plutonium (PuO2). Ces pastilles, de forme cylindrique à section
écoulements dans le cœur d’un réacteur à eau pressurisée compte
circulaire, sont placées dans des gaines en Zircaloy de diamètre
tenu de la complexité de ce type d’essais, des moyens de mesures
intérieur légèrement supérieur au diamètre extérieur des pastilles.
associés encore insuffisants et de la géométrie en tube très éloi-
À la fabrication, la gaine est remplie d’hélium qui vient occuper
gnée des conditions réelles.
l’espace compris entre pastilles et gaine. Au cours de l’irradiation,
des produits de fission gazeux migrent dans cet espace.
4.5.5 Transfert de chaleur par condensation L’ensemble de l’empilement des pastilles constitue le crayon,
d’une hauteur de 3,66 ou 4,27 m selon le type de réacteur. Les
La vapeur au contact d’une paroi froide va se condenser. Dans crayons sont disposés de façon uniforme pour constituer l’assem-
ce cas, pour les applications nucléaires, la contrainte n’est pas sur blage de combustible, avec les tubes-guides de grappes de com-
la température car l’écoulement cède de la chaleur à l’extérieur. mande et d’instrumentation. Dans les réacteurs occidentaux, la
On cherche à estimer le taux de condensation, notamment dans la structure de l’assemblage est carrée avec un ensemble de 14 à 18
phase de renoyage (cf. § 5.2.2.3). La vapeur se condense sur les crayons ou tubes-guides par coté, en France tous les assemblages
structures du plénum supérieur, l’eau liquide pouvant ensuite par sont des assemblages 17x17. Dans les REP de conception sovié-
gravité rejoindre le cœur et participer au refroidissement du com- tique dits VVER (cf. [BN3101]), la structure est hexagonale et les
bustible. crayons sont disposés en pas triangulaire.
La quantité de vapeur condensée est liée directement au trans- On peut également trouver des pastilles creuses, ce qui réserve
fert de chaleur au travers de la puissance échangée relative à la un volume d’expansion pour les produits de fission et permet de
chaleur latente selon : diminuer la température maximale de la pastille.
Dans le crayon combustible, le principal mode de transmission
de la chaleur est la conduction, au sein de la pastille, dans le jeu
pastille/gaine et dans la gaine. Les autres modes de transmission
ne sont pas considérés ici (cf. [BN3060]).
Sur des structures horizontales cylindriques (diamètre D), pour
de la condensation en film laminaire, le nombre de Nusselt inté- La conduction dans un solide suit la loi de Fourier selon
gré sur tout le diamètre et la longueur peut s’écrire [5] selon : laquelle :
(125)
Le gradient thermique étant donné (différence de température
entre deux points du matériau), c’est la conduction thermique λ
avec qui limite la transmission de la chaleur. À différence de tempéra-
ture donnée entre deux surfaces, le flux est bien plus important si
la conductivité thermique du matériau est élevée (matériau
conducteur) que si le matériau est isolant. La dissipation ther-
mique n’est fonction que du solide considéré au travers de sa
conductivité qui peut varier en fonction de la température, tous
Les propriétés du liquide sont prises à la température de film les autres paramètres (pression, composition du crayon en entier,
(moyenne de la température en paroi et la température de satura- vitesse externe…) ont une influence qui peut être négligée. Le

BN 3 050v2 – 29
tableau 19 récapitule les principales caractéristiques physiques rant de la conduction radiale pure sans effets azimutaux et à une
des composants du crayon combustible. hauteur donnée)
4.5.6.1 Pastille de combustible

Nous allons considérer le problème général d’une pastille ayant
la géométrie de la figure 12. Le centre de la pastille est vide sur ce qui intégré donne :
un rayon r0. Pour les pastilles pleines, les plus courantes, r0 est
donc nul.
Environ 92 % de la puissance thermique du réacteur provient de
l’énergie de fission générée dans la pastille (les 8 % restant sont En r = r0, le flux de chaleur est nul et la loi de Fourier impose alors :
générés directement dans le modérateur et les matériaux de struc-
ture). Nous supposerons toutefois que les pastilles sont le siège
de la totalité de la puissance générée. Considérons un élément de
crayon de hauteur unitaire. En régime permanent, production ce qui conduit à la détermination de la constante selon :
(puissance volumique q’’’) et dissipation sont égales. La dissipa-
tion peut, en conduction radiale pure (absence d’effets azimu-
taux), s’exprimer par la divergence du flux de chaleur. En
considérant la divergence en coordonnées cylindriques on écrit :
La production q’’’ est égale à la dissipation (divergence du flux), et

la loi de Fourier (relation (125)) permet d’écrire (toujours considé-
Si r0 = 0 et pour une conductivité constante, alors
Tableau 19 – Ordre de grandeur des principales

caractéristiques physiques des composants La différence maximale de température dans la pastille est obtenu
du crayon combustible (neuf) au bord où r = rp.
Le produit de la puissance volumique à la section donne la puis-
Zircaloy-4 UO2 Oxyde mixte sance linéique (proportion de puissance générée par le crayon par
Matériau
Copy [57] [58] unité de longueur), ainsi :
Masse 6 600 11 000 11 500

volumique à
20 °C [kg/m3] Le maximum de température est fonction uniquement de la puis-
sance linéique et donc indépendant du rayon de la pastille à puis-
Température de 1 850 2 850 2 425 sance de réacteur donné (la puissance linéique est calculée
fusion [°C] suivant la hauteur du crayon, cf. exemple 2).
Conductivité 100 °C : 14,0 500 °C : 4,2 500 °C : 3,5 Lorsque les variations de λ avec la température sont significatives,
thermique ce qui est le cas dans la pastille, pour résoudre l’équation en régime
[W/(m · K)] 430 °C : 17,0 2 000 °C : 2,1 2 000 °C : 2,0 permanent avec λ = f (T), on doit faire appel à l’intégrale de conduc-
tivité. Une estimation est possible en utilisant la figure 13.
Capacité 100 °C : 0,30 730 °C : 0,32 730 °C : 0,34
thermique
[kJ/(kg · K)] 430 °C : 0,33 2 230 °C : 2 230 °C :
0,50 0,39 Intégrale
Intégrale [kW/m] de conductivité
λ [W / (m ∙ k)] [kW/m] à partir de 0 °C
Conductivité
8
thermique
[W / (m ∙ K)]
6 Intégrale
de conductivité
[kW/m] à partir
4 de 500 °C
r0
2
r1
0
rp 0 1 000 2 000 3 000
Température [°C]
Figure 13 – Conductivité thermique et intégrale de conductivité de

Figure 12 – Géométrie d’une pastille de combustible l’UO2 (95 % densité)

BN 3 050v2 – 30
4.5.6.2 Jeu pastille-gaine et gaine

Il n’y a pas de terme source ni dans la gaine ni dans le jeu pas- 3 000
tille-gaine. On peut appliquer directement la loi de Fourier (rela- UO2 fusion Jeu
TEMPÉRATURE (°C)
tion (125)). Alors,
2 500 Gaine
Protections 2 000
(590 W ∙ cm–1)
et par définition de la puissance linéique et du flux de chaleur (à
puissance linéique donnée, le flux diminue à mesure que la sur- 1 500
Point chaud
face au travers de laquelle passe le flux augmente) :
(420 W ∙ cm–1)
1 000
Moyenne
(186 W ∙ cm–1) 500
Alors la résolution de l’équation différentielle donne : Tgaine = 340 °C
Pastille
0
0,0
0,5
1,0
1,5
2,0
2,5
3,0
3,5
4,0
4,5
5,0
Rayon (mm)
Dans l’épaisseur du jeu, la conductivité thermique peut être prise
constante et vaut λjeu, la valeur de r2 est égale au rayon de la pas-
tille rp et de l’épaisseur e du jeu et r1 est égal à rp, alors : Figure 14 – Profils de température typique dans un crayon
plus élevée que le réacteur. Certaines conceptions récentes de

réacteurs [30] reposent sur le principe du thermosiphon. Soit
majoritairement dans certains types de réacteurs modulaires de
(on considère que e/rp << 1 et le logarithme peut se simplifier). petites puissance (SMR pour small modular reactor), soit comme
L’hélium présent dans le jeu contribue à une meilleure conduc- système de sûreté additionnel comme dans le réacteur AP1000 de
tion que l’air et évite la corrosion. À noter que le jeu est trop faible Westinghouse ou l’AES 2006 de Rosatom (système SPOT passif
pour qu’un mouvement de convection naturelle s’établisse. d’évacuation de la puissance résiduelle par les générateurs de
Dans la gaine d’épaisseur E et de rayon externe R : vapeur, voir [BN3101]).
La valeur du débit se stabilise à l’équilibre entre la force motrice
crée par le différentiel de masse volumique (compte tenu des tem-
pératures) entre le cœur (chaud) et les générateurs de vapeur
(froid) et la perte de charge totale du circuit, frottement et singula-
rités. D’une part
(on considère que E/(R-E) << 1)
Ces relations sont des relations approchées mais qui donnent
de bons ordres de grandeur.
et, d’autre part avec la relation de Darcy-Weisbach (90) dans
4.5.6.3 Transfert thermique entre le crayon et le fluide et laquelle on met dans le coefficient de perte de charge global K
profil général l’ensemble des pertes de charge par frottement et dues aux singu-
larités, toutes rapportées à la masse volumique de référence ρ0 et
Le transfert thermique entre la gaine du crayon et le fluide est
à la vitesse massique moyenne :
en fonctionnement normal principalement en convection forcée
simple phase (§ 4.5.1).
Le profil général (figure 14, cf. [BN 3060]) est fortement dépen-
dant de la valeur de la puissance linéique, le gradient de tempéra-
ture principal étant au sein de la pastille. La perte de charge résistante est égale à la force motrice et, en
utilisant la relation de Boussinesq qui linéarise la variation de
masse volumique :
4.5.7 Circulation naturelle
La circulation naturelle se met en place dès qu’un fluide est sou-
mis à des différences de température telles que le fluide chaud se
retrouve à une altitude plus basse que le fluide froid. Dans un cir- alors
cuit fermé, le fluide s’écoule alors et l’on parle de thermosiphon.
Dans un espace ouvert, l’écoulement du fluide va s’organiser sous
forme de mouvements de rotation et l’on parle de convection
naturelle. Les outils d’analyse pour ces deux situations sont détail-
lés dans la suite. La vitesse massique moyenne à l’équilibre pour le thermosiphon
est proportionnelle à la racine de la différence de température
entre les points chaud et froid :
4.5.7.1 Thermosiphon
La situation de thermosiphon se rencontre typiquement en cas
de perte des pompes primaires. Le réacteur est arrêté mais la
puissance résiduelle continue de chauffer le fluide dans le cœur.
L’eau chauffée s’élève alors, passe dans les branches chaudes et De plus, l’utilisation du bilan enthalpique (relation (5)) et l’approxi-
se refroidit dans les générateurs de vapeur situés à une altitude mation de l’élévation d’enthalpie par le produit de la différence de

BN 3 050v2 – 31
température et du cp (3) permet de relier le ΔT et la puissance Π L’approche en convection naturelle est semblable à celle de la
dissipée dans le fluide et on obtient : convection forcée. Elle repose sur l’établissement de corrélations
reliant les nombres adimensionnels typiques de Nusselt, Prandtl
et Grashof. Le nombre de Grashof a le même rôle que le nombre
de Reynolds pour la convection forcée. Le nombre de Grashof
représente le rapport des forces d’Archimède (variation de la
masse volumique) aux forces dues à la viscosité et s’écrit :
Le débit massique à l’équilibre pour le thermosiphon est propor-
tionnel à la racine cubique de la puissance résiduelle dans le
cœur. (127)
La grandeur β, coefficient de dilatation, utilisée dans la relation
de Boussinesq, est calculée selon : Si l’on est en écoulement diphasique (liquide dans la partie froide
et vapeur pour la partie chaude) le produit βΔT peut s’écrire Δρ/ρ
et :
(126)
pour laquelle quelques valeurs relatives à l’eau liquide sont ras-

semblées dans le tableau 20. L est une longueur de référence qui devient plus complexe à
Lorsque l’évacuation de la chaleur au générateur de vapeur est caractériser dans les grands volumes et varie selon les configura-
insuffisante, le fluide peut devenir diphasique Les relations sont tions étudiées. Cela peut être la hauteur sur laquelle le fluide est
plus complexes à établir mais reposent sur les mêmes hypo- chauffé par exemple. Les relations entre Nusselt, Grashof et
thèses. Prandtl présentées sont de fait liées à la définition de cette lon-
Cette circulation de l’eau permet l’évacuation de la puissance gueur de référence.

résiduelle après l’arrêt automatique du réacteur, les pompes pri- Le nombre de Nusselt est lié au nombre de Grashof suivant :
maires étant stoppées, tant que le GV peut refroidir le fluide pri-
maire (s’il n’est donc pas asséché) et tant que l’inventaire en eau (128)
du primaire n’est pas trop dégradé, le niveau d’eau devant être
suffisant pour permettre la circulation du fluide. Le produit Gr · Pr forme le nombre de Rayleigh :
4.5.7.2 Corrélations d’échange thermique en convection

naturelle
Les valeurs de la constante C et de l’exposant m dépendent du
De manière générale, les corrélations présentées au para- type d’écoulement et de la configuration [40] (tableau 21). Les
graphe 4.5.1 peuvent être appliquées. La corrélation de Dittus- propriétés du fluide sont estimées à la température de film
Boelter (116) a un grand domaine de validité et peut être utilisée (moyenne de la température de paroi et de la température « à
mais requiert l’évaluation de la vitesse moyenne. Dans le cas de l’infini »). Ces relations sont valables pour des conditions aux
très faibles vitesses, si l’écoulement est laminaire, les équations limites imposées sur les températures. Si le flux est imposé sur la
peuvent être résolues analytiquement et dans les cas d’un écoule- plaque, alors on prendra ces relations en considérant la tempéra-
ment en conduite circulaire, les relations exactes proposées en ture de film (pour évaluer les propriétés du fluide) au milieu de la
4.5.1 sont applicables. plaque. Des itérations sont nécessaires pour converger vers la
Dans des grands volumes, lorsque l’eau est chauffée en partie valeur de la température de film.
basse et refroidie en partie haute, des mouvements de circulation En utilisant la température de film et avec la hauteur de la
naturelle s’établissent sous forme de macrostructures en rotation. plaque comme longueur caractéristique, on peut utiliser la rela-
tion suivante plus précise :
Tableau 20 – Coefficients de dilatation pour

quelques valeurs relatives à l’eau liquide
(129)
P [bar] 1,0 P [bar] 150,0
T [°C] β [1/K] T [°C] β [1/K]
5,0 0,000016 25,0 0,000274

Tableau 21 – Convection naturelle : valeurs
10,0 0,000088 75,0 0,000599 des constantes dans l’équation (128)
15,0 0,000151 125,0 0,000849 Configuration C m Condition
20,0 0,000207 175,0 0,001119 Plaque verticale 0,59 1/4 104 < Gr · Pr < 109
25,0 0,000257 225,0 0,001496 0,10 1/3 109 < Gr · Pr < 1013
50,0 0,000457 275,0 0,002183 Plaque 0,54 1/4 104 < Gr · Pr < 107
horizontale
75,0 0,000613 300,0 0,002865 0,15 1/3 107 < Gr · Pr < 1011

BN 3 050v2 – 32
Pour un cylindre horizontal, en utilisant la température de film On souhaite maintenant estimer la température sur la gaine des
et avec le diamètre comme longueur caractéristique, on peut utili- crayons de l’assemblage, ou plus précisément estimer si la tempéra-
ser la relation suivante : ture de vaporisation peut être atteinte. Si l’on suppose qu’en sortie
d’assemblage la température moyenne de l’eau est de 60 °C, on
impose alors une différence de température maximum de 40 °C entre
la température de paroi et la température de l’eau.
D Avec cet écart maximum de température (paroi-liquide), on utilise
(130) les mêmes valeurs des propriétés physiques de l’eau que ci-dessus
(50 °C) en première approximation :
Exemple 11
Un assemblage de cœur de réacteur à eau sous pression est placé Pr = 3,57 et Nu = 7046 – relation (129).
en piscine de désactivation, un mois après l’arrêt du réacteur.
La puissance d’un crayon au régime nominal est de 78,5 kW
(exemple 2). La longueur du crayon est L = 4,27 m. Le diamètre
hydraulique considéré est Dh = 10,86 mm (exemple 1). La section de
passage de l’assemblage (exemple 1) est A = 0,002404 m². Le flux qui peut être évacué par convection naturelle avec un ΔT de
40 °C vaut :
La puissance du crayon un mois après l’arrêt est voisine de 2 % de
la puissance nominale soit :
La surface totale d’un crayon combustible valant : 4,27 × 9,5 × 4,27

= 0,127 m2, la puissance évacuée vaut 42 279 × 0,127 = 5 388 W,
ce qui donne pour les 264 crayons d’un assemblage : soit 5,4 kW.
Cette puissance est nettement supérieure à la puissance réelle du
Le taux d’expansion β de l’eau à la pression atmosphérique à 50 °C barreau 1,57 kW. En conséquence, on peut dire que le combustible
est 0,000457 K-1 (relation (126)) : est refroidi par convection naturelle à l’intérieur de la piscine de stoc-
kage en fonctionnement normal.
De façon identique à la démarche adoptée pour le thermosiphon,
nous considérons ici l’écoulement naturel du fluide dans l’assem-
blage. Dans la piscine, au-dessus des racks d’assemblage, des mou-
vements libres vont se mettre en place dans le large volume.
Dans l’assemblage, la perte de charge est égale à la force
5. Sûreté des REP
d’Archimède :
L’objectif fondamental de sûreté est de protéger les personnes
et l’environnement contre les effets nocifs des rayonnements ioni-
sants [34]. Il est donc essentiel de confiner les matières radioac-
tives dans l’installation. Ceci est réalisé au moyen de trois
barrières indépendantes et successives ; la gaine contenant le
et donc
combustible nucléaire, le circuit primaire et l’enceinte du bâtiment
réacteur. La démonstration de la sûreté du réacteur repose notam-
ment sur l’évaluation de la tenue de ces trois barrières. Cette
démonstration doit prendre en compte non seulement tous les
états du réacteur, du fonctionnement à pleine puissance à l’arrêt à
froid pour rechargement mais aussi les conditions qui peuvent
résulter de tout dysfonctionnement ou événement, du possible au
très improbable.
ΔT étant la différence de température du fluide entre l’entrée et la Les critères de sûreté s’appuient sur la nécessité de préserver
sortie de l’assemblage. en fonctionnement normal, mais aussi dans les cas d’incidents de
Le coefficient de perte de charge par frottement f est estimé, après fonctionnement les plus fréquents (au sens de la probabilité
itérations, à 0,039 d’après la formule (68) : d’occurrence), l’intégrité de la première barrière de confinement
(la gaine du combustible) ; du point de vue thermohydraulique,
c’est le risque de crise d’ébullition qu’il faut évaluer (§ 5.1). Pour
ce qui concerne les accidents (statistiquement très improbables),
ce sont les accidents de perte de réfrigérant primaire (APRP) suite
à une brèche dans le circuit primaire, qui font appel à un grand
nombre de phénomènes thermohydrauliques détaillées dans la
suite (§ 5.2). On trouvera dans [35] une description détaillée des
analyses permettant d’évaluer la maitrise des situations acciden-
Le débit massique vaut 5,19 kg/s, la vitesse massique 216 kg/m².s et la telles du système réacteur.
vitesse d’écoulement de l’eau dans l’assemblage vaut 0,219 m/s.
On vérifie a posteriori le calcul de f avec Re = 4292 et f = 0,039.
5.1 Crise d’ébullition
L’échauffement de l’eau entre l’entrée et la sortie de l’assemblage
est : La crise d’ébullition caractérise une soudaine diminution du
transfert de chaleur entre la gaine du combustible et le fluide.
Π Cette réduction se produit lorsque le régime d’écoulement est
modifié. La gaine se recouvre d’une poche de vapeur qui l’isole du
fluide. Il s’ensuit une augmentation brutale de la température de

BN 3 050v2 – 33
paroi. L’ébullition nucléée qui jusqu’à cette limite améliorait le

transfert de chaleur en convection forcée devient alors ‹ critique ›.
Le déclenchement du phénomène est appelé crise d’ébullition ou Assèchement Caléfaction
Departure from Nucleate Boiling (DNB). Pour une valeur trop éle-
vée du flux de chaleur surfacique, le taux de production de vapeur
est trop important pour être évacué par l’écoulement depuis la
paroi chauffante vers le cœur de l’écoulement. La crise d’ébulli-
tion (ou DNB) apparaît donc pour une valeur de flux maximale
autorisée : le flux critique. On préserve alors l’intégrité de la gaine
de combustible en évitant d’atteindre cette valeur critique.
La valeur de flux critique est, par essence même du phénomène
physique, liée aux conditions de l’écoulement. Le flux critique n’a
pas la même valeur en écoulement fortement sous-saturé ou en
régime saturé, et dépend aussi, entre autres, de la vitesse de
l’écoulement. En l’absence de mouvement du fluide dans des
grands volumes, on parle de crise d’ébullition en vase ou en bain
(cf. [BE8235]) et les phénomènes les plus influents sont la créa-
tion, l’expansion et les tailles de bulles ainsi que les mouvements Figure 15 – Atteinte de la crise d’ébullition par assèchement ou par
caléfaction
de convection naturelle. Dans les réacteurs à eau pressurisée la
crise d’ébullition en vase peut se rencontrer en cas d’accident
grave pour les réacteurs de conception récente (AP1000) qui privi- À géométries identiques, et pression et vitesse massique
légient la rétention en cuve du corium en cas de fonte significative fixées, de manière générale le flux critique est une fonction
du cœur. La crise d’ébullition peut apparaître dans l’eau lors du décroissante, voire quasi-linéaire, de la sous-saturation. De ce
noyage du puits de cuve. Pour ce qui concerne la crise d’ébullition fait, le flux critique est une fonction décroissante du titre.
en convection, historiquement les études ont porté sur la crise Cependant, on observe deux comportements assez différents
d’ébullition dans la géométrie la plus simple, le tube de section lors de la réalisation des essais. Considérons un protocole
circulaire. Ultérieurement, des essais ont été réalisés dans des d’essai qui consiste à augmenter progressivement et lentement
géométries plus proches de la configuration des cœurs (figure 1) la puissance aux bornes de la section d’essais, tout en mainte-
avec des simulateurs électriques des barreaux combustible. nant constantes la vitesse massique, la pression et la tempéra-
ture en entrée. La crise d’ébullition est détectée lorsque l’on
Différents auteurs ont abordé le problème de l’évaluation du observe une augmentation rapide de la température sur la paroi
flux critique par des approches analytiques avec comme hypo- chauffante. Dans les conditions fortement sous-saturées (donc
thèse soit la présence d’une couche de vapeur en paroi, soit l’éva- forte puissance et faible titre), l’augmentation de la température
poration de la sous-couche liquide en contact avec la paroi. Ces est brutale, rapide et importante (plusieurs dizaines de degrés).
modèles ne sont pas encore assez développés et précis pour être À mesure qu’en entrée de la section d’essais, on approche de la
appliqués dans l’industrie nucléaire. saturation, le flux critique est plus faible. Cela correspond aussi
à un titre plus important. Dans ce cas, la crise d’ébullition est
5.1.1 Évaluation du flux critique plus difficile à identifier. La température en paroi augmente len-
tement. Dans le premier cas, la crise d’ébullition est obtenue par
On présentera des approches expérimentales relativement caléfaction, création locale – limitée à l’environnement proche
faciles à mettre en œuvre. Ces relations, ou corrélations de flux de la paroi chauffante – de vapeur, le reste de l’écoulement est
critique, relient une estimation du flux critique mesuré lorsque la liquide. Dans le second cas, le taux de vapeur dans l’écoulement
crise d’ébullition est détectée, à divers paramètres thermohydrau- est important et le liquide à la paroi ‹ s’assèche › progressive-
liques ou géométriques. Le flux critique peut s’exprimer en fonc- ment jusqu’à la crise d’ébullition (cf. figure 15). En consé-
tion des conditions globales mesurées lors des essais, vitesse quence, il sera difficile d’obtenir une régression qui puisse
massique ou débit, température ou enthalpie d’entrée, pression couvrir avec une précision identique ces deux types de crise
en entrée ou en sortie de la section d’essais. Dans ce cas, des d’ébullition. Par ailleurs, les conditions thermohydrauliques et
paramètres géométriques, et notamment la hauteur chauffante, les valeurs locales de flux correspondant aux plages de fonc-
sont indispensables pour évaluer le flux critique. L’autre approche tionnement des réacteurs à eau pressurisée amènent majoritai-
consiste à relier le flux critique aux grandeurs thermohydrauliques rement au risque de crise d’ébullition par caléfaction.
locales, principalement la vitesse massique, la pression et le titre L’assèchement relève plus du domaine de fonctionnement des
thermodynamique à l’équilibre (relation (6)). De fait les grandeurs réacteurs à eau bouillante. Dans certaines situations, la crise
locales ne sont pas mesurées mais reconstituées à partir des d’ébullition est obtenue par atteinte d’une valeur critique de la
conditions en entrée ou sortie de la section d’essais et de la géo- qualité. On s’intéresse dans cet article aux relations évaluant le
métrie des essais. En tube, cette reconstitution est immédiate. La flux critique par atteinte de la crise d’ébullition en caléfaction.
vitesse massique est constante, la pression estimée à l’aide des
pertes de charge (cf. § 4.3) et le titre thermodynamique est calculé 5.1.1.1 Corrélations en tube
à l’aide d’un bilan enthalpique. Les incertitudes sur l’enthalpie
En tube, le titre thermodynamique χ, en sortie d’un tube de lon-
d’entrée et sur la puissance se retrouvent donc combinées dans
gueur L et diamètre D, se calcule très simplement à partir des
l’évaluation du titre. Pour cette raison, les relations exprimées en
conditions d’essai pression P, de la vitesse massique G, de la tem-
fonction des paramètres locaux ont une dispersion plus impor-
pérature d’entrée Te et du flux critique φc selon :
tante que les relations à partir des conditions d’entrée. À l’aide
de régressions mathématiques, le flux critique est exprimé à
partir des conditions d’entrée ou locales. Des tables ont aussi
été proposées, ce qui supprime l’étape de la régression et
l’erreur mathématique entre la valeur issue de la régression et
la valeur mesurée. Cependant cette approche impose d’intro-
duire l’ensemble des valeurs tabulées et d’interpoler entre les Suivant l’approche choisie, les relations suivantes sont éta-
valeurs de la table. blies. En condition d’entrée ou en condi-

BN 3 050v2 – 34
tions locales . L’introduction de la longueur et du – pour Pr > 1 (MPa)

diamètre du tube est discutable. Pour des essais relatifs à une
longueur et un diamètre, alors une corrélation en condition
d’entrée peut à l’évidence ne pas introduire ces paramètres
dans son expression. Son application sera restreinte à cette
unique géométrie. A contrario, en conditions locales, certains
auteurs ont observé à partir des essais des effets de la longueur
chauffante et/ou du diamètre qu’ils ont introduits en consé-
quence dans leurs corrélations.
■ Corrélation de Biasi
Le titre étant calculé à l’aide du flux critique, une régression en
conditions locales peut sembler mathématique- La corrélation de Biasi [3] est valable pour des géométries en tube :
– de diamètre D compris entre 0,003 et 0,037 m ;
ment mal posée. La fonctionnelle va lier la variable ‹ dépendante ›
– de longueur L comprise entre 0,2 et 6 m.
φc aux variables ‹ indépendantes › P, G et χ. Il faut dans ce cas
considérer le titre non pas comme le résultat d’un calcul mais Les paramètres thermohydrauliques suivants sont couverts par
comme une grandeur mesurée, indépendamment de la pression, cette corrélation :
de la vitesse massique et du flux critique. Malgré cette lacune, – pression P entre 0,27 et 14 MPa ;
nous présentons ici des relations en conditions locales car plus – vitesse massique G entre 100 et 6 000 kg/m².s ;
générales que les corrélations en conditions d’entrée. Signalons 1
– titre entre et 1.
cependant le cas particulier de la corrélation W3 qui utilise à la ρ
fois l’enthalpie d’entrée et le titre local. Dans ce cas, l’enthalpie et 1 + ρl
v
le titre d’entrée sont forcément liés et ne peuvent être indépen-
Le flux critique est le maximum des deux fonctionnelles suivantes :

dants. La fonctionnelle est alors mal posée.
Pour des écoulements en tube (donc assez éloignés de l’écoule-
ment dans le cœur des REP), on peut citer les corrélations sui-
vantes.
■ Corrélation de Bowring
La corrélation de Bowring [3] est valable pour des géométries
en tube :
– de diamètre D compris entre 0,002 et 0,045 m ; La constante n en puissance du diamètre du tube est la suivante :
– de longueur L comprise entre 0,15 et 3,7 m. • n = 0,6 pour D < 0,01 m ;
Les paramètres thermohydrauliques suivants sont couverts par • n = 0,4 pour D ≥ 0,01 m.
cette corrélation : Les fonctions f(P) et h(P) de la pression (P en MPa) s’expriment
– pression P entre 0,2 et 19 MPa ; suivant :
– vitesse massique G entre 136 et 18 600 kg/m².s ;
– pas de restriction en titre.
Comme le facteur B est directement proportionnel au diamètre,
au-delà d’une certaine valeur du diamètre, le flux critique devient
négatif aux fortes vitesses massiques et à forte qualité. Il est donc
■ Corrélation W3
recommandé de valider les prédictions obtenues avec la corréla-
tion de Bowring avec d’autres corrélations : La corrélation de W3 [3] a été proposée par Westinghouse et
établie par L. S. Tong. Cette corrélation a été développée à partir
d’essais en tube :
3 – de diamètre D compris entre 5 et 17,8 mm ;
– de longueur L comprise entre 0,254 et 3,65 m.
avec Les paramètres thermohydrauliques suivants sont couverts par
cette corrélation :
2,317
– pression P entre 6,9 et 15,9 MPa ;
– vitesse massique G entre 1356 et 6780 kg/m².s ;
– titre thermodynamique à l’équilibre χ entre –0,15 et +0,15.
Le flux critique est le produit des fonctions :
Les fonctions F1 à F4 sont différentes suivant la valeur de Pr :

– pour Pr < 1 (MPa)
avec

BN 3 050v2 – 35
On remarque que cette corrélation est fonction à la fois de dans la corrélation qui est déterminée de façon empirique. Les
conditions locales (titre thermodynamique) et en entrée (enthal- approches plus analytiques et à des échelles plus fines que le
pie). Le titre est lui-même fonction de l’enthalpie d’entrée et du sous-canal ne sont pas encore suffisamment abouties pour être
flux (intégré entre l’entrée et la hauteur considérée). Cette non- intégrées dans les analyses de sûreté.
indépendance des paramètres explicatifs (P, G, χ, he) utilisés dans Les caractéristiques de conception du combustible nucléaire, en
la régression conduit à une réduction artificielle de la dispersion particulier les grilles, la structure (géométrie) de l’assemblage, ont
de la valeur expliquée φ [41]. De plus, il n’est pas possible d’éva- une influence importante sur le flux critique. En favorisant la créa-
luer avec confiance les situations où la combinaison du flux cal- tion de vitesses transverses ou l’intensité de la turbulence ou
culé par la corrélation et de l’enthalpie d’entrée considérée ne encore l’éclatement des poches de vapeur, on repousse la valeur
correspond pas au titre introduit dans la corrélation. du flux critique. Les corrélations de flux critique sont dépendantes
■ Tables de Groeneveld du combustible et sont confidentielles, chaque fabricant établit sa
corrélation de flux critique applicable sur ses produits. Les autori-
À partir de données en tubes, initialement compilées par V.E. tés de sûreté vérifient le bien-fondé de cette corrélation. Pour
Doroshchuk en 1975, des valeurs du flux critique ont été propo- obtenir un ordre de grandeur du flux critique en grappe, des
sées par D.C. Groeneveld [52] [53] [54], pour des valeurs fixées de auteurs proposent des facteurs correctifs à introduire dans des
pression, vitesse massique et titre et uniformisées à un diamètre corrélations établies en tube [27]. Cependant, un grand nombre de
interne de tube de 8 mm. campagnes d’essais réalisées dans les années 1980 à l’université
L’utilisation de ces tables [54] est valable pour des géométries Columbia de New-York ont été mises dans le domaine public [44]
en tube : et peuvent être utilisées librement. Une de ces campagnes a été
– de diamètre D compris entre 0,003 et 0,025 m avec un facteur exploitée et une corrélation simple a été proposée [45]. Les para-
correctif multiplicatif à appliquer pour un diamètre différent de mètres thermohydrauliques suivants sont couverts par cette cor-
8 mm ; rélation :
– de rapport longueur sur diamètre supérieur à 100. – pression P entre 10 et 17 MPa ;

Les paramètres thermohydrauliques suivants sont couverts par – vitesse massique G entre 1 300 et 4 700 kg/m².s ;
ces tables [54] : – titre thermodynamique à l’équilibre χ entre 0,05 et 0,30.
– pression P entre 0,1 et 21 MPa ;
– vitesse massique G entre 0 et 8 000 kg/m².s ;
– titre inférieur à 1.
Divers facteurs correctifs multiplicatifs sont proposés [52] et [3]
pour prendre en compte les effets de géométries différentes du Cette corrélation contient exclusivement les paramètres locaux
tube ou l’effet de grilles. et en particulier n’introduit pas de paramètres géométriques rela-
Ces tables permettent d’avoir un ordre de grandeur satisfaisant tifs à l’assemblage. Cependant les grilles de mélange ont une
pour des évaluations préliminaires. L’introduction de toutes les influence significative sur les performances en flux critique et une
valeurs tabulées et la mise en place de moyens d’interpolation grille optimale va permettre d’atteindre des valeurs de flux cri-
pour des conditions en dehors de valeurs tabulées rendent cette tique élevées, à conditions locales identiques. La distance entre le
méthode peu pratique. La correction de diamètre hydraulique est point d’évaluation et la position de la grille en amont est donc un
aussi limitée car uniquement liée au rapport des diamètres sans paramètre premier dans les corrélations de flux critique en
prendre en compte les conditions thermohydrauliques locales. grappe. La corrélation présentée calcule le flux critique en consi-
dérant une grille en amont suffisamment loin pour que l’effet de
5.1.1.2 Corrélations en grappe cette distance ne soit plus significatif.
Le terme générique de corrélation en grappe est utilisé pour Pour des écoulement dans des grappes de crayons disposés en
décrire les moyens de calculs applicables aux géométries des pas triangulaire, caractéristique des réacteurs à eau pressurisée
cœurs de réacteur. Étant données les géométries très particulières de conception russe (voir [BN3101]), la corrélation de Bezrukov ou
rencontrées, une approche expérimentale est indispensable pour OKB Gidopress de 1976 est assez simple à mettre en œuvre et
évaluer les valeurs du flux critique. À partir des données expéri- donne des résultats satisfaisants y compris pour des évaluations
mentales, de la même façon que pour les géométries en tube (cf. en grappe à pas carré pour des grilles peu mélangeuses.
§ 5.1.1.1), les régressions sont établies à partir de conditions glo- Les paramètres thermohydrauliques suivants sont couverts par
bales telles que mesurées dans les essais ou de conditions locales cette corrélation :
calculées. Contrairement à une géométrie en tube, ces données
– pression P entre 7,45 et 16,7 MPa ;
locales, en particulier vitesse massique et titre thermodynamique,
ne peuvent se calculer simplement. La vitesse massique peut – vitesse massique G entre 700 et 3800 kg/m².s ;
varier selon sa position dans la section d’essai compte tenu des – titre thermodynamique à l’équilibre χ entre –0,07 et 0,40.
redistributions de débits. Le titre thermodynamique varie signifi-
cativement selon les sous-canaux (figure 1). Il faut alors recourir à
des logiciels de calcul adaptés, les codes d’analyse par sous-
canaux, tels que FLICA [42] ou COBRA [43]. Les cœurs des réac-
teurs à eau pressurisée sont ouverts, à la différence des cœurs
des réacteurs à eau bouillante où chaque assemblage de combus- Exemple 12
tible est isolé de ses voisins, en outre un assemblage de REP com- Les valeurs de flux critique pour différentes corrélations et condi-
porte plusieurs centaines de crayons. Des essais en vraie échelle tions locales sont données dans le tableau 22. Pour les corrélations
sont alors irréalisables. L’approche en conditions locales est alors en tube, un diamètre de 8 mm est utilisé. Les corrélations en grappe
indispensable pour des applications de la corrélation de flux cri- n’ont pas de terme de diamètre mais le diamètre hydraulique corres-
tique à des calculs de cœur de REP. pondant est d’environ 11 mm.
Par approche en conditions locales on entend ici des conditions Les valeurs sur fond grisé (tableau 22) sont hors du domaine sur
calculées tout au long de l’écoulement mais à l’échelle du sous- lequel les modèles correspondants ont été établis. Pour la corrélation
canal (figure 1). Les phénomènes existants à une échelle plus fine W3, une enthalpie en entrée identique a été introduite pour chaque
(par exemple dans la couche limite thermique entre la paroi et le valeur de la pression à savoir à 15 MPa 1284 kJ/kg (soit 290 °C) et à
fluide) ne sont pas modélisés et sont pris en compte globalement 7 MPa 991 kJ/kg (230 °C).

BN 3 050v2 – 36
On remarque une assez grande dispersion des résultats et des (ζ = 0). Le facteur K est établi à partir de données expérimentales
variations avec la vitesse massique et la pression assez différentes suivant le modèle suivant :
selon les corrélations. Les extrapolations sont acceptables (sur les
domaines présentés). Une grande prudence est recommandée dans
les estimations des valeurs de flux critique qui peuvent varier selon
les sources de +/–25 % et pour lesquelles l’incertitude à associer
peut atteindre +/–50 %. Rosal [3] propose les valeurs suivantes pour a, b et c, pour une
intégration à partir de l’ébullition locale z0 : a = 185,6 ; b = 4,31 ;
5.1.1.3 Variation axiale du flux de chauffage c = 0,478. K est en m–1 si G est en kg/(m2.s).
Il est observé (dans des géométries en tube ou en grappe) que le Établies à l’origine pour des unités anglo-saxonnes (in–1 et
flux critique atteint pour des essais avec un chauffage uniforme est Mlb.h ft² ; a = 0,15 et c = 0,478), le coefficient a est seul modifié et
significativement différent (de quelques pourcents à plusieurs se calcule à l’aide des facteurs de conversions (cf. § 8) :
dizaines de pourcents) de celui mesuré si le chauffage est variable
axialement. Notamment, le flux critique en chauffage uniforme est
plus élevé qu’en flux variable dans la zone de flux décroissant. Étant
donné que dans le réacteur le flux de conception considéré est sou-
vent un flux en forme de cosinus et que pour de nombreuses appli-
cations, le flux est variable, cet aspect doit être pris en compte. 5.1.2 Évaluation du risque de crise d’ébullition
L’effet « d’histoire » est introduit par plusieurs auteurs en support au
modèle qu’ils proposent. Plus directement, on peut relier cet effet L’évaluation du risque de crise d’ébullition est indispensable à
d’histoire au constat que même si les paramètres thermohydrau- l’examen de l’intégrité de la première barrière de confinement (cf.
liques sont identiques pour un chauffage uniforme ou non en § 5.1). Des calculs sont réalisés dans tout le cœur à l’aide des
moyenne à l’échelle « macroscopique », la répartition locale du taux codes d’analyse par sous-canaux et ceci pour les conditions nor-
de vide et du titre, entre proche paroi et centre du sous-canal, peut males et accidentelles couvrant l’ensemble des situations et
être sensiblement différente selon le mode de chauffage axial. modes de fonctionnement du réacteur. À partir de ces conditions,
les valeurs locales des paramètres thermohydrauliques ainsi cal-
Le modèle le plus courant consiste à appliquer un facteur cor- culées permettent de déterminer un flux critique qui est comparé
rectif multiplicatif pour tenir compte d’une variation axiale de la au flux local considéré. On forme alors le rapport de flux ther-
forme du flux. Il faut diviser la valeur du flux critique par la valeur mique critique RFTC :
du facteur de flux non uniforme FNU qui a été modélisé suivant la
relation :
Ces informations sont synthétisées dans la figure 16. Les para-

mètres importants sont tracés en fonction de la hauteur dans le
La distance d’intégration zc est comptée, selon les auteurs, à partir cœur. Le flux local considéré est issu d’une forme de conception
du début de l’ébullition z0 ou du début de la longueur chauffante (en cosinus) et relatif à la puissance totale du réacteur et à la dis-
Tableau 22 – Valeurs de flux critiques pour différentes corrélations et conditions locales
PGX Groeneveld
Bowring Biasi W3 RR1 BZK
MPa kg/m² 1986 1995 2006
15 2000 –0,10 3,425 4,009 3,372 3,244 3,351 3,350 2,233 3,266
15 2000 0,00 2,572 3,118 2,584 2,391 2,464 2,468 1,824 2,327
15 2000 0,10 1,719 2,228 1,829 1,673 1,843 1,687 1,415 1,640
15 4000 –0,10 3,888 3,268 4,395 3,906 4,502 4,317 2,842 3,791
15 4000 0,00 2,525 2,475 3,034 2,953 2,995 2,874 2,230 2,644
15 4000 0,10 1,161 1,682 1,804 1,906 2,270 2,093 1,617 1,824
7 2000 –0,10 5,930 6,275 4,965 6,044 6,142 6,769 2,900 3,632
7 2000 0,00 5,132 5,384 4,028 4,728 4,952 5,327 2,600 2,804
7 2000 0,10 4,335 4,494 3,438 3,685 3,785 4,020 2,301 2,159
7 4000 –0,10 6,456 5,067 6,471 6,334 7,398 7,118 3,508 4,216
7 4000 0,00 5,096 4,274 4,730 4,638 4,372 4,226 3,006 3,185
7 4000 0,10 3,736 3,480 3,391 2,929 2,928 2,855 2,504 2,401

BN 3 050v2 – 37
(gaussienne) de ces écarts et l’on utilise les coefficients d’Owen

[46]. Wilks [47] propose aussi l’établissement d’un critère à partir
RFTC RFTC χ d’une distribution non paramétrique des écarts.
4 0,3
χ (–)
Flux critique φcr
5.2 Accident de perte de réfrigérant
0,2
du circuit primaire (APRP)
Dans ces accidents, on postule une brèche du circuit primaire
3 (seconde barrière) et une dégradation limitée de la première bar-
0,1 rière (gaine du combustible) peut éventuellement être acceptée. Il
est indispensable non seulement de conserver l’intégrité de
l’enceinte de confinement mais aussi de pouvoir refroidir le cœur
Flux local du réacteur, d’éviter la fusion du combustible et les réactions
0,0 entre l’eau et le Zircaloy. Des systèmes de secours (injections de
2 sécurité moyenne ou haute pression, accumulateurs passifs et
injections de sécurité basse pression, aspersion enceinte) sont
conçus pour limiter les conséquences de tels accidents. L’analyse
–0,1
de ces accidents a pour objectif de vérifier le bon dimensionne-
Valeurs numériques ment de ces systèmes, avec leurs réserves d’eau associées, et la
indicatives… maîtrise des conséquences de tels accidents.
1
0 1,00 2,00 3,00 4,00 5.2.1 Méthode d’étude

z : hauteur dans le cœur [m]
La complexité des phénomènes thermohydrauliques consécu-
tifs à une brèche du circuit primaire requiert le recours à des logi-
Figure 16 – Évolution du RFTC et des grandeurs thermohydrauliques ciels de calcul [48] pour évaluer les valeurs des paramètres cibles,
dans un sous-canal en particulier la température de la gaine du combustible et la
pression dans l’enceinte. Les valeurs calculées sont comparées
aux valeurs limites règlementaires et valident (ou non) le dimen-
tribution radiale de cette puissance au sein des assemblages et sionnement des systèmes de sauvegarde et les limites de fonc-
des crayons des assemblages les plus chauds. À partir de la tem- tionnement du réacteur. On rappelle que ces analyses sont
pérature d’entrée et de la vitesse locale (non représentée), le titre menées non seulement pour les APRP mais aussi pour l’ensemble
thermodynamique est calculé. Il est donc croissant tout au long de des incidents et accidents résultant de tous les évènements initia-
la hauteur du cœur. Le flux critique est calculé sur toute la hauteur teurs imaginables [BN3820].
et décroît (vitesse massique et pression sont relativement
Lors d’un APRP, l’onde de décompression associée à l’ouverture
constants et le titre augmente). Les petites ruptures de pente
de la brèche se propage dans le circuit primaire et soumet les struc-
indiquent que la corrélation utilisée comporte un terme de dis-
tures internes de la cuve, les assemblages combustibles et les com-
tance à la grille. Lorsque cette distance est réinitialisée (à chaque
posants primaires à des efforts hydrauliques importants. Les
hauteur comportant une grille), le flux critique croît légèrement.
contraintes mécaniques imposées aux internes de cuve et aux
Le rapport du flux critique au flux local (RFTC) amplifie ces rup-
assemblages combustibles doivent rester limitées pour que la géo-
tures de pente au passage d’une grille. L’évolution du RFTC n’est
métrie du cœur et des assemblages de combustible permette le
pas monotone. Le RFTC décroît au départ compte tenu de la
refroidissement de ces derniers par les moyens de sauvegarde.
décroissance du flux critique puis augmente lorsque le flux local
Cependant, on s’intéresse ici uniquement à l’aspect thermohydrau-
est faible.
lique des APRP, les conséquences mécaniques, la démarche de la
Le RFTC présente donc un minimum. À la hauteur correspon- démonstration de sûreté, la détermination et le respect des critères,
dante, le risque de crise d’ébullition est le plus élevé. Si la valeur en particulier l’introduction de conservatismes ainsi que la prise en
minimum du RFTC est inférieure à l’unité, le flux local dépasse la compte des incertitudes, dépassant le cadre de cet article.
valeur du flux critique et la crise d’ébullition ne peut être évitée.
En France, l’étude thermohydraulique précise de cette catégorie
Cependant, l’évaluation du flux critique repose sur des essais d’accident a nécessité le développement de l’outil de calcul scien-
(entachés d’incertitudes expérimentales) et sur une fonctionnelle tifique CATHARE (Code Avancé de THermohydraulique pour les
mathématique (la corrélation) qui n’ajuste pas exactement les Accidents des Réacteurs à Eau) [29], réalisé par le CEA, Frama-
valeurs mesurées. De plus, les essais ne constituent qu’un échan- tome, EDF et l’IRSN à partir de 1979. CATHARE simule tout le
tillon de la « population » du flux critique. Il faut définir un inter- spectre des grandes aux petites brèches soit du circuit primaire,
valle de tolérance statistique permettant de prendre en compte soit du circuit secondaire, soit des circuits auxiliaires (RRA), ainsi
ces différents facteurs : que tout autre transitoire accidentel (autre que les brèches) ou des
– les incertitudes expérimentales et l’erreur de la fonctionnelle transitoires d’exploitation. La description de ces transitoires
mathématique sont transcrites au travers de la moyenne et de dépasse le cadre de cet article. Les transitoires d’APRP mettent en
l’écart type associé à la corrélation et aux essais ; lumière de nombreux phénomènes thermohydrauliques et sont
– la taille de l’échantillon, le taux de confiance et le pourcentage décrits dans la suite. CATHARE constitue un outil d’analyse de
de la population déterminent un coefficient statistique qui permet sûreté et de mise au point des procédures de conduite en situa-
de calculer un seuil unilatéral critique de tolérance appelé critère tion accidentelles et post-accidentelles. Les capacités de calcul
de la corrélation. actuelles permettent son intégration dans les simulateurs de cen-
L’absence de risque de crise d’ébullition est couvert avec une trale pour calculer en temps réel la thermohydraulique des circuits
confiance de 95 % et pour une proportion de 95 % si le RFTC est primaire, secondaire et certains circuits auxiliaires. C’est par
supérieur au critère « 95/95 » de la corrélation. Ce critère peut être exemple le cas du simulateur SOFIA (Simulateur d’Observation du
établi à partir de la distribution des écarts entre les mesures expé- Fonctionnement Incidentel et Accidentel) [49].
rimentales et les calculs établis avec la corrélation. Le critère peut CATHARE est un code à deux fluides et six équations (cf. § 3.4).
être établi à l’aide de l’hypothèse d’une distribution normale Les équations de bilan de masse, d’énergie et de quantité de mou-

BN 3 050v2 – 38
vement pour les phases vapeur et liquide du fluide sont résolues diées séparément du fait des scénarios très différents qui les
dans un circuit modélisé par des éléments à une dimension, caractérisent.
comme les conduites ou des éléments 0D (volume). Un module
3D est également disponible pour simuler les écoulements possi-
blement tridimensionnels dans le cœur et la cuve du réacteur. Les 5.2.2 Accident de perte de réfrigérant primaire :
dernières versions du code offrent également la possibilité de grosse brèche
modéliser le comportement de l’enceinte de confinement en situa- Un transitoire APRP grosse brèche résulte d’une brèche d’un
tion accidentelle (thermohydraulique avec forte présence d’un diamètre équivalent supérieur à 25-35 cm, le cas extrême étant la
incondensable) ou de simuler des réacteurs de génération IV rupture franche (rupture guillotine doublement débattue) de la
(réacteur rapide au sodium, réacteur à haute température…). La tuyauterie, le fluide s’écoulant alors par les deux tronçons
version 2.5 du code CATHARE est actuellement l’outil de référence ouverts. Lors de la conception d’un nouveau réacteur, ce cas
des analyses de sûreté en France. Son successeur CATHARE 3 limite de l’APRP dimensionne les accumulateurs et l’injection de
[31] intègre notamment un modèle à deux fluides et trois champs sécurité basse pression.
afin de pallier le manque de précisions du modèle à six équations
pour certaines configurations d’écoulement. Ce sont les brèches localisées en branche froide (entre une
pompe primaire et la cuve) qui sont, du point de vue thermohy-
Les circuits primaire et secondaire sont découpés en volumes draulique, les plus pénalisantes pour la tenue des crayons com-
élémentaires (mailles) dans chacun desquels les équations du bustible. En effet, cette situation conduit à postuler que l’eau des
modèle à deux fluides sont résolues. La figure 17 en donne un systèmes de sauvegarde (injections de sécurité haute pression et
exemple pour le circuit primaire. une partie des accumulateurs) de la boucle rompue ne sert pas à
À l’international et en particulier aux États-Unis [50], les logi- refroidir le cœur car elle est directement perdue à la brèche.
ciels TRACE [51], RELAP [33] proposent des fonctionnalités simi- L’APRP grosse brèche est un transitoire rapide de l’ordre de
laires à CATHARE. Ces logiciels représentent d’une manière aussi 200 secondes, pendant lesquelles une brusque chute de la pression
réaliste que possible (« Best Estimate ») l’ensemble des phéno- primaire, un dénoyage total du cœur puis un renoyage de ce dernier
mènes physiques intervenant dans la thermohydraulique de l’acci- sont observés. D’un point de vue thermohydraulique le scénario de
dent. Cette approche impose une prise en compte détaillée des cet accident peut alors se décomposer suivant trois grandes phases, à
incertitudes et des lacunes dans les modèles physiques et la dis- savoir la dépressurisation, le remplissage du plenum inférieur de la
crétisation géométrique adoptée. cuve et le renoyage oscillant puis stabilisé du cœur. La figure 18 posi-
Ces logiciels requièrent une validation rigoureuse et doivent tionne ces différentes phases en regard de l’évolution des tempéra-
montrer que l’effet d’échelle entre le support expérimental et les tures de la pastille de combustible et de la gaine.
applications en réacteur est correctement pris en compte. Les L’énergie stockée dans le combustible (initialement et au travers
applications de ces logiciels sont bien entendu limitées au de la loi de décroissance de la puissance résiduelle) et les trans-
domaine couvert par la validation. ferts de chaleur avec le fluide déterminent les valeurs de tempéra-
ture atteintes. De nombreux phénomènes thermohydrauliques
Pour ce qui concerne les calculs APRP, l’emplacement et la taille
complexes, décrits ci-dessous, gouvernent le processus de trans-
de la brèche déterminent le déroulement et la physique de l’acci-
fert de chaleur qui peut rencontrer les divers régimes présentés
dent. Pour les petites brèches d’un diamètre équivalent compris
au paragraphe 4.5.
entre 9,5 et 25 mm, l’étude de sûreté a pour objectif de démontrer
le non-découvrement du cœur. Les brèches intermédiaires (d’un
diamètre compris entre 2,5 et 25 à 35 cm) et les grosses brèches 5.2.2.1 Dépressurisation
(d’un diamètre équivalent supérieur à 25-35 cm) conduisent à un La phase de dépressurisation débute à l’ouverture de la brèche
découvrement plus ou moins prononcé des crayons combustible et s’achève lorsque la pression primaire s’équilibre avec la pres-
et font intervenir des phénomènes physiques complexes. Les
grosses et les petites brèches ou brèches intermédiaires sont étu-
Température de pastille [°C] Température maximale de gaine [°C]
Tpastille = de 430 à 1 300 °C

1 100
1 000 T3
T2
T1
Tgaine = 360-400 °C
10 20 40 80 150 Temps [s]

Phase de décompression
Renoyage
oscillant Renoyage stabilisé
Phase de remplissage
du plénum inférieur
Figure 17 – Exemple de schéma de discrétisation du circuit primaire Figure 18 – Évolution de la température de gaine dans un transitoire
dans un code de calcul type CATHARE grosse brèche

BN 3 050v2 – 39
sion de l’enceinte (après 40 s environ). Dès l’ouverture de la teur annulaire qui doit vaincre les pertes de charge de l’écoule-
brèche, le circuit primaire se vide rapidement, le fluide s’écoule ment de la vapeur et des gouttelettes d’eau vers le haut du cœur,
des deux côtés dans la boucle rompue. Un point de stagnation (où le plénum supérieur, les générateurs de vapeur et les pompes. Au
la vitesse du fluide est nulle) progresse dans le circuit primaire au cours de cette phase, les crayons sont à une température telle que
cours du temps. S’il est situé dans le cœur, les échanges ther- l’eau liquide ne peut venir directement au contact des gaines
miques par convection sont alors très fortement réduits. Lorsque parce qu’elle y est immédiatement vaporisée (point de Leidenfrost
la pression atteint la valeur à saturation correspondant à la tempé- voir [BE8235]). Il y a établissement d’un écoulement de type annu-
rature de l’eau (autour de 330 °C soit 13 MPa dans les zones les laire inverse, la vapeur entoure les gaines et l’eau liquide est reje-
plus chaudes), la vaporisation (« flashing ») très forte entraîne, par tée au centre du sous-canal fluide. La vaporisation en partie basse
effet neutronique, la chute de la réactivité. La pression continue du cœur crée un fort courant ascensionnel de vapeur dont la
de baisser et la vaporisation du fluide augmente, de ce fait la vitesse est de l’ordre de plusieurs dizaines de mètres par seconde.
baisse de pression est ralentie. Tout au long de cette phase, le Les gaines atteignent un nouveau pic de température lorsque le
transfert thermique est limité et les gaines s’assèchent rapide- niveau d’eau atteint la côte du flux maximal de puissance du
ment, leur température s’élève entre 800 et 1 000 C, selon l’éner- cœur, pic T3 (figure 18), dit de contournement (turnaround) puis
gie stockée dans le combustible (et donc très dépendant de la leurs températures baissent sensiblement. Avec la contribution de
conductivité et de la chaleur spécifique de celui-ci). C’est le pic T1 la conduction axiale, il y a progression du front de trempe et
de la figure 18. Puis une circulation temporaire de fluide à l’état retour à saturation. L’accident est terminé pour sa phase la plus
diphasique peut se rétablir dans la cuve, d’où une baisse possible importante, la phase à court terme.
de température après ce pic. La phase de dépressurisation
Les phénomènes thermohydrauliques prépondérants sont liés au
s’achève lorsque l’équilibre des pressions est atteint entre le cir-
contre-courant en sortie du cœur, aux échanges thermiques entre la
cuit primaire (ouvert) et l’enceinte.
gaine et le fluide, à l’entraînement des gouttes d’eau par la vapeur
Les phénomènes thermohydrauliques prépondérants sont tout qui peuvent être piégées par les structures du plénum supérieur (et
d’abord le débit critique à la brèche (§ 4.4), les transferts de cha-
ensuite condensées cf. § 4.5.5) ou vaporisées dans le générateur de

leur (avec notamment la détermination de l’instant de la crise vapeur (steam binding) qui influence la vitesse de montée du front
d’ébullition au-delà duquel les transferts de chaleurs sont forte- de trempe au travers des différentiels de pression.
ment réduits), les pertes de charge dans le circuit (notamment
Dans le GV, le transfert de chaleur est inversé, le secondaire
pour la vidange des autres boucles vers la brèche en particulier
dans la zone haute de l’espace annulaire où l’eau provenant des chauffe le primaire lorsque la pression du secondaire est supé-
boucles froides saines va préférentiellement aller vers la brèche et rieure à celle du primaire.
non vers le cœur) et les phénomènes de séparation des phases
dans les zones à faible vitesse (proche du point de stagnation). 5.2.3 Accident de perte de réfrigérant primaire :
petite brèche et brèche intermédiaire
5.2.2.2 Remplissage
Lors d’un accident de perte de réfrigérant primaire par petite
Cette phase débute avant la fin de la dépressurisation (après 20
brèche ou brèche intermédiaire, on postule la rupture partielle
s environ), dès que la pression dans le circuit primaire est suffi-
d’une grosse tuyauterie primaire ou la rupture d’une petite tuyau-
samment basse pour permettre l’ouverture des clapets des accu-
terie branchée au circuit primaire ou encore une soupape de
mulateurs (4 MPa), ce qui conduit à une arrivée d’eau froide et
sûreté du pressuriseur restée bloquée ouverte après fonctionne-
borée dans les boucles. Une partie en est perdue à la brèche (on
ment (cas de l’accident de Three Mile Island). La rupture d’un tube
postule de plus que celle provenant de l’accumulateur connecté à
de générateur de vapeur (RTGV) est un cas très particulier d’APRP
la branche rompue est perdue en totalité). En effet, la vapeur pro-
du fait de son déroulement totalement différent où les interactions
duite empêche cette eau de descendre dans le collecteur annu-
entre les pressions au primaire et secondaire sont prédominantes.
laire, une autre partie est vaporisée par les structures métalliques
On s’intéresse ici aux aspects purement thermohydrauliques pour
chaudes (branches froides et entrée cuve). L’eau froide peut aussi
lesquels la RTGV ne présente pas des phénomènes spécifiques, ce
condenser la vapeur et participe ainsi à la décompression dans le
type d’APRP ne sera pas exposé ici.
circuit. L’eau parvient finalement à remplir l’espace annulaire et le
plénum inférieur de la cuve. Durant cette phase, le cœur est Le déroulement d’un APRP due à une petite brèche ou une
encore dénoyé et, n’étant pas refroidies, les gaines s’échauffent brèche intermédiaire est plus lent que dans le cas d’une grosse
de façon quasi-adiabatique (car les échanges de chaleur avec la brèche ; l’écoulement est moins homogène (les phases liquides et
vapeur sont négligeables) jusqu’au pic T2 (figure 18). Une fois le vapeur pouvant être « séparées ») et les effets gravitaires plus
plénum inférieur rempli d’eau, l’eau atteint le bas du cœur (début importants déterminent les phases du transitoire ; le circuit pri-
de la phase de renoyage). Cette eau, au contact des crayons com- maire se dépressurise plus lentement, ce qui laisse au système de
bustible très chauds, se vaporise rapidement ce qui provoque sauvegarde (injection de sécurité haute pression notamment) le
l’apparition d’un fort débit vapeur dans le cœur et donc son refroi- temps d’agir avant la vidange du cœur. En outre les actions de
dissement (pic T2). Ce phénomène provoque au début un l’opérateur peuvent modifier le déroulement de l’accident. La taille
renoyage dit oscillant (succession de phases de pénétration d’eau de la brèche, sa position, l’état initial du réacteur sont aussi des
froide puis de vaporisation). facteurs prédominants. Il en résulte une multitude de scenarios
Les phénomènes thermohydrauliques prépondérants sont le possibles qui dépassent le cadre de cet article et dont on trouvera
contre-courant et la condensation de la vapeur au contact de l’eau des descriptions dans [36] et [37]. On peut néanmoins distinguer
des accumulateurs et le destockage d’énergie des masses métal- plusieurs phases parmi lesquelles une phase de dépressurisation,
liques. une autre de baisse du niveau de l’eau dans la cuve (pour les
brèches suffisamment importantes) et enfin la phase d’injection
des accumulateurs avec remplissage de la cuve et renoyage du
5.2.2.3 Renoyage
cœur (le cas échéant). Le pressuriseur se vidange au début de la
Après 40 s, le système d’injection de sûreté basse pression phase de dépressurisation. L’inventaire en eau diminue en fonction
(environ 2 MPa) vient injecter de l’eau dans les boucles, en com- du débit critique à la brèche. La dépressurisation s’achève lorsque
plément des accumulateurs avant leur vidange complète. Cette la pression dans le circuit atteint la pression de saturation de la
eau va continuer à remplir la cuve et permettre un renoyage température de l’eau dans les parties chaudes. Pour les situations
stable de la partie basse du cœur. Le renoyage est un phénomène limites, la dégradation des conditions locales (apparition de vapeur
gravitaire : la progression du niveau d’eau dans le cœur se produit due à la baisse de pression) peut conduire à la crise d'ébullition
sous l’effet du poids de la colonne d’eau présente dans le collec- dans les zones du cœur aux plus fortes puissances.

BN 3 050v2 – 40
Lors de la dépressurisation du circuit primaire consécutive à la tionnel, dans la partie montante des tubes de GV. La fraction
brèche, le débit à la brèche est liquide. Un premier palier va être condensée, eau « claire » sans bore, redescend à contre-courant
atteint avec l’équilibre entre les pressions au primaire et au secon- vers la cuve, mais ne passe plus dans la partie descendante des
daire (figure 19), la création de vapeur dans le cœur est compen- tubes.
sée par les échanges aux générateurs de vapeur qui condensent L’inventaire en eau diminuant et la circulation de l’eau dans
la vapeur d'eau du primaire. La puissance résiduelle est évacuée l’ensemble du circuit primaire n’étant plus assurée (caloduc), le
principalement par le secondaire. niveau d’eau baisse dans la branche en U et dans la cuve (baisse
Tant que les pompes primaires continuent d’opérer, une circula- manométrique). La partie dénoyée du cœur (la plus haute) n’est
tion diphasique s’établit dans le cœur. L’écoulement est relative- plus suffisamment refroidie. La baisse de niveau se termine
ment homogène (niveau « gonflé » dans le cœur) ce qui favorise lorsque le liquide à la brèche passe à l’état vapeur une fois le
le refroidissement du cœur mais accélère la perte de l'inventaire « bouchon liquide » poussé de la branche en U par la vapeur
en eau dans le circuit primaire. venant de la sortie GV. L’inventaire en eau dépend alors du débit
à la brèche et du débit des injections de sécurité. La pression
Lorsque la production de vapeur dans le cœur est insuffisante,
continuant de diminuer, pour les brèches de plus grande taille, il y
la pression baisse alors au primaire (figure 19) et se retrouve en
a vaporisation dans la cuve et le downcomer, le cœur pouvant se
dessous de la pression au secondaire. Les échanges de chaleur
vidanger totalement. Finalement, le débit des injections de sécu-
sont alors inversés, le GV réchauffant le primaire. L’eau perdue à
rité peut amener au renoyage (ralenti cependant pour les plus
la brèche reste liquide tant que l’eau liquide présente dans la
grosses brèches par la vapeur qui s’évacue à la brèche à contre-
branche intermédiaire en amont de la brèche n’est pas expulsée à
courant du liquide dans le cœur et dans l’espace annulaire.
la fois par la vapeur et par la vidange à la brèche. Dès le passage
de la brèche en vapeur, (figure 19) la dépressurisation est signifi- Vis-à-vis de la thermohydraulique, les brèches en branche
cative et peut mener au découvrement du cœur (baisse manomé- froide sont plus pénalisantes que celles en branche chaude, l’eau
trique) compensé cependant par l’eau injectée des accumulateurs. injectée à la brèche étant perdue dans la boucle rompue, et la
perte en masse plus importante du fait d’un passage en phase

En cas d’arrêt des pompes primaires, le niveau de fluide dans le vapeur à la brèche plus tardif. On se place donc dans la suite dans
cœur se tasse, le fluide s’échauffe dans le haut du cœur et la l’hypothèse d’une brèche en branche froide.
vapeur se porte en sortie cœur, dans le plénum supérieur, dans
les branches chaudes et dans la partie montante des tubes des Outre la convection naturelle (thermosiphon et caloduc) en
générateurs de vapeur. Une circulation en thermosiphon simple phase liquide et en diphasique, les phénomènes thermohy-
(§ 4.5.7.1) simple phase puis diphasique s’établit, le GV étant drauliques prépondérants sont le contre-courant et le traitement
placé une petite dizaine de mètres au-dessus du cœur. La vidange des incondensables qui peuvent interrompre la circulation natu-
se poursuivant, de la vapeur passe en partie descendante des relle. La répartition du liquide et de la vapeur dans le circuit pri-
tubes de GV (générateur de vapeur) et, les pertes de charges aug- maire (séparation des phases) influence notablement le
mentant, le débit va en diminuant jusqu’à s’annuler lorsque déroulement de l’accident.
l’inventaire en eau devient insuffisant pour maintenir la circulation En cas d’APRP grosse brèche ou brèche intermédiaire, le phéno-
naturelle. mène de « steam binding » a de l’importance sur les délais de
Pour les plus grosses brèches intermédiaires, un point de sta- renoyage. Dans le GV (principalement mais aussi dans le plénum
gnation peut apparaître amenant possiblement à une crise d’ébul- supérieur ou encore les branches chaudes), l’eau liquide entraînée
lition localisée due à des vitesses locales faibles. Selon les par la vapeur peut se vaporiser et la perte de charge est accrue du
configurations, le primaire soit continue d'évacuer partiellement la fait de l’accélération de l’écoulement, le coefficient de frottement
puissance résiduelle au GV (température primaire supérieure à peut aussi être plus important.
celle du secondaire) soit la pression au primaire continue de bais- Enfin les phénomènes d’entraînement (ou de « dé-entraînement »)
ser de façon indépendante des GV, et les conditions d’injection des gouttes de liquide par la vapeur contribuent aux délais de
des accumulateurs puis de l’ISBP (entre 1 et 2 MPa) sont atteintes. renoyage et de refroidissement. En effet, lorsque les sections de
Le transport d’énergie aux GV se fait alors en thermosiphon passage varient, la vitesse de la vapeur (et sa capacité à emporter
local double phase apparenté à un caloduc, la vapeur sortant de la les gouttes de liquide) varie. Par exemple le liquide entraîné par la
cuve allant se condenser, si le refroidissement au GV est opéra- vapeur dans le cœur va rester dans le plénum supérieur où la
vitesse de la vapeur diminue fortement (section de passage bien
plus grande), ce liquide va pouvoir ensuite retomber en périphérie
du cœur. En outre, dans la phase de renoyage, l’abaissement de la
pression en haut du cœur (par condensation de la vapeur avec les
Pression gouttes de liquide entraînées) va abaisser la pression et permettre
la remontée du liquide venant du downcomer.
155 Enfin la montée en température des crayons peut conduire à
130 AAR l’oxydation et la déformation (ballonnement) pouvant s’accompa-
120 IS gner d’une relocalisation du combustible.
5.2.4 Accident grave

86 Les scenarii présentés précédemment concernent la phase court
72 terme de l’APRP. À plus long terme, les opérateurs de la centrale
ont un certain nombre d’actions à réaliser afin de maintenir le
réacteur dans un état sûr. Si l’on suppose que les opérateurs
n’arrivent pas à intervenir (comme dans des situations telles que
l’accident de Fukushima), le refroidissement du cœur du réacteur
Pression primaire Temps
peut à nouveau s’interrompre. À la fin de la phase court terme de
Pression secondaire
l’APRP, le refroidissement du cœur est assuré par le système
d’injection de secours, ce dernier puisant l’eau dans une bâche
Figure 19 – Évolution des pressions primaire et secondaire dans le (réservoir) qui a une contenance limitée. L’eau perdue à la brèche
cas d’un APRP brèche intermédiaire se retrouve cependant au sol du bâtiment réacteur. Cette eau peut

BN 3 050v2 – 41
alors être pompée et le refroidissement peut se poursuivre. des phénomènes est multidimensionnelle et constitue une étape
Cependant, en cas de dysfonctionnement, l’absence de refroidis- essentielle pour parvenir à un haut niveau de modélisation. Cela
sement conduit à l’augmentation de la température du cœur. La demande des moyens d’analyse sophistiqués, moyens de calcul et
puissance résiduelle produite par les crayons combustibles reste moyens expérimentaux. Bien que les capacités de calcul aient for-
effective sur le très long terme. Il s’ensuit alors une dégradation tement progressé lors des dernières décennies, la base de
du cœur avec possible fusion des gaines et des pastilles de com- connaissance s’appuyant sur des interprétations d’essais requiert
bustibles (corium fondu dans le cœur). On parle d’accident grave encore des développements et fait l’objet de programmes en pro-
lorsque la dégradation du cœur commence. jet ou en cours de réalisation. Les besoins du nucléaire sont nom-
Ce processus de dégradation du cœur est caractérisé par des breux et doivent plus particulièrement appuyer les applications
phénomènes complexes s’accompagnant de changement de géo- liées aux accidents graves et aux systèmes innovants notamment
métrie dans le cœur et dans la cuve du réacteur (relocalisation et de sûreté passifs.
fusion du combustible). Par ailleurs, lors de tels accidents, plu-
sieurs réactions chimiques entre la gaine et la vapeur d’eau ou
l’air sont à prendre en considération car elles ont un impact non
négligeable sur le déroulement de l’accident. Par exemple, la réac- 7. Glossaire
tion d’oxydation de la gaine par l’eau conduit à de multiples
conséquences, à savoir une fragilisation de la gaine, la formation
AAR
d’hydrogène (gaz explosif) et une accélération de la montée en
température (réaction exothermique). En outre, on peut craindre Arrêt automatique du réacteur
l’explosion vapeur suite à la vaporisation immédiate du liquide au APRP
contact du combustible fondu. Un code de thermohydraulique
seul pour prendre en compte tous ces phénomènes n’est donc Accident par perte de réfrigérant primaire
plus suffisant, c’est pourquoi des couplages entre plusieurs codes DNB
(notamment avec des codes de thermomécanique) sont réalisés. Departure from nucleate boiling
Il existe deux phases d’étude successives dans l’accident grave, GRS
la première dédiée au comportement du combustible jusqu’à sa
fusion et la seconde dédiée au comportement du combustible Gesellschaft für Anlagen- und Reaktorsicherheit mbH
fondu. Lors de la première phase, les phénomènes physiques GV
importants concernent la montée en température des crayons,
Générateur de vapeur
l’oxydation et la déformation (ballonnement) puis la rupture des
gaines du combustible s’accompagnant d’une relocalisation dans HEM
la gaine et dans la cuve de fragments de pastilles combustibles. Homegenous equilibrium model
Lors de la seconde phase de l’accident, les phénomènes phy-
siques à considérer sont relatifs au comportement du corium et IRSN
ses interactions avec les composants non fondus dans le cœur, Institut de radioprotection et de sûreté nucléaire
lors de son déplacement et son accumulation dans le fond de IS
cuve, pouvant entraîner une rupture de la cuve. Le code ASTEC
(Accident Source Term Évaluation Code) co-développé par l’IRSN Injection de sécurité
et GRS [55] est en particulier dédié à l’étude des accidents graves, ISBP, ISMP, ISHP
la référence [56] présente un état des lieux détaillé des approches Injection de sécurité basse pression, moyenne pression, haute
et des outils dédiés aux accidents graves. pression
ONB
Onset of nucleate boiling
6. Conclusion OSV
Onset of significant void
La thermohydraulique doit relever les défis de la mécanique des
REP
fluides à savoir parvenir à des relations de fermeture modélisant
au mieux les phénomènes. Cependant, la coexistence de plusieurs Réacteur à eau pressurisée
phases dans l’écoulement rend la tâche hautement complexe RFTC
compte tenu des interactions multiples entre les phases et des
interfaces évoluant sans cesse. En l’état actuel des connaissances, Rapport de flux thermique critique
suivant une approche moyennée à une dimension, des modèles
satisfaisants ont pu être établis et intégrés aux outils de calculs
scientifiques. Ces outils de calculs ont atteint un niveau de
confiance et de robustesse qui reproduisent assez fidèlement les 8. Sigles, notations
très nombreuses expériences réalisées en support au développe-
ment des centrales nucléaires. Cependant, la compréhension fine et symboles
Dimension Ordre
Unité
Notation Définition Unité SI de grandeur Équation Facteur de conversion (2)
anglaise
(M L Θ t) (1) (SI)
A Surface de passage m2 ft2 L2 1ft2 = 0,092903 m2
cp Capacité calorifique à kJ/kg–°C Btu/lbm–°F L2t–2Θ–1 4,5 1 Btu/lb°F = 4,1868 kJ/kg°C

pression constante

BN 3 050v2 – 42
Dimension Ordre
Unité
anglaise (1)
(M L Θ t) (SI)
c Vitesse du son m/s ft/s L t–1 tableau 12 1ft/s = 0,3048 m/s
C Paramètre d’écoulement (–) (–) (–) rela- C = α/β (–)

tion (54)
C0 Constante du modèle à (–) (–) (–) tableau 3 (–)

dérive de flux
Cx Coefficient de Fanning (–) (–) (–) 0,0005 (–)
D Diamètre des crayons m in L 0,0095 1 in = 0,0254 m

chauffants
Dh Diamètre (hydraulique) m ft L 0,0010 Dh = 4A/ 1ft/s = 0,3048 m/s

–1
d Vitesse de dérive m/s ft/s Lt d = Uv – j 1ft/s = 0,3048 m/s
drift velocity
dfv Dérive de flux m/s ft/s L t–1 dfv = αd = jv – αj 1ft/s = 0,3048 m/s
« drift flux »
f Coefficient de perte de (–) (–) (–) 0,02 f = 4Cx (–)

charge
g Accélération pesanteur m2/s ft/s2 L2t–1 9,81 1 ft/s2 = 0,3048 m/s2
G Vitesse massique kg/m2 – s 106lb/ft2 – hr ML–2t–1 3 500 G = W/A = ρU 1 Mlb/h – ft2 = 1356,23 kg/m2 s
h Enthalpie kJ/kg Btu/lbm L2t–2 1 200 Δh≈cpΔT 1 Btu/lb = 2,326 kJ/kg
hc Coefficient de transfert W/m2–°C Btu/h ft2 °F M Θ–1 t–3 50 000 h = φ"/ΔT 1 Btu/h–ft2 °F = 5,678 W/m2
thermique par °C
convection
hlv Chaleur latente kJ/kg Btu/lbm L2t–2 1 000 1 Btu/lb = 2,326 kJ/kg
j Flux volumétrique ou m/s ft/s L t–1 j = Q/A 1ft/s = 0,3048 m/s

vitesse superficielle
jv Flux volumétrique m/s ft/s L t–1 jv = Gx/ ρv 1ft/s = 0,3048 m/s

vapeur
jl Flux volumétrique m/s ft/s L t–1 jl = Ql/ A 1ft/s = 0,3048 m/s

liquide
m Masse kg lbm M t–1 1 lb = 0,4535924 kg
Débit massique kg/s lbm/h M t–1 1 lb/h = 0,000126 kg/s
P Pression Pa = N/m2 psia = lbf/in² M L–1 t–2 15,5 106 1 psi = 6894,76 Pa
pas Pas entre crayons m in L 0,0126 1 in = 0,0254 m
P, Π Puissance thermique W = J/s Btu/hr ML2t–3 1 Btu/hr = 0,2931 W
Périmètre mouillé m in L 1 in = 0,0254 m
q’ Puissance linéique W/m Btu/h ft MLt–3 20 000 1 Btu/h ft = 0,9615 W/m
q’’ φ Flux de chaleur W/m² Btu/h ft2 Mt–3 600 000 1 Btu/h ft2 = 3,1546 W/m2

BN 3 050v2 – 43
Dimension Ordre
Unité
anglaise (1)
(M L Θ t) (SI)
Q Débit volumique m3/s ft3/h (3) L3 t–1 1 USGPM = 0,2271247 m3/h
Q Chaleur W = J/s Btu/hr ML2 t–3 1 Btu/hr = 0,2931 W
R Rayon des crayons m in L 1 in = 0,0254 m

chauffants
S Glissement (–) (–) (–) 1 à 10 S = Uv/Ul (–)
T Température °C °F Θ 300 0°C = 32 °F
°F = x1,8°C+32
Vitesse m/s ft/s L t–1 5,35 U = Q / A 1ft/s = 0,3048 m/s

U,
Ud Vitesse débitante m/s ft/s L t–1 1ft/s = 0,3048 m/s
ud Vitesse de dérive m/s ft/s L t–1 tableau 4 vd = <αd>/<α> 1ft/s = 0,3048 m/s

moyenne
x Qualité (–) (–) (–) 0,0 x = Wv/Wtotal (–)
V Volume m3 US Gal =231 L3 1 US gal = 0,0037854 m3

in.3
v Volume spécifique m3/kg ft3 /lb M L–3 v = 1/ρ 1 ft3/lb = 0,062428 m3/kg
W Travail W = J/s Btu/hr ML2 t–3 1 Btu/hr = 0,2931 W
α Taux de vide (–) (–) (–) 0à1 α = Vv/Vtotal (–)
β Taux de débit volumé- (–) (–) (–) β = Qv/Qtotal (–)

trique vapeur
β Coefficient de dilatation (–) (–) (–) tableau 20 (–)
λ Conductivité thermique W/m–°C Btu/h ft °F ML Θ–1 t–3 0,6 1 Btu/h ft °F = 1,730735 W/m

°C
τ Frottement N/m2 ML–1t–2 τ = μdU/dx

–1 –1 –3
μ Viscosité dynamique kg/m – s lbm/ft – s ML t 10 = 10 poise 1 lbm/ft – s = 1,48816 Pa s
(poi-
seuille)
ν Viscosité cinématique m2/s ft2/s L2t–1 10–6 = μ/ρ 1 ft2/s = 0,0929 m2/s
(myriasto
kes)
κ Diffusivité thermique m2/s ft2/s L2t–1 κ = λ/ρ cp 1 ft2/s = 0,0929 m2/s
η Rendement (–) (–) (–)
χ Titre thermodynamique (–) (–) (–)
ρ Masse volumique kg/m3 lbm/ft3 M L–3 700 1 lb/ft3 = 16,018 kg/m3
σ Tension superficielle N/m lbf/in M t–2 0,05 1 lbf/in = 175,127 N/m
φ ou q'' Flux surfacique W/m2 Btu/h – ft2 Mt–3 600 000 1 Btu/h ft2 = 3,1546 W/m2
BO Nombre d’ébullition (–) (–) (–) φ/(Ghlv) (–)

BN 3 050v2 – 44
Dimension Ordre
Unité
anglaise (1)
(M L Θ t) (SI)
Fr Nombre de Froude (–) (–) (–) U2/gD (–)
Gr Nombre de Grashof (–) (–) (–) gL3βΔT/ν2 (–)
Ja Nombre de Jakob (–) (–) (–) cpΔT/hlv (–)
Nu Nombre de Nusselt (–) (–) (–) hD/λ (–)
Pe Nombre de Peclet = Re Pr (–) (–) ρcpUD/λ (–)
Pr Nombre de Prandtl (–) (–) (–) 1 ν/κ = μcp/λ (–)

3
Ra Nombre de Raleigh = Gr Pr (–) (–) (gL βΔT)/(νκ) (–)
Re Nombre de Reynolds (–) (–) (–) 500000 ρUD/μ (–)
St Nombre de Stanton = Nu/ Re (–) (–) h/ρUcp (–)

Pr
(1) M : masse L : longueur Θ : température t : temps.

L’unité de force (newton) est en ML t–2. L’unité d’énergie (joule = Nm) est en ML2 t–2. L’énergie spécifique (par unité de masse) est en (L/t)2.
(2) On obtient les unités SI en multipliant les valeurs en unités anglaises par le facteur indiqué
1 ft = 12" = 12* 2,54/100 = 0,3048 m.
(3) On trouve aux États-Unis l’US gal/min soit 6,309 10–5 m3/s.
Indices et exposants Indices et exposants
l Liquide he Hors équilibre
e Équilibre ou condition en entrée sat Saturation
v Vapeur cr Critique

BN 3 050v2 – 45
P
O
U
Thermohydraulique des réacteurs R
à eau sous pression
E
N
Ingénieur chargé d’expertises
Institut de Radioprotection et de Sûreté Nucléaire, Fontenay-aux-Roses, France
S
Avec la participation de Tony Glantz (IRSN) pour la rédaction des sous-sections 3.4 et 5.2.
A
V
Sources bibliographiques
O
[1] COPPOLANI (P.) et al. – La chaudière des ré- http://www.osti.gov/scitech/servlets/purl/ [29] GEFFRAYE (G.) et al. – CATHARE 2 V2.5 2 : A
acteurs à eau sous pression. EDP Sciences
(2004). [15]
6871478/
CODDINGTON (P.) et MACIAN (R.). – A study
single version for various applications
Nuclear Engineering and Design 241 4456-
4463 (2011).
I
[2] WAGNER (W.) et KRETZSCHMAR (H.-J.). – In- of the performance of void fraction correla-
ternational Steam Tables – Properties of Wa-
ter and Steam based on the Industrial Formu-
tions used in the context of drift-flux two-
phase flow models. Nuclear Engineering and
[30] AIEA. – Passive Safety Systems and Natural
Circulation in Water Cooled Nuclear Power
R
lation IAPWS-IF97. Springer-Verlag (2008). Design 215, 199-216 (2002). Plants TECDOC-1624 (2012).
[3] TODREAS (N.) et KAZIMI (M.). – Nuclear Sys- [16] BESTION (D.). – The physical closure laws in [31] EMONOT (P.), SOUYRI (A.), GANDRILLE
tems : Thermal Hydraulic Fundamentals, Se- the CATHARE code. Nuclear Engineering and (J.L.) et BARRÉ (F.). – CATHARE 3 : a new sys-
[4]
cond Edition, Volume 1. CRC Press (2011).
DELHAYE (J.-M.). – Thermohydraulique des
réacteurs. EDP Sciences (2013).
[17]
Design 124, 229-245 (1990).
WOLDESEMAYAT (M.A.) et GHAJA (A.J.). –
tem code for thermal-hydraulic in the context
of the NEPTUNE project. Nuclear Enginee-
ring and Design 241 4476-4481 (2011).
P
Comparison of void fraction correlations for
[5] COLLIER (J.G.) et THOME (J.R.). – Convective
Boiling and Condensation. Clarendon Press
different flow patterns in horizontal and
upward inclined pipes. International Journal
[32] LEE (S.C.), BANKOFF (S.G.). – A comparison
of predictive models for the onset of signifi-
L
(1996). of Multiphase Flow 33, 347-370 (2007). cant void at low pressures in-forced-convec-
[6] THOME (J.R.) et al. – Encyclopedia of Two-
Phase Heat Transfer and Flow. World Scien-
[18] XU (Y.) et FANG (X.) et al. – Evaluation of fric-
tional pressure drop correlations for two-
tion subcooled boiling. KSME Int. 5. 12.504-
513 -1998)
U
tific Publishing (2015). phase flow in pipes Nuclear Engineering and [33]
[7] GHIAASIAAN (S.M.). – Two-Phase Flow, Boi-
ling, and Condensation In Conventional and [19]
Design 253, 86-97 (2012).
IDEL’CIK (I.E.). – Mémento des pertes de
RELAP5-3D Code Development Team. – Code
Manual Volume I : Code Structure, System,
Models, And Solution Methods INEEL-EXT-
S
Miniature Systems. Cambridge University charges Eyrolles (1984). 98-00834 rev2.4 (2005).
Press (2011). [20] IDEL’CIK (I.E.). – Handbook of Hydraulic Re- [34] AIEA. – Principes fondamentaux de sûreté.
[8] LEVY (S.). – Two-Phase Flow in Complex Sys- sistance Begell (2008). SF-1 (2007). http://www-pub.iaea.org/MTCD/
tems. Wiley (1999). [21] HOLMAN (J.). – Heat Transfer Mac Graw Hill publications/PDF/P1273_F_web.pdf
[9] ZUBER (F.). – The effects of non-uniform flow (2009). [35] TARRIDE (B.). – Physique, fonctionnement et
and concentration distributions and the ef- [22] BAILLY (H.), MENESSIER (D.) et PRUNIER sûreté des REP Maîtrise des situations acci-
fect of the local relative velocity on the ave- (C.). – Le combustible nucléaire des réacteurs dentelles du système réacteur. EDP Sciences
rage volumetric concentration in two-phase à eau sous pression et de réacteurs à neu- (2013).
flow. GEAP-4542, General Electric Co. (1964). trons rapides Eyrolles Collection du Commis- [36] ABRAMSON (P.B.). – Guidebook to light wa-
[10] ISHII (H.). – Thermo-Fluid Dynamics of Two- sariat à l’énergie atomique (1996). ter reactor safety analysis Hemisphere Publi-
Phase Flow. Springer (2011). shing Corporation (1985).
[23] THOM (J.R.S.). – Prediction of pressure drop
[11] CHEXAL, HOROWITZ et LELLOUCHE. – An during forced convection boiling of water Int. [37] BURCHILL (W.E.). – Physical phenomena of
Assessment of Eight Void Fraction Models J. Heat Mass Transfer 7,709-724 (1964). small-break loss-of-coolant accident in PWR
for Vertical Flows. NSAC-107 EPRI (1986) Nucl. Safety, 23 (No. 5) (1982).
[24] ZOHURI (B.) et FATHI (N.). – Thermal-Hydrau-
http://www.epri.com/abstracts/Pages/Pro-
ductAbstract.aspx?ProductId=NSAC-107
lic Analysis of Nuclear Reactors Springer [38] WINTERTON (R.H.S.). – Where did the Dittus
(2015). and Boelter equation come from? Int. J. Heat
[12] CHEXAL et LELLOUCHE. – A Full-Range Drift-
[25] WHALLEY (P.B.). – Boiling, Condensation, Mass Transfer Vol. 41, Nos 4-5, pp. 809-810
Flux Correlation for Vertical Flows (Revision (1998).
and Gas-Liquid Flow Clarendon Press Oxford
1). NP-3989-SR EPRI (1986)
Engineering Science Series (1990). [39] CHEN (J.C.), SUNDARAM (R.K.) et OZKAY-
http://www.epri.com/abstracts/Pages/Pro-
ductAbstract.aspx?ProductId=NP-3989-SR [26] GROENEVELD (D.C.). – An investigation of NAK (F.T.). – Phenomenological Correlation
heat transfer in the liquid deficient regime for Post-CHF Heat Transfer USNRC Report.
[13] CHEXAL et al. – Void Fraction Technology for
Chalk River AECL-3281 (1969). NUREG-0237 (1977).
Design and Analysis. TR-106326 EPRI (1997) :
http://www.epri.com/abstracts/Pages/Pro- [27] AIEA. – Thermohydraulic relationships for [40] CENGEL (Y.A.). – Introduction to Thermody-
ductAbstract.aspx?ProductId=TR-106326 advanced water cooled reactors TECDOC- namics and Heat Transfer Mc Graw Hill
[14] ISHII et MAMORU. – One-dimensional drift- 1203 (2001). (2007).
flux model and constitutive equations for re- [28] AWAD (M.) et MUZYCHKA (Y.). – Effective [41] AYED (R.) et BOURTEELE (J.-P.). – Theoreti-
lative motion between phases in various property models for homogeneous two- cal Evaluation of Currently Used Critical Heat
two-phase flow regimes. Argone National phase flows Exp. Therm. Fluid Sci. 33, 106- Flux Correlations Trans. ANS vol.59, Atlanta
Laboratory ANL-77-47 (1977) 113 (2008). (1989).

Doc. BN 3 050v2 – 1
P THERMOHYDRAULIQUE DES RÉACTEURS À EAU SOUS PRESSION ____________________________________________________________________________

O
U [42] FILLION (P.) et al. – FLICA-OVAP : A new plat- [48] D’AURIA (F.). – Thermal-Hydraulics of Water [54] GROENEVELD (D.C.) et al. – The 2006 CHF
R form for core thermal – hydraulic studies

Nuclear Engineering and Design 241 4348-
4358 (2011).
Cooled Nuclear Reactors Woodhead Publi-
shing (2017).
look-up table Nuclear Engineering and De-
sign 237 1909-1922 (2007).
[49] CAPPON (H.), DUFEIL (P.) et VEILLEUX (L.). – [55] Collectif. – ASTEC, un code intégré de simu-
[43] GABRIEL (H.) et GLÜCK (M.). – Prediction of Integration of the CATHARE 2 Safety Analy- lation d’accident grave dans les réacteurs à
void fraction for PWR and BWR conditions sis Code on the SOFIA Engineering Simula- eau légère Rapport scientifique et technique
with the sub-channel code F-COBRA-TF tors Proceedings of 2010 PowerPlantSim de l’IRSN (2002).
E Nuclear Engineering and Design 241 3952-
3966 (2011).
Conference, San Diego (USA), 21-26, Fe-
bruary 2010.
http://www.irsn.fr/FR/Larecherche/publica-
tions-documentation/aktis-lettre-dossiers-
[50] SWINDLEHURST (G.B.). – Data Sources for thematiques/RST/RST-2002/Documents/
N [44] FIGHETTI (C.F.) et REDDY (D.G.). – Parametric
Study of CHF Data, Volume 3, Part 1 : Critical
Heat Flux Data NP-2609-V3P1, EPRI (1982)
Capability Assessments of Next Generation
Safety Analysis Codes EPRI, Palo Alto [56]
Chap01_art5.pdf
AIEA. – Approaches and Tools for Severe
http://www.epri.com/abstracts/ Pages/Pro- 3002003110 (2014) Accident Analysis for Nuclear Power Plants
http://www.epri.com/abstracts/Pages/Pro- Safety Report Series-56 (2008)
ductAbstract.aspx?ProductId=NP-2609-V3P1
ductAbstract.aspx?Produc- http://www-pub.iaea.org/MTCD/publications/
[45] HERER (C.). – Conservatisms, provisions and
S margins in DNB analyses, Proceedings of
2016 ICAPP Conference, San Francisco
[51]
tId=000000003002003110
Collectif. – TRACE V5.0 Theory Manual
ML120060218 US Nuclear Regulatory Com-
[57]
PDF/Pub1327_web.pdf
AIEA. – Thermophysical properties database
of materials for light water reactors and
(USA), 17-20, April 2016.
A [46] OWEN (D.B.). – Factors for one-sided tole-
rance limits and for variable sampling plans
mission (2008)
http://www.nrc.gov/docs/ML1200/
ML120060403.html
heavy water reactors IAEA-TECDOC-1496
(2006)
[58] CARBAJO (J. J.), et al. – A review of the ther-
V SCR-607 Sandia Corporation (1963)
http://www.nrc.gov/docs/ML1403/
ML14031A495.pdf
[52] GROENEVELD (D.C.) et al. – 1986 AECL-UO
Critical Heat Flux Lookup Table Heat Transfer
mophysical properties of MOX and UO2
fuels Journal of Nuclear Materials 299, 3
Engineering Volume 7, Issue 1-2 (1986). (2001)
O [47] [53]
WILKS (S.S.). – Determination of Sample GROENEVELD (D.C.) et al. – The 1995 look-
Sizes for Setting Tolerance Limits Ann. Math. up table for critical heat flux in tubes Nuclear
Statist. 12 (1941). Engineering and Design 163 1-23 (1996).
I
R À lire également dans nos bases
LALLEMAND (A.). – Écoulement des fluides – Écou- LALLEMAND (M.). – Transferts en changement de PANDEY (P.). – Transfert en changement de
lements en conduites. Réseaux. [BE8161], phase – Ébullition libre. [BE8235], Transferts phase – Condensation sur des surfaces lisses.
P Fluides et combustion (2015).

LALLEMAND (A.). – Écoulement des fluides – Dy-
namique des fluides réels. [BE8157], Fluides et
thermiques (2005).
LALLEMAND (M.). – Transferts en changement de
phase – Ébullition convective. [BE8236], Trans-
[BE8238], Transferts thermiques (2006).
BORZOV, DIMITROV et HERER. – Réacteurs à eau
sous pression de conception russe : les VVER.
L combustion (2015).
LALLEMAND (A.). – Thermodynamique appli-
ferts thermiques (2006).
PADET (J.). – Transferts en changement de phase –
[BN3101], Génie nucléaire (2016).
BLANPAIN, GARNIER, IMBERT et SALAUN. – Ther-
quée – Premier principe. Énergie. Enthalpie. Ébullition con Convection thermique et mas- momécanique du combustible des réacteurs à
U [BE8005], Thermodynamique et énergétique
(2015).
sique – Principes généraux. [BE8205], Trans-
eau sous pression. [BN3060], Génie nucléaire
(2011).
GOSSE (J.). – Mécanique des fluides. [A1870], PADET (J.). – Convection thermique et massique – LIBMANN (J.). – Sûreté des réacteurs français en
S Sciences fondamentales (1995). Nombre de Nusselt : partie 1. [BE8206], Trans-
service. [BN3820], Génie nucléaire (2000).
Sites Internet
Cours en accès libre du MIT (2007) Cours en accès libre de l’AIEA (2008).
http://ocw.mit.edu/courses/nuclear-engineering/22-313j-thermal-hydraulics- Fundamentals and Basic Professional Training Courses : Thermal Hydrau-
in-power-technology-spring-2007/ lics, Heat Transfer
Buongiorno, Jacopo. 22.313J Thermal Hydraulics in Power Technology. http://nucleus.iaea.org/Pages/thermal-hydraulics-heat-transfer.aspx
Spring 2007. Massachusetts Institute of Technology : MIT OpenCourseWare
http://ocw.mit.edu.
Événements
Conférence NURETH (NUclear REactor THermalhydraulics) tous les deux ans.
Réglementation
Loi n° 2006-686 du 13 juin 2006 relative à la Transparence et à la sécurité Guide de l’ASN n° 22 : Conception des réacteurs à eau sans pression
en matière Nucléaire (TSN). (2017).
Arrêté du 7 février 2012 fixant les règles générales relatives aux installa- Guide de l’ASN n° 28 : Qualification des outils de calcul scientifiques utili-
tions nucléaires de base (INB). sés dans la démonstration de sûreté nucléaire – première barrière (2017).
Loi n° 2015-992 du 17 août 2015 relative à la transition énergétique pour la
croissance verte.

Doc. BN 3 050v2 – 2
GAGNEZ DU TEMPS ET SÉCURISEZ VOS PROJETS
EN UTILISANT UNE SOURCE ACTUALISÉE ET FIABLE
Techniques de l’Ingénieur propose la plus importante

collection documentaire technique et scientifique
en français !
Grâce à vos droits d’accès, retrouvez l’ensemble
des articles et fiches pratiques de votre offre,
leurs compléments et mises à jour,
et bénéficiez des services inclus.
   
RÉDIGÉE ET VALIDÉE MISE À JOUR 100 % COMPATIBLE SERVICES INCLUS
PAR DES EXPERTS PERMANENTE SUR TOUS SUPPORTS DANS CHAQUE OFFRE
NUMÉRIQUES
 + de 350 000 utilisateurs

 + de 10 000 articles de référence
 + de 80 offres
 15 domaines d’expertise
Automatique - Robotique Innovation
Biomédical - Pharma Matériaux
Construction et travaux publics Mécanique
Électronique - Photonique Mesures - Analyses
Énergies Procédés chimie - Bio - Agro
Environnement - Sécurité Sciences fondamentales
Génie industriel Technologies de l’information
Ingénierie des transports
Pour des offres toujours plus adaptées à votre métier,

découvrez les offres dédiées à votre secteur d’activité
Depuis plus de 70 ans, Techniques de l’Ingénieur est la source

d’informations de référence des bureaux d’études,
de la R&D et de l’innovation.
www.techniques-ingenieur.fr
CONTACT : Tél. : + 33 (0)1 53 35 20 20 - Fax : +33 (0)1 53 26 79 18 - E-mail : [email protected]
LES AVANTAGES ET SERVICES
compris dans les offres Techniques de l’Ingénieur
  
ACCÈS
Accès illimité Téléchargement des articles Consultation sur tous

aux articles en HTML au format PDF les supports numériques
Enrichis et mis à jour pendant Pour un usage en toute liberté Des contenus optimisés
toute la durée de la souscription pour ordinateurs, tablettes et mobiles
 
SERVICES ET OUTILS PRATIQUES
Questions aux experts* Articles Découverte Dictionnaire technique multilingue

Les meilleurs experts techniques La possibilité de consulter des articles 45 000 termes en français, anglais,
et scientifiques vous répondent en dehors de votre offre espagnol et allemand
 
Archives Impression à la demande Alertes actualisations
Technologies anciennes et versions Commandez les éditions papier Recevez par email toutes les nouveautés
antérieures des articles de vos ressources documentaires de vos ressources documentaires
*Questions aux experts est un service réservé aux entreprises, non proposé dans les offres écoles, universités ou pour tout autre organisme de formation.
ILS NOUS FONT CONFIANCE
www.techniques-ingenieur.fr
CONTACT : Tél. : + 33 (0)1 53 35 20 20 - Fax : +33 (0)1 53 26 79 18 - E-mail : [email protected]

bn3050

Transféré par

Droits d'auteur :

Formats disponibles

bn3050

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

bn3050

Transféré par

Droits d'auteur :

Formats disponibles

Réf.

Cet article est issu de : Énergies | Génie nucléaire

par Christophe HERER

Pour toute question :

Thermohydraulique des réacteurs

1. Notions élémentaires .............................................................................. BN 3 050v2 - 2

2. Caractérisation de l’état du fluide .......................................................... — 3

L a thermohydraulique est l’étude du transfert de la chaleur couplé à

Copyright © – Techniques de l’Ingénieur – Tous droits réservés

THERMOHYDRAULIQUE DES RÉACTEURS À EAU SOUS PRESSION ____________________________________________________________________________

Une des spécificités de l’industrie nucléaire et du réacteur à eau

Si le concept du diamètre hydraulique est acceptable, son appli-

Copyright © – Techniques de l’Ingénieur – Tous droits réservés

_____________________________________________________________________________ THERMOHYDRAULIQUE DES RÉACTEURS À EAU SOUS PRESSION

les tubes-guide des assemblages combustible. Du fait de la loi

Copyright © – Techniques de l’Ingénieur – Tous droits réservés.

THERMOHYDRAULIQUE DES RÉACTEURS À EAU SOUS PRESSION ____________________________________________________________________________

Tableau 1 – Grandeurs thermodynamiques de l’eau liquide

15,5 292,2 1 296 2,585 742 91,6

15,5 329,6 1 514 4,774 653 76,3

7,3 229,5 989 1,622 832 117

Tableau 2 – Grandeurs thermodynamiques de l’eau à saturation

Température Masse Masse Tension de

2,0 212,4 850 10,0 908,6 2 798 0,03483

7,3 288,7 734 38,3 1 283 2 769 0,01697

10,0 311,0 688 55,5 1 408 2 725 0,01186

15,5 344,8 594 102 1 630 2 596 0,004669

Les tables d’eau [2] sont maintenant disponibles sous de nom-

À noter que s’il n’y a pas de changement de phase et pour des

La quantité de chaleur par unité de temps correspond à la puis-

Copyright © – Techniques de l’Ingénieur – Tous droits réservés

_____________________________________________________________________________ THERMOHYDRAULIQUE DES RÉACTEURS À EAU SOUS PRESSION

centrale produisant une puissance de 1 400 MW avec un rendement

– fluide à l’état de vapeur surchauffée : (9)

(6) 2.2.1.1 Débits et vitesses

La vitesse moyenne de l’écoulement est le rapport du débit

Copyright © – Techniques de l’Ingénieur – Tous droits réservés.

THERMOHYDRAULIQUE DES RÉACTEURS À EAU SOUS PRESSION ____________________________________________________________________________

Témpérature moyenne du fluide

Monophasique liquide Écoulement

Figure 3 – Écoulement dans un tube chauffé

Le débit volumique étant rapporté à la surface totale, la relation (22)

Copyright © – Techniques de l’Ingénieur – Tous droits réservés

_____________________________________________________________________________ THERMOHYDRAULIQUE DES RÉACTEURS À EAU SOUS PRESSION

2.3 Régimes d’écoulements

Copyright © – Techniques de l’Ingénieur – Tous droits réservés.

THERMOHYDRAULIQUE DES RÉACTEURS À EAU SOUS PRESSION ____________________________________________________________________________

plus probable, puis à vérifier la cohérence des résultats obtenus et

0,60 – dans des écoulements verticaux, lorsque la phase vapeur

7 MPa S = 1 15,5 MPa S = 1 soit , la corrélation de Wallis éta-

avec m = 1 et C = 1 lorsque l’effet de la gravité est prépondérant

Les grandeurs C et m de la relation de Wallis ci-dessus peuvent

Les propriétés physiques de l’eau à saturation sont les suivantes

Les grandeurs thermohydrauliques suivantes sont calculées :

La relation avec C = 1 et m = 1 conduit à :

Comme indiqué au paragraphe 2.1.3, l’équilibre thermodyna-

Copyright © – Techniques de l’Ingénieur – Tous droits réservés

_____________________________________________________________________________ THERMOHYDRAULIQUE DES RÉACTEURS À EAU SOUS PRESSION

3. Modélisations Quelles que soient les corrélations, lois de fermeture ou