Analyse Thermohydraulique Transitoire Dans Un Canal Du Coeur D'un Réacteur Nucléaire de Recherche de Type MTR

RÉ PUBLIqUE ALGÉ RIEnnE DÉ MOCRATIqUE ET POPULAIRE
mInISTÈ RE DE L’EnSEIGnEMEnT SUPÉ RIEUR ET DE LA RECHERCHE SCIEnTIfIqUE
École Nationale Polytechnique
Département d’Hydraulique
Laboratoire de Recherches Sciences de l’Eau
Projet de fin d’études

Mémoire de projet de fin d’études pour l’obtention du diplôme d’ingénieur d’état en
Hydraulique
Analyse thermohydraulique
transitoire dans un canal du coeur
d’un réacteur nucléaire de recherche
de type MTR
Réalisé par :
BENSEMMANE Yasmine Salima
TIHAL Ouarda
Soutenu publiquement le 11/07/2021 devant le jury composé de :
Président Dr TACHI Salah Eddine MCA ENP
Encadreur Dr SIDI ALI Kamel Chercheur Expert CRND
Co-Encadreur Dr BENMAMAR Saâdia Pr ENP
Examinateur Dr BELGACEM Ibtissem MCA ENP
Mr. BELHADJ Omar MAA ENP
Mr. OUCHENE Samir Doctorant G.M. ENP
ENP 2021
Annexe D. Étapes de simulation FLUENT en régime transitoire
Mémoire de PFE TIHAL et BENSEMMANE 150

RÉ PUBLIqUE ALGÉ RIEnnE DÉ MOCRATIqUE ET POPULAIRE
mInISTÈ RE DE L’EnSEIGnEMEnT SUPÉ RIEUR ET DE LA RECHERCHE SCIEnTIfIqUE
École Nationale Polytechnique
Département d’Hydraulique
Laboratoire de Recherches Sciences de l’Eau
Projet de fin d’études

Mémoire de projet de fin d’études pour l’obtention du diplôme d’ingénieur d’état en
Hydraulique
Analyse thermohydraulique
transitoire dans un canal du coeur
d’un réacteur nucléaire de recherche
de type MTR
Réalisé par :
BENSEMMANE Yasmine Salima
TIHAL Ouarda
Soutenu publiquement le 11/07/2021 devant le jury composé de :
Président Dr TACHI Salah Eddine MCA ENP
Encadreur Dr SIDI ALI Kamel Chercheur Expert CRND
Co-Encadreur Dr BENMAMAR Saâdia Pr ENP
Examinateur Dr BELGACEM Ibtissem MCA ENP
Mr. BELHADJ Omar MAA ENP
Mr. OUCHENE Samir Doctorant G.M. ENP
ENP 2021
Abstract
In order to ensure the safety of nuclear research reactors, the determination of their
thermal-hydraulic characteristics requires a permanent monitoring of some physical
quantities such as pressure, temperature and heat fluxes which prevail in their core .
Thus, to have a correct evolution of these physical quantities a computer program in 1D
and CFD simulations in steady state and unsteady state in 2D and 3D were established.
The core channel of the nuclear research reactor is studied and the temperatures of the
coolant, the cladding and the flesh are calculated. The transient simulations reproduce a
LOFA (Loss Of Flow Accident) in order to properly design the redundant pumps of the
protection system. The ONB (onset of nucleate boiling) is verified in all simulations.
Keywords : Thermohydraulics, Nuclear reactor, plate fuel, CFD, LOFA
Résumé
Afin d’assurer la sûreté des réacteurs nucléaires de recherche, la détermination de leurs
caractéristiques thermohydrauliques nécessite un suivi permanent de certaines grandeurs
physiques telles que la pression, la température et les flux de chaleur qui règnent
dans leurs cœur. Ainsi, pour avoir une évolution correcte de ces grandeurs physiques
un programme informatique en 1D et des simulations CFD en régime stationnaire et
transitoire en 2D et 3D ont été établi. Le canal du cœur du réacteur nucléaire de recherche
est étudié et les températures du caloporteur, de la gaine et de la chair sont calculées.
Les simulations transitoires reproduise un accident LOFA (Loss Of Flow Accident) dans
le but de bien dimensionner les pompe redondantes des système de protection. L’ONB
(onset of nucleate boiling) est vérifié dans toutes les simulations.
Mots clés : Thermohydraulique, réacteur nucléaire, combustible plaque, CFD, LOFA.
Remerciements
On remercie Dieu le tout-puissant de nous avoir donné le courage, la santé

et de nous avoir accordé son soutien durant les périodes les plus difficiles.
On remercie vivement notre cher encadreur, docteur Kamel SIDI ALI de

nous avoir aidé et orienté dans notre travail et surtout d’avoir cru en nos
compétences.
On remercie également notre chère professeur et directrice de laboratoire

de recherche, madame Benmamar Sadiaa, pour sa disponibilité ainsi que
ses judicieux conseils.
On tient à adresser nos sincères remerciements à Pr.Smaili Arezki et

monsieur Ouchene Samir pour le temps qu’ils ont consacré pour nous aider
durant notre projet de fin d’étude.
Nos vifs remerciements vont également aux membres du jury pour l’intérêt
qu’ils ont porté à notre travail et d’avoir bien voulu l’examiner et l’enrichir
de leurs propositions.
On remercie chaleureusement nos familles et tous nos amis pour leur

soutien et leurs encouragements.
On souhaite aussi remercier tous les enseignants de l’École Nationale

Polytechnique, et en particulier nos professeurs d’hydraulique qui nous ont
encadré auparavant et tous nos enseignants pour les connaissances qu’ils
nous ont transmises, leur disponibilité et leurs efforts.
Que tous ceux qui ont contribué de près ou de loin à la réalisation de ce

modeste travail trouvent ici l’expression de notre sincère gratitude.
Merci à toutes et à tous.

Dédicaces
À mes parents, merci d’avoir toujours cru en moi. Merci pour votre amour,
vos prières, vos encouragements, votre soutien et votre présence tout au
long de mes études.
Merci de m’avoir donné la chance de devenir celle que je suis aujourd’hui.
J’espère que vous êtes fiers de moi.
À mon frères Hamoud et ma grand-mère qui m’ont soutenu durant toute

ma vie et qui ont toujours su me remonter le morale quand il le fallait.
À ma binôme Ouarda, pour tous les bons moments passés ensemble, pour
cette expérience et son soutien tout au long de nos années d’études
À tout mes chers amis Anes, Nesrine, Assia, Asma, Lilya et Fazila qui
m’ont encouragés jusqu’à la fin
À tous ceux qui ont eu confiance en mes capacités pendant mon parcours.
À tous ceux qui me sont chers, je dédie ce modeste travail.
Yasmine
Dédicaces
A ma chère Maman
Aucune dédicace ne saurait exprimer mon amour, ma gratitude et ma

reconnaissance pour la force et le soutien que tu m’as apporté et les
sacrifices que tu as fait pour nous chaque jour pour qu’on ne manque de
rien et pour que je devienne la personne que je suis aujourd’hui. Que ce
modeste travail soit l’expression des vœux que tu n’as cessé de formuler
durant tes prières.
A mon cher Papa
A qui je dois ma vie, ma réussite et tout mon respect. Tu as toujours été là

pour m’encourager et me soutenir tout au long de mon parcours. Que ce
travail soit l’expression de mon amour et de ma profonde reconnaissance.
A mes chères sœurs et mon frère
Qui ont partagé avec moi tous les moments d’émotions lors de la réalisation
de ce travail. Aucun mot ne pourrait exprimer mon amour et ma profonde
reconnaissance pour tout ce que vous faites pour moi au quotidien.
A mes chères Tantes et ma cousine
Zou, Mimi et Lila que j’aime, je vous considère comme ma deuxième

famille, merci pour votre amour, votre soutien et pour vos conseils
judicieux que vous ne cessassiez de me donner.
A mon cher Papi
L’homme qui m’a fait grandir, qui as toujours été là et qui n’a jamais dit
non à mes demandes, merci pour tout.
A ma chère Binôme Yasmine
Avec qui j’ai partagé cette expérience enrichissante, merci pour ta patience
et pour ton engagement dans ce travail.
Enfin, je dédie ce travail à tous mes amis proches avec qui j’ai passé les
meilleurs moments durant mon parcours.
Ouarda
Table des matières
Table des figures 9

Liste des tableaux 1
Liste des symboles 4
Liste des abréviations 1

5
Introduction générale 19
1 Généralités sur les réacteurs nucléaires de recherche ;

technologie, sûreté nucléaire et mode de fonctionnements 21
1.1 Introduction.......................................................................................22
1.2 Composants d’un réacteur nucléaire de recherche de type MTR
(Material Testing Reactor)...............................................................22
1.2.1 Le coeur du réacteur nucléaire de recherche..........................24
1.2.2 Le système de refroidissement du coeur.................................28
1.3 Types d’accidents pouvant survenir dans un réacteur nucléaire
de recherche de type MTR...............................................................32
1.3.1 LOCA (Loss Of Coolant Accident) Accident de perte du
fluide de refroidissement.........................................................33
1.3.2 LOFA (Loss Of Flow Accident) accident de perte de débit 34
1.3.3 RIA (Reactivity Insertion Accident) accident d’insertion
de réactivité............................................................................35
1.4 Modes de fonctionnement d’un réacteur nucléaire de recherche
de type MTR................................................................................35
1.4.1 Fonctionnement normal..........................................................35
1.4.2 Fonctionnement accidentel.....................................................36
1.5 Conclusions.............................................................................................37
2 La thermohydraulique d’un canal du coeur d’un réacteur

nucléaire de recherche en monodimensionnel 38
2.1 Introduction.......................................................................................39

Table des
2.2 Caractéristiques physiques et de design des réacteurs de

recherche nucléaire de type MTR.....................................................39
2.3 Mise en équation thermohydraulique d’un canal du coeur d’un
réacteur nucléaire de recherche.........................................................41
2.3.1 Étude thermohydraulique asymptotique................................41
2.3.2 Analyse de sûreté nucléaire....................................................47
2.4 Évolution des quantités physiques dans trois réacteurs de
différentes puissances.......................................................................48
2.4.1 Résultats obtenus pour une application au réacteur
nucléaire de recherche de 2MW (AIEA)................................48
nucléaire de recherche de 10MW (AIEA)..............................56
nucléaire de recherche de 1MW (NUR).................................63
2.5 Conclusions.............................................................................................68
3 Modélisation thermohydraulique en 2D d’un canal du

coeur d’un réacteur nucléaire de recherche utilisant un
combustible à flux de chaleur uniforme et constant 70
3.1 Introduction.......................................................................................71
3.2 Configuration du canal thermohydraulique en 2D...........................71
3.3 Mise en équation thermohydraulique...............................................72
3.3.1 Les équations de conservation................................................72
3.3.2 Modélisation de la turbulence................................................73
3.4 Analyse thermohydraulique à l’état stationnaire............................76
3.4.1 Mise en équation.....................................................................76
3.4.2 Conditions aux limites............................................................77
3.5 Analyse thermohydraulique à l’état transitoire...............................78
3.5.1 Mise en équation.....................................................................78
3.5.2 Conditions aux limites............................................................79
3.6 Présentation du logiciel Comsol.......................................................80
3.7 Maillage du domaine numérique.......................................................80
3.8 Méthode des éléments finis...................................................................82
3.9 Résultats obtenus pour l’écoulement stationnaire............................83
3.10 Résultats obtenus pour l’écoulement transitoire..............................86
3.10.1 Évolution en fonction du temps.............................................86
3.10.2 Évolution temporelle de la température maximale................88
3.11 Conclusion..............................................................................................90
Mémoire de PFE TIHAL et 7

Table des
4 Analyse thermohydraulique, en 3D à l’état stationnaire et
transitoire, d’un canal du coeur d’un réacteur nucléaire de
recherche 91
4.1 Introduction.......................................................................................92
4.2 Présentation du Code FLUENT.......................................................92
4.3 Mise en équation thermohydraulique du canal en 3D.....................93
4.3.1 Les équation de conservation................................................93
4.3.2 Le modèle de turbulence........................................................95
4.4 Conditions aux limites......................................................................95
4.5 Maillage du domaine numérique à étudier.......................................96
4.6 Méthode des volumes finis....................................................................99
4.7 Analyse thermohydraulique à l’état stationnaire...........................100
4.8 Résultats et interprétations de la simulation stationnaire...........102
4.8.1 Pour une vitesse constante...................................................102
4.8.2 Analyse d’un LOFA constant en mode stationnaire pour
plusieurs pertes de débits......................................................104
4.9 Analyse thermohydraulique à l’état transitoire..............................107
4.10 Résultats et interprétations de la simulation transitoire..............108
4.11 Conclusion............................................................................................116
Conclusion générale 118
Bibliographie 121
Annexes 124
A Structure du programme Matlab pour le calcul des

températures en monodimensionnel 125
B Étapes de simulation sur Comsol en régime stationnaire et

transitoire 130
C Étapes de simulation FLUENT en régime stationnaire 139
D Étapes de simulation FLUENT en régime transitoire 147

Table des figures
1.1 Schéma d’une réaction de fission.........................................................23

1.2 Composants d’un réacteur nucléaire de recherche...........................24
1.3 Schéma simplifié du coeur d’un réacteur nucléaire de recherche
MTR type piscine.............................................................................25
1.4 Assemblage combustible de type plaque..........................................26
1.5 Les barres de contrôle d’un réacteur nucléaire de recherche MTR. 27
1.6 Coupe transversale du coeur d’un réacteur nucléaire de recherche. 28
1.7 Échangeur de chaleur de type plaque...............................................29
1.8 Échangeur de chaleur de type tube et calandre...............................30
1.9 Tour de refroidissement humide.......................................................31
1.10 Tour de refroidissement sèche..........................................................31
1.11 Tour de refroidissement hybride.......................................................32
1.12 Fonctionnement d’un réacteur nucléaire de recherche.....................36
2.1 Canal du coeur d’un réacteur nucléaire de recherche......................42
2.2 Encombrements d’un canal du coeur...............................................44
2.3 Distribution de la température en fonction de la longueur
adimensionnelle du canal pour un réacteur nucléaire de
recherche 2MW.................................................................................49
2.4 Comparaison de la distribution des températures avec des
travaux publiés en fonction de la longueur adimensionnelle du
canal pour un réacteur nucléaire de recherche 2MW.......................50
adimensionnelle du canal..................................................................52
2.6 Comparaison des profils des températures des deux modèles
en fonction de la longueur adimensionnelle du canal pour un
réacteur nucléaire de recherche de 2MW.........................................53
2.7 Flux de chaleur critique ONB en fonction de la longueur
recherche 2MW.................................................................................55
2.8 Rapport du flux de chaleur critique ONB le long du canal pour
un réacteur nucléaire de recherche 2MW.........................................56

Table des
2.9 Distribution des températures en fonction de la longueur

recherche 10MW...............................................................................57
2.10 Comparaison des distributions de température en fonction de la
longueur adimensionnelle du canal pour un réacteur nucléaire
de recherche 10MW...............................................................................58
recherche 10MW...............................................................................59
2.12 Comparaison de la distribution des températures des deux
modèles en fonction de la longueur adimensionnelle du canal
pour un réacteur nucléaire de recherche 10MW...............................60
recherche de 10MW..........................................................................62
un réacteur nucléaire de recherche de 10MW..................................63
recherche de 1MW............................................................................64
recherche 1MW.................................................................................65
2.17 Comparaison de la distribution des températures des deux
modèles en fonction de la longueur adimensionnelle du canal
pour un réacteur nucléaire de recherche de 1MW............................66
recherche 1MW.................................................................................67
un réacteur nucléaire de recherche de 1MW....................................68
3.1 Géométrie du canal en 2D................................................................72
3.2 interface Comsol multiphysics..........................................................80
3.3 Type d’éléments de maillage............................................................81
3.5 discrétisation d’un domaine par éléments finis.................................83
3.6 Distribution de la température du fluide de refroidissement le
long de l’axe y........................................................................................84

Table des
3.7 Distribution de la température de la gaine le long de l’axe y. . 85

3.8 Distribution de la température de la chair le long de l’axe y. . . 85
3.9 Évolution de la température le long du canal à t=1s....................86
3.12 Évolution de la température à y=0,6 m en fonction du temps.......89
4.1 schéma du canal thermohydraulique en 3D.....................................93
4.3 Géométrie du canal 3D.....................................................................97
4.4 Maillage du canal selon les trois plans de l’espace............................98
4.5 Température en fonction du nombre de mailles...............................99
4.6 Cellule discrétisée à deux dimensions.............................................100
4.7 Fonction UDF du flux de chaleur...................................................102
4.8 Contour de la température du caloporteur et de la gaine sur le
plan xz en stationnaire....................................................................103
4.9 Profils des températures en fonction de la longueur
adimensionnelle..............................................................................103
4.10 Contour des températures de la gaine et du caloporteur sur le
plan xz pour 5%, 10 %, 15% et 20%..............................................104
4.11 Profil de la température pour plusieurs débits en fonction de la
longueur adimensionnelle...............................................................105
4.12 Évolution au point le plus chaud de la gaine de la température
en fonction du pourcentage de perte de débit................................106
4.13 Fonction UDF de la variation de vitesse........................................108
4.14 Contour des températures du caloporteur sur le plan xz à
chaque pas de temps de 6s à 10s....................................................109
4.15 Profil de la température en fonction du temps..............................109
4.16 Contour des températures de la gaine sur le plan xz à chaque
pas de temps de 6s à 10s................................................................111
4.17 Profil de température en fonction du temps...................................112
4.19 Contour des températures du caloporteur sur le plan xz à
chaque pas de temps 6s à 15s.........................................................114
4.20 Contour des températures de la gaine sur le plan xz de contact
avec le caloporteur à chaque pas de temps 6s à 15s.......................115
A.1 Structure du code Matlab . . . . . . . . . . . . . . . . . . . . 127
A.2 Code Matlab de la ligne 1 à 31 . . . . . . . . . . . . . . . . . . 128

Table des
A.3 Code Matlab de la ligne 32 à 62.....................................................128

A.4 Code Matlab de la ligne 63 à 93.....................................................129
A.5 Code Matlab de la ligne 92 à 120...................................................129
B.1 Création de la géométrie sur COMSOL.........................................131
B.2 Fenêtre Géométrie (COMSOL)......................................................131
B.3 Choix des matériaux des composant du canal (COMSOL)...........132
B.4 Choix du matériau pour le caloporteur..........................................132
B.5 Choix du matériau pour la gaine....................................................133
B.6 Choix du matériau pour la chair....................................................133
B.7 Choix du modèle de turbulence......................................................134
B.8 Fenêtre de l’écoulement turbulent..................................................134
B.9 Indication de la paroi......................................................................134
B.10Vitesse à l’entrée du canal..............................................................135
B.11Vecteur force de gravité..................................................................135
B.12Indication de la pression à la sortie du canal.................................135
B.13Axe de symétrie sur la géométrie...................................................136
B.14Fenêtre de la physique : transfert de chaleur entre solide et fluide.136
B.15Température d’entrée du caloporteur.............................................136
B.16Script de calcul de la moyenne à l’état stationnaire (ligne 1-31) 137
B.17Script de calcul de la moyenne à l’état stationnaire (ligne 31-59).137
B.18Script de calcul de la moyenne à l’état transitoire à 1s.................138
C.1 Tableau Fluent sur Workbench.......................................................140
C.2 Configuration générale........................................................................140
C.3 Fenêtre "models".............................................................................141
C.4 Matériaux........................................................................................141
C.5 Propriétés de la chair......................................................................141
C.6 Propriétés de la gaine......................................................................142
C.7 Propriétés de l’eau..........................................................................142
C.8 Fenêtre cell zone..................................................................................142
C.9 Conditions aux limites....................................................................143
C.10Vitesse à l’entrée du canal..............................................................143
C.11Pression à la sortie du canal...........................................................144
C.12Source de chaleur.................................................................................144
C.13Fonction UDF de la source de chaleur...........................................144
C.14Methods.................................................................................................145
C.15Residual.................................................................................................145

Table des
C.16Initialisation..........................................................................................146
C.17Calculs et itération.........................................................................146
D.1 Fonction UDF de la vitesse............................................................148
D.2 Calculs et itérations........................................................................149
D.3 Script du code matlab de calcul de la moyenne en stationnaire. 149

Liste des tableaux
1.1 Répartition de l’énergie provenant d’une réaction de fission . . 23

2.1 Données principales d’un réacteur de 1 MW...................................40
2.2 Données principales d’un réacteur 2 MW........................................40
2.3 Données principales d’un réacteur 10 MW......................................41
2.4 Les corrélations pour le calcul du nombre de Nusselt.....................45
2.5 Flux de chaleur pour un réacteur 2 MW..........................................54
2.6 Flux de chaleur d’un réacteur 10 MW.............................................61
4.1 Paramètre pour le calcul de l’épaisseur de la couche limite..............98

LIsT DEs sYMBoLEs
E
Am Section transversale de la chair (m2)

Az Surface de la section d’écoulement (m2)
Cp Chaleur spécifique du caloporteur(J/kg.C°)
Dh Diamètre hydraulique (m)
Dω Diffusion croisée ( −)
E Énergie totale J
ec Épaisseur de la gaine (m)
ech Épaisseur du canal (m)
em Épaisseur de la chair (m)
F Facteur de puissance ( −)
F1, F2 Fonction de mélange ( −)
Fax Facteur de puissance axiale ( −)
Fd Facteur multiplicatif ( −)
F NUC Facteur de puissance nucléaire (−)
Frad Facteur de puissance radiale ( −)
FP T N Nombre total de plaques combustibles (−)
G Distance résiduelle à la paroi (m)
GK Générations de l’énergie cinétique turbulente (J)
Gm Vitesse massique (kg/s)
h Coefficient de transfert de chaleur par convection (W/m2.C)
hf Enthalpie de l’écoulement (J)
hj Coefficient de transfert de chaleur par convection pour une espèce j
K Énergie cinétique turbulente
k Conductivité thermique (W/m2.C)
kc Conductivité thermique de la gaine (W/m2.C)
km Conductivité thermique de la chair (W/m2.C)
Lch Longueur totale du canal (m)

Le Longueur extrapolée (m)

