CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    气液两相 没有纯气体

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

      DZ44GBIONA}7`WR})4VI(2V.png
      GU0962NTFA1OXF5XX{B5.png
      上面是加热壁面,下面全都是多孔区域(孔隙率0.73),然后我产生的液相分布是这样的~ 几乎没有纯气体区域~甚至在固液界面处的液相含量反而高~我已经调过不同的毛细力模型,发现不是毛细力的影响~
      现在还可以做的是增大相变系数~但是我相变系数一增大,就会在开始相变的时候发散~

      有没有老师可以分析一下 可能是哪里导致的问题~

      solvers
      {
          "alpha.liquid.*"
          {         
              cAlpha          1;
              nAlphaCorr      2;
              nAlphaSubCycles 2;
              MULESCorr       yes;
              nLimiterIter    5;       
              solver          smoothSolver;
              smoother        symGaussSeidel;
              tolerance       1e-7;
              relTol          0;
          }
          rho
          {
              solver           diagonal;
              tolerance        1e-7;
              relTol           0.1;
          }
          rhoFinal
          {
              $rho;
              tolerance        1e-7;
              relTol           0;
          }
          mDotSmearFinal
          {
              solver           PCG;
              tolerance        1e-7;
              preconditioner   DIC;
              relTol           0.00;
              smoother         DIC;
          }
          "(U|T)"
          {
              solver           smoothSolver;
              smoother         symGaussSeidel;
              tolerance        1e-7;
              relTol           0.;
              minIter          15;
              maxIter          50;
          };  
          "(U|T)Final"
          {
              $U;
              tolerance        1e-7;
              relTol           0;
              maxIter          50;
          }
          p
          {
              tolerance        1e-8;
              solver           PCG;
              preconditioner   DIC;
              relTol           0.001;
              smoother         DIC;
          }
          pFinal
          {
              $p;
              tolerance       1e-07;
              relTol          0;
          }    
          "pcorr.*"
          {
              $p_rgh;
              tolerance        1e-9;
              relTol           0;
          };
      }
      PIMPLE
      {
          momentumPredictor   yes;
          nCorrectors         5;
          nNonOrthogonalCorrectors 0;
          pRefCell            0;
          pRefValue           1e5;
      }
      relaxationFactors
      {
          alpha.liquid   0.2;
          U              0.5;
          p              0.5;
      }
      
      ddtSchemes
      {
          default         Euler;
      }
      gradSchemes
      {
          default         Gauss linear;
      }
      divSchemes
      {
          div(phi,alpha)  Gauss vanLeer01;
          div(phirb,alpha) Gauss vanLeer01;
          //div(phirb,alpha) Gauss interfaceCompression;
          //div(rhoPhi,U)      Gauss vanLeerV;
          //div(rhoPhi,U)      Gauss limitedLinearV 1;
          div(rhoPhi,U)      Gauss linearUpwind;
          div(rhoPhiByEps,U)  Gauss linearUpwind;
          div((mu*dev2(grad(U).T()))) Gauss linear;
          div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
          div(((rho*nu)*dev2(T(grad(U))))) Gauss linear;
          div(rhoCpPhi,T) Gauss vanLeer;
          div((interpolate(cp)*rhoPhi),T) Gauss vanLeer;
      }
      laplacianSchemes
      {
          default         Gauss linear corrected;
      }
      interpolationSchemes
      {
          default         linear;
      }
      snGradSchemes
      {
          default         corrected;
      }
      
      1 条回复 最后回复 回复 引用
      • First post
        Last post