面向对象气候模拟器两相连水体温度平均计算方案问询

问题根源 #

• 数值稳定性约束：你当前采用的显式欧拉差分格式仅在小时间步长下稳定，当time_scale超过导热系数对应的阈值时，会出现数值振荡，导致温度计算结果违背能量守恒，出现高于理论平均值的异常结果。
• 重复计算逻辑漏洞：如果你的所有物体都会遍历自身邻居执行热传递计算，那么A与B这对邻居会被分别计算两次热传递，等效于时间步长被翻倍，进一步放大了数值不稳定的问题。

基础物理模型 #

忽略对流、辐射等其他热传递形式，仅考虑接触导热的两个相同质量、相同比热容物体，温度演变符合如下规律：

1. 最终平衡温度固定为两物体初始温度的平均值，满足能量守恒：T_eq = (T1 + T2) / 2
2. 温度差随时间指数衰减：ΔT(t) = ΔT0 * e^(-kt)，其中ΔT0为初始温差，k为和导热率、接触面积、热容、质量相关的自定义导热系数。

修正方案 #

逻辑层面：避免重复计算 #

不要在遍历邻居时实时修改温度，改为先统一计算所有物体的温度增量，再批量更新温度，或者通过约定计算规则（如仅高ID物体计算与低ID邻居的热传递）避免重复计算。

数值层面：解决大时间步长不稳定问题 #

增加传递比例钳位，保证单次步长内的热传递最多抹平当前温差，不会出现超调：

# 第一步：所有物体预计算温度增量，不修改当前温度
def pre_calculate_delta(self, time_scale=0.001, heat_conductivity=1.0):
    self.temp_delta = 0.0
    for n in self._neighbours:
        temp_diff = n.temperature - self.temperature
        # 限制传递比例最大为1，保证不会出现温度超调
        transfer_ratio = min(time_scale * heat_conductivity, 1.0)
        self.temp_delta += transfer_ratio * temp_diff / 2

# 第二步：所有物体统一应用温度增量
def apply_temp_update(self):
    self.temperature += self.temp_delta

每次仿真帧更新时，先遍历所有仿真对象执行pre_calculate_delta，再遍历所有对象执行apply_temp_update即可。该方案下无论time_scale多大，最终温度都会收敛到理论平均值，不会出现异常升高的问题。

内容的提问来源于stack exchange，提问作者Nathan Marotte