Liste des
Lh Longueur active de la plaque (m)

lch Largeur du canal (m)
lm Largeur de la chair (m)
mch Débit massique du fluide (kg/s)
Nu Nombre de Nusselt ( −)
n Vecteur normal extérieur ( −)
Ph Périmètre chauffé (m)
Pr Nombre de Prandlt ( −)
PU Puissance du réacteur nucléaire (W )
p Pression (bar)
pcr pression critique (bar)
Qp Source de chaleur ponctuelle (J)
Qsv Débit standard (m3/s)
Qted Amortissement thermo-élastique (s−1)
Qvd Dissipation visqueuse
q′ Flux de chaleur (W/m2)
q′ Densité de puissance axiale (W/m3)
3
q′ a Densité de puissance moyenne (W/m )
3
qc Densité de puissance maximale (W/m )
q0 Flux de chaleur entrant (W/m2)
qav Flux de chaleur moyen (W/m2)
qcr Flux de chaleur critique (W/m2)
Rce Rayon à l’entrée du canal (m)
Rcs Rayon à la sortie du canal (m)
S Terme source (J)
Sφ Source par unité de volume (J/m3)
SV R Rapport volume/surface de la chair (m)
T Température (°C)
Tc Distribution de la température axiale de la gaine (°C)
Tf Distribution de la température axiale du fluide de refroidissement
Tf1 Température d’entrée du caloporteur (°C)
Tm Distribution de la température axiale de la chair (°C)
Tsat Température de saturation (°C)
t Temps (s)
U Vitesse du fluide (m/s)
U ∗ Vitesse de friction locale
u Vitesse selon x (m/s)
ui Composante de vitesse dans la direction correspondante (m/s)
u′i Vitesse fluctuante (m/s)

Liste des
u¯i Vitesse moyenne (m/s)

Vch Flux de chaleur moyen (W/m2)
Vm Volume de la chair (m3)
v Vitesse selon y (m/s)
WT Débit total (m3/s)
w Vitesse selon z (m/s)
y Direction axiale
y+ position sans dimension de y
A→ Vecteur de surface (m2 )

→g Vecteur gravité (m/s2)
J→j Vitesse d’une espèce j (m/s)
U→ Vecteur vitesse
Γ frontières du domaine
Γe frontière d’un élément
Γφ Coefficient de diffusion
δij Tenseur d’unité ( −)
ε Taux de dissipation d’énergie turbulente (m2/s3)
λ Conductivité thermique du fluide (W/m.°C)
µ Viscosité dynamique du fluide (kg/m.s)
µt Viscosité turbulente (kg/m.s)
ν Viscosité cinématique du fluide (m2/s)
νt Viscosité statique turbulente (m2/s)
Ωij Tenseur du taux de rotation moyen
Ωe élément
ω Taux de dissipation spécifique (m2/s3)
ρ Masse volumique du fluide (kg/m3)
ρst Masse volumique standard (kg/m3)
τ¯ Tenseur des contraintes visqueuses (kg/m.s2)
∇φ Gradient de φ

LIsT DEs abRÉvIaTIons
E
AIEA Agence Internationale de l’Energie Atomique

CAO Conception Assistée par Ordinateur
CATIA Conception Assistée Tridimensionnelle Interactive Appliquée
CFD Computationnal Fluid Dynamics
DNB Departure of Nucleate Boilling
DNS Direct Numerical Simulation
I −DEAS Integrated Design Engineering and Analysis Software
LES Large Eddy Simulation
LOCA Loss Of Coolant Accident
LOFA Loss Of Flow Accident
MTR Material Testing Reactor
NUR Nuclear Uranium Reactor
OFI Onset of Flux Instability
ONB Onset Nucleate Boiling
RANS Reynolds Averaged Navier
Stokes RIA Reactivity Insertion
Accident
RNG Re-Normalization Group
SST Shear Stress Transport
UDF User Define Function
InTRodUcTIon généRaLE
Dans le domaine nucléaire, Les réacteurs nucléaires de recherche jouent

un rôle clé dans son évolution et son développement. Ils permettent des
expériences nucléaires dans diverses disciplines on cite par exemple la
neutrographie ou encore les effets des transmutations tels que le dopage du
silicium.
Les réacteurs nucléaires de recherche sont des installations de puissance

modéré qui nécessitent moins de combustible comparé aux réacteurs
nucléaires de production d’énergie électrique et ont la particularité de
fonctionner à des températures plus basses. Ils comprennent une multitude
de design dont les plus communs sont les réacteurs nucléaires de recherche à
coeur ouvert et fermé.
Les réacteurs nucléaires de recherche à coeur ouvert utilisent les

combustibles nucléaire de type plaque et sont refroidit dans le but de faire
baisser la température qui se dégage du coeur car un quelconque
dépassement de celle ci peut provoquer des dégâts colossaux des
composants du coeur comme par exemple la gaine du combustible. A cet
effet, une étude thermohydraulique lors d’un fonctionnement normal est
primordiale pour surveiller la distribution des températures du fluide
caloporteur, de la gaine et de la chair et faire en sorte que les critères de
sûreté soient respectés. Il existe plusieurs manières pour procéder à cette
étude. Le cas le plus simple serait de faire une analyse monodimensionnelle
mais il est possible de faire l’analyse bidimensionnelle et tridimensionnelle
en utilisant les codes CFD.
Notre projet de fin d’étude est organisé en quatre chapitres présentées

comme suit :

Introduction générale
Le premier chapitre est consacré à la description des réacteurs nucléaires de

recherche et à la présentation de la problématiques et des notions de base
nécessaire dans notre travail.
Dans le chapitre II, une analyse thermohydraulique stationnaire est effectuée

en monodimensionnel et appliquée sur différents réacteurs nucléaires
de recherche de puissance distinct et cela en utilisant deux modèles
mathématiques valables. On a par la suite confirmer la sûreté du réacteur
en vérifiant certains critères de sûreté comme le flux de chaleur critique et
moyen ONB.
Dans le chapitre III on a reproduit l’évaluation expérimentale dans un

réacteur nucléaire de recherche avant insertion du combustible (Uranium)
à l’aide du logiciel COMSOL. Deux simulations sont réalisées à cet effet,
la première en mode stationnaire et la seconde en mode transitoire. Le
but recherché est de connaître la température que peut atteindre le fluide
caloporteur lors d’un fonctionnement normal ou accidentel.
On a présenté dans le dernier chapitre à une simulation tridimensionnelle de

l’écoulement du fluide caloporteur dans un canal du réacteur nucléaire de
recherche à l’aide du code de calcul Fluent. De même que pour le chapitre
III, on commence par une simulation stationnaire puis une simulation
transitoire qui reproduit un accident LOFA (Loss Of Flow Accident).

CHAPITRE 1
GénéRaLITÉs sUR LEs RÉacTEURs

nUcLÉaIREs de REcheRche ;
TEchnoLogIE, sûRETÉ nUcLÉaIRE ET
mode de foncTIonnemenTs

Chapitre 1. Généralités sur les réacteurs nucléaires de recherche ;
technologie, sûreté nucléaire et mode de
1.1 Introduction
Les réacteurs nucléaires de recherche sont des installations de puissance
modérées. Leurs puissances peuvent varier de quelques KW à des dizaines
de MW. Ils permettent des expériences nucléaires dans diverses disciplines.
On citera l’exemple de la neutrographie qui permet de prendre des clichés
d’équipements plus précis que ceux des rayons x, ou encore les effets des
transmubtions tels que le dopage en silicium.
Dans ce chapitre, on va présenter la conception classique d’un réacteur

nucléaire de recherche de type MTR utilisant un combustible à plaque. Puis
on donnera les accidents qui peuvent se produire dans ce type de réacteurs.
En dernier, deux mises en équation thermohydrauliques d’un canal d’un
réacteur nucléaire de recherche sont présentées. La première concerne l’état
stationnaire et la seconde l’état transitoire.
1.2 Composants d’un réacteur nucléaire de recherche

de type MTR (Material Testing Reactor)
Il existe plusieurs types de réacteurs nucléaire de recherche. Ces derniers
diffèrent de part leurs design et leur conception. Dans ce travail, l’intérêt
sera porté sur un réacteur nucléaire de recherche type piscine utilisant un
combustible à plaques.
Un réacteur nucléaire de recherche type piscine fonctionne à pression

atmosphérique. Il permet un accès facile aux emplacements d’irradiations
pour l’expérimentation dans le coeur ou à l’extérieur du coeur car il utilise
des faisceaux intenses de neutrons. Il donne aussi une vision directe sur le
fonctionnement du réacteur et la réaction qui a lieu.
Dans ce réacteur, on attribue à l’eau légère ( H2O) le rôle de modérateur

et de fluide de refroidissement. L’eau est aussi utilisée comme moyen de
protection efficace contre le rayonnement généré par la fission. L’eau permet
donc de maintenir des niveaux acceptables de radioactivité autour de la
piscine. Une réaction de fission (*, 2008) est définie par la scission d’un
noyau atomique lourd dit "fissile" en deux fragments lorsqu’il y a collision
avec un neutron. Il en résulte un dégagement d’une quantité importante
d’énergie accompagnée par l’émission de neutrons et de rayonnement
ionisant comme


montré sur la figure (1.1).
FIGURE 1.1 – Schéma d’une réaction de fission.
L’énergie libérée qui provient de la réaction de fission d’un atome est

d’environ 200 MeV. Elle est répartie sur les produits de fission. On donne
sur le tableau (1.1) un exemple de répartition de l’énergie entre les produits
de fission de l’atome U235 (Stacey, 2018).
Tableau 1.1 – Répartition de l’énergie provenant d’une réaction de fission
Produits de fission de U235 Énergie libérée (MeV)

Énergie cinétique des fragments de fission 165
Énergie cinétique des photons gamma 7
Énergie cinétique des neutrons 5
Radioactivité gamma et bêta 11
Énergie (invisible) des neutrinos 13
Total 200
Pour ce qui est de la conception d’un réacteur nucléaire de recherche de

type MTR, on décrit sur la figure (1.2) les principaux composants dont est
constitué ce réacteur.

FIGURE 1.2 – Composants d’un réacteur nucléaire de recherche.
1 : Coeur du réacteur nucléaire.

2 : Pompe primaire.
3 : Échangeurs de chaleur.
4 : Pompe secondaire.
5 : Tours de refroidissement.
I : Circuit primaire.
II : Circuit secondaire.
Le composant (1) représente le coeur du réacteur nucléaire qui est déposé

dans une cuve remplie d’eau légère. Cette eau sert à refroidir les éléments
combustibles et à la modération des neutrons. Cette dernière circule dans le
circuit primaire (I) où l’écoulement est assuré par une pompe primaire (2).
Les deux circuits ( primaire (I) et secondaire (II)) sont reliés entre eux par
des échangeurs de chaleur (3). La circulation de l’eau légère du circuit (II)
se fait entre les échangeurs de chaleur et les tours de refroidissement (5).
L’écoulement dans le circuit (II) se fait à l’aide de pompes secondaires (4)
(Sidi Ali, 2002).
1.2.1 Le coeur du réacteur nucléaire de recherche
Le coeur du réacteur nucléaire de recherche est positionné au fond de

la cuve comme le montre la figure (1.3). Il est constitué d’un arrangement
d’éléments combustibles insérés dans une plaque matricielle fabriquée en
aluminium.

FIGURE 1.3 – Schéma simplifié du coeur d’un réacteur nucléaire de recherche MTR type
piscine.
Sur la figure (1.3), (A) représente l’eau légère, (B) le coeur (combustible
plaque), (C) la cuve du réacteur, (b1) la chair du combustible et (b2) la
gaine du combustible.
En plus des éléments combustibles, le coeur comporte des éléments de

contrôle et de sécurité. La plaque matricielle permet la circulation de l’eau
de refroidissement entre les surfaces extérieures des éléments combustibles.
En règle générale, le cœur du réacteur est protégé par une épaisseur
d’eau d’environ 8 m de hauteur et, latéralement, par 1,2 m d’eau plus
2,4 m de béton baryté ( béton lourd qui offre une barrière efficace contre
les rayons Gamma et les autres émanations radioactives) (Michel et al.,
2015).
En plus du fluide caloporteur et modérateur, quatre des principaux

composants du coeur d’un réacteur nucléaire de recherche type MTR sont
présentés ci-après :
1.2.1.1 L’ élément combustible
Il est formé de plusieurs plaques combustibles incurvées ou plates. Dans

une plaque combustible, il y a ce qu’on appelle une gaine en aluminium qui
sert à confiner la chair qui contient la matière fissile. Cette configuration

de "type plaque" permet d’isoler la matière fissile de l’environnement du
réacteur. De plus, Les combustibles plaques sont disposés en parallèle pour
former un assemblage combustible comme le montre la figure (1.4) (*,
2008).
FIGURE 1.4 – Assemblage combustible de type plaque.
1.2.1.2 Les éléments de contrôle
Il s’agit de barres de contrôle comme montré sur la figure (1.5) posées au

dessus du combustible contenant des matériaux absorbants dont le rôle est
de maintenir un équilibre entre les neutrons qui naissent ou les neutrons qui
meurent, c’est ce qu’on appelle la criticité. leurs rôle est aussi d’arrêter ou de
ralentir la réaction de fission par absorption de neutrons qui se trouvent dans
le réacteur. Ces barres ont trois fonctions principales : la compensation de
la réactivité pendant un cycle de fonctionnement, le contrôle du démarrage
et de l’arrêt pendant le fonctionnement normal, et l’arrêt rapide dans des
conditions anormales.

FIGURE 1.5 – Les barres de contrôle d’un réacteur nucléaire de recherche MTR.
1.2.1.3 Un réflecteur
Un réflecteur sert à minimiser la dispersion et l’échappement des

neutrons du coeur du réacteur et garder le régime de la réaction de fission.
1.2.1.4 Les canaux expérimentaux
Ce sont des canaux horizontaux qui traversent la cuve du réacteur

nucléaire de recherche pour atteindre le coeur (figure 1.6). Il servent à
effectuer des expériences au coeur du réacteur nucléaire.

FIGURE 1.6 – Coupe transversale du coeur d’un réacteur nucléaire de recherche.
1.2.2 Le système de refroidissement du coeur
La réaction de fission dans le coeur du réacteur est une réaction

exothermique, c’est à dire qu’elle produit de la chaleur. Pour faire diminuer
cette chaleur on utilise un système de refroidissement.
Le système de refroidissement du réacteur est réparti en deux circuits,

on a un circuit primaire et un circuit secondaire. Ces deux circuits sont reliés
entre eux à l’aide d’un échangeur de chaleur. L’objectif de ces circuits est
de faire évacuer la chaleur produite par réaction de fission dans le coeur du
réacteur (Michel et al., 2015).
a) Le circuit primaire
Le circuit primaire est un circuit fermé qui passe par le coeur. Il

assure le refroidissement du combustible en faisant passer l’eau à
travers les plaques. Ainsi, chaque deux plaques adjacentes forment ce
qu’on appelle un canal. La circulation de l’eau légère se fait grâce aux
pompes hydrauliques qui aspirent l’eau du coeur et la font passer par
un échangeur de chaleur et retourner au coeur. L’écoulement peut se
faire de façon ascendante ou descendante selon le design du réacteur
(Michel et al., 2015).

b) Le circuit secondaire
Le circuit secondaire récupère la chaleur transmise du circuit primaire

par l’intermédiaire d’un échangeur de chaleur puis l’eau du circuit
secondaire passe par des tours de refroidissement afin de dissiper
l’énergie thermique dans l’atmosphère (Michel et al., 2015).
Le système de refroidissement (circuit primaire et circuit secondaire )
comprend des équipements de refroidissement :
a) Échangeur de chaleur
Il assure le transfert de l’énergie thermique sous forme de chaleur

du circuit primaire au circuit secondaire sans que les eaux des deux
circuits ne se mélangent. Il fait partie des équipements de conversion
d’énergie. Il existe plusieurs types, on cite comme exemples les
échangeurs de chaleur type plaque et type tube (Lewzisjr et al., 1964).
a.1) Échangeurs à Plaques
Il est formé de plusieurs plaques comme cela est montré sur la

figure (1.7). Il est très compact ce qui le rend très efficace. Le transfert
de chaleur se fait à travers le métal des plaque.
FIGURE 1.7 – Échangeur de chaleur de type plaque.
a.2) Échangeurs à tubes et calandre

Il convient aux applications à haute pression comme cela est montré

sur la figure (1.8). Ce type d’échangeurs peut utiliser plusieurs passes
pour un bon transfert de chaleur.
FIGURE 1.8 – Échangeur de chaleur de type tube et calandre.
b) Les tours de refroidissement
Une tour de refroidissement est généralement nécessaire pour évacuer

la chaleur provenant des échangeurs de chaleur figure (1.2). Cette
chaleur est transmise à l’environnement via des tours de
refroidissement (Sidi-Ali, 2018). Il existe différents type de tours de
refroidissement. On cite trois d’entre elles :
b.1) Tour de refroidissement humide (figure 1.9)
Le contact est direct entre l’eau et l’air et le transfert de chaleur

est augmenté par un transfert de masse ce qui le rend très performant.

FIGURE 1.9 – Tour de refroidissement humide.
b.2) Tour de refroidissement sèche (figure 1.10)
Dans une tour de refroidissement sèche, il y a des tubes qui assurent le

transfert de chaleur entre l’eau à l’intérieur des tubes et l’air.
FIGURE 1.10 – Tour de refroidissement sèche.
b.3) Tour de refroidissement hybride (figure 1.11)
Pour les tours de refroidissement hybrides, on prend le même

écoulement que dans les tours de refroidissement sèches, auquel
on ajoute un système de pulvérisation d’eau comme celui utilisé dans
les tours de refroidissement humides. Le refroidissement est encore
meilleur que pour les deux autres types de tours de refroidissement.

FIGURE 1.11 – Tour de refroidissement hybride.
c) Pompes hydrauliques
Les pompes hydrauliques sont indispensables pour l’écoulement de

l’eau dans le système de refroidissement. On en trouve dans les deux
circuits primaire et secondaire. Elles peuvent être disposées soit en
parallèle pour augmenter le débit soit en série pour augmenter la
pression.
1.3 Types d’accidents pouvant survenir dans un

réacteur nucléaire de recherche de type MTR
La sûreté occupe une place très importante dans le domaine nucléaire.
En effet peu importe le degré de protection et les dispositions mises en
place pour assurer un fonctionnement normal d’un réacteur nucléaire. Il y a
toujours un risque qu’un accident se produise.
Parmi les accidents les plus courants dans un réacteur de recherche,

on cite le LOCA (Loss Of Coolant Accident) qui désigne une perte en
quantité de fluide de refroidissement. Le LOFA (Loss Of Flow Accident) où
il s’agit d’une perte de débit d’écoulement. Le RIA (Reactivity Insertion
Accident) qui veut dire accident d’insertion de réactivité. Puis vient les
accidents de Coupure du courant électrique, de manipulation erronée
ou défaillance d’un équipement ou d’un composant, évènements internes
particuliers, évènements externes et erreurs humaines.
Les accidents qui concernent le coeur seulement sont les trois premiers

c’est à dire le LOCA, le LOFA et le RIA. Lorsque ces accidents ont lieu,
le fonctionnement du réacteur nucléaire doit rester dans les limites exigées
pour la sûreté nucléaire (GONNIER, 2008).
1.3.1 LOCA (Loss Of Coolant Accident) Accident de perte du

fluide de refroidissement
1.3.1.1 Causes
Cet accident se produit suite à :

— Un endommagement de la piscine
— La vidange de la piscine
— La défaillance des canaux expérimentaux.
On parle de perte de fluide lorsque la circulation naturelle interne est
supprimée. On obtient comme résultat prévisible, un échauffement excessif
du cœur ce qui peut potentiellement endommager le combustible si on
atteint une certaine puissance (deLimaReis et al., 2019).
1.3.1.2 Conséquences
On doit examiner toutes les conséquences afin de déterminer la

probabilité d’endommager le combustible. Même si cette probabilité s’avère
être très faible, il faut tenir compte du risque d’irradiation du personnel
d’exploitation à cause de la perte de la gaine. Par exemple un LOCA dont la
conséquence est une exposition partielle du cœur, doit être analysé pour
déterminer son potentiel d’endommagement du combustible.
Un autre exemple : Lorsque le niveau d’eau dans le cœur du réacteur

a baissé de manière significative, les plaques de combustible transférent
l’énergie dans un environnement gazeux seulement. La fusion du cœur dans
ce cas peut se produire, en commençant par les régions où la densité de
puissance est la plus élevée.
La simulation numérique de cet accident a pour but d’étudier le

comportement du coeur du réacteur nucléaire (gaine et chair du combustible)
et l’évolution de la température, de la pression et du niveau d’eau en
fonction du temps avant pendant et après l’accident c’est à dire lorsque le
réacteur revient à son état initial de fonctionnement normal.

1.3.2 LOFA (Loss Of Flow Accident) accident de perte de débit
La réduction du débit du fluide de refroidissement survient à la suite

de divers événements tels que la défaillance d’une pompe ou d’une vanne,
l’obstruction d’un canal ou la redistribution du débit. Ceci peut entraîner la
rupture des gaines par surchauffe (Qing et al., 2009).
Pour certains réacteurs nucléaire de recherche, ces accidents sont considérés
comme des accidents de dimensionnement, et certains facteurs de sécurité
technique sont intégrés pour faire face à leurs conséquences.
1.3.2.1 Causes
Comme pour un accident LOCA, le fonctionnement d’un réacteur à

pleine puissance lors d’un accident LOFA conduit à une réduction
significative de la capacité de transfert de chaleur. On cite les causes pour
lesquelles un accident LOFA peut se produire :
— Perte de puissance de pompage.
— Obstruction du canal d’entrée du cœur.
— perte de pression dans le système de refroidissement
1.3.2.2 Conséquences
La privation d’écoulement dans plusieurs canaux exigerait que la chaleur

soit transférée au caloporteur dans les canaux voisins non affectés. Si
la chaleur générée dans les canaux affectés ne peut pas être dissipée, le
combustible peut s’échauffer jusqu’à fondre. Lors de la fusion, la formation
de mousse dans le combustible (dont la quantité dépend du taux de
combustion) peut entraîner le contact des plaques entre elles, et donc créer
une propagation des dommages dans l’assemblage combustible.
Dans certains réacteurs nucléaires, il est indésirable que le liquide de

