Movatterモバイル変換


[0]ホーム

URL:


Aller au contenu
Wikipédial'encyclopédie libre
Rechercher

Empoisonnement au xénon

Un article de Wikipédia, l'encyclopédie libre.

Dans unréacteur nucléaire, l'empoisonnement au xénon est le phénomène de production et d'accumulation duxénon 135, un puissantabsorbeur de neutron, étouffant laréaction nucléaire d'un réacteur à l'arrêt, ou provoquant des oscillations de puissance dans les réacteurs de grande taille.

Ce phénomène est un desfacteurs ayant conduit à l'accident de Tchernobyl.

La chaîne Xe

[modifier |modifier le code]

Production du xénon 135

[modifier |modifier le code]
Chaine xénon.

Lexénon 135 est unradioisotope (demi-vie : 9,2 h[1]), qui apparaît dans lesproduits de fission de lafission nucléaire d'unematière fissile (matière employée pour le fonctionnement d'unréacteur nucléaire,uranium 235 ouplutonium pour les réacteurs de troisième génération).

  • Le135Xe se forme surtout indirectement, parfiliation radioactive, principalement à partir de l'iode 135 de période 6,7 h[1] (lequel peut lui-même être produit par la désintégration de l'éphémèreTellure 135, voire en amont par des parents à durée de vie encore plus brève, sans intérêt pour le fonctionnement d'unréacteur nucléaire).
  • Pour mémoire, le135Xe peut apparaître aussi de manière marginale comme produit de fission direct, avec unrendement de fission de 0,4 % pour l'uranium 235.

Dans la fission de l'uranium 235, lerendement de fission pour le poids atomique de 135 est de 6,6 % d'atomes produits par fission d'isotope fissile pour des neutrons thermiques (6,3 % pour des neutrons rapides)[2].

Disparition du xénon 135

[modifier |modifier le code]

Le135Xe peut donner ducésium 135, pardécroissance radioactive de période 9,2 h, mais peut également donner enréacteur nucléaire du xénon 136 parcapture d'un neutron, quand il est soumis à unflux neutronique. La proportion relative de ces deux consommations dépend du flux de neutrons (donc de la puissance de fonctionnement et du type de réacteur) : plus le flux est important, et moins le xénon aura le temps de se désintégrer en césium avant de capturer un neutron supplémentaire.

Par la suite, ces deux noyaux136Xe et135Cs sont pratiquement stables et de section efficace négligeable. Pour mémoire, le135Cs est unproduit de fission faiblement radioactif (d'une demi-vie de 2,6 millions d'années), ce qui pose d'autres problèmes pour le gérer en tant quedéchet radioactif, mais il n'a pas d'incidence sur le fonctionnement d'unréacteur nucléaire.

Les conditions de fonctionnement d'un réacteur nucléaire sont proches du point où le rapport136Xe/135Cs bascule : pour unflux neutronique en réacteur de l'ordre de 1 × 1014 n cm−2 s−1 (un peu élevé pour un réacteur de puissance, mais en-dessous d'un réacteur d'irradiation), la section efficace du135Xe, de l'ordre de σ=2,65 × 10−18 cm2, conduit à une probabilité de capture de 2,65 × 10−4 s−1, c'est-à-dire une demi-vie par capture de l'ordre de l'heure, neuf fois plus que la demi-vie par désintégration. Les demi-vies ayant une incidence exponentielle, la capture l'emporte alors dans une proportion de 29 contre un, soit à 99,8 % : en fonctionnement dans ces conditions, la quasi-totalité du135Xe fait l'objet d'une capture neutronique. Inversement, si le flux neutronique est dix fois moins intense (comme c'est le cas pour des réacteurscandu à pleine puissance), les deux demi-vies seront du mêmeordre de grandeur, et le rapport136Xe/135Cs sera donc proche de 50 %.

Dynamique du xénon 135

[modifier |modifier le code]

Dans unréacteur nucléaire en fonctionnement stabilisé, la quantité de xénon produit (surtout par la désintégration de l'iode 135) est asymptotiquement égale à la quantité de xénon qui disparaît (par capture neutronique, ou par désintégration) ; mais cet équilibre dépend du flux neutronique, et est très sensible à ses variations, donc aux variations de puissance du réacteur.

La dynamique du xénon en réacteur est gouvernée par les particularités suivantes :

L'empoisonnement ausamarium 149 (section efficace de capture = 41 140 barns) dont le père est le prométhium 149 (période = 53,08 h) suit un mécanisme similaire, si ce n'est que cet élément n'est pas radioactif : après son accumulation à l'arrêt, il ne disparaît pas ; et un minimum d'excès de réactivité est nécessaire par rapport à la situation d'équilibre pour redémarrer le réacteur, même après un arrêt prolongé. L'excès de samarium 149 par rapport à sa valeur d'équilibre disparaît lentement après redémarrage, « mangé par le flux ». Les phénomènes se déroulent beaucoup plus lentement que dans le cas du xénon et les effets en réactivité sont nettement moindres. Bien que leprométhium 149 soit présent à l'équilibre en plus grande quantité que l'iode 135 dans le réacteur[Note 1], la section efficace du samarium 149 est très inférieure à celle du xénon 135 et l'anti-réactivité du samarium à l'équilibre est environ 3 fois moindre que celle du xénon 135.

Effet du xénon sur le pilotage d'un réacteur

[modifier |modifier le code]

Étouffement au démarrage

[modifier |modifier le code]

Au premier démarrage d'un réacteur nucléaire thermique, la concentration en135Xe est nulle. La concentration en xénon va progressivement augmenter (si le flux neutronique est constant) jusqu'à atteindre une valeur d'équilibre au bout d'un jour ou deux. Cet empoisonnement en fonctionnement est de l'ordre de 3 000 pcm pour un réacteurREP de 900 MW.

De ce fait, après la divergence initiale, un réacteur perd progressivement saréactivité, ce qui doit être compensé par le pilotage (soit pour un réacteur REP: levée desbarres de contrôle et/ou dilution de la concentration enacide borique) faute de quoi le réacteur s'arrêtera de lui-même.

Cet effet fut découvert de manière imprévue dans leRéacteur B, premier réacteur de puissance industrielle, réalisé dans lecomplexe nucléaire de Hanford. Après sa divergence initiale, les opérateurs eurent la surprise de voir la puissance du réacteur chuter, et la réaction s'arrêter, le réacteur restant impossible à démarrer pendant quelques jours[4]. Après analyse du phénomène, le problème put être corrigé, les ingénieurs concepteurs du réacteur ayant prévu des marges suffisantes pour pouvoir rajouter des éléments combustibles par rapport au plan initial, ce qui permit de restaurer la réactivité de l'ensemble en régime de fonctionnement permanent.

Pic xénon

[modifier |modifier le code]
Perte de réactivité due à l'augmentation transitoire de la concentration en xénon après l'arrêt d'un réacteur nucléaire.

« Toute réduction de puissance cause un pic de xénon transitoire. Plus la réduction de puissance est faible, plus le pic est faible et plus il se produit rapidement »[5]. Le passage brutal à un faible régime peut conduire à l'étouffement du réacteur. Le phénomène est connu en anglais sous le nom deiodine pit, dû à l'aspect de la courbe de réactivité du réacteur après un arrêt (courbe en rouge).

Le pic xénon est particulièrement gênant dans le cas desréacteurs de propulsion navale : une baisse importante de puissance ou un arrêt au port ou a fortiori survenant à la mer, ne devant pas conduire à rendre le réacteur donc le bâtiment indisponible. De ce fait, les réacteurs de propulsion navale doivent conserver une réserve de réactivité beaucoup plus importante qu'un réacteur civil : leurcombustible est d'un enrichissement supérieur, et il ne peut pas atteindre lestaux de combustion réalisés couramment pour les réacteurs civils.

Description du phénomène

[modifier |modifier le code]

À l'arrêt (ou lors d'une réduction de charge) d'un réacteur, la production d'iode 135 s'arrête (ou diminue), ainsi que la consommation du135Xe. Le « réservoir d'iode 135 » se transforme alors encore en135Xe (avec une demi-vie de 6,7 h). La demi-vie du xénon 135 (avec une demi-vie de 9,2 h) étant supérieure à celle de son père, il y a un effet d'embouteillage : le xénon 135 s'élimine nettement moins vite qu'il n'est formé. Il commence alors par s'accumuler (augmentant ainsi l'anti-réactivité associée), jusqu'à ce que la source d'iode 135 se tarisse et que la désintégration propre du xénon 135 soit suffisante pour l'éliminer progressivement. Le niveau du réservoir xénon passe ainsi par un maximum (au bout d'une dizaine d'heures lors d'un arrêt, ou moins lors d'une réduction de charge), puis tend plus lentement vers zéro.Le135Xe étant unpoison neutronique de premier ordre, cette augmentation transitoire de la teneur en xénon va correspondre à une perte très importante en réactivité du réacteur, pouvant aller jusqu'à interdire le redémarrage du réacteur. Ainsi, à l'arrêt d'un réacteur, quand l'anti-réactivité apportée par le135Xe au « pic xénon » est trop importante et que les réserves de réactivité sont insuffisantes (par exemple, en fin de vie du réacteur), le réacteur ne peut plus rediverger pendant quelques heures.

