simpleFoam中的 rAtU

lllwonderliquid

前辈们好！本人刚刚接触openfoam，在阅读simpleFoam求解器程序时对于其中rAtU的意义始终不能理解......
附上simpleFoam的pEqn.H文件中的一段代码

 if (simple.consistent())
    {
        rAtU = 1.0/(1.0/rAU - UEqn.H1());
        phiHbyA +=
            fvc::interpolate(rAtU() - rAU)*fvc::snGrad(p)*mesh.magSf();
        HbyA -= (rAU - rAtU())*fvc::grad(p);
    }

其中rAtU这个系数在算法中对应的含义让我困惑，如果我没理解错的话，rAU应当是1/Ap,那rAtU的含义是？希望能得到大神们的解答！

李东岳

那个是SIMPLEC算法哦，在1/Ap之后减去了临点的影响。

好像在2015年我发了个simplec的文章，不过后来删掉了，最近事情比较多，可以再9月份之后写个simpleFoam解析 :cheeky:

lllwonderliquid

@李东岳感谢东岳前辈的指点，恍然大悟:cheeky:

程迪

@李东岳一直看不懂所谓的"临点"是什么，对应的全称或者英文术语是啥呀？

李东岳

@程迪
邻点：相邻的节点。参考《数值传热学》里面用了大量的这个称呼。比如网格点P，西边的邻点W，东边的邻点E

程迪

@李东岳
........我就说这个称呼好诡异，FVM一般用单元和面，不用结点的。一直搞不清楚邻点是哪个点。。。

李东岳

@lllwonderliquid
链接文本
看方程25，是SIMPLE和PISO中略去临点思想的本源，SIMPLEC需要重新写个网页，以后有时间写写:cheeky:
simpleFoam和pisoFoam，icoFoam，pimpleFoam变动其实不超过10%，理解理解icoFoam有助于理解其他各种Foam。

李东岳

@lllwonderliquid 已更新simpleFoam解析 http://dyfluid.com/simplefoam.html