Empoisonnement au samarium

Effet en réactivité d'un arrêt et d'un redémarrage au « pic » xénon + samarium. Le xénon est prépondérant

L'empoisonnement par le samarium 149 est une réduction de la réactivité d'un réacteur nucléaire, provoquée par la capture de neutrons par le samarium 149, un produit de fission stable qui est un poison nucléaire. Il s'agit d'un processus similaire à l'empoisonnement par le xénon 135.

La fission produit du néodyme 149 qui interagit peu avec les neutrons et qui se dégrade avec une demi-vie de 1,728 h en prométhium 149 qui se désintègre à son tour en samarium 149 avec une période de 53,08 h. Le samarium ainsi formé est un poison neutronique qui présente une section efficace de capture très élevée (41 140 barns).

Description qualitative

Arrêt et redémarrage au « pic » samarium

$\text{Fission}\underset{1,09\mathrm{\%}}{\overset{}{\to }}{}^{149}\mathrm{N}\mathrm{d}\underset{1,728h}{\overset{\beta }{\to }}{}^{\mathrm{1}\mathrm{4}\mathrm{9}}\mathrm{P}\mathrm{m}\underset{53,08h}{\overset{\beta }{\to }}{}^{\mathrm{1}\mathrm{4}\mathrm{9}}\mathrm{S}\mathrm{m}\underset{\sigma =40140b}{\overset{+n}{\to }}{}^{\mathrm{1}\mathrm{5}\mathrm{0}}\mathrm{S}\mathrm{m}$

En fonctionnement stable prolongé, production et dégradation de ces corps s'équilibrent et une concentration stable s'établit à terme :

• le néodyme 149 est produit par les fissions (donc proportionnellement au flux) avec un rendement de fission de 1,09 %; il disparaît par décroissance radioactive (émission d'un électron de façon à réduire son excès de neutrons) en donnant le prométhium 149
• le prométhium 149 est produit par la transmutation bêta moins du néodyme 149 ; il disparaît par radioactivité bêta moins en donnant le samarium 149
• le samarium 149 (stable) n'est produit que par la désintégration du prométhium 149 ; il n'est pas radioactif et ne disparaît que par capture neutronique sous flux.

En cas d'arrêt du réacteur :

• le néodyme 149 présent dans le cœur cesse d'être produit et se transforme assez rapidement en prométhium
• de même le prométhium se désintègre en samarium 149
• le samarium 149 qui est stable cesse d'être consommé par le flux et en une semaine, la quasi-totalité du néodyme 149 et du prométhium 149 s'est désintégrée en samarium 149 augmentant ainsi l'empoisonnement neutronique du réacteur.

Le processus conduisant à l'empoisonnement du cœur par le samarium 149 est ainsi largement similaire à l'empoisonnement par le xénon 135. Des différences existent cependant :

• Le samarium 149 est stable et ne disparaît que par capture neutronique sous flux ;
• Les vitesses d'évolution des concentrations sont beaucoup plus lentes que dans le cas de l'iode 135 et du xénon 135 ;
• Les quantités de prométhium 149 et de samarium 149 présentes dans le cœur à l'équilibre sont plus importantes que celles de l'iode 135 et le xénon 135. En effet, bien que le rendement de fission du tellure 135 soit supérieur à celui du néodyme 149, les périodes radioactives du néodyme 149 et du prométhium 149 sont plus élevées que celles du tellure 135 et de l'iode 135 ce qui compense et au-delà le moindre rendement de fission. (1,09 % × 53,08 pour le prométhium 149) contre (6,4 % × 6,58 pour l'iode 135). Les concentrations à l'équilibre sont respectivement dans le même rapport.
• Les effets en réactivité nucléaire sont moindres car la section efficace de capture du samarium 149 (40 140 barns) est très inférieure celle du xénon 135 (2 650 000 barns), ceci malgré les quantités plus importantes de samarium 149 et de prométhium 149 que de xénon 135 et d'iode 135 présentes dans le réacteur dans le cours du fonctionnement.

Quelques équations

On se propose de retrouver les ordres de grandeur de l'empoisonnement en samarium dans le cas d'un réacteur de caractéristiques proches d'un REP 900 MWe en résolvant les équations du phénomène. Les relations modélisant les concentrations et effet en réactivité sont formellement les mêmes que dans le cas du couple iode 135/xénon 135, elles sont toutefois un peu plus simples.

On évalue en un premier temps les quantités de prométhium et de samarium 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 samarium.

