@Amadeus 好啊好啊！加个微信吧，我微信号：wx-200115qsw:146:

@李东岳 感谢东岳老师解答

@李东岳 好的。以前在CFD-DEM的时候要几千个颗粒，DNS的话应该得缩短一下管道长度，算少一点颗粒量。

@李东岳 谢谢您

还可以考虑GPU运行，不使用cpu运行，具体可以看一下这个WX公众号的内容：TFluid多物理场仿真

@学流体的小明 好的，谢谢大佬，准备用starccm+试试画多面体网格，感觉六面体这么搞没得玩

@李东岳 好的好的，谢谢老师回复。

各位大佬，我要模拟一个圆柱体，上表面逐步加载，上面表面没有流入和流出。水从下表面流出，算出来的结果很奇怪，要求计算10S，前3S都算正常，水向下流，但算到5S就停止了，最后一步水反向流出了。以下是流速的计算结果，帮忙看看，边界条件设置是否有问题，为什么结果怪怪的，谢谢。

circle.jpg

circle2.jpg

压力边界条件
inlet
{
type fixedValue;
value uniform 0.01;
}

pipe
{
type zeroGradient;
}
outlet
{
type zeroGradient;

}

流速边界条件
inlet
{
type zeroGradient;

}

pipe
{
type noSlip;
}
outlet
{
type inletOutlet;
inletValue uniform (0 0 0.1);
value uniform (0 0 0.1);
}

我有几个问题啊，inlet上表面压力为什么是固定值啊? 我要是跑着跑着里面的压力变得大于入口压力了，那岂不是反流了?
此外，速度边界条件的Outlet为什么给uniform (0 0 0.1);这个值啊?

wallShearStress可以计算tau_w。然后用来算y+就好啦。

感谢🙏，虚拟机总共给的4个，CSDN也查到说是并行太多，单J无后缀就是拉满，等会儿用J2重新编

@东方白杨 不清楚什么原因，你看在程序最后再重新计算一遍保存，看看结果有无区别

GradCenPot = - fvc::grad(CenPot); GradCenPot.write();

@Voynich 目前看到的处理都是把k忽略掉，这个帖子OpenFOAM 不可压缩湍流模型的 divDevReff 函数 | Giskard's CFD Learning Tricks中的2个观点，其实本质都是k被忽略掉，只是不知清楚k忽略对结果影响具体是多大，我也没具体搜过有无人做过这方面的比较

这个是自己写的求解器，所以有可能差别确实很大

https://www.cfd-china.com/topic/7301/waves2foam的一键安装解决方案?_=1738765131135

应该是偏算法吧，线下课应该也是以算法为主吧？那我还是学C++吧，到时候看赶不上这个三月底的线下课了

@李东岳 感谢大佬的回复！！
您的意思是我需要创建cellZones来划分这些区域吗，那么关于这些区域的控制方程以及其他属性我需要如何进行分配，之前是在求解器内部进行布置的。

我尝试过了kOmegaSSTIDDES在OpenFOAM-v2212版本中可以使用，但是你上面字典里IDDESDelta全写成
IDDEESDelta了，多了个E会报错，你的字典去掉E就可以正常运行了，0文件内包含P、U、k、nut、nuTilda、omega

@李东岳 一个普通的打工研究生罢了。。。

div(phi,U) Gauss limitedLinear 0;

我还第一次见有人这样用。我下周找个小网格测试一下。

div(phi,U) Gauss linear;

另外你这个结果看起来不太对，有震荡