你可以改一下。不太费事。

这个是全系列的虚拟机么？

这个就是高斯定理。没啥别的了。

边界面上的

也就是说这里的tSU.ref()遍历了PtrList中的所有粒子团的指针然后把每一个粒子团所属类的parcelCloud.SU函数的返回值相加了得到了一个源项的场，可以这么理解吗？

你的理解是对的

http://dyfluid.com/EL.html

@李东岳 在 icoFoam中压力泊松方程的求解 中说：

@AppleKiller http://dyfluid.com/piso.html

谢谢李老师

33ae34d9-4a5a-4ebd-8603-c18a65d1b2a3-image.png 按照经验来看，在高约化速度下，圆柱振幅应该降低的，结果振幅一直下不去，有没有大佬遇见过类似的问题？用的是圆形计算域

已知某点坐标，获取该点的速度。单核可行，多核并行就会出现由于分块找不到cellID的情况。代码如下：
point b (1,1,1);
label cellID = mesh.findCell(b);
vector intU = U[cellID];
怎么解决呐？可以用reduce函数？但是没看明白。或者有没有其它办法获得该点的速度值，使其在并行的情况下有效。
鞠躬感谢！

不知道楼主还看不看得到，我是想把翼型3维主动变形和overset结合起来。但现在第一步，用codefixedvalue来定义displacement边界，就出现了不收敛问题，不知道能否分享下您当时这部分的codes，我学习一下。

不知道楼主还看不看得到，我是想把翼型3维主动变形和overset结合起来。但现在第一步，用codefixedvalue来定义displacement边界，就出现了不收敛问题，不知道能否分享下您当时这部分的codes，我学习一下。

感谢李老师的回答！我刚接触of不久，请问contour是什么边界条件吗，我找了好久没有找到，of的算例里面有类似的吗

感谢大佬

想请教一下大家openfoam中有涉及长波辐射模拟吗，如何实现呢？

@wangfei9088 不好意思，才看见消息。
非常感谢您的回复，或许我还是应该使用FLUENT去实习自己的需求。

@青山 好的谢谢！

我算358

给大家分享一下当采用自定义的求解器对波浪进行数值仿真的时候采用wave2Foam waveGaugeProbes对波高进行监测
当前对波高进行监测主要有wave2Foam中的waveGaugeProbes，IHFoam和olaFlow都有自己的波高监测仪，OpenFOAM中也有interfaceHeight函数对波高进行监测。
但是有一个问题是，当采用interfaceHeight对alpha.water进行监测的时候，会导致OpenFOAM运行速度慢，而采用自定义的求解器对波高进行监测的时候，那只能采用sampleDict对alpha.water进行监测，这样的结果是在postProcessing会生成一系列的时间文件，需要采用其他方法对这一些列的时间文件中的波高重新提取，也是很不方便。
因此我这边就是结合wave2Foam中的waveGaugeProbes和interFoam，这样在interFOAM求解器运行的时候，可以采用wave2Foam中的波高仪进行监测，这样的效率会得到很大的提高。

wave2foam waveGuages.png
1.把interFoam求解器拷贝到wave2Foam路径的waves2Foam/applications/solvers/solvers2206_PLUS（当然这个求解器solvers2206_PLUS对应的是OpenFOAM相应的版本）
2.在interFoam文件夹中修改Make/files中修改EXE = $(FOAM_USER_APPBIN)/DinterFoam，使得编译后的程序在platforms中
3.在这个拷贝的interFoam文件夹中的Make/Options中添加如上图所示2和3步的关于wave2foam的代码
4.对wave2Foam进行编译就行
使用方法
1.在constant文件中添加probeDefinitions文件，这个文件是在wave2Foam中就有的，主要是用来设置波高仪
2.采用wave2foam的命令waveGaugesNProbes，就可以生成waveGaugesNProbes文件
3.在controlDict中添加waveGaugesNProbes文件中的 自定义波高仪名称_controlDict
4.OKOK
完美实现
地下这个是放在github的文件，和上面也是一模一样
https://github.com/byChen47/waveHeightGauge

@tidedrinker 你好，我发现of9中可以在force文件中单独定义相名称来提取相分数和对应的密度，输出的力和力矩如下：

forces forces write: sum of forces: pressure : (-10.813 0 -5.80003) viscous : (-0.000457782 -0.000707157 0.000145593) porous : (0 0 0) sum of moments: pressure : (-0.00441956 4.53963 -0.00703327) viscous : (-3.21179e-05 7.86609e-05 -2.42533e-05) porous : (0 0 0)

我的搅拌轴是y轴，所以我的扭矩是否就是4.53963N·M，但是这和实验值相差甚远，实验值大概只有0.06N·M，这是为什么？是不是还需要对这个力矩进行下一步处理才可以得到扭矩？需要除以密度嘛？