CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    InterFoam相方程代码疑惑

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

      各位朋友大家好!

      本人最近在看InterFoam里面alphaEqn.H代码(4.1版本),有一段代码不是很明白,希望大家给一些指点

      if (alphaApplyPrevCorr && MULESCorr)
          {
              talphaPhiCorr0 = alphaPhi - talphaPhiCorr0;
          }
      
          if
          (
              word(mesh.ddtScheme("ddt(rho,U)"))
           == fv::EulerDdtScheme<vector>::typeName
          )
          {
              rhoPhi = alphaPhi*(rho1 - rho2) + phiCN*rho2;
          }
          else
          {
              if (ocCoeff > 0)
              {
                  // Calculate the end-of-time-step alpha flux
                  alphaPhi = (alphaPhi - (1.0 - cnCoeff)*alphaPhi.oldTime())/cnCoeff;
              }
      
              // Calculate the end-of-time-step mass flux
              rhoPhi = alphaPhi*(rho1 - rho2) + phi*rho2;
          }
      

      上面这段代码 是利用MULES算法求解出修正的总通量alphaPhi 之后的 一段代码,

      现在不太明白这段代码的作用,,同时希望大家推荐 有关这些代码的 一些参考公式推导 或者公式推导的思路

      恳请好心人解答!同时也随时欢迎交流一下有关InterFoam里面求解相方程的问题

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

        去年我完成了MULES-FCT的相关研究,提交给Computers Fluids,审了半年只找到一个审稿人,被拒稿,后来就没有再投,也不想投了,我计划6月份之前把这部分内容直接挂在网上,MULES-FCT在网上现存一些讨论,你可以自己找找,包括一些sci,但我总觉得写的不详细不通透。这部分内容在论坛里面基本讨论不清楚,不是3 5个方程就能解释清楚的。如果可以等我后续会把详细推导发布。

        这部分内容非常复杂,目前我建议你直接将其理解为对$\alpha$的显性离散后,求解即可。

        CFD高性能服务器 http://dyfluid.com/servers.html
        2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

        1 条回复 最后回复 回复 引用
        • 小
          小考拉 最后由 编辑

          好的,谢谢东岳老师指点!

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