Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新
    1. Home
    2. Caijinjin
    C
    • Profile
    • Following 0
    • Followers 0
    • Topics 8
    • Posts 29
    • Groups 0

    Caijinjin

    @Caijinjin

    8
    Profile views
    29
    Posts
    0
    Followers
    0
    Following
    Joined Last Online

    Caijinjin Unfollow Follow

    Latest posts made by Caijinjin

    • RE: 关于multiphaseEulerFoam运算分散相体积分数问题

      @bestucan 谢谢老师,这个问题我再去认真了解一下并尝试尝试。

      posted in OpenFOAM
      C
      Caijinjin
    • RE: 关于multiphaseEulerFoam运算分散相体积分数问题

      @caijinjin 老师,这是我的离散格式和求解器设置

      ddtSchemes
      {
          default         Euler;
      }
      
      gradSchemes
      {
          default         Gauss linear;
          limited         cellLimited Gauss linear 1;
      }
      
      divSchemes
      {
          default                         none;
      
          "div\(phi,alpha.*\)"            Gauss vanLeer;
          "div\(phir,alpha.*\)"           Gauss vanLeer;
      
          "div\(alphaRhoPhi.*,U.*\)"      Gauss linearUpwind limited;
          "div\(phi.*,U.*\)"              Gauss linearUpwind limited;
      
          //"div\(alphaRhoPhi.*,Yi\)"       Gauss linearUpwind limited;
          "div\(alphaRhoPhi.*,(h|e).*\)"  Gauss linearUpwind limited;
          "div\(alphaPhi.*,f.*\)"         Gauss linearUpwind limited;
          "div\(alphaRhoPhi.*,K.*\)"      Gauss linearUpwind limited;
          "div\(alphaRhoPhi.*,\(p\|thermo:rho.*\)\)" Gauss linearUpwind limited;
          "div\(alphaRhoPhi.*,(k|epsilon|omega).*\)"  Gauss linearUpwind limited;
      
          "div\(\(\(\(alpha.*\*thermo:rho.*\)\*nuEff.*\)\*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear;
      }
      
      laplacianSchemes
      {
          default         Gauss linear corrected;
      }
      
      interpolationSchemes
      {
          default         linear;
      }
      
      snGradSchemes
      {
          default         corrected;
      }
      
      fluxRequired
      {
          default         no;
      }
      
      wallDist
      {
          method          meshWave;
          nRequired       yes;
      }
      

      求解器设置

      solvers
      {
          "alpha.*"
          {
              nAlphaCorr      1;
              nAlphaSubCycles 1;
          }
      
          droplets
          {
              nCorr                   1;
              tolerance               1e-4;
              scale                   true;
              solveOnFinalIterOnly    true;
              sourceUpdateInterval    20;
          }
      
          p_rgh
          {
              solver          GAMG;
              smoother        DIC;
              tolerance       1e-8;
              relTol          0.01;
              maxIter         20;
              minIter         2;
          }
      
          p_rghFinal
          {
              $p_rgh;
              relTol          0;
          }
      
          "U.*"
          {
              solver          smoothSolver;
              smoother        symGaussSeidel;
              tolerance       1e-7;
              relTol          0;
              minIter         1;
          }
      
          "(e|h).*"
          {
              solver          smoothSolver;
              smoother        symGaussSeidel;
              tolerance       1;
              relTol          0;
              minIter         0;
              maxIter         0;
          }
      
          "(k|epsilon|omega).*"
          {
              solver          smoothSolver;
              smoother        symGaussSeidel;
              tolerance       1e-7;
              relTol          0;
              minIter         1;
          }
      
          "f.*"
          {
              solver          PBiCGStab;
              preconditioner  DILU;
              tolerance       1e-6;
              relTol          0;
          }
      }
      
      PIMPLE
      {
          nOuterCorrectors    3;
          nCorrectors         1;
          nNonOrthogonalCorrectors 0;
          pRefPoint (0.0276 0 0);
          pRefValue 101325;
      }
      
      relaxationFactors
      {
          equations
          {
              ".*"            0.7;
          }
      }
      
      posted in OpenFOAM
      C
      Caijinjin
    • RE: 关于multiphaseEulerFoam运算分散相体积分数问题

      @李东岳 老师,我这个网格是三维的,图上展示了二维,用的结构网格,of9计算。目前发现我只要把流速设的大一点,就不会有这种问题出现,但是流速变小就不行了。另外,网格壁面不加密且数量少的时候它也不报错,自己也还在找原因中。

      posted in OpenFOAM
      C
      Caijinjin
    • 关于multiphaseEulerFoam运算分散相体积分数问题

      目前本人在用multiphaseEulerFoam求解器求解孔板中的油水两相流动,孔板架构如图
      46841618-0476-4681-8e39-5570171d6e3c-image.png

      模拟过程中水是分散相且体积分数设为了0.005,运算时发现有时候运算会报错,计算过程显示Mules算法求解相分数时,不管连续相还是分散相,最小相分数都为0,偶尔为负。请问各位老师和同学们有没有了解这一问题的?不吝赐教。
      4bd22744b0c828cfb40585cf0042cc4.jpg

      posted in OpenFOAM
      C
      Caijinjin
    • RE: 离散法求解群体平衡方程中的程序含义

      我自己暂时的理解是,程序中fjfj.phase()/fj.x()和fkfk.phase()/fk.x()指代的是j和k的颗粒数目,也就是右边公式中的Nj和Nk。但是其中的fj.phase和fk.phase是什么含义就有点不太明白了。还有就是聚并速率的量纲是m3/s的,程序里面为什么还要再乘以i颗粒体积fi.x()?这样子的话量纲就不是每秒了。

      posted in OpenFOAM
      C
      Caijinjin
    • 离散法求解群体平衡方程中的程序含义

      在离散法求解群体平衡方程中,程序文件中颗粒聚并生成项(birthByCoalescence)和文献中的给的式子形式不太一样,请问哪位大佬可以帮忙解释以下程序中为什么还要多乘以颗粒的直径?
      C.jpg coalescence.jpg
      方程里面η在程序里面应该对应gamma,指的是颗粒的分配体积,c对应程序里的coalescenceRate

      posted in OpenFOAM
      C
      Caijinjin
    • RE: 并行计算速度没有变快

      @李东岳 另外请问分块方式对计算速度有影响吗,我都是用simple方法进行分块的

      posted in OpenFOAM
      C
      Caijinjin
    • RE: 并行计算速度没有变快

      @李东岳 我台式是AMD 3600 6核, 工作站AMD EPYC 7502 32核,平时用都快不少,就这个算例基本没有加速效果,可能和网格有关系吗?

      posted in OpenFOAM
      C
      Caijinjin
    • 并行计算速度没有变快

      最近画了个釜式反应器,然后用的是结构化网格,一共50多万个六面体,使用自己台式机分四块并行计算和工作站分24块计算作比较,发现计算速度是一样的,请问这是什么原因?按理说速度最起码要提升一点啊,一点没提升是什么问题?网格是这样比较规整的b2d4d3ed-8136-4b3b-97d7-0e0b067505b4-image.png

      posted in OpenFOAM
      C
      Caijinjin
    • RE: 循环边界条件的设置

      @tens 谢谢,这种我也考虑过,但是在进口凸起部分应该会出现涡流,流场结果可能会有偏差。我在尝试用mappeed边界条件来实现。

      posted in OpenFOAM
      C
      Caijinjin