q-DNS计算槽道流遇到了一些问题，求大神们指点
-
-
9月CFD算法编程课： http://dyfluid.com/class.html
需要帮助debug算例的看这个 https://cfd-china.com/topic/8018
-
@李东岳 李老师，我看很多文献里，槽道湍流DNS计算时，除了时间平均，还可以做展向和流向的空间平均，在OpenFOAM自带LES Ret=395的算例中，有一个postChannel脚本可以改来用，个人觉得用了这个效果挺好的hhh~
-
@xiexiaoyang 有道理
-
@XieXiaoyang 您好，我想请教一下，像这种 DNS / LES 槽道流的结果，我想画速度剖面，进行正则化时候的摩擦速度${u_{\tau}}$和粘性长度${\delta_{\nu}=\nu/u_{\tau}}$是怎样计算的？
我当前的做法是：使用postChannel（就像channel395那个算例一样）将空间时间都做平均之后，可以得到${u(y)}$这样一个分布，然后计算壁面上的速度曲线的斜率，计算壁面摩擦应力，计算摩擦速度：
\begin{equation}
{u_\tau } = \sqrt {\frac{{{\tau _w}}}{\rho }} {\rm{ = }}\sqrt {\nu \frac{{\partial u}}{{\partial y}}}
\end{equation}
但算出来的摩擦速度和预设的相差很大。所以最后的曲线是这样的
算了很久，对数率那里也无法和理论解对上。计算结果的整个曲线也与理论解的形状不符合。
感觉我使用的方法不太对，想请问一下您的做法。 -
学 学流体的小明 被引用 于这个主题
-
-
-
@SSSSK
两种方法估计摩擦应力：
第一种是经验公式：
$${{\mathop {\rm Re} \nolimits} _\tau } \approx 0.09 {{\mathop {\rm Re} \nolimits} ^ {0.88}}$$
$${\mathop{\rm Re}\nolimits} = {U_m}\left( {2h} \right)/\nu $$
$h$为槽道半高，$U_m$为槽道平均速度。
第二种是完全发展槽道流中的计算：
$$ \frac{{\partial p}}{{\partial x}} \times h = {\tau _w} $$ -
