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

Zhoujq

仿照DTCHull的例子计算半潜体流场，设置为2D仿真，里面把k-w模型改为k-e模型
仿真计算域及标注：其中的RECTANGULAR为半潜体

计算后发现流场完全不符合物理现象：
半潜体浮到了空中，水位一直下降

0文件夹下文件：

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，无法上传，请教如何把文件上传

对于下图中所示代码

在此编辑的时候是下图所示，请问如何处理呢

李东岳

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

后面编辑问题我看下

Zhoujq

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