LES模拟两射流撞击
-
目前，我正在使用LES模拟两射流撞击，因为新学LES没多久，所以对相关设置不是很熟悉，希望能有大佬帮忙解答：
入口处湍射流的生成，按我目前的理解，可以有这几种方法：
1、先使用RANS计算进口的雷诺应力，然后把这些雷诺应力和附加随机扰动添加在进口。目前我有RANS的结果，但我不太清楚这个附加的随机扰动怎样生成？又怎么添加呢？
2、cfd-online上面的perturbU 小程序好像是在一个平均速度场上面添加扰动，
这里面的Ubar应该就是设置平均速度场，但我的是两股方向相反的入口射流，所以这个平均场应该怎么设置呢？
如果有大神帮忙解答疑惑，不胜感激！ -
我又来了，在使用LES模拟单射流撞击壁面时，边界条件和初始条件全部按LES的pitzdaily算例的设置，湍流模型使用动态Smagorinsky模型，时间项使用C-N格式离散，对流项TVD格式，求解器pisoFoam加一个温度输运方程，
文献中给出了沿底部壁面的nu数分布情况。
下图是我是用RANS和LES的模拟结果，可以看到使用RANS时结果吻合还可以，但使用LES时偏差非常大，疑惑。。。
-
我来还愿了。
这是新的模拟结果，LES结果可能还会调，看看会不会更好。
在修改亚格子模型，入口边界条件，离散格式等结果还是没有好转，我快要崩溃的时候，突然发现求解器中温度输运方程中，有效热扩散率alphaEff=alphat+alpha，其中alphat=nut/prt，这里面的nut在RANS中是湍流粘度，可是在LES中是亚格子粘度，会非常小，这就会导致大涡模拟计算alphat非常小，所以我直接把alphat赋值，不进行alphat=nut/prt的计算，也就是将RANS中计算得到的alphat结果直接放在0文件夹计算，就得到了上面与实验比较吻合的结果。 -
@jinlinna 在 LES模拟两射流撞击 中说：
有效热扩散率alphaEff=alphat+alpha，其中alphat=nut/prt，这里面的nut在RANS中是湍流粘度，可是在LES中是亚格子粘度，会非常小，这就会导致大涡模拟计算alphat非常小，
但这句话确实应该是这样的。
9/24