Instabilité d'un redémarrage au pic xénon

[modifier |modifier le code]
Croquis figurant l'évolution de l'empoisonnement xénon lors d'un redémarrage au « pic ».

Lors d'un redémarrage au moment du « pic xénon », la disparition du xénon par interaction avec le flux provoque un apport de réactivité important, d'autant plus rapide que la puissance est grande et que la quantité de xénon restante est importante. La montée en puissance du réacteur dans ces conditions doit donc être suffisamment lente pour que la variation de réactivité puisse en permanence être compensée par l'action sur les barres de contrôle, faute de quoi le réacteur pourrait s'emballer.

Un non-respect des consignes de sécurité sur l'extraction des barres de contrôle neutronique et une mauvaise préparation de l'essai comptent parmi les évènements ayant provoqué l'explosion de ce réacteur le. Cet étouffement de la réaction nucléaire s'était produit sur le réacteur 4 de lacentrale nucléaire de Tchernobyl à la suite d'une réduction volontaire de puissance. Mais ensuite, afin de maintenir la puissance du réacteur dans le but de réaliser un essai sur la partie non nucléaire de l'installation, les opérateurs ont relevé les barres de contrôle au-delà de la limite autorisée ce qui a fait entrer le réacteur dans une zone de fonctionnement instable (coefficient de vide positif) et provoqué son auto-emballement[6].

Instabilités aux changements de régime

[modifier |modifier le code]

Après un temps de fonctionnement suffisant à puissance constante, la concentration en xénon atteindra un premier équilibre dans un réacteur. À ce stade, si la puissance demandée au réacteur augmente d'un palier, l'augmentation de puissance correspond à une augmentation du flux neutronique, qui augmentera à terme tous les flux et toutes les concentrations des éléments de la chaîne xénon (en proportion, la seule différence induite par la puissance de fonctionnement étant le rapport entre les productions de136Xe et135Cs). Le point d'équilibre du xénon se déplacera donc vers une plus grande concentration, puisque l'augmentation du flux doit permettre de détruire plus de xénon. Cependant, ce nouveau point d'équilibre asymptotique ne sera pas atteint directement.

Dans un premier temps, l'augmentation de puissance ne fait augmenter le réservoir d'iode que linéairement par rapport au temps. Comme 95 % du xénon 135 provient de la désintégration de l'iode 135, la quantité supplémentaire de xénon apportée par la plus grande puissance varie donc initialement ent2. En revanche, la vitesse de dégradation du xénon est proportionnelle à la puissance du réacteur : le déficit en xénon résultant de cette destruction est donc initialement proportionnel au temps. L'effet linéaire étant prédominant par rapport à l'effet quadratique, le bilan sera initialement une destruction nette de xénon.

La concentration en xénon suit donc une dynamique paradoxale : dans un premier temps la hausse de puissance entraîne (à l'échelle de l'heure) une diminution transitoire de cette concentration, qui ne rejoint son niveau d'équilibre plus élevé que dans un deuxième temps (à l'échelle de la journée). Par ailleurs, comme le xénon 135 est un poison à neutrons, sa baisse de concentration augmente la réactivité du cœur, et donc la puissance du réacteur, ce qui est compensé par l’insertion des barres de contrôle.

Donc, quand la puissance du réacteur doit être augmentée d'un palier, on verra d’abord les barres de contrôle se lever pour augmenter la puissance, puis se baisser (pour compenser la consommation de xénon) puis repasser au-dessus du niveau initial quand l'équilibre est atteint (pour compenser la quantité de xénon à présent plus importante). Ces variations de réactivité sont suffisamment lentes pour ne pas poser en soi de problème de pilotage, à condition que les réserves de réactivité et d'anti-réactivité soient suffisantes pour maintenir en permanence le réacteur dans une plage de fonctionnement où il peut être piloté, ce qui peut imposer des limites sur la vitesse de variation de la puissance du réacteur.

Lorsque la puissance du réacteur doit être diminuée, le processus ci-dessus est inversé[7]. Toutes choses égales par ailleurs, la baisse de puissance provoque une accumulation de xénon 135 neutrophage, qui tend à diminuer encore plus la puissance. Si le réacteur n'a pas de réserve de réactivité suffisante, cette accumulation peut rendre impossible une stabilisation du réacteur en fonctionnement à faible puissance, après un passage à forte puissance, et provoquer un arrêt du réacteur, par étouffement sous l'effet du pic de xénon non contrôlé.


Oscillations dans la distribution de puissance

[modifier |modifier le code]

L'anti-réactivité locale introduite par la concentration en xénon peut provoquer desoscillations dues au xénon, radiales ou axiales[8],[9],[10]. Le mécanisme de cetterétroaction passe par quatre étapes[4] :

