Dans une industrie où la sécurité est primordiale, les modèles de calcul offrent un moyen rentable et sûr d'explorer une vaste gamme de conditions d'accident sans exposer le public, l'environnement ou des infrastructures coûteuses à des dommages. Au cours des dernières décennies, les progrès de la puissance de calcul et des méthodes numériques ont élargi le rôle de la modélisation, qui consiste à passer de simples approximations unidimensionnelles à des simulations multiphysiques très détaillées qui couplent l'hydraulique thermique, la neutronique et la mécanique structurelle en temps quasi réel. Cet article explore l'état actuel de la modélisation des scénarios d'accident de la PWR, les types de modèles utilisés, leurs applications et les orientations futures qui promettent d'améliorer davantage la sécurité des réacteurs.

Comprendre les RPP et la nécessité de modéliser

Dans un réacteur à eau pressurisée, l'eau est maintenue à haute pression (environ 15,5 MPa ou 2250 psi) à l'intérieur de la boucle primaire de refroidissement, ce qui l'empêche de bouillir même à des températures supérieures à 300 °C. Cette eau pressurisée absorbe la chaleur des assemblages de combustible du cœur du réacteur et la transporte vers un générateur de vapeur, où il transfère de la chaleur à une boucle secondaire, produisant de la vapeur pour conduire une turbine. La boucle primaire contient le réacteur, les pompes de refroidissement, le pressuriseur et le faisceau de tubes de générateur de vapeur. La boucle secondaire est séparée physiquement, assurant que les matériaux radioactifs de la boucle primaire n'atteignent pas la turbine ou l'environnement.

La complexité des systèmes de RAP introduit de nombreux scénarios d'accidents potentiels, à la fois des accidents de base de conception (ACD) et des accidents hors de la conception (ACD), qui doivent être compris pour assurer un fonctionnement sûr. Les accidents comme un accident de perte de coloc (LOCA), un accident de réactivité ou un blackout de station impliquent une cascade de processus physiques interdépendants : débit de fluide par des géométries complexes, transfert de chaleur du combustible au liquide de refroidissement, cinétique des neutrons dans le cœur et comportement du matériau sous une température et une pression extrêmes.

En outre, la modélisation joue un rôle central dans l'élaboration d'évaluations probabilistes des risques (APR) et de lignes directrices sur la gestion des accidents graves. Sans modèles de calcul précis, l'industrie nucléaire ne disposerait pas de la base quantitative nécessaire pour démontrer qu'une installation peut résister en toute sécurité à toute une gamme d'événements d'initiation postulés.

Types de modèles informatiques utilisés

L'analyse de sûreté nucléaire repose sur un ensemble diversifié de modèles informatiques, chacun conçu pour capter des phénomènes physiques spécifiques à des échelles appropriées.Les modèles sont généralement classés en codes de systèmes thermohydrauliques, en dynamique des fluides informatiques (CFD), en codes cinétiques des neutrons et en logiciels d'analyse structurelle.

Codes du système thermique-hydraulique

Les codes du système thermique-hydraulique sont les chevaux de travail de l'analyse des accidents par les câbles de refroidissement, simulant le système global de refroidissement du réacteur (SRC) et sa réponse aux transitoires et aux accidents en utilisant une représentation unidimensionnelle ou quasi-dimensionnelle des voies de circulation (p. ex., tuyaux, pompes, récipients). Les codes tels que RELAP5 (développé par la Commission de réglementation nucléaire des États-Unis), TRACE[ et ATHLET résolvent les équations de régulation de la masse, de l'impulsion et de la conservation de l'énergie pour un fluide biphasé (eau et vapeur).

Dynamique des fluides informatiques (CFD)

