请教DPMFoam和MPPICFoam中颗粒大小和网格大小之间的关系

[bubbleEL]

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

[李东岳]

在E-L方法里颗粒大小是必须在三个维度上小于网格大小么？还是可以只在至少一个维度上小于网格尺度？

三个维度

壁面附近的边界层网格在y+方向上尺度是很小的，如果颗粒的半径要超过边界层网格y+方向的尺度，那E-L方法还能算吗？

你要模拟的什么工况？如果颗粒太大会出现连续相体积分数为0的真空状态，我觉得不能算。你可以直接放弃y+的严格限定。

[bubbleEL]

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

[李东岳]

之前我也看过气泡减阻的文章，确实会出现这个问题，你模拟的是哪一篇数据？

[bubbleEL]

@东岳
是这一篇文献，没找到比较合适的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.

[李东岳]

我主要参考的是你的文献中的参考文献15Numerical simulation of micro-bubble drag reduction using population balance model，模型类似。双欧拉这面假定没有那么精细，我看到很多文献双欧拉网格比粒子小，双欧拉的相分数不是从直径计算出来的。E-L双向耦合中相分数的计算严格从粒子直径算过来的：$\alpha =k_v{d}^3/\mathrm{\Delta }V$，粒子直径太大会导致粒子的$\alpha =1$，流体的$\alpha =0$。这个倒不是主要问题。关键是流体和粒子的力交互不好计算。

我建议你多找找欧拉拉格朗日模拟气泡减阻的SCI，看看他们怎么评价这方面的问题？希望得到反馈！你也可以写道你的文章里面

[金石为开]

@bubbleEL 你好啊 我也在做气泡减阻的 想找你咱们交流下啊 你qq是啥啊 我的1277909926

[李东岳]

@bubbleEL 有什么进展么

[bubbleEL]

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

[李东岳]

@bubbleEL 好犀利，静候大作，学习学习

[李东岳]

https://www.sciencedirect.com/science/article/pii/S0021999112007462

一个网格可以小于颗粒直径的算法

[Ant]

@东岳 李老师，您好，用MPPIC计算大尺度的流化床，是不是也需要（最大）颗粒/parcel直径小于（最小）网格体积？如果计算颗粒尺度很大，是不是只能增大（最小）网格尺寸和减小parcel中的颗粒数，也就是增加parcel的数目？目前在测算算例，parcel包含颗粒较大，会导致发散，如果一直增加parcel数也不太现实，一直也找不到合适解决办法。同时有点不理解为什么需要三个维度上小于网格尺度？

[星星星星晴]

http://www.cfd-china.com/topic/1152/多相流模拟的数学模型/2

[Ant]

@星星星星晴 多谢

[Ant]

@bubbleEL 您好，请问您德高斯分布函数已经实现了吗？最近也在实现这种办法，不知道能否指点一下？另外一直在思考超出网格之后发散的数学原理是什么？不知道您这边是不是已经弄清楚了？