Une augmentation locale du flux neutronique dans le cœur (induite, par exemple, par un retrait des barres de contrôle servant à ajuster la puissance du réacteur) va provoquer un écart (radial ou axial) par rapport à la distribution de puissance.

Comme expliqué précédemment, la section efficace du xénon rend instable la distribution de puissance. Dans les zones où le flux a augmenté, le135Xe est consommé plus rapidement, ce qui augmente la réactivité et donc la puissance locale.

L'ensemble accentue le déséquilibre initial. Le déséquilibre s'accentue jusqu'à ce que l'effet thermique soit suffisant pour compenser la différence de réactivité.

Comme le réacteur est piloté pour fonctionner à puissance globale constante, une augmentation locale de puissance impose une diminution de puissance dans les autres parties du cœur.

Comme expliqué précédemment, la section efficace du xénon rend instable la distribution de puissance.

Dans les zones où la puissance a baissé, l'accumulation locale de135Xe entraîne une baisse de puissance locale. L'ensemble accentue le déséquilibre initial. Le déséquilibre s'accentue jusqu'à ce que l'effet thermique soit suffisant pour compenser la différence de réactivité.

Dans les zones de haut flux, la concentration en135I produit par les fissions augmente alors progressivement. L'augmentation de la concentration en135I finit par être suffisante pour rétablir la production de135Xe, et faire remonter sa concentration au niveau d'équilibre correspondant à la puissance locale, ce qui ramène à une valeur d'équilibre la contre-réactivité due au135Xe.Dans les zones de flux moindre, la concentration en135I produit par les fissions diminue progressivement. La diminution de la concentration en135I ralentit la production de135Xe, et fait diminuer sa concentration, diminuant la contre-réactivité due au135Xe.
Quand l'équilibre de réactivité est atteint, la concentration en135I est au-dessus de sa valeur d'équilibre, et continue à produire plus de135Xe que le flux local ne peut en consommer.

La réactivité de la zone diminue au-delà de sa valeur d'équilibre, conduisant à une baisse locale du flux neutronique dans le cœur, plus ou moins amortie par la baisse de température.

Dans les zones de faible flux, la concentration en135I tombe en-dessous de sa valeur d'équilibre, ce qui entraîne un déficit dans la production de135Xe.

La réactivité de la zone augmente au-delà de sa valeur d'équilibre.

L'instabilité au changement de régime joue alors dans l'autre sens, accentuant la baisse de puissance locale dans la zone précédemment en excès.Le réacteur étant piloté pour fonctionner à puissance constante, cette baisse de puissance locale est compensée par la hausse concomitante de puissance dans les zones où le flux avait été réduit : le déséquilibre (radial ou axial) s'inverse.
Dans la zone à présent de bas flux, étouffées par un excès de135Xe, la production en135I passe en-dessous de sa valeur d'équilibre, et malgré une consommation moindre par capture neutronique, la concentration en135Xe diminue lentement par décroissance radioactive. Au bout d'un certain temps, la contre-réactivité s'élimine et retrouve une valeur d'équilibre, mais dans une zone à présent déficitaire en135I, ce qui provoquera à terme un déficit en135Xe, une augmentation locale du flux, et enclenchera un nouveau cycle.La zone précédemment à bas flux devient à haut flux : la dissymétrie a été inversée par rapport à la situation initiale.

La répétition du cycle déplace les déséquilibres de distribution de puissance dans le cœur, avec une périodicité de l'ordre de 15 heures[11].

Ces oscillations peuvent induire des déséquilibres dans la distribution de puissance d'un facteur trois voire plus[11]. Un tel déséquilibre peut créer un point chaud dans le réacteur et mettre en danger les éléments combustibles.

Dans un réacteur présentant un coefficient de température fortement négatif, ces oscillations sont réduites assez rapidement[11].

Afin de contrôler le déséquilibre de puissance axial, l'instrumentation desREP informe en permanence l'exploitant de sa valeur (axial-offset (Dpmax)).

À l'apparition d'une oscillation xénon, l'exploitant du réacteur se doit de la contenir par un « contre-pilotage » adapté (ajustement des barres de réglage et de la concentration enbore), avant d'atteindre les limites de fonctionnement du réacteur nécessitant une baisse de puissance, et éventuellement, son arrêt automatique par action des systèmes de protection.

Cette instabilité ne peut apparaître que dans les grands réacteurs, la diffusion des neutrons étant suffisante pour gommer les écarts de concentration à petite échelle.Les petits réacteurs non soumis à ces oscillations xénon sont dits :réacteurs couplés[réf. nécessaire].

Quelques équations

[modifier |modifier le code]

On donne les équations du phénomène de l'empoisonnement au xénon et les ordres de grandeurs dans le cas d'un réacteur de caractéristiques proches d'unREP 900 MWe.

On évalue en un premier temps les quantités d'iode et de xénon formées en fonction du flux (donc de la puissance du réacteur), on apprécie ensuite l'effet négatif sur la réactivité du cœur apporté par le xénon.

Les données et notations sont celles de l'articleflux neutronique

Données et notations principales
Grandeur physiqueNotationValeurUnitéGrandeur physiqueNotationValeurUnité
Puissance thermique
du cœur
W{\displaystyle W}2 768MWTempérature moyenne de
l'eau primaire dans le cœur
Tm304,5°C
Masse d’uraniumMU72 367kgFlux neutronique
thermique
ϕ{\displaystyle \phi }n/cm2/s
Enrichissement moyen
en uranium 235
2,433 %sans dimFlux neutronique[Note 2]
thermique à l'équilibre
ϕo{\displaystyle \phi _{o}}3,2 × 1013n/cm2/s
Volume du cœur26,571m3Concentration en iode 135I{\displaystyle I}at/cm3
Concentration en iode 135
à l'équilibre
Io{\displaystyle I_{o}}at/cm3
Concentration en xénon 135X{\displaystyle X}at/cm3
Concentration en xénon 135
à l'équilibre
Xo{\displaystyle X_{o}}at/cm3
Autres données et notations
Grandeur physiqueNotationValeurUnitéGrandeur physiqueNotationValeurUnité
Rendement de fission
de l'iode 135
γi{\displaystyle \gamma _{i}}0,064sans dimMasse d’uranium 235
= 72 367 × 2,433 %
MU51 760,93kg
Rendement de fission
du xénon 135
γx{\displaystyle \gamma _{x}}0,004sans dimMasse d'eau primaire
dans le cœur
11 073,8kg
Période
de l'iode 135

