Énergie nucléaireÉnergie

Réacteur nucléaire

La neutronique d’un réacteur à partir des premiers principes — cinétique ponctuelle à six groupes avec rétroaction Doppler et xénon-135 : le saut prompt, la falaise de criticité prompte, les excursions de puissance auto-limitées, le puits d’iode post-arrêt, le suivi de charge, et le SCRAM.

DjiniousLabModélisation et simulation

Un graphique en échelle logarithmique de la puissance du réacteur après une insertion de réactivité — s’emballant sans rétroaction, mais s’auto-limitant à un niveau borné une fois la rétroaction en température Doppler incluse

cinétique des neutrons retardés: 6 groupescinétique des neutrons retardés
limité par effet Doppler (contre un emballement ×28 830): 2.4×limité par effet Doppler (contre un emballement ×28 830)
puits xénon, 8 h après l’arrêt: −4 459 pcmpuits xénon, 8 h après l’arrêt
période pour une efficacité de barre de +0,2 $: 37 spériode pour une efficacité de barre de +0,2 $
puissance résiduelle (1 s → 1 h): 6 % → 1 %puissance résiduelle (1 s → 1 h)
notebooks de conception: 10notebooks de conception

Pourquoi un réacteur est pilotable — et pourquoi il s’arrête de lui-même.

Un réacteur nucléaire vit sur le fil du rasoir : la population de neutrons peut doubler en une fraction de seconde, et pourtant les opérateurs le pilotent à la main. Deux choses rendent cela possible — les neutrons retardés, qui ramènent la réponse à des échelles de temps humaines, et la rétroaction en température, qui fait que le cœur combat sa propre montée en puissance. Ce programme modélise les deux, depuis la cinétique ponctuelle à six groupes jusqu’à l’auto-limitation par effet Doppler, en passant par le poison xénon qui peut empêcher le redémarrage d’un réacteur, le pilotage en suivi de charge, et l’arrêt d’urgence.

Les neutrons retardés achètent le temps nécessaire.

Si un réacteur ne fonctionnait qu’avec des neutrons prompts, une petite insertion de réactivité ferait exploser la puissance en quelques millisecondes — incontrôlable. Mais une fraction de pour cent des neutrons apparaît des secondes à des minutes plus tard, issue de la désintégration des produits de fission, et ce sont ces neutrons retardés qui donnent le rythme : insérer un peu de réactivité fait monter la puissance sur une période de quelques dizaines de secondes, assez lentement pour être pilotée. Le modèle capture le saut prompt (un échelon de +0,5 dollar fait bondir la puissance à 2,07×, conformément à la théorie), la relation inhour entre réactivité et période, et la falaise au dollar où les neutrons retardés cessent de compter et le réacteur devient prompt-critique.

Une famille de réponses de puissance du réacteur à des échelons de réactivité en échelle logarithmique, la pente se raidissant à mesure que la réactivité insérée augmente — La puissance du réacteur après des échelons de réactivité, en échelle logarithmique : chaque échelon se stabilise sur une montée exponentielle dont la période se raccourcit à mesure que la réactivité insérée augmente. C’est la réponse contrôlée par les neutrons retardés qui rend un réacteur pilotable — jusqu’à ce que l’insertion approche le dollar.

La courbe inhour — période du réacteur en fonction de la réactivité insérée en dollars, la période s’effondrant à mesure que la réactivité approche le dollar — La courbe inhour : comment la période du réacteur dépend de la réactivité insérée. Un dixième de dollar donne une période confortable de 98 secondes ; un demi-dollar, 5,7 secondes ; approcher un dollar plein et la période s’effondre vers zéro — la falaise de criticité prompte que tout système de pilotage de réacteur est conçu pour éviter.

Le cœur qui combat sa propre montée en puissance.

La réactivité n’est pas fixe — elle dépend de la température. À mesure que le combustible chauffe, l’élargissement Doppler des résonances d’absorption capture davantage de neutrons, ce qui fait chuter la réactivité. Cette rétroaction négative est ce qui fait qu’une excursion de puissance s’auto-limite au lieu de s’emballer. Le modèle le montre de façon frappante : une insertion de +0,6 dollar qui explose à près de trente mille fois la puissance sans rétroaction atteint à la place un maximum de 2,4× et se stabilise, une fois l’effet Doppler intégré à la boucle. C’est la même physique qui rend un réacteur bien conçu intrinsèquement stable.

La puissance du réacteur après une insertion de réactivité en échelle logarithmique : un emballement exponentiel sans rétroaction, comparé à une réponse bornée et auto-limitée avec rétroaction Doppler — Auto-limitation par effet Doppler : la même insertion de +0,6 dollar s’emballe jusqu’à environ 30 000 fois la puissance sans rétroaction en température (rouge), mais avec la rétroaction Doppler (vert), la montée en température du combustible fait rechuter la réactivité et l’excursion s’auto-limite à quelques fois la puissance. La stabilité intrinsèque, illustrée sur un seul graphique.

La réactivité xénon après l’arrêt d’un réacteur, plongeant plus négativement jusqu’à un puits vers 8 heures avant de remonter sur le jour suivant — Le puits d’iode : après un arrêt, le xénon-135 continue de s’accumuler à partir de la désintégration de l’iode alors même que sa combustion s’arrête, si bien que sa réactivité négative s’approfondit pendant environ huit heures — jusqu’à −4 459 pcm ici — avant de se dissiper. Un réacteur qui ne peut pas surmonter ce puits ne peut pas redémarrer tant qu’il n’est pas passé ; le transitoire qui façonne, on le sait bien, l’exploitation.

La puissance résiduelle après un SCRAM du réacteur, chutant d’environ 6 % de la pleine puissance à une seconde vers 1 % après une heure — Après le SCRAM : même une fois la réaction en chaîne arrêtée, la désintégration des produits de fission continue de faire produire de la chaleur au cœur — environ 6 % de la pleine puissance une seconde après le déclenchement, encore 1 % une heure plus tard. Évacuer cette puissance résiduelle est le problème de sûreté qui définit le refroidissement d’arrêt du réacteur.

Chaque valeur est une valeur que vous pouvez recalculer.

Le notebook de validation finale re-dérive chaque exigence à partir du même modèle de cinétique ponctuelle que construit le programme.

Résultat

Saut prompt (échelon +0,5 $) : 2,07×
Excursion limitée par effet Doppler : maximum 2,4×
Puits d’iode xénon : −4 459 pcm à 8 h
Erreur de poursuite en suivi de charge : 2,4 %
Exigences vérifiées : 6 / 6

Exigence

Saut prompt (échelon +0,5 $) : théorie 2,0
Excursion limitée par effet Doppler : bornée, pas d’emballement
Puits d’iode xénon : post-arrêt
Erreur de poursuite en suivi de charge : 100→60→100 %
Exigences vérifiées : RÉUSSI

Neutronique 0D, un modèle pédagogique et de régulation.

Le réacteur est modélisé en cinétique ponctuelle 0D — la population de neutrons réduite à un seul nombre, sans forme spatiale du flux, donc sans mode d’oscillation xénon spatial (un développement identifié pour la suite). La rétroaction thermique est localisée, le xénon est à nœud unique, et les sections efficaces sont illustratives. C’est un modèle pédagogique et de conception de régulation, pas un code de transport neutronique de qualité licenciement. C’est exactement la fidélité qui construit l’intuition et soutient la logique de régulation — le saut prompt, la relation inhour, les rétroactions de sûreté, le transitoire xénon, le suivi de charge et le SCRAM — sur des chiffres reproductibles, avant un solveur de diffusion spatiale ou de transport Monte-Carlo.