在Fluent或者其他软件中用电势法求解磁流体（MHD）控制方程组时，各个控制方程的计算过程，欢迎赐教与讨论。

[zhaoqijin]

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

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

[bestucan]

fluent 的可以看《ANSYS Fluent Advanced Add-On Modules》，Part V 讲的就是 MHD，其中 chapter 3: Implementation 中开头提到，MHD 计算依托于 UDS 和 UDM。而 UDS 和 UDM 的计算流程《ANSYS Fluent UDF Manual》中有，在其中搜 “solution procedure for” 就能搜到不同求解器的求解流程图。UDS 和 UDM 都在其中。

OpenFOAM 的可以看这个的 3.3.6 节。也可以搜类似的文章。