有限体积法温度迭代循环异常问题求助

问题概述 #

用有限体积法做瞬态温度场计算，通过时间-位置嵌套循环，基于j时刻的节点温度推导j+1时刻温度。目前出现两个异常：

• 第7号节点的conduction_out计算值异常
• 第2号节点的conduction_in疑似有误

需求是仅更新内部节点温度，边界节点保持边界条件，总节点数为9，因此设置位置循环范围为1到i-1，但怀疑Python索引逻辑出错，无法定位具体问题。

核心排查点（针对Python索引特性） #

1. 节点编号与0索引的匹配错误 #

Python的数组/列表默认从0开始索引，若你的节点编号是1_{9（共9个），对应的数组索引应为0}8：

• 边界节点：索引0（对应原1号）、索引8（对应原9号）
• 内部节点：索引1_{7（对应原2}8号）

你当前设置的循环范围1到i-1（i=9）即1~8，会把原9号边界节点（索引8）也纳入计算，既违反了“保留边界条件”的要求，还可能导致该节点的conduction_out访问越界索引（9），直接引发异常值。

2. 导热项邻接节点的索引逻辑错误 #

对于内部节点k（数组索引），conduction_in依赖左邻节点k-1，conduction_out依赖右邻节点k+1：

• 若循环误包含索引8（原9号节点），计算conduction_out时会访问索引9，超出数组范围，要么报错要么取到垃圾值
• 若原2号节点对应索引1，计算conduction_in时本应取索引0（左边界节点），如果代码中误写为k+1，就会错误取到索引2的温度值，导致conduction_in异常

3. 修正后的循环逻辑示例 #

调整内部节点的循环范围为1到i-2（i=9，即1~7），确保只处理真正的内部节点，同时校验导热项的索引：

# 假设T是存储当前时刻温度的数组，长度为9（索引0-8）
total_nodes = 9
for time_step in range(total_time_steps):
    # 复制当前温度，避免迭代中覆盖原值
    T_next = T.copy()
    # 仅遍历内部节点：索引1到7（对应原2-8号节点）
    for k in range(1, total_nodes - 1):
        # 左邻导热项：来自k-1（左节点）到k（当前节点）
        conduction_in = (T[k-1] - T[k]) * thermal_coeff
        # 右邻导热项：来自k（当前节点）到k+1（右节点）
        conduction_out = (T[k] - T[k+1]) * thermal_coeff
        # 更新下一时刻温度，边界节点保持不变
        T_next[k] = T[k] + dt * (conduction_in + conduction_out + other_source_terms)
    # 更新温度数组
    T = T_next

进一步排查需要的信息 #

为了精准定位问题，请补充：

• 代码中循环范围定义和导热项计算的具体片段
• 当前结果与预期结果的关键差异（比如第7号节点的异常值是偏大/偏小、是否出现NaN等）

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