refroidissement s’écoule vers le bas dans le cœur. Par conséquent, tout
débris logé à l’entrée du cœur aura tendance à rester piégé. De plus, lors
d’accidents de perte de capacité de pompage, l’inversion de l’écoulement du
caloporteur, du mode de convection forcée au mode de convection naturelle,
entraîne des périodes de stagnation et éventuellement l’apparition de
dommages dans le cœur.

1.3.3 RIA (Reactivity Insertion Accident) accident d’insertion

de réactivité
L’accident d’insertion de réactivité RIA résulte, comme son nom

l’indique, d’une insertion de réactivité dans le coeur. À la différence du
LOCA, sa cinétique est très rapide, ce qui le rend très contraignant pour le
dimensionnement des réacteurs. Il concerne plus les réacteurs à eau
pressurisés REP (GONNIER, 2008).
1.4 Modes de fonctionnement d’un réacteur nucléaire

de recherche de type MTR
Dans le domaine nucléaire, le coeur d’un réacteur doit être soumis à un
refroidissement même après qu’une réaction de fission se soit terminer. On
décrit alors deux types de fonctionnement, normal et accidentel où chacun
d’entre eux a sa propre mise en équation et ses propres conditions aux
limites.
1.4.1 Fonctionnement normal
Le fonctionnement normal d’un réacteur nucléaire de recherche traduit

un état stationnaire pour lequel les variables qui le décrivent ne soient pas
fonction du temps. En d’autres termes les dérivés par rapport au temps dans
les équations de conservation de la masse, de la quantité de mouvement et
de l’énergie sont nuls.
Lors d’un fonctionnement normal l’ensemble des équipements fonctionne

normalement comme cela est schématisé dans la figure (1.12).

FIGURE 1.12 – Fonctionnement d’un réacteur nucléaire de recherche.
En thermohydraulique, il existe plusieurs façons pour procéder au calcul de

la température en monodimensionnel, bidimensionnel ou tridimensionnel. Il
suffit pour cela de poser les équations de base afin de les résoudre et obtenir
de bons résultats.
Dans le cas monodimensionnel, beaucoup de modèles mathématiques

ont été développés. On cite par exemple le modèle de El-Wakil (1971) ainsi
que le modèle de Todreas et al. (1989) dont les équations sont
asymptotiques. Pour le cas bidimensionnel et tridimensionnel, le calcul
se fait à travers des simulations CFD (Computationnal Fluid Dynamics) en
raison de la complexité de la géométrie et des équations.
1.4.2 Fonctionnement accidentel
Un fonctionnement accidentel est un fonctionnement pour lequel on note

un changement qui va rendre le coeur en état d’instabilité. Une situation
accidentelle signifie un fonctionnement à l’état instationnaire qu’on appelle
aussi transitoire.
Dans un réacteur nucléaire de recherche, un transitoire est mis en œuvre à

des fins expérimentale, il soumet le combustible à des conditions d’accident
sans affecter l’ensemble du réacteur, qui est protégé par son système de
sûreté contre les évolutions préjudiciables (**, 2017). La défaillance d’une
pompe, par exemple, mène à la réduction du débit d’eau de refroidissement

de 20%.
L’étude qui suivra devra répondre à deux comportements. Le premier

dans le cas du fonctionnement normal en vérifiant tous les paramètres
thermohydrauliques et le second dans le cas accidentel et vérifier les
problèmes de sûreté nucléaire qui lui sont liés. A savoir de début de
l’ébullition de l’eau et la fusion de la gaine du combustible. Le temps dans
ce dernier cas est un facteur déterminant pour l’activation des systèmes
redondants devant le SCRAM du réacteur et rester opérationnel en toute
sécurité.
On considère dans notre cas la défaillance de la pompe P1 sans que la
pression du système ne soit touchée, donc en considérant un LOFA de
quelques secondes.
En transitoire, la mise en équations est d’autant plus complexe. Pour cela,
les simulations se font qu’avec des codes de calculs de CFD.
1.5 Conclusions
Dans ce chapitre, la description d’un réacteur nucléaire de recherche
type MTR a été faite. Les composants des circuits primaire et secondaire
sont expliqués. Le processus de fission et l’ensemble de la technologie du
combustible sont décrits avec ample détails.
Une partie a été consacrée aux différents types d’accidents nucléaires

qui peuvent avoir lieu dans un réacteur nucléaire. Ceci a été fait dans le but
de bien présenter un des trois accidents qui va faire l’objet d’une partie de
ce travail.
Ce chapitre a été clos par la présentation des deux types de fonctionnement

possible dans un réacteur nucléaire de recherche qui sont normal et
accidentel.
Il est très important de noter que le coeur d’un réacteur nucléaire doit
être refroidi quelque soit la situation, soit en fonctionnement normal ou
accidentel.

CHAPITRE 2
La THERMohYDRaULIqUE D’Un canaL

dU coeUR D’Un RÉacTEUR nUcLÉaIRE
de REcheRche en monodIMEnsIonneL

Chapitre 2. La thermohydraulique d’un canal du coeur d’un
réacteur
2.1 Introduction
Les combustibles nucléaire de type plaque sont utilisés dans les réacteurs
nucléaires de recherche à coeur ouvert fonctionnant à de faibles puissances
et sont refroidis dans le but de diminuer et réguler la température du coeur
car un quelconque dépassement de celle ci provoque une dégradation de
la gaine du combustible. A cet effet, une étude thermohydraulique à l’état
stationnaire est primordiale pour surveiller le distribution des températures
du fluide de refroidissement, de la gaine et la chair et faire en sorte que les
critères de sûreté soient respectés.
Dans ce travail, on va effectuer une analyse thermohydraulique en régime

stationnaire sur différents réacteurs nucléaires de recherche de puissance
distinct et cela en usant de deux modèles mathématiques valables en
monodimensionnel. On essayera par la suite de vérifier que les critères de
sûreté nucléaire sont assurés.
2.2 Caractéristiques physiques et de design des

réacteurs de recherche nucléaire de type MTR
Chaque réacteur nucléaire de recherche a un design et des caractéristiques
qui lui sont propres. Le réacteur de recherche qu’on va étudier pour faire
une simulation est un réacteur de recherche MTR à coeur ouvert. Pour avoir
une meilleure analyse, on choisit de prendre trois réacteurs de ce type avec
des puissances différentes.
Les données relatives au réacteur nucléaire de recherche 1MW sont

présentées sur le tableau (2.1) tiré de l’article (Meftah et al., 2006).

réacteur
Tableau 2.1 – Données principales d’un réacteur de 1 MW.
Propriétés Valeurs
Nombre de plaque de combustible (SFE) 19
Nombre d’éléments combustible(SFE) 12
Épaisseur de la chair 0,070 cm
Épaisseur de la gaine 0,040 cm
Épaisseur du canal 0,270 cm
Largeur de la chair (SFE) 6,0 cm
Largeur du canal 6,6 cm
Hauteur de la chair (ou hauteur active) 61,5 cm
Conductivité thermique de la chair 0,15 W/cm.C
Conductivité thermique de la gaine 1,8 W/cm.C
Facteur de puissance axial 1,31
Facteur de puissance radial 2,03
Facteur de puissance total 2,9
Nombre de plaque de contrôle (SCE) 14
Nombre d’éléments de contrôle 5
Dans ce qui suit (tableau 2.2), sont présentées les principales données
pour un réacteur de 2 MW (Labani et al., 2014).
Tableau 2.2 – Données principales d’un réacteur 2 MW.
Nombre des éléments combustibles standard 19
Nombre des éléments combustibles de contrôle 4
Largeur du canal lch 6,64 cm
Largeur de la chair lm 6,3 cm
Longueur active de la plaque Lh 60 cm
Température d’entrée Tf1 38 C◦
Pression de sortie P 1,961 bar
Conductivité thermique de la gaine kc 180 W/cm.K
Conductivité thermique de la chair km 53,6 W/cm.K
Nombre de plaque des éléments Combustibles standard 19
Nombre de plaque des éléments Combustible de contrôle 15
Longueur totale de la plaque (canal) LT 62,5 cm
Épaisseur du canal ech 0,2916 cm
Épaisseur de la gaine ec 0,0381 cm
Épaisseur de la chair em 0,051 cm
Facteur de puissance radial Frad 2
Facteur axial de puissance Fax 1,58
Taux du débit total wT 300 m3/h
Vitesse d’écoulement dans le canal U 0,94 m/s
Longueur extrapolée Le 75,5 cm
Le tableau (2.3) ci-dessous donne les principales données pour un réacteur

de 10 MW (Labani et al., 2014)

réacteur
Tableau 2.3 – Données principales d’un réacteur 10 MW.
Nombre des éléments combustibles standard 23
Nombre des éléments combustibles de contrôle 5
Largeur du canal 6,65 cm
Largeur du combustible 5,1
Longueur du chauffage 60 cm
Température d’entrée 38 C◦
Pression de sortie
Conductivité thermique de la gaine 1,80 W/cm.K
Conductivité thermique de la chair 0,5 W/cm.K
Nombre de plaque des éléments Combustibles standard 23
Nombre de plaque des éléments Combustible de contrôle 17
Longueur totale de la plaque (canal) 62,5 cm
Épaisseur du canal 0,219 cm
Épaisseur de la gaine 0,038 cm
Épaisseur de la chair 0,051 cm
Facteur de puissance radial 1,4
Facteur axial de puissance 1,5
Taux du débit total 1000 m3/h
Vitesse d’écoulement dans le canal 2,97 m/s
2.3 Mise en équation thermohydraulique d’un canal du

coeur d’un réacteur nucléaire de recherche
La mise en équation concerne l’analyse thermohydraulique et l’analyse
de sûreté dans un canal chaud du coeur du réacteur nucléaire de recherche.
L’analyse thermohydraulique est défini comme étant une étude de transfert
de chaleur couplé à l’hydraulique des écoulements. Elle s’intéresse à l’état
du fluide, à son mode de circulation et aux échanges thermiques dans le
coeur du réacteur nucléaire.
2.3.1 Étude thermohydraulique asymptotique
On va utiliser deux modèles thermohydrauliques. Le premier modèle est

celui d’El-Wakil (1971) et le second celui de Todreas et al. (1989).
Le calcul se fait par rapport à l’axe z comme sur la figure (2.1) ci-dessous :

réacteur
FIGURE 2.1 – Canal du coeur d’un réacteur nucléaire de recherche.
2.3.1.1 Modèle d’ El-Wakil (1971)
2.3.1.1.1 Distribution de puissance dans le canal du coeur
Les équations de El-Wakil (1971) sont des équation asymptotiques et

donnent directement la distribution de la température en fonction de la
longueur axiale du canal. Avant d’introduire ces équations, on donne la
forme de la distribution de puissance axiale dans le canal du coeur. Elle est
de forme sinusoïdale et est donnée par l’équation (2.1).
πz
q′(z) = q′ cos (2.1)
c
Le
Où, q (z) est la densité de puissance axiale, Le est la longueur extrapolée, qc′
′
est la densité de puissance maximale dans le canal à la moitié de la longueur

de la plaque combustible.
La densité de puissance maximale s’écrit comme suit équation (2.2) :
′ ′
qc = F q a (2.2)
Où, qa′ est la densité de puissance moyenne et F est le facteur de puissance.
Dans un canal moyen, F et est égal au facteur de puissance axial Fax
seulement alors que dans un canal chaud F est égal au facteur de
puissance nucléaire F NUC qui est le produit du facteur de puissance
axial Fax et du facteur de puissance radial Frad (Khedr, 2008).
La densité de puissance moyenne est donnée par (2.3) :
′ PU
qa = (2.3)
Vm
Où, PU est la puissance du réacteur nucléaire et Vm le volume de la chair.

réacteur
La distribution axiale du flux de chaleur dans le canal en fonction de

la densité de puissance axiale est donnée par l’équation (2.4) :
q′(z)
q(z) = (2.4)
SV R
Où, SVR est le rapport surface/volume de la chair du combustible
2.3.1.1.2 Distribution de température du fluide caloporteur, de la gaine et

de la chair
Les distributions de température du fluide caloporteur, de la gaine et

de la chair sont calculées selon l’axe z à partir d équations de El-Wakil
(1971).
La distribution de la température du fluide caloporteur (Tf ) dans le

canal axial est donnée par :
′
cq
πz πLh
Tf (z) = Tf1 A + sin + (2.5)
πC
Lmpmche Le sin 2Le
+
La distribution de la température axiale de la gaine (Tc) dans le canal est

évaluée par :
′
Tc(z) = Tf (z) q cem πz

cos (2.6)
+ 2h Le
La distribution axiale de la température de la chair (Tm), dans le canal s’écrit :

′
Tm(z) = Tc(z) + q cos

πz e2 emec
m + (2.7)
c Le 8km 2kc
Où, Am est la section transversale de la chair (=emlm), lm la largeur de
la chair, em l’épaisseur de la chair, Cp la chaleur spécifique du fluide de
refroidissement, Lh la longueur active de la plaque, Le est la longueur
extrapolée, mch le débit massique du fluide, ec l’épaisseur de la gaine, z
la longueur axiale du canal qui est égale à zéro au centre du canal, Tf1
la température d’entrée du fluide de refroidissement, km la conductivité
thermique de la chair, kc la conductivité thermique de la gaine et h le
coefficient de transfert de chaleur par convection.
Les encombrements du canal sont donnés sur la figure (2.2) où (2.2a)

réacteur
représente la section transversale et (2.2b) la section le long du canal selon
l’axe z.

réacteur
(a) Section transversale du canal (b) Coupe en long du canal suivant z
FIGURE 2.2 – Encombrements d’un canal du coeur.
Le débit massique du fluide de refroidissement dans le canal étudié est

donné par l’équation (2.8) :
mch ρFdWT (2.8)
= FP T N
Où, FPTN est le nombre total de plaques combustibles, WT est le

débit total, Fd est un facteur multiplicatif inférieur à l’unité (Fd=0,9) et ρ
la masse volumique du fluide.
2.3.1.1.3 Calcul du coefficient de transfert de chaleur par convection dans le

canal
Afin de calculer la température de la gaine (Tc), il est impératif de faire

le calcul du coefficient de transfert de chaleur par convection h. Il est donné
par l’équation (2.9) :
Nuλ (2.9)
h = Dh
Où, Nu est le nombre de Nusselt, Dh le diamètre hydraulique et λ la

conductivité thermique du fluide .
Le nombre de Nusselt traduit la qualité de l’échange thermique : une

augmentation de ce nombre traduit une contribution importante de
l’écoulement sur l’échange de chaleur avec la paroi. Pour le calculer, il
existe plusieurs corrélations, parmi elles on cite la corrélation de Dittus-
Boelter donnée par Glasstone et al. (1967), la corrélation de Sieder-Tate
utilisée par El-Wakil (1971), la corrélation de Colburn citée par Glasstone

réacteur
et al. (1967)

réacteur
et la corrélation de Petukhov dans Holman (1989).
Ces corrélations sont utilisées dans le cas d’un écoulement turbulent :

cas stationnaire d’un réacteur nucléaire de recherche (Khedr, 2008). Le
tableau (2.4) donne les principales formules des corrélations pour le calcul
du nombre de Nusselt :
Tableau 2.4 – Les corrélations pour le calcul du nombre de Nusselt
Corrélation formule de Nu domaine d’application

LT
Dittus-Boelter 0, 023Re0,8Pr 0,4 0, 7 ≤ Pr ≤ 102 Re ≥ 104 ≥ 10
Dh
LT
Sieder-Tate 0, 023Re 0,8Pr 0,4 µw 0,14 0, 7 ≤ Pr ≤ 16, 7 Re ≥ 104 ≥ 10
µ Dh
LT
Colburn 0,8 0,33 0, 7 < Pr < 160 104 < Re < 1, 2.105 > 60
0, 023Re Pr Dh
LT
Petukhov f
R P µ
0,11 0, 5 ≤ Pr ≤ 2.103 3.103 < Re < 5.106 ≥ 10
8 e r b Dh
√f 2
1,07+127 8 (P 3
−1) µw
r
1
Où f =
(1,82log10Re−1,64)2
UDh
Où, Re est le nombre de Reynolds donné par Re = , Pr le nombre de
µ
νCp
Prandtl donné par Pr = , ech l’épaisseur du canal, lch la largeur du canal,
λ
U la vitesse du fluide, ν la viscosité cinématique du fluide et µ sa viscosité
dynamique.
2.3.1.2 Modèle mathématique de Todreas et al. (1989)
Dans ce modèle, la distribution axiale de la température est liée à la

distribution radiale de la température. En régime permanent, l’équation
d’énergie dans un canal avec une surface d’écoulement axiale s’écrit :
d q Ph Gm dp Gm|Gm|
G m
h f= + +f (2.10)
A ρ d 2Dh
Où Gm est la Vitesse massique, hf l’enthalpie de l’écoulement, Ph le
périmètre chauffé, Dh le diamètre hydraulique, ρ la masse volumique du
fluide et Az la surface de la section d’écoulement.
En négligeant les termes d’énergie dûs au gradient de pression et à la

réacteur
dissipation par frottement, l’équation (2.10) s’écrit :

d
GmAz hf = q Ph (2.11)
dz
Ou
mch d hf = q′(z) (2.12)
dz
′
Où, mch est le débit massique et q (z) la distribution de puissance axiale.
2.3.1.2.1 Température du caloporteur
En intégrant l’équation (2.12) sur la longueur axiale après avoir remplacé

′
le terme q (z) par l’équation (2.1), l’équation (2.12) devient :
hf ∫ z πz
mch ∫ (z) dhf = qc′ cos dz (2.13)
−L/2 Le
he
L’équation (2.13) pour un écoulement monophasique est exprimée par :

Tf ∫ z πz
mchCp ∫ (z) dT = qc′ cos dz (2.14)
−L/2 Le
Tf 1
Le résultat de l’intégration de l’équation (2.14) est :
Tf (z) − Tf1 = ′ πz πLh

qc Le sin + (2.15)
L 2 Le
mch cp π e
sin
2.3.1.2.2 Température de la gaine
La variation axiale de la température extérieure de la gaine peut être

déterminée en considérant le flux de chaleur à la surface extérieure de la
gaine :
′
h[Tc(z) − Tf (z)] = q ′ q
(z) = (z) (2.16)
Ph
Où, Ph = 2(2ec + em + lch) et h le coefficient de transfert de chaleur par
convection.

réacteur
En injectant l’équation (2.1) dans l’équation (2.16) :
′
qc πz
Tc(z) = Tf (z) c
2(2ec + em + lch)
L

réacteur
2.3.1.2.3 Température au centre de la chair
Il est possible d’extrapoler cette approche pour déterminer la température

maximale de la chair. L’équation de la distribution de la température de la
chair s’écrit : h i
1 R cs 1 1 πz
T (z) = T (z) + q ′ ln + + cos
m c
c 2πkc Rce ((2ech) + (2lch))h 4πkf Le
(2.18)
Où, km et kc sont les conductivités thermiques de la chair et de la gaine,
respectivement.
2.3.2 Analyse de sûreté nucléaire
Le bon fonctionnement d’un réacteur nucléaire de recherche est vérifié

seulement en faisant une analyse de sûreté nucléaire. Cette dernière est
définie comme étant l’ensemble des dispositifs mis en place pour veiller à
ce que la conception, la construction, et le fonctionnement, à l’arrêt et au
démarrage des réacteurs nucléaire de recherche, se fassent avec un faible
risque d’accidents.
Parmi les phénomènes critiques capables de provoquer un

dysfonctionnement au niveau du coeur, il y a le départ de l’ébullition
nucléée DNB (Departure of Nucleate Boilling ), le début de l’instabilité de
l’écoulement OFI (Onset of Flux Instability) et le début de l’ébullition
nucléé ONB (Onset of Nucleate Boilling) (SIDI ALI et al., 2012).
On traite, ici, le phénomène ONB car il s’agit d’une contrainte conservatrice.

On l’évalue en calculant les flux de chaleur moyen et critique à partir de la
corrélation de Bergles et Rohsenow. Par conséquent, le flux de chaleur
moyen au niveau de l’ONB, est calculé comme suit :
p0,0284
2
5 9
T s — f + ,23 Fax Frad
9 qav
5 TFrad
qav = 2.10 Lhlm
ρ U ech Cp lch +
Frad Fax
h (2.19)
Le flux de chaleur critique est donné par la corrélation de Mirshak tiré de
(***, 1980) :
qcr = 151(1 + 0, 01112 U )[1 + 0, 00914(Tsat − Tf1)](1 + 0, 0009 p)

27,9 Lh lm (2.20)
1− (1 + 0, 01112 U )(1 + 0, 0009 p)
ρ U ech C F d

