请教divSheme中设置对溃坝自由液面形态的影响, RNG k-e湍流模型
-
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湍流模型时,自由液面出现这样的现象,应该做何调整,来使得自由液面尽量的光滑贴近实际呢? 有没有推荐的资料,我好好研究看看这部分。有点无头绪,万分感谢!!! -
把你要对比的文献也发上来
-
我做过三维溃坝跟实验结果的对比,效果还是很好的
-
@波流力 嗯嗯,好的,算例在这里damBreakRNGkepsilonmodified.zip
-
@cccrrryyy 谢谢建议!是1m长的tank,这样的长度网格多少就可以了呢?也在搜索一些信息和帖子看如何能调节好,看到有这样的建议,我就调节对流项了。不知道您还有其他的建议或是学习资料吗,感谢!
-
damBreakRNGkepsilonmodified.zip
自己加密网格算算 -
@波流力 您好!我找了下,基本都是在YouTube上的视频,好多是带obstacle的,paper里的我找了下,是不同时刻的image,能看出来特征,主要是翻卷和溅落,供您参考。链接:http://www.cjcp.org.cn/fileup/HTML/2017-6-641.htm
以及这篇paper里的图片:
、https://www.sciencedirect.com/science/article/abs/pii/S0029801818304049
还有这个:https://www.mdpi.com/2073-4441/11/5/1110/htm
9/36