Full Text:   <2180>

Summary:  <13>

Suppl. Mater.: 

CLC number: 

On-line Access: 2026-03-25

Received: 2025-05-21

Revision Accepted: 2025-09-09

Crosschecked: 2026-03-25

Cited: 0

Clicked: 1988

Citations:  Bibtex RefMan EndNote GB/T7714

 ORCID:

Daosheng LING

https://orcid.org/0000-0002-0604-1175

-   Go to

Article info.
Open peer comments

Journal of Zhejiang University SCIENCE A 2026 Vol.27 No.3 P.231-245

http://doi.org/10.1631/jzus.A2500203


Monte Carlo simulation method for estimating the fine rattler fraction in large-ratio binary mixtures


Author(s):  Jing WANG, Changyu SHI, Daosheng LING

Affiliation(s):  1. MOE Key Laboratory of Soft Soils and Geoenvironmental Engineering, Institute of Geotechnical Engineering, Zhejiang University, Hangzhou 310058, China more

Corresponding email(s):   dsling@zju.edu.cn

Key Words:  Granular materials, Rattlers, Fines content, Size ratio, Discrete element method (DEM), Monte Carlo


Jing WANG, Changyu SHI, Daosheng LING. Monte Carlo simulation method for estimating the fine rattler fraction in large-ratio binary mixtures[J]. Journal of Zhejiang University Science A, 2026, 27(3): 231-245.

@article{title="Monte Carlo simulation method for estimating the fine rattler fraction in large-ratio binary mixtures",
author="Jing WANG, Changyu SHI, Daosheng LING",
journal="Journal of Zhejiang University Science A",
volume="27",
number="3",
pages="231-245",
year="2026",
publisher="Zhejiang University Press & Springer",
doi="10.1631/jzus.A2500203"
}

%0 Journal Article
%T Monte Carlo simulation method for estimating the fine rattler fraction in large-ratio binary mixtures
%A Jing WANG
%A Changyu SHI
%A Daosheng LING
%J Journal of Zhejiang University SCIENCE A
%V 27
%N 3
%P 231-245
%@ 1673-565X
%D 2026
%I Zhejiang University Press & Springer
%DOI 10.1631/jzus.A2500203

TY - JOUR
T1 - Monte Carlo simulation method for estimating the fine rattler fraction in large-ratio binary mixtures
A1 - Jing WANG
A1 - Changyu SHI
A1 - Daosheng LING
J0 - Journal of Zhejiang University Science A
VL - 27
IS - 3
SP - 231
EP - 245
%@ 1673-565X
Y1 - 2026
PB - Zhejiang University Press & Springer
ER -
DOI - 10.1631/jzus.A2500203


Abstract: 
Reliable estimation of the fine particle rattler fraction is crucial for understanding the structural and mechanical responses of binary granular systems with large size ratios. However, such estimation is challenged by the general inability to obtain interparticle contact force information directly from experimental images and by the lower accuracy of positional and size identification of fine particles compared with coarse particles. To address these challenges, in this study, we focused on 2D bidisperse granular assemblies with large size ratios (α=7, 9, 12, and 16) and developed an approach based on monte Carlo simulation (MCS) that relies solely on the size and positional information of coarse particles, avoiding the need for force-resolved computations. The performance of the method was evaluated against experimental measurements and discrete element method (DEM) simulations. The MCS-based predictions show close agreement with experimental results, with a slight overall overestimation. At low fines content, the approach tends to overestimate the fine particle rattler fraction relative to DEM results, whereas at higher fines content, it underestimates the rattler fraction. Overall, the proposed MCS-based approach enables robust and relatively accurate estimation of the fine particle rattler fraction. This study provides a practical framework for predicting the rattler fraction, contributes to advancing both experimental analysis and theoretical modeling in granular physics, and demonstrates the conceptual extendibility of the MCS framework to more complex 3D packings.

高粒径比二元颗粒堆积体系中悬浮细颗粒分数的蒙特卡洛统计分析

