dynamicKEqn模型更换滤波器奇怪报错

流散然

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是如何做到区分两者的？

hoversoar

您好，关注到你也在做高雷诺数下的圆柱绕流。看到你之前的Cd出现负值的情况，不知道解决了没有？
可否加个联系方式探讨一下？
十分感谢！

东岳

@流散然在 dynamicKEqn模型更换滤波器奇怪报错中说：

Smagorinsky模型更改滤波器就能正常运行。

Smagorinsky模型并不需要显性的定义滤波函数

hoversoar

请问这个滤波为什么会报错有什么结果吗? 我现在也碰到报错了，但是其他mesh又可以用，在怀疑是不是mesh的问题？