CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    计算出现负相率和负湍动能

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

      1b12af59-4375-480a-bb3b-de5d017006aa-image.png
      如图,是一个搅拌釜内的液固两相流,用的twoPhaseEulerFoam,计算过程中逐渐发散,出现了很多负值,请问这个可能是什么原因导致的?对网格进行了质量检查,网格为非结构网格,这个非正交性的值是偏大还是正常范围?
      d5838f8d-b15f-4118-8771-7b8350ed29e5-image.png
      fvSolution设置为:(为了不求解温度,我将h/e的迭代数设为0,是否对计算有影响?)

      solvers
      {
          "alpha.*"
          {
              nAlphaCorr      1;
              nAlphaSubCycles 2;
      
              smoothLimiter   0.1;
      
              implicitPhasePressure yes;
              solver          smoothSolver;
              smoother        symGaussSeidel;
              tolerance       1e-9;
              relTol          0;
              minIter         1;
          }
      
          p_rgh
          {
              solver          GAMG;
              smoother        DIC;
              tolerance       1e-8;
              relTol          0;
          }
      
          p_rghFinal
          {
              $p_rgh;
              relTol          0;
          }
      
          "U.*"
          {
              solver          smoothSolver;
              smoother        symGaussSeidel;
              tolerance       1e-5;
              relTol          0;
              minIter         1;
          }
      
          "(h|e).*"
          {
              solver          smoothSolver;
              smoother        symGaussSeidel;
              // tolerance       1e-6;
              relTol          0;
              minIter         0;
              maxIter         0;
          }
      
          "Theta.*"
          {
              solver          smoothSolver;
              smoother        symGaussSeidel;
              tolerance       1e-6;
              relTol          0;
              minIter         1;
          }
      
          "(k|epsilon).*"
          {
              solver          smoothSolver;
              smoother        symGaussSeidel;
              tolerance       1e-5;
              relTol          0;
              minIter         1;
          }
      }
      
      PIMPLE
      {
          nOuterCorrectors 3;
          nCorrectors      1;
          nNonOrthogonalCorrectors 0;
      }
      
      relaxationFactors
      {
          equations
          {
              ".*"            1;
          }
      }
      

      希望各位大佬帮忙解答,感谢!

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

        另外附上我的离散格式,出现负值是否因为离散格式设置错误?未包含有界限制?

        ddtSchemes
        {
            default     Euler;
        }
        
        gradSchemes
        {
            default     Gauss linear;
        }
        
        divSchemes
        {
            default                         none;
        
            "div\(phi,alpha.*\)"            Gauss vanLeer;
            "div\(phir,alpha.*\)"           Gauss vanLeer;
        
            "div\(alphaRhoPhi.*,U.*\)"      Gauss limitedLinearV 1;
            "div\(phi.*,U.*\)"              Gauss limitedLinearV 1;
        
            "div\(alphaRhoPhi.*,(h|e).*\)"  Gauss limitedLinear 1;
            "div\(alphaRhoPhi.*,K.*\)"      Gauss limitedLinear 1;
            "div\(alphaPhi.*,p\)"           Gauss limitedLinear 1;
        
            div(alphaRhoPhi.particles,Theta.particles) Gauss limitedLinear 1;
        
            "div\(alphaRhoPhi.*,(k|epsilon).*\)"  Gauss limitedLinear 1;
        
            div((((alpha.water*thermo:rho.water)*nuEff.water)*dev2(T(grad(U.water))))) Gauss linear;
        
            div((((thermo:rho.particles*nut.particles)*dev2(T(grad(U.particles))))+(((thermo:rho.particles*lambda.particles)*div(phi.particles))*I)))  Gauss linear;
        }
        
        laplacianSchemes
        {
            default     Gauss linear uncorrected;
            bounded     Gauss linear uncorrected;
        }
        
        interpolationSchemes
        {
            default     linear;
        }
        
        snGradSchemes
        {
            default     uncorrected;
            bounded     uncorrected;
        }
        
        1 条回复 最后回复 回复 引用
        • 李东岳
          李东岳 管理员 最后由 编辑

          openfoam几,用新的openfoam比较好
          这个求解器对网格之类比较敏感,最好先上迎风格式

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

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

            @李东岳 用的of7,那请问我这个网格的非正交性算大还是正常,如果网格质量不好那使用reactingTwoPhaseEulerFoam来计算是否会网格的敏感性有所改善?还是用最新的of9里面的multiphaseEulerFoam会好一点?

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

              @caijinjin multiphaseEulerFoam > reactingTwoPhaseEulerFoam > twoPhaseEulerFoam

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

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

                @李东岳 好的,谢谢东岳老师

                C 1 条回复 最后回复 回复 引用
                • C
                  cyberk @Caijinjin 最后由 编辑

                  @caijinjin 您好,我也遇到了类似的问题 想问一下您解决了吗?

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

                    负相率和负湍动能

                    这个很难处理。一般只能调整设置。这是纯算法上带来的问题。

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

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

                      @李东岳 李老师您好,除了上面提到1.更换结构网格 2.采用新版本的openfoam 之外,您说的调整设置具体是指fvscheme里的离散格式和fvsolution中的残差、矩阵求解器、残差等吗?您能给一个大概的调整方向的建议吗?我看到您之前发的帖子中有提到twoPhaseEulerFoam等pimple类算法的的亚松弛因子就算给1也会有保证矩阵对角占优,我观察pimple系列的求解器的tut中亚松弛因子也基本都是给的1,我如果降低亚松弛因子会有助于提高twoPhaseEulerFoam计算的稳定性吗?

                      1 条回复 最后回复 回复 引用
                      • 田畔的风
                        田畔的风 讲师 最后由 田畔的风 编辑

                        离散格式用vanLeer01和linear01之类的格式可以把标量强制限制在0和1之间。

                        1 条回复 最后回复 回复 引用
                        • C
                          Caijinjin @cyberk 最后由 编辑

                          @cyberk 我这边是因为缺少了createBaffles字典,如果案例的几何边界中存在零厚度面就需要对其进行设置,我这边的搅拌釜的搅拌桨就是零厚度面

                          C 1 条回复 最后回复 回复 引用
                          • C
                            cyberk @Caijinjin 最后由 编辑

                            @caijinjin 你好,请问你是用ICEM画的网格吗?我最近也在尝试用ICEM画搅拌槽的网格,可以看我之前的发帖记录。但是我画的网格在计算的时候搅拌桨面上的epsilon值存在一些问题,初步排查原因似乎是由于周向划分数量导致桨叶处第一层网格的y+过小或者过大?想问一下你遇到过这种情况吗?可以加微信交流一下吗?我的微信是wangyu1895

                            C 1 条回复 最后回复 回复 引用
                            • C
                              Caijinjin @cyberk 最后由 编辑

                              @cyberk 抱歉,这个问题我没有遇到过,但我看了你之前用ICEM画的网格,这个六叶Ruston桨的搅拌槽的网格不太复杂的呀,这个你可以试试只画1/6个区域,然后对称过去,画完整体之后再去openfoam中用topoSet划分搅拌的区域

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