= 6,58 × 3 600
24 1710sConcentration des atomes
d'uranium 235
dans le cœur
U5{\displaystyle U_{5}}at/cm3
Constante radioactive
de l'iode 135
= ln(2)/(6,58 × 3 600)
λi{\displaystyle \lambda _{i}}2,926 × 10−5s−1Concentration des atomes
d'hydrogène dans le cœur
H{\displaystyle H}at/cm3
Période
du xénon 135
= 9,17 × 3 600
33 012s Masse de zirconium
dans le cœur
19 096kg
Constante radioactive
du xénon 135 = ln(2)/ 33 012
λx{\displaystyle \lambda _{x}}2,099 6 × 10−5s−1Concentration des atomes
d'oxygène du modérateur
dans le cœur
Om{\displaystyle O_{m}}at/cm3
Section efficace microscopique
de capture d'un neutron
par l'iode 135 = 7 barn
σi{\displaystyle \sigma _{i}}7 × 10−24cm2Facteur de fissions rapidesε{\displaystyle \varepsilon }1,07sans dim
Section efficace microscopique
de capture d'un neutron
par le xénon 135 = 2,65 × 106 barn
σx{\displaystyle \sigma _{x}}2,65 × 10−18cm2Section efficace microscopique
de fission thermique
de l'uranium 235 = 579,5 barn
σf{\displaystyle \sigma _{f}}5,795 × 10−22cm2
Section efficace microscopique
de capture d'un neutron
par le xénon 135 = 2,65 × 106 barn
σx{\displaystyle \sigma _{x}}2,65 × 10−18cm2Section efficace macroscopique
de fission thermique
Σf{\displaystyle \Sigma _{f}}0,098 40cm−1
Section efficace macroscopique
de capture d'un neutron
par le xénon 135
Σx{\displaystyle \Sigma _{x}}cm−1Section efficace macroscopique
d'absorption d'un neutron thermique
dans le combustible
Σu{\displaystyle \Sigma _{u}}cm−1
Section efficace microscopique
de fission de l'uranium 235 =
σf{\displaystyle \sigma _{f}}579,5barn Concentration des atomes
de zirconium dans le cœur
Z{\displaystyle Z}at/cm3
Section efficace microscopique
d'absorption[Note 3] par l'uranium 235 =
σa{\displaystyle \sigma _{a}}679,9barn Facteur d'utilisation thermique
avant empoisonnement xénon
f{\displaystyle f}sans dim
Section efficace microscopique
de capture d'un neutron thermique
par l'oxygène = 0,267 × 10−3 barn
σo{\displaystyle \sigma _{o}}0,267 × 10−27cm2 Facteur d'utilisation thermique
après empoisonnement xénon
f{\displaystyle f'}sans dim
Section efficace microscopique
de capture d'un neutron thermique
par le zirconium = 0,182 barn
σz{\displaystyle \sigma _{z}}0,182 × 10−24cm2Réactivité du cœur avant
empoisonnement xénon
ρ{\displaystyle \rho }pcm
Section efficace macroscopique
de capture d'un neutron thermique
dans le modérateur
Σm{\displaystyle \Sigma _{m}}cm−1Réactivité du cœur après
empoisonnement xénon
ρ{\displaystyle \rho '}pcm
Section efficace microscopique
de capture d'un neutron thermique
par l'hydrogène = 0,332 barn
σh{\displaystyle \sigma _{h}}0,332 × 10−24cm2Concentration des atomes
d'oxygène du combustible
dans le cœur
Ou{\displaystyle O_{u}}at/cm3
Section efficace microscopique
de capture d'un neutron thermique
par l'uranium 238 = 2,72 barn
σ8{\displaystyle \sigma _{8}}2,72 × 10−24cm2

Évolution de la concentration en iode 135

[modifier |modifier le code]
  • Lors d'un démarrage et d'une montée en puissance de quelques dizaines de minutes[Note 4], suivie d'un palier, à partir d'une situation cœur vierge sans iode ni xénon, l'iode atteint sa valeur d'équilibre au bout d'un temps de l'ordre de :1λi{\displaystyle {1 \over \lambda _{i}}} = 9,5 h
  • La concentration à l'équilibre s'écrit :
Io=γi×Σfλi+σi×Φo×Φoγi×Σfλi×Φo=6,886×1015{\displaystyle I_{o}={\gamma _{i}\times \Sigma _{f} \over \lambda _{i}+\sigma _{i}\times \Phi _{o}}\times \Phi _{o}\approxeq {\gamma _{i}\times \Sigma _{f} \over \lambda _{i}}\times \Phi _{o}=6{,}886\times 10^{15}}at/cm3
  • On peut voir que la concentration à l'équilibre est proportionnelle à la valeur du flux ; plus le flux est élevé plus il y a d'iode.
  • La fonction générale donnant I(t) à partir d'un démarrage sans iode ni xénon est la suivante:
I(t)=Io×(1eλi×t){\displaystyle I(t)=I_{o}\times (1-e^{-\lambda _{i}\times t})}[Note 5]
I(t)=γiλi×Σf×Φo×(1eλi×t){\displaystyle I(t)={\gamma _{i} \over \lambda _{i}}\times \Sigma _{f}\times \Phi _{o}\times (1-e^{-\lambda _{i}\times t})}
Si le flux est variable, on remplace Φo par Φ(t)
I(t)=γiλi×Σf×Φ(t)×(1eλi×t){\displaystyle I(t)={\gamma _{i} \over \lambda _{i}}\times \Sigma _{f}\times \Phi (t)\times (1-e^{-\lambda _{i}\times t})}
Évaluation de la quantité d'iode 135 dans le cœur
  • Liminaire : On néglige la présence du tellure 135 dont la période est faible devant celle de l'iode 135 et qui n'est pas capturant des neutrons et on fait l'hypothèse simplificatrice que l'iode 135 est produite directement par les fissions avec le rendement de fission du tellure 135 soit : 0,064
  • Variation de la concentration en iode 135 = dI/dt

dIdt=[Production par fission] - [Disparition par décroissance radioactive] - [Disparition par capture d'un neutron]{\displaystyle {\mathrm {d} I \over \mathrm {d} t}={\text{[Production par fission] - [Disparition par décroissance radioactive] - [Disparition par capture d'un neutron]}}}

dIdt=γi×Σf×Φ(t)λi×Iσi×I×Φ(t){\displaystyle {\mathrm {d} I \over \mathrm {d} t}=\gamma _{i}\times \Sigma _{f}\times \Phi (t)-\lambda _{i}\times I-\sigma _{i}\times I\times \Phi (t)}

0=γi×Σf×Φoλi×Iσi×I×Φo{\displaystyle 0=\gamma _{i}\times \Sigma _{f}\times \Phi _{o}-\lambda _{i}\times I-\sigma _{i}\times I\times \Phi _{o}}

Io=γi×Σfλi+σi×Φo×Φo{\displaystyle I_{o}={\gamma _{i}\times \Sigma _{f} \over \lambda _{i}+\sigma _{i}\times \Phi _{o}}\times \Phi _{o}}
  • On simplifie l'équation différentielle en remarquant que le terme σi × Φ est négligeable devant λi, ce qui revient à dire que l'iode n'est pas un poison, en effet :

σi×Φo=7×1024×3,2×1013=2,24×1010{\displaystyle \sigma _{i}\times \Phi _{o}=7\times 10^{-24}\times 3{,}2\times 10^{13}=2{,}24\times 10^{-10}} captures par seconde à comparer à :

λi=2,926×1005{\displaystyle \lambda _{i}=2{,}926\times 10^{-05}} désintégrations par seconde

d'où l'équation simplifiée :

