@李东岳 感谢东岳老师提供的思路，我在代码中把ShihQuadraticKE改了个名，目前在源文件中用Allwmake编译成功了。
不过孩子还是很不明白怎么单独拿出来编译这类非线性湍流模型，看来暂时只能通过这种方式代替了。

@bestucan 好的，谢谢老师

@wangfei9088 好的，谢谢，可以了

@chen_hao
解决了，改成了无量纲，就避免出问题：

volSymmTensorField Sij(twoSymm(tgradU())/omega_); volTensorField Wij(2.0*skew(tgradU())/omega_);

@cfdngu 收费也没办法，我之前找了一圈也没有找到比这个好用的

@tens 感谢！这正是我想要的

如果已知流场内某一点的坐标（a，b，c），想通过插值获得这一点的流场物理量，可以使用interpolate（）函数，这一函数的定义如下： template<class Type> Type Foam::AveragingMethods::Basic<Type>::interpolate ( const barycentric& coordinates, const tetIndices& tetIs ) const { return data_[tetIs.cell()]; }

代码中coordinates、tetIs分别表示什么呢？
如果已知点（a，b，c）如何获取coordinates、tetIs呢？
谢谢！

@李东岳 十分感谢

@李东岳 正在学习，谢谢李老师！

@李东岳 好的李老师，里边的数据已经拷出来了。多谢。

@cccrrryyy 感谢感谢！给我指了一个好的修改方向！我做过RANS模型的对比，确实影响不大，但我看文献里有的是Smagorinsky算的很好，比如文献1，有的就是算的误差比较大，比如文献2，一头雾水:136: 还是得自己测试呀！

@sibo 我靠，前辈，我也是做结冰的

不仅是64以上不线性。10多核以上就不线性了。

@bestucan谢谢老师不吝赐教[抱拳]

[General] Export full spreadsheet from Paraview -- CFD Online Discussion Forums (11_11_2022 8_57_13 PM).html

@bestucan 谢谢这位老师了，我的这个问题的确是李老师说的那种情况。

https://www.cfd-online.com/Forums/openfoam-solving/147135-how-switch-off-combustion-reaction-reactingfoam.html

https://github.com/vitst/dissolFoam

上来就自定义，完了还要跨版本号编译代码。把门远远甩在身后的人干起来也是横垄地里拉车，一步一个坎。:chigua2:

先照着视频教程之类的，实现自定义函数什么的。

提供的信息很难判断问题在哪

s.component(0) = 123 s.component(1) = 456

看看这样行不，这样行的话，那就是ux、vx的值有问题了