Python中abfluss函数返回空列表或内存错误问题排查求助

• 循环变量未更新：V_geblieben和neu_tiefe仅在进入while循环前计算了一次，循环内部未对这两个变量做更新。如果初始计算的V_geblieben >= 0，循环条件会永久成立，程序会无限向列表中追加元素最终触发内存溢出；如果初始计算的V_geblieben < 0，循环不会执行，直接返回空列表，完全匹配两类异常现象。
• 流量计算逻辑错误：剩余体积计算中存在无物理依据的/V0操作，流出体积的正确计算逻辑为单位时间流出量 * 时间步长，不需要除以初始体积。
• 时间统计逻辑错误：原代码直接修改时间增量dt作为时间点，不符合业务需求中“存储各时间增量、最后一个元素对应排空时刻”的要求。
• 缺少边界处理：未对体积为0的排空状态做特殊处理。

• 调整变量更新逻辑：在循环内部逐次计算当前步的剩余体积、湖泊深度、累计时间，保证循环条件会随着迭代逐步触发终止。
• 修正物理公式：删除冗余的/V0操作，采用正确的流出体积计算公式。
• 拆分时间变量：将时间步长和累计时间拆分为两个独立变量，固定时间步长，逐次累加得到各时间节点。
• 增加边界处理：当剩余体积小于等于0时，手动添加排空状态的时间和深度后终止循环。

from numpy import pi
p = 0.036
def tiefe(V):
    a = (3*V*p**2)/pi
    z = a**(1/3)
    return z

U_unten = 35.56
A = 2                                                                       
def abfluss(V0, dt_step):
    t_list = []
    Z_list = []
    current_t = 0
    current_V = V0
    # 写入初始状态
    t_list.append(current_t)
    Z_list.append(tiefe(current_V))
    
    while True:
        outflow = U_unten * A * dt_step
        current_V -= outflow
        current_t += dt_step
        # 排空后写入终止状态
        if current_V <= 0:
            t_list.append(current_t)
            Z_list.append(0.0)
            break
        # 写入当前步状态
        current_z = tiefe(current_V)
        t_list.append(current_t)
        Z_list.append(current_z)
     
    return t_list, Z_list

# 测试调用：初始体积50000，时间步长3600秒
print(abfluss(50000,3600))

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