Modélisation des phénomènes de dissolution lors des phases précoces et avancées d'un accident grave de réacteur nucléaire

Julien Belloni

Description

Sous la direction de Benoît Goyeau
Thèse soutenue le 09 février 2009: Ecole centrale Paris
Cette thèse porte sur la modélisation des phénomènes de dissolution lors de la phase précoce d'un accident grave de réacteur nucléaire. L'étude s'intéresse à la dissolution de céramiques solides (ZrO2 et UO2) par un métal liquide (Zr). En cas d'hypothétique accident grave dans un Réacteur à Eau Pressurisée, les phénomènes de dissolution jouent un rôle primordial dans l'aggravation de la dégradation et la fusion des oxydes à des températures parfois largement inférieures à leur température de fusion normale. Cela concerne en particulier les principaux constituants des crayons combustibles : les pastilles d'UO2 ainsi que la couche de ZrO2 formée sur la surface externe de la gaine qui vont subir une dissolution par le zircaloy des gaines à partir de 2100 K (la température de fusion de ces oxydes étant supérieure à 2800 K). Dans l'état actuel des connaissances, on peut supposer que les phénomènes de dissolution sont responsables, d'une part, de l'effondrement prématuré des crayons combustibles dans le cœur et, d'autre part, de la formation rapide d'un bain liquide en fond de cuve si des oxydes de fer sont présents. De nombreuses études expérimentales ont été menées sur ce sujet mais la modélisation n'est pas encore satisfaisante à ce jour. Les modèles actuels sont essentiellement des modèles 1D qui ne prennent pas en compte de façon explicite la convection naturelle ni la présence d'une zone de transition diphasique au niveau de l'interface métal / céramique. Un modèle 2D, décrivant les transferts de quantité de mouvement, de chaleur et de masse, établi par prise de moyenne volumique des équations microscopiques de transport a été développé. Ce modèle est basé sur des hypothèses d'équilibre thermique local et de non-équilibre massique local et s'inspire d'un modèle de solidification de mélanges binaires (P. Bousquet-Mélou, 2000). Sur le plan théorique, l'approche est identique. Cependant, la dissolution introduit des contraintes physico-numériques supplémentaires qu'il a fallu prendre en compte. La prise en compte des deux aspects mentionnés (convection et zone diphasique) constitue une nouveauté significative par rapport aux modèles existants. Le modèle a été ensuite étendu au cas ternaire grâce à des développements nouveaux permettant de résoudre certaines difficultés supplémentaires (ajout d'une équation d'espèce, relations d'équilibre non bijectives, plusieurs coefficients de diffusion à déterminer, indétermination sur les flux à l'interface, densité du solide non constante...) afin d'étudier la dissolution UO2 / Zr. Une validation des modèles obtenus est faite à partir de résultats expérimentaux de dissolution de creusets en ZrO2 et UO2 par le zircaloy liquide et montrent la bonne qualité prédictive de notre modèle.
-dissolution
-milieux poreux
-modèle 2D
This thesis focuses on the modeling of dissolution phenomena during the early phase of a serious accident of nuclear reactor. The study focuses on the dissolution of solid ceramic (ZrO2 and UO2) by a liquid metal (Zr). In case of hypothetical severe accident in a pressurized water reactor, the phenomena of dissolution play a role in the further decay and melting oxides at temperatures sometimes significantly below the normal melting temperature. This concerns in particular the main constituents of fuel rods : pellets of UO2 and ZrO2 layer formed on the outer surface of the sheath that will suffer by the dissolution of zircaloy sheaths from 2100 K (the melting temperature of these oxides being greater than 2800 K). In the present state of knowledge, we can assume that dissolution phenomena are responsible for, on the one hand, the premature collapse of the fuel rods in the heart and, on the other hand, the rapid formation of a liquid bath in the bottom of the tank where iron oxides are present. Many experimental studies have been conducted on this subject but modeling is not yet satisfactory today. Current models are essentially 1D models that do not explicitly account for natural convection nore the presence of a biphasic transition zone at the interface metal / ceramic. A 2D model describing the transfer of momentum, heat and mass, determined by taking average density of microscopic transport equations has been developed. This model is based on assumptions of local thermal equilibrium and non-local mass balance and it is based on a model of solidification of binary mixtures (P. Bousquet-Mélou, 2000). In theory, the approach is identical. However, the dissolution introduced physico-digital extra constraint that had to be taken into account. Taking into account the two aspects mentioned (convection and two-phase zone) is a novelty in comparison to existing models. The model was then extended to the ternary case with new developments to resolve certain difficulties (addition of an equation of this case, non-bijective equilibrium relations, several diffusion coefficients to be determined, uncertainty on flows in interface, density of the solid non-constant ...) to study the dissolution UO2 / Zr. Validation of models is obtained from experimental results of dissolution of crucibles of UO2 and ZrO2 by liquid zircaloy and show the quality of our predictive model.
Source: http://www.theses.fr/2009ECAP0008/document