dIdt=γi×Σf×Φλi×I{\displaystyle {\mathrm {d} I \over \mathrm {d} t}=\gamma _{i}\times \Sigma _{f}\times \Phi -\lambda _{i}\times I}

et la valeur à l'équilibre approximée :

γi×Σfλi×Φo=0,064×0,09840×3,2×10132,926×105=6,886×1015{\displaystyle {\gamma _{i}\times \Sigma _{f} \over \lambda _{i}}\times \Phi _{o}={0{,}064\times 0{,}098\,40\times 3{,}2\times 10^{13} \over 2{,}926\times 10^{-5}}=6{,}886\times 10^{15}}at/cm3[Note 6]

dIdt×eλi×t+λi×I×eλi×t=γi×Σf×Φ(t)×eλi×t{\displaystyle {\mathrm {d} I \over \mathrm {d} t}\times e^{\lambda _{i}\times t}+\lambda _{i}\times I\times e^{\lambda _{i}\times t}=\gamma _{i}\times \Sigma _{f}\times \Phi (t)\times e^{\lambda _{i}\times t}}

d(I×eλi×t)dt=γi×Σf×Φ(t)×eλi×t{\displaystyle {\mathrm {d} (I\times e^{\lambda _{i}\times t}) \over \mathrm {d} t}=\gamma _{i}\times \Sigma _{f}\times \Phi (t)\times e^{\lambda _{i}\times t}}

I(t)=I(o)+(γi×Σf×otΦ(u)×eλi×u×du)×eλi×t{\displaystyle I(t)=I(o)+\left(\gamma _{i}\times \Sigma _{f}\times \int _{o}^{t}\Phi (u)\times e^{\lambda _{i}\times u}\times \mathrm {d} u\right)\times e^{-\lambda _{i}\times t}}

  • On fait l'hypothèse d'un démarrage du réacteur vierge (sans iodes ni xénon), d'où I(o) = 0, et d'une croissance du flux pour atteindre la valeur Φo en quelques dizaines de minutes durée faible devant le rythme d'évolution de l'iode. On peut alors simplifier la fonction donnant I(t) en remplaçant Φ (t) par la valeur à terme du flux Φo[Note 7]

I(t)=γi×Σf×Φo×[oteλi×u×du]×eλi×t{\displaystyle I(t)=\gamma _{i}\times \Sigma _{f}\times \Phi _{o}\times \left[\int _{o}^{t}e^{\lambda _{i}\times u}\times \mathrm {d} u\right]\times e^{-\lambda _{i}\times t}}

I(t)=γi×Σf×Φo×eλi×t1λi×eλi×t{\displaystyle I(t)=\gamma _{i}\times \Sigma _{f}\times \Phi _{o}\times {e^{\lambda _{i}\times t}-1 \over \lambda _{i}}\times e^{-\lambda _{i}\times t}}

I(t)=γi×Σfλi×Φo×(1eλi×t)=Io×(1eλi×t){\displaystyle I(t)={\gamma _{i}\times \Sigma _{f} \over \lambda _{i}}\times \Phi _{o}\times (1-e^{-\lambda _{i}\times t})=I_{o}\times (1-e^{-\lambda _{i}\times t})}

 

Évolution de la concentration en xénon 135

[modifier |modifier le code]
  • Lors d'un démarrage et d'une montée en puissance de quelques dizaines de minutes suivie d'un palier, à partir d'une situation cœur vierge sans iode ni xénon, le xénon atteint sa valeur d'équilibre au bout d'un temps supérieur à celui de l'iode ; dit de façon simpliste, il faut attendre que l'iode ait atteint un régime stable pour que le xénon soit à l'équilibre.
  • La concentration à l'équilibre s'écrit :
Xo=(γx+γi)×Σf×Φoλx+σx×Φo=2,024×1015{\displaystyle X_{o}={(\gamma _{x}+\gamma _{i})\times \Sigma _{f}\times \Phi _{o} \over \lambda _{x}+\sigma _{x}\times \Phi _{o}}=2{,}024\times 10^{15}}at/cm3
  • On peut voir que si le flux augmente la concentration en xénon à l'équilibre plafonne vers la valeur :

Xomax=γx+γi×Σfσx=2,525×1015{\displaystyle X_{omax}={\frac {\gamma _{x}+\gamma _{i}\times \Sigma _{f}}{\sigma _{x}}}=2{,}525\times 10^{15}}at/cm3

qui ne dépend pas du flux. Ce résultat est à rapprocher de celui concernant l'iode dont la concentration est proportionnelle au flux.

Dans le cas d'un réacteur à flux élevé il y a une réserve d'iode d'autant plus importante en valeur relative devant le xénon présent dans le cœur à l'équilibre que le flux est élevé. Le "pic de xénon" suivant un arrêt du réacteur en sera d'autant plus important.

IoXo=γi×Σf×Φoλi(γi+γx)×Σf×Φoλx+σx×Φo=γiγi+γx×λx+σx×Φoλi{\displaystyle {I_{o} \over X_{o}}={{\gamma _{i}\times \Sigma _{f}\times \Phi _{o} \over \lambda _{i}} \over {(\gamma _{i}+\gamma _{x})\times \Sigma _{f}\times \Phi _{o} \over \lambda _{x}+\sigma _{x}\times \Phi _{o}}}={\gamma _{i} \over \gamma _{i}+\gamma _{x}}\times {\lambda _{x}+\sigma _{x}\times \Phi _{o} \over \lambda _{i}}}

Le rapport le rapportIoXo{\displaystyle {I_{o} \over X_{o}}} augmente avec la valeur du flux.. γi γx σx λx et λi .{\displaystyle .~\gamma _{i}~\gamma _{x}~\sigma _{x}~\lambda _{x}~et~\lambda _{i}~.} sont des constantes physiques, le rapportIoXo{\displaystyle {I_{o} \over X_{o}}} dépend du dessin du réacteur considéré uniquement par la valeur du flux lequel est inversement proportionnel à la concentration moyenne de matière fissile présente dans le cœur. Plus l'usure du combustible est importante plus la concentration du xénon à l'équilibre est faible en valeur relative devant la quantité d'iode "stockée" conduisant à un effet de "pic" accru.

La fonction générale donnant X(t) en cas de démarrage cœur vierge sans iode ni xénon est assez lourde :

X(t)=Xo×[1+1γi+γx×(γi×λiλxλi+σx×Φoγx)×e(λx+σx×Φo)×tγiγi+γx×(λx+σx×Φoλxλi+σx×Φo)×eλi×t]{\displaystyle X(t)=X_{o}\times \left[1+{1 \over \gamma _{i}+\gamma _{x}}\times \left({\gamma _{i}\times \lambda _{i} \over \lambda _{x}-\lambda _{i}+\sigma _{x}\times \Phi _{o}}-\gamma _{x}\right)\times e^{-(\lambda _{x}+\sigma _{x}\times \Phi _{o})\times t}-{\gamma _{i} \over \gamma _{i}+\gamma _{x}}\times \left({\lambda _{x}+\sigma _{x}\times \Phi _{o} \over \lambda _{x}-\lambda _{i}+\sigma _{x}\times \Phi _{o}}\right)\times e^{-\lambda _{i}\times t}\right]}
avec Xo = valeur à l'équilibre