réacteur
Ainsi que par la corrélation de Labuntsov mentionnée dans l’article (***,

1980) et qui s’écrit :
#
!0,25 "
15, 1 C (T − T )
2, 5 U 2 p s f
qcr = 145, 5 1 1 (2.2
θ( λw
!
p
Avec θ(p) = 0, 99531p0.33 1 1,33
pcr
(2.22)
−
Où, lm est la largeur du canal chauffé, ech l’épaisseur du canal d’eau, pcr la
pression critique de l’eau et λw la chaleur latente de l’eau.
Lorsqu’il y a un phénomène d’ONB il y a formation de bulles d’air . Cela

donne naissance à un écoulement diphasique dans le canal du coeur. Pour
pouvoir bien vérifier la sûreté du coeur du réacteur il faut calculer le rapport
du flux thermique critique (RFTC) qui est le rapport du flux de chaleur
critique sur la puissance axiale.
qcr
RFTC = . (2.23)
q′(z)
Le rapport de flux thermique critique doit être supérieur à 1,3 pour garantir
un bon fonctionnement du réacteur nucléaire de recherche et éviter une
ébullition nucléée comme cela a été précisé dans l’article de Boudali et al.
(2011).
2.4 Évolution des quantités physiques dans trois

réacteurs de différentes puissances
Après la mise en équation établie précédemment, on évalue pour chacun
des réacteurs nucléaires de recherche (1MW, 2MW et 10MW) la
température de la chair Tm(z), celle de la gaine Tc(z) et celle du caloporteur Tf
(z) avec les modèles cités en plus des flux de chaleur q(z). Pour se faire, un
programme informatique est établi et exploité ( voir annexe A)

nucléaire de recherche de 2MW (AIEA)
L’ensemble des résultats obtenus par les deux modèle sont présentés ci-
après.

réacteur
2.4.1.1 Résultats en utilisant le modèle d’El-Wakil (1971)
Les résultats obtenus pour un réacteur de 2 MW sont présentés sur la

figure 2.3. Cette figure regroupe les trois profils des température.
FIGURE 2.3 – Distribution de la température en fonction de la longueur adimensionnelle

du canal pour un réacteur nucléaire de recherche 2MW.
La figure (2.3) montre que la température du caloporteur augmente le long

du canal en commençant par la valeur donnée initialement qui est de 38 °C
jusqu’à ce qu’elle atteigne la valeur de 53,59°C à la sortie du canal.
On remarque que le profil de la température de la gaine a une allure

en forme de cloche. A l’entrée, elle est évaluée à 48,15 °C puis elle
augmente jusqu’à une valeur maximale de 81,8°C à une longueur
adimensionnelle de 0,56. Après ce pic, la température de la gaine diminue
pour arriver à une valeur 63,16 °C à la sortie du canal.
Quant à la température de la chair, elle vaut 48,38 °C à l’entrée. La

température de la chair augmente et suit une allure similaire à celle de la
gaine. Sa valeur maximale se situe près du centre du canal c’est à dire à
0,56 et est égale à 82,68°C. Après ce pic, la température de la chair diminue
jusqu’à la fin du canal avec une valeur de 63,42 °C.

réacteur
Pour s’assurer que les courbes de température trouvées sont cohérentes,

elles sont comparées aux résultats obtenus dans d’autres travaux. Pour ce
réacteur nucléaire de recherche on prend comme référence (Labani et al.,
2014) et les courbes du code TERMIC (1998). Ces courbes sont données sur
la figure (2.4).
FIGURE 2.4 – Comparaison de la distribution des températures avec des travaux publiés en
fonction de la longueur adimensionnelle du canal pour un réacteur nucléaire de recherche
2MW.
D’après la figure (2.4), on constate que les courbes de la température de la

gaine ont une allure assez similaire avec des températures à l’entrée et la
sortie du canal qui diffèrent. A l’entrée, la température de la gaine trouvée
par Labani et al. (2014) est de 46,34 °C, il s’agit de la température d’entrée
la plus faible avec une valeur proche de celle du présent travail avec une
différence relative de 3,75%. La température du présent travail est de 48,15
°C puis vient celle du code TERMIC qui est la température d’entrée la plus
forte parmi ces trois températures avec une valeur de 51,73 °C sa différence
relative par rapport au présent travail est de 6,92 % . A la sortie, on trouve
que la température de la gaine de Labani et al. (2014) est égale à 60°C avec
une différence relative de 5 % et celle du présent travail égale à 63,16°C
puis en dernier celle du code TERMIC avec une valeur de 64,66°C donc une
différence relative de 2,31 %.
Sachant que ces courbes sont en forme de cloche elles passent par un

réacteur
maximum près du milieu du canal. Le maximum le plus élevé est celui de

Labani et al. (2014) avec une valeur de 82,43 °C et une différence relative
de 0,76 % lorsque la valeur de l’abscisse est égale à 0,54, le deuxième pic
est celui du présent travail avec une valeur de 81,8 °C à 0,56 de la longueur
totale du canal puis vient le pic du travail effectué avec le code TERMIC
avec une valeur de 80,8 °C et une différence relative de 1,22 % à 0,54 de la
longueur adimensionnelle du canal.
De même pour les profils de température du fluide de refroidissement qui

débutent avec des températures égales à 38 ◦ C puis évoluent à la sortie telle
que la température de sortie la plus élevée est celle de (Labani et al., 2014)
avec une valeur de 59,4 °C et dont la différence relative est égale à 9,7 %
puis vient celle du présent travail avec une valeur de 53,59 °C et en dernier
celle du code TERMIC (Doval, 1998) avec une valeur de 52,17 °C et une
différence relative de 2,65 %.
On déduit de cette comparaison que les résultats obtenus avec le présent

travail sont satisfaisants puisque les différences relatives calculées pour toute
les températures significatives sont inférieurs à 5 % à part la température
de sortie du caloporteur obtenue par Labani et al. (2014).
2.4.1.2 Résultats en utilisant le modèle de Todreas et al. (1989) :
On donne les mêmes quantités physiques avec le modèle de (Todreas et

al., 1989) et on présente les résultats obtenu sur la figure (2.5).

réacteur

du canal.
Sur la figure (2.5), la température du caloporteur à l’entrée du canal est de

38 °C. Cette température va augmenter le long du canal jusqu’à atteindre la
valeur de 53,59°C à la sortie.
On remarque que la température de la gaine à l’entrée est de 47,41

°C. Cette température augmente au fur et à mesure qu’on s’approche du
centre pour ainsi atteindre une valeur maximale de 79°C à 0,56 de la
longueur adimensionnelle du canal. Après ce pic, la température de la gaine
diminue jusqu’à la sortie du canal avec une valeur de 62,5 °C.
Quant à la température de la chair, elle est égale à l’entrée à 47,5 °C.

La température de la chair augmente jusqu’à une valeur maximale. Cette
valeur maximale est de 79,35°C à 0,56 de la longueur adimensionnelle du
canal. Après ce pic, la température de la chair diminue jusqu’à sa sortie avec
une valeur de 62,6 °C.
Étant donné qu’on a vérifié les profils de température de la méthode

de (El-Wakil, 1971) aux autres travaux. Il suffit, pour ce modèle , de
comparer ces courbes avec celle du modèle (El-Wakil, 1971). On présente
sur la figure (2.6) la comparaison des profils de température trouvés avec le
modèle (El-Wakil, 1971) et le modèle (Todreas et al., 1989).

réacteur
FIGURE 2.6 – Comparaison des profils des températures des deux modèles en fonction de
la longueur adimensionnelle du canal pour un réacteur nucléaire de recherche de 2MW.
On remarque que pour le profil de température du fluide de refroidissement

les courbes sont les mêmes. Elles commencent par 38 ◦C et évoluent jusqu’à
56,62 °C à la sortie.
Contrairement aux températures de la gaine et de la chair où il y a

une légère différence telle que celle du modèle El-Wakil (1971) est
supérieure à celle du modèle de (Todreas et al., 1989). En effet, à l’entrée la
température de la gaine du modèle de Todreas et al. (1989) est égale à
47,41 °C et celle de la chair à 47,51 °C, c’est à dire qu’elle est
respectivement inférieure de 1,5 % et 1,7 % à celle du modèle de El-Wakil
(1971) qui est égale à 48,15 °C pour la gaine et 48,38 °C pour la chair.
Au pic, la température de la gaine atteint la valeur de 78,99 °C et

celle de la chair 79,35 °C tout deux à 0,58 de la longueur adimensionnelle
du canal ce qui est inférieur respectivement à 3,43 % et 4,01 % de la
température du modèle d’El-Wakil (1971) qui est égale à 81,8 °C pour la
gaine et 82,67 °C pour la chair.
A la sortie, la température de la gaine obtenue par le modèle de Todreas

réacteur
et al. (1989) est égale à 63,48 °C et celle du modèle El-Wakil (1971) égale
à 63,16 °C avec une différence de 0,5 %. Pour la chair, à la sortie, la
température obtenue par le modèle de Todreas et al. (1989) est égale à 62,59
°C et celle du modèle El-Wakil (1971) à 63,42 °C avec une différence
évaluée à 1,3 %.
Puisque les différences relatives sont inférieures à 5 % on en déduit

que le deuxième modèle présente des résultats très acceptables.
2.4.1.3 Résultats de l’analyse de sûreté nucléaire
Pour éviter qu’il y ait un phénomène ONB, on calcule le flux de

chaleur moyen qav et le flux de chaleur critique qcr et le rapport du flux
thermique critique RFTC en fonction de z.
On calcule le flux de chaleur moyen à l’aide de la corrélation de Bergles

et Rohsenow (***, 1980), qui donne lieu à une équation non linéaire qui
a été résolue en utilisant la méthode du point fixe. On trouve le flux de
chaleur moyen égal à qav = 15, 57W/cm2. On calcule également le flux
de chaleur critique avec la corrélation de Mirshak (***, 1980), on trouve
qcr = 239, 73W/cm2 à T=Tf 1 =38 ◦ C.
Le tableau (2.5) affiche les flux de chaleur lors d’un ONB trouvé par (Khedr,
2008) et ceux trouvés dans le présent travail pour T = Tf1.
Tableau 2.5 – Flux de chaleur pour un réacteur 2 MW.
Travail Flux de chaleur moyen ONB Flux de chaleur critique ONB

[W/cm2] [W/cm2]
Présent travail 15,57 239,73
référence Khedr (2008) 13,58 230
On remarque que les deux valeurs de flux de chaleur critique et moyen

sont proches de leurs valeurs de référence tel que le flux de chaleur moyen
trouvé dans le présent travail diffère du flux trouvé dans (Khedr, 2008) de
12,8% contrairement au flux critique dont la valeur s’approche de la valeur
de référence tiré du travail de (Khedr, 2008) avec un pourcentage de 4,06 %.
Par conséquent les résultats sont satisfaisants.

réacteur
Dans la figure (2.7), on donne la variation du flux de chaleur critique

le long du canal.
FIGURE 2.7 – Flux de chaleur critique ONB en fonction de la longueur adimensionnelle

Le flux de chaleur critique à l’entrée est égal à qcr = 239, 73W/cm2. En

allant vers la sortie du canal le flux de chaleur diminue jusqu’à une valeur de
qcr = 218W/cm2. On déduit que le flux de chaleur critique est inversement
proportionnel à la température du caloporteur.
On présent sur la figure 2.8 le rapport du flux de chaleur critique en chaque

point du canal.

réacteur
FIGURE 2.8 – Rapport du flux de chaleur critique ONB le long du canal pour un réacteur
nucléaire de recherche 2MW.
Ce qui intéresse la sûreté nucléaire sur la figure (2.8) c’est la plus faible
valeur atteinte. Dans notre cas la valeur minimale se situe près du centre du
canal et est est égale à 3. Cette valeur est supérieure à 1,3 ce qui permet de
déduire que le phénomène d’ébullition (ONB) est évité.

nucléaire de recherche de 10MW (AIEA)
Étant donné que l’analyse effectuée pour ce réacteur est la même que
pour le réacteur de 2MW, les résultats sont présentés de la même façon. On
présente ci-après les résultats des deux modèles.
2.4.2.1 Résultats en utilisant le modèle El-Wakil (1971)
Les résultats obtenues pour un réacteur de puissance 10 MW en utilisant

le modèle d’El-Wakil (1971) :
Les trois profils de température le long du canal sont présentés sur la figure
(2.9).

réacteur
FIGURE 2.9 – Distribution des températures en fonction de la longueur adimensionnelle

Dans la figure (2.9) la température du caloporteur à l’entrée du canal est

égale à 38 °C. Cette température va augmenter le long du canal jusqu’à
atteindre la valeur de 53,55°C à la sortie.
A l’entrée, la température de la gaine vaut 48,88 °C. Cette température

augmente jusqu’à une valeur maximale de 83,96°C à 0,56 de la longueur
adimensionnelle du canal. Après ce pic, elle diminue jusqu’à la sortie du
canal avec une valeur de 63,65 °C. Quant à la température de la chair, elle
vaut à l’entrée 49,57 °C. La température de la chair augmente jusqu’à 86,54
°C à 0,54 de la longueur adimensionnelle du canal puis elle diminue jusqu’à
sa sortie pour atteindre 64,39 °C.
Pour vérifier nos résultats, trois références seront utilisées Labani et al.
(2014), TERMIC (1998) et (Lu et al., 2009). Toutes les courbes sont
présentées sur la figure (2.10).

réacteur
FIGURE 2.10 – Comparaison des distributions de température en fonction de la longueur

adimensionnelle du canal pour un réacteur nucléaire de recherche 10MW.
En comparant d’abord les profils de température du caloporteur, on voit

qu’à l’entrée du canal la température est égale à 38 °C puis à la sortie
la température donnée par TERMIC (1998) est égale à 54,52 °C, celle
du présent travail vaut 53,55°C, celle de Labani et al. (2014) est égale à
50,87°C puis celle de Lu et al.(2009) égale à 49,62°C. On remarque que la
température de sortie la plus élevée est celle du code TERMIC (1998), il
s’agit notamment de la courbe qui se rapproche le plus de la courbe du
présent travail avec un pourcentage de 1,7 % de différences relatives.
Pour le profil de la température de la gaine, la valeur à l’entrée du

présent travail est égale à 48,88 °C alors que celle de Labani et al.(2014) est
égale à 38°C, celle de Lu et al.(2009) égale à 38,07°C et celle de TERMIC
(1998) égale à 49,87°C. La température maximale de la gaine vaut 83,96°C
pour le présent travail à 0,56 de la longueur adimensionnelle du canal, elle
est de 83,42°C pour Labani et al.(2014) à 0,67 de la longueur totale du canal
avec une différence de 0,64 % par rapport au présent travail, pour Lu et
al. (2009) elle est égale à 77,3°C à 0,64 de la longueur adimensionnelle du
canal et présente une différence relative par rapport au présent travail de
7,9 % puis pour TERMIC (1998) elle est de 79,97°C à 0,59 de la longueur
adimensionnelle du canal avec une différence relative de 4,75 % par rapport

réacteur
au présent travail. A la sortie les valeurs diminuent, pour le présent travail on

trouve une valeur de 63,65°C, pour Labani et al (2014) 68,88°C, pour Lu et
al.(2009) 65,53°C puis pour TERMIC (1998) on trouve une valeur de 65,53°C.
Au vu des différents profils et différentes différences relatives on peut dire

que les valeurs obtenues sont très acceptables et la méthode utilisée est
fiable.
2.4.2.2 Modèle de Todreas et Kazimi (1990)
Les profils de température calculés avec le modèle de Todreas et Kazimi

(1990) sont donnés sur la figure (2.11) :

Sur la figure (2.11), la température du caloporteur à l’entrée du canal vaut

38 °C. Cette température va augmenter le long du canal pour atteindre
la valeur de 56,36°C. La température de la gaine à l’entrée est égale à
53,08 °C par la suite elle augmente au fur et à mesure qu’on s’approche
du centre à 0,54 de la longueur adimensionnelle du canal pour atteindre
une valeur maximale de 98,7°C. Après ce pic, la température de la gaine
diminue jusqu’à la sortie du canal à 73,57 °C. Concernant la température de
la chair, à l’entrée elle vaut 53,53 °C. La température de la chair augmente
jusqu’à une valeur maximale. Cette valeur maximale est de 100°C à 0,54 de

réacteur
la longueur adimentionnelle du canal. Après ce pic, la température de la

chair diminue jusqu’à la sortie du canal et sa valeur descend à 74,18 °C.
Étant donné qu’on a vérifié les profils de température obtenus par la méthode
d’El-Wakil(1971) avec d’autres travaux. Il suffit, pour le modèle de Todoreas
et Kazimi(1940) de le comparer ses courbes avec celles du modèle d’El-Wakil.
Les résultats obtenus par les deux modèles sont présentés sur la figure (2.12) :
FIGURE 2.12 – Comparaison de la distribution des températures des deux modèles en

10MW.
On remarque pour le profil de température du fluide de refroidissement que

les allures des courbes sont identiques. Elles commencent par une valeur
de 38 ◦C puis à la sortie elles ont une valeur de 56,62 ◦ C. Concernant
les températures de la gaine et de la chair il y a une grande différence en
amplitude telle que celle du modèle de Todreas et Kazimi (1990) arrivée au
pic à 0,54 elle devient égale à 98,7°C pour la gaine et 100°C pour la chair et
est supérieure à celle du modèle de El-Wakil (1971) qui atteint la valeur de

réacteur
83,96°C pour la gaine et 86,54°C pour la chair à 0,56 de la longueur totale

du canal.
A l’entrée, la température de la gaine du modèle de Todreas et Kazimi

(1990) est égale à 53,08°C et à la sortie elle à 73,57°C comparé à la
température trouvée avec le modèle El-Wakil(1971) qui, à l’entrée, est égale
à 53,55°C et à la sortie égale à 63,65°C.
Pour la température de la chair à l’entrée elle est égale à 53,53°C et à

la sortie elle est égale à 74,18°C pour le modèle de Todreas et Kazimi
(1990). Pour le modèle de El-Wakil la température de la chair à l’entrée est
égale à 49,57°C et à la sortie égale à 64,39°C.
2.4.2.3 Résultats de l’analyse de sûreté nucléaire
La même analyse faite au paragraphe 2.3.1.2 est appliquée ici.
On a calculé le flux de chaleur moyen à l’aide de la corrélation de

Bergles et Rohsenow (***, 1980), qui donne lieu à une équation non
linéaire qui a été résolu en utilisant la méthode du point fixe. On trouve
le flux de chaleur moyen égal à qav = 44, 7W/cm2. On calcule également
le flux de chaleur critique avec la corrélation de Mirshak ***, 1980, on
trouve qcr = 245, 1W/cm2 à T=Tf 1 =38 ◦ C.
Le tableau (2.6) affiche les flux de chaleur lors d’un ONB trouvés par
(Khedr, 2008) et ceux trouvé dans le présent travail pour T = Tf1.
Tableau 2.6 – Flux de chaleur d’un réacteur 10 MW.
Travail Flux de chaleur moyen ONB Flux de chaleur critique ONB

[W/cm2] [W/cm2]
Présent travail 44,7 245,1
référence (Khedr, 2008) 41,97 265,6
On remarque que les deux valeurs de flux de chaleur critique et moyen

sont proches de leurs valeurs de référence tel que le flux de chaleur moyen
trouvé dans le présent travail diffère du flux trouvé dans Khedr, 2008 de
6,11%. Concernant le flux critique dont la valeur s’approche de la valeur de
référence tiré du travail de (Khedr, 2008) avec un pourcentage de 7,71 %.
Par conséquent les résultats sont acceptables.

réacteur
Dans la figure (2.13), on donne la variation flux de chaleur critique.

du canal pour un réacteur nucléaire de recherche de 10MW.
Le flux de chaleur critique à l’entrée est égal à qcr = 245, 1W/cm2. En allant
vers la sortie du canal le flux de chaleur diminue jusqu’à une valeur de
qcr = 212W/cm2. On déduit que le flux de chaleur critique est inversement
proportionnel à la température du fluide de refroidissement.
On présente sur la figure (2.14) le rapport du flux de chaleur critique en

chaque point du canal.

réacteur
nucléaire de recherche de 10MW.
Sur la figure (2.14) la plus faible valeur du rapport du flux thermique

critique est égale à 1,6. Dans ce cas le phénomène d’ébullition (ONB) ne se
produira pas (RFTC> 1,3).

nucléaire de recherche de 1MW (NUR)
On présente les résultats des deux modèles mathématique El-Wakil

(1971) et Todreas et al. (1989).
2.4.3.1 Résultats en utilisant le modèle El-Wakil (1971)
Les résultats obtenus pour un réacteur de 1 MW en utilisant le modèle

d’El-Wakil (1971) sont donnés sur la figure (2.15).

réacteur

du canal pour un réacteur nucléaire de recherche de 1MW.
Dans la figure (2.15) la température du fluide de refroidissement à l’entrée

du canal est égale à 38°C . Cette température va augmenter le long du canal
jusqu’à atteindre la valeur de 48,9°C à la sortie.
On remarque que la température de la gaine à l’entrée est de 44,58°C.

La température de la gaine augmente jusqu’à une valeur maximale de
66,36°C à 0,58 de la longueur adimensionnelle du canal. Après ce pic, la
température de la gaine diminue jusqu’à la sortie du canal à une valeur de
54,86°C.
Quant à la température de la chair, elle est égale à l’entrée à 44,89°C.

Elle augmente jusqu’à 68°C à 0,56 de la longueur adimensionnelle du canal
puis elle diminue jusqu’à la sortie pour atteindre 55,33°C.
2.4.3.2 Résultats en utilisant le modèle Todreas et al. (1989)
Les profils de température calculés avec le modèle de (Todreas et al.,

1989) sont présentés sur la figure (2.16) :

réacteur

On remarque sur la figure (2.16) que la température du caloporteur à l’entrée

du canal est de 38°C. Cette température augmente le long du canal jusqu’à
atteindre la valeur de 48,98°C à la sortie.
On remarque que la température de la gaine à l’entrée est de 43,44°C

. Cette température augmente au fur et à mesure qu’on s’approche du centre
pour ainsi atteindre une valeur maximale de 61,9°C à 0,56 de la longueur
adimensionnelle du canal. Après ce pic, la température de la gaine diminue
jusqu’à la sortie du canal avec une valeur de 53,6°C.
Quant à la température de la chair, elle est égale à l’entrée à 43,4°C.

La température de la chair augmente jusqu’à une valeur maximale de 62,6°C
à 0,56 de la longueur adimensionnelle du canal. Après ce pic, la température
de la chair diminue jusqu’à la sortie du canal avec une valeur qui est de 53,8°C.
Comme pour les deux cas précédents à 2MW et à 10MW, on a aussi

comparer ces résultats avec ceux obtenus par le modèle d’El-Wakil (1971).

réacteur
FIGURE 2.17 – Comparaison de la distribution des températures des deux modèles en

de 1MW.
On remarque pour le profil de température du fluide de refroidissement que

les courbes sont identiques. Elles commencent par une valeur de 38◦C puis
à la sortie elles ont une valeur de 48,98°C.
Concernant les températures de la gaine et de la chair il y a une différence

en amplitude comme celle du modèle de Todreaset Kazimi (1990) arrivée
au pic à 0,6 elle est égale à 61,9°C pour la gaine et 62,6°C pour la chair et
elle est inférieure à celle du modèle d’El-Wakil (1971) qui atteint la valeur
de 63,76°C pour la gaine et 65,43°C pour la chair à 0,6 de la longueur
adimensionnelle du canal.
A l’entrée, la température de la gaine du modèle de Todreas et Kazimi(1990)

est égale à 43,44°C et à la sortie elle vaut 53,6°C comparé à la température
trouvée avec le modèle d’El-Wakil(1971) qui, à l’entrée, est égale à 44,5°C
et à la sortie 50,43°C.
Pour la température de la chair à l’entrée elle est égale à 44,89°C et à

réacteur
la sortie elle est égale à 50,89°C pour le modèle d’El-Wakil . Pour le modèle
de Todreas et Kazimi(1990) la température de la chair à l’entrée est égale à
43,4°C et à la sortie égale à 53,8°C.
2.4.3.3 Résultats de l’analyse de sûreté nucléaire :
La même analyse faite pour les cas de 2MW et 10MW est appliquée
ici. On trouve le flux de chaleur moyen égal à une valeur de qav = 35,
65W/cm2. On calcule également le flux de chaleur critique avec la
corrélation de Mirshak (***, 1980), on trouve qcr = 263, 24W/cm2 à T=Tf 1
=38 ◦ C.
Dans la figure (2.18), on trace la variation du flux de chaleur critique.

Le flux de chaleur critique à l’entrée est égal à qcr = 263, 1W/cm2. A la

sortie le flux diminue jusqu’à une valeur de qcr = 254W/cm2.
On donne sur la figure (2.19) le rapport du flux de chaleur critique en chaque

point du canal.

réacteur
nucléaire de recherche de 1MW.
Pour vérifier la sûreté nucléaire il faut que le rapport RFTC soit supérieur
à la valeur 1,3. Sur la figure (2.19) la plus faible valeur est égale à 7,5. Dans
ce cas le phénomène d’ébullition (ONB) ne se produit pas.
2.5 Conclusions
Dans ce chapitre une analyse thermohydraulique en utilisant le modèle
d’El-Wakil (1971) et le modèle de Todreas et al. (1989) a été faite. Ces deux
modèles proposent des solutions à l’équation de l’énergie. La différence
entre les deux modèles se situe dans les termes et les constantes dans
l’expression des équations donnant les différentes températures.
Les résultats obtenus à partir de ces deux modèles ont été vérifiés en
se basant sur des travaux valorisés, on cite les travaux de Labani et al.
(2014) (Labani et al., 2014), Lu et al. (2009) (Lu et al., 2009) et TERMIC
(1998)(Doval, 1998) et après vérification on peut dire que les résultats sont
satisfaisants puisque les différences relatives entre les résultats du présent
travail et ceux des travaux précédents n’excèdent pas les 5 %. De plus, on
a appliqué ces deux modèles sur trois réacteurs nucléaire de recherche de
puissance distincts 1MW, 2MW et 10MW.

réacteur
Par la suite, on a examiné la sûreté nucléaire des trois réacteurs par le

calcul du flux de chaleur moyen et critique et le rapport de flux de chaleur
critique en utilisant la corrélation de Bergles et Rohsenow et la corrélation
de Mirshak. On a déduit que les réacteurs nucléaire étudiés fonctionnent en
toute sûreté car les résultats de calculs ne dépassent pas la limite critique .
Cette analyse de sûreté nucléaire est faite dans le but de vérifier et d’éviter
qu’il n’y ait un phénomène d’ébullition qu’on appelle début d’ébullition
nuclée (ONB). Ce travail a été effectué au moyen d’un programme Matlab,
qu’on a mis au point et qu’on peut adapter à n’importe quel réacteur
nucléaire de recherche utilisant un combustible de type plaque, il suffit de
changer les données d’entrée spécifique au réacteur nucléaire étudié.
On peut finir par dire que ces méthodes asymptotiques sont très utiles
pour la compréhension de la thermohydraulique des réacteurs nucléaires de
recherche et une bonne base pour aboutir à des résultats exploitables en des
temps limités.

CHAPITRE 3
ModéLIsaTIon THERMohYDRaULIqUE
en 2D d’Un canaL DU coeUR D’Un
RÉacTEUR nUcLÉaIRE DE REcheRche
UTILIsanT Un combUsTIBLE à fLUx de
chaLEUR UnIFoRM ET consTanT
E
Chapitre 3. Modélisation thermohydraulique en 2D d’un canal du
coeur d’un réacteur nucléaire de recherche utilisant un combustible à
flux de
3.1 Introduction
Dans le présent chapitre, On va s’intéresser au cas bidimentionnel de
l’écoulement du fluide caloporteur dans un canal du réacteur nucléaire
de recherche. Deux mises en équation seront établis, la première en mode
stationnaire et la seconde en mode transitoire.
En fait, on va reproduire l’évaluation expérimentale dans un réacteur

nucléaire de recherche avant insertion du combustible (Uranium) à l’aide
du logiciel COMSOL. Le combustible utilisé pour évaluer les performances
expérimentales avant le premier démarrage du réacteur est un combustible
à base de résistances électriques. Donc une source de chaleur uniforme et
constante. Le but recherché est de savoir jusqu’à quelle température du
fluide caloporteur on peut arriver pour un fonctionnement normal ou
accidentel.
3.2 Configuration du canal thermohydraulique en 2D

L’écoulement dans le canal se fait suivant l’axe y de haut en bas. Compte
tenu du fait qu’il ait une forme symétrique, il est possible de configurer une
moitié seulement comme cela est montré sur la figure (3.1). L’axe de
symétrie coupe la géométrie au niveau du fluide caloporteur.

flux de

flux de
FIGURE 3.1 – Géométrie du canal en 2D.
3.3 Mise en équation thermohydraulique

Afin de bien modéliser le problème, on introduit les équations qui le
gouvernent et qui sont respectivement l’équation de conservation de la
masse, de la quantité de mouvement et de l’énergie. On rajoute un modèle
de turbulence pour considérer les écoulements turbulents.
3.3.1 Les équations de conservation
Les équations de conservation sont énoncé ci-après (Sidi Ali ; Djaber

et al., 2017) :
Équation de conservation de la masse
L’équation de conservation de la masse s’écrit :

∂ρ
+ . (ρU→ ) = 0 (3.1)
∂t ∇

flux de
Où, ρ est la masse volumique, U→ le vecteur vitesse et t le temps.
Équation de conservation de quantité de mouvement
L’équation de conservation de la quantité de mouvement est donnée

par :
∂
(ρU→ ) + ∇ . (ρU→ U→ ) = −∇p + ∇ . + (∇U→ )T )] + ρ→g (3.2)
∂ [µ(∇U→
Où, p est la pression, µ la viscosité dynamique du fluide et →g le

vecteur gravité.
Équation de conservation de l’énergie
L’équation de conservation de l’énergie est comme suit :

h i h Σ i
∂ρE
→ = ∇. (k∇T ) − →
∂t + ∇. (ρE + p)U J + (τ¯.U ) + S (3.3)
→j
hj
j
Où, E est l’énergie, k la conductivité thermique, T la température, hj

le coefficient de transfert de chaleur par convection d’une espèce j , J→j le
vecteur vitesse d’une espèce j, (τ¯ le tenseur des contraintes visqueuses et S
la source
∂
A noter qu’en stationnaire les dérivés par rapport au temps sont
∂
nulles.
3.3.2 Modélisation de la turbulence
L’écoulement considéré est un écoulement turbulent étant donné que le

nombre de Reynolds est supérieur à 3000.Ainsi pour prédire les écoulements
turbulents, trois approches peuvent être utilisées : Simulation des grands
tourbillons (LES : Large Eddy Simulation), la simulation numérique directe
(DNS : Direct Numerical Simulation) et le modèle Navier-Stokes moyenné
par Reynolds (RANS : Reynolds-Averaged Navier-Stokes) tels que
présentés par (Chen et al., 2007).
flux de
Afin de mieux comprendre l’effet de la turbulence sur l’écoulement on

commence par décomposer la vitesse en deux termes, un terme moyen et

flux de
un terme fluctuant comme expliqué par (Wilcox, 1993) :

′
ui = ui + u i (3.4)
Où, u¯i est la vitesse moyenne et u′i la vitesse fluctuante(i=1, 2, 3).
En remplaçant l’équation (3.4) dans l’équation (3.2) et après calcul des

moyennes, on obtient l’équation (3.5) :
∂ ∂ ∂p ∂ h i
(ρu )+ (ρu iu j) = − + ∂ui ∂uj 2 ij∂ul
∂ ′ ′
i ∂ + ∂ (−ρu
µ ∂ + ∂ − 3δ ∂x u)
∂ ∂ ∂ l
i j
(3.5)
La sommation des indices (convention d’Einstein) est considérée et pour
simplifier les expressions, les barres sur les vitesses sont supprimées.
L’équation (3.5) est appelée équation de Navier-Stokes moyennée par

′ ′
Reynolds (RANS), −ρ u ui jest appelé les contraintes de Reynolds, et doit
être modélisée. Hinze et al. (1975) a relié les contraintes de Reynolds au
gradient de vitesse moyen, on obtient l’équation suivante :
′ ′
— ρu u ∂ui ∂ui 2 ∂ul (3.6)
i j= µ + − ρK + µ δ
t t i
∂ ∂ 3
Où, µt est la viscosité turbulente, K l’énergie cinétique turbulente
La famille RANS contient plusieurs modèles de turbulence. On cite le

modèle
— K ε standard (Launder et al., 1974), modèle K ε Re-
Normalization Group (RNG) (Yakhot et al., − 1992), modèle K ε
Réalisable (Shih et al., 1995), modèle
− standard de K ω (Wilcox, 1993)
et le modèle choisit dans − nos simulations K ω SST (Shear-Stress
Transport) (Menter, 1994) développé ci après.
Modèle K − ω SST (Shear-Stress Transport)(Menter,1994)
Le modèle Shear Stress Transport K −ω, proposé par Menter en 1994

est un modèle hybride, c’est à dire utilisant des modèles de turbulence
différents pour la région proche de la paroi et dans le reste de l’écoulement.
La transition d’un modèle vers l’autre se fait progressivement à partir d’un
mélange. L’approche K −ω est utilisée pour prédire l’écoulement dans la
région affectée par la viscosité, et possède l’avantage de ne pas nécessiter
−
flux de
de fonction d’amortissement, alors que l’approche K ε est utilisée pour
l’autre partie de l’écoulement. Une fonction permet d’assurer la transition

flux de
d’un modèle à l’autre. Les équations de transport du modèle SST K − ω

sont similaires aux équations utilisées par le modèle K − ω standard.
La viscosité turbulente et les équations de transport appliquées dans le

modèle K − ω SST sont comme suit :
ρK 1
µt = h √ (3.7)
ω 1 2Ωij Ωij i
max F2
α∗ α1ω
Où Ωij est le tenseur du taux de rotation moyen et α1 = 0, 31. F2 est une

fonction de mélange Menter 1994. Les équations de transport de ce modèle
sont :
h i
∂ ∂ µt ∂K
(ρK) + ∂ µ + (3.8)
(ρKu ) = — β ρωK +
∗
∂ ∂
G
i K
∂ ∂ ∂ σ ∂
(ρ ω) + ∂ h
∂ ∂ (ρ ω u ) =
γ 2 µt ∂ω i
i
G K— β ρ ω + µ+ + (3.9)
ν
D
ω
∂ GK due
La génération de l’énergie cinétique turbulente σ aux
∂ gradients de
vitesse moyenne est donnée par :
′ ′ ∂uj
G = −ρu u
(3.10)
K i j
∂xi
La quantité Dω est le terme de diffusion croisé et est donné par :
ρ ∂K ∂ω
D = 2(1 − F σω,2 ω ∂xj ∂xj
)ω 1
Les quantités σk et σω sont données par :
1
σK = F1 σK,1 +(1−F1 )σK,2 avec σK,1 = 0, 85, σK,2 = 1, 0
1
Et σω = avec σω,1 = 0, 5, σω,2 = 0, 856 et β∗ = Cµ =
F1σω,1+(1−F1)σω,2
0, 09
Où, µt est la viscosité turbulente, σk, σω, β , ∗β. La viscosité statique
turbulente est donnée par :
µt flux de
νt = (3.11)
ρ
Enfin, F1 apparaissant dans l’expression qui donne σK est une fonction de
mélange comme F2. De plus amples explications sur la modélisation des
différentes quantités peuvent être trouvées dans Menter (1994).

flux de
3.4 Analyse thermohydraulique à l’état stationnaire
3.4.1 Mise en équation
Étant donné que la modélisation concerne un écoulement turbulent où un

échange de chaleur est effectué entre le fluide et les plaques de combustible
qu’on appelle convection forcée, alors les équations qui rentrent en jeu sont
les équations de conservation de la masse, de la quantité de mouvement et
de l’énergie (COMSOL, 2020).
3.4.1.1 Les équations de conservation
a) L’équation de conservation de la masse s’écrit comme suit :
ρ∇.U→ = 0 (3.12)
b) L’équation de conservation de la quantité de mouvement est donnée par

(3.13) :
ρ(U→ .∇)U→ = ∇.[−p| + τ¯] + F→ + ρ→g (3.13)
Où, U→ est le vecteur vitesse du fluide, p la pression, τ¯ tenseur des
contraintes visqueuses, ρ la masse volumique, le vecteur force de
F→ volume et →g le vecteur gravité.
et :
τ¯ = (µ + µt) + (∇U→ )T ) (3.14)
(∇U→
Où, µ est la viscosité dynamique du fluide et µt est la viscosité
turbulente.
c) L’équation de conservation de l’énergie est donnée par (3.15) :
ρ Cp U→ .∇ T + ∇. →q = S + q0 (3.15)
Où, S est la source de chaleur, →q est le vecteur du flux de chaleur, Cp
la chaleur spécifique du caloporteur, q0 flux de chaleur entrant.
Le flux de chaleur est donné dans l’équation (3.16).
q = −k ∇ T (3.16)

