请教divSheme中设置对溃坝自由液面形态的影响, RNG k-e湍流模型

Sloan

dam break算例，使用RNG k-epsilon湍流模型，最开始fvScheme中divScheme如下：

divSchemes
{
    div(rhoPhi,U)   Gauss linearUpwind grad(U);
    div(phi,alpha)  Gauss vanLeer;
    div(phirb,alpha) Gauss linear;
    div(phi,k)      Gauss upwind;
    div(phi,epsilon) Gauss upwind;
    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

自由液面在撞击到右侧壁上升过程出现裂隙现象，随后的跌落翻卷也不尽如意，如下图：

使用limitedLinear, 对divSchem做如下改动：

divSchemes
{
    div(rhoPhi,U)   Gauss linearUpwind grad(U);
    div(phi,alpha)  Gauss limitedLinear01 1;
    div(phirb,alpha) Gauss interfaceCompression;
    div(phi,k)      Gauss upwind;
    div(phi,epsilon) Gauss upwind;
    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

液面中裂隙现象减弱了些，但依然存在，液面翻卷跌落不是很好，如下图所示：

如果把divSchemes中第一项也是用linearlimited:

divSchemes
{
    div(rhoPhi,U)   Gauss limitedLinearV 1;
    div(phi,alpha)  Gauss limitedLinear01 1;
    div(phirb,alpha) Gauss interfaceCompression;
    div(phi,k)      Gauss upwind;
    div(phi,epsilon) Gauss upwind;
    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

液面在上升过程中形态还可以，但随后的跌落，翻卷过程显得有些不自然，如下图：

所以向大家请教，在使用RNG k-epsilon湍流模型时，自由液面出现这样的现象，应该做何调整，来使得自由液面尽量的光滑贴近实际呢? 有没有推荐的资料，我好好研究看看这部分。有点无头绪，万分感谢！！！

波流力

@Sloan 算例发上来看看，我也在做这方面的研究，可以帮你看看

波流力

把你要对比的文献也发上来

波流力

我做过三维溃坝跟实验结果的对比，效果还是很好的

Sloan

@波流力嗯嗯，好的，算例在这里damBreakRNGkepsilonmodified.zip

Sloan

@波流力现在做的这个没有对比的文献，但是有好多dam break的粒子，不管是paper还是视频，跟我算出来的不一样，我的计算里自由液面不对劲。。。

CFD中文网

请教divSheme中设置对溃坝自由液面形态的影响, RNG k-e湍流模型