还可以考虑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();

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

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

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

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

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

div(phi,U) Gauss limitedLinear 0;

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

div(phi,U) Gauss linear;

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

你网格变密了，库朗数也得跟着改吧，是自适应时间步长吗? 可以考虑把库朗数限制在0.5以下?

@Amadeus 是的，一个是Cp一个是Cv，但是目前看的各种版本推导的温度方程，有的用Cp有的用Cv不太一致。。。。

@Amadeus @aurora16 感谢各位回复，最后发现是由于对流-扩散方程中，离心力项太大了(高达几百量级)，而其它项还在10的量级，离散求解方程时，量级差异太大了，导致发散。计算还是要保证方程中各项量级差异不能太大

不需要potentialFoam

我之前算过高亚声速叶栅，用rhoSimpleFoam求解器收敛不了，换成LUSGS就很好收敛

如题，现在已经计算完成了，在paraview中也可以查看云图，怎么绘制最大值随时间变化曲线