LES模拟求助
-
如题,我正在算一个自由射流,LES模型采用的Smagorinsky,同时计算了RANS,几何模型就是一个三维圆柱。LES的计算结果显示,在过渡区,射流轴向中心线的密度值总是比实验值高,如下图所示。麻烦各位大佬,可以帮忙分析一下可能的原因吗?
-
@strive_dut LES相比RANS的nut、alphat都要小,普适性来讲是这样的
-
@strive_dut 可以分享下结果吗?调残差管用吗?
-
@strive_dut 两个点需要关注下,一是网格是否足够小,二是入口的条件是否给的合适。我之前也用OpenFOAM做过LES,有点难搞,不知道你这个结果怎么样了,可以分享下么?
-
东岳老师,各位同行,我近期也一直在修改参数重新计算,我抽空整理一下,放到论坛,还望各位大佬多指点。
-
最近做了一些算例测试,更新一下结果,欢迎各位提出宝贵意见!
入口脉动通过InflowGenerator产生,主要修改了入口的湍流强度和湍流普朗特数。有关残差问题,我尝试缩小至1e-9,结果并没有变化。
另外,图中的红色曲线是两篇不同文献的结果,同样使用Smag模型,看起来差距很大。这里附上参考文献2的一些模拟结果,感觉是Smagorinsky模型不太合适呢。
240万网格模拟结果:
330万网格模拟结果:
参考文献2的模拟结果(SM:Smagorinsky;WL:Wale;VR:Verman;SMD:Dynamic Smagorinsky):
-
@strive_dut 你这个就是射流,然后物性剧烈变化么?有传热的过程么?看330万网格的结果,好像湍流强度影响很大,蓝色那条线其实和实验对的还可以呀。话说湍流强度在inflowGenerator里面怎么体现的,我记得它只能给雷诺应力的6个分量吧。
-
@cccrrryyy 您好!我算的是可压缩射流,物性变化比较剧烈,具有高密度梯度。射流与环境存在传热过程。
蓝色的线,射流的液核长度比较短,跟实验差别比较大,其他区域还算理想。我正在测试一个Prt=1.0的Case,看看结果。
您说的对!是给的雷诺应力,雷诺应力计算的时候是(UI)^2,这个I应该指的湍流强度吧?