Les données et notations sont celles des articles flux neutronique et empoisonnement au xénon

Données et notations principales

Grandeur physique	Notation	Valeur	Unité	Grandeur physique	Notation	Valeur	Unité
Puissance thermique du cœur	W	2 785	MW	Température moyenne de l'eau primaire dans le cœur	Tm	304,5	°C
Masse d’uranium	MU	72 367	kg	Flux neutronique thermique	$\mathrm{\Phi }$		n/cmModèle:2/s
Enrichissement moyen en uranium 235		2,433 %	ss dim	Flux neutronique thermique à l'équilibre	Φo	[[Flux neutronique|Modèle:Nb]]	n/cmModèle:2/s
Volume du cœur		26,571	mModèle:3	Concentration en prométhium 149	P		Modèle:Nb
				Concentration en en prométhium 149 à l'équilibre =	$P_o=\frac{{\gamma }_p\times {\mathrm{\Sigma }}_f}{{\lambda }_p}\times {\mathrm{\Phi }}_o$	Modèle:Nb	Modèle:Nb
				Concentration en samarium 149	S		Modèle:Nb
				Concentration en samarium 149 à l'équilibre	So		Modèle:Nb

Autres données et notations

Grandeur physique	Notation	Valeur	Unité	Grandeur physique	Notation	Valeur	Unité
Rendement de fission du néodyme 149	${\gamma }_p$	0,010 9	ss dim	Masse d’uranium 235 = 72 367 × 2,433 %	MU5	1 760,93	kg
				Masse d'eau primaire dans le cœur		11 073,8	kg
Période du prométhium 149 = 53,08 h		191 088	s	Concentration des atomes d'uranium 235 dans le cœur	U5		Modèle:Nb
Constante radioactive du prométhium 149 = ln(2)/(53,08 × 3 600)	${\lambda }_p$	Modèle:Nb	sModèle:-1	Concentration des atomes d'hydrogène dans le cœur	H		Modèle:Nb
				Masse de zirconium dans le cœur		19 096	kg
				Concentration des atomes d'oxygène du modérateur dans le cœur	Om		Modèle:Nb
Section efficace microscopique de capture d'un neutron par le prométhium 149 = 1 350 barns	${\sigma }_p$	Modèle:Nb	cmModèle:2	Facteur de fissions rapides	$\epsilon$	1,07	ss dim
Section efficace microscopique de capture d'un neutron par le samarium 149	${\sigma }_s$	Modèle:Nb	cmModèle:2	Section efficace microscopique de fission thermique de l'uranium 235 = 579,5 barns	${\sigma }_f$	Modèle:Nb	cmModèle:2
		'	cmModèle:2	Section efficace macroscopique de fission thermique	${\mathrm{\Sigma }}_f$	0,098 40	cmModèle:-1
Section efficace macroscopique de capture d'un neutron par le samarium 149	${\mathrm{\Sigma }}_s$		cmModèle:-1	Section efficace macroscopique d'absorption d'un neutron thermique dans le combustible	${\mathrm{\Sigma }}_u$		cmModèle:-1
Section efficace microscopique de fission de l'uranium 235 =	${\sigma }_u$	579,5	barn	Concentration des atomes de zirconium dans le cœur	Zr		Modèle:Nb
Section efficace microscopique d'absorption par l'uranium 235 =	${\sigma }_a$	679,9	barn	Facteur d'utilisation thermique avant empoisonnement samarium	f		ss dim
Section efficace microscopique de capture d'un neutron thermique par l'oxygène = Modèle:Nb	${\sigma }_o$	Modèle:Nb	cmModèle:2	Facteur d'utilisation thermique après empoisonnement samarium	f'		ss dim
Section efficace microscopique de capture d'un neutron thermique par le zirconium = 0,182 barn	${\sigma }_z$	Modèle:Nb	cmModèle:2	Réactivité du cœur avant empoisonnement samarium 149	$\rho$		pcm
Section efficace macroscopique de capture d'un neutron thermique dans le modérateur	${\mathrm{\Sigma }}_m$		cmModèle:-1	Réactivité du cœur après empoisonnement samarium	${\rho }^{\prime }$		pcm
Section efficace microscopique de capture d'un neutron thermique par l'hydrogène = 0,332 barn	${\sigma }_h$	Modèle:Nb	cmModèle:2	Concentration des atomes d'oxygène du combustible dans le cœur	Ou		Modèle:Nb
Section efficace microscopique de capture d'un neutron thermique par l'uranium 238 = 2,72 barns	${\sigma }_8$	Modèle:Nb	cmModèle:2

