请教一个问题，关于SIMPLEC算法的改进！

[航中魁]

我想在Of中实现上述算法，但是不知道红色的部分该如何解决？这个算法是SIMPLEC算法更一般的形式，我想测试一下它的性能，各位大神有没有什么好的想法？非常感谢

[李东岳]

你贴的图就是SIMPLEC，OpenFOAM已经植入进去了

[航中魁]

@李东岳 你好，我上面贴的是SIMPLEC算法的一个改进，of中植入的SIMPLEC算法de为de=Ae/(ae-anb)，上面这个算法中，de中多了两项，分别是Unb/Ue，这两项我不知道应该在原始的SIMPLEC算法上面怎么改？求指导，非常感谢

[李东岳]

SIMPLEC算法：
\begin{equation}
A_e u_e' - A_N u_e' =A_N u_N' - A_N u_e' -\Delta p
\end{equation}
确认一下，你的方程16是：
\begin{equation}
A_e u_e' - A_N u_N/u_e u_e' = -\Delta p
\end{equation}
?

[航中魁]

@李东岳 你好，对的，你的（2）式就是方程（16），把Ue提出来就是方程（17），所以李老师，有没有什么办法，解决这个问题呢，在每个迭代层次上，如何处理Un/Ue
，这是我想要表达的，非常感谢！

[李东岳]

新版本中的SIMPIEC处理的更成熟，在SIMPLE下定义关键词就可以，新的求解器代码我还没看。

看老的OpenFOAM求解器，SIMPLE算法通过H()植入，SIMPLEC算法通过H1()植入，就是这个：

HbyA = rAU*UEqn.H();//公式(15), SIMPLE
HbyA = rAU*UEqn.H1();//SIMPLEC

或许你可以看看把他植入到fvMatrix下面的H1()函数下，比如写一个H2()函数 :cheeky:

[航中魁]

@李东岳 多谢李老师！