CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    关于compressibleinterFoam中的非线性源项dgdt的请教

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

      李老师对compressibleinterfoam中链接文本具体对体积分数运输方程的推导及后续处理。
      其中的一部分推导我用图片截出f9ac8b0f4278529f323234476336b59.png
      这是该求解器中关于dgdt的代码。
      6a4af49cd2d4eafee83c44ef62cb63a.png
      后续为了保持相分界性的特殊处理。
      65ae80079f8fe1025e2a4e387de74bb.png
      能找到的关于dgdt创建
      3476a004f8911af510721f1c38ce412.png
      我的问题:
      1、dgdt(alpha1*fvc::div(phi))和f9a17a7c7685e4105a4309ecb633278.png

      能对的上吗,我现在想的是我是不是没有推导出来。
      2、dimensionedScalar(dimless/dimTime,0)中的dimless/dimTime是个什么量纲
      谢谢大家<抱拳>

      T 1 条回复 最后回复 回复 引用
      • T
        Tens 讲师 @LYT 最后由 编辑

        @lyt
        1.

        dgdt(alpha1*fvc::div(phi))
        

        这段是在createFields.H文件里的,只是初始化用,而真正计算dgdt是在pEqn.H里

        dgdt =
                    (
                        alpha1*(p_rghEqnComp2 & p_rgh)    
                      - alpha2*(p_rghEqnComp1 & p_rgh)    
                    );
        

        相当于19c274dc-0102-42b3-ad3e-7df1351f0ffe-image.png
        2. dimless量纲为1,dimTime量纲为时间单位s,所以dimless/dimTime就是1/s

        L 1 条回复 最后回复 回复 引用
        • L
          LYT @Tens 最后由 编辑

          @tens 谢谢老师,找到了。

          杨 1 条回复 最后回复 回复 引用
          • 杨
            杨英狄 @LYT 最后由 编辑

            @LYT 我有个问题要请教你,在更新相分数场的时候就要用到dgdt,此时dgdt的表达式就是setfields中的形式,但是这样的话跟公式是对不上的,这个是什么原因呢,谢谢!

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