@散漫守望2016 感谢回复！我找出问题了，这个颗粒直径应该设置成particle diameter，即单个小颗粒的直径，而不是parcel equivalent diameter，即假设的parcel等效直径。这个发散的原因可能是直径设置错误当把nParticle 设置成100时颗粒相体积分数太大了，所以就没法算了。
解释一下供后面的新手参考

@孙智勇 格拉茨大学他们搞的那个CFD-DEM， 好像是ligggts，用LAMMPS处理颗粒跟踪的，据说可以搞不同的颗粒形状。你可以看看他们最近还有没有继续搞

昨天在一个大佬的文章里面看到这么一句 “At an even coarser scale, it is common to adopt an Euler–Euler (EE) approach consisting of separate Eulerian transport equations for the locally averaged variables of both the gas and solid phases (6, 16–18). This approach, commonly referred to as a two-fluid model (TFM), ......”
出自 Sundaresan, Sankaran, Ali Ozel, and Jari Kolehmainen. "Toward constitutive models for momentum, species, and energy transport in gas–particle flows." Annual review of chemical and biomolecular engineering 9 (2018): 61-81. https://www.annualreviews.org/doi/10.1146/annurev-chembioeng-060817-084025

@ibelief 可能还没自由到投资CFD 软件的钱等于早饭钱

@李东岳 那只能祝李博早日财富自由，就可以自主投资了！

@五好青年 sorry 回错人了……

@五好青年 哎，我看2月28号不是还在招人吗……怎么突然就黄了？

@下里巴人 看到这个帖子想起我之前做那个流化床模拟，加了跟管子就死活不收敛了，估计和你这个一样的问题。
话说有尺度差异的时候出来的这个网格好呆……明明圆形的楞给画成多边形了。
我想起Gmsh他们组最近有个博士答辩的东西就是做curve mesh的，不知道能不能用在计算上了。。

好像这货是24C 48T的……
说起来windows的内存管理我一直很不爽，因为他老是不让我用完内存……

@韬智tz 看了下你的回复，非常建议你去系统的学习一下tutorials，至少学习下OpenFOAM或者Fluent的基础操作再开始做计算

@hxy970122 扫了一眼Fluent的组分输运方程，对于各个组分就没有重力的影响，所以如果组分分离就肯定是因为扩散率不同。但是在binary mixture里面，似乎两个组分的扩散系数是相同的？

重力的影响是在质量方程和动量方程里面，所以就是所有组分受到的重力作用是平均重力……