横向圆柱绕流湍流计算不收敛问题

OpenFOAM
1

image
边界也是根据表格中的设定的，设置内容都是计算或文章中得来，不知道哪里有问题。

2

看不到你的模型，边界，没法做评价啊

3 
  object      epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -3 0 0 0 0];

internalField   uniform 0.02547;

boundaryField
{
   
    INLET
    {
         type            inletOutlet;
         inletValue      uniform 0.02547;
         value           uniform 0;
    }
    OUTLET
    {
         type            inletOutlet;
         inletValue      uniform 0.02547;
         value           uniform 0;
    }
    FRONT_AND_BACK
    {
        type  zeroGradient;
       
  
    }
    LEFT
    {
       type  zeroGradient;
  
      
    }
    RIGHT
    {
      type  zeroGradient;
   
    
    }
    INTERFACE_INTERFACE
    {
       type  epsilonWallFunction;
       value uniform 0.02547;
    }
   
object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0.03375;

boundaryField
{
   
     INLET
    {
        type            inletOutlet;
        inletValue      uniform 0.03375;
        value           uniform 0;
    }
    OUTLET
    {
         type            inletOutlet;
         inletValue      uniform 0.03375;
         value           uniform 0;
    }
    FRONT_AND_BACK
    {
        type zeroGradient;
        value uniform 0;
  
    }
    LEFT
    {
        type zeroGradient;
        value uniform 0;
      
    }
    RIGHT
    {
       type zeroGradient;
        value uniform 0;
    
    }
    INTERFACE_INTERFACE
    {
        type kqRWallFunction;
        value uniform 0.03375;
    }
  object      nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    
   INLET
    {
        type           calculated;
        //value          uniform 0;
    }
    OUTLET
    {
        type            calculated;
        // inletValue      uniform 0;
         //value           uniform 0;
    }
    FRONT_AND_BACK
    {
        type           zeroGradient;
        value           uniform 0;
  
    }
    LEFT
    {
       type           zeroGradient;
       value           uniform 0;
      
    }
    RIGHT
    {
      type            zeroGradient;
      value           uniform 0;
    
    }
    INTERFACE_INTERFACE
    {
      type            nutkWallFunction;
      value           uniform 0;
    }
  object      nuTilda;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    INLET
    {
       type   zeroGradient;
    }
    OUTLET
    {
        type   zeroGradient;
    }
    FRONT_AND_BACK
    {
         type   zeroGradient;
    }
    LEFT
    {
        type   zeroGradient;
    }
    RIGHT
    {
       type   zeroGradient;
    }
    INTERFACE_INTERFACE
    {
       type   zeroGradient;
    }
  object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
   
    INLET
    {
        type           zeroGradient;
        value          uniform 0;
    }
    OUTLET
    {
         type           fixedValue;
         value          uniform 5.409e3;
       
    }
    FRONT_AND_BACK
    {
        type           zeroGradient;
  
    }
    LEFT
    {
        type            zeroGradient;
      
    }
    RIGHT
    {
        type            zeroGradient;
    
    }
    INTERFACE_INTERFACE
    {
        type           zeroGradient;
        value          uniform 0;
    }
 object      pointMotionU;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    INLET
    {
        type           fixedValue;
        value          uniform (0 0 0);
    }
    OUTLET
    {
         type           zeroGradient;
         value          uniform (0 0 0);
       
    }
    FRONT_AND_BACK
    {
        type           fixedValue;
         value          uniform (0 0 0);
  
    }
    LEFT
    {
        type           fixedValue;
        value          uniform (0 0 0);
      
    }
    RIGHT
    {
        type           fixedValue;
         value         uniform (0 0 0);
    
    }
    INTERFACE_INTERFACE
    {
        type           fixedValue;
        value          uniform (0 0 0);
    }
 object      R;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform (0 0 0 0 0 0);

boundaryField
{
   
     INLET
    {
        type kqRWallFunction;
        value uniform ( 0 0 0 0 0 0 );
    }
    OUTLET
    {
         type kqRWallFunction;
        value uniform ( 0 0 0 0 0 0 );
    }
    FRONT_AND_BACK
    {
        type kqRWallFunction;
        value uniform ( 0 0 0 0 0 0 );
    }
    LEFT
    {
       type kqRWallFunction;
        value uniform ( 0 0 0 0 0 0 );
    }
    RIGHT
    {
        type kqRWallFunction;
        value uniform ( 0 0 0 0 0 0 );
    
    }
    INTERFACE_INTERFACE
    {
         type kqRWallFunction;
        value uniform ( 0 0 0 0 0 0 );
    }
 object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 3 0);

boundaryField
{
     INLET
    {
        type           fixedValue;
        value          uniform (0 3 0);
    }
    OUTLET
    {
         type           zeroGradient;
         value          uniform (0 0 0);
       
    }
    FRONT_AND_BACK
    {
        type           fixedValue;
         value         uniform (0 0 0);
  
    }
    LEFT
    {
        type           fixedValue;
        value          uniform (0 0 0);
      
    }
    RIGHT
    {
        type           fixedValue;
         value         uniform (0 0 0);
    
    }
    INTERFACE_INTERFACE
    {
       type        fixedValue;
       value       uniform (0 0 0);
    }
4

image
另外 该模型在虚拟机啊下运算买虚拟机会自动重启

5

我觉得你的边界有点混淆，没有弄清楚。

仔细检查下边界吧，最好能列个表，看看边界有没有问题？

在这里列表，可以参考这个格式 http://xianbai.me/learn-md/article/extension/table.html

这里有个示例和效果：

你 copy 过来就可以修改下

6

边界跟这个表格是一一致的，除了u下面的wall我用的是固定值；
另外一个案例的内侧流体和管道的计算速度还不错，我把u中wall的固定值改成了movingwallvelocity就可以算了

7

那就有问题啦

case 都不同，怎么能用一样的边界？你列一下你的边界条件吧。列出来，你就能看出问题在哪了。

8 
Time = 0.002, iteration: 3
Maximal accumulated displacement of interface points: 0.098614664
GAMG:  Solving for cellMotionUx, Initial residual = 0.95370137, Final residual = 0.00045599299, No Iterations 1
GAMG:  Solving for cellMotionUy, Initial residual = 0.96348561, Final residual = 0.00014276369, No Iterations 3
GAMG:  Solving for cellMotionUz, Initial residual = 0.99999652, Final residual = 0.00015547833, No Iterations 3
GAMG:  Solving for cellMotionUx, Initial residual = 0.0007328271, Final residual = 4.7399663e-07, No Iterations 6
GAMG:  Solving for cellMotionUy, Initial residual = 0.00012143445, Final residual = 6.0776397e-07, No Iterations 4
GAMG:  Solving for cellMotionUz, Initial residual = 0.00012980668, Final residual = 7.3704718e-07, No Iterations 4
Evolving flow model
Courant Number mean: 0.076610348 max: 0.67708466 velocity magnitude: 5.8490433
DILUPBiCG:  Solving for Ux, Initial residual = 0.00095234971, Final residual = 9.7218489e-08, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.0039401773, Final residual = 4.2292268e-07, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.00031208718, Final residual = 2.3436977e-08, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.033030465, Final residual = 4.3804072e-07, No Iterations 12
GAMG:  Solving for p, Initial residual = 0.017281581, Final residual = 5.5453573e-07, No Iterations 9
time step continuity errors : sum local = 1.7850388e-09, global = 2.4516785e-10
GAMG:  Solving for p, Initial residual = 0.011900619, Final residual = 4.2335046e-07, No Iterations 10
GAMG:  Solving for p, Initial residual = 0.0022941043, Final residual = 7.7017365e-07, No Iterations 7
time step continuity errors : sum local = 2.4450769e-09, global = -8.8047757e-10
GAMG:  Solving for p, Initial residual = 0.0013032211, Final residual = 8.3179127e-07, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.00017690006, Final residual = 7.5331011e-07, No Iterations 4
time step continuity errors : sum local = 2.3982512e-09, global = -3.3739974e-10
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00021256906, Final residual = 1.7128835e-08, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.00035567439, Final residual = 5.5818143e-08, No Iterations 2
Setting traction on solid patch
Total force (fluid) = (-89.990217 -1.8384034 -0.0028548367)
Total force (solid) = (90.238228 1.8243108 0.0028836326)
Evolving stress model
DICPCG: Solving for D, Initial residula = 0.00042239694, Final residual = 0.00084485725, No outer iterations = 25
Max relative residual = 0.0078576518, Relative residual = 5.6923678e-06, enforceLinear = 0
Current fsi relative residual norm: 1
[1]+  Segmentation fault      (core dumped) fsiFoam > log

这是管内无水的fsi，库朗数正确怎么也报错，segmentation fault？
在另一个算例中，我把管当成实心的做fsi却没有报错。

9 
Evolving stress model
DICPCG: Solving for D, Initial residula = 0.00042239694, Final residual = 0.00084485725, No outer iterations = 25

初始的残差才 0.0004，算完了残差还增加到 0.0008，这算不算有问题呢？