blastFoam这个求解器里面的算法还挺有意思的。求解器的作者也很厉害。我一直在关注。

@李东岳 效果非常好，非常感谢李老师:140:

自问自答，我粗略的整出来了，如图1，要注意画红圈地方的选择:chitang:
746340cd-5bc6-47e8-a92f-b507dd49d9a8-image.png
图1

@wangfei9088 你好，哈哈昨天晚上搞得太晚神智有点不清了，刚改了编译通过了，多谢

@杨英狄 就是说想以这样的方式注入拉格朗日粒子：每个时间步注入流场的粒子的数量、位置都与当前时间步的流场有关，这样的注入方式怎么实现呢？

@yfclark 你好，你这个功能实现了吗，可以请教一下如何实现的吗？

大家好，目前我可以输出温度场的三维数据，比如
QQ截图20221222151315.png
我想垂直于z轴做100个等间距垂直切面，在每个切面上对温度求积分并取平均值。
\begin{equation}
\int_{A} T dx dy /A
\end{equation}
然后沿着z轴画出一条温度平均值曲线。

谢谢大家！

@李东岳 李老师您好，D()是一个volVectorField。主要是发现建立两个fvMatrix时，前面一个fvMatrix会对后面的fvMatrix产生影响，主要体现在源项上。按理说这是两个独立的fvMatrix对象，不应该互相影响啊

你先备份一下。

然后你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: