我在环形系统中有一个带有 Robin 边界条件的 PDE 方程，我应该通过有限体积法求解它：

$$\begin{align} \frac{\partial T_f}{\partial t} - k \left(\frac{\partial^2 T_f}{\partial r^2} + \frac{1}{r} \frac{\partial T_f}{\partial r} \right) = -Q \end{align}$$

通过有限体积法对内点进行离散化没有任何问题，但是在$r_{ID}$和$r_{OD}$，它有 robin 边界条件，我不知道如何离散化它。

$$\begin{align} \left.\frac{\partial T_f}{\partial r}\right|_{r=r_{ID}} &= h(T_f-T_w ) \\ \left.\frac{\partial T_f}{\partial r}\right|_{r=r_{OD}} &= h(T_f-T_g ) \end{align}$$

（h&k 是常数）
如果你能帮助我离散化这个边界条件，我将不胜感激。