求助：Openfoam中的LES/DNS如何输出湍流耗散率

[李东岳]

加一个LES湍流动能的计算，这个相对简单。OpenFOAM中可以在计算平均的过程中计算$\overline{\bfU}$，同时，这个量：$\overline{\bfU'\bfU'}$也可以计算出来。那解析的湍流动能$k_{res}$就是$0.5\cdot\mathrm{tr}(\overline{\bfU'\bfU'})$。模化的湍流动能$k_{mod}$一般在LES模型里面计算。

[XieXiaoyang]

@李东岳 !李老师，就是这样：

求得脉动速度后，直接做双内积，然后把这个标量场存下来，再用fieldAverage做时间平均。

[李东岳]

@xiexiaoyang 我求的是这个：$\epsilon=2\nu \overline{|\bfS'|^2}$，不是$\epsilon=\nu \overline{\dfrac{\partial u_{i}^{\prime}}{\partial x_{k}} \dfrac{\partial u_{i}^{\prime}}{\partial x_{k}}}$。这两个有点区别，但是结果应该差不多。你试试channel395，看看结果如何？

[李东岳]

@李东岳 Wilcox说这两个的区别在2%左右..

[Dingcy]

@李东岳 这个式子里的delta是大涡里过滤的那个delta吧，与网格尺寸有关，能不能提示一下输出它的方法呀

[快乐的山炮]

目前我也在做相同的工作，我困惑的是这里面的粘度项是否需要进行平均，如果看定义的话是不需要的，但是李老师的这个做法是不是把粘性的平均也考虑进去了

[李东岳]

哪个粘度项？$\nu$么？这个是常数

[快乐的山炮]

@李东岳 李老师，想问一下这里的计算应该用层流粘性还是湍流粘性

[快乐的山炮]

@快乐的山炮 我算的是LES

[快乐的山炮]

@李东岳 李老师好，还是没太弄清楚，按照定义来说我用的运动粘度进行计算，运动粘度用sutherland定律来计算，也就是在温度高的地方粘性大，相乘的话也就是在温度高的地方耗散率相对比较大，这样是合理的吗？如果是用湍流粘度来计算则完全不同 对此非常疑惑，请李老师能否有时间解答一下，另外我做的的LES燃烧

[快乐的山炮]

@快乐的山炮 并且按照定义只对了 做时间平均 没有对粘性项做平均

[李东岳]

@快乐的山炮 你有没有看到你们领域那面的epsilon怎么定义，把公式发来看看

[快乐的山炮]

@李东岳 老师好，我暂时还没找到领域内有定义的书籍，主要还是参考poinsot 2005和wilcox 2006这两本书

[李东岳]

运动粘度用sutherland定律来计算，也就是在温度高的地方粘性大，相乘的话也就是在温度高的地方耗散率相对比较大，这样是合理的吗？

我重新核查了一下$\epsilon$的定义，这个看起来就是层流粘度。如果这样说的话，看起来就是跟粘度相关的变量。这个我目前也不能确定，你最好查一下相关的文献，是不是这样的。也告诉我一下

https://www.cfd-online.com/Wiki/Turbulence_dissipation_rate

你这个确实是个好问题，需要严肃的研究一下