请教DPMFoam中粒子与壁面碰撞的问题

bubbleEL

最近做了一个气泡上升与上面平板发生碰撞的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函数里是怎么得到变形量的？为什么这么小？

李东岳

气泡上升与上面平板发生碰撞的case

从物理层面考虑，气泡和壁面会发生碰撞回弹么？还是贴附？

关于变形量，首先确定下：你的网格是不是足够大了？

bubbleEL

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

bubbleEL

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

李东岳

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

你这是在拉格朗日时间步的输出么？还是欧拉时间步

你增加网格尺寸试试

bubbleEL

@东岳感谢李老师的意见，我增加网格尺寸试一下。
我是在拉格朗日时间步输出的，每个流体时间步力输出了很多次，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是什么意思，什么作用，是如何得到的，李老师可否指导一下，非常感谢！

李东岳

这个site是什么意思

目前我也不清楚知道了可以告诉我一下

CFD中文网

请教DPMFoam中粒子与壁面碰撞的问题