flux de
Où, k est la conductivité thermique.

flux de
On adapte l’équation (3.15) sur la gaine, la chair et le caloporteur.

Ainsi, pour la chair et la gaine l’équation se réduit à (3.15) se réduit à
:
ρ Cp U→ . ∇ T + ∇ . →q = S + q0 + Qted (3.17)
Où, Qted est l’amortissement thermo-élastique.
Pour le caloporteur, l’équation (3.15) s’écrit alors :
ρ Cp U→ .∇ T + ∇ . →q = S + q0 + Qp + Qvd (3.18)
Où Qp est la source de chaleur ponctuelle et Qvd la dissipation visqueuse.
3.4.2 Conditions aux limites
On introduit la vitesse d’entrée au canal qui est de 0.94 m/s, puis on

pose les conditions aux limites à la paroi, à l’entrée et à la sortie. Puisque
la géométrie est symétrique, on rajoute la condition de symétrie. On donne
aussi les conditions liées au transfert de chaleur entre la chair, la gaine et
le fluide où le combustible expérimental représente la source de chaleur
(COMSOL, 2020).
a) A la paroi
U |lw =0 = 0 (3.19)
h1
l = (3.20)
w
2
∇k.→n = 0, (3.21)
Où, lw est ,→n le vecteur normal
b) A l’entrée du canal :
On a :
ρ → M
∫ (U .→n)dA = , ρst = n (3.22)
ρst Vm
— Qsv
∂ 1
3 2 k2
, ω= , ∇G.→n = 0 (3.23)
k = (Uref IT ) ∗1 4
(β0 )
2
LT
Où, Qsv est le débit standard du caloporteur, Vm le volume molaire
standard du caloporteur, Mn la masse molaire moyenne du caloporteur,
flux de
IT , LT , β0∗ , G la distance résiduelle à la paroi, ρst la masse
volumique standard du caloporteur.

flux de
c) A la sortie du canal
On a :
∧ 1
p = p − ρ|U→ |
2
hydro (3.24)
+p 2
∧
p ≤ p0 , U→ . t = 0 (3.25)
∇k.→n = 0 , ∇ω.→n = 0 , ∇G.→n = 0 (3.26)
∧
Où, p0 est la pression totale, p la pression adimentionnelle, phydro
d) Consernant la symétrie on a
U→ .→n = 0 (3.27)
k→n − (k→n.→n)→n = 0 , (3.28)
∇k.→n = 0 , ∇ω.→n = 0 , ∇G.→n = 0 (3.29)
e) Température du fluide à l’entrée au canal
La température du fluide est fixée à 38°C à l’entrée du canal.
f) Source de chaleur
La source de chaleur est donnée sous forme de puissance thermique

uniforme et constant.
q′ = 10, 38.104W/m3 (3.30)
3.5 Analyse thermohydraulique à l’état transitoire

On passe à l’état transitoire où la variable du temps rentre en jeu . On
réécrit les équations (3.1) et (3.14) .
3.5.1 Mise en équation
Les équations qui décrivent l’écoulement cité ci-après ont été fourni du
guide COMSOL (COMSOL, 2020).

flux de
3.5.1.1 Les équations de conservation
a) L’équation de conservation de la masse s’écrit :

∂ρ
+ρ =0
.
Tel que :
U→
∂t ∇ + (∇U→ )T )
τ¯ = (µ + µT )
(∇U→
Où, µ est la viscosité dynamique et µT est la viscosité de turbulence.
b) L’équation de conservation de la quantité de mouvement devient :

∂
U + ρ(U→ = ∇.[−p| + τ¯] + + ρ→g (3.31)
→
.∇)U→ F→
ρ
∂t
On développe l’équation (3.31) par rapport à l’axe x et y. Ainsi, on

donne l’équation (3.32) qui est l’équation de conservation de la quantité
de mouvement suivant x :
∂u ∂u ∂u ∂p (3.32)
ρ +u +v =− + τ¯ + F→ + x
∂ ∂ ∂ ∂
ρg
L’équation de quantité de mouvement suivant y s’écrit :
∂v ∂v
ρ +u ∂v ∂p
+v =− + τ¯ + + y
(3.33)
∂ ∂ → ρg
F
∂ ∂
c) L’équation de conservation de l’énergie est donnée par l’équation (3.34) :
∂T
ρ Cp + ρ Cp U.∇T + ∇. q = S + q0 (3.34)
∂t
3.5.2 Conditions aux limites

flux de
Les conditions aux limites restent exactement les mêmes qu’en situation
stationnaire. Excepté pour la condition à l’entrée du canal où la vitesse
diminue de 20% de sa valeur initiale et devient donc égale à :
U = 0, 752m/s (3.35)

flux de
3.6 Présentation du logiciel Comsol

COMSOL Multiphysics est un logiciel avancé de simulation polyvalent
pour la modélisation de designs, de dispositifs et de procédés dans tous les
domaines de l’ingénierie, de la production et de la recherche scientifique et
pour tout processus multiphysique . Il a remplacé FEMLAB qui était pour le
coup un logiciel de modélisation très puissant, basé à l’origine sur le langage
de programmation MATLAB et l’environnement de développement intégré .
Comsol comprend deux interfaces qui sont Comsol multiphysics et Comsol

script. La première concerne l’interface graphique et la deuxième montre
les lignes de commande. Sur la figure (3.2) on montre l’interface Comsol
multiphysics lors de la modélisation (Comsol, 2020).
FIGURE 3.2 – interface Comsol multiphysics.
3.7 Maillage du domaine numérique

Un Maillage est un ensemble d’éléments de volumes unitaires utilisés
pour discrétiser un domaine de calcul ou une géométrie. Il convient de
distinguer plusieurs types de maillages selon leurs caractéristiques
géométriques. Un maillage sert à découper une géométrie 2D ou 3D en
plusieurs éléments pour réaliser des études numériques et calculer en tous
points de la géométrie les variables dont on veut connaître les valeurs. Ces
calculs sont effectués à l’aide de méthodes numériques qui dépendent du
choix du maillage.
Il convient de distinguer plusieurs types de maillages selon leurs
flux de

flux de
caractéristiques géométriques. Il existe deux types de maillages, structuré ou

non structuré. Il peut être simple, de forme triangulaire ou tétraèdre ou bien
mixte en combinant plusieurs type d’éléments structurés et non structurés.
Chaque code de calcul utilise un maillage spécifique. Ainsi les codes de
calcul d’éléments structurés utilisent des quadrilatères en (2D) et des
hexaèdres en (3D). Les solveurs d’éléments non structurés utilisent souvent
des triangles (2D) ou des tétraèdres (3D), mais les codes plus récents
peuvent utiliser des polyèdres arbitraires. Les représentations de éléments de
maillage sont donnés sur la figure (3.3) (Sadrehaghighi, 2021).
FIGURE 3.3 – Type d’éléments de maillage.
Dans notre cas, le type de maillage choisi est un mélange de mailles

triangulaires 2D non structuré et de mailles quadrilatères structuré. Les
mailles triangulaires concernent la gaine, le combustible et le fluide loin de
la paroi alors que les mailles quadrilatères concerne la couche limite
laminaire. La figure (3.4) montre une partie du maillage généré par
COMSOL avec les différents encombrements.

flux de
FIGURE 3.4 – Maillage du domaine numérique.
3.8 Méthode des éléments finis

Comsol utilise la méthode des éléments finis pour la résolution des
équations différentielles et intégrales qui se posent dans divers domaines de
l’ingénierie et des sciences appliquées. La méthode est une généralisation
des méthodes variationnelles classiques ( la méthode de Ritz) et des
méthodes résiduelles pondérées (par exemple, moindres carrés), qui reposent
sur l’idée que la solution u d’une équation différentielle peut être représentée
comme une combinaison linéaire de paramètres inconnus cj et de
fonctions φj choisies de manière appropriée dans le domaine entier du
problème. Les paramètres cj sont alors déterminés de telle sorte que
l’équation différentielle soit satisfaite, généralement dans un sens intégral
pondéré. Les fonctions φj , appelées fonctions d’approximation, sont
choisies de telle sorte qu’elles satisfassent les conditions aux limites du
problème.
L’idée de base de la méthode des éléments finis est de considérer un

domaine donné comme un assemblage de formes géométriques simples,
appelées éléments finis pour lesquelles il est possible de générer
systématiquement les fonctions d’approximation nécessaires à la résolution
d’équations différentielles. La capacité de représenter des domaines à
géométrie irrégulière par un ensemble d’éléments finis fait de cette méthode
un outil pratique pour la résolution (Reddy et al., 2010).

flux de
La figure (3.5) montre comment est discrétisé un domaine par la méthode

