q-DNS计算槽道流遇到了一些问题,求大神们指点
-
@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} $$