超音速下液膜演化问题--可压缩多相流求解器

[EricLiu]

在网上找到了一个叫hybridCentralFoam的可压缩两相流求解器，用于模拟激波下液滴破碎，默湍流模型为laminar，我打算用该求解器来计算超音速气流下液膜演化的问题。第一张是用求解器默认的laminar湍流模型计算的液相体积分数结果，并没有出现液膜在超音速来流下扭转破碎的现象

因此我将湍流模型换成了k-epsilon，得到了下面的结果

液膜出现了一定程度的破碎，但奇怪的是变薄了很多，感觉跟层流差异特别大。
我是按照公式算的k,nut,的大概初值，求助各位大佬

[李东岳]

你看一下你壁面附近的湍流粘度云图

[EricLiu]

@李东岳

感谢东岳老师回复，这是湍流粘度云图，液膜射流入口给的是calculated条件，然后射流下游的壁面给的是壁面函数，这个分布我看的不是很明白，您能帮忙分析一下问题吗，十分感谢东岳老师！
下面是湍流粘度的边界条件

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

internalField   uniform 0.0001;

boundaryField
{
    UPPER-WALL
    {
type            fixedValue;
      value           uniform 0.0001;
        
}

    F-INLET
    {
     type            calculated;
  value           uniform 0.0001;
         
    }

    OUTLET
    {
       type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }

   BOTTOM-WALL
    {
      type            nutkWallFunction;
        value           uniform 0.0001;
    }
    C-INLET
    {
       type            calculated;
        value           uniform 0;

    }
    frontAndBackPlanes
    {
        type            empty;
    }
}

[李东岳]

2D还是3D? 液相部分VOF处理的么

[EricLiu]

@李东岳 东岳老师，这个case是2D的，液相是VOF模型，这个求解器之前用来求解激波下的液滴演化破碎，默认湍流模型是层流，我现在换成k-epsilon，但是一直算的不合理，调了各种参数，但是没有动过fvScheme和fcSolution这两个文件，因为了解的不够，还请您指点

[李东岳]

VOF模型2D跟3D差距非常大。如果你要看真实的破碎，要换3D的。

[EricLiu]

明白了，谢谢东岳老师，我接下来试试3D模型。老师那您觉得是用层流还是湍流比较好呢