CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    fvc::DDt是什么意思?

    OpenFOAM
    6
    12
    2546
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • Y
      yfclark 讲师 最后由 编辑

      我看到有人算
      DDtnew.png
      用了如下的表达式

      fvc::DDt(phi,p)
      

      我试着写入代码里,可以运行,但是在网上搜索了半天也没有发现这部分代码的定义。有老哥看到过这个运算的定义吗?
      Openfoam 2.3.0

      1 条回复 最后回复 回复 引用
      • zym604
        zym604 最后由 编辑

        看上去像是ddt+div
        查了一下代码:
        https://github.com/OpenFOAM/OpenFOAM-3.0.x/blob/master/src/finiteVolume/finiteVolume/fvc/fvcDDt.C

        template<class Type>
        tmp<GeometricField<Type, fvPatchField, volMesh> >
        DDt
        (
            const surfaceScalarField& phi,
            const GeometricField<Type, fvPatchField, volMesh>& psi
        )
        {
            tmp<GeometricField<Type, fvPatchField, volMesh> > ddtDivPhiPsi
                = fvc::ddt(psi) + fvc::div(phi, psi);
        
            if (phi.mesh().moving())
            {
                return ddtDivPhiPsi - fvc::div(phi + phi.mesh().phi())*psi;
            }
            else
            {
                return ddtDivPhiPsi - fvc::div(phi)*psi;
            }
        }
        

        这是3.0的版本,不确定2.3是不是也是如此。从这个版本里看,就是ddt+div。

        至于定义,应该是随体导数(Material derivative)。
        c665ddfc-c0f7-4e01-b907-e970d6dce5fe-image.png

        strive_DUT 1 条回复 最后回复 回复 引用
        • bestucan
          bestucan 版主 副教授 最后由 编辑

          6.0版的定义。估计之前的难理解,又加了个U。
          zzz.png

          滚来滚去……~(~o ̄▽ ̄)~o 滚来滚去都不能让大家看出来我不是老师么 O_o

          异步沟通方式(《posting style》from wiki)(下载后打开):
          https://www.jianguoyun.com/p/Dc52X2sQsLv2BRiqnKYD
          提问的智慧(github在gitee的镜像):
          https://gitee.com/bestucan/How-To-Ask-Questions-The-Smart-Way/blob/master/README-zh_CN.md

          1 条回复 最后回复 回复 引用
          • 李东岳
            李东岳 管理员 最后由 编辑

            感谢诸位大佬分享,赠送66元漏油奔驰代金券

            :xiexie: :xiexie:

            线上CFD课程开始报名:http://www.dyfluid.com/class.html

            CFD高性能服务器 http://dyfluid.com/servers.html

            1 条回复 最后回复 回复 引用
            • Y
              yfclark 讲师 最后由 编辑

              谢谢各位大佬

              1 条回复 最后回复 回复 引用
              • strive_DUT
                strive_DUT @zym604 最后由 编辑

                @zym604 随体导数最后不是梯度项吗?源代码中两个值作差都是散度呀?这个该怎么理解?

                X 1 条回复 最后回复 回复 引用
                • X
                  xpqiu 教授 @strive_DUT 最后由 编辑

                  @strive_DUT
                  根据张量等式
                  $$
                  \nabla \cdot ( \mathbf{U} \phi) = \mathbf{U} \cdot \nabla \phi + \phi (\nabla \cdot \mathbf{U})
                  $$
                  所以,
                  $$
                  \frac{\partial \phi}{\partial t} + \mathbf{U} \cdot \nabla \phi = \frac{\partial \phi}{\partial t} + \nabla \cdot ( \mathbf{U} \phi) - \phi (\nabla \cdot \mathbf{U})
                  $$
                  这里等式右边的减号就是二楼贴出来代码里面减号的来源了。

                  另外,

                  @bestucan 在 fvc::DDt是什么意思? 中说:

                  6.0版的定义。估计之前的难理解,又加了个U。
                  zzz.png

                  这个是来搞笑的么?

                  strive_DUT bestucan 2 条回复 最后回复 回复 引用
                  • strive_DUT
                    strive_DUT @xpqiu 最后由 编辑

                    @xpqiu 多谢多谢!这个张量等式我还真不了解!长知识了:xiexie:

                    1 条回复 最后回复 回复 引用
                    • bestucan
                      bestucan 版主 副教授 @xpqiu 最后由 编辑

                      @xpqiu 是啊,我的签名就声明了我就是来搞笑的:chitang: 滚来滚去,滚去滚来,薛定谔的滚蛋状态:mihu: 布朗滚蛋运动:mihu:

                      像您这样的能答上来的人没有出现,只有先搞笑一下,不至于别人问的问题没人搭理。还能刺激一下大家的回答欲。忘了在哪看的研究,一般不愿意回答的问题要是有个错误回答,就会引来很多纠错。因为纠错的成就感更强一些。

                      :146:

                      滚来滚去……~(~o ̄▽ ̄)~o 滚来滚去都不能让大家看出来我不是老师么 O_o

                      异步沟通方式(《posting style》from wiki)(下载后打开):
                      https://www.jianguoyun.com/p/Dc52X2sQsLv2BRiqnKYD
                      提问的智慧(github在gitee的镜像):
                      https://gitee.com/bestucan/How-To-Ask-Questions-The-Smart-Way/blob/master/README-zh_CN.md

                      X 1 条回复 最后回复 回复 引用
                      • 李东岳
                        李东岳 管理员 最后由 编辑

                        小弟只能在这里多谢各位捧场了,最近事情太多了,早五晚八的,哎...每天就是在瞎忙..

                        :mihu:

                        线上CFD课程开始报名:http://www.dyfluid.com/class.html

                        CFD高性能服务器 http://dyfluid.com/servers.html

                        bestucan 1 条回复 最后回复 回复 引用
                        • bestucan
                          bestucan 版主 副教授 @李东岳 最后由 编辑

                          @东岳 老师言重了,我从东岳流体上学到很多很多,现在来论坛逛也是回报,不算付出。虽然东岳流体上的内容不收费,那是老师的好,我白占便宜我受不了。:huahua:

                          对,我是来论坛打工换东岳流体上的内容的:high:

                          滚来滚去……~(~o ̄▽ ̄)~o 滚来滚去都不能让大家看出来我不是老师么 O_o

                          异步沟通方式(《posting style》from wiki)(下载后打开):
                          https://www.jianguoyun.com/p/Dc52X2sQsLv2BRiqnKYD
                          提问的智慧(github在gitee的镜像):
                          https://gitee.com/bestucan/How-To-Ask-Questions-The-Smart-Way/blob/master/README-zh_CN.md

                          1 条回复 最后回复 回复 引用
                          • X
                            xpqiu 教授 @bestucan 最后由 编辑

                            @bestucan 在 fvc::DDt是什么意思? 中说:

                            @xpqiu 是啊,我的签名就声明了我就是来搞笑的:chitang: 滚来滚去,滚去滚来,薛定谔的滚蛋状态:mihu: 布朗滚蛋运动:mihu:

                            像您这样的能答上来的人没有出现,只有先搞笑一下,不至于别人问的问题没人搭理。还能刺激一下大家的回答欲。忘了在哪看的研究,一般不愿意回答的问题要是有个错误回答,就会引来很多纠错。因为纠错的成就感更强一些。

                            :146:

                            :146: :146:

                            1 条回复 最后回复 回复 引用
                            • First post
                              Last post