CFD中文网

@液固两相流 我最近没有在做这方面的工作了，之前我遇到的问题在更新版本的OpenFOAM（v7、v8等）中由于官方优化得到了一定的改善，你也可以试试。

请教各位老师，在DPMFoam的文件夹中还有一个叫DPMTurbulenceModels的文件夹，里面只有一个C文件，内容如下图，在编译的时候好像是把要用的湍流模型又编译了一下，不太明白这样做的必要性，在createFields的时候直接用原本的湍流类来创建turbulence对象不可以吗？

@东岳 我在计算alpha的时候用了高斯分布函数把气泡体积光滑到周围的网格当中，分布范围就是气泡直径，这样当气泡直径没超过网格大小的时候，还是像原来那样减，如果像我之前说的，比边界层网格要高了，气泡实际占据的网格每一个网格都会分一部分气泡的体积。这样比较真实，稳定性也好

@东岳
昨晚粒子可以进入旋转域了，但是在里面的粒子还是有问题，就是某一个时间步突然一个网格速度猛增（如下图1），然后库朗数就迅速上升到几百，然后下一步就不能算了，报错是说粒子追踪的时候没法对目标面做三角分割了，我这个是在1706用老的追踪方法算的，我也试过v6和1812这种新的质心追踪计算，不过旋转域里面的粒子还是有问题。并且所有版本的OF好像都不能在粒子追踪求解器里面把AMI动域并行分块到多个处理器中，这是另一个致命打击- -

如题请教一下，我现在尝试计算的时候颗粒无法穿过AMI面，一接触就会发散。怀疑是不是patchInteractionModel给的有问题，我现在对AMI面是给的none，不过我看好像也没有什么合适的interaction形式。请教各位，感谢。

@东岳
是这一篇文献，没找到比较合适的E-L方法模拟的，这个是有试验和E-E方法模拟对比的。
Qin S , Chu N , Yao Y , et al. Stream-wise distribution of skin-friction drag reduction on a flat plate with bubble injection[J]. Physics of Fluids, 2017, 29(3):037103.

我是在模拟气泡减阻的问题，参照文献中的试验，流速为4m/s，气泡直径大约1mm，但是平板第一层网格厚度大概在0.2mm以下才能把摩擦阻力计算的准确。在这样的网格条件下，我通气量小一点还可以算，通气量过大就会导致流体计算发散。很可能是您说的真空现象导致的。思考了好几天也没有什么好方法，老师您有没有什么意见？

请教各位一个问题，在E-L方法里颗粒大小是必须在三个维度上小于网格大小么？还是可以只在至少一个维度上小于网格尺度？壁面附近的边界层网格在y+方向上尺度是很小的，如果颗粒的半径要超过边界层网格y+方向的尺度，那E-L方法还能算吗？

@东岳 感谢李老师的意见，我增加网格尺寸试一下。
我是在拉格朗日时间步输出的，每个流体时间步力输出了很多次，normaloverlapMag的值没什么差别，基本都是2e-6。
这个normalOverlap的确定是跟这个函数参数site相关的。如第一行：

  vector r_PW = p.position() - site;

这个函数在lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/WallModel/WallSpringSliderDashpot/WallSpringSliderDashpot.C中
这个site的来源好像是PairCollision.C中的wallInteraction函数
我没太明白在这个函数里这个site是什么意思，什么作用，是如何得到的，李老师可否指导一下，非常感谢！

@东岳 网格尺寸是稍大于气泡直径的，1.2倍气泡直径

@东岳 有试验结果，气泡会稍微回弹一下，然后贴附在上平板上，维持一定的变形量，在我这个工况下大概是0.2倍的半径。

最近做了一个气泡上升与上面平板发生碰撞的case，原本DPMFoam中弹簧阻尼模型中的刚度是用杨氏模量确定的，我将刚度改成了气泡碰撞模型常用的刚度。case模拟的过程如下图所示，模拟结果想要看看气泡在与上方平板发生碰撞后的变形量是否合理。



运行结束后，我查看最后时刻的粒子position，质心y方向是0.0081，上方平板的高度是0.0084，粒子半径是0.0006，那变形量就应该是0.0003，这个量相对半径很大，我想弹性力计算是不是有问题，检查代码的过程中我想在evaluateWall函数中输出变形量

vector r_PW = p.position() - site;

vector U_PW = p.U() - data.wallData();

scalar U_PW_mag = mag(U_PW);

scalar r_PW_mag = mag(r_PW);

scalar normalOverlapMag = max(pREff - r_PW_mag, 0.0);

Info << "normal overlap = " << normalOverlapMag <<nl;

输出的结果却是2e-6这样一个很小的变形量，我不太确定变形量到底应该是哪个？还有就是PairCollision.C中wallInteraction函数里是怎么得到变形量的？为什么这么小？