des éléments finis.
FIGURE 3.5 – discrétisation d’un domaine par éléments finis.
Sur la figure (3.5), Γ représente les frontières du domaine, Ωe un élément et

Γe ses frontières.
3.9 Résultats obtenus pour l’écoulement stationnaire

On tire de la simulation les résultats des températures en fonction de l’axe
x et de l’axe y. Puisque la géométrie est en bidimensionnel, il est difficile
de juger correctement l’évolution de la température dans le canal. Donc
pour bien la visualiser, il suffit de moyenner les valeurs par rapport à l’axe
x pour obtenir la courbe de distribution de température en fonction de l’axe
y.
Pour moyenner les valeurs des températures dans un domaine 2D on le

divise en plusieurs sections de dimensions égales. Tel que pour une section,
on effectue une moyenne de toutes les valeurs des noeuds se trouvant dans le
domaine. Cette moyenne se calcule en intégrant la surface de la section entre
les bornes a et b. Un code Matlab a été élaboré à cet effet (voir le script dans
l’annexe B)
∫
b
f (x, y) = 1 f (x, y)dxdy (3.36)
b−a a

flux de
Les courbes de températures sont fonction de l’axe y dans l’intervalle de la

flux de
longueur du combustible. La courbe de distributions de température du

fluide de refroidissement est donnée sur la figure 3.6.
FIGURE 3.6 – Distribution de la température du fluide de refroidissement le long de l’axe

y.
On remarque sur la figure (3.6) une augmentation de la température de 40,1

°C à l’entrée du canal jusqu’à 47 °C à sa sortie. Le profil est presque linéaire.
La courbe de distribution de températures de la gaine est donnée sur la

figure 3.7

flux de
FIGURE 3.7 – Distribution de la température de la gaine le long de l’axe y.
On remarque sur la figure (3.7) une augmentation de la température linéaire.

La température à l’entrée est égale à 40,19°C et la température à la sortie
est égale à 51,45 °C.
La courbe de distribution de températures de la chair est donnée sur

la figure (3.8).
FIGURE 3.8 – Distribution de la température de la chair le long de l’axe y.

flux de
On évalue la température de la chair à partir de y égale à 0,05m où on

trouve une température de 50,22 °C puis à y égale à 0,57 m la température
augmente jusqu’à 57,9 °C.
3.10 Résultats obtenus pour l’écoulement transitoire
3.10.1 Évolution en fonction du temps
On prend une période de 4 secondes de simulation avec un pas de temps

de 0,1 secondes. Ce qui importe le plus dans cette simulation c’est de
connaître l’évolution de la température au niveau de la chair , de la gaine
et du fluide. Puisqu’il s’agit d’un transitoire de 3 secondes, on tracera les
profils de température à 1,2 et 3 secondes.
La figure (3.9) donne le profil de la température de la chair, de la gaine et

du fluide à t= 1 seconde pour un canal équivalent à la longueur de la
chair.
FIGURE 3.9 – Évolution de la température le long du canal à t=1s.
On remarque sur la figure (3.9) une augmentation pratiquement linéaire

de la température de fluide, de la gaine et de la chair car la source de
chaleur est une source électrique qui fournie une puissance thermique
linéaire en fonction de la température. En effet à y=0,013m la
température du fluide est égale à 38,9 °C puis elle augmente jusqu’à la
valeur de 54,26°C à y=0,6m. La température de la gaine augmente de

flux de
81,07°C jusqu’à 90,21°C et la

flux de
température de la chair augmente de 82,65°C à 90,69°C. On remarque que la

différence de température à l’entrée et la sortie est de 8,04°C pour la chaire
et 9,14°C pour la gaine alors qu’elle est de 15,36°C pour le fluide.
La figure (3.10) donne le profil de la température de la chair, la gaine

et le fluide à t= 2 seconde.
De même pour t=2s l’évolution est linéaire. Cependant la température

devient plus élevée. A y=0,013m la température de la chair est égale à
87,04°C, celle de la gaine est égale à 82,45°C et celle du fluide 41,49°C. A
y=0,6 m, la température de la chair passe à 104,34°C, celle de la gaine
103,86°C et celle du fluide 61,92°C. On remarque que la différence de
température à l’entrée et la sortie est de 17,3°C pour la chaire et 21,41°C
pour la gaine alors qu’elle est de 20,43°C pour le fluide.
La figure (3.11) donne le profil de température de la chair, de la gaine

et du fluide à t= 3 seconde.

flux de
A t=3s. La température s’élève d’autant plus. A y=0,013m la température

de la chair est égale à 91,14°C, celle de la gaine est égale à 87,05°C et celle
du fluide elle est égale à 44,89°C. A y=0,6 m, la température de la chair est
égale à 109,14°C, celle de la gaine est égale à 108,44°C et celle du fluide
elle est égale à 64,5°C. On remarque que la différence de température à
l’entrée et la sortie est de 17,3°C pour la chaire et 21,41°C pour la gaine
alors qu’elle est de 19,61°C pour le fluide. On remarque que la différence de
température à l’entrée et la sortie est de 18°C pour la chaire et 21,39°C
pour la gaine alors qu’elle est de 19,61°C pour le fluide.
3.10.2 Évolution temporelle de la température maximale
Après avoir tracé les courbes de température sur toute la longueur de la

chair pour trois pas de temps différents, on trace pour une position
déterminée de la longueur de la chair l’évolution de la température en
fonction du temps. La figure (3.12) illustre le profil de température de la
chair, la gaine et le fluide de refroidissement à la position y= 0,6 m en
fonction du temps.

flux de
FIGURE 3.12 – Évolution de la température à y=0,6 m en fonction du temps.
On voit sur la figure (3.12) trois zones différentes. La première zone de zéro
à 3s est la partie de la simulation de l’écoulement transitoire où la vitesse
diminue de 20 % de sa valeur initiale qui est de 0,94 m/s. La deuxième partie
est la partie où il y a rétablissement du débit initial. Enfin, la troisième
partie "Évolution stationnaire" où le débit redevient constant.
Pour la zone "Évolution instationnaire", la température enregistre une

augmentation significative en trois secondes. Tel que, pour la chair la
température s’élève de 38°C à 109,15 °C. Pour la gaine, la température passe
de 38°C à 108,43°C et enfin la température du fluide de refroidissement
augmente de 38°C à 64,5 °C.
Après rétablissement du débit, L’évolution de la température attendue

est représentée dans la zone "Rétablissement du débit" où les trois courbes

flux de
de température diminuent avec une tendance linéaire jusqu’à ce qu’elles

atteignent une température qui correspond au régime stationnaire de
l’écoulement.
3.11 Conclusion
Dans ce travail, on a étudié une simulation expérimentale en utilisant
une expérimentation numérique. L’approche utilisée a permis de déterminer
l’évolution de la température du fluide pour une période de 3 secondes. La
température du fluide a atteint la valeur de 64.5°C. A 64.5°C il n’y a pas
d’ONB (Onset Of Nucleate Boiling) initiation de l’ébullition nuclée. Ceci
est suffisant pour l’étude de sûreté car la pompe qui assure la redondance
de la pompe défaillante met environ 3 secondes pour démarrer (Soit
automatiquement sans Scram soit avec Scram) donc la sûreté de
l’installation est pleinement assurée. Comme l’expérimentation 2D a donnée
de bons résultats avec le code COMSOL en 2D, on va essayer de les
confirmer avec un autre code de calcul et en 3 dimensions de l’espace.

CHAPITRE 4
AnaLYse THERMohYDRaULIqUE, en 3D
à sTaTIonnaIR ET TRansIToIRE,
L’éTaT E
d’Un canaL DU coeUR D’Un RÉacTEUR
nUcLÉaIRE DE REcheRche
Chapitre 4. Analyse thermohydraulique, en 3D à l’état stationnaire
et transitoire, d’un canal du coeur d’un réacteur nucléaire de
4.1 Introduction
Dans ce chapitre, on va se pencher sur la cas tridimensionnel de
l’écoulement du fluide caloporteur dans un canal du réacteur nucléaire de
recherche. La simulation se fera sous le code de calcul FLUENT. Les mises
en équations qu’on établira concernent l’analyse thermohydraulique en
mode stationnaire et en mode transitoire.
Le but principal recherché dans le présent travail est de suivre l’évolution de

la température des composants du coeur d’un réacteur nucléaire de recherche
lors d’un fonctionnement normal ainsi qu’en situation accidentelle de perte
de débit afin de vérifier sa sûreté nucléaire.
4.2 Présentation du Code FLUENT

Le code FLUENT est un code de calcul conçu pour la simulation
des écoulements de fluide et du transfert de chaleur dans des géométries
complexes, il présente une grande flexibilité des résultats et une adaptation
avec n’importe quel type de maillage, il utilise des maillages non-structurés
2D (à maille triangulaire / quadrilatérale) ou 3D (à maille tétraédrique /
hexaédrique / pyramide).
Ce produit inclut également un véritable environnement CAO et un

mailleur paramétrique de dernière génération, permettant ainsi de mettre
en place, rapidement les modèles numériques ou de s’intégrer aux outils de
conception déjà existants comme SOLIDWORKS, CATIA, I-DEAS etc.
Ce code utilise l’approche des volumes finis pour discrétiser et résoudre les
équations de la dynamique des fluides (équation de conservation de la
masse, de la quantité de mouvement et de l’énergie) et les équations de
modélisation de la turbulence. Il est écrit en langage C, et emploie toute
la flexibilité et la puissance qu’offre ce langage.

4.3 Mise en équation thermohydraulique du canal en

3D
Avant de poser les équations qui décrivent l’écoulement dans le canal il
est fondamental d’introduire sa géométrie qui vient sous la forme de deux
plaques de combustible parallèles comme cela a été mentionné dans les
chapitres précédents. L’écoulement s’effectue verticalement de haut en bas
suivant l’axe x comme cela est montré sur la figure (4.1) :
FIGURE 4.1 – schéma du canal thermohydraulique en 3D.
4.3.1 Les équation de conservation
Équation de conservation de la masse
En développant le terme de divergence en tridimensionnel dans l’équation

(3.1) du chapitre précédent elle devient :
∂ρ ∂ρ ∂ρ ∂ρ ∂u ∂v ∂ω
+u +v +ω +ρ + + =0 (4.1)
∂ ∂ ∂ ∂ ∂ ∂ ∂
On considère que l’écoulement est incompressible, ρ = cste
On trouve de ce fait :
∂ρ ∂ρ ∂ρ ∂ρ
= = = =0 (4.2)
∂t ∂x ∂y ∂z

Par conséquent l’équation de conservation de la masse s’écrit ainsi :

∂u ∂v ∂ω
+ + =0 (4.3)
∂x ∂y ∂z
Équation de conservation de la quantité de mouvement
En développant le terme de divergence en tridimensionnel de l’équation

(3.2) du chapitre précèdent elle devient :
∂ ∂U ∂ ∂ 1
→ = gr→adP +
U U→ U→ ν∆U→ ρ
+ ρ→g (4.4)
→ +u +v +ω
∂x ∂y ∂z
∂t −
On procède à une projection de l’équation (4.4) dans le chapitre précédent
sur les deux axe x, y et z :
L’écoulement se fait selon la direction x donc :

ρgy = ρgz = 0.
La projection selon x :
∂u ∂u ∂u ∂u 1 ∂p ∂2u ∂2u ∂2u
+u +v +w =− + + + + ρgx (4.5)
∂t ∂x ∂y ∂z ρ ∂x ∂x2 ∂y2 ∂z2
ν
La projection selon y :
∂v ∂v ∂v ∂v
1 ∂p 2 ∂2v ∂2v
∂v
+u +v +w =− + + + + ρgy (4.6)
∂t ∂x ∂y ∂z ρ ∂y ∂x2 ∂y2 ∂z2
ν
La projection selon z :
∂w ∂w ∂w ∂w 1 ∂p ∂ 2w ∂2w ∂2w
+u +v +w =− + + ρgz
∂t ∂x ∂y ∂z ρ ∂z ∂x2 ∂y2 ∂
Où u est la vitesse du caloporteur selon x, v est sa vitesse selon y et w sa
vitesse selon z.
En moyennant les équations de conservation de mouvement par rapport au

temps on obtient l’équation précédemment cité (3.5) :
∂ ∂ h i
(ρu )+ ∂p ∂ ∂ui ∂uj 2 ∂ul ∂ ′ ′
i
(ρu iu j) = − + µ + − δ + (−ρu u )
∂ ∂ ∂

Équation de conservation de l’énergie
L’équation de conservation de l’énergie a été déjà mentionnée dans le

chapitre précédent On va réécrire l’équation (3.3) :
∂ρE h i h Σ i
∂t + ∇. (ρE + p)U→ = ∇. (k∇T ) − J
+ (τ¯.U→ ) + S
→j
hj
j
4.3.2 Le modèle de turbulence
Le modèle de turbulence choisit pour l’analyse thermohydraulique

tridimensionnelle est le même que celui qu’on a choisit dans la simulation
bidimentionnel du chapitre précédant. Il s’agit du Modèle K− ω SST
(Shear-Stress Transport) (Menter, 1994) car il correspond parfaitement à
l’écoulement turbulent établi dans notre cas et réunit les meilleures
propriétés
— de K ε (valeurs indépendantes en champ lointain) etK ω de Wilcox
(1993) (robustesse, comportement asymptotique correct près des murs).
4.4 Conditions aux limites

Le modèle des deux couches est utilisé pour le calcul numérique de
l’écoulement turbulent dans la région située près de la paroi. Dans cette
approche, l’ensemble du domaine est subdivisé en une région affectée par la
viscosité et une région entièrement turbulente (Sidi Ali ; Gatignol, 2010).
La démarcation des deux régions est déterminée par un nombre de

Reynolds tur√bulent, Rey , basé sur la distance à la paroi dont
l’expression
est Rey = ρy k/µ où y est la distance du point considéré par rapport à la
paroi, tel que présenté par (Wilcox, 1993).
La première grille numérique près de la paroi est située suffisamment loin

de d’elle pour que le nombre de Reynolds turbulent soit supérieur à l’unité,
ce qui signifie que les forces de turbulence sont plus importantes que les
forces visqueuses. Les deux expressions de la vitesse adimensionnelle
donnée dans l’équation (4.8), et utilisées dans ce travail, ont été introduites
par Prandtl (1927, 1932) et sont des fonctions de la paroi :

et transitoire, d’un canal
(du coeur d’un réacteur nucléaire de
U y+ y+ ≤ y+
= 1 +
ν (4.8)
ln(E ≥
U k
+ y+
)
ν

Où, la vitesse de friction locale U ∗ est calculée grâce à U ∗ = 1/2

C1/4 µ k et
où y+ = ρU ∗y/µ est la position sans dimension de y, et y+ est l’épaisseur
ν
sans dimension de la sous-couche visqueuse égale à 11,225. La constante E
est égale à 9,793 et la constante de Von Karman kvk est égale à 0,4187. Les
flux de diffusion pour toutes les quantités dans la direction de la sortie sont
égaux à zéro comme proposé par Chung (2002) :
∂Ui ∂k ∂ε
=0; = 0 et =0 (4.9)
∂n ∂n ∂n
avec (∂/∂n) la dérivée par rapport à la direction normale à la la surface de
sortie.
La vitesse près de la paroi est supposée nulle :
Upr = 0 (4.10)
À la paroi, le modèle des deux couches est utilisé pour le calcul numérique
de l’écoulement turbulent dans les régions situées à proximité de la paroi on
utilise l’expression de la vitesse donnée dans (4.8).
4.5 Maillage du domaine numérique à étudier

La géométrie tridimensionnelle a été modélisée par le logiciel Solid
Works. Elle comporte la chair, la gaine et le caloporteur. Sa création a été
réalisée de sorte à ce que les éléments du maillage soient de la forme d’un
hexaèdre régulier. Pour se faire, la gaine a été découpée en 17
parallélépipèdes. Sur la figure (4.2), le découpage de la géométrie selon les
trois plans XY, YZ, XZ sur Ansys.

FIGURE 4.2 – Maillage du domaine numérique
FIGURE 4.3 – Géométrie du canal 3D.
L’importance du choix du maillage n’est pas négligeable. En effet, il affecte

grandement la cohérence des résultats obtenues par la suite. Dans ce travail,
le maillage a été effectué via Ansys workbench. Grâce au découpage de la
géométrie montrée sur la figure (4.2), on a pu élaborer un maillage conforme
et uniforme. Lorsqu’il s’agit d’un solide comme pour le cas de la chair et de
la gaine, le nombre des noeuds et la taille des éléments ne pose généralement
pas de problème quant aux résultats de simulation. Néanmoins, le nombre
choisi doit être suffisamment grand pour obtenir des résultats en tout point
du solide.
Pour un fluide en écoulement turbulent, comme dans notre cas, on

doit tenir compte de la couche limite en créant de petites subdivisions dans
le maillage à son niveau car à un nombre de Reynolds élevé, la couche

limite est si mince qu’il est difficile d’utiliser suffisamment de points de

grille pour la résoudre. Pour y parvenir, il suffit de calculer son épaisseur
ainsi qu’un paramètre qu’on appelle fonction de paroi Y Plus (Y +). On
prend Y Plus (Y +) égal à 20. Sur le tableau (4.1) les paramètres entrant dans
le calcul de la couche limite sont présentés.
Tableau 4.1 – Paramètre pour le calcul de l’épaisseur de la couche limite
Paramètre Valeurs
Vitesse d’écoulement 0, 94 m/s
Masse volumique de l’eau 993 kg/m3
Viscosité dynamique 6, 78 . 10−4 kg/m.s
Diamètre hydraulique 0, 0058 m
YPlus (Y+) 20
Épaisseur de la première cellule prés de la paroi 2, 5 . 10−5 m
Puisque l’épaisseur de la cellule près de la paroi est égale à 2, 5.10−5m et que

le facteur d’augmentation est de 1, 2. Les cellules loin de la paroi et celle de
la gaine et de la chair sont de l’ordre de 4 . 10−4 m. Le nombre de mailles
total est de 1273927 mailles. Sur la figure (4.4) est donnée le maillage de la
géométrie selon les trois plans.
FIGURE 4.4 – Maillage du canal selon les trois plans de l’espace.
Pour savoir si notre maillage est correct, il faut que les résultats obtenus à
savoir la température soient indépendants du nombre de mailles. On
procède, par conséquent, à une évaluation de l’influence de la résolution du
maillage sur la température au centre de la gaine. A cet effet, plusieurs
simulations ont été réalisés pour l’étude de ce paramètre avec différentes
résolutions de maillage. Sur la figure (4.5) est présenté la courbe donnant la
température trouvée au centre de la gaine en fonction du nombre de
noeuds.
FIGURE 4.5 – Température en fonction du nombre de mailles.
On remarque sur le graphe de la figure (4.5) que pour un même point on

obtient différentes valeurs de température quand le nombre de mailles est
inférieur à 943121 mailles. On en déduit donc que le résultat dépend du
nombre de mailles.
Lorsqu’on dépasse 943121 mailles, le graphe se stabilise et la température

obtenue au même point devient constante et égale à 79,75 °C quelque soit le
maillage. On en conclut que le maillage choisi de 1273927 mailles convient
à notre géométrie. La sensibilité au maillage est ainsi vérifiée.
4.6 Méthode des volumes finis

La résolution se fait avec la méthode des volumes finis qui est une
méthode de calcul numérique utilisée pour résoudre les équations aux
dérivées partielles qu’elles soient de forme elliptiques, paraboliques ou bien
hyperboliques. Cette méthode est adaptée aux équations de la conservation
dans leur forme intégrale et concerne un large domaine d’application comme
par exemple la mécanique des fluides, transfert de masse et transfert de
chaleur. Ce qui la caractérise des autres méthodes de résolution c’est sa
capacité à conserver localement une quantité d’une cellule à une autre. Cette
caractéristique est due au fait qu’elle soit basée sur une approche de balance.
Ainsi une balance locale, aussi appelé volume de contrôle, est présente dans
chaque cellule de discrétisation ce qui permets la modélisation des
problèmes ou le flux est très important.

La résolution des équations de l’écoulement est faite par la méthode des

volumes finis. La discrétisation des équations est classique et peut être
illustrée par l’équation de conservation à l’équilibre pour le transport d’une
quantité scalaire φ. Cette équation, sous forme intégrale pour un volume de
contrôle arbitraire V , est donnée par l’équation (4.11)
I I ∫
ρ φ U→ d→A = Γφ ∇φ . d→A + Sφ dV (4.11)

v
Où, A→ est le vecteur de surface, Γφ le coefficient de diffusion,∇φ est le gradient
de φ et Sφ la source par unité de volume.
L’équation (4.11) est appliquée à chaque volume de contrôle, ou cellule dans
le domaine de calcul et est discrétisée. On donne l’équation (4.12) :
Nfaces
Nfaces
Σ ρ φ → → Σ
f f U f . A f Γφf ∇φf . A→f + Sφ V (4.12)
f
= f
Où, N[faces est le nombre de faces entourant la cellule, φf la valeur de φ

convexe à travers la surface f, ρf φf U→f . A→f le flux de φ à travers la face f,
A→f le vecteur de la surface de la face f,∇ φf le gradient de φ à la face f et
V le volume de la cellule. L’équation (4.12) conduit à une équation
discrétisée de la forme (Sidi Ali ; Gatignol, 2010) :
ap φ f = Σ
af φf + S (4.13)
f
La figure (4.6) montre une cellule discrétisé par la méthode des volumes
finis :
FIGURE 4.6 – Cellule discrétisée à deux dimensions.
4.7 Analyse thermohydraulique à l’état stationnaire

Dans cette analyse thermohydraulique, on effectue plusieurs simulations

stationnaire. La première reprend le fonctionnement normal d’un réacteur

nucléaire de recherche de 2MW. Puis on passe à des débits réduit de 5%,

10%, 15% et 20% toujours à l’état stationnaire. Les conditions aux limites
du domaine concernent la vitesse et la température à l’entrée du canal, la
pression à sa sortie et la condition limite aux parois de contacte entre le
fluide caloporteur et la gaine.
a) A l’entrée du canal
La vitesse du caloporteur à l’entrée est constante et est égale à 0,94 m/s

en fonctionnement normal. Elle est égale à 0,893 m/s lorsque le débit est
réduit à 5%, à 0,846 m/s lorsqu’il est réduit de 10 %, à 0,799 m/s lorsqu’il
est réduit à 15 % puis à 0,752 % lorsque le débit est réduit à 20 %.
Pour ce qui est de la température du fluide caloporteur, sa valeur est de
38°C à l’entrée du canal.
b) A la sortie du canal
La pression de sortie est égale à zéro et diffère de la pression à l’entrée afin

