q-DNS计算槽道流遇到了一些问题，求大神们指点

[悬铃神木]

@东岳 是的，求解本身没有用到平均速度 ，也没有规定平均速度，给它一个压力梯度（实际相当于增加了一个体积力），让摩擦力自己跟压力梯度平衡去，平均速度会在某个值附近上下波动。

[李东岳]

我试试你这种方法

[XieXiaoyang]

@李东岳 李老师，我看很多文献里，槽道湍流DNS计算时，除了时间平均，还可以做展向和流向的空间平均，在OpenFOAM自带LES Ret=395的算例中，有一个postChannel脚本可以改来用，个人觉得用了这个效果挺好的hhh~

[李东岳]

@xiexiaoyang 有道理

[vbcwl]

@李东岳 东岳老师您好，我之前了解到meanVelocityForce适合槽道这种几何，然后不可压缩，如果我想做圆环柱计算域的传热，然后采用Boussinesq的话，这种计算域以及轻微的热膨胀还可以使用meanVelocityForce么？

[vbcwl]

@悬铃神木 请问你的rms拟合效果怎么样呢？

[学流体的小明]

@XieXiaoyang 您好，我想请教一下，像这种 DNS / LES 槽道流的结果，我想画速度剖面，进行正则化时候的摩擦速度$u_{\tau }$和粘性长度${\delta }_{\nu }=\nu /u_{\tau }$是怎样计算的？
我当前的做法是：使用postChannel（就像channel395那个算例一样）将空间时间都做平均之后，可以得到$u(y)$这样一个分布，然后计算壁面上的速度曲线的斜率，计算壁面摩擦应力，计算摩擦速度：
$$\begin{array}{}\text{(1)}& u_{\tau }=\sqrt{\frac{{\tau }_w}{\rho }}=\sqrt{\nu \frac{\mathrm{\partial }u}{\mathrm{\partial }y}}\end{array}$$
但算出来的摩擦速度和预设的相差很大。所以最后的曲线是这样的

算了很久，对数率那里也无法和理论解对上。计算结果的整个曲线也与理论解的形状不符合。
感觉我使用的方法不太对，想请问一下您的做法。

[李东岳]

你得壁面函数用的哪个壁面函数？标准壁面函数还是Spalding那个

[学流体的小明]

nut 在壁面上使用壁面函数 nutUSpaldingWallFunction，会出现上面那样的速度分布图
但另外一个$R{e}_{\tau }$也等于1000的算例，得到的速度分布也差不多
我等等重开一个帖子说明下问题

[SSSSK]

@悬铃神木 您好，我想请教一下，我也是用了这个工具进行打乱，然后给了一个固定的压力梯度，但他随着时间推移速度还是在不断变小。想问一下您这个也是先变小然后再增大吗？大概算了多久才能和文献对上，压力梯度是给的多少啊

[学流体的小明]

@SSSSK 建议看下我这个帖子。我自己使用的时候也经常估计不好压力梯度数值和具体的平均速度值，流场速度一直掉，其实就是压力梯度给的动力小于流场的耗散。所以我先用设定平均速度的方法，先让程序算，稳定之后的那个压力梯度就可以知道了。之后再直接用固定压力梯度算。
https://www.cfd-china.com/topic/6248/les直流槽道边界层模拟-如何得到正则化速度u-以及正则化坐标y/3?_=1698989899271

[SSSSK]

@学流体的小明 感谢回复，还有一点我想问的就是固定梯度和固定速度这两个方法最终算出来的结果是一样的吗？如果给定固定梯度我该怎么确定我算的摩擦雷诺数是我想要的那个值呢，是在计算后提取摩擦速度来重新计算然后根据算出来的来调整吗

[学流体的小明]

@SSSSK
两种方法估计摩擦应力：
第一种是经验公式：
$${\mathrm{Re}}_{\tau }\approx 0.09{\mathrm{Re}}^{0.88}$$
$$\mathrm{Re}=U_m\left(2h\right)/\nu$$
$h$为槽道半高，$U_m$为槽道平均速度。
第二种是完全发展槽道流中的计算：
$$\frac{\mathrm{\partial }p}{\mathrm{\partial }x}\times h={\tau }_w$$

[SSSSK]

@学流体的小明 感谢大佬，我是按第二种给的压力梯度为1，结果一直掉，可能是我算的时间太短了，想问一下如果是用pertubU给了一个初始场之后，用icofoam算准dns，大概得算多少秒才能得到完全发展得湍流？

[学流体的小明]

@SSSSK
我用给定槽道平均速度的方法算摩擦雷诺数550的槽道，根据这个摩擦雷诺数的变化曲线，感觉得到10s以后吧，你也可以画这样一条曲线看看。

[SSSSK]

@学流体的小明 感谢，我这算了三天才算了不到六秒，感觉算的有问题就停了，我再继续算算看