q-DNS计算槽道流遇到了一些问题，求大神们指点
-
@东岳 谢谢东岳老师
这个meanVelocityForce的源项我用过，但是和我想要的不太一样，meanVelocityForce也是通过添加一个压力梯度来控制平均速度的，但它给添加的压力梯度项是随时间步不断计算调整的，我想添加一个简单的常数的压力梯度。不过我确实可以通过修改meanVelocityForce来达到这个目的，今后我会试一下这个。
另外，我那个平均速度的代码就是从meanVelocityForce中copy过来的，稍微改了一下下让他能在求解器中运行，现在来看运行的还可以。而且这个平均速度只是为了在log中输出一个当前流场的平均速度，方便我检测求解状况，不会影响迭代求解的，所以还是不太清楚结果中的问题出在了哪里。ORZ -
-
-
@李东岳 李老师，我看很多文献里，槽道湍流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} $$ -