Évolution de la quantité de prométhium

• 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 prométhium ni samarium, le prométhium atteint sa valeur d'équilibre au bout d'un temps de l'ordre de :

$\frac{p\acute{e}riode}{\ln (2)}=\frac{1}{{\lambda }_p}=76,6\mathrm{h}$

• La concentration à l'équilibre s'écrit :

$P_o=\frac{{\gamma }_p\times {\mathrm{\Sigma }}_f}{{\lambda }_p+{\sigma }_p\times {\mathrm{\Phi }}_o}\times {\mathrm{\Phi }}_o\approxeq \frac{{\gamma }_p\times {\mathrm{\Sigma }}_f}{{\lambda }_p}\times {\mathrm{\Phi }}_o=9,461\times 10^{15}\mathrm{a}\mathrm{t}/\mathrm{c}{\mathrm{m}}^{\mathrm{3}}$

• On peut voir que la concentration à l'équilibre est proportionnelle à la valeur du flux ; plus le flux est élevé plus il y a de prométhium.
• La fonction générale donnant P(t) est la suivante :

$P(t)=P_o\times (1-{\mathrm{e}}^{-{\lambda }_p\times t})$

$P(t)=\frac{{\gamma }_p}{{\lambda }_p}\times {\mathrm{\Sigma }}_f\times {\mathrm{\Phi }}_o\times (1-{\mathrm{e}}^{-{\lambda }_p\times t})$

Si le flux est variable, on remplace Φ_o par Φ(t)

$P(t)=\frac{{\gamma }_p}{{\lambda }_p}\times {\mathrm{\Sigma }}_f\times \mathrm{\Phi }(t)\times (1-{\mathrm{e}}^{-{\lambda }_p\times t})$

• À noter que la concentration à l'équilibre du prométhium 149 est supérieure à celle de l'iode 135, en effet :

Concentration à l'équilibre en iode 135 : $I_o=\frac{{\gamma }_i}{{\lambda }_i}\times {\mathrm{\Sigma }}_f\times {\mathrm{\Phi }}_o$

Concentration à l'équilibre en prométhium 149 : $P_o=\frac{{\gamma }_p}{{\lambda }_p}\times {\mathrm{\Sigma }}_f\times {\mathrm{\Phi }}_o$

Rapport : $\frac{P_o}{I_o}=\frac{{\gamma }_p}{{\gamma }_i}\times \frac{{\lambda }_i}{{\lambda }_p}=\frac{0,0109}{0,064}\times \frac{2,926\times 10^{-5}}{3,627\times 10^{-6}}=1,374$

Il y a 1,4 fois plus de prométhium 149 que d'iode 135 dans le cœur à l'équilibre

• Liminaire : On néglige la présence du néodyme 149 dont la période est petite devant celle du prométhium 149 et qui n'est pas capturant des neutrons et on fait l'hypothèse simplificatrice que le prométhium est produit directement par les fissions avec le rendement de fission du néodyme 149 soit : 0,010 9
• Variation de la concentration en prométhium 149 = dP/dt

$\frac{\mathrm{d}P}{\mathrm{d}t}=+[𝖯𝗋𝗈𝖽𝗎𝖼𝗍𝗂𝗈𝗇𝗉𝖺𝗋𝖿𝗂𝗌𝗌𝗂𝗈𝗇]-[𝖣𝗂𝗌𝗉𝖺𝗋𝗂𝗍𝗂𝗈𝗇𝗉𝖺𝗋𝖽𝖾𝖼𝗋𝗈𝗂𝗌𝗌𝖺𝗇𝖼𝖾𝗋𝖺𝖽𝗂𝗈𝖺𝖼𝗍𝗂𝗏𝖾]-[𝖣𝗂𝗌𝗉𝖺𝗋𝗂𝗍𝗂𝗈𝗇𝗉𝖺𝗋𝖼𝖺𝗉𝗍𝗎𝗋𝖾{𝖽}^{\prime }𝗎𝗇𝗇𝖾𝗎𝗍𝗋𝗈𝗇]$

$\frac{\mathrm{d}P}{\mathrm{d}t}={\gamma }_p\times {\mathrm{\Sigma }}_f\times \mathrm{\Phi }(t)-{\lambda }_p\times P-{\sigma }_p\times P\times \mathrm{\Phi }(t)$

