风沙跃移运动的数值模拟.pdfVIP

学兔兔 第1期(总第188期) 机械 工 程 与 自动 化 No．1 2015年O2月 MECH ANICAI ENGINEERING AUTOMATION Feb． 文章编号：1672—6413(2015)01—0004～03 风沙跃移运动的数值模拟鬻 刘 博，黄新成，王旭峰，蒋建云，何欢欢 (塔里木大学机械电气化工程学院，新疆 阿拉尔 843300) 摘要：根据纳维一斯托克斯方程，分析风沙流的特征和风沙运动的基本规律，推导出流体相控制方程。采用 离散单元法和硬球模型对单颗沙粒碰撞规律进行分析。微观上通过建立风沙运动颗粒碰撞的数学模型，如跃 移颗粒在空间上的碰撞等，得出颗粒的运动规律，同时模拟出颗粒的受力状况、沙粒速度随高度的变化规律 以及输沙量随高度的分布状况。 关键词：风沙运动；硬球模型；输沙量；颗粒碰撞；数值模拟 中图分类号：TP15 文献标识码：A 0 引言 由于沙粒粒径较大，且较为干燥，沙粒间的液膜黏 2O世纪3O年代，美国和加拿大西部大草原的沙 结力和静电力可忽略不记，沙粒近似为无黏性流体，又 尘暴灾害受到全世界的高度重视。基于这次沙尘暴灾 因为跃移层中的沙粒拖拽力可以改变风速分布，其对 害，Bagnold¨】 将风沙运动作为一个空气动力学来研 风速的作用可以近似认为是纳维一斯托克斯方程式添 究，于 1941年发表了关于风沙物理学的著名论著。他 加一个修正项F ，即： 通过利比亚沙漠中的野外观测和实验室风洞试验确定 。 l0，[ +(，· )，]一P，g— 夕+ r，+Fdra ．(3) 了沙粒运动的力学机制，指出了沙粒的运动主要发生 L， 在0 m～1 ITI的高度范围内，在10 cm～2O cm的高度 1 Fd 一‘ ，dra ． (4 范围内较为集中 ]。在他的著作中，深刻地描述了未 i= 1 起沙地表和起沙地表上方的风速剖面、跃移沙粒的运 其中：△ 为控制体的体积；Fd 为单颗粒所受拖拽力； 动轨迹、沙粒运动的临界风力、风沙的单宽输沙率等物 ， ， 为控制体内第i个颗粒对沙粒的拖拽力；”为控 理学问题，对理性认识风沙运动具有重要意义。 制体内的颗粒数。 本文采用硬球模型建立了单颗沙粒的运动模型， 1．2 颗粒碰撞模型 并充分考虑颗粒间的碰撞作用，模拟了沙粒速度随高 颗粒间碰撞的瞬时性是硬球模型的基本思想，颗 度的变化以及输沙量随高度的分布状况。 粒本身在碰撞过程中是不会发生变形的。根据动量守 1 沙粒跃移运动模型 恒原理得到颗粒碰撞后的运动状态模型的突出特点 1．1 流体相控制方程 是：通过给定摩擦因数和恢复系数就能够求出碰撞后 近地层风沙运动满足纳维一斯托克斯方程，矢量 的颗粒运动速度，此模型适合低浓度颗粒流的计算。 式为： 沙粒在运动过程中不断地与其他沙粒发生碰撞， pr[ +(，· )，]=