@coolhhh 感谢！

@xjwang

CFDRFG的附录code还有个问题，求解$K$，不是按照论文中的原理计算的。他先用$PQ$没有代入$S(f)$的值，算了一组$K$，然后所有坐标点都用了同一组$K$。而根据理论，应该是要根据每个坐标点的非均匀$PQ$，求对应该点的$K$。DSRFG如果生成非均匀，理论上也要这么计算。
739858e5-ac4f-4736-94bd-9f4ec9f07a9d-image.png

关于CDRFG方法，他们的code里是对每一个点都生成一组PQK

最近看到一篇JFM文章《Guo, H., et al., An efficient and low-divergence method for generating inhomogeneous and anisotropic turbulence with arbitrary spectra. Journal of Fluid Mechanics, 2023. 970: p. A2.》，用的就是三维能谱，不是频谱，生成非均匀湍流场也是low-divergence。

@luofq-sysu 这个靠谱 差不多 嗯

通过vorticity求解涡量； 写个程序求解$\nabla w$
2.1 在这个程序里面，$\nabla w$的边界场$\cdot\bfn$再乘以粘性系数

说简单也简单，只不过繁琐了一些，好像没有几个命令就能出来的东西

可以尝试用functionObject中的surfaces提取物体表面的压力

1693360823669.png还需要添加其他限制条件，具体可以查看这个网站的说明
https://openfoamwiki.net/index.php/Parameter_Definitions_-_dynamicMotionSolverFvMesh

@linhan-ge 我觉得也应该是这样，但我还没理解李老师说的，您理解了吗？

如果你的d是个dimensionedScalar，那We是个dimensionedScalar

好的，我试一下1912版本的试一下

@Banbor 温度是个场。你要修改哪部分温度数据

我看到了一个bug修复，可以用新版看看还有没有负数

https://github.com/OpenFOAM/OpenFOAM-dev/commit/1ae52882748922f4e02c92ba52301c416954e426

@李东岳 好的老师，谢谢您

@李东岳 多谢李老师，可惜我现在人在国外向日葵貌似用不了。但是我跟您说下最终测试结果。我们可以使用./Allwmake -j -k把那一步跳过去。这样编译出来的还是能用的。从日志上来看，貌似是这个FoamtoFireMesh这里
3b01921f-39c2-4d31-a36d-b8310fcd497a-image.png
而且，如果我不下载代码包，直接用官网的命令在线安装是没有任何报错的。
我把报错的日志发给您看下，您看看是否有头绪。

链接：https://pan.baidu.com/s/1RcbLj7qFiBwl3yK4eEhDjA
提取码：wo84

OpenFOAM热解反应中某一种反应物的质量随时间的变化情况，应该怎样去计算呢？还有controlDict下时间步长adjustableRunTime，库朗数从0.5修改到0.8，结果却发生很大变化，这是为什么呢？

@yhdthu 大佬别闹，求赐教:140:

请问大家:

为什么执行foamToVTK指令之后，流场的cell总数、points总数都发生了变化？其中的原理是什么？

使用paraview加载的.vtk数据是否为体心变量？谢谢！

15c3027b-3697-4f52-abf4-80866407053f-image.png

@李东岳 谢谢李老师的帮助，虽然目前没能解决，但是找到原因了，目前能单核计算也挺好的至少可以用，等之后有能力优化代码再分享解决办法

@李东岳 好的，谢谢老师

@李东岳 谢谢李老师，interFoam里粘性应力计算我是按上面那几行公式算的，所有的场都输出来验证过，最后得到的devRhoReff和rho()*turb.devReff是一样的；受力的计算方法不统一的话我试试别的公式，找个差别小的就用了。

@xpc um, 我做surfaceLIC还没遇到过这个问题，也就不知道咋整了，看看别的大佬怎么说