I(t)=γi×Σfλi×Φo×eλi×t{\displaystyle I(t)={\gamma _{i}\times \Sigma _{f} \over \lambda _{i}}\times \Phi _{o}\times e^{-\lambda _{i}\times t}}

X(t)=Σf×Φo×[γiλxλi×eλi×t+(γi+γxλx+σx×Φoγiλxλi)×eλx×t]{\displaystyle X(t)=\Sigma _{f}\times \Phi _{o}\times \left[{\gamma _{i} \over \lambda _{x}-\lambda _{i}}\times e^{-\lambda _{i}\times t}+\left({\gamma _{i}+\gamma _{x} \over \lambda _{x}+\sigma _{x}\times \Phi _{o}}-{\gamma _{i} \over \lambda _{x}-\lambda _{i}}\right)\times e^{-\lambda _{x}\times t}\right]}

On en déduit l'instant du « pic » suivant un arrêt rapide après une marche continue à pleine puissance :

tpic=1λxλi×ln[λx×(λxλi)λi×γi×(γiλxλiγi+γxλx+σx×Φo)]{\displaystyle t_{pic}={1 \over \lambda _{x}-\lambda _{i}}\times \ln \left[{\lambda _{x}\times (\lambda _{x}-\lambda _{i}) \over \lambda _{i}\times \gamma _{i}}\times \left({\gamma _{i} \over \lambda _{x}-\lambda _{i}}-{\gamma _{i}+\gamma _{x} \over \lambda _{x}+\sigma _{x}\times \Phi _{o}}\right)\right]}

Dans l'exemple type choisi on trouvetpic = 3,051 × 104 s = 8,47 h. Plus le flux est élevé plus le pic est tardif et ample.

Évaluation de la concentration du xénon 135 dans le cœur

Variation de la concentration en xénon 135 =dXdt{\displaystyle {\frac {\mathrm {d} X}{\mathrm {d} t}}}

dXdt=+[Production par decroissance de liode]+[Production par fission][Disparition par decroissance radioactive][Disparition par capture d un neutron]{\displaystyle {\mathrm {d} X \over \mathrm {d} t}=\scriptstyle +[{\mathsf {Production~par~decroissance~de~l'iode}}]+[{\mathsf {Production~par~fission}}]-[{\mathsf {Disparition~par~decroissance~radioactive}}]-[{\mathsf {Disparition~par~capture~d'~un~neutron}}]}

dXdt=γx×Σf×Φ(t)+λi×Iλx×Xσx×X×Φ(t){\displaystyle {\mathrm {d} X \over \mathrm {d} t}=\gamma _{x}\times \Sigma _{f}\times \Phi (t)+\lambda _{i}\times I-\lambda _{x}\times X-\sigma _{x}\times X\times \Phi (t)}

0=γx×Σf×Φo+λi×Ioσx×Xo×Φoλx×Xo{\displaystyle 0=\gamma _{x}\times \Sigma _{f}\times \Phi _{o}+\lambda _{i}\times I_{o}-\sigma _{x}\times X_{o}\times \Phi _{o}-\lambda _{x}\times X_{o}}

Xo=γx×Σf×Φo+λi×Ioλx+σx×Φo{\displaystyle X_{o}={\gamma _{x}\times \Sigma _{f}\times \Phi _{o}+\lambda _{i}\times I_{o} \over \lambda _{x}+\sigma _{x}\times \Phi _{o}}}

Io=γiλi×Σf×Φo{\displaystyle I_{o}={\gamma _{i} \over \lambda _{i}}\times \Sigma _{f}\times \Phi _{o}}

Xo=γx×Σf×Φo+λi×(γiλi×Σf×Φo)λx+σx×Φo=(γx+γi)×Σf×Φoλx+σx×Φo{\displaystyle X_{o}={\gamma _{x}\times \Sigma _{f}\times \Phi _{o}+\lambda _{i}\times ({\gamma _{i} \over \lambda _{i}}\times \Sigma _{f}\times \Phi _{o}) \over \lambda _{x}+\sigma _{x}\times \Phi _{o}}={(\gamma _{x}+\gamma _{i})\times \Sigma _{f}\times \Phi _{o} \over \lambda _{x}+\sigma _{x}\times \Phi _{o}}}

Xo=(0,004+0,064)×0,09840×3,2×1013 2,0997×105+2,65×1018×3,2×1013=2,024×1015{\displaystyle X_{o}={(0{,}004+0{,}064)\times 0{,}098\,40\times 3{,}2\times 10^{13} \over \ 2{,}0997\times 10^{-5}+2{,}65\times 10^{-18}\times 3{,}2\times 10^{13}}=2{,}024\times 10^{15}}at/cm3

Xomax=(0,064+0,004)×0,098402,65×1018at/cm3=2,525×1015{\displaystyle Xomax={(0{,}064+0{,}004)\times 0{,}098\,40 \over 2{,}65\times 10^{-18}}\;{\rm {at/cm^{3}=2{,}525\times 10^{15}}}}

La valeur limite Xomax ne dépend pas du flux.

Production de xénon par désintégration de l'iode =λi×Io=2,926×105×6,886×1015=2,0148×1011{\displaystyle \lambda _{i}\times I_{o}=2{,}926\times 10^{-5}\times 6{,}886\times 10^{15}=2{,}0148\times 10^{11}}at/cm3/s ; soit 94,1 % du total produit

Production de xénon directement par fission =γx×Σf×Φo=0,04×0,09840×3,2×1013=1,2593×1010{\displaystyle \gamma _{x}\times \Sigma _{f}\times \Phi _{o}=0{,}04\times 0{,}098\,40\times 3{,}2\times 10^{13}=1{,}2593\times 10^{10}}at/cm3/s ; soit 5,9 % du total produit[Note 8]

avec

ψ(u)=λx+σx×Φ(u){\displaystyle \psi (u)=\lambda _{x}+\sigma _{x}\times \Phi (u)}

Il vient :

dXdt×eotψ(u)×du+ψ(t)×X(t)×eotψ(u)×du=[λi×I(t)+γx×Σf×Φ(t)]×eotψ(u)×du{\displaystyle {\mathrm {d} X \over \mathrm {d} t}\times e^{\int _{o}^{t}\psi (u)\times \mathrm {d} u}+\psi (t)\times X(t)\times e^{\int _{o}^{t}\psi (u)\times \mathrm {d} u}=\left[\lambda _{i}\times I(t)+\gamma _{x}\times \Sigma _{f}\times \Phi (t)\right]\times e^{\int _{o}^{t}\psi (u)\times \mathrm {d} u}}