• La concentration en prométhium à l'équilibre $(\frac{\mathrm{d}P}{\mathrm{d}t}=0)$ est solution de l'équation :

$0={\gamma }_p\times {\mathrm{\Sigma }}_f\times {\mathrm{\Phi }}_o-{\lambda }_p\times P-{\sigma }_p\times P\times {\mathrm{\Phi }}_o$

$P_o=\frac{{\gamma }_p\times {\mathrm{\Sigma }}_f}{{\lambda }_p+{\sigma }_p\times {\mathrm{\Phi }}_o}\times {\mathrm{\Phi }}_o$

• On simplifie l'équation différentielle en remarquant que le terme σ_p × Φ est négligeable devant λ_p, ce qui revient à dire que le prométhium n'est pas un poison, en effet :

${\sigma }_p\times {\mathrm{\Phi }}_o=1350\times 10^{-24}\times 3,2\times 10^{13}=4,346\times 10^{-8}\text{ captures par seconde}$ à comparer à :

${\lambda }_p=3,627\times 10^{-6}\text{désintégrations par seconde}$

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

$\frac{\mathrm{d}P}{\mathrm{d}t}={\gamma }_p\times {\mathrm{\Sigma }}_f\times \mathrm{\Phi }-{\lambda }_p\times P$

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

$\frac{|}{{\gamma }_p\times {\mathrm{\Sigma }}_f}{\lambda }_p\times {\mathrm{\Phi }}_o=\frac{0,0109\times 0,09840\times 3,2\times 10^{13}}{3,627\times 10^{-6}}=9,461\times 10^{15}\mathrm{a}\mathrm{t}/\mathrm{c}{\mathrm{m}}^{\mathrm{3}}$

• L'équation générale est résolue en multipliant l'équation différentielle par ${\mathrm{e}}^{{\lambda }_p\times t}$

$\frac{\mathrm{d}P}{\mathrm{d}t}\times {\mathrm{e}}^{{\lambda }_p\times t}+{\lambda }_p\times P\times {\mathrm{e}}^{{\lambda }_p\times t}={\gamma }_p\times {\mathrm{\Sigma }}_f\times \mathrm{\Phi }(t)\times {\mathrm{e}}^{{\lambda }_p\times t}$

$\frac{\mathrm{d}(P\times {\mathrm{e}}^{{\lambda }_p\times t})}{\mathrm{d}t}={\gamma }_p\times {\mathrm{\Sigma }}_f\times \mathrm{\Phi }(t)\times {\mathrm{e}}^{{\lambda }_p\times t}$

$P(t)=P(o)+({\gamma }_p\times {\mathrm{\Sigma }}_f\times \underset{o}{\overset{t}{\int }}\mathrm{\Phi }(u)\times {\mathrm{e}}^{{\lambda }_p\times u}\mathrm{d}u)\times {\mathrm{e}}^{-{\lambda }_p\times t}$

• On fait l'hypothèse d'un démarrage du réacteur vierge (sans prométhium ni samarium), d'où P(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 du prométhium. On peut alors simplifier la fonction donnant P(t) en remplaçant Φ(t) par la valeur à terme du flux Φ_o

$P(t)={\gamma }_p\times {\mathrm{\Sigma }}_f\times {\mathrm{\Phi }}_o\times \left[\underset{o}{\overset{t}{\int }}{\mathrm{e}}^{{\lambda }_p\times u}\mathrm{d}u\right]\times {\mathrm{e}}^{-{\lambda }_p\times t}$

$P(t)={\gamma }_p\times {\mathrm{\Sigma }}_f\times {\mathrm{\Phi }}_o\times \frac{{\mathrm{e}}^{{\lambda }_p\times t}-1}{{\lambda }_p}\times {\mathrm{e}}^{-{\lambda }_p\times t}$

$P(t)=\frac{{\gamma }_p\times {\mathrm{\Sigma }}_f}{{\lambda }_p}\times {\mathrm{\Phi }}_o\times (1-{\mathrm{e}}^{-{\lambda }_p\times t})=P_o\times (1-{\mathrm{e}}^{-{\lambda }_p\times t})$ Modèle:Boîte déroulante/fin

Évolution de la quantité de samarium 149

• 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 prométhium ni samarium, le samarium atteint sa valeur d'équilibre au bout d'un temps supérieur à celui du prométhium ; dit de façon simpliste, il faut attendre que le prométhium ait atteint un régime stable pour que le samarium soit à l'équilibre.
• La concentration à l'équilibre s'écrit :

