温度场，热流值

模拟40年都是正确的（图1），但模拟到第50年时，出现下部温度波动呢（图2），这是为什么呢？
边界条件：上部温度边界；下部热流值。左右边界为wall。