请教一个问题,关于SIMPLEC算法的改进!



  • 0_1506508984259_TIM图片20170927183330.png
    我想在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:



  • @李东岳 多谢李老师!


登录后回复
 

与 CFD中文网 的连接断开,我们正在尝试重连,请耐心等待