simpleFoam中的 rAtU
-
前辈们好!本人刚刚接触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:
-
@李东岳 感谢东岳前辈的指点,恍然大悟: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
论坛登录问题反馈可联系 li.dy@dyfluid.com