Temperatur in einem nuklearen Brennstab

Voraussetzungen

  • mathematische Grundlagen von Randwertproblemen

  • Ableitungsregeln in mehrdimensionalen ODEs höherer Ordnung

  • Reduktion von ODEs höherer Ordnung

Lerninhalte

  • Anwendung der mathematischen Grundlagen von ODEs höherer Ordnung auf thermodynamische Modellgleichungen

Temperatur in einem nuklearen Brennstab#

Der Brennstab eines Nuklearreaktors ist eine brennstoffgefüllte zylindrische Hülle. Innerhalb des Zylinders und auch in der Hülle entsteht Hitze durch Kernreaktionen. Die äußere Oberfläche des Zylinders wird von fließendem Wasser auf \(T_\infty = 473\;\textrm{K}\) heruntergekühlt.

../_images/brennstab.png

Der Wärmeübergangskoeffizient ist \(h =10^4 \; \textrm{W}/(\textrm{m}^2\mathrm{K})\). Die thermische Leitfähigkeit des Hüllenmaterials ist \(k = 16.75 \; \textrm{W}/(\textrm{mK})\). Der innere Radius des Stabes beträgt \(R=1.5 \cdot 10^{-2}\;\textrm{m}\) und die Dicke der Hülle um den Brennstab \(w = 3 \cdot 10^{-3}\;\textrm{m}\).

Konstante

Symbol

Wert

Einheit

Wärmeübergangskoeffizient

\(h\)

\(10^4\)

\(\frac{\textrm{W}}{\textrm{m}^2\textrm{K}}\)

Thermische Leitfähigkeit

\(k\)

\(16.75\)

\(\frac{\textrm{W}}{\textrm{mK}}\)

Wassertemperatur

\(T_\infty\)

\(473\)

\(\textrm{K}\)

Radius des Stabes

\(R\)

\(1.5 \cdot 10^{-2}\)

\(\textrm{m}\)

Hülle des Stabes

\(w\)

\(3 \cdot 10^{-3}\)

\(\textrm{m}\)

Die Temperaturverteilung in der Hülle kann durch Lösen der folgenden Differentialgleichung beschrieben werden:

\[ \frac{\mathrm{d}}{\mathrm{d}r} \left( r k \frac{\mathrm{d}T}{\mathrm{d}r} \right) = -10^8 \cdot e^{-\frac{r}{R}} \]

für \(R \leq r \leq R + w\). Die Randbedingungen sind

\[\begin{split} \begin{align} &\frac{\mathrm{d}T}{\mathrm{d}r}(r=R) = - \frac{6.32 \cdot 10^5 \mathrm{W}\mathrm{m}^{-2}}{k} \notag \\ &\frac{\mathrm{d}T}{\mathrm{d}r}(r=R+w) = - \frac{h}{k}(T(r=R+w) - T_\infty) \notag \end{align} \end{split}\]

Lösen Sie die Differentialgleichung mit der Matlab-Routine bvp4c.