LES定义入口速度的问题(DSRFG方法)
-
@霜染丹枫 换了台电脑回复你^-^
能不能贴一下你不同位置的速度的紊流度以及功率谱的图片?
还有一个问题是,比如你的计算域顺风向长度是15,你目标结构应该放在5-8左右,这样紊流度等衰减都不是很大。
我之前犯的错是,入口的参数搞错了,紊流积分尺度这些设置错了。
以下是在公众号下边的留言:
抱歉在论坛上提问问题后没有及时更新后续。
很惭愧,这个问题后来我发现是参数设置错了,生成风场的参数是按照缩尺的参数选取的,而DSRFG这个方法,本来是对网格尺寸不敏感的,但因为参数选错了,所以导致了在入口处引起了divergence-free问题。至于后续问的那个沿顺风向紊流度减小的问题,确实是会减小,不过不是太过于明显,因为我错误选取了参数,导致顺风向的衰减程度特别大。之所以会问在两个面定义成inlet,是因为看过一篇文章,Kim, Y., Castro, I. P., & Xie, Z. T. (2013). Divergence-free turbulence inflow conditions for large-eddy simulations with incompressible flow solvers. Computers & Fluids, 84, 56-68.
看楼上有位同学或老师回复过这个问题,他用的CDRFG方法,这个方法是产生inhomogeneous的inflow的,所以在入口处不满足divergence-free,而且文章中感觉有一处错误(也许是我理解不够吧),所以就没有用这个方法。
在入口处加尖劈或者粗糙元,我个人觉得不太好控制想要的风场参数(如紊流积分尺度等),而且由此导致的网格加密会增加计算负担,虽然DSRFG方法也是很耗费计算资源。
-
@xjwang 感谢回复,我的计算域如图所示,几何尺寸:X=1.2m、Y=0.2m、Z=0.4m。展向和轴向为均匀网格。空气沿x轴正向进入,平均速度是1.2m/s。入口采用速度入口边界条件,出口为压力出口边界条件,展向是周期性边界条件(即Z方向),上部和底部为壁面,采用无滑移边界条件。三个方向网格节点:64X64X128,垂直壁面方向采用非均匀网格,首层网格高度△y=0.0001m,保证Yplus<1。目前程序中湍流积分尺度设置为Ls=0.014(大致估算的,目前看来是有问题的)。Huang文章说这个参数很重要,给了三种设置的方法,但是我一直没引起重视。。。。您这个参数是如何来设置的呢?谢谢!