完全解析模拟的数值公式提出了湍流中自由运动的刚性颗粒。这项工作基于虚拟域的快速计算方法Sharma&Patankar([1]参考J.Compt。物理。,(205), 2005). 该方法避免了显式计算分布式拉格朗日乘子用于施加刚体运动并减少粒子相导致的计算开销。实施基于有限体积的能量共存的数值算法在结构化笛卡尔网格上保留分数步格式是提出了。数值方法首先验证了不同雷诺数下的固定球和自由流动重力作用下的落球。网格和时间步长收敛研究用于评估方法的准确性。最后,模拟衰减各向同性湍流中125立方粒子的研究湍流模拟的可行性自由移动、任意形状的刚性粒子的存在。