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

[李东岳]

在一些文献里面我看有人这么用，更简单：
$$
\varepsilon=\frac{\mu_{eff}^3}{\rho_c^3(C_s\Delta)^4}
$$

Simulation of coalescence breakup and mass transfer in polydisperse multiphase

[李东岳]

@李东岳" class=" img-fluid img-markdown">

没个老铁来讨论讨论咩

[XieXiaoyang]

@李东岳 李老师您好~
求得时均速度之后，就可以用时均速度和当前时间步的脉动速度求瞬时的脉动速度应变率张量S’，可以直接对S’求双内积（&&）嗷，为什么要用这一步呢？

[李东岳]

@xiexiaoyang 嗯？没懂

[李东岳]

加一个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

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