在控制固体粒子相分数时,当粒子超过"极限",那么粒子怎么被分配呢?

[zhe]

大家好,我想问一下,在研究MPPIC这个方法时,大家应该都遇到了"固体粒子相分数极限"这个项.也就是控制网格内固体粒子的占比最大值.但是我的问题是,如果超过这个占比的粒子想要移动到这个网格里,是被阻挡了吗?是通过压力吗?那么被阻挡下来的粒子又是怎么安置的呢?是阻挡在了附近网格,还是在一个时间步长内,把原本装满粒子的网格挤出相对数量的粒子,来保证质量守恒?希望大家能多发言讨论,哪怕是猜测的,只要有道理,欢迎!!!

[李东岳]

单独从理论来看，这是通过固相应力（粒子应力）实现的，跟双流体模型类似，看Andrew and O'Rourke大牛1996年的文章，里面所有stress $\tau$的地方，固相应力是一个矢量，这个量用于修正接近极限的粒子速度

https://www.jianguoyun.com/p/DRX1P0wQ9s3ZBhiKwKoD

[李东岳]

坚果云竟然挂了？？？？

[马乔]

当然是固相应力呀！Harris & Crighton模型

[zhe]

@马乔 这个是固相应力以压力方式排斥粒子，不让它进入该网格。但是由于用的是normal stress，那么粒子返回的路径就是个问题了。我做了几个小tests，粒子都是直接被弹走，并且不能预测它去哪里。还有就是，被弹走的粒子也对原来饱和的网格内的粒子有反作用力，使得它们分散看来。不知道这个现象有没有什么办法可以规避？希望您多指导一下。

[马乔]

你要精确预测轨迹建议还是用DEM吧。

[zhe]

@马乔 您好,我并不是想要追踪精确的粒子运动轨迹,只是当网格饱和时,这个方法的处理方式确实会引起结果的不确定性.我再看看是否可以运用其他的方法,谢谢您!

[李东岳]

@zhe 你在模拟稀疏流还是稠密的？

[zhe]

@东岳 稠密的.sediment transportation.