Sujets

• Dissolution
• Porosité

Informations

• Publié par: Thesee
• Nombre de lectures: 45
• Langue: Français
• Poids de l'ouvrage: 4 Mo

Extrait

ÉCOLE CENTRALE DES ARTS
ET MANUFACTURES
« ÉCOLE CENTRALE PARIS »


THÈSE
présentée par

Julien BELLONI

pour l’obtention du

GRADE DE DOCTEUR

Spécialité : Energétique

Laboratoire d’accueil : EM2C

SUJET : MODELISATION DES PHENOMENES DE DISSOLUTION LORS
DE LA PHASE PRECOCE D’UN ACCIDENT GRAVE DE REACTEUR
NUCLEAIRE




Soutenue le 9 Février 2009

devant un jury composé de :

Florian FICHOT IRSN-DPAM, Cadarache (Ingénieur) Examinateur
Dominique GOBIN FAST-Paris 6, Orsay (D.R.) Président
Benoît GOYEAU EM2C-ECP, Châtenay-Malabry Examinateur
(Professeur, co-directeur de Thèse)
Christophe JOURNEAU CEA-DEN, Cadarache (Ingénieur) Rapporteur
Michel QUINTARD IMFT-INPT, Toulouse Examinateur
(D.R.,co-directeur de Thèse)
Gérard VIGNOLES LCTS-UB1, Bordeaux (Professeur) Rapporteur


2009-ECAP0008



tel-00453295, version 1 - 4 Feb 2010
tel-00453295, version 1 - 4 Feb 2010ma
ts,
famille.
Sophie
A
et
paren
mes
A
tel-00453295, version 1 - 4 Feb 2010tel-00453295, version 1 - 4 Feb 2010he.
sincères
fut
ts
é
Je
bres
tiens
d'acquérir
tout
hel
d'ab
disp
ord
ail
à
thèse.
remercier
à
M.
le
Dominique
de
Gobin
Fluides
d'a
thèse
v
à
oir
b

LESAM.
de
soutien
présider

mon
de
jury

de

thèse
n
ainsi
Benoît
que
P
Christophe
l'Institut
Journeau
(
et
on
Gérard

Vignoles
b
p
ail
our
à
en

a

v
je
oir
oratoire
été
trois
les

rapp
années.
orteurs.
réel
Leurs
v
remarques

m'on
p
t
eaucoup
été
t
très
du
utiles
J'adresse
p
ts
our
y
p
Cen

et

tard
mémoire.

Cette
T
thèse
andmer
a
).
été
dirigé
réalisée
leurs
au
que
sein
on
du
app
Lab
tra
oratoire

d'Etude
je
et
que
de
apprécié
Sim
tra
ulation
que
des
règne
A


les
ts

Ma
our
jeurs
t
(LESAM)
de
du
au
Service
de
d'Etude
trois
et
Ce
de
un
Mo
plaisir
délisation
tra
du
ailler
Com
ses
bustible
et
en
m'a
situation
ermis
A
b

de

notammen
telle
dans
de
domaine
l'IRSN
monde
à
ucléaire.

mes
he.
remerciemen
Je
à
remercie
Go
tout
eau
particulièremen
l'Ecole
t
trale
Jo
aris
ëlle

Fleurot
Quin
qui
de
m'a
de

de
dans
de
son
oulouse
lab
gr
oratoire.

Je
!
tiens
Ils
à
t
remercier

