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



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


  • 管理员

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

    三个维度

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

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



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


  • 管理员

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



  • @东岳
    是这一篇文献,没找到比较合适的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_vd^3/\Delta V$,粒子直径太大会导致粒子的$\alpha=1$,流体的$\alpha=0$。这个倒不是主要问题。关键是流体和粒子的力交互不好计算。

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



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


  • 管理员

    @bubbleEL 有什么进展么



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


  • 管理员

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


  • 管理员

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

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



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







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


Log in to reply
 


CFD中文网 | 东岳流体学术 | 东岳流体商业 | 吉ICP备20003622号-1