$S_o=\frac{{\gamma }_p\times {\mathrm{\Sigma }}_f}{{\sigma }_s}=2,6542\times 10^{16}\mathrm{a}\mathrm{t}/\mathrm{c}{\mathrm{m}}^{\mathrm{3}}$

• La concentration en samarium à l'équilibre ne dépend pas du flux. Ce résultat est à rapprocher de celui concernant le prométhium dont la concentration est proportionnelle au flux. Plus le flux est élevé plus la quantité de prométhium est élevée en valeur relative devant celle de samarium conduisant à un effet d'empoisonnement accru lors de l'arrêt du réacteur
• Le rapport des concentrations en prométhium et samarium s'écrit :

$\frac{P_o}{S_o}=\frac{\frac{{\gamma }_p\times {\mathrm{\Sigma }}_f\times {\mathrm{\Phi }}_o}{{\lambda }_p}}{\frac{{\gamma }_p\times {\mathrm{\Sigma }}_f}{{\sigma }_s}}=\frac{{\sigma }_s}{{\lambda }_p}\times {\mathrm{\Phi }}_o$

σ_s, λ_p sont des constantes physiques le rapport P_o / S_o dépend du dessin du réacteur considéré uniquement par la valeur du flux lequel est inversement proportionnel à la quantité de matière fissile présente dans le cœur. Plus l'usure du combustible est importante plus la quantité de prométhium à l'équilibre est importante par rapport au samarium présent dans le cœur conduisant à un effet d'empoisonnement post arrêt du réacteur augmenté.

• Il est également intéressant de comparer la concentration en samarium avec celle du xénon :

La concentration en xénon à l'équilibre s'écrit :

$X_o=\frac{({\gamma }_i+{\gamma }_x)\times {\mathrm{\Sigma }}_f\times {\mathrm{\Phi }}_o}{{\lambda }_x+{\sigma }_x\times {\mathrm{\Phi }}_o}$

Si le flux augmente la concentration en xénon à l'équilibre plafonne vers la valeur :

$X_{omax}=\frac{({\gamma }_i+{\gamma }_x)\times {\mathrm{\Sigma }}_f}{{\sigma }_x}$

qui ne dépend pas du flux.

Le rapport S_o / X_omax vaut :

$\frac{So}{X_{omax}}=\frac{\frac{{\gamma }_p\times {\mathrm{\Sigma }}_f}{{\sigma }_s}}{\frac{({\gamma }_i+{\gamma }_x)\times {\mathrm{\Sigma }}_f}{{\sigma }_x}}=\frac{{\gamma }_p}{{\gamma }_i+{\gamma }_x}\times \frac{{\sigma }_x}{{\sigma }_s}=\frac{0,0109}{0,004+0,064}\times \frac{2,65\times 10^6}{40140}=10,582$

Dans un réacteur à l'équilibre il y a près de 10 fois plus de samarium 149 que de xénon 135, en dépit de cela l'anti-réactivité du xénon est 2,85 fois plus importante.

• La fonction générale donnant S(t) s'écrit :

$S(t)=S_o\times [(1+\frac{{\lambda }_p}{{\sigma }_s\times {\mathrm{\Phi }}_o-{\lambda }_p})\times {\mathrm{e}}^{-{\sigma }_s\times {\mathrm{\Phi }}_o\times t}-\frac{{\sigma }_s\times {\mathrm{\Phi }}_o}{{\sigma }_s\times {\mathrm{\Phi }}_o-{\lambda }_p}\times {\mathrm{e}}^{-{\lambda }_p\times t}]$

dans laquelle S_o est la valeur à l'équilibre.

Si le flux est variable, on remplace Φ_o par Φ(t)

Si Φ(t) reste constant, le deuxième terme de l'équation converge vers la valeur 1.

• En pratique dans les simulateurs d'entrainement du personnel d'exploitation des réacteurs on retient une modélisation simplifiée dans laquelle seules les concentrations à l'équilibre du prométhium et du samarium en fonction du flux (donc de la puissance thermique) sont modélisées.

