请教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湍流模型时,自由液面出现这样的现象,应该做何调整,来使得自由液面尽量的光滑贴近实际呢? 有没有推荐的资料,我好好研究看看这部分。有点无头绪,万分感谢!!! -
@波流力 您好!我找了下,基本都是在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