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}} $) 就是：

计算 $\overline{u_x}$
计算 $u_x^{\prime}$
计算 $\overline{u_y}$
计算 $u_y^{\prime}$
计算每个时间点的 $ {u_x^{\prime}u_y^{\prime}} $
计算 $ \overline{u_x^{\prime}u_y^{\prime}} $

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

兔子安东尼

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