FLuent使用两方程模型如何统计流场的雷诺应力分布

[solarzeng]

各位好！我最近在模拟树冠对风场的模拟时，在两方程模型中添加了一些动量和湍流的附加源项，想在在后处理的时候，想统计一下某一个截面的湍流应力，比如U'V'和U'W'，不知道大家有什么方法。已经尝试过在源项UDF中添加一个用户自定义内存（C_UDMI）显示湍流应力，但是失败了，特求助大家。先谢谢各位了。

[李东岳]

不懂Fluent友情帮顶
不过你做的很有工业意义！

[solarzeng]

@东岳 谢谢东岳老师，在论文中看到了有相关的统计分析，可是困惑于不知如何在商业软件中实现。

[李东岳]

在两方程模型中添加了一些动量和湍流的附加源项，

只研究过森林的这部分东西，后处理那面没玩过了~

[solarzeng]

@东岳 请问您对带有城市植被的风场模拟有什么高见么？最近在对比不同的湍流模型和不同形式的源汇项对风场模拟的影响。

[李东岳]

只是在壁面跟你一样植入了湍流动能和动量的耗散，增加了湍流耗散。计算没问题。就没有再玩了。之前是个工业咨询项目。所以说你这个很有工业前景！具体的验证你们来做吧

[ibelief]

report 可以么？我没有尝试过

[solarzeng]

@ibelief hi,report里面没有找到相应的变量

[random_ran]

只要得到了速度场随时间变化的数据，你自己就可以求出雷诺应力。

雷诺分解就是把瞬时速度重写成两项：时均速度与速度的波动(瞬时速度减去时均速度)。用公式写出来就是：
$$\begin{array}{}\text{(1)}& u=\overline{u}+u^{\mathrm{\prime }},v=\overline{v}+v^{\mathrm{\prime }},w=\overline{w}+w^{\mathrm{\prime }}.\end{array}$$

把雷诺应力的定义及其展开形式：
$${\tau }_{ij}^{\mathrm{\prime }}\equiv \rho \left[\begin{array}{ccc}\overline{u_x^{\mathrm{\prime }}{u}_x^{\mathrm{\prime }}}& \overline{u_x^{\mathrm{\prime }}{u}_y^{\mathrm{\prime }}}& \overline{u_x^{\mathrm{\prime }}{u}_z^{\mathrm{\prime }}}\\ \overline{u_y^{\mathrm{\prime }}{u}_x^{\mathrm{\prime }}}& \overline{u_y^{\mathrm{\prime }}{u}_y^{\mathrm{\prime }}}& \overline{u_y^{\mathrm{\prime }}{u}_z^{\mathrm{\prime }}}\\ \overline{u_z^{\mathrm{\prime }}{u}_x^{\mathrm{\prime }}}& \overline{u_z^{\mathrm{\prime }}{u}_y^{\mathrm{\prime }}}& \overline{u_z^{\mathrm{\prime }}{u}_z^{\mathrm{\prime }}}\end{array}\right],$$

其中 (举个例子)，$\overline{u_x^{\mathrm{\prime }}{u}_y^{\mathrm{\prime }}}$的意思是 $u_x^{\mathrm{\prime }}{u}_y^{\mathrm{\prime }}$ 的时均值。

具体的算法 ($\overline{u_x^{\mathrm{\prime }}{u}_y^{\mathrm{\prime }}}$) 就是：

1. 计算 $\overline{u_x}$
2. 计算 $u_x^{\mathrm{\prime }}$
3. 计算 $\overline{u_y}$
4. 计算 $u_y^{\mathrm{\prime }}$
5. 计算每个时间点的 $u_x^{\mathrm{\prime }}{u}_y^{\mathrm{\prime }}$
6. 计算 $\overline{u_x^{\mathrm{\prime }}{u}_y^{\mathrm{\prime }}}$

将最后的结果乘以流体的密度，就得到了雷诺应力。 如果有错，还请大家指正。

[兔子安东尼]

@random_ran 您好，我在fluent用LES算湍流流动，也需要求雷诺应力，想详细请教您，QQ965968052，诚心请教