CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    请教interFoam求解器求解半潜体问题

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

      仿照DTCHull的例子计算半潜体流场,设置为2D仿真,里面把k-w模型改为k-e模型
      仿真计算域及标注:其中的RECTANGULAR为半潜体
      c2c93b9c-ec74-4f85-a7c7-b790e4312272-image.png
      计算后发现流场完全不符合物理现象:
      半潜体浮到了空中,水位一直下降

      a2e16801-fe8b-4e8d-8b6c-14e5ea525871-image.png
      0文件夹下文件:
      4fd1cea6-cc9d-493f-a34e-15925ec9ceca-image.png
      alpha.water文件:

      boundaryField
      {
         INLET
          {
              type            fixedValue;
              value           uniform 0;
          }
          OUTLET
          {
               type            variableHeightFlowRate; 
               lowerBound      0;
               upperBound      1;
               value           uniform 0;
           }
          TOP
          {
               type            inletOutlet;
               inletValue      uniform 0;
               value           uniform 0;
           }
          BOTTOM
           {
               type            zeroGradient;
           }
          RECTANGULAR
           {
               type            zeroGradient; 
           }
          frontAndBackPlanes
           {
               type            empty;
           }
      }
      

      k文件:

      dimensions      [0 2 -2 0 0 0 0];
      internalField   uniform 2e-04;   
      boundaryField 
      {
         INLET
          {
              type            fixedValue; 
              value           $internalField;
          }
         OUTLET
          {
             type            inletOutlet;
             inletValue      2e-04;  
             value           $internalField;
          }
         TOP
          {
             type            inletOutlet;
             inletValue      2e-04;  
             value           $internalField;
          }
         BOTTOM
          {
             type            kqRWallFunction;
             value           $internalField;
          }
         RECTANGULAR
          {
             type            kqRWallFunction;
             value           $internalField;
          }
         frontAndBackPlanes
          {
             type            empty;
          }
      }
      

      epsilon文件:

         dimensions      [0 2 -3 0 0 0 0];
         internalField   uniform 3.6e-03;
         boundaryField
          {
             INLET
              {
                  type            fixedValue; 
                  value           $internalField;
              }
              OUTLET
               {
                  type            inletOutlet;
                  inletValue      3.6e-03;
                  value           $internalField;
               }
              TOP
               {
                  type            inletOutlet;
                  inletValue      3.6e-03 ;
                  value           $internalField;
               }
              BOTTOM
               {
                  type            epsilonWallFunction;
                  value           $internalField;
               }
              RECTANGULAR
               {
                  type            epsilonWallFunction;
                  value           $internalField;
                }
              frontAndBackPlanes
               {
                  type            empty;
               }
          }
      

      nut文件:

      dimensions      [0 2 -1 0 0 0 0];
      internalField   uniform 9.917e-07;
      boundaryField
      {
        INLET
          {
              type            fixedValue;          
              value           $internalField;
          }
         OUTLET
          {
              type            zeroGradient;
          }
         TOP
          {
              type            zeroGradient;
          }
         BOTTOM
          {
               type            nutkWallFunction;
               value           $internalField;
          }
         RECTANGULAR
          {
               type            nutkWallFunction;
               value           $internalField;
          }
         frontAndBackPlanes
          {
              type            empty;
          }
      }
      

      p_rgh文件:

      dimensions      [1 -1 -2 0 0 0 0];
      internalField   uniform 0;
      boundaryField
      {
         INLET
          {
              type           fixedFluxPressure; 
              value          uniform 0;
          }
         OUTLET
          {
              type            zeroGradient;
          }
         TOP
          {
              type            totalPressure;
              p0              uniform 0;
          }
         BOTTOM
           {
              type            fixedFluxPressure;
              value           $internalField;
          }
         RECTANGULAR
          {
              type            fixedFluxPressure;
              value           $internalField;
          }
         frontAndBackPlanes
          {
              type            empty;
          }
      }
      

      U文件:

      dimensions      [0 1 -1 0 0 0 0];
      internalField   uniform (0.2 0 0); 
      boundaryField
       {
          INLET
           {
              type            fixedValue;          
              value           $internalField;
           }
          OUTLET
           {
               type            outletPhaseMeanVelocity;
               alpha           alpha.water;
               UnMean          0.2; 
               value           $internalField;
           }
          TOP
           {
               type             pressureInletOutletVelocity;
                value           uniform (0 0 0);
           }
          BOTTOM
           {
               type             fixedValue;
               value            uniform (0 0 0);
           }
          RECTANGULAR
           {
              type            fixedValue;
              value           uniform (0 0 0);
            }
          frontAndBackPlanes
           {
              type             empty;
           }
       }
      

      谢谢各位大神帮我解答

      向上传文件的,发现超过了500k,无法上传,请教如何把文件上传

      对于下图中所示代码
      6c74b1e4-1a7c-400c-8bff-7e1929e6c309-image.png

      在此编辑的时候是下图所示,请问如何处理呢
      6dfbf4da-c959-428f-9fe1-c2395352afc8-image.png

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

        你先把动网格关掉,看看流场是不是正确,然后再加动网格,你的alpha看起来有问题

        后面编辑问题我看下

        CFD高性能服务器 http://dyfluid.com/servers.html
        2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

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

          @东岳 谢谢老师,确实是这样,setFields从新设置就好了。太感谢老师了。

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