qu’un écoulement puisse s’établir.
c) Condition aux paroi
La vitesse du fluide caloporteur aux parois est nulle.
d) Source de chaleur
Puisqu’il y a transfère de chaleur entre les différents composants du canal

(chair, gaine et caloporteur), on désigne la chair comme étant la source de
chaleur et on impose un flux sinusoïdale dont l’équation (2.1) est donnée :
πx
q′(z) = q′ cos
c
Le
Cette équation est introduite à l’aide d’une fonction UDF (User Define
Function) dont le script est présenté sur la figure (4.7).

FIGURE 4.7 – Fonction UDF du flux de chaleur.
4.8 Résultats et interprétations de la simulation

stationnaire
Les profils des températures obtenus sont données dans ci-après.
4.8.1 Pour une vitesse constante
Les résultats obtenus pour une vitesse de 0,94 m/s sont traduits sous
forme de graphes et de contours. La variable dont on désire connaître
l’évolution est la température. Donc pour mieux interpréter les résultats,
on établit un code Matlab qui sert à moyenner les températures à chaque
pas de x et par conséquent avoir un seul profil de température le long du
canal. Pour le fluide caloporteur, on prend les données au niveau du plan
de symétrie et pour la gaine, on prend la surface qui est en contacte avec le
fluide caloporteur.
Sur la figure (4.8), (1) et (2) représentent respectivement la distribution de

la température de la gaine et du fluide caloporteur dans le plan xz.

FIGURE 4.8 – Contour de la température du caloporteur et de la gaine sur le plan xz en

stationnaire
Sur la figure (4.9), le profil de température de la gaine et du caloporteur en

fonction de la longueur adimensionnelle du canal.
FIGURE 4.9 – Profils des températures en fonction de la longueur adimensionnelle.
Les courbes de la figure (4.9), donne la distribution de la température du

caloporteur et de la gaine du début de la chair à 0,0209 jusqu’à la fin de la
chair à 0,97.

La température du caloporteur à 0,0209 est égale à 38°C et augmente jusqu’à

52,30 °C à 0,97. On peut voir sur la figure (4.8) (2) la distribution de la
température du caloporteur dans le plan xz.
La température de la gaine est de 57,08 °C à 0,0209 et augmente jusqu’à

un pic de 79,9 °C à 0,45 puis diminue et devient égale à la température du
caloporteur de 52,30 °C à 0,97.
4.8.2 Analyse d’un LOFA constant en mode stationnaire pour

plusieurs pertes de débits
La diminution du débit de façon progressive de 5% à 20 % en stationnaire

est étudié. On a calculé en stationnaire pour chaque débit les différentes
évolution de température qui nous intéressent.
La figure (4.10) représente les Contours des températures du caloporteur et

de la gaine sur le plan pour une diminution de débit de 5%, 10 %, 15% et
20%.
FIGURE 4.10 – Contour des températures de la gaine et du caloporteur sur le plan xz

pour 5%, 10 %, 15% et 20%.
(a1) : Contour de température de la gaine à 5 % de réduction de débit

(a2) : Contour de température du caloporteur à 5 % de réduction de débit
(b1) : Contour de température de la gaine à 10 % de réduction de débit

(b2) : Contour de température du caloporteur à 10 % de réduction de débit

(c1) : Contour de température de la gaine à 15 % de réduction de débit
(c2) : Contour de température du caloporteur à 15 % de réduction de débit
(d1) : Contour de température de la gaine à 20 % de réduction de débit
(d2) : Contour de température du caloporteur à 20 % de réduction de débit
Les températures en fonction de la longueur adimensionnelle du canal sont
présentées sur la figure (4.11).
FIGURE 4.11 – Profil de la température pour plusieurs débits en fonction de la longueur

adimensionnelle.
La figure (4.11) englobe les profils des températures du caloporteur et de la

gaine de 5% , 10% , 15% et 20 % de perte de débit pour la partie en contact
direct entre la chair, la gaine et le caloporteur. En d’autre termes, pour une
longueur adimensionnelle du canal allant de 0,0209 à 0,97.
Pour le caloporteur à 0,0209, la température à l’entrée vaut 38°C puis elle

augmente et diffère à 0,97 de la longueur adimensionnelle du canal. En effet,
la température du caloporteur à 5 % de perte de débit est égale à 53,05 °C.
Celle avec une perte de 10 % est égale à 53,87 °C. Celle avec une perte de
15 % est égale à 54,79°C et celle avec une perte de 20% égale à 55,83°C. On

remarque une augmentation progressive de la température à chaque fois que

le débit diminue et cela est d’autant plus remarquable à la sortie du fluide
caloporteur.
La gaine à 0,0209 a une température de 57,76°C lorsque la perte de débit

est de 5 %. Elle est de 58,48°C à 10 % de perte de débit, 59,28°C à 15 % de
perte de débit et 60,15°C à 20 % de perte de débit.
La gaine atteint une température maximale à une distance adimensionnelle

de 0,45 quelque soit le débit. Lorsque la perte de débit est de 5 % elle est
égale à de 81,75°C . Elle atteint les 83,88°C lorsque la perte de débit est de
10
%. Puis à 15 % de perte la température maximale de la gaine devient égale
à 86,25°C et à 20 % de perte elle devient égale à 88,89 °C.
A une distance adimensionnelle de 0,97 la température de la gaine devient
égale à la température du caloporteur dont les valeurs ont été données
précédemment.
FIGURE 4.12 – Évolution au point le plus chaud de la gaine de la température en fonction

du pourcentage de perte de débit.

La figure (4.12) donne un aperçu de l’augmentation de la température

lorsqu’on diminue le débit de 20 %. Tel que pour la gaine au pic
de température c’est à dire à une distance adimensionnelle de 0,45 la
température évolue linéairement avec une pente de 0,45 et va de 79,91°C
à 88,89°C. La température de la gaine à une distance adimensionnelle de
0,0209 va de 57,08°C à 60,15 % de façon linéaire avec une pente de 0,15.
Pour le caloporteur à une distance adimensionnelle de 0,97 l’augmentation
est aussi à tendance linéaire avec une pente de 0,18 et va de 52,33°C à 55,87°C.
Ce qu’on peut en déduire c’est que l’augmentation de température la

plus rapide est celle de la gaine puisqu’elle se trouve en contact direct avec
la source puis vient celle du fluide caloporteur étant donné qu’il est en
mouvement et qu’il absorbe la chaleur transmise par la gaine.
La simulation d’un LOFA stationnaire a donné des résultats où on voit

une augmentation de température à chaque fois que le débit du caloporteur
est réduit ce qui semble corrects. Cependant le but de la simulation d’un
accident LOFA est de voir à quel moment la diminution du débit va affecter
la sûreté nucléaire du coeur du réacteur pour pouvoir dimensionner un
système de protection adapté. Une simulation transitoire est donc plus
appropriée puisque le temps n’est pas négligé.
4.9 Analyse thermohydraulique à l’état transitoire

Dans cette analyse thermohydraulique, on établit la simulation transitoire
de l’écoulement dans le cas d’un fonctionnement accidentel dans un réacteur
nucléaire de recherche de 2MW. Dans cette dernière on va poser les
même conditions aux limites de l’état stationnaire et on diminue le débit
du caloporteur en diminuant sa vitesse, la section de passage du fluide
caloporteur étant constante.
La vitesse est d’abord égale à 0,94 m/s pour un intervalle de temps

de trois secondes puis diminue de 20 % pendant les deux prochaines
secondes puis le calcul transitoire débute. Cette variation est implémentée à
l’aide d’une fonction UDF qu’on donne sur la figure (4.13).

FIGURE 4.13 – Fonction UDF de la variation de vitesse.
On essayera dans cette simulation d’évaluer à quel moment l’ébullition du

caloporteur commence et donc le temps nécessaire pour mettre en marche la
pompe redondante.
4.10 Résultats et interprétations de la simulation

transitoire
Pour la simulation transitoire l’évolution des températures en fonction du
temps est données sur les figures suivantes. Tel que le transitoire commence
à partir de 6 secondes.
La figure (4.14) représente les Contours des températures du caloporteur sur
le plan à chaque pas de temps de 6s à 10s.

FIGURE 4.14 – Contour des températures du caloporteur sur le plan xz à chaque pas de
temps de 6s à 10s.
La figure (4.15) donne la distribution de la température du caloporteur le

long du canal à différents pas de temps :
FIGURE 4.15 – Profil de la température en fonction du temps.

On remarque sur la figure (4.15) qu’en chaque point du canal la température

du caloporteur augmente en fonction du temps. A 0,0209 de distance
adimensionnelle la température à 6s est égale à 38°C, à 7s elle est égale à
38,98°C, à 8s elle est égale à 41,71°C, à 9s elle est égale à 46,34°C et à 10s
elle est égale à 52,16°C.
Concernant le pic de température, il est atteint pour chaque pas de temps en

différents points. A 6s la température maximale est de 54,35°C à une
distance adimentionnelle de 0,97. A 7s la température maximale est de
54,68°C à une distance adimentionnelle de0,87. A 8s la température
maximale est de 56,54°C à une distance adimentionnelle de 0,73. A 9s la
température maximale est de 62,28°C à une distance adimentionnelle de
0,58 et enfin à 10s la température maximale est de 71,05°C à une distance
adimentionnelle de 0,5.
A 0,97 de distance adimensionnelle toutes les courbes de température

s’approche en diminuant vers une valeur de 54°C.
On observe un changement de la forme de la courbe de distribution de la

température du caloporteur de 0,0209 à 0,97 de la distance adimensionnelle
en fonction du temps. Tel que plus le temps augmente de 6 à 10 secondes
plus la courbe épouse la forme d’une cloche.
La figure (4.16) représente les Contours des températures du caloporteur sur

le plan à chaque pas de temps de 6s à 10s.

FIGURE 4.16 – Contour des températures de la gaine sur le plan xz à chaque pas de temps
de 6s à 10s.
La figure (4.17) donne la distribution de la température de la gaine en

fonction de la longueur adimensionnelle du canal.

FIGURE 4.17 – Profil de température en fonction du temps.
On remarque sur la figure (4.17) que plus on augmente le temps, plus

la température de la gaine augmente. A une distance adimentionnelle de
0,0209, pour un temps de 6s la température de la gaine est égale à 59,86°C.
Pour un temps de 7s, elle est égale à 71,62°C. Pour un temps de 8s, elle est
égale à 83,34°C. Pour un temps de 9s, elle est égale à 94,23°C et puis pour
un temps de 10s elle est égale à 104,46°C.
Le pic de température à 6s est égal à 86,028°C à une distance

adimentionnelle de 0,43. A 7s il est égal à 103,38°C à une distance
adimentionnelle de 0,41. A 8s il est égal à 122,32°C à une distance
adimentionnelle de 0,39. A 9s il est égale à 140,62°C à une distance
adimentionnelle de 0,389. 10s il est égal à 158,27°C à une distance
adimentionnelle de 0,379. Les courbe de température se rejoignent à une
distance adimentionnelle de 0,97 avec une valeur qui s’approche de 53°C.
On peut observer que plus le temps avance plus la pente de la courbe
parabolique augmente. Cela est du au fait que le flux de chaleur généré est
sinusoïdale.
La figure (4.18) représente l’évolution de la température en fonction du temps

de 6s à 10s.
FIGURE 4.18 – Profil de température en fonction du temps.
Les températures ont été prise à un point donné de la longueur

adimensionnelle du canal pour évaluer le taux d’augmentation de la
température en fonction du temps. On choisit de prendre le point où le pic
est atteint à 10s c’est à dire à une distance adimentionnelle de 0,5 pour le
caloporteur et à une distance adimentionnelle de 0,38 pour la gaine.
On remarque d’abord qu’en stationnaire de 0 à 3s la température ne

varie pas elle est égale à 48°C pour le caloporteur et 75°C pour la gaine.
Après avoir diminuer le débit de 20 % de 3 à 6s il y a une légère
augmentation de température car il y a eu une transition d’un état
stationnaire à une autre état stationnaire. C’est à partir de la 6 ème
seconde que le transitoire se lance. On observe alors que l’augmentation de
la température de la gaine est plus rapide que celle du fluide avec une
tendance linéaire. Donc en 4 secondes de transitoire le caloporteur
enregistre une différence de 23,84°C de 47,21°C à 71,05°C alors que la
gaine augmente de 72,6°C de 85,65°C à 158,26°C.
On remarque que le caloporteur n’atteint pas la barre des 100°C après

quatre secondes de transitoire. Suite à ça, on a décidé de lancer une
deuxième simulation avec un temps de calcul plus long.

La figure (4.19) donne les contours des températures du caloporteur lors du

transitoire de la deuxième simulation de 6s à 15s.
FIGURE 4.19 – Contour des températures du caloporteur sur le plan xz à chaque pas de
temps 6s à 15s.
La figure (4.19) montre les contours des températures de la gaine lors du

transitoire de la deuxième simulation de 6s à 15s.

FIGURE 4.20 – Contour des températures de la gaine sur le plan xz de contact avec le
caloporteur à chaque pas de temps 6s à 15s.
La figure (4.21) donne l’évolution de la température en fonction du temps

avant et après le lancement du transitoire qui commence à la 6ème secondes.
FIGURE 4.21 – Profil de température en fonction du temps
De 0 à 3 secondes la simulation est à l’état stationnaire et la température de

la gaine à une distance de 0,38 et du caloporteur à une distance
adimensionnelle de 0,5 sont constantes tel que le caloporteur a une
température de 48°C et
la gaine a une valeur de 75°C.
De 3 à 6 secondes, on fixe la vitesse à −20% de la vitesse initiale.

Dans ce cas là, on enregistre une légère augmentation.
Le lancement du transitoire se fait alors à partir de 6 secondes. On

remarque pour cela une augmentation linéaire très rapide de la température
de la gaine et du caloporteur. A 15 secondes de la simulation, la température
de la gaine est égale à 240,30 °C et celle du caloporteur est égale à 130 °C.
Les résultats de la deuxième simulations ont permit d’observer qu’à

12,5 secondes les 100°C sont dépassés pour le caloporteur. Le phénomène
d’ONB ne se manifeste donc qu’à 12,5 secondes. L’apparition des bulles
d’air peut détériorer le transfert de chaleur entre la gaine et le fluide et peut
ainsi accélérer la fusion de cette gaine. Il est donc impératif qu’une pompe
redondante comblant le débit manquant puisse être démarré dans les 10
premières secondes du début de l’accident.
4.11 Conclusion
Dans ce travail, on a étudié des simulations stationnaires et transitoires
tridimensionnelle de l’écoulement à travers un canal du coeur d’un réacteur
nucléaire de recherche. En stationnaire, on a reproduit l’écoulement du
caloporteur lors d’un fonctionnement normal du réacteur. Les résultats
obtenus sont satisfaisants. Par la suite, on a lancé des simulations
stationnaires où le débit, a été diminué à chaque fois de 5, 10, 15 et
20%, le but étant de faire un LOFA. Après l’analyse de l’évolution des
profils de température, on a constaté que malgré la diminution du débit, la
température d’ébullition n’a pas été atteinte.
En transitoire, on a simulé un accident de LOFA de 4s où une diminution de

débit de 20 % a été imposée. Les résultats de cette simulation ont démontré
qu’après 4 s de transitoire il n y a pas de phénomène d’ONB. Par conséquent
le temps a été rallongé de 9s. Cette intervalle de temps a été suffisant pour
observer le moment pour lequel il y a apparition d’une ébullition puisque
le caloporteur arrive à une température de 100°C après 6,5 secondes de
transitoire.
En bref, l’analyse thermohydraulique de sûreté nucléaire de la première

barrière de sûreté a été vérifiée et les limites à ne jamais franchir ont été
bien définies.

ConcLUsIon généRaLE
Dan le présent travail, plusieurs outils pour l’étude thermohydraulique

pour vérifier la sûreté nucléaire d’un canal du coeur d’un réacteur nucléaire
de recherche ont été utilisés en 1D, 2D et 3D.
Dans le premier chapitre, la description d’un réacteur nucléaire de recherche

type MTR a été faite. Les composants des circuits primaire et secondaire
ont été expliqués. Le processus de fission et l’ensemble de la technologie
du combustible sont décrits avec force détails. Une partie a été consacrée
aux différents types d’accidents nucléaires qui peuvent avoir lieu dans un
réacteur nucléaire. Ceci a été fait dans le but de bien présenter un des trois
accidents qui va faire l’objet d’une partie de ce travail. Il est très important
de noter que le refroidissement du coeur d’un réacteur nucléaire doit être
refroidi quelque soit la situation, soit en fonctionnement normal ou
accidentel.
Dans le second chapitre, une analyse thermohydraulique en utilisant le

modèle d’El-Wakil (1971) et le modèle de Todreas et al. (1989) a été faite.
Les résultats obtenus ont été vérifiés et validés en se basant sur des travaux
valorisés. De plus, on a appliqué ces deux modèles sur trois réacteurs
nucléaire de recherche de puissance distincts 1MW, 2MW et 10MW. Par la
suite on a examiné leurs sûreté nucléaire par le calcul du flux de chaleur
moyen et critique et le rapport de flux de chaleur critique et on en a déduit
qu’ils fonctionnent en toute sûreté. Le but de l’analyse de sûreté est de
vérifier et d’éviter qu’il n’y ait un phénomène d’ébullition (ONB). Ce travail
a été effectué au moyen d’un programme Matlab qu’on peut adapter à
n’importe quel réacteur nucléaire de recherche utilisant un combustible de
type plaque. On peut finir par dire que ces méthodes asymptotiques sont
très utiles pour la compréhension de la thermohydraulique des réacteurs
nucléaires de recherche et une bonne base pour aboutir à des résultats

Conclusio
exploitables en des temps limités.
Dans le troisième chapitre, on a étudié une simulation expérimentale en

utilisant une expérimentation numérique. L’approche utilisée a permis de
déterminer l’évolution de la température du fluide pour une période de 3
secondes. La température du fluide a atteint la valeur de 64.5°C. A 64.5°C
il n’y a pas d’ONB (Onset Of Nucleate Boiling) initiation de l’ébullition
nuclée. Ceci est suffisant pour l’étude de sûreté car la pompe qui assure la
redondance de la pompe défaillante met environ 3 secondes pour démarrer
(Soit automatiquement sans Scram soit avec Scram) donc la sûreté de
l’installation est pleinement assurée. Comme l’expérimentation 2D a donnée
de bons résultats avec le code COMSOL en 2D, on va essayer de les confirmer
avec un autre code de calcul et en 3 dimensions de l’espace.
Enfin, on a étudié des simulations stationnaires et transitoires

tridimensionnelle de l’écoulement à travers un canal du coeur d’un réacteur
nucléaire de recherche. En stationnaire, on a reproduit l’écoulement du
caloporteur lors d’un fonctionnement normal du réacteur. Les résultats
obtenus sont satisfaisants. Par la suite, on a lancé d’autre simulations
stationnaire où le débit a été diminué à chaque fois de 5, 10, 15 et 20%.
Le but étant de faire un LOFA. Après analyse de l’évolution des profils de
température, on a constaté que malgré la diminution de débit, la température
d’ébullition n’est pas atteinte.
En transitoire, on a simulé un accident de LOFA de 4s où une diminution de

débit de 20 % a été imposée. Les résultats de cette simulation ont démontré
qu’après 4 s de transitoire il n y a pas de phénomène d’ONB. Par conséquent
le temps a été rallongé de 9s. Cette intervalle de temps a été suffisant pour
observer le moment pour lequel il y a apparition d’une ébullition puisque
le caloporteur arrive à une température de 100°C après 6,5 secondes de
transitoire.
Donc pour assurer la sûreté nucléaire du réacteur, il faudrait dimensionner

une pompe redondante capable de démarrer avant que le LOFA n’atteint les
6,5 secondes si la perte de débit est de 20 %.
En dernier, l’objectif de ce travail a été pleinement atteint. Des programmes

de calcul ont été établis et validés, des simulations en 2D et 3D
d’écoulements en mode stationnaire et transitoire ont été réalisées. Tous les
résultats
Conclusio
obtenus versent dans le respect des critères de sûreté nucléaire qui doivent
être respectés.

BIBLIogRaPHIE
*, 2008. Les combustibles nucléaires. Jean François Parisot. "monographie EDEN".

**, 2017. Research reactors for the development of materials and fuels for innovative
nuclear energy systems. the IAEA Publishing Section. "nuclear energy series".
***, 1980. Research reactor core conversion from the use of highly enriched uranium
to the use of low enriched enriched iranium fuels guidebook. PHYSICS SECTION
INTERNATIONAL ATOMIC ENERGY AGENCY. "IAEA guide book".
BOUDALI, Z ; SALHI, R, 2011. Contribution à l’étude thermique d’un réacteur nucléaire
embarqué et ses systèmes annexes. PFE Report.
CHEN, Y ; ZHAI, J ; WANG, L, 2007. Computer modeling of multiscale fluid flowand heat
and mass transfer in engineered spaces. Chemical Engineering Science 62, p. 3580-3588.
COMSOL, 2020. CFD Module User’s Guide. comsol licence. "Comsol guides".
COMSOL, 2020. Introduction to Comsol Multiphysics. Comsol. "Nuclear Systems I".
DELIMAREIS, P ; LOMBARDICOSTA, A ; PEREIRA, C ; AUXILIADORA, M
;
HERREROMIRÓ, R ; GUMERSINDO, V, 2019. Evaluation of the neutronic feedback
effects in loss of coolant accident simulation of the IPR R1 TRIGA REACTOR.
Computational Thermal Sciences. T. 11, no 3, p. 243-254.
DOVAL, A, 1998. Validation and verification of the MTR-PC thermohydraulique
package. International Meeting on Reduced Enrichment for Research and Test
Reactors Sao Paulo, Brazil. T. 1, no 1, p. 1-9.
GLASSTONE, S ; SASONSKE, A, 1967. Nuclear Reactor Engineering. Van Nostrand
Reinhold Company New York.
GONNIER, C, 2008. Étudier les situations accidentelles. le moniteur. "CEA".
HINZE, Je ; HINZE, JO, 1975. Turbulence (McGraw-Hill Mechanical Engineering).
McGraw-Hill Companies. "COMPUTATIONAL MECHANICS and APPLIED
ANALYSIS".
HOLMAN, P, 1989. Heat transfer. Mc Graw-Hill Book Company New York.
KHEDR, A, 2008. Thermal-hydraulic fortran program for steady state calculation of plate
type fuel research reactors. research gate. T. 1, no 23, p. 1-12.
LABANI, O ; SIDI ALI, K ; SAIM, B ; OUKIL, K, 2014. Thermal-Hydraulic Analysis
of a Nuclear Research Reactor Core Channel. Journal of Science Research. T. 27, no
7, p. 34-42.
Bibliographi
LAUNDER, BE ; SPALDING, DB, 1974. The numerical computation of turbulent flows.