sp
et
écialemen

t
ainsi
Florian
leur

onibilité
hot,
t
ingénieur
eaucoup
de
orté


herc
v
he
de
et
herc
exp
Enn,
ert
tiens
à
dire
l'IRSN,
j'ai
qui
eaucoup
a
les

de

v
tra
ainsi
v
l'am
ail
qui
de
au
thèse
A
et
titre
p
remercie
our
mem
tout
de
le
lab
sa
p
v
leur
oir
duran
qu'il

a
années
pu
orter
m'app
Remerciemen
tel-00453295, version 1 - 4 Feb 2010tel-00453295, version 1 - 4 Feb 2010
le
mouv
Cette
dèle
thèse
v
p
est
orte
La
sur
èce,
la
Zircalo
mo
métal
délisation
v
des

phénomènes
duit
de

dissolution
eaux
lors
les
de
résultats
la
dèle.
phase
transition
préco
an

masse,
d'un
a

thermique
t
de
gra
tique.
v
supplémen
e
men
de
ort
réacteur
grâce
n
dicultés
ucléaire.

L'étude
non
s'in
dèles
téresse
en
à
b
la
explicite
dissolution

de
eau

mo
solides
de
(ZrO

ternaire
de
et

UO
é.
binaire,
h
)
de
par
mo
un
.
métal

liquide
la
(Zr).
tes
En
fallu

des
d'h
ection
yp

othétique
existan

étendu
den
elopp
t
de
gra
jout
v

e
à
dans
terface,
un
d'étudier
Réacteur
v
à
faite
Eau
de
Pressurisée,
UO
les
mon
phénomènes
e
de
dissolution,
dissolution

jouen
ni
t
zone
un
au
rôle
l'in
primordial

dans
2D,
l'aggra
les
v
tité
ation
t,
de
et
la
par
dégradation
y
et
des
la
transp
fusion
dév
des
mo-
o
sur
xydes
othèses
à

des
massique
temp
s'inspire
é-
de
ratures
binaires
parfois
2000).
largemen
théorique,
t
est
inférieures
endan
à
in
leur

temp
ysico-n
érature
qu'il
de
en
fusion
en
normale.
asp
Cela
(con

zone
en
nouv
particulier
par
les
mo

Le

été
ts

des
des

ts
y
ermettan
ons


taires
bustibles
équation
:
d'équilibre
les
plusieurs
pastilles
de
d'UO
indétermination
non-équilibre,
à
ainsi
du
que
te...)
la
dissolution

Zr.
he
des
de
us
ZrO
partir
ection,
érimen
formée
de
sur
diphasique,
la
par
surface
liquide
externe
t
de
qualité
la
notre
gaine

qui
p
v
la
on
v
t
naturelle
subir
la
une
d'une
dissolution
de
par
diphasique
le
niv
zircalo
de
y
terface
des
/
gaines
Un
à
dèle
partir

de
t
2100
transferts
K
quan
(la
de
temp
emen
érature
de
de
haleur
fusion
de
de
établi

prise
o
mo
xydes
enne
étan
olumique
t
équations
sup
de
érieure
ort
à
été
2800
elopp
K).
Ce
Dans
dèle
l'état
basé
actuel
des
des
yp

d'équilibre
on
lo
p
et
eut
non-équilibre
sup-
lo
p
et
oser
d'un
que
dèle
les
solidication
phénomènes
mélanges
de
(P
dissolution
Bousquet-Mélou,
son
Sur
t
plan
resp
l'appro
onsables,
he
d'une
iden
part,
Cep
de
t,
l'eondremen
dissolution
t
tro
prématuré
des
des
train

ph
y
umériques
ons
taires

a
bustibles
prendre
dans

le
prise


eur
deux
et,
ects
d'autre
tionnés
part,
v
de
et
la
diphasique)
formation
une
ra-
eauté
pide
e
d'un
rapp
bain
aux
liquide
dèles
en
ts.
fond
mo
de
a

ensuite
e
au
si
ternaire
des
à
o
dév
xydes
emen
de
nouv
fer
p
son
t
t
résoudre
présen
taines
ts.
supplémen
De
(a
nom
d'une
breuses
d'esp