其实也不难，你用grep大概能猜出个所以然来。
grep -rn "validate"
其他都是空函数只有一个函数有点意思

template<class BasicTurbulenceModel> void Foam::ReynoldsStress<BasicTurbulencemModel>::validate() { correctNut(); }

然后我又grep 了correctNut()
常用的是kOmega 我看了下

template<class BasicTurbulenceModel> void kOmega<BasicTurbulenceModel>::correctNut() { this->nut_ = k_/omega_; correctBC; BasicTurbulenceModel::correctNut(): }

我虽然没有一行行debug， 大概猜出来就是湍流初始化。最开始给定k omega 初始值计算出第一步的nut 初始值给k omega方程，即previous time step，欢迎指正。

遇到过，忘记怎么解决的了，换个paraview低版本？

@史浩 按照老师说的用alpha识别相界面，不管我设置的那个横截面是否把网格切分，比如我设置z轴0.0002m一下为液体区域，当z轴方向50个网格时，初始计算他的液相含量是0.2；当我z轴方向51个网格时，初始计算他的液相含量变成了0.196078，导致0文件中alpha非0即1，因此并没有识别出相界面。

@chengan-wang 对

2004里面g9已经有了，如果装2012，换回g9就可以，任意切换使用

那可能是你的组分里根本就没有定义N2，你把inertSpecie N2换成inertSpecie H2O试试。

@bestucan 我去试试，感谢大佬::140:

@李东岳 东岳老师，我在使用multiphaseEulerFoam（OF8）的时候出现了这种情况，有没有可能是局部固相分数太小造成的这种发散？

謝謝老哥，後來用的歪門邪道搞定的，你説的這個方法不錯 @bestucan

@李东岳 在 关于cht的网格划分问题 中说：

writeZones

万分感谢李老师的解答，转化格式（写入命令语句）为
fluentMeshToFoam xxx.msh -writeZones，
本人OpenFoam小白，望老师多多包涵:chouchou:

可以啊 这个跟几何没关系啊

@bestucan 万分感谢！

@李东岳 对的，大概有十几公里，那谢谢东岳老师，我再使劲努努力钻研一下。个人觉得根据文献植入算法太难了，虽然能看懂离散和求解过程，但是想在openfoam里面实现觉得难度挺大的，当时老师您带我们植入的单相可压缩的，目前这是两相可压缩，我努力努力，再次感谢东岳老师的持续关心:xinxin:

问题已解决，勾选paraView中的“Use VYkPolyhedron”选项后交错就会消失

@bestucan 我看了目前用到这个 -I. 的地方，包括reactingParcelFoam跟sprayFoam，都有子文件夹，比如sprayFoam下的egineFoam，应该跟compile子求解器有关

@李东岳 东岳老师是想诱惑我再报一次您的辅导班:141:

@cresendo 我也是看了亿遍~~~~~

在opemfoam中有一个内聚力，可以使粒子之间相互黏结。它是通过粒子之间的碰撞的能量转换实现的。其中通过一个叫cohesionEnergyDensity的系数来决定它的强度。这个密度越大，黏结越明显。希望我说的这个是你想问的。

@史浩 VOF和LES都是高度网格依赖类求解器。主要是LES只能捕获大于网格尺度的涡，因此，网格越小，LES越能捕获到更小的涡。类似的，VOF捕获到的界面只能大于网格尺度，如果计算域存在气泡，如果网格越小，VOF可以捕获到的气泡越小。

@李东岳 实在是有点困难~请问老师如何让两个场（如图片中的a和vDotAlpha）对应网格的值相乘，这不是内积也不是外积，用forAll循环的话，会报assignment of read-only location错误。

forAll(a,celli) { vDotcAlphal[celli] = a[celli]*vDotcAlphal0[celli]; vDotvAlphal[celli] = a[celli]*vDotvAlphal0[celli]; }

图片2.png

麻烦帮忙改正一下。想要规定thetaBoundary为一个常数值，然后对theta（相分数）进行对比，总是实现不了。谢谢。

const Field<scalar>& theta();
const Field<scalar>& thetaBoundary();
thetaBoundary = {0.5};
if (theta < thetaBoundary)