作者:王景1,2,施昌宇1,凌道盛1,2
机构:1浙江大学,软弱土与环境土工教育部重点实验室,中国杭州,310058;2浙江大学,超重力研究中心,中国杭州,310058
目的:在高粒径比二元颗粒体系中,准确统计悬浮细颗粒分数对理解体系的结构与力学响应至关重要。然而,试验中通常难以直接获得颗粒间接触力信息,同时细颗粒的位置与尺寸识别精度有限。本文提出了一种基于蒙特卡洛模拟(MCS)的统计方法,该方法仅依赖于粗颗粒的尺寸与空间位置信息,通过几何分析实现对悬浮细颗粒分数的有效估算,为相关的统计分析提供一种高效可行的替代途径。
创新点:1.克服离散元模拟需要进行动力学循环迭代所带来的计算量偏大等局限;2.实现以小规模样本的重复随机抽样来近似大样本的统计特性。
方法:1.利用Delaunay三角剖分提取粗颗粒体系的孔隙分布,并建立细颗粒填充的MCS框架(图3);2.在多次蒙特卡洛采样中,根据孔隙容量判定细颗粒是否为悬浮颗粒,并通过统计得到平均悬浮细颗粒比例(图6~8);3.将MCS预测结果与试验观测及离散元方法(DEM)模拟进行对比,评估方法的精度与适用性(图15)。
结论:1.所提出的MCS方法能够在无需动力学求解条件下可对悬浮细颗粒分数进行稳健、高效且较为准确的统计;2.研究为颗粒材料的试验分析和理论建模提供了实用框架,并展示了方法向三维复杂体系扩展的潜力。

关键词:颗粒材料;悬浮颗粒;细颗粒含量;粒径比;离散元方法(DEM);蒙特卡洛

Darkslateblue:Affiliate; Royal Blue:Author; Turquoise:Article

Reference

[1]Al-RaoushR, ThompsonK, WillsonCS, 2003. Comparison of network generation techniques for unconsolidated porous media. Soil Science Society of America Journal, 67(6):1687-1700.

[2]BabuV, VinuthaHA, BiDP, et al., 2023. Discontinuous rigidity transition associated with shear jamming in granular simulations. Soft Matter, 19(48):9399-9404.

[3]BagiK, 2005. An algorithm to generate random dense arrangements for discrete element simulations of granular assemblies. Granular Matter, 7(1):31-43.

[4]BehringerRP, ChakrabortyB, 2019. The physics of jamming for granular materials: a review. Reports on Progress in Physics, 82(1):012601.

[5]BerziD, VescoviD, 2021. Cooling after shearing: three possible fates for dense granular materials. Granular Matter, 23(2):47.

[6]BiazzoI, CaltagironeF, ParisiG, et al., 2009. Theory of amorphous packings of binary mixtures of hard spheres. Physical Review Letters, 102(19):195701.

[7]BryantSL, KingPR, MellorDW, 1993. Network model evaluation of permeability and spatial correlation in a real random sphere packing. Transport in Porous Media, 11(1):53-70.

[8]ChengZ, WangJF, 2018. Experimental investigation of inter-particle contact evolution of sheared granular materials using X-ray micro-tomography. Soils and Foundations, 58(6):1492-1510.

[9]CGAL (The Computational Geometry Algorithms Library), 2024. CGAL 6.1–Manual. https://doc.cgal.org/latest/Manual/packages.html

[10]CundallPA, StrackODL, 1979. A discrete numerical model for granular assemblies. Géotechnique, 29(1):47-65.

[11]DuvergerS, DuriezJ, PhilippeP, et al., 2021. Rattlers’ involvement for possibly looser critical states under higher mean stress. EPJ Web of Conferences, 249:11002.

[12]ElishakoffI, 2003. Notes on philosophy of the Monte Carlo method. International Applied Mechanics, 39(7):753-762.

[13]FengYT, HanK, OwenDRJ, 2003. Filling domains with disks: an advancing front approach. International Journal for Numerical Methods in Engineering, 56(5):699-713.

[14]GaoGH, LiDQ, CaoZJ, 2023. Identification of the design point based on Monte Carlo simulation. Computers and Geotechnics, 159:105438.

[15]HuTT, LingDS, TangY, et al., 2024. Evolution of rattlers under different confining pressures and its influence on shear strength. Transportation Geotechnics, 44:101176.

[16]HuY, WangY, 2020. Probabilistic soil classification and stratification in a vertical cross-section from limited cone penetration tests using random field and Monte Carlo simulation. Computers and Geotechnics, 124:103634.

[17]HuangX, HanleyKJ, O’SullivanC, et al., 2014. Exploring the influence of interparticle friction on critical state behaviour using DEM. International Journal for Numerical and Analytical Methods in Geomechanics, 38(12):1276-1297.

[18]HurleyRC, LimKW, RavichandranG, et al., 2016a. Dynamic inter-particle force inference in granular materials: method and application. Experimental Mechanics, 56(2):217-229.

[19]HurleyRC, HallSA, AndradeJE, et al., 2016b. Quantifying interparticle forces and heterogeneity in 3D granular materials. Physical Review Letters, 117(9):098005.

