用电势法求解瞬态磁流体（MHD）控制方程组时，离散后的方程求解顺序问题。

[zhaoqijin]

刚在Fluent区发了一下，然后感觉Openform区搞CFD算法的比较多，所以特再发一次讨教。用电势法求解瞬态磁流体（MHD）控制方程组时（图1），离散后的控制方程计算的流程给出两种可能如下（图2和图3）。问题如下：
①图2和图3中关于时间步推进的表示是否是对的？
②电荷守恒方程应该在什么地方进行求解？
③电磁力和焦耳热的计算是在判断收敛前，还是收敛后呢？如果是Fluent中自带的MHD模块，图2和图3的求解流程哪个是是正确的，或者是都有错误？
④如果用udf编程计算电磁力和焦耳热，是否下图的两种求解流程是否都可以实现，且都是对的？

图1 控制方程

4b374628-1dab-4275-a48f-

图2 流程1

图3 流程2

每个问题都可以讨论，然后不是磁流体研究领域的也欢迎讨论。希望大家不吝赐教

[李东岳]

你确定Fluent那面求解的算法用的是速度场与电流场的强耦合么

[zhaoqijin]

@李东岳 强耦合应该不算吧，他是先假定速度场，根据磁场和速度场计算电势，然后再计算电流、电磁力和焦耳热，将计算得到的电磁力和焦耳热作为动量源项和能量源项带入下一次的迭代计算。

[zhaoqijin]

@李东从这个过程而言我觉得流程图算是合理的吧

[李东岳]

图2看起来更合理

[zhaoqijin]

@李东岳 :确是这样，每一个迭代步内理应是都要把相关方程进行求解的。

[zhouxu]

@zhaoqijin 我也遇到了这样的问题，我是用MHD方程求解等离子体。没看出来流程1和2有什么本质的区别？我认为求解顺序不重要吧