d(X×eotψ(u)×du)dt=[λi×I(t)+γx×Σf×Φ(t)]×eotψ(u)×du{\displaystyle {\mathrm {d} (X\times e^{\int _{o}^{t}\psi (u)\times \mathrm {d} u}) \over \mathrm {d} t}=\left[\lambda _{i}\times I(t)+\gamma _{x}\times \Sigma _{f}\times \Phi (t)\right]\times e^{\int _{o}^{t}\psi (u)\times \mathrm {d} u}}

  • Dans l'hypothèse d'un démarrage du réacteur vierge (sans iode ni xénon), d'où I(o) = X(o) = 0, et d'une croissance du flux pour atteindre la valeur Φo en quelques dizaines de minutes durée faible devant l'évolution de l'iode et du xénon, on peut simplifier la fonction donnant X(t) en remplaçant Φ (t) par la valeur à terme Φo ;

Par ailleurs, on a vu que :

I(t)=γi×Σfλi×Φo×(1eλi×t){\displaystyle I(t)={\gamma _{i}\times \Sigma _{f} \over \lambda _{i}}\times \Phi _{o}\times (1-e^{-\lambda _{i}\times t})}X(t)=Xo×[1+1γi+γx×(γi×λiλxλi+σx×Φoγx)×e(λx+σx×Φo)×tγiγi+γx×(λx+σx×Φoλxλi+σx×Φo)×eλi×t]{\displaystyle X(t)=X_{o}\times \left[1+{1 \over \gamma _{i}+\gamma _{x}}\times \left({\gamma _{i}\times \lambda _{i} \over \lambda _{x}-\lambda _{i}+\sigma _{x}\times \Phi _{o}}-\gamma _{x}\right)\times e^{-(\lambda _{x}+\sigma _{x}\times \Phi _{o})\times t}-{\gamma _{i} \over \gamma _{i}+\gamma _{x}}\times \left({\lambda _{x}+\sigma _{x}\times \Phi _{o} \over \lambda _{x}-\lambda _{i}+\sigma _{x}\times \Phi _{o}}\right)\times e^{-\lambda _{i}\times t}\right]}

avec Xo = valeur à l'équilibre

  • Dans le cas d'un arrêt rapide du réacteur à partir d'un régime permanent initial établi les équations s'allègent un peu :

I(t)=γi×Σfλi×Φo×eλi×t{\displaystyle I(t)={\gamma _{i}\times \Sigma _{f} \over \lambda _{i}}\times \Phi _{o}\times e^{-\lambda _{i}\times t}}

X(t)=Σf×Φo×[γiλxλi×eλi×t+(γi+γxλx+σx×Φoγiλxλi)×eλx×t]{\displaystyle X(t)=\Sigma _{f}\times \Phi _{o}\times \left[{\gamma _{i} \over \lambda _{x}-\lambda _{i}}\times e^{-\lambda _{i}\times t}+\left({\gamma _{i}+\gamma _{x} \over \lambda _{x}+\sigma _{x}\times \Phi _{o}}-{\gamma _{i} \over \lambda _{x}-\lambda _{i}}\right)\times e^{-\lambda _{x}\times t}\right]}

La fonction donnant X(t) est aisément dérivable :

