Calcul De L Apparition Du R Gime D Bullition Nucl E

Ce calculateur estime l’apparition de l’ébullition nuclée pour l’eau à partir de la pression locale, de la température du fluide, de la température de paroi et du coefficient de convection. L’approche combine une température de saturation issue de l’équation d’Antoine et une corrélation simplifiée de surchauffe de paroi à l’ONB pour une première vérification d’ingénierie.

Calculateur ONB

Guide expert du calcul de l’apparition du régime d’ébullition nuclée

Le calcul de l’apparition du régime d’ébullition nuclée, souvent abrégé ONB pour onset of nucleate boiling, est une question centrale en thermique avancée, en ingénierie nucléaire, dans les échangeurs compacts et dans l’analyse de sûreté des systèmes à haute densité de flux. Lorsque la paroi chauffée dépasse localement une combinaison critique entre la température de saturation et la surchauffe minimale nécessaire à l’activation de microcavités, les premières bulles stables apparaissent. Cette transition modifie profondément le coefficient d’échange, le champ de température, les marges vis-à-vis du CHF et parfois les vibrations thermo-hydrauliques.

1. Que signifie réellement l’apparition de l’ébullition nuclée ?

Dans un écoulement monophasique chauffé, la chaleur est d’abord évacuée par convection forcée. À mesure que la température de paroi augmente, certaines aspérités ou microcavités de la surface peuvent devenir des sites de nucléation. L’ONB ne signifie pas forcément une ebullition massive dans tout le canal. Il s’agit d’un événement local, généralement observé d’abord sur la paroi, au voisinage du pic de flux thermique, dans une région où la température de fluide est la plus proche de la saturation.

Le phénomène dépend de plusieurs familles de paramètres :

• la pression locale, qui fixe la température de saturation de l’eau,
• la température du fluide au voisinage de la paroi,
• la température réelle de la paroi,
• la rugosité, la mouillabilité et l’état métallurgique de la surface,
• la vitesse et le régime d’écoulement, qui influencent le coefficient de convection,
• le sous-refroidissement local, c’est-à-dire l’écart entre la température de saturation et la température du liquide.

En pratique, l’ONB apparaît lorsque la surchauffe de paroi par rapport à la saturation devient suffisante pour maintenir des bulles stables malgré la condensation rapide dans un liquide encore sous-refroidi.

2. Pourquoi ce calcul est-il si important en génie nucléaire ?

Dans les cœurs de réacteurs, la densité de flux au niveau des crayons combustibles peut être élevée alors même que la marge à la crise d’ébullition doit rester rigoureusement maîtrisée. L’apparition de l’ébullition nuclée n’est pas forcément un événement défavorable. Dans de nombreux cas, elle améliore même temporairement le transfert thermique grâce à l’agitation induite par le départ de bulles. En revanche, une estimation erronée de son apparition peut fausser :

1. l’évaluation des températures de gaine,
2. le bilan de sous-refroidissement le long du canal,
3. la prévision des régimes diphasiques locaux,
4. les marges au DNB ou au CHF selon le type de réacteur,
5. les résultats d’un dimensionnement d’échangeur à haute performance.

Dans un REP, par exemple, l’eau reste généralement sous-refroidie en cœur à haute pression, si bien que l’ONB peut survenir localement avant que le liquide en vrac n’atteigne la saturation. Dans un REB, l’ébullition dans les canaux est au contraire une partie normale du fonctionnement. La valeur précise du seuil conditionne alors le taux de vide, la stabilité et la marge thermique.

3. Le principe physique utilisé par ce calculateur

Le calculateur ci-dessus procède en trois étapes simples, adaptées à un pré-diagnostic :

1. conversion de la pression d’entrée en bar absolus,
2. calcul de la température de saturation de l’eau par inversion de l’équation d’Antoine,
3. estimation d’une surchauffe minimale de paroi à l’ONB via une corrélation simplifiée décroissante avec la pression.

La logique d’ingénierie est la suivante :

• on calcule la surchauffe réelle de paroi, soit ΔT réelle = T paroi – T saturation,
• on estime la surchauffe minimale requise, soit ΔT ONB,
• si ΔT réelle dépasse ΔT ONB + marge, alors l’apparition du régime d’ébullition nuclée est jugée probable.

Le calculateur estime aussi le flux thermique associé par la relation convective classique q” = h × (T paroi – T fluide). Cette quantité est utile pour situer votre point de fonctionnement, mais il faut se rappeler qu’en présence d’ONB, le coefficient d’échange effectif peut évoluer rapidement. En d’autres termes, le modèle présenté est une porte d’entrée opérationnelle, pas un code système à deux fluides ni une corrélation de qualification réglementaire.

4. Pression et température de saturation de l’eau : données de référence

Le premier facteur gouvernant l’ONB est la pression. Plus la pression augmente, plus la température de saturation s’élève. Cela modifie la marge disponible entre le liquide sous-refroidi et la surface chauffante. Le tableau ci-dessous rappelle quelques valeurs de saturation couramment utilisées en thermique de l’eau.

Pression absolue	Température de saturation de l’eau	Contexte d’ingénierie typique	Commentaire thermique
1 bar	99,6 °C	Essais de laboratoire, circuits ouverts	Ébullition très accessible, faible sous-refroidissement requis pour l’ONB.
10 bar	179,9 °C	Boucles d’essai pressurisées	La surchauffe de paroi nécessaire reste modérée mais le fluide peut rester largement sous-refroidi.
70 bar	285,8 °C	Ordre de grandeur REB	Les premières bulles peuvent apparaître sur paroi alors que le liquide en vrac n’est pas saturé.
155 bar	344,8 °C	Ordre de grandeur REP	La haute pression décale la saturation très haut et rend le sous-refroidissement du cœur déterminant.

