老铁们，一般报这种错怎么解决？
a01bec9b-dbe4-484b-a741-3611893c01d3-image.png

在本地可以多算几步，但是超算就出问题了

@cresendo 在 计算钝体涡激振动时发散（pimpleDyMFoam） 中说：

songtao.chen@sjtu.edu.cn

谢谢！非常感谢您的建议！

@李东岳 好的，谢谢李老师

谢谢李老师的建议！问题已解决，更新一下！
出现上述未定义的引用（undefined reference）错误提示，是因为我在编译自己的库specie时没有改名，编译后出现libspecie.o，可能是与官方的重复了，在我手动删除该文件后，求解器可以正常使用。
参考bestucan 版主的回复！链接文本
不过我有一个疑问，我的编译路径是FOAM_USER_LIBBIN，为什么还会导致这样的错误呢？

@zousiyu 如果能够构造出一个量，它既能表征不同相，又是连续的，本身还守恒（没有源项），你这个问题是不是就好解决了？

各位前辈好！

如图1 ，我现在流场中有一些气泡，想将每一个气泡的内部压力映射到底部壁面上，效果图如图2，我应该如何在paraview中进行映射/投影呢？我自己试验了pointDatasetInterpolator好像不太对（图3）。

哪位前辈能够指点一个思路呢？谢谢

图1：
1111111111.png

图2：

2222222.png

图3：
3333333.png

@李东岳

@李东岳 在 OF中在有限的计算域中怎么消去边界对气泡的影响 中说：

既用小的计算域同时又能消除边界的强影响，从而达到使用cyclic边界条件

算法上倒是有一些。比如那些advection边界条件。但是扩大计算域是最稳妥的方法。你这个是空炮么？可以设置非均匀网格处理。外围网格非常非常稀疏，里面足够细，也不会太多的网格。

锯齿网格可能开始会有，但随着计算，表面张力的作用会变平

李老师你好，算例的确是计算空泡的；

扩大计算域的话，可能没办法满足我的需求，原因如下：

气泡的间距是算例的一个主要变量，很需要近距离的情况 如果足够远的话，循环边界条件可能就没有作用了 非均匀网格的话，一方面无法确保在流体中的压力、速度的正确捕捉；其次在网格密-疏位置处，是不是会使得结果误差比较大呢？ 如果一开始就是锯齿网格的话，那后续的计算精度感觉不能保证。原因是在高的内压驱动下，初始泡的大小决定之后的每步计算结果，所以一开始的边界捕捉及其重要

所以还是想咨询下，能不能一开始在0文件中就能很好的捕捉相分数场，粗的背景网格+细气泡边界网格，之后就用自适应网格计算。

其次还问下advection边界条件的使用方法或相关链接

麻烦李老师了~

@李东岳 您好，李老师，请问"c0"是我所设定的cellZone的名字吗，然后"s"是所设定的源的值。我本来是想先通过toposet选定一个cellZone，然后在fvoption里面对所选cellZone进行源的赋值，但是不清楚浓度源应该选择哪种type。

十分感谢两位大佬:140:

PISO算法中的对流项fvm::div(phi,U)是如何迭代求解的呢？

@李东岳 李老师，对于周期边界下的压力场，piso循环中有pEqn.setReference(pRefCell,pRefValue);用来确定压力场。那在计算周期边界下的被动标量时，应该用什么来对温度场进行修正呢？
:135:

ff12dfb1-d78b-4f7d-812d-db36e6749698-图片.png
请问各位老师，根据参差图 --Initial residual , Final residual
这个能判定收敛吗？
初始残差没到设定的值，且计算到后边开始波动，红色方框所示；
Final residual虽然可以到设定值，但是在一部分计算中波动很大，红色圆圈所示。
请问是什么原因呢。
期待大佬回答解惑。

@五好青年 可以，就是easycae的网盘，你注册个用户应该就有，

@李东岳 哈哈哈，物理隔绝了:mihu:

图一梯度算子的分量$\frac{\p }{\p x}$，不是某个物理量的梯度分量

$\frac{\p }{\p x}$不是一个函数么，你要写一个函数？

@李东岳 收到，谢谢李老师。

液滴尺寸比网格大，蒸发源项(组分、能量)加到网格里面会不会出问题，有没有老师知道OpenFoam这方面是怎么处理这个问题的

@李东岳 谢谢李老师。

网格流动性是什么..

在WallSpringSliderDashpot.C中，有这样几句：

scalar etaN = alpha_*sqrt(p.mass()*kN)*pow025(normalOverlapMag); vector fN_PW = rHat_PW *(kN*pow(normalOverlapMag, b_) - etaN*(U_PW & rHat_PW)); ... scalar etaT = etaN; ... (if条件：颗粒不滑移的条件下) fT_PW = - kT*tangentialOverlap_PW - etaT*USlip_PW;

这里的颗粒碰撞中的切向阻尼系数，etaN，等于法向值，然后分别应用到求碰撞力的切向和法向分量，fN, fT，中。是不是意味着，我在颗粒配置文件中通过alpha设置的回弹系数值e, 就是同时设定了切向和法向的回弹系数，eT, eN, 都等于e？

请问这个理解对不对？