[20]IndraratnaB, HaqS, RujikiatkamjornC, et al., 2022. Microscale boundaries of internally stable and unstable soils. Acta Geotechnica, 17(5):2037-2046.

[21]Itasca Consulting Group, Inc., 2014. PFC—Particle Flow Code, Ver. 5.0. Itasca, Minneapolis, USA.

[22]JiangHR, KawamotoR, MatsushimaT, 2024. Exploring the combined effect of particle shape and friction on cell structures in granular packings via LS-DEM modeling. IOP Conference Series: Earth and Environmental Science, 1330(1):012047.

[23]KumarN, MagnanimoV, RamaioliM, et al., 2016. Tuning the bulk properties of bidisperse granular mixtures by small amount of fines. Powder Technology, 293:94-112.

[24]KwanAKH, ChanKW, WongV, 2013. A 3-parameter particle packing model incorporating the wedging effect. Powder Technology, 237:172-179.

[25]LadeP, LiggioC, YamamuroJ, 1998. Effects of non-plastic fines on minimum and maximum void ratios of sand. Geotechnical Testing Journal, 21(4):336-347.

[26]LeeCL, BililignE, AzémaE, et al., 2025. Loading-dependent microscale measures control bulk properties in granular material: an experimental test of the stress-force-fabric relation. Physical Review E, 112(3):035401.

[27]LeeK, HurleyRC, 2025. Interplay between forces, particle rearrangements, and macroscopic stress fluctuations in sheared two-dimensional granular media. Physical Review E, 112(1):015413.

[28]LiJ, WangD, 2023. Comparison of PDEM and MCS: accuracy and efficiency. Probabilistic Engineering Mechanics, 71:103382.

[29]LiRD, YinZY, HeSH, 2025. 3D reconstruction of arbitrary granular media utilizing vision foundation model. Applied Soft Computing, 169:112599.

[30]LinMC, WangGQ, ZhouJ, et al., 2023. Discrete element method study of hysteretic behavior and deformation characteristics of rockfill material under cyclic loading. Journal of Zhejiang University-SCIENCE A, 24(4):350-365.

[31]LingDS, HuTT, WangJ, et al., 2024. Numerical study on the shear strength of granular materials under the low confining pressure. Soils and Foundations, 64(3):101447.

[32]MaherCE, TorquatoS, 2024. Hyperuniformity scaling of maximally random jammed packings of two-dimensional binary disks. Physical Review E, 110(6):064605.

[33]MajmudarTS, SperlM, LudingS, et al., 2007. Jamming transition in granular systems. Physical Review Letters, 98(5):058001.

[34]MeshkinghalamH, Emami TabriziM, ChenaghlouMR, 2025. Influence of the contact model on energy fluctuations in non-cohesive granular materials subjected to confinement axial cyclic loading using DEM. Granular Matter, 27(1):15.

[35]NavarroHA, de Souza BraunMP, 2013. Determination of the normal spring stiffness coefficient in the linear spring-dashpot contact model of discrete element method. Powder Technology, 246:707-722.

[36]NewhallKA, 2025. Sampling the large-dimensional energy landscape of a 2D granular system with the hydra string method. Soft Matter, 21(24):4839-4848.

[37]NguyenNS, TahaH, MarotD, 2021. A new Delaunay triangulation-based approach to characterize the pore network in granular materials. Acta Geotechnica, 16(7):2111-2129.

[38]Otto von Guericke University Magdeburg, 2025a. The Best Known Packings of Equal Circles in a Circle. http://hydra.nat.uni-magdeburg.de/packing/cci

[39]Otto von Guericke University Magdeburg, 2025b. The Best Known Packings of Equal Spheres in a Sphere. http://hydra.nat.uni-magdeburg.de/packing/ssp/ssp.html

[40]PetitJC, SperlM, 2025. Additional jamming transition in two-dimensional bidisperse granular packings. Physical Review Research, 7(2):L022073.

[41]PrasadI, SantangeloC, GrasonG, 2017. Subjamming transition in binary sphere mixtures. Physical Review E, 96(5):052905.

[42]ReboulN, VincensE, CambouB, 2008. A statistical analysis of void size distribution in a simulated narrowly graded packing of spheres. Granular Matter, 10(6):457-468.

[43]ReboulN, VincensE, CambouB, 2010. A computational procedure to assess the distribution of constriction sizes for an assembly of spheres. Computers and Geotechnics, 37(1-2):195-206.

[44]RedaelliI, di PriscoC, 2019. Three dimensional steady-state locus for dry monodisperse granular materials: DEM numerical results and theoretical modelling. International Journal for Numerical and Analytical Methods in Geomechanics, 43(16):2525-2550.