Les codes système offrent une vue au niveau du système, mais les modèles CFD captent des modèles de débits tridimensionnels transitoires et des effets locaux critiques pour certains accidents. Des outils CDF tels que ANSYS Fluent[, STAR-CCM+[ et OpenFOAM[ résolvent les équations Navier-Stokes avec des modèles de turbulence (p. ex., k-ε, k-----SST) et des modèles multiphasés (p. ex., eulérien-eurien, volume de fluide).

  • Stratification thermique dans le système de refroidissement du réacteur pendant la circulation naturelle
  • Mélange d'eau de refroidissement de l'eau de refroidissement du cœur de secours avec le liquide de refroidissement primaire
  • Bouillant de piscine et flux thermique critique (CHF) sur les surfaces des barres de combustible
  • Distribution de l'hydrogène dans le confinement après un accident grave
  • Performance de la buse de pulvérisation dans le refroidissement de confinement

Le coût élevé du calcul du CFD, qui nécessite souvent des milliers d'heures de CPU pour un seul passage, limite son utilisation aux régions localisées ou à de courtes échéances. Cependant, sa résolution s'est révélée inestimable pour valider et améliorer les modèles de code système, en particulier dans les zones où les approximations unidimensionnelles échouent, comme le mélange du plénum ou le contournement du plénum.

Simulations de transport de neutrons

Les codes cinétiques de Neutron calculent la distribution spatiale et énergétique des neutrons (flux) couplés à la température et à la densité de rétroaction du liquide de refroidissement et du carburant. Deux grandes classes existent:

  • [Codes déterministes (p. ex., PARCS, SIMULATE-3): Résolvez l'équation de diffusion des neutrons ou l'équation de transport en modes stationnaire et transitoire. Ils sont couramment utilisés pour la conception de cœur et le support opérationnel.
  • [p. ex., MCNP[, Serpent[): Utiliser un échantillonnage aléatoire pour simuler les antécédents individuels de neutrons, fournissant des distributions de flux et des vitesses de réaction très précises.

Les codes thermohydrauliques/neutroniques couplés permettent aux analystes de simuler la rétroaction de réactivité lors d'accidents (p. ex., élargissement de Doppler dans le carburant, écoulement du liquide de refroidissement, rétroaction de la température du modérateur).

Intégrité structurelle et modèles chimiques

Les codes d'analyse des éléments finis (AFE) tels que ANSYS Mécanique et ABAQUS[ évaluent les marges de contrainte, de contrainte et de défaillance du réservoir de pression du réacteur, des tuyauteries et du confinement sous des charges thermiques et sous pression. Les modèles de réaction chimique, souvent intégrés dans des codes d'accident graves comme MAAP ou MELCOR[, prédisent des phénomènes tels que :

  • oxydation du Zircaloy et génération d'hydrogène
  • Interaction entre le claquage des combustibles
  • Délivrance et transport de produits de fission (iode, césium, etc.)
  • Interaction cœur-béton fondu lors de scénarios de fusion du noyau

Applications dans l'analyse des scénarios d'accident

Les normes de sécurité de l'AIEA exigent que les analyses de sécurité déterministes soient effectuées en utilisant des codes validés pour tous les accidents de la base de conception. Voici les principales catégories de scénarios et les approches de modélisation utilisées.

Accidents liés à la perte de coloc (LOCA)

Pour le LBLOCA, les codes du système thermique-hydraulique simulent la phase de soufflage, le remplissage et la remise en état du cœur par des systèmes de refroidissement du cœur de secours. Les principaux phénomènes comprennent le débit critique en deux phases à travers la rupture, l'annulation induite par la dépressurisation et la propagation frontale d'extinction à mesure que l'eau monte du fond pour remouvoir les tiges de combustible. Des modèles CFL fiables (limite du débit contrecourant) sont essentiels pour prédire si l'eau peut atteindre le cœur. Le CDF est souvent utilisé pour analyser l'efficacité de l'injection descendante. Pour le SBLOCA, le comportement naturel de circulation et de joint de boucle devient dominant, exigeant des codes système avec nodalisation détaillée des boucles primaires et secondaires.

Accidents réactifs (RIA)

Les événements de l'EIR, comme l'éjection de la tige de commande ou un faisceau de tiges larguées, provoquent une insertion soudaine de la réactivité dans le noyau, entraînant une excursion rapide en puissance. La température du carburant peut s'accentuer, provoquant une interaction mécanique avec les boules de plomb (PCMI) et, dans de graves cas, une fragmentation et une dispersion du carburant. Ces accidents sont habituellement simulés par l'accouplement d'un code cinétique des neutrons avec un modèle thermique-hydraulique qui tient compte de la rétroaction du carburant et du liquide de refroidissement.

Arrêt de la station et accidents graves

Une panne de station (SBO) entraîne la perte de toute la puissance AC, entraînant la défaillance des systèmes de refroidissement actifs. L'usine compte sur des systèmes passifs (ou des générateurs diesel de secours) et, éventuellement, le système de refroidissement du réacteur chauffe et pressurise. Si le refroidissement n'est pas rétabli, le noyau peut découvrir, surchauffer et fondre.Les codes d'accident sévère comme MELCOR[ ou [MAAP sont utilisés pour modéliser la progression complète de la dépressurisation du système vers les dommages de cœur, la défaillance du navire et la réponse au confinement.

Autres scénarios importants

  • Transient anticipé sans éclaboussure (ATWS) : Une défaillance des tiges de contrôle à insérer pendant un événement prévu (p. ex., écoulement de turbine, perte d'eau d'alimentation). Les codes neutroniques/thermiques-hydrauliques couplés simulent l'insertion de réactivité à partir de la dilution du bore ou de l'écoulement du liquide de refroidissement et évaluent si la plante peut survivre sans brouillage par rétroaction inhérente.
  • Cassage de la ligne principale de vapeur (MSLB)[: Une rupture du côté secondaire qui augmente la charge du générateur de vapeur, potentiellement surrefroidissant le système primaire et insérant une réactivité positive.
  • Réaction de confinement et sécurité de l'hydrogène[: Après un accident grave, l'hydrogène provenant de l'oxydation du zirconium peut s'accumuler et poser un risque de déflagration ou de détonation.

Avantages et limites de la modélisation par calcul

Avantages

La modélisation computationnelle offre des avantages indéniables pour l'analyse de la sécurité des câbles de protection individuelle. Premièrement, elle réduit la nécessité de procéder à des essais d'effet intégral coûteux et à grande échelle (ETI) et à des essais d'effet séparés (ETI), bien qu'ils demeurent essentiels pour la validation. Les modèles peuvent explorer beaucoup plus de variations de scénarios — par exemple, la taille des ruptures, l'emplacement, les actions de l'opérateur, la configuration du système — que les expériences le permettent.

Limites et défis

Malgré leur puissance, les modèles calculateurs ont des limites inhérentes. Un défi principal est la quantification d'incertitude[. Les paramètres d'entrée – comme les conditions initiales, les conditions limites et les propriétés matérielles – sont incertains et les structures des modèles simplifient la réalité. Les utilisateurs du code doivent appliquer les méthodes de la meilleure estimation plus l'incertitude (BEPU), comme le préconisent des règlements comme le 10 CFR 50.46 de l'USNRC pour la LOCA. De plus, le couplage multiphysique souffre d'erreurs numériques et de problèmes de synchronisation en échelon de temps.

Orientations futures

L'industrie nucléaire fait avancer activement la modélisation informatique afin de tirer parti des capacités de calcul à haute performance (HPC), de l'apprentissage automatique et des approches plus intégrées.

  • Couplage multiphysique de haute fidélité: Les efforts tels que le CASL (Consortium for Advanced Simulation of Light Water Reactors) visent à relier les performances neutronique, thermique et énergétique à une résolution sans précédent, en utilisant des CFD 3D et des neutroniques Monte Carlo sur des ordinateurs exascales. Ces modèles détaillés permettront de réduire le conservatisme dans les marges de sécurité.
  • Modélisation assistée par l'IA et conduite de données: L'apprentissage automatique est appliqué pour accélérer la propagation de l'incertitude, développer des modèles de substitution pour les simulations d'accidents et améliorer la surveillance de l'état des installations.
  • Digital Twins[: Des jumeaux numériques en temps réel de PWR, combinant les données de capteurs et la simulation en direct, pourraient fournir un avertissement précoce des conditions anormales et soutenir la prise de décision de l'opérateur en cas d'accident.
  • Modélisation améliorée des accidents graves[: Des modèles améliorés pour le comportement du corium, l'interaction combustible-concret fondu et la physique des aérosols sont en cours de développement pour soutenir la prochaine génération d'analyses de sécurité, en particulier pour les usines du groupe III+.

À mesure que la puissance de calcul continue de croître et que les données provenant d'installations expérimentales (p. ex., les tests ATLAS, LOBI et PKL deviennent plus accessibles, la modélisation computationnelle jouera un rôle de plus en plus central dans le fonctionnement des RFP avec les normes de sécurité et de fiabilité les plus élevées.

En conclusion, la modélisation computationnelle a transformé la façon dont l'analyse des scénarios d'accident est effectuée pour les réacteurs à eau pressurisée. En fournissant un banc d'essai virtuel pour explorer un vaste espace de défaillances potentielles, les modèles sont devenus essentiels pour améliorer la conception, le développement de cas de sûreté et la conformité réglementaire.