辐射防护 ›› 2026, Vol. 46 ›› Issue (3): 224-232.doi: 10.27045/j.1000-8187.202603005

• 辐射防护评价 • 上一篇    下一篇

基于DEM-CFD耦合的放射性颗粒再悬浮起动机制与阈值模拟研究

郝毅杰1,2, 康晶1,2, 崔锦江1,2, 陈海龙1,2, 武翡翡1,2, 廉冰1,2   

  1. 1.中国辐射防护研究院,太原 030006;
    2.生态环境部辐射环境与健康重点实验室,太原 030006
  • 收稿日期:2026-03-03 出版日期:2026-05-20 发布日期:2026-06-18
  • 通讯作者: 廉冰。E-mail:lianbing@cirp.org.cn
  • 作者简介:郝毅杰(2001—),男,2023年毕业于清华大学工程物理系,现为中国辐射防护研究院核环境科学研究所辐射防护与环境保护专业在读硕士研究生。E-mail:haoyijie@cirp.org.cn

Microscopic simulation of resuspension initiation mechanisms and thresholds for radioactive particles on rough granular beds using DEM-CFD coupling

HAO Yijie1,2, KANG Jing1,2, CUI Jinjiang1,2, CHEN Hailong1,2, WU Feifei1,2, LIAN Bing1,2   

  1. 1. China Institute for Radiation Protection,Taiyuan 030006;
    2. Key Laboratory of Radiation Environment & Health of the Ministry of Ecology and Environment,Taiyuan 030006
  • Received:2026-03-03 Online:2026-05-20 Published:2026-06-18

摘要: 放射性污染土壤颗粒的再悬浮途径是辐射环境安全评估中环境排放途径之一。本文采用离散元法与计算流体力学(DEM-CFD)双向耦合方法对10 μm的普通沙尘(SiO2)和典型放射性颗粒(U3O8、PuO2)在粗糙颗粒床层表面的起动过程进行微观解析和研究。结果表明:密铺粗糙床层的几何形态对边界层流动有显著影响,粗糙床层相较理想光滑壁面显著抬升边界层结构,并增大近壁垂直速度分量;颗粒再悬浮起动是原位微振、滚动爬升与最终脱附的复杂时空演化过程。模拟结果显示:10 μm的普通沙尘的临界起动风速为1.60 m/s,而典型10 μm的放射性颗粒(U3O8和PuO2)分别提升至2.05 m/s和2.21 m/s,增幅约28%和38%,明显表征出由于放射性核素在密度方面的加持增强了重力与强表面粘附的协同阻碍作用。本文能够为放射性颗粒再悬浮阈值参数化及污染场地环境风险评估模型提供微观机理支撑和阈值参考。

关键词: 放射性颗粒, 再悬浮, 离散元法, 计算流体力学

Abstract: Resuspension of radioactively contaminated soil particles represents one of the key environmental release pathways and critical source terms in radiological environmental safety assessment. This study employs a bidirectional DEM-CFD coupling approach to perform microscopic analysis of the incipient motion process of 10 μm ordinary dust particles (SiO2) and typical radioactive particles (U3O8 and PuO2) on rough particle bed surfaces. The results reveal that the geometric morphology of densely packed rough beds exerts a significant influence on boundary layer flow; compared to an ideal smooth wall, the rough bed markedly displaces the boundary layer structure upward and increases the near-wall vertical velocity component. Particle resuspension initiation emerges as a complex spatiotemporal evolutionary process involving in-situ micro-vibration, rolling ascent, and eventual detachment. Simulation results indicate that the critical threshold wind speed for ordinary dust is 1.60 m/s, whereas for the typical radioactive particles U3O8 and PuO2 it rises to 2.05 m/s and 2.21 m/s, respectively—an increase of approximately 28% and 38% respectively. This clearly demonstrates the synergistic inhibitory effect of enhanced gravity and stronger surface adhesion resulting from the higher density imparted by radioactive nuclides. The present work provides microscopic mechanistic insight and more accurate threshold references for parameterizing radioactive particle resuspension and for environmental risk assessment models at contaminated sites.

Key words: radioactive contamination, particle resuspension, discrete element method, computational fluid dynamics

中图分类号: 

  • X591