你先备份一下。

然后你wclean一下，看一下这些东西会不会消失。消失的话再编译一下

@李东岳 好滴，感谢李老师，我编一下

@李东岳 在 VOF方法中如何准确计算自由表面曲率？ 中说：

这个方法在函数calculateK()里面，可直接调用。计算结果在3232网格上能得到光滑的结果，但不对。在6464网格上出现了数值振荡。

VOF结果在网格比较细的时候结果不对？这个结论很有挑战性

李老师，我用interIsoFoam算了一下，还是不对。在32$\times$32的网格上可以得到光滑的结果，但是和解析解不一样；在64$\times$64的网格上结果出现了振荡。interfoam也跑了一下，结果和interisofoam不一样，也不对。。。（这个算例的设置参照和解析解在主贴的参考文献[1]中）

算例在这里，大家有兴趣可以跑一下。要是算例设置有问题也可以指出来。另外，我是在v2006环境下跑的(这个链接可能需要注册坚果云，我没找到其他的上传方法:xinlei: )
capillaryWave

@李东岳 并行文件夹里面的2e-08文件夹是有pointProcAddressing文件的，但是后面算出来的2e-02，4e-02就没有这个文件了。所以reconstructPar说找不到这个文件，报错了。

@exthan 您好，我刚刚也下载了OF7版本的makeAxialmesh,但好像没找到编译文件，没有Allmake文件或者wmake文件，请问您是怎么编译的呢d27bf862-45db-4120-8ba1-050e18d46128-image.png

@bestucan 谢谢老师，我试了下，还是不行~

他是由于none量 他把它删了，看来我的问题就是那个omega，只要用

const volScalarField& omega_ = ( this->mesh_.objectRegistry::template lookupObject<volScalarField>("omega") );

就会找不到，估计是结构问题，难道只能重写了。

@李东岳 李老师您好，感谢您提出的意见，现在只是算的有些慢但是似乎问题已经解决了，谢谢。33.png

@jacob-z

mpirun --oversubscribe -np 20 xxFoam -parallel

应该是进程超限锁死，需要添加申请允许超限，就可以了

还有种方法，通过一个工具将OF中颗粒坐标转换成旧版本的格式，然后就能被paraview和tecplot等读取了
https://github.com/blueCFD/lagrangianExtraFunctionObjects

@李东岳 多谢李老师，我的物理问题中速度是一定会求解的。所以放在泊松方程中能够更加稳定，避免震荡。
明白了，多谢您

你那个77行的ttc，是一种简单的声明，换一种比较全的写法，就是那种IOobejct都写上的那种，然后命名一下，估计可以

求问各位，以往的OpenFOAM Workshop 的论文集在哪可以找到？在这个网站里没有呢？https://openfoamworkshop.org/previous-workshops/
我想找一个OpenFOAM Workshop 2016年的一些报告，请问有没有知晓如何获取的？:quwan:

@李东岳
针对 ncompressible::turbulenceModel 模板特化后的 ShihQuadraticKE 模型 ，仔细研究了一下其怎么实现非线性项的实现，发现其是在 NonlinearEddyViscosity.H 中定义 nonlinearStress_ 项，然后用 ShihQuadraticKE 作为子类继承父类 NonlinearEddyViscosity 以便在 ShihQuadraticKE.C 中直接引用omega_。

但由于我修改的是SSTkomegabase.H( 继承于TurbulenceModel类 )与 linearViscousStrss.H ( 继承于BasicturbulenceModel原则类 ) 没有特例化，不能简单类似上面的ShihQuadraticKE 模型。需要从新仿造框架，是一项大工程！

我的策略是：在现有的文件上直接硬植入，简单很多，毕竟涉及到湍流模型的类错综复杂，针对找不到omega，按照我的分析应该是对的，运行求解器时linearViscousStrss.C先运行此时并未运行SSTkomegabase.C 所以没有omega场注册，找不到他，所以我在creatfield.H中IOobject一下omega场，目前可以算了，结果还未检查，不知道我的理解有没有问题。

Data-Create Zone-Rectangular提取出想截取的位置
d2469623-64c0-45c3-952d-a6339a3b560b-image.png
但是提取出来的Zone没有数据，Data-Interpolate根据原数据对创建的Zone进行插值就可以
abd5cf18-caaf-4a49-a528-3b9afd696f35-image.png
这样就得到一个名为Rectangular zone的Zone
758ce629-1262-40bf-83b2-1477d15a474d-image.png

tecplot应该可以吧？
先提取出来流线，然后对流线积分

cfd151b6-cb6a-4521-a36a-a3fb27fe3b41-image.png

a4a297dc-0192-46eb-8194-342bf20663cd-image.png

Info << patchCentre.size() << endl; // 若注释掉上一条Info语句，在这里输出patchCentre的size得到的值变为了0

换成这个呢？

Info << patch.faceCentres();

可以在 OpenFOAM cavity 算例里执行 blockMesh 后，查看 constant/plymMesh 文件夹。就是这些东西。
不过记录形式向会有区别。但都是典型的“非结构网格”。

网格系统，自己开发的通用性不强，处理起来全都得自己手写一套。当然好处就是不用去学成熟的网格系统。便于上手，但越往后，想要的功能越多，越麻烦。

如果是 OpenFOAM 的用的网格系统，paraview 就可以看，也可以用 blender 结合 blockMesh 生成。snapHexMesh 也可以。对网格、流场的各种处理也可以调用 VTK 库直接在程序里进行。

看了一下你上传的数据，没有像样的文件头，估计是刚接手一套祖传代码:136:

@dxl 噢，那可以在相应的位置加一些测试语句试试。

会不会是质量有问题？sum变成了0.5就不对了吧

更新一下
发现了问题的源头，是网格的问题。
我使用ICEM画网格，或者OpenFOAM的blockMesh工具画网格，都是先画背景网格，然后对近壁面进行多次加密。加密是八叉树类型的加密，所以在网格尺寸变化的地方，大网格的尺寸直接是小网格的二倍。这样的网格确实不好，可能一些拓扑结构上的不好导致了singleGraph程序进行插值时找不到插值的基点。
解决办法：
举例说距离壁面0~1的网格尺寸是1,1~2的网格尺寸是2，那么就使用singleGraph在0~1这个范围内提取一条线，在1~2再提取一条线，这时的type可以是uniform，也可以是midPoint或者midPointAndFace。一句话，就是分段提取信息。
我是重新画了网格， 老师说之前画的这网格根本不能用:mihu: