Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

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

    OpenFOAM
    2
    10
    8645
    Loading More Posts
    • Oldest to Newest
    • Newest to Oldest
    • Most Votes
    Reply
    • Reply as topic
    Log in to reply
    This topic has been deleted. Only users with topic management privileges can see it.
    • 金石为开
      金石为开 last edited by

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

      1 Reply Last reply Reply Quote
      • 李东岳
        李东岳 管理员 last edited by

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

        线上CFD课程 7月1日报名截止 http://dyfluid.com/class.html
        CFD高性能服务器 http://dyfluid.com/servers.html

        1 Reply Last reply Reply Quote
        • 李东岳
          李东岳 管理员 last edited by

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

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

          线上CFD课程 7月1日报名截止 http://dyfluid.com/class.html
          CFD高性能服务器 http://dyfluid.com/servers.html

          金石为开 1 Reply Last reply Reply Quote
          • 金石为开
            金石为开 @李东岳 last edited by

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

            李东岳 1 Reply Last reply Reply Quote
            • 李东岳
              李东岳 管理员 @金石为开 last edited by

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

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

              线上CFD课程 7月1日报名截止 http://dyfluid.com/class.html
              CFD高性能服务器 http://dyfluid.com/servers.html

              金石为开 1 Reply Last reply Reply Quote
              • 金石为开
                金石为开 @李东岳 last edited by

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

                李东岳 1 Reply Last reply Reply Quote
                • 李东岳
                  李东岳 管理员 @金石为开 last edited by

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

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

                  线上CFD课程 7月1日报名截止 http://dyfluid.com/class.html
                  CFD高性能服务器 http://dyfluid.com/servers.html

                  金石为开 1 Reply Last reply Reply Quote
                  • 金石为开
                    金石为开 @李东岳 last edited by

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

                    李东岳 1 Reply Last reply Reply Quote
                    • 李东岳
                      李东岳 管理员 @金石为开 last edited by

                      @金石为开

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

                      线上CFD课程 7月1日报名截止 http://dyfluid.com/class.html
                      CFD高性能服务器 http://dyfluid.com/servers.html

                      金石为开 2 Replies Last reply Reply Quote
                      • 金石为开
                        金石为开 @李东岳 last edited by

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

                        1 Reply Last reply Reply Quote
                        • First post
                          Last post

                        CFD中文网 | 东岳流体 | 京ICP备15017992号-2
                        论坛登录问题反馈可联系 li.dy@dyfluid.com