q-DNS计算槽道流遇到了一些问题,求大神们指点
-
@东岳 谢谢东岳老师
这个meanVelocityForce的源项我用过,但是和我想要的不太一样,meanVelocityForce也是通过添加一个压力梯度来控制平均速度的,但它给添加的压力梯度项是随时间步不断计算调整的,我想添加一个简单的常数的压力梯度。不过我确实可以通过修改meanVelocityForce来达到这个目的,今后我会试一下这个。
另外,我那个平均速度的代码就是从meanVelocityForce中copy过来的,稍微改了一下下让他能在求解器中运行,现在来看运行的还可以。而且这个平均速度只是为了在log中输出一个当前流场的平均速度,方便我检测求解状况,不会影响迭代求解的,所以还是不太清楚结果中的问题出在了哪里。ORZ -
我想用q-DNS模拟一个Re_tau=180的槽道湍流,网格与这篇链接文本文献中基本相同。
如果你要实现这个q-DNS模拟槽道湍流,我觉得你可以试试meanVelocityForce,之前我也是用你这种植入代码的方式计算,后来发现meanVelocityForce更简单并且两者效果相同。
这个平均速度只是为了在log中输出一个当前流场的平均速度
我当时植入的代码平均速度是给定的,因为你完全可以不植入任何代码直接用icoFoam做直接模拟,但是速度回往下掉, 因此给定一个固定的平均速度。如果你只是为了输出,那应该没问题吧
-
-
@李东岳 李老师,我看很多文献里,槽道湍流DNS计算时,除了时间平均,还可以做展向和流向的空间平均,在OpenFOAM自带LES Ret=395的算例中,有一个postChannel脚本可以改来用,个人觉得用了这个效果挺好的hhh~
-
@xiexiaoyang 有道理
-
@XieXiaoyang 您好,我想请教一下,像这种 DNS / LES 槽道流的结果,我想画速度剖面,进行正则化时候的摩擦速度
和粘性长度 是怎样计算的?
我当前的做法是:使用postChannel(就像channel395那个算例一样)将空间时间都做平均之后,可以得到 这样一个分布,然后计算壁面上的速度曲线的斜率,计算壁面摩擦应力,计算摩擦速度:
但算出来的摩擦速度和预设的相差很大。所以最后的曲线是这样的
算了很久,对数率那里也无法和理论解对上。计算结果的整个曲线也与理论解的形状不符合。
感觉我使用的方法不太对,想请问一下您的做法。 -
-
@SSSSK 建议看下我这个帖子。我自己使用的时候也经常估计不好压力梯度数值和具体的平均速度值,流场速度一直掉,其实就是压力梯度给的动力小于流场的耗散。所以我先用设定平均速度的方法,先让程序算,稳定之后的那个压力梯度就可以知道了。之后再直接用固定压力梯度算。
https://www.cfd-china.com/topic/6248/les直流槽道边界层模拟-如何得到正则化速度u-以及正则化坐标y/3?_=1698989899271
12/23