Introduction au CFD en génie nucléaire

Les réacteurs nucléaires fonctionnent en exploitant la chaleur générée par les réactions de fission contrôlées. Les barres de combustible qui contiennent le combustible nucléaire sont les sources de chaleur primaires, et l'efficacité et la sûreté de l'ensemble du réacteur dépendent de l'élimination fiable de cette chaleur tout en maintenant la température du combustible dans des limites sûres.

Dans une tige de combustible nucléaire, la chaleur est effectuée à travers la boule de combustible, à travers l'écart (si présent), à travers le revêtement, et enfin dans le flux de liquide de refroidissement. Le transfert de chaleur conjugué – la combinaison de conduction solide et de convection fluide – et le mélange turbulent dans le canal de liquide de refroidissement nécessitent des techniques de simulation sophistiquées.

Par exemple, les simulations à haute fidélité aident à prédire le flux critique de chaleur, le départ de l'ébullition du nucléate et les températures centrales du combustible, tous les paramètres clés pour les marges de sûreté du réacteur. Les progrès récents dans le calcul à haute performance ont permis de réaliser des modèles à grande échelle de sous-canaux, réduisant ainsi le besoin de boucles expérimentales coûteuses.

Équations de régulation pour le transfert de chaleur et le débit de fluide

Pour l'analyse des tiges de combustible nucléaire, les équations qui régissent les équations doivent tenir compte des effets de compressibilité (souvent négligeables pour les flux liquides à basse vitesse), des propriétés dépendantes de la température et, éventuellement, des flux entraînés par la flottabilité.

  • Équation de continuité:[ Assure la conservation de la masse dans le domaine du liquide de refroidissement.
  • équations de Navier-Stokes: Décrire le mouvement du fluide, y compris les contraintes visqueuses et les tourbillons turbulents.
  • Équation énergétique:[ Modéliser le transfert de chaleur dans les régions fluides et solides, y compris la production de chaleur volumétrique dans le combustible.

Dans le cas de l'ANSYS Fluent, ces équations sont discrétées à l'aide d'une méthode de volume fini. Les ingénieurs doivent choisir des modèles de turbulence et des schémas de solution appropriés pour équilibrer la précision et le coût de calcul. Pour les applications nucléaires, le modèle réalisable k-ε ou Shear Stress Transport (SST) k---- est couramment utilisé pour saisir avec précision le transfert de chaleur à proximité du mur.

Un aspect important est la modélisation de la production de chaleur à l'intérieur de la boule de combustible. La distribution de puissance de fission est souvent axiale et radiale non uniforme. INSYS Fluent permet à l'utilisateur de définir une source de chaleur volumétrique en fonction de sa position ou par l'intermédiaire d'une fonction définie par l'utilisateur (UDF).

Configuration de la simulation dans Ansys Fluent

Géométrie et génération de mesh

La première étape consiste à créer un domaine de calcul qui comprend la tige de combustible (solide) et le canal de refroidissement environnant (fluide). Dans la plupart des analyses pratiques, la symétrie est exploitée pour réduire la taille du modèle. Par exemple, une seule tige de combustible dans un réseau carré ou hexagonal est modélisée en utilisant des conditions de limite périodiques pour représenter l'ensemble complet. La géométrie consiste généralement en:

  • La granule de combustible (cylindrique, éventuellement avec un trou central pour certains modèles).
  • L'écart entre le gaz rempli et le revêtement (parfois modélisé comme une interface mince avec une résistance thermique spécifiée).
  • Le revêtement métallique (par exemple, Zircaloy-4 ou M5®).
  • Le domaine liquide du liquide de refroidissement (eau, sodium liquide ou hélium selon le type de réacteur).

Un filet hexaédral structuré est préféré pour les géométries cylindriques car il s'aligne sur la direction du flux et réduit la diffusion numérique. INSYS Meshing ou ICEM CFD peut générer de tels mesh avec des améliorations de la couche limite près de la paroi de revêtement pour capturer les gradients de température raides dans la sous-couche visqueuse. Une bonne pratique est de s'assurer que la distance de paroi sans dimension (y+) est d'environ 1 pour les modèles de turbulence qui résolvent la sous-couche visqueuse. Pour le transfert de chaleur conjugué, le mesh doit être continu à l'interface fluide-solide, ou un couplage conformal est nécessaire pour assurer la conservation de l'énergie.

Propriétés du matériel et conditions de délimitation

Dans les barres de combustible nucléaire, la conductivité thermique, la densité et la chaleur spécifique du combustible (dioxyde d'uranium ou MOX) varient fortement en fonction de la température et de la combustion. ANSYS Fluent permet de définir les propriétés dépendantes de la température par des ajustements polynomiaux ou des fonctions linéaires à la pièce. Pour le revêtement, les propriétés isotropes ou anisotropes peuvent être utilisées. La densité et la viscosité du liquide de refroidissement varient également en fonction de la température et pour les scénarios d'ébullition, des modèles multiphasés deviennent nécessaires.

Les conditions limites typiques comprennent :

  • Inlet: Vitesse ou débit massique, température, intensité de turbulence et diamètre hydraulique.
  • Outlet: Sortie de pression avec contre-pression ambiante ou système.
  • Tailles: Pas de glissement pour le liquide, souvent avec des conditions adiabatiques ou symétriques pour les limites éloignées. Les surfaces de la tige de combustible peuvent être définies comme des parois couplées pour le transfert de chaleur conjugué.
  • Limites périodiques: Pour simuler un motif de treillis répété.
  • Source de chaleur: Production de chaleur volumétrique dans la région des granulés de combustible (W/m3).

Pour les simulations transitoires, les conditions initiales doivent être spécifiées (p. ex., un champ de température uniforme ou une solution à l'état d'équilibre convergent). Les fonctions définies par l'utilisateur (FDU) peuvent être écrites en C pour mettre en œuvre des conditions limites personnalisées, comme les coefficients de transfert de chaleur qui dépendent des régimes d'écoulement ou des profils de puissance axiale.

Modèles de physique et paramètres du solvant

Pour un liquide de refroidissement monophasé (la plupart des systèmes de refroidissement PWR et BWR en fonctionnement normal), les choix suivants sont typiques:

  • Équation énergétique: Active, avec un chauffage visqueux souvent négligeable.
  • Modèle de turbulence:[ SST k---- pour une bonne performance près du mur ou standard k-ε avec traitement mural amélioré.
  • Radiation:[ Généralement, les barres de combustible ne sont pas considérées comme étant à l'intérieur du réservoir en raison de matériaux opaques, mais peuvent être pertinentes dans les réacteurs refroidis au gaz à très haute température.
  • Construction: Permete si la circulation naturelle est importante; utiliser l'approximation de Boussinesq pour des changements de densité modérée.

Pour le solveur, l'algorithme couplé à la pression est recommandé pour les flux à l'état d'équilibre car il accélère la convergence. Un système de vent ascendant de deuxième ordre pour l'impulsion et l'énergie fournit une précision suffisante. Les facteurs de sous-relaxation peuvent nécessiter un ajustement pour une convergence stable, en particulier avec une production de chaleur forte et des propriétés dépendantes de la température.

Transfert de chaleur conjugué dans les barres de combustible

Dans les barres de combustible nucléaire, le CHT est essentiel parce que le revêtement est un solide qui sépare le combustible du liquide de refroidissement. ANSYS Fluent gère le CHT en résolvant simultanément l'équation d'énergie dans les zones solides et fluides, les conditions d'interface assurant automatiquement la continuité de la température et du flux thermique.

Deux approches de THC sont possibles chez Fluent :

  • Mesures non-conformelles:[ Les mailles solides et fluides sont indépendantes et appariées à l'interface en utilisant une condition de limite d'interface de maille. Cela permet différentes densités de mailles dans chaque région.
  • Mesures de forme:[ Le même maillage s'étend à la fois solide et fluide, avec des zones cellulaires attribuées à différents matériaux.

Pour les barres de combustible, la présence de l'écart de gaz introduit une résistance thermique supplémentaire. Souvent, l'écart est modélisé comme une paroi mince avec une résistance thermique définie par l'utilisateur (conductivité) qui explique la conductivité de l'écart (fonction de gonflement du carburant, composition du gaz et pression de contact).

La température maximale de la ligne médiane du carburant est un paramètre de sécurité critique qui doit rester en dessous du point de fusion (p. ex. ~2850°C pour UO2). La température de la ligne médiane du carburant permet d'identifier les changements dans le débit du liquide de refroidissement, le niveau de puissance ou l'épaisseur du revêtement.

Considérations relatives à la modélisation de la turbulence

La prédiction précise du transfert de chaleur du revêtement au liquide de refroidissement repose fortement sur la modélisation du flux turbulent dans le sous-canal. La turbulence à proximité du mur contrôle le coefficient convectif de transfert de chaleur. Le choix du modèle de turbulence peut conduire à des différences significatives dans la température de paroi prédite et le flux de chaleur.

Pour les réacteurs refroidis à l'eau, le nombre de Reynolds dans les sous-canaux se situe généralement entre 104 et 105, ce qui indique un débit turbulent. Le modèle standard k-ε avec traitement mural amélioré est un choix courant car il est robuste et efficace sur le plan calculateur. Cependant, il peut sous-estimer le transfert de chaleur dans les débits avec une forte courbure ou séparation. Le modèle SST k---- combine la formulation k------ dans la région proche du mur avec k-ε dans le champ lointain, donnant souvent de meilleurs résultats pour le transfert de chaleur dans les débits internes.

Il est important de valider la sélection des modèles de turbulence à partir de données expérimentales (p. ex., expériences de mélange de sous-canaux PWR ou corrélations de transfert de chaleur de faisceau de tiges).

Analyser les résultats et optimiser la conception

Distribution de la température et identification des points chauds

Une fois la simulation CFD convergente, le post-traitement commence. Les outils de visualisation intégrés ANSYS Fluents affichent des courbes de température sur les surfaces de la tige de combustible et le volume du liquide de refroidissement. Les profils de température axial le long de la ligne centrale de la tige et de la surface extérieure de revêtement sont extraits pour évaluer les marges de sécurité.

Par exemple, une simulation pourrait montrer que la température maximale de la surface de revêtement se produit près du sommet de la tige en raison du chauffage du liquide de refroidissement le long de la longueur. En ajustant le profil de puissance axiale ou en augmentant le débit du liquide de refroidissement, les ingénieurs peuvent réduire la température de pointe.

Les modèles de débit de refroidissement et l'efficacité de l'élimination de la chaleur

Dans les géométries de faisceau, le mélange de flux croisé entre les sous-canaux est important pour l'élimination uniforme de la chaleur. CFD peut quantifier les taux de mélange et identifier les régions de faible échange de liquide de refroidissement, comme près des grilles d'espacement ou dans les réseaux serrés. Les ingénieurs peuvent alors modifier la géométrie du réseau ou utiliser des caractéristiques de renforcement du flux.

Les coefficients de transfert de chaleur (CCT) sont calculés à partir de la simulation en divisant le flux thermique mural par la différence entre la température de paroi et la température locale du fluide en vrac. Ces CCT sont comparés à des corrélations (p. ex. Dittus–Boelter, Gnielinski) pour valider le modèle.

Les simulations de CFD transitoires, comme un voyage de pompe ou un accident de perte de débit, sont particulièrement utiles pour évaluer la réponse temporelle de la tige de carburant à la température.

Avantages de l'utilisation du CFD dans la conception des réacteurs nucléaires

L'intégration du CFD dans la conception des tiges de combustible nucléaire offre des avantages tangibles au-delà des codes de systèmes unidimensionnels traditionnels (p. ex. RELAP5, TRACE).

  • Haute résolution spatiale: CFD fournit des champs de température et de débit en trois dimensions, captant des effets tels que des blocages de débit ou des variations de température azimutales qui sont perdues dans les modèles de paramètres de masse.
  • Coûts expérimentaux réduits :[ De nombreuses itérations de conception peuvent être testées pratiquement avant de construire des maquettes physiques ou d'effectuer des expériences coûteuses hors du pieu.
  • Une analyse de sécurité améliorée :[ CFD aide à évaluer des scénarios difficiles ou dangereux à reproduire expérimentalement, comme des blocages partiels du débit ou des accidents de perte de liquide de refroidissement.
  • Optimisation de la conception:[ Les études paramétriques (puissance variable, débit, géométrie) peuvent être automatisées en utilisant AnsYS Workbench pour trouver des configurations optimales qui maximisent l'efficacité thermique tout en maintenant des marges de sécurité.
  • Soutien à l'octroi de licences:[ Les autorités de réglementation acceptent de plus en plus les résultats validés des CFD dans le cadre du cas de sûreté, surtout lorsque les données expérimentales sont rares.

Défis pratiques et solutions

Malgré sa puissance, CFD pour les barres de combustible nucléaire présente plusieurs défis. Le coût informatique reste élevé pour les modèles à cœur complet ou même à assemblage complet. Pour y remédier, les ingénieurs utilisent souvent des modèles représentatifs à échelle sous-canale avec des conditions de limites périodiques, réduisant le nombre de cellules de centaines de millions à quelques millions. ]La qualité de la maille[ près de l'interface entre l'écart de calage et le refroidisseur de revêtement nécessite une attention particulière – les cellules à haut rapport d'aspect peuvent causer des instabilités numériques.

L'incertitude de propriété matérielle est un autre problème. La conductivité thermique du combustible se dégrade avec la combustion, et la dégradation n'est pas toujours bien caractérisée. Les études de sensibilité utilisant des gammes de propriétés peuvent lier l'incertitude.

La validation[ est essentielle pour la crédibilité.Les résultats de CFD doivent être comparés aux données expérimentales, comme ceux des expériences de transfert de chaleur de l'OCDE/NEA . Sans validation, la simulation reste une prédiction non vérifiée. La meilleure pratique consiste à comparer les profils de température, la chute de pression et les coefficients de transfert de chaleur à de multiples conditions de fonctionnement.

Tendances futures du DFC pour le génie nucléaire

Le rôle de la FCD dans l'analyse des barres de combustible nucléaire se développe rapidement. Des méthodes de haute fidélité comme La grande simulation de Eddy (LES) deviennent réalisables pour les domaines sous-canaux à mesure que la puissance de calcul augmente.Les ERP capturent des structures turbulentes transitoires avec plus de précision, améliorant les prédictions des pics de mélange et de transfert de chaleur. La FCD multiphase pour les débits bouillants (p. ex., modèles eulériens à deux fluides avec bouillie de paroi) est en cours de développement pour simuler directement les conditions DNB et post-DNB, réduisant ainsi la dépendance à l'égard des corrélations empiriques.

L'apprentissage automatique entre également dans le champ : les réseaux neuronaux formés sur les bases de données CFD peuvent servir de substituts rapides pour le contrôle en temps réel ou la quantification de l'incertitude. INSYS Fluent soutient désormais l'intégration avec les outils d'optimisation pilotés par l'IA.

Les plateformes open-source comme OpenFOAM gagnent en traction dans le monde universitaire, mais AnsYS Fluent reste le standard de l'industrie en raison de sa validation complète, interface conviviale et support robuste pour les UDF. L'énergie nucléaire continue de jouer un rôle dans la production d'énergie à faible teneur en carbone, la demande d'analyse thermique précise et à haute résolution ne fera que croître.

Conclusion

La capacité de modéliser le transfert de chaleur, la turbulence et les géométries complexes fournit aux ingénieurs les idées nécessaires pour concevoir des réacteurs plus sûrs et plus efficaces. En suivant les pratiques exemplaires établies pour la création de géométrie, le maillage, la configuration physique et la validation, les analystes peuvent produire des prédictions fiables qui appuient à la fois l'optimisation de la conception et l'évaluation de la sûreté.

Pour les ingénieurs qui sont nouveaux dans ce domaine, il est crucial d'investir du temps dans la compréhension de la physique sous-jacente et la validation des données expérimentales.