求助:Openfoam中的LES/DNS如何输出湍流耗散率
-
目前在用icoFoam算圆柱绕流,算例和 http://www.dyfluid.com/icoFoam.html 中设置基本相同。
现在想生成无量纲化平均应变张量
OpenFOAM利用fieldAverage可以直接生成平均速度与雷诺应力
。湍动能 即平均速度梯度可以后处理“postProcess -func “grad(UMean)””得到。
那耗散率在openfoam里面可以怎么求呢?
看定义是可以通过改编fieldAverage求。但是试了一下没有成功。。请问一下能通过什么方法求得呢。 -
最近在备课,准备了很多湍流的东西,
这个定义不是真正的epsilon的定义,真正的定义应该是 。其中所以在写代码的时候,大体就是:
- 通过求解器求解速度
U
, - 通过fieldAverage获得
UMean
- 二者相减是
U'
- 通过
U'
获得 - 求解瞬态量
,对瞬态量 做平均。这个也可以在paraview里面通过滤镜来实现:Filter-temporal-statistics
大体上,需要做两次平均过程。第一次平均获得uMean,第二次获得
定义里面的平均。这样可以获得解析的 。同时还要加上 。可以参考Wilcox那本湍流书的109页。
我算了算
的算例,左边是湍流动能耗散率,右边是速度。看起来是符合物理的。就是没有数据可以对。 - 通过求解器求解速度
-
@李东岳 李老师您好~
求得时均速度之后,就可以用时均速度和当前时间步的脉动速度求瞬时的脉动速度应变率张量S’,可以直接对S’求双内积(&&)嗷,为什么要用这一步呢? -
@xiexiaoyang 嗯?没懂
-
@李东岳 !李老师,就是这样:
求得脉动速度后,直接做双内积,然后把这个标量场存下来,再用fieldAverage做时间平均。 -
@xiexiaoyang 我求的是这个:
,不是 。这两个有点区别,但是结果应该差不多。你试试channel395,看看结果如何? -
目前我也在做相同的工作,我困惑的是这里面的粘度项是否需要进行平均,如果看定义的话是不需要的,但是李老师的这个做法是不是把粘性的平均也考虑进去了
-
运动粘度用sutherland定律来计算,也就是在温度高的地方粘性大,相乘的话也就是在温度高的地方耗散率相对比较大,这样是合理的吗?
我重新核查了一下
的定义,这个看起来就是层流粘度。如果这样说的话,看起来就是跟粘度相关的变量。这个我目前也不能确定,你最好查一下相关的文献,是不是这样的。也告诉我一下https://www.cfd-online.com/Wiki/Turbulence_dissipation_rate
你这个确实是个好问题,需要严肃的研究一下
发布 20 的 12