CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    代码实现的问题

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

      大家好,我最近在修改OpenFoam里的雷诺应力模型。因为在对比OpenFOAM和Fluent之间的湍流模型的区别时,发现Fluent里面的快速压力应变项比Openfoam里多了由对流引起的一项C。
      phi12.png
      其中CodeCogsEqn-6.gif
      CodeCogsEqn-7.gif
      所以想添加CodeCogsEqn-8.gif
      这一项到代码里。但是改了很多种表达都通不过编译,比如说定义
      volTensorField C(fvm:: div(alphaRhoPhi,R));
      然后在方程里使用 dev(C)
      都不行,于是想问问大家有什么代码思路。编程小白,真诚求助

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

        啊,我现在改成这样了:
        外层定义volSymmTensorField devR(dev(R));
        R方程里面用 fvm::div(alphaRhoPhi,devR)
        编译可以通过,而且貌似是可以代表上述表达的,不知道这样对不对?

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

          更新:使用fvm的时候可以编译但计算不行,改成fvc 可以计算接下来就是看能不能收敛了

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