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

[solarzeng]

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

[李东岳]

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

[solarzeng]

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

[李东岳]

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

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

[solarzeng]

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

[李东岳]

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

[ibelief]

report 可以么？我没有尝试过

[solarzeng]

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

[random_ran]

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

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

把雷诺应力的定义及其展开形式：
$$ \tau_{ij}^{\prime} \equiv \rho
\begin{bmatrix}\overline{u_x^{\prime}u_x^{\prime}}&\overline{u_x^{\prime}u_y^{\prime}} &\overline{u_x^{\prime}u_z^{\prime}} \\ \overline{u_y^{\prime}u_x^{\prime}}&\overline{u_y^{\prime}u_y^{\prime}} &\overline{u_y^{\prime}u_z^{\prime}} \\ \overline{u_z^{\prime}u_x^{\prime}}&\overline{u_z^{\prime}u_y^{\prime}} &\overline{u_z^{\prime}u_z^{\prime}}
\end{bmatrix},
$$

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

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

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

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

[兔子安东尼]

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