@kanjunao 您好，我是用的动网格

@rhzhang 谢谢您

本人小白，正在复刻国外一个大型水槽实验，原实验对于某双色波的组次共造了18个波，总模拟时长200s，由于两种波具有不同周期，故两种波应该在不同的时间点结束造波，结束造波后波高的时间序列应该会出现自然衰减的效果。

请问各位大佬，如何修改waveproperties文件能使两个波在不同的时间点结束造波？

对于waveProperties文件的设置，一开始就简单设置为

waveType bichromaticFirst; relaxationNames ( inlet ); initializationName inlet; pName p_rgh; inletCoeffs { waveType bichromaticFirst; period1 5.26; period2 6.45; ......

后来又试过另一种造波类型

waveType combinedWaves; relaxationNames ( inlet ); ... inletCoeffs { waveType combinedWaves; combinedWaveNames ( wave1 wave2 ); relaxationZone { ... } } wave1Coeffs { waveType stokesFirst; period 5.26; ... startTime 0; endTime 94.68; } wave2Coeffs { waveType stokesFirst; period 6.45; startTime 0; endTime 116.1; }

但是从计算后输出的log文件和结果来看，这个starttime和endtime并没有正确运行（可能是被直接忽略了），不知道是根本就没有这个功能呢还是我设置的不对？
有做过类似案例的老师吗:135:

另外一个就是把rapidcfd重写，目前这个适用于of23，cuda8，但是因为我不懂cuda，所以这个活我也干不了。

https://exafoam.eu/benchmarks/

请大佬关注

（暂且不谈显性LES）

最近闲着没事跑了个LES，发现隐性LES能算出湍流，隐性滤波LES湍流直接磨平了。顺着这个思路，想到一个问题：隐性LES（一些文献里面说是no model LES），在不需要LES湍流粘度项的时候，就可以预测比较好的结果。那么为什么还需要在隐性LES湍流模型上，添加显性的湍流粘度项？

类似下图，no model LES与Smagorinsky已经都很好了，那添加Smagorinsky模型的意义是？

JCP 345 (2017) 565–595

捕获.JPG

@a602657563 您好，麻烦问一下您提到的这个算例的文件能麻烦发一份给我么，我装的openfoam11中没有这个案例了，十分感激

我搜了一下，还没看到谁把OpenFOAM和HOS-NWT耦合到一起并开源代码

@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();