球装问题 Random Close Packing

做DEM模拟的第一步,寻找合适的初始构型。于是找到随机球装问题(RCP, random close packing),这个问题本身就很有意思。简单而言就是如何在有限空间内存放更多小球。估计是一个几千年的问题了,不过用计算机模拟,可以发现很多的问题,例如很容易在局部产生拥塞结构(jammed structure)。
下面是一个满足周期性边界条件的1000个等尺寸球装构型,相对密度62.8%。更高密度的结构,需要指数增长的CPU时间;但同时又可能发生局部规则结构的形成(结晶),而失去随机结构的意义。
可视化使用Mathematica,数学软件,做可视化就慢多了。

Number of Particles: 1000
Contraction Rate: 0.10000000
Radius: 0.05312100
Density: 0.62789550
File: P_N1000_0.dat (14.03.2007)

发布在科学工作. 将该链接存入书签发表评论或留个互链:互链地址.

2 条评论