沙粒跃移运动的数值模拟论文.pdfVIP

第25卷增刊 计算力学学报 V01．25，Sup． 2008 2008年12月 Medl锄ics December Ch戤Jo,,rnalofCc删pI|伽i伽al 沙粒跃移运动的数值模拟 马高生， 郑晓静’ (兰州大学西部灾害与环境力学教育部重点实验室，兰州730000) 摘要：沙粒在风场作用下形成的风沙运动，直接导致一系列环境问题．现有的数值模拟并没有揭示风沙流的 二维性质．本文提出一种合理有效的模拟方法来揭示风沙流时空的变化．气相运动采用大涡模型，而固相【沙 粒)处理为高散体(DEM)，并利用已有的沙粒与床面碰撞的激溅函数，以获得整个风沙系统运动特性．结果 表明：其发展过程与以往的试验结果相比较符合：风沙流达到稳定后，不同截面输沙率沿高度分布与平均输沙 宰相比差别很大，并大约在9cm高处，平均输沙率出现极值，单宽输沙率沿风向分布同样具有非均匀性． 关键词：流体体积分数；相问阻力；风沙耦合；输沙率 中图分类号：TP273 文献标识码：A 沙粒跃移高度在0．02m以下，因此较大的计算网 1引言 格可能导致在对沙粒跃移和近地面风场的模拟方 面不够精细，如风沙流稳定后，近地面摩阻风速 对沙粒在风场作用下形成的风沙流动，由于 为0．32m／s，大于计算模型中的起动风速0．21m／s， 直接导致一系列环境问题，而得到学者们的关注。 不少学者通过野外观测和风洞实验以及理论建 采用k一占湍流模型处理湍流风场。在对沙粒相的 模，试图对近地表的风沙流的规律有清楚的认识 处理方面，他们追踪单颗沙粒与床面的碰撞和跃 并能予以准确的预测。然而，由于野外环境的多 样性、风洞实验的局限性以及风沙两相流动问题 移运动，得到了输沙通量随高度的分布和风沙流 发展时间等结果。这种处理的优点在于：可以得 的复杂性，现有的实验和理论研究所给出的结果 主要集中在考虑风场【1．2】和输沙强度，如：浓度和 到每个沙粒运动信息，如沙粒运动速度分布和粒 输沙率等‘】，沿地表高度的变化，以及电场‘叼 床碰撞。但是，由于计算量较大，使之限制了所 能处理的沙粒数，由此导致计算结果与试验有较 方面。 大出入。例如，风速沿高度的分布规律与实测结 最早采用三维N-S方程研究风沙运动的sh∞ 果…J有较明显的差距。 等【8】，摒弃了仅适合一维的混合长理论，他们对瞬 尽管文献8，10处理的是二维和三维的风沙 态的N-S方程，采用大涡模拟，次网格湍流应力 流，但他们并没有揭示风沙流的不均匀性。我们 项采用k—s模型，在沿风向的水平方向上取网格 知道：即使不考虑实际风速的波动和地表性质的 长度为0．1m，垂向网格长度O．02m。计算了区域 差异，风速沿风向的分布也不是均匀的。这是因 为(