CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    求解器编译求助

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

      浓度方程如图:0_1511939161106_$NYO@ZHM)G9T5CH)YASEKJX.png
      对应的代码该怎么写?写了一些编译都不通过,其中之一为:![0_1511939317349_DD[S]})BF$JL}7$5{$H3]}B.png](/assets/uploads/files/1511939317270-dd-s-bf-jl-7-5-h3-b-resized.png)
      求助大神帮忙修改一下:crying:

      厚积薄发,贵在坚持!

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

        surfaceScalarField nu("nu",mesh.magSf()/mesh.magSf()*nu)
        solve
        (
        fvm::ddt(c)
        +fvm::div(nu, c)
        -D/(3*phi)*fvc::laplacian(c),
        mesh.solver(“c”)
        );
        

        厚积薄发,贵在坚持!

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

          你这个方程写的不对,\nabla c是矢量,c是标量

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

          txwwbz 1 条回复 最后回复 回复 引用
          • txwwbz
            txwwbz @李东岳 最后由 编辑

            @李东岳 0_1512028978469_8}10SS`XQM(_[%{2]8Q12BG.png
            原文中的方程如图,应该是作者的问题,当做散度的话,这个方程该怎么编译?

            厚积薄发,贵在坚持!

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

              surfaceScalarField cPhi("cPhi", c & mesh.sf());
              
              fvScalarMatrix c
              (
                  fvm::ddt(c)
                  +  fvc::div(cPhi)
                  == source.
              );
              

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

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