Computer Methods in Applied Mechanics and Engineering, p. 269-289.
LEWZISJR, G ; ZIZZA, M ; DERIENZO, P, 1964. Heat exchangers in nuclear power
plants. Advances in Nuclear Science and Technology. T. 2, no 1, p. 41-106.
LU, Q ; QIU, S ; SU, G, 2009. Development of a thermal–hydraulic analysiscode for
research reactors with plate fuels. Annals of nuclear Energy. T. 36, no 1, p. 433-447.
MEFTAH, B ; ZIDI, T ; BOUSBIA-SALAH, A, 2006. Neutron flux optimization in
irradiation channels at NUR research reactor. Annals of Nuclear Energy, p. 1164-1175.
MENTER, F, 1994. Two Equation Eddy-Viscosity Turbulence Modeling for Engineering
Applications. AIAA Journal, p. 1598-1605.
MICHEL, A ; HUMBERTO GRACHER, R, 2015. Procedures for Manufacturing Nuclear
Research Reactor Fuel Elements. LAMBERT academic publishing. "monographie
DEN". ISBn 978 3 659 47963 2.
QING, L ; SUIZHENG, Q ; SU, G, 2009. Development of a thermal–hydraulic analysis
code for research reactors with plate fuels. Elsevier. T. Annals of nuclear energy, no
36, p. 433-447.
REDDY, JN ; GARTLING, DK, 2010. The finite Element Method in Heat transfer
and Fluid Dynamics. Taylor et Francis. "COMPUTATIONAL MECHANICS and
APPLIED ANALYSIS".
SADREHAGHIGHI, I, 2021. Structure Meshing for CFD. ANNAPOLIS. "CFD Open
Series".
SHIH, T ; LIOU, W ; SHABBIR, A ; YANG, Z ; ZHU, J, 1995. A new k-ε eddy viscosity
model for high reynolds number turbulent flows. Computers and Fluids, p. 227-238.
SIDI ALI, K, 2002. A COMPUTER PROGRAM PREDICTING STEADY -
STATE PERFORMANCE OF A NUCLEAR RESEARCH REACTOR’S COOLING
SYSTEM. 10th International Conference on Nuclear Engineering, p. 1-6.
SIDI ALI, K ; BOUDALI, Z ; SALHI, R, 2012. Thermal-Hydraulic behavior of physical
quatities at critical velocities in a nuclear research reactor core channel using plate
type fuel. Nuclear Technology and Radiation Protection. T. 27, no 3, p. 229-238.
SIDI ALI, K ; DJABER, A ; ELMOUNDIR, M ; BELMRABET, T, 2017. Thermal
Hydraulic Modeling of a Nuclear Reactor Core Channel Using CFD. Springer.
T. Computational Methods and Experimental Testing In Mechanical Engineering, n o
2, p. 9-16.
SIDI ALI, K ; GATIGNOL, R, 2010. Interfacial friction factor determination using CFD
simulations in a horizontal stratified two-phase flow. Chemical Engineering Science,
p. 5160-5169.
SIDI-ALI, T.Rahem A.Aissati K., 2018. Comparative study between three types of
cooling towers : Dry, wet and hybrid. ResearchGate. T. Congrès National de la
Physique et de ses Applications, no 12, p. 1-10.
STACEY, W, 2018. Nuclear Reactor Physics. "wiley VCH".
TODREAS, N ; KAZIMI, M, 1989. Thermal hydraulic fundamentals. Hemisphere
Corporation. "Nuclear Systems I".

Bibliographi
EL-WAKIL, M, 1971. Nuclear Heat Transport. Cranton Penn USA. "International

Textbook Company".
WILCOX, Dc, 1993. Turbulence Modeling for CFD. DCW Inductries Inc. P. 1-3.
YAKHOT, V ; ORSZAG, SA, 1992. Development of Turbulence Models for Shear Flows
by a Double Expansion technique. Physics of Fluids A Fluid Dynamics, p. 4-15.

Annexes

AnnEXE A
STRUcTURE DU PRogRamme MaTLab

PoUR LE caLcUL DEs TEMPÉRaTUREs en
monodIMEnsIonneL

Annexe A. Structure du programme Matlab pour le calcul
des
Organigramme Matlab
L’organigramme donne les commandes et les étapes du programme Matlab.

Il permet en effet d’afficher les courbes des distributions de température du
fluide caloporteur, de la gaine et de la chair à l’intérieur d’un seul canal
chaud du coeur pour un réacteur de recherche d’une puissance donnée.
Quant aux quantités nucléaire liées à la sûreté nucléaire il permet de calculer
le flux de chaleur critique et moyen pour un ONB tout le long du canal. Le
programme est conçu à partir d’un script principale et un script fonction
qui donne la méthode du point fixe qui sert à calculer le flux de chaleur moyen.
Les températures dans le canal sont évaluées à l’aide du modèle

mathématique El-Wakil (1971). Pour se faire, on crée une boucle dont
le nombre d’itération est égal à la taille du vecteur z qui donne les points
pour lesquelles on calcule les quantités physique souhaitées. On écrit les
équations de ces quantités physique à l’intérieur de la boucle de façon à ce
qu’à chaque itération on attribue à un point du canal la valeur de celle ci. En
plus de ces équations, on ajoute dans la boucle les équations des propriétés
physiques du fluide en fonction de la température puisqu’on les retrouve
dans les équations du modèle mathématique choisi.
Le flux de chaleur moyen est calculé en dehors de la boucle car il est donné
dans une équation non linéaire qu’on résout à l’aide d’une fonction créée au
préalable et qui donne la méthode du point fixe.
Le deuxième modèle mathématique suit les mêmes étapes du programme

Matlab. La différence se situe dans l’expression des équations qui donnent
les températures du caloporteur, de la gaine et de la chair.

des
Programme
Principal
Entrer les données de

dimensionnement des plaques de
combustible et leurs données
Initialiser les données qui varient
en fonction de la
Créer un vecteur
i=1
NON OUI
i<=lenght
Tf(i)=Tf1+(((qc*Am*Le)/(pi*Cpi*((Fd*Wt*rhoi)/(F PTN))))*(sin(pi*z(i)/Le)
ρ(i)=1003.1-(0.1532*Ti)-(0.003*(Ti^2))
landai=0.5551+(0.00219*Ti)-(9.32401*10^(-6)*(Ti^2))
Fonction pointFixe
µ(i)=0.0015-(3.1313*10^(-5)*Ti)+(3.01339*10^(-7)*Ti^2)- (1.10053*10^(-9)*Ti^(3
Cp(i)=(4186)-((0.52402)*Ti)+((0.0080)*(Ti)^2)
Tc(i)=Tf(i)+((qc*em)/(2*hi))*(cos(pi*z(i)/Le))
Tm(i)=Tc(i)+(qc*((em^2/(8*km))+((em*ec)/(2*kc)))*cos(pi*z(i)/Le))
g(Qav)=10*(Tsat- Ti+((5/9)*((9.23*Fax*Frad*Qav)/(P^1.156))^((P^0.0284)/(216))))
/((2*10^5*Frad*LH*lm)/(rho*U*ech*Cp*lch))+(100*Fax*Frad/h))
[x,N]=pointFixe(g,x0,epsilon)
QcrM(i)=(151*(1+0.01112*U)*(1+0.00914*(Tsat-Ti))*(1+0.0009*P))/(1-
(((27.6*LH*lm)/(rho*U*ech*Cp*Fd*Wt))*(1+0.01112*U)*(1+0.0009*P)));
Dessiner les graphes Tf ,Tc et Tm, Qcr et RFTC

des
Script Matlab pour un réacteur nucléaire de recherche 2MW
Les figures A.2, A.3, A.4 et A.5 montrent le script Matlab qui a servi
dans le calcul thermohydraulique en monodimensionnel dans le chapitre II.
FIGURE A.2 – Code Matlab de la ligne 1 à 31
FIGURE A.3 – Code Matlab de la ligne 32 à 62.

des

AnnEXE B
ÉTaPEs de sIMULaTIon sUR ComsoL En

RÉgIME sTaTIonnaIRE ET TRansIToIRE

Annexe B. Étapes de simulation sur Comsol en régime stationnaire
et
Les étapes de simulation stationnaire et transitoire sur Comsol sont les

mêmes excepté pour la dernière étape "study" où ça diffère. En stationnaire
on choisi le solveur stationnaire alors qu’en transitoire on choisit deux
solveurs "initiation à la paroi" et "temporel" car le modèle de turbulence
choisi est le modèle SST K −ω. Les étapes de simulation sont décrites sur les
figures qui suivent :
a) Géométrie
La première étape consiste à créer la géométrie bidimensionnelle. Pour
cela, il suffit de créer deux rectangles pour la chair et le caloporteur et
un polygone pour la gaine ensuite on appuis sur "Union" ( figure B.1).
FIGURE B.1 – Création de la géométrie sur COMSOL
Puis on construit la séquence (figure B.2).
FIGURE B.2 – Fenêtre Géométrie (COMSOL).
b) Matériaux
La prochaine étape consiste à choisir le matériaux de chaque composant

comme cela est montré sur la figure B.3.

et
FIGURE B.3 – Choix des matériaux des composant du canal (COMSOL)
Pour le caloporteur, on choisit l’eau et on sélectionne la géométrie

appropriée.
(b) Géométrie du fluide caloporteur
(a) Propriétés du fluide caloporteur.
FIGURE B.4 – Choix du matériau pour le caloporteur.
Pour la gaine, on choisit l’aluminium 6061 et on sélectionne la

géométrie appropriée.

et
(b) Géométrie de la gaine
(a) Propriétés de la gaine
FIGURE B.5 – Choix du matériau pour la gaine
Pour la chair, on choisit l’Uranium et on sélectionne la géométrie

appropriée
(b) Géométrie de la chair
(a) Propriétés de la chair.
FIGURE B.6 – Choix du matériau pour la chair.
c) Modèle de turbulence

et
Le modèle de turbulence choisi pour le canal est le modèle K −ω SST

de la famille RANS (figure B.7).
FIGURE B.7 – Choix du modèle de turbulence.
d) Conditions aux limites
On a des conditions à la paroi, à l’entrée, à la sortie et puis on a une

condition de symétrie car la géométrie a été modélisée à moitié (figure
B.8). On fait rentrer la force de gravité.
FIGURE B.8 – Fenêtre de l’écoulement turbulent.
A la paroi : On indique la paroi qui sépare le fluide du solide sur la

géométrie (figure B.9)
FIGURE B.9 – Indication de la paroi.
A l’entrée : On indique la vitesse initiale (figure B.10).

et
FIGURE B.10 – Vitesse à l’entrée du canal.
La force appliquée On active la force de gravité car il s’agit d’un

écoulement verticale (figure B.11).
FIGURE B.11 – Vecteur force de gravité.
A la sortie On indique la valeur de la pression à la sortie (figure B.12).
FIGURE B.12 – Indication de la pression à la sortie du canal.
Symétrie On indique la ligne qui représente la symétrie en bleue

(figure B.13).

et
FIGURE B.13 – Axe de symétrie sur la géométrie.
e) Transfert de chaleur
Il y a transfert de chaleur entre fluide et solide. Par conséquent, on

indique les solides et les fluides qui rentrent en jeu. Dans notre cas il
s’agit de la gaine, la chair et le caloporteur et on impose un flux de
chaleur linéaire provenant de la chair (figure B.14).
FIGURE B.14 – Fenêtre de la physique : transfert de chaleur entre solide et fluide.
Pour établir un transfert de chaleur il faut entrer la température initiale

du caloporteur (figure B.15).
FIGURE B.15 – Température d’entrée du caloporteur.
Code pour moyenner les températures
Pour obtenir les profils des température du caloporteur et de la gaine,

on a élaboré un code Matlab pour moyenner les températures des résultats
des deux simulations stationnaire et transitoire. Les figures (B.16) et (B.17)
présentent le script du code Matlab de calcul de la moyenne en stationnaire.

et
FIGURE B.16 – Script de calcul de la moyenne à l’état stationnaire (ligne 1-31)
FIGURE B.17 – Script de calcul de la moyenne à l’état stationnaire (ligne 31-59).
Les figures (B.18), (B.19) et (B.20) présentent le script du code Matlab de

calcul de la moyenne en transitoire à 1s, 2s et 3s.

et
FIGURE B.18 – Script de calcul de la moyenne à l’état transitoire à 1s
FIGURE B.19 – Script de calcul de la moyenne à l’état transitoire à 2s.
FIGURE B.20 – Script de calcul de la moyenne à l’état transitoire à 3s.

AnnEXE C
ÉTaPEs de sIMULaTIon FLUENT en

RÉgIME sTaTIonnaIRE

Annexe C. Étapes de simulation FLUENT en régime
Les premières étapes d’une simulation CFD consiste à créer la géométrie

et à générer son maillage. Sur Ansys Workbench, il est possible de le faire
en utilisant spaceclaim, design modler ou en important la géométrie d’un
autre logiciel CAO comme dans le cas de cette simulation, la géométrie a été
réalisée sur Solid Works. Sur la figure C.1 le tableau du solver fluent.
FIGURE C.1 – Tableau Fluent sur Workbench.
Après l’étape du maillage vient l’étape de la configuration. On commence

par préciser qu’il s’agit d’une simulation stationnaire sans variation dans le
temps puis on introduit la force de pesanteur dans la direction de l’axe x
(figure C.2) puisque l’écoulement se fait verticalement de haut en bas.
FIGURE C.2 – Configuration générale.
Dans la fenêtre Models on précise qu’il y a un transfère de chaleur en

cochant la case "Energy" puis on détermine le type de l’écoulement qui est
turbulent en choisissant le modèle de turbulence k − ω SST.

FIGURE C.3 – Fenêtre "models".
Il est aussi important de préciser le type de matériaux de chaque solide et

fluide comme sur la figure C.4. On a défini la chair comme étant de
l’uranium
, la gaine est fabriquée en aluminium 6061 et le caloporteur utilisé n’est
autre que de l’eau légère.
FIGURE C.4 – Matériaux
Sur la figure C.5, les propriétés de l’Uranium attribué à la chair.
FIGURE C.5 – Propriétés de la chair.
Sur la figure C.6, les propriétés de l’aluminium attribué à la gaine.

FIGURE C.6 – Propriétés de la gaine.
Sur la figure C.7, les propriétés de l’eau attribuées au caloporteur.
FIGURE C.7 – Propriétés de l’eau.
Dans la fenêtre Cell zone sont listés les différents composant du canal (C.8)
FIGURE C.8 – Fenêtre cell zone.
En ouvrant la fenêtre des conditions aux limites, on peut déterminer les

conditions aux limites énoncées précédemment (figure C.9).

FIGURE C.9 – Conditions aux limites.
Les conditions aux limites sont la vitesse d’entrée, la température d’entrée,

la pression de sortie et la source de chaleur qui est sinusoïdale. On fait donc
entrer la vitesse comme montré sur la figure C.10.
FIGURE C.10 – Vitesse à l’entrée du canal.
La pression à la sortie est entrée comme suit (figure C.11) :

FIGURE C.11 – Pression à la sortie du canal.
La source de chaleur est la chair (figure C.12). Puisqu’elle est sinusoïdale il

faut faire entrer une fonction UDF qu’on donne sur la figure C.13.
FIGURE C.12 – Source de chaleur
La fonction UDF s’écrit sous langage C. Il faut alors l’interpréter pour

pouvoir l’utiliser.
FIGURE C.13 – Fonction UDF de la source de chaleur.

La détermination de la méthode de résolution est faites en choisissant le

schéma SIMPLE (figure C.14).
FIGURE C.14 – Methods.
La convergence est assurée lorsque les erreurs (risiduals) deviennent inférieurs

à 10−6 pour l’énergie et 10−4 pour le reste (figure C.15).
Erreurs residuelles
FIGURE C.15 – Residual.
On précède par la suite à l’initialisation en choisissant une initialisation

hybride. Cette initialisation permet de définir des valeurs initiales pour les
variables d’écoulement et d’initialiser la solution à l’aide de ces valeurs.

(figure C.16).
Initialisation
FIGURE C.16 – Initialisation.
Pour commencer le calcul on donne un nombre d’itération assez grand mais

le calcul s’arrête dès qu’il y a convergence (figure C.17).
FIGURE C.17 – Calculs et itération.

AnnEXE D
ÉTaPEs de sIMULaTIon FLUENT en

RÉgIME TRansIToIRE

Annexe D. Étapes de simulation FLUENT en régime
Les étapes de simulation en transitoire sont les mêmes qu’en stationnaire

hormis quelques modifications qu’on expliquera dans les étapes suivantes :
Sur la fenêtre Configuration générale (figure C.2), il faut sélectionner

Transient au lieu de steady. Puis s’en suit les mêmes étapes jusqu’à l’étape
des conditions aux limites où la vitesse d’entrée dans notre cas doit diminuer
de 20 % après un certain temps. Sa fonction UDF est présentée sur la figure
D.1.
FIGURE D.1 – Fonction UDF de la vitesse.
Sur la fenetre "Run calculation" on choisit 30000 itérations avec un pas de

temps de 0,001 secondes comme cela est montré sur la figure D.2

Annexe D. Étapes de simulation FLUENT en régime
FIGURE D.2 – Calculs et itérations.
La figure D.3 montre le script du code Malab servant à moyenner pour

obtenir les profils de températures sur Fluent
FIGURE D.3 – Script du code matlab de calcul de la moyenne en stationnaire.

Analyse Thermohydraulique Transitoire Dans Un Canal Du Coeur D'un Réacteur Nucléaire de Recherche de Type MTR

Transféré par

Droits d'auteur :

Formats disponibles

Analyse Thermohydraulique Transitoire Dans Un Canal Du Coeur D'un Réacteur Nucléaire de Recherche de Type MTR

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

Analyse Thermohydraulique Transitoire Dans Un Canal Du Coeur D'un Réacteur Nucléaire de Recherche de Type MTR

Transféré par

Droits d'auteur :

Formats disponibles

RÉ PUBLIqUE ALGÉ RIEnnE DÉ MOCRATIqUE ET POPULAIRE

mInISTÈ RE DE L’EnSEIGnEMEnT SUPÉ RIEUR ET DE LA RECHERCHE SCIEnTIfIqUE

École Nationale Polytechnique

Laboratoire de Recherches Sciences de l’Eau

Projet de fin d’études

Mémoire de PFE TIHAL et BENSEMMANE 150

École Nationale Polytechnique

Laboratoire de Recherches Sciences de l’Eau

Projet de fin d’études

On remercie Dieu le tout-puissant de nous avoir donné le courage, la santé

On remercie vivement notre cher encadreur, docteur Kamel SIDI ALI de

On remercie également notre chère professeur et directrice de laboratoire

On tient à adresser nos sincères remerciements à Pr.Smaili Arezki et

On remercie chaleureusement nos familles et tous nos amis pour leur

On souhaite aussi remercier tous les enseignants de l’École Nationale

Que tous ceux qui ont contribué de près ou de loin à la réalisation de ce

Merci à toutes et à tous.

À mon frères Hamoud et ma grand-mère qui m’ont soutenu durant toute

À tous ceux qui me sont chers, je dédie ce modeste travail.

Aucune dédicace ne saurait exprimer mon amour, ma gratitude et ma

A mon cher Papa

A qui je dois ma vie, ma réussite et tout mon respect. Tu as toujours été là

A mes chères sœurs et mon frère

A mes chères Tantes et ma cousine

Zou, Mimi et Lila que j’aime, je vous considère comme ma deuxième

A mon cher Papi

A ma chère Binôme Yasmine

Table des figures 9

Liste des abréviations 1

1 Généralités sur les réacteurs nucléaires de recherche ;

2 La thermohydraulique d’un canal du coeur d’un réacteur

Mémoire de PFE TIHAL et BENSEMMANE 6

2.2 Caractéristiques physiques et de design des réacteurs de

3 Modélisation thermohydraulique en 2D d’un canal du

Mémoire de PFE TIHAL et 7

Conclusion générale 118

A Structure du programme Matlab pour le calcul des

B Étapes de simulation sur Comsol en régime stationnaire et

C Étapes de simulation FLUENT en régime stationnaire 139

D Étapes de simulation FLUENT en régime transitoire 147

Mémoire de PFE TIHAL et 8

1.1 Schéma d’une réaction de fission.........................................................23

Mémoire de PFE TIHAL et 9

2.9 Distribution des températures en fonction de la longueur

Mémoire de PFE TIHAL et 1

3.7 Distribution de la température de la gaine le long de l’axe y. . 85

Mémoire de PFE TIHAL et 1

A.3 Code Matlab de la ligne 32 à 62.....................................................128

Mémoire de PFE TIHAL et 1

Mémoire de PFE TIHAL et 1

1.1 Répartition de l’énergie provenant d’une réaction de fission . . 23

Mémoire de PFE TIHAL et 1

Am Section transversale de la chair (m2)

Mémoire de PFE TIHAL et 1

Mémoire de PFE TIHAL et 1

Lh Longueur active de la plaque (m)

Mémoire de PFE TIHAL et 1

u¯i Vitesse moyenne (m/s)

A→ Vecteur de surface (m2 )