CFD中文网

流

1 最近在尝试在dynamicKEqn模型中使用不同的滤波器，参照CFD-online上的帖子和LESfilter.H的设置更改了滤波器。
按理来说，根据dynamicKEqn头文件的说明如此可以更改用于亚格子动能的滤波器
但是出现了如下报错：

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#5  Foam::LESModels::dynamicKEqn<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::Ck(Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
#6  Foam::LESModels::dynamicKEqn<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correctNut(Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#7  Foam::LESModels::dynamicKEqn<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correctNut() at ??:?
#8  ? at ??:?
#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10  ? at ??:?
Floating point exception (core dumped)

令人感到奇怪的是对于Smagorinsky模型更改滤波器就能正常运行。

2 另一处令我感到困惑的地方在于OpenFOAM中盒式滤波和高斯滤波居然使用同一参数进行开启，而这两者是完全不同的滤波器，OpenFOAM是如何做到区分两者的？

流

@random_ran
感谢您这么长时间以来的耐心解答与帮助，问题我已解决，所有的症结都出在边界条件上，这一方法适用于矩形网格和圆形网格。具体的细节我不太方便透露，如果您对此有兴趣的话，我们可以加下联系方式（邮箱或微信等皆可）具体的聊下
最后再次对您一直以来的帮助致以诚挚的感谢，如果日后我因为算力不足转向DES，到时还请您多加指教！

流

@random_ran
再次对您详细认真的回复表示感谢，我现在在使用二维的圆形网格计算（虽然根据我使用dynamic kEqn进行二维计算的众多经验来看，这不可能有什么好结果），因此，暂时没有新的计算结果与您讨论，仅就以下三点和您进行讨论：

1. 您是否做过turbulentInlet对脉动值的影响的测试？尤其是对于Cd’
1. 您能简略的说下one equation eddy-viscosity model和dynmaic kEqn的相同与不同之处么？
1. 您说您的计算资源不足以进行LES计算，但是我看到您有在本帖的计算中说过您曾经使用144核和96核进行计算，这里指的是物理核心，还是线程（逻辑核心）？我在您的论文种看到您也做了Re10e5左右的算例，但都使用的是DES，一百余核也不足以支撑在Re10e5左右的LES计算么？

流

@random_ran
这么短感谢您在这么短的时间内，给出如此详细的答复，但是我还有几点想要和您详细讨论下：

1，更正，上次提问中我误将“计算时长”打为“计算步长”，我实际想要问的是，我总共只计算了25个周期，负Cd是否与计算时间长短有关。
2，在您的论文中，我注意到您在近壁使用了Spalding的壁面函数，我想问下您的壁面函数是在nut文件中设置的么？使用壁面函数后，近壁y+为多少，周向节点与展向节点有是否需要相应的做出调整？
3，我想和您确认下圆形网格的划分，您是边界条件是这样设置的么？
4，使用dynmaic kEqn模型的时候k的初始值应当怎样设定？对计算结果有何影响？我虽然知道k是亚格子动能，但是不知道如何在inlet处设置，我现在在出口处暂时使用fixedValue格式的湍动能作为替代。
5，我注意到您在论坛中发过如下帖子

我查找了4.1的dynamickEqn的头文件

通过对比论文我确认了OpenFOAM中源代码与头文件中声明的模型一致
但是您在论文中使用的LES为

您在毕业论文中更换湍流模型的原因是什么？是因为计算结果的好坏么？还是其他的原因？

参考文献：Kim, W.-W., & Menon, S. (1995). A new dynamic one-equation subgrid-scale model for large eddy simulations. In 33rd Aerospace Sciences Meeting and Exhibit.

流

@random_ran 你好，我初步尝试使用dynamic K Eqn对Re=2e5情况下的D=1m三维圆柱绕流进行了计算，计算了25个周期，但是升阻力系数曲线始终没有出现任何明显的周期现象，甚至Cd出现了负值，Cd，Cl频谱曲线也十分混乱，

经过初步检查，我认为有以下几种可能：
1，计算时长过短
2，离散格式设置有误
3，Z方向节点过少
请问您在调试过程中是否遇到同样的问题，是怎样解决的？如果没有，您对解决这一问题有何建议与意见？先行谢过！
计算域如下

近壁核心区是一个4D*4D的矩形，进口到中央矩形的距离为5D，出口到中心矩形的距离为11D，两侧距离中心矩形距离都是3D。
近壁核心区的布置如下图

离散格式如下

边界条件如下

CFD中文网

流散然

帖子