关于李东岳老师icoFoam求解器解析细节问题



  • 关于李东岳老师icoFoam求解器解析细节问题,我在网站上看到他的解析,看到细节有点不明白,想跟网友确认一下是笔误还是自己理解的错误呢?0_1463915940740_upload-d533e1b6-418e-4d0e-a7cd-a14befd1e7dd
    还有后面的:加入压力梯度,我看一开始公式是有压力除以密度rho的,后面离散加的时候没有rho,这是咋回事?是漏写了吗还是就是这样子的?谢谢



  • 非常感谢!
    已经将公式(13)和公式(14)更新。



  • 另外,在你发帖的时候,我应该把公式的速度重新换了一个符号(粗体),看截图还是之前的?是另存的pdf?

    在前天更新的icoFoam解析我把速度矢量$\mathbf{U}$和其他文档的速度分量$u$进行了区分。并且强调了PISO循环。



  • @李东岳 您的意思是变成了粗体的U来表示矢量吗?我的图就是截图您之前的内容。另外非常感谢你做的工作,让我们这些小研究生学习openfoam过程中少走弯路,真的很感谢。另外您以后会讲讲动网格,重叠网格相关东西可以吗?比如interDyMFoam的东西,我会用到动模型,动网格等。谢谢~~



  • @金石为开
    可以参考这里和这里
    interDyMFoam就是interFoam和DyMFoam的结合。

    重叠网格OpenFOAM目前没有,我也没做过植入,不过extend版本的浸没边界法很有意思。



  • @李东岳 我也第一次看这个piso循环。之前问题差不多明白了,想问下关于循环的。for (int corr=0; corr<nCorr; corr++),for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++),if (nonOrth == nNonOrthCorr)等关于些循环步分别是什么意思呢?PISO循环的截止条件是什么?nCorr,nNonOrthCorr指什么呢?这个循环和非正交网格有关系?问题有点多哈,不过确实不理解。。。



  • @金石为开
    for (int corr=0; corr<nCorr; corr++)为PISO的内迭代,for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)为非正交修正。PISO循环无截止条件。不过Issa在1985年的论文里面已经指出2-3次即可。在我们使用中,需要强制定义2次或者3次。

    有关非正交修正,需要码公式。在我处理完LES文档后说一下。也欢迎交流。



  • @李东岳 恩恩,非常感谢,也就是这里nCorr的值也就是2-3的意思吧?



  • @金石为开

    是的,为2或者3或者其他。



  • @李东岳 东岳老师我想说一下自己的理解,不知道对不对:最后速度更新的时候,U = HbyA - rAUfvc::grad(p)中HbyA是通过HbyA = rAUUEqn.H();//公式(15)得出的吧,也就是Ur,预测的速度而不是n+1步的速度对不?我还没看SIMPLE,不过按照我的理解是不是HbyA要随着速度变化不断更新?那内迭代的最终截止条件是什么呢?非常感谢!


Log in to reply