基于GEKKO的乙醇生产补料分批生物反应器DAE模拟收敛问题

我正尝试使用GEKKO求解用于乙醇生产的补料分批生物反应器DAE系统，以便后续更便捷地进行优化以最大化乙醇产量。该问题此前已在MATLAB中求解并得到对应结果图，但在GEKKO中使用相同常量与变量值时，无法复现相同结果，仅在短时间范围（如m.time= np.linspace(0,1,11)）下能收敛。现附上原系统方程与GEKKO代码，询问代码存在何种问题？

内容的提问来源于stack exchange，提问作者Chemical Engineering Student