5圆柱坐标系二维稳态控制方程有限差分Python Solver问题？

我正在开发一个固定床反应器中的热量和质量传递模型，这个反应器中存在吸热反应，稳态。 如下所示的是一个二维图，从左侧有恒定温度和浓度的入口，顶部有恒定温度的墙壁，底部是一个中心线（对称线），温度和浓度没有变化，出口（右侧）也一样，这些都定义了初始条件和边界条件。[图片]三个方程，即质量平衡，能量平衡和反应速率，两个控制方程是通过反应速率耦合在一起的。[图片]我计划使用Python中的有限差分...