非稳态沙波纹流场的数值模拟研究.pdfVIP

第25卷增刊 计算力学学报 V01．25，Sup． 2008年12月 Med恤缸 December2008 0fC‘唧岫如nal Cl抵,eJo．mal 非稳态沙波纹流场的数值模拟 佟鼎’， 蒋红，黄宁 (兰州大学西部灾害与环境力学教育部重点实验室，兰州730000) 摘要：风成沙波纹是沙质地表在风力作用下形成的最基本的地貌类型．本文研究了沙波纹在形成过程中，风 场对波纹的形态变化的响应．采用目前比较流行的有限体积法，通过求解非稳态情况下的NavieT-Stokcs方程 来计算风场：地表形态采用阶梯网格进行逼近，为了对地表形态变化进行有效的控制，根据相关参数采用正弦 函数作为沙波纹的基本形状：最后根据计算结果详细讨论了沙波纹从形成到稳定的发展过程中地表流场的行为 变化特征． 关键词：沙波纹；simpler算法；非稳态；有限体积法 中国分类号：TP273 文献标识码：A 行模拟研究，使得现有的风沙跃移运动理论与数 1 引 言 值模型还远未达到能够准确预测野外实际情况下 的风沙流运动的程度。基于此，本文研究了沙波 风成沙纹是风力作用下地表形成的最基本的 纹在形成过程中，风场对波纹的形态变化的响应。 地貌类型，它既与沙粒输移过程紧密联系，又是 理解其他类型风沙地貌的基础。风成沙纹虽然具 2数学模型 有规则的几何图案且一般表现为迎风坡缓，背风 坡陡，但由于自然条件的千变万化，其形态又决 2．1流体控制方程 不是整齐划一的。风与不同粒级颗粒的相互作用 直角坐标系下流体二维基本微分方程组的统 也使沙纹的外部形状和内部构造变得复杂多样。 一形式为： 风成沙纹是风气流与沙质地表相互作用的直接表 达，其运动是风沙运动的一个重要的组成部分， 0r 对风成沙纹的研究有助于风沙动力学【14】的探讨研 其展开形式为 究。在风成沙波纹的研究中，虽然实验观测和测 塑+塑堕+塑鲤：旦(r型1 量一直是最为重要和常见的研究手段[5-s]，但随着 a 缸 勿 反、玉’ (2) 计算机技术的发展，数值模拟方法逐渐被重视。 +旦仃翌)+s 砂。砂 目前，对沙波纹形成的模拟研究取得了一定 式l中，方程左边第一项为瞬态项，第二项为对 的成绩。但是在大部分的研究方法中，缺乏对地 流项，方程右边第一项为扩散项，第二项为源项。 表流场有效信息的描述【9-11】。同时，目前对近地层 当矽取不同的值，方程1可表示为连续性方程、动 风沙流的理论与数值模拟研究都是在理想情况， 量方程和能量方程f1舶。 如定常风速、平坦沙面等给定条件下对风沙流进 2．2边界条件及数值算法 边界条件：入口采用速度进口，上边界为对 收稿日期：2008