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



  • 仿照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



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

    后面编辑问题我看下



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


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2