• Variation de la concentration en samarium 149 = dS/dt = (Production par décroissance du prométhium) - (Disparition par capture d'un neutron)

$\frac{\mathrm{d}S}{\mathrm{d}t}={\lambda }_p\times P(t)-{\sigma }_s\times S(t)\times \mathrm{\Phi }(t)$

• La concentration en samarium à l'équilibre (dS/dt = 0) est solution de l'équation :

$0={\lambda }_p\times P_o-{\sigma }_s\times S_o\times {\mathrm{\Phi }}_o$

$S_o=\frac{{\lambda }_p\times P_o}{{\sigma }_s\times {\mathrm{\Phi }}_o}$

$S_o=\frac{{\gamma }_p\times {\mathrm{\Sigma }}_f}{{\sigma }_s}$

en remplaçant P_o par sa valeur

$S_o=\frac{0,0109\times 0,09840\mathrm{c}{\mathrm{m}}^{\mathrm{-}\mathrm{1}}}{4,014\times 10^{-20}\mathrm{c}{\mathrm{m}}^{\mathrm{-}\mathrm{2}}}=2,654\times 10^{16}\mathrm{a}\mathrm{t}/\mathrm{c}{\mathrm{m}}^{\mathrm{3}}$

La concentration en samarium 149 à l'équilibre ne dépend pas du flux.

• L'équation différentielle générale est résolue en multipliant l'équation par

${\mathrm{e}}^{\underset{o}{\overset{t}{\int }}\psi (u)\mathrm{d}u}$

avec

$\psi (u)={\sigma }_s\times \mathrm{\Phi }(u)$

Il vient :

$\frac{\mathrm{d}S}{\mathrm{d}t}\times {\mathrm{e}}^{\underset{o}{\overset{t}{\int }}\psi (u)\mathrm{d}u}+\psi (t)\times S(t)\times {\mathrm{e}}^{\underset{o}{\overset{t}{\int }}\psi (u)\mathrm{d}u}=({\lambda }_p\times P(t))\times {\mathrm{e}}^{\underset{o}{\overset{t}{\int }}\psi (u)\mathrm{d}u}$

$\frac{\mathrm{d}(S(t)\times {\mathrm{e}}^{\underset{o}{\overset{t}{\int }}\psi (u)\mathrm{d}u})}{\mathrm{d}t}={\lambda }_p\times P(t)\times {\mathrm{e}}^{\underset{o}{\overset{t}{\int }}\psi (u)\mathrm{d}u}$

Finalement :

$S(t)={\mathrm{e}}^{-\underset{o}{\overset{t}{\int }}\psi (u)\mathrm{d}u}\times [\underset{o}{\overset{t}{\int }}{\mathrm{e}}^{\underset{o}{\overset{u}{\int }}\psi (v)\mathrm{d}v}\times {\lambda }_p\times P(u)\mathrm{d}u+S(o)]$

• Dans l'hypothèse d'un démarrage du réacteur vierge (sans prométhium ni samarium), d'où P(o) = 0, et d'une croissance du flux pour atteindre la valeur Φ_o en quelques dizaines de minutes durée faible devant l'évolution du prométhium et du samarium. On peut simplifier la fonction donnant S(t) en remplaçant Φ (t) par la valeur à terme Φ_o;

Par ailleurs on a vu que :

$P(t)=\frac{{\gamma }_p\times {\mathrm{\Sigma }}_f}{{\lambda }_p}\times {\mathrm{\Phi }}_o\times (1-{\mathrm{e}}^{-{\lambda }_p\times t})$

$S(t)=S_o\times [(1+\frac{{\lambda }_p}{{\sigma }_s\times {\mathrm{\Phi }}_o-{\lambda }_p})\times {\mathrm{e}}^{-{\sigma }_s\times {\mathrm{\Phi }}_o\times t}-\left(\frac{{\sigma }_s\times {\mathrm{\Phi }}_o}{{\sigma }_s\times {\mathrm{\Phi }}_o-{\lambda }_p}\right)\times {\mathrm{e}}^{-{\lambda }_p\times t}]$

avec S_o = valeur à l'équilibre.

• Dans une modélisation simplifiée on remarque que, dès lors que le temps après démarrage est supérieur à une dizaine de jours, le terme ${\sigma }_s\times {\mathrm{\Phi }}_o\times {\mathrm{e}}^{-{\lambda }_p\times t}$ devient faible devant la valeur de ${\lambda }_p\times {\mathrm{e}}^{-{\sigma }_s\times {\mathrm{\Phi }}_o\times t}$ en effet :

${\sigma }_s\times {\mathrm{\Phi }}_o\times {\mathrm{e}}^{-{\lambda }_p\times 864000}=4,014\times 10^{-20}\times 3,2\times 10^{13}\times {\mathrm{e}}^{-3,627\times 10^{-6}\times 864000}=5,592\times 10^{-8}{\mathrm{s}}^{\mathrm{-}\mathrm{1}}$

${\lambda }_p\times {\mathrm{e}}^{-{\sigma }_s\times {\mathrm{\Phi }}_o\times 864000}=3,627\times 10^{-6}\times {\mathrm{e}}^{-4,014\times 10^{-20}\times 3,2\times 10^{13}\times 864000}=1,196\times 10^{-6}{\mathrm{s}}^{\mathrm{-}\mathrm{1}}$

Il est alors possible de ramener l'évolution du samarium à long terme à l'expression plus simple suivante :

$S(t)=S_o\times (1+\frac{{\lambda }_p}{{\sigma }_s\times {\mathrm{\Phi }}_o-{\lambda }_p})\times {\mathrm{e}}^{-{\sigma }_s\times {\mathrm{\Phi }}_o\times t}$

avec S_o = valeur à l'équilibre. Modèle:Boîte déroulante/fin

Empoisonnement du réacteur par le samarium - Effet en réactivité

• Pour un flux initial de Modèle:Nb typique d'un réacteur à eau pressurisée, l'empoisonnement dû au samarium après un fonctionnement prolongé en puissance stable vaut environ 1 300 pcm. L'anti-réactivité supplémentaire apporté par la désintégration du prométhium 149 après un arrêt est de l'ordre de 460 pcm ; cette perte est proportionnelle au flux, et peut avoir des valeurs plus élevées dans le cas d'un réacteur à haut flux.
• Le réacteur doit être conçu pour disposer d'une marge en réactivité suffisante (retrait des barres de contrôle ou dilution du poison soluble) pour pouvoir être redémarré sans problème après un arrêt qui peut toujours être nécessaire de façon inopinée. Le délai procuré par la décroissance du prométhium 149 (plus de 72 heures) peut permettre un redémarrage intermédiaire, toutefois, la décroissance du xénon 135 qui intervient entretemps au bout de 24 heures environ procure en pratique un gain de réactivité plus important que la perte due au samarium 149.
• Parmi les 4 facteurs de la formule $k_{\mathrm{\infty }}=\epsilon \times \eta \times p\times f$, c'est le facteur f qui est concerné par la présence du samarium 149 formé

Cœur vierge avant empoisonnement : $f=\frac{{\mathrm{\Sigma }}_u}{{\mathrm{\Sigma }}_u+{\mathrm{\Sigma }}_m}$

Après empoisonnement $f^{\prime }=\frac{{\mathrm{\Sigma }}_u}{{\mathrm{\Sigma }}_u+{\mathrm{\Sigma }}_m+{\mathrm{\Sigma }}_s}$

$\frac{f^{\prime }}{f}=0,99260=\frac{k_{ef{f}^{\prime }}}{k_{eff}}=\frac{{\mathrm{e}}^{{\rho }^{\prime }}}{{\mathrm{e}}^{\rho }}={\mathrm{e}}^{{\rho }^{\prime }-\rho }$

${\rho }^{\prime }-\rho =100000\times \ln (0,99260)=-742\mathrm{p}\mathrm{c}\mathrm{m}$

soit donc un effet environ 4,9 fois moindre que dans le cas du xénon

• Après arrêt du réacteur la concentration en samarium se trouve augmentée de la valeur de la concentration en prométhium à l'équilibre ; l'empoisonnement en samarium augmente de 264 pcm. On peut voir que le surcroît d'empoisonnement en samarium est inférieur à l'effet de « pic » xénon.

• L'empoisonnement par le samarium affecte le facteur f de la formule des 4 facteurs

f avant empoisonnement = $f=\frac{{\mathrm{\Sigma }}_u}{{\mathrm{\Sigma }}_u+{\mathrm{\Sigma }}_m}$

f' après empoisonnement = $f^{\prime }=\frac{{\mathrm{\Sigma }}_u}{{\mathrm{\Sigma }}_u+{\mathrm{\Sigma }}_m+{\mathrm{\Sigma }}_s}$

Le rapport f / f' exprime la variation de k_eff

$\frac{f}{f^{\prime }}=1+\frac{{\mathrm{\Sigma }}_s}{{\mathrm{\Sigma }}_u+{\mathrm{\Sigma }}_m}$

• Évaluation de ${\mathrm{\Sigma }}_u={\sigma }_5\times U_5+{\sigma }_8\times U_8+{\sigma }_o\times O_u+{\sigma }_z\times Zr$

${\mathrm{\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}\mathrm{c}{\mathrm{m}}^{\mathrm{-}\mathrm{1}}$

Concentration des atomes d'uranium 235 = $U_5=\frac{1760,93}{0,235044}\times \frac{N_A}{26,571\times 10^6}=1,6980\times 10^{20}\mathrm{a}\mathrm{t}/\mathrm{c}{\mathrm{m}}^{\mathrm{3}}$

Concentration des atomes d'uranium 238 = $U_8=\frac{72367-1760,93}{0,238051}\times \frac{N_A}{26,571\times 10^6}=6,722\times 10^{21}\mathrm{a}\mathrm{t}/\mathrm{c}{\mathrm{m}}^{\mathrm{3}}$

Concentration des atomes d'oxygène du combustible = $O_u=2\times (U_5+U_8)=1,3784\times 10^{22}\mathrm{a}\mathrm{t}/\mathrm{c}{\mathrm{m}}^{\mathrm{3}}$

Concentration des atomes de zirconium = $Zr=\frac{19096}{0,091224}\times \frac{N_A}{26,571\times 10^6}=4,744\times 10^{21}\mathrm{a}\mathrm{t}/\mathrm{c}{\mathrm{m}}^{\mathrm{3}}$

D'où :${\mathrm{\Sigma }}_u=0,13460\mathrm{c}{\mathrm{m}}^{\mathrm{-}\mathrm{1}}$

• Évaluation de ${\mathrm{\Sigma }}_m={\sigma }_h\times H+{\sigma }_o\times O_m$

${\mathrm{\Sigma }}_m={\sigma }_h\times H+{\sigma }_o\times O_m$

${\mathrm{\Sigma }}_m=(0,332\times H+0,267\times 10^{-3}\times O_m)\times 10^{-24}\mathrm{c}{\mathrm{m}}^{\mathrm{-}\mathrm{1}}$

Nombre de moles d'eau dans le cœur = $\frac{11073,8\times 1000}{15,9994+2\times 1,00794}=614690$

Concentration des atomes d'hydrogène = $\frac{614690\times 2\times N_A}{26,571\times 10^6}=2,786\times 10^{22}\mathrm{a}\mathrm{t}/\mathrm{c}{\mathrm{m}}^{\mathrm{3}}$

Concentration des atomes d'oxygène du modérateur =$\frac{2,786\times 10^{22}}{2}=1,3932\times 10^{22}\mathrm{a}\mathrm{t}/\mathrm{c}{\mathrm{m}}^{\mathrm{3}}$

${\mathrm{\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,0092543\mathrm{c}{\mathrm{m}}^{\mathrm{-}\mathrm{1}}$

• Évaluation de Σ_s à l'équilibre = ${\mathrm{\Sigma }}_s={\sigma }_s\times S_o=4,014\times 10^{-20}\times 2,654\times 10^{16}=0,0010724\mathrm{c}{\mathrm{m}}^{\mathrm{-}\mathrm{1}}$
• $\frac{f}{f^{\prime }}=1+\frac{{\mathrm{\Sigma }}_s}{{\mathrm{\Sigma }}_u+{\mathrm{\Sigma }}_m}$

$\frac{f}{f^{\prime }}=1+\frac{0,00107294}{0,13460+0,0092543}=1,00745$

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

$\frac{k_{e^{\prime }ff}}{k_{eff}}=\frac{f^{\prime }}{f}=\frac{1}{1,01300}=0,99260$

L'effet en réactivité correspondant vaut :$\frac{k_{ef{f}^{\prime }}}{k_{eff}}=\frac{e^{{\rho }^{\prime }}}{e^{\rho }}=e^{{\rho }^{\prime }-\rho }$

exprimé en pcm

${\rho }^{\prime }-\rho =100000\times \ln (0,98717)=-742\mathrm{p}\mathrm{c}\mathrm{m}$

• Après arrêt du réacteur la concentration en samarium se trouve augmentée de la valeur de la concentration en prométhium soit donc S_arrêt = S_o + P_o = Modèle:Nb + Modèle:Nb = Modèle:Nb soit un facteur 1,356

L'empoisonnement en samarium passe de - 742 pcm à -1 006 pcm soit un empoisonnement accru de - 264 pcm qu'il faudra compenser pour pouvoir redémarrer. Modèle:Boîte déroulante/fin