这个问题解决了，主要是缺少了一个参数，导致和别的文件参数的声明不匹配

addBoundaryDiag(boundaryDiagCmpt, cmpt);

请参考：

void fvMatrix<Type>::addCmptAvBoundaryDiag(scalarField& diag) const { forAll(internalCoeffs_, patchI) //对每一个边界，都处理这个边界的internalCoeffs_ { addToInternalField//假如针对某个零法向梯度速度边界 ( lduAddr().patchAddr(patchI), // 这个边界的网格单元的面对应的内部网格单元编号 cmptAv(internalCoeffs_[patchI]), // 对这个patch上的internalCoeffs_，进行三个方向的平均，形成一个。 //比如原来internalCoeffs_可能是(1,2,0), (1,2,0), (1,2,0)，做完了之后成为了(1.5), (1.5), (1.5) diag//原本的对角线元素的贡献。加上考虑内部网格单元编号，与cmptAv的操作，形成一个新的internalCoeffs_的贡献 //也就是原本对角线元素 + 平均后internalCoeffs_的贡献 ); } }

boundaryDiagCmpt.negate();

H操作符是把矩阵系数移到右边。下面的上下对角线移到右边都取负。边界导致的对角线移到右边，也要取负。

for (register label face=0; face<l.size(); face++) { Hphi[u[face]] -= Lower[face]*sf[l[face]];//取负 Hphi[l[face]] -= Upper[face]*sf[u[face]];//取负 }

@Tens 都试了不行，减少并行块数量有点用，但是偶尔也犯病，这玩意我现在看玄学了

谢谢您！问题已解决，路径问题,,文件放错位置了

@李东岳 老师您好，用户手册编译库那里看了，目前任务不是编译库，只是植入求解器，问题已解决，感谢老师!

你把流速设置为0，看看还飘不。

不飘就说明不是计算出错了，是A、B都太轻了。多小的流速，都直接带走

@李东岳 老师，你好，就是最顶上的那个图片，我想在一定的时间段内以um的速度流入，一定的时间段内速度为零，达到脉冲射流的效果。

:146: :146: :146:

@星星星星晴 老师您好！我在将别人写好的算例导入我自己的OpenFOAM,编译库的时候出现了问题，请教一下您，库是如何编译的呀？链接是我编译过程中遇到的问题，感谢您的回复！Openfoam植入算例，编译库报错

@lyc 好嘞，谢谢哈

已解决，是因为FENE-P模型的数值格式不稳定导致的。我之前将槽道宽为1m算例的入口速度设置为1000m/s能够收敛是我将壁面条件设置为了slip滑移，所以可以收敛。当改为无滑移边界时，就发现不收敛了。本质上就是该数值格式可能在高雷诺数下是不稳定的。

@李东岳 现在明白啦

@李东岳 打扰李老师了，我现在有RANS的KEpsilon的DES模型，不过这个模型还是会将非选定的壁面进行KEpsilon处理，李老师您的思路是通过KEpsilon的DES模型里面通过划分区域来调用两个不同的湍流模型嘛？我在想能不能通过改变壁面边界条件让KEpsilon的DES模型识别不到选定的壁面，从而使得选定的壁面用LES模型进行计算。

@李东岳 谢谢东岳老师

@学流体的小明 好的，谢谢您

@李东岳 这回好使了，谢谢李老师 ：）

@gxf 您好，请问您计算孤立液滴蒸发的case使用的是哪个求解器呢？是否进行了重新编译？我之前试算过，结果误差很大。
image.png

我想问一下这个overInterDyMFoam求解器中移动网格的边界条件可以设置为其他的吗（不是墙壁），就是其中有一个边界在往外部喷气体，气体也可以进入背景网格当中，形成气液接触面。移动网格受到反作用力的效果，它与背景网格产生相对运动。

看起来是除以0了。 可能komega里面你设的初始条件有点问题。

这个植入在了湍流模型里面。这个turbulence->divDevRhoReff(rho, U)