Your browser does not seem to support JavaScript. As a result, your viewing experience will be diminished, and you have been placed in read-only mode.
Please download a browser that supports JavaScript, or enable it if it's disabled (i.e. NoScript).
我又来了，在使用LES模拟单射流撞击壁面时，边界条件和初始条件全部按LES的pitzdaily算例的设置，湍流模型使用动态Smagorinsky模型，时间项使用C-N格式离散，对流项TVD格式，求解器pisoFoam加一个温度输运方程， 文献中给出了沿底部壁面的nu数分布情况。 下图是我是用RANS和LES的模拟结果，可以看到使用RANS时结果吻合还可以，但使用LES时偏差非常大，疑惑。。。
为什么不在用buoyantPimpleFoam？
@李东岳 求解器这个是因为，作者的动量方程没有考虑温度的影响，而且在做RANS的时候，我是用的ｓｉｍｐｌｅＦｏａｍ添加温度输运方程，看结果吻合的还可以，所以做LES就直接用了ｐｉｓｏＦｏａｍ添加温度输运方程，
@李东岳我待会试一下ｂｕｏｙａｎｔｐｉｍｐｌｅＦｏａｍ看结果会不会有改善，请问李老师除了求解器有问题，还可能是哪里的原因导致的呢？
我来还愿了。 这是新的模拟结果，LES结果可能还会调，看看会不会更好。 在修改亚格子模型，入口边界条件，离散格式等结果还是没有好转，我快要崩溃的时候，突然发现求解器中温度输运方程中，有效热扩散率alphaEff=alphat+alpha，其中alphat=nut/prt，这里面的nut在RANS中是湍流粘度，可是在LES中是亚格子粘度，会非常小，这就会导致大涡模拟计算alphat非常小，所以我直接把alphat赋值，不进行alphat=nut/prt的计算，也就是将RANS中计算得到的alphat结果直接放在0文件夹计算，就得到了上面与实验比较吻合的结果。
你开始给的alphat是一个场？
@李东岳 都是定义了一个标量场，但原本是
这样子计算的，没搞清楚RANS中nut和LES中nusgs的区别，还是跟RANS中一样让入口处的alphat直接用calculated边界条件计算，就导致热扩散率一直很小
@jinlinna 在 LES模拟两射流撞击 中说：
有效热扩散率alphaEff=alphat+alpha，其中alphat=nut/prt，这里面的nut在RANS中是湍流粘度，可是在LES中是亚格子粘度，会非常小，这就会导致大涡模拟计算alphat非常小，
但这句话确实应该是这样的。
@李东岳 不好意思李老师，我没明白您的意思，您是说不应该修改温度输运方程吗？alphat应该用nusgs/prt计算吗？ 我现在把温度输运方程修改为这样，然后alphat直接使用RANS的结果计算，这样做是不对的吗？
alphat应该用nusgs/prt计算吗？
是的。你为何不用buoyantPimpleFoam呢，已经处理好了，并且温度T本身不是守恒量，只有beta版求解器才求解温度，求解能量更好
然后你这个算力好像没有考虑密度，不需要考虑温度引起密度变化的浮力？
@李东岳 我这个新的使用的是buoyantBoussinesqPimpleFoam求解器，这个求解器就是有一个温度输运方程，我直接修改的它
那你不用修改了，直接用就可以了。alphat用壁面函数了么
@李东岳 没有用，壁面处用的zeroGradient，网格画的很密，壁面处yplus小于1
@李东岳 李老师，但是我直接使用buoyantBoussinesqPimpleFoam求解器计算的话，撞击区的Nu数就一直很小，无论我改变什么条件，如果说我的那种做法不可行的话，请问我可以上传一下我的算例麻烦您看一下吗？不知道您是否方便？
网格数多少？
因为是二维的，所以网格数好像才十万左右，挺小的
用LES算二维？有没有想过用3D？
二维只是为了复现这篇文献，后面自己研究的话会做三维的
@jinlinna 您好，想问下为什么alphat初场会对计算结果有明显影响呢，在0时刻后，alphat不还是用nutSGS/prt计算么