X(t)÷λi×γiλxλi×eλi×t+λx×(γiλxλiγi+γxλx+σx×Φo)×eλx×t{\displaystyle X'(t)\div -\lambda _{i}\times {\gamma _{i} \over \lambda _{x}-\lambda _{i}}\times e^{-\lambda _{i}\times t}+\lambda _{x}\times \left({\gamma _{i} \over \lambda _{x}-\lambda _{i}}-{\gamma _{i}+\gamma _{x} \over \lambda _{x}+\sigma _{x}\times \Phi _{o}}\right)\times e^{-\lambda _{x}\times t}}

L'instant du "pic xénon" est évaluable analytiquement.

λi×γiλxλi×e(λxλi)×tpic=λx×(γiλxλiγi+γxλx+σx×Φo){\displaystyle \lambda _{i}\times {\gamma _{i} \over \lambda _{x}-\lambda _{i}}\times e^{(\lambda _{x}-\lambda _{i})\times t_{pic}}=\lambda _{x}\times \left({\gamma _{i} \over \lambda _{x}-\lambda _{i}}-{\gamma _{i}+\gamma _{x} \over \lambda _{x}+\sigma _{x}\times \Phi _{o}}\right)}

Pour finir :

tpic=1λxλi×ln[λx×(λxλi)λi×γi×(γiλxλiγi+γxλx+σx×Φo)]{\displaystyle t_{pic}={1 \over \lambda _{x}-\lambda _{i}}\times \ln \left[{\lambda _{x}\times (\lambda _{x}-\lambda _{i}) \over \lambda _{i}\times \gamma _{i}}\times \left({\gamma _{i} \over \lambda _{x}-\lambda _{i}}-{\gamma _{i}+\gamma _{x} \over \lambda _{x}+\sigma _{x}\times \Phi _{o}}\right)\right]}

Dans l'exemple type choisi du REP 900 MWe on trouve tpic = 3,051 × 104 s = 8,47 h post arrêt. Plus le flux est élevé plus le pic est tardif et ample.

 

Empoisonnement xénon - Effet en réactivité

[modifier |modifier le code]

Cœur vierge avant empoisonnement :f=ΣuΣu+Σm{\displaystyle f={\frac {\Sigma _{u}}{\Sigma _{u}+\Sigma _{m}}}}[Note 9]

Après empoisonnementf=ΣuΣu+Σm+Σx{\displaystyle f'={\frac {\Sigma _{u}}{\Sigma _{u}+\Sigma _{m}+\Sigma _{x}}}}

ff=0,96406=keffkeff=eρeρ=eρρ{\displaystyle {f' \over f}=0{,}964\,06={k_{eff}' \over k_{eff}}={e^{\rho '} \over e^{\rho }}=e^{\rho '-\rho }}

ρρ=100000×ln(0,96406)={\displaystyle \rho '-\rho =100\,000\times \ln(0{,}964\,06)=} -3 660 pcm

Évaluation de l'effet en réactivité du xénon 135
  • L'empoisonnement par le xénon affecte le facteurf de la formule des 4 facteurs

f avant empoisonnement =f=ΣuΣu+Σm{\displaystyle f={\Sigma _{u} \over \Sigma _{u}+\Sigma _{m}}}

f' après empoisonnement =f=ΣuΣu+Σm+Σx{\displaystyle f'={\Sigma _{u} \over \Sigma _{u}+\Sigma _{m}+\Sigma _{x}}}

Le rapportff{\displaystyle {\frac {f}{f'}}} exprime la variation dekeff

ff=1+ΣxΣu+Σm{\displaystyle {f \over f'}=1+{\Sigma _{x} \over \Sigma _{u}+\Sigma _{m}}}

Σu=(679,9×U5+2,72×U8+0,267×103×Ou+0,182×Zr)×1024{\displaystyle \Sigma _{u}=(679{,}9\times U_{5}+2{,}72\times U_{8}+0{,}267\times 10^{-3}\times O_{u}+0{,}182\times Zr)\times 10^{-24}}cm−1

Concentration des atomes d'uranium 235 =U5=1760,93/0,235044×NA/26,571×106=1,6980×1020{\displaystyle U_{5}=1760{,}93/0{,}235\,044\times N_{A}/26{,}571\times 10^{6}=1{,}6980\times 10^{20}}at/cm3

Concentration des atomes d'uranium 238 =U8=723671760,93)0,238051×NA/26,571×106=6,722×1021{\displaystyle U_{8}={\frac {72\,367-1760{,}93)}{0{,}238\,051\times N_{A}/26{,}571\times 10^{6}}}=6{,}722\times 10^{21}}at/cm3

Concentration des atomes d'oxygène du combustible =Ou=2×(U5+U8)=1,3784×1022{\displaystyle O_{u}=2\times (U_{5}+U_{8})=1{,}3784\times 10^{22}}at/cm3

Concentration des atomes de zirconium =Zr=19096/0,091224×NA/26,571×106=4,744×1021{\displaystyle Zr=19\,096/0{,}091\,224\times N_{A}/26{,}571\times 10^{6}=4{,}744\times 10^{21}}at/cm3

D'où :Σu=0,13460{\displaystyle \Sigma _{u}=0{,}134\,60}cm−1

Σm=σh×H+σo×Om{\displaystyle \Sigma _{m}=\sigma _{h}\times H+\sigma _{o}\times O_{m}}

Σm=(0,332×H+0,267×103×Om)×1024{\displaystyle \Sigma _{m}=(0{,}332\times H+0{,}267\times 10^{-3}\times O_{m})\times 10^{-24}}cm−1

Nombre de moles d'eau dans le cœur =11073,8×100015,9994+2×1,00794=614690{\displaystyle {\frac {11\,073{,}8\times 1000}{15{,}9994+2\times 1{,}007\,94}}=614\,690}

Concentration des atomes d'hydrogène =614690×2×NA26,571×106=2,786×1022{\displaystyle {\frac {614\,690\times 2\times N_{A}}{26{,}571\times 10^{6}}}=2{,}786\times 10^{22}}at/cm3

Concentration des atomes d'oxygène du modérateur =2,786×1022/2=1,3932×1022{\displaystyle 2{,}786\times 10^{22}/2=1{,}3932\times 10^{22}}at/cm3

Σm=2,786×1022×0,332+1,3932×1022×0,267×103×1024=0,0092543{\displaystyle \Sigma _{m}=2{,}786\times 10^{22}\times 0{,}332+1{,}3932\times 10^{22}\times 0{,}267\times 10^{-3}\times 10^{-24}=0{,}009\,2543}cm−1

ff=1+0,00536300,13460+0,0092543=1,03728{\displaystyle {f \over f'}=1+{0{,}005\,3630 \over 0{,}134\,60+0{,}009\,2543}=1{,}037\,28}

La variation relative de keff, si le cœur était juste critique avant empoisonnement, s'écrit:

keffkeff=ff=11,03728=0,96406{\displaystyle {k_{eff}' \over k_{eff}}={f' \over f}={\frac {1}{1{,}037\,28}}=0{,}964\,06}

L'effet en réactivité correspondant vaut :keffkeff=eρeρ=eρρ{\displaystyle {k_{eff}' \over k_{eff}}={e^{\rho '} \over e^{\rho }}=e^{\rho '-\rho }}

exprimé en pcm
ρρ=100000×ln(0,96406)={\displaystyle \rho '-\rho =100\,000\times \ln(0{,}964\,06)=} -3 660 pcm

ff=1+σx×(γi+γx)×Σf×Φ(t)λx+σx×Φ(t)Σu+Σm{\displaystyle {f \over f'}=1+{\sigma _{x}\times {(\gamma _{i}+\gamma _{x})\times \Sigma _{f}\times \Phi (t) \over \lambda _{x}+\sigma _{x}\times \Phi (t)} \over \Sigma _{u}+\Sigma _{m}}}

ff=1+1,2326×1019×Φ(t)2,0997×105+2,65×1018×Φ(t){\displaystyle {f \over f'}=1+{1{,}2326\times 10^{-19}\times \Phi (t) \over 2{,}0997\times 10^{-5}+2{,}65\times 10^{-18}\times \Phi (t)}}

 

Références et liens

[modifier |modifier le code]
  1. a etbHandbook of Chemistry and Physics 58th edition
  2. (en)Chain Fission Yields, AIEA.
  3. Fission product data: Thermal neutron cross sections, resonance integrals, and Westcott factors, AIEA.
  4. a etbThe Influence of Xenon-135 on Reactor Operation, Paul L. Roggenkamp.
  5. Pic xénon transitoire, banque de données terminologiques et linguistiques du gouvernement du Canada, consulté le 18 avril 2022
  6. « Tchernobyl, 25 ans après : un tel scénario peut-il se répéter ?» réponses de Michel Chouha, physicien, représentant de l'IRSN en Europe de l'Est,Sciences et avenir, 26 avril 2011.
  7. DOE Handbook, pp. 35–42.
  8. CONTRÔLE ET PILOTAGE DES RÉACTEURS À EAU SOUS-PRESSION - Clefs CEA N°45 Automne 2001 - pp 39,40 - Anne Nicolas (Direction de l'énergie nucléaire CEA/Saclay)
  9. Le pilotage des réacteurs - Réunion SFEN du 15 mars 2007 - François BOUTEILLE ST6 « Physique des Réacteurs (AREVA) »
  10. La physique des réacteurs nucléaires (Coll EDF R&D) - p1058 & suivantes - Auteur : MARGUET Serge - Date de parution : 04-2011 - Éditions Lavoisier
  11. ab etcReactor Theory (Nuclear Parameters) DOE-HDBK-1019/2-93 XENON

Notes

[modifier |modifier le code]
  1. Le produit : (période de décroissance du prométhium 149 = 53,08 h + 1,728 h) × (rendement de fission = 1,09 %) est supérieur au même produit pour l'iode 135 (6,58 h × 6,4 %) donc la quantité de prométhium à l'équilibre est plus élevée.
  2. En toute rigueur, la valeur du flux thermique donnée dans l'article flux neutronique serait à augmenter de 7 % pour tenir compte de ce que le flux thermique est évalué dans cet article en tenant compte d'un flux rapide responsable de 7 % des fissions
  3. On entend par absorption le total fission + capture
  4. Le rythme maximal de montée en allure d'un REP est de 6 % par minute ; le rythme normal est de 3 % par minute
  5. Relation similaire à une décroissance radioactive simple
  6. Les valeurs du flux et de Σf sont calculées dans l'articleFlux neutronique
  7. Ceci revient à considérer un "échelon" de puissance pour le réacteur
  8. La quantité de xénon produite par fission est un peu plus élevée car il faut tenir compte des fissions rapides
  9. L'incertitude principale du calcul effectué dans cet article vient de l'estimation de la capture dans le combustible : Σu
Ce document provient de « https://fr.wikipedia.org/w/index.php?title=Empoisonnement_au_xénon&oldid=223326852 ».
Catégorie :
Catégories cachées :

[8]ページ先頭

©2009-2025 Movatter.jp