其中,$\alpha$为热扩散系数,$q_{i,j}$表示该小格子的热源(如果有)。 由于本问题为二维问题,我们需要同时考虑 $x$ 方向和 $y$ 方向的温度变化。假设小格子的长度为 $h$,时间步长为 $\Delta t$,则有: $\frac{T_{i,j}^{n+1} - T_{i,j}^{n}}{\Delta t} = \alpha \frac{T_{i+1,j}...