Ces valeurs sont cohérentes avec les tables de vapeur de référence. Elles montrent bien qu’une discussion sur l’ONB sans pression locale fiable reste incomplète. Une erreur de quelques bars peut suffire à déplacer la température de saturation de plusieurs degrés, ce qui modifie immédiatement la marge à l’apparition des premières bulles.

5. Comparaison de conditions thermo-hydrauliques représentatives

Le tableau suivant donne des ordres de grandeur utiles pour comparer des contextes réels d’utilisation de l’eau pressurisée. Ces chiffres sont représentatifs d’applications industrielles et d’analyses ouvertes publiées par des organismes techniques, mais ils ne doivent pas être utilisés comme données de qualification de matériel.

Système	Pression typique	Température fluide typique	Flux surfacique typique	Lecture vis-à-vis de l’ONB
Canal chauffé de laboratoire à eau	1 à 20 bar	90 à 210 °C	0,05 à 0,8 MW/m²	L’ONB est souvent facilement observable et instrumentable.
Réacteur à eau bouillante	65 à 75 bar	270 à 290 °C	0,3 à 1,2 MW/m²	L’ébullition fait partie du fonctionnement normal dans les canaux.
Réacteur à eau pressurisée	150 à 160 bar	290 à 325 °C	0,5 à 1,5 MW/m²	Le liquide reste globalement pressurisé et sous-refroidi, mais l’ONB locale peut rester pertinente.

Ce tableau montre qu’il ne faut jamais juger l’ONB à partir de la seule température de fluide en vrac. Dans un REP, par exemple, la température moyenne du liquide peut être inférieure à la saturation à 155 bar, tout en laissant subsister des zones de paroi nettement plus chaudes à cause du flux radial et de la résistance thermique de la gaine.

6. Méthode pratique de calcul pas à pas

Pour une vérification rapide, on peut suivre la séquence suivante :

1. Déterminer la pression absolue locale. Évitez les confusions entre pression relative et absolue.
2. Calculer la température de saturation. Une table vapeur ou une équation type Antoine convient pour une estimation.
3. Mesurer ou estimer la température de paroi. Elle peut provenir d’une thermocouple, d’un bilan 1D à travers paroi ou d’un calcul CFD conjugué.
4. Évaluer le sous-refroidissement du liquide. C’est la différence entre saturation et température du fluide local.
5. Estimer la surchauffe requise pour l’ONB. On utilise une corrélation adaptée au domaine de pression, au fluide et à la géométrie.
6. Comparer. Si la surchauffe réelle de paroi dépasse la surchauffe requise, l’ONB est probable.
7. Contrôler le flux thermique. Vérifiez que le flux estimé reste cohérent avec les ordres de grandeur de l’installation.

Cette séquence est particulièrement utile lors de campagnes d’essais, de revues de conception thermique ou de contrôles de cohérence entre instrumentation et modèle. Elle permet aussi de hiérarchiser les incertitudes. Très souvent, les plus fortes dispersions viennent non pas de la formule elle-même, mais de la température de paroi locale et de l’état réel de surface.

7. Sources principales d’erreur

• Pression mal définie : une pression mesurée en amont peut être différente de la pression locale près de la zone chauffée.
• Température de paroi indirecte : une thermocouple noyée dans la paroi ne voit pas exactement la température de surface.
• Rugosité et conditionnement de surface : oxydation, dépôts, usure et nettoyage modifient les sites de nucléation.
• Écoulement non uniforme : les zones de recirculation et les pics de vitesse changent le coefficient d’échange local.
• Simplification excessive de la corrélation : une formule unique ne remplace pas les corrélations validées pour une géométrie donnée.

Pour des études de sûreté, de qualification combustible ou de dimensionnement final, on emploie donc des outils plus complets : corrélations expérimentales dédiées, calculs sous-canal, CFD multiphase, ou codes système thermo-hydrauliques validés sur une base d’essais. Le présent calculateur doit être compris comme un outil d’avant-projet, de formation ou d’analyse rapide.

8. Comment interpréter le résultat du calculateur

Trois cas de figure se présentent généralement :

• Résultat négatif net : la température de paroi reste en dessous de la température requise pour l’ONB. L’écoulement est probablement en convection monophasique.
• Résultat proche du seuil : l’installation est dans une zone sensible. Une variation mineure de flux, de débit, de pression ou d’état de surface peut déclencher les premières bulles.
• Résultat positif franc : l’apparition du régime d’ébullition nuclée est probable. Il faut alors examiner l’impact sur le coefficient d’échange, le glissement des bulles, le taux de vide local et la marge à une transition plus sévère.

Une bonne pratique consiste à appliquer une marge conservatrice, comme le permet ce calculateur. En exploitation réelle, cette marge couvre partiellement les fluctuations instrumentales, la dispersion des surfaces et l’incertitude de corrélation.

9. Références et ressources d’autorité

Pour approfondir les bases physiques, les données d’eau saturée et le contexte nucléaire, consultez aussi :

• NIST Chemistry WebBook pour les propriétés thermophysiques et les constantes utiles.
• MIT OpenCourseWare pour des supports universitaires de thermique et de thermo-hydraulique.
• U.S. Nuclear Regulatory Commission pour le cadre technique et réglementaire des réacteurs à eau.

Ces sources complètent utilement les manuels de transfert de chaleur et les bases d’essais. Si votre application concerne un composant nucléaire classé, un échangeur critique ou un essai à haute puissance, il est recommandé de recourir à des corrélations validées dans votre domaine exact de pression, géométrie, état de surface et qualité de mesures.