Smagorinsky模型系数问题
-
@李东岳 李老师, 这两个链接在讨论
of
中的Smagorinsky
模型时,也类似讨论of中为何植入的是 ,相比理论公式看似缺少 系数。但其实是一致的,这是因为of中把 加到了前面的系数中。下面截图中的理论公式符号代表: ,但OF公式代表: ,在CFD-online上的帖子因为把符号写成一样导致了混淆。
Smagorinsky model details -- CFD Online Discussion Forums
[转载]有关OpenFoam中LES smagorinsky模型的推导_馨予_新浪博客
湍流模型我不太懂,在想
dynSmagorinsky
是否类似于CFD-online的讨论,是否也是把系数 加到了前面的系数中?也就是cD(D)
可能比dynSmagorinsky
理论公式多了 系数?void dynSmagorinsky::updateSubGridScaleFields(const volSymmTensorField& D) { nuSgs_ = cD(D)*sqr(delta())*sqrt(magSqr(D)); nuSgs_.correctBoundaryConditions(); }
-
@李东岳 这个公式结论是对的,不用加
。CFD-online讨论帖子初始问题混淆点是把OF计算的 当做 ,而理论公式的 表示的是 。OF是单独计算了 ,为与理论对应, 在前面的系数里计算再相乘。Smagorinsky SGS model in OpenFOAM | CFD WITH A MISSION这个帖子推导就没有混淆符号,统一用
表示
-
-
-
@李东岳 在 Smagorinsky模型系数问题 中说:
你的意思就是,OpenFOAM里面植入的应该是
,其中而不是
是这样么
https://bugs.openfoam.org/view.php?id=816
这个确实是个bug,不过OpenFOAM新版已经把动态Smagorinsky删掉了,所以无关紧要了。
-
我按照
的方式植入进入了。从云图来看没啥区别。希望有大佬来详细验证一下是否均衡。下面这个文件,可以在OpenFOAM-10下进行编译,步骤:wmake cd pitzDaily ./Allrun
注意算例中的
的值忘记改了,应该改为0.065
个人疑问是植入:
那么式 (1d) 中的k是怎么计算的?这会不会就是两个方法的主要差别?@coolhhh 应该是从
反推出来的, -
@cfdem小白 可以看3#的帖子,这是《Moin, P. and J. Kim, Numerical investigation of turbulent channel flow. Journal of fluid mechanics, 1982. 118: p. 341-377.》建议取
Cs=0.065
,然后根据of2206中planeChannel
算例的设置,LES采用Smagorinsky,提到设置Ce=1.048,Ck=0.0265463553,能够Updated to give Cs = 0.065
。我也用这个系数计算过Channel Flow,结果跟实验吻合还比较好
19/47