[45]ReyesSC, IglesiaE, 1991. Monte Carlo simulations of structural properties of packed beds. Chemical Engineering Science, 46(4):1089-1099.

[46]RobleeLHS, BairdRM, TierneyJW, 1958. Radial porosity variations in packed beds. AIChE Journal, 4(4):460-464.

[47]SantosAP, BolintineanuDS, GrestGS, et al., 2020. Granular packings with sliding, rolling, and twisting friction. Physical Review E, 102(3):032903.

[48]ShireT, O’SullivanC, 2016. Constriction size distributions of granular filters: a numerical study. Géotechnique, 66(10):826-839.

[49]ShireT, O’SullivanC, BarretoD, et al., 2013. Quantifying stress-induced anisotropy using inter-void constrictions. Géotechnique, 63(1):85-91.

[50]ShireT, O’SullivanC, HanleyKJ, et al., 2014. Fabric and effective stress distribution in internally unstable soils. Journal of Geotechnical and Geoenvironmental Engineering, 140(12):04014072.

[51]SrivastavaI, RobertsSA, ClemmerJT, et al., 2021. Jamming of bidisperse frictional spheres. Physical Review Research, 3(3):L032042.

[52]SufianA, RussellAR, WhittleAJ, et al., 2015. Pore shapes, volume distribution and orientations in monodisperse granular assemblies. Granular Matter, 17(6):727-742.

[53]TejadaIG, ArroyoM, CiantiaM, 2025. Stress field on intermediate length scales during triaxial compression. Géotechnique, 75(6):716-731.

[54]ThakurMM, PenumaduD, 2021. Influence of friction and particle morphology on triaxial shearing of granular materials. Journal of Geotechnical and Geoenvironmental Engineering, 147(11):04021118.

[55]ThorntonC, 2000. Numerical simulations of deviatoric shear deformation of granular media. Géotechnique, 50(1):43-53.

[56]TingJM, KhwajaM, MeachumLR, et al., 1993. An ellipse-based discrete element model for granular materials. International Journal for Numerical and Analytical Methods in Geomechanics, 17(9):603-623.

[57]WangJ, HuTT, ShiCY, et al., 2026. Partial coordination number and contact proportion of densely packed binary random mixtures of quasi-two-dimensional granular materials. Powder Technology, 467:121533.

[58]WangY, 2011. Reliability-based design of spread foundations by Monte Carlo simulations. Géotechnique, 61(8):677-685.

[59]WangY, CaoZJ, AuSK, 2010. Efficient Monte Carlo simulation of parameter sensitivity in probabilistic slope stability analysis. Computers and Geotechnics, 37(7-8):1015-1022.

[60]WautierA, BonelliS, NicotF, 2019. Rattlers’ contribution to granular plasticity and mechanical stability. International Journal of Plasticity, 112:172-193.

[61]XiaCJ, LiJD, KouBQ, et al., 2017. Origin of noncubic scaling law in disordered granular packing. Physical Review Letters, 118(23):238002.

[62]XuMQ, GuoN, YangZX, 2021. Particle shape effects on the shear behaviors of granular assemblies: irregularity and elongation. Granular Matter, 23(2):25.

[63]YangDZ, ChuXH, 2024. Influence of particle shape on creep and stress relaxation behaviors of granular materials based on DEM. Computers and Geotechnics, 166:105941.

[64]YangS, LacasseS, SandvenR, 2006. Determination of the transitional fines content of mixtures of sand and non-plastic fines. Geotechnical Testing Journal, 29(2):102-107.

[65]YoshiiK, OtsukiM, 2025. Mechanical and geometrical properties of jammed weakly cohesive granular materials. Journal of the Physical Society of Japan, 94(8):084801.

[66]ZhangJ, MajmudarTS, TordesillasA, et al., 2010. Statistical properties of a 2D granular material subjected to cyclic shear. Granular Matter, 12(2):159-172.

[67]ZhangL, GongWP, LiXX, et al., 2022. A comparison study between 2D and 3D slope stability analyses considering spatial soil variability. Journal of Zhejiang University-SCIENCE A, 23(3):208-224.

Open peer comments: Debate/Discuss/Question/Opinion

<1>

Please provide your name, email address and a comment





Journal of Zhejiang University-SCIENCE, 38 Zheda Road, Hangzhou 310027, China
Tel: +86-571-87952783; E-mail: cjzhang@zju.edu.cn
Copyright © 2000 - 2026 Journal of Zhejiang University-SCIENCE