风沙运动的DEM数值模拟.pdfVIP

中国工程热物理学会 多相流 学术会议论文 编号：056075 风沙运动的D酬数值模拟 亢力强 郭烈锦“ (蘸安交通人学动力工程多相流国家重点实验室，西安，710049) (nl：02E—mail：kJq@sohu．com) I翼一采用离教单元法(DEM)对风沙气河随相流动进行了数值计算。在该模型中，采用体平均的 Navic卜stokes方程米描述气相的运动，对离散颗粒采用拉格朗日方法模拟，求解颗粒运动方程。采用 硬球模型描述颗粒问碰撞作用。成功地模拟了风沙运动的动态过程。计算结果表明颗粒的平均水1r速 度随高度的变化町以表示为按对数函数或幂函数规律增加，输沙通量随高度按指数规律衰减，这与已 有文献的结果一致。本文的模拟宙助于研究风沙运动规律。 哭蝴tDEM模拟{硬球模型；风沙运动 1前言 ‘ 风畋颗粒运动可以导致一些严重的环境问题，例如沙尘暴和荒漠化ll“j。风沙运动 的研究受到人{f』越来越多的关注，一些数值模型也已报道Ij4】。在这些模型中．风沙运 动通常被分解为四个子过程：空气动力学卷吸颗粒、颗粒运动轨迹的计算、颗粒床面碰 撞以及颗粒和气流动量交换引起的风速修正，通过关联这些子过程，建立起负反馈机制 的非稳态过程，其描述颗粒床面碰撞的溅射函数是通过对单颗粒碰撞床面的统计计算而 来。McEwan和晰lle似n99I。19931f6|7】也描述一个跃移颗粒群的计算模型，利用高速摄 移过程的模拟，分析了不稳定风场对颗粒跃移行为的影响，发现输沙率的非稳态响斑与 风波动的幅度和频率有关。 上述这些模型中，颗粒床面作用通过溅射函数p4】或实验数据f6州进行描述，而不考 虑颗粒在运动过程中颗粒闻碰撞作用。 风沙运动是一个复杂的气同两相流动，在颗粒间蚍及气流和颗粒间都存在着复杂的 的沙床形态和颗粒运动轨迹。这表明了颗粒碰撞在风沙模拟中有重要作用。sorcnsen和 地会修改颗粒运动行为，降低颞粒床面碰撞强度。因此，颗粒床面碰撞以及气流中碰撞 在风沙运动模拟中不应该忽略。 另～方面，在气同流化床领域应用的离散颗粒方法…。川能在单颗粒层次上给出颗粒 ’目家重点书础研究发腥计划(973计划)G2000048702项目资助内容 “通讯作者，E—mail：U—guo@mamxjtu．cdu．c11 413 运动的详细信息，可以描述风沙运动中颗粒的随机运动。 本文建立了描述风沙运动的离散颗粒模型，采用硬球模型描述颗粒间碰撞作用。本 模型不仅考虑颗粒间碰撞佟用，俪且还考虑了流体颗粒闻的相互耦合作用。本工作主要 目的是基丁cFD的观点，在单颗粒层次上揭示风沙运动的动态变化规律和非均匀特征。 2数学模型 2．1流体相控镧方程 流体相连续方程和动量守恒方程分别表示如’F： 喜@，乃)+V·(口，乃u，)=o (1) (2) 昙(艺乃u，)+V·(町p，u，u，)=～，印+V·位，t，)+口，乃g—k 式中口，为流体体积分数，k为体平均的颗粒流体阻力·分别表示为 a，=l一∑％，△矿， ‰=古喜‰ 式中△矿和％分别为计算控制体的体积和颗粒i在该控制体中的体积·F．慨为单颗粒 所受的流体粘性阻力。，z为控制体中的颗粒数。对于二维流动计算，控制体体积取为 △矿=缸匈砬p，△x和缈分别为控匍体在x和y方向上的长度，dP为颗粒直径a 由于气流流动处于湍流状态，故流体相采用标准七一s湍流模型描述。 2．2囊粒运动方程 单颗粒的运动采用牛顿第二运动定律描述。颗粒运动方程为 脚，誓=郴+‰ (3) 采用DiFelice修正的流体阻力表达式计算单颗粒所受的流体粘