双流体方法颗粒堆积模拟

[youminghao]

各位大佬好，本人用双流体方法尝试了一个颗粒堆积情况，如下图所示，选定一个区域初始填充为颗粒床，然后在重力作用下下落，最后颗粒堆积没有形成实际颗粒堆积的堆积角，请问一下是我设置有问题还是双流体方法不能够实现？

[李东岳]

我在想控制方程里面有哪一项可以导致这种情况发生。 @xpqiu 大佬怎么看？

DEM那面是哪一项实现的堆积角？摩擦应力比较大导致粒子不变？

[youminghao]

@李东岳 李老师，您好，DEM中实现颗粒堆积角的应该是自身的内摩擦角吧，您说的摩擦应力大导致粒子不变是什么意思？密集堆积的话颗粒动理论的作用基本可以忽略，主要还是摩擦应力模型计算占主导，那是不是关于颗粒粘度的计算会影响颗粒堆积？

[xpqiu]

@李东岳
对，摩擦应力(frictional stress)很重要。加上摩擦应力才能在双流体模型中模拟出来堆积角。

[youminghao]

@xpqiu 老师，您好，我是采用了JohnsonJacksonSchaeffer模型来计算摩擦应力，但是最终没有发现堆积角，您认为可能是哪方面的问题？是应力模型本身的问题还是模型参数的影响？

[xpqiu]

@youminghao

可能是J-J 摩擦应力模型的问题。

上面图里面第一个是不考虑摩擦应力的结果，下面两个是考虑摩擦应力的结果。摩擦应力模型用的是 Laux，对应参考文献是：

A. Nikolopoulos, N. Nikolopoulos, N. Varveris, S. Karellas, P. Grammelis, and E.
Kakaras, Investigation of proper modeling of very dense granular flows in the
recirculation system of CFBs. Particuology, 2012, 10(6):699-709

[youminghao]

@xpqiu 好的，谢谢老师，我学习一下