风沙运动问题的sph-fvm耦合方法数值模拟研究 - 物理学报.pdfVIP

物理学报 Acta Phys. Sin. Vol. 63, No. 13 (2014) 130202 风沙运动问题的SPH-FVM耦合方法 数值模拟研究∗ 陈福振 强洪夫 高巍然 (第二炮兵工程大学601 室, 西安 710025) ( 2014 年1 月11 日收到; 2014 年3 月26 日收到修改稿) 针对离散颗粒模型和欧拉-欧拉双流体模型在求解气粒两相流动问题中存在的不足, 提出了一种新方法 ——SPH-FVM 耦合方法, 并应用于风沙运动过程的数值模拟计算. 新方法基于拟流体模型, 采用光滑粒子流 体动力学方法(SPH) 对离散颗粒相进行求解, 追踪单颗粒运动轨迹, 采用有限体积法(FVM) 求解连续气体 相, 捕捉流场特性, 两相间通过曳力、压力梯度、体积分数等参量进行耦合, 建立了两种不同坐标系下方法间的 耦合框架. 对SPH 粒子所承载的物质属性进行了重定义, 改造成了适用于离散颗粒相求解的光滑离散颗粒流 体动力学模型(SDPH), 阐述了SPH 粒子与离散相中颗粒之间的关系, 推导得到了拟流体SPH 离散方程组. 模拟了风沙运动中沙粒跃移过程和自由来流风速作用下沙粒的运动过程以及沙丘在风力作用下缓慢向前蠕 移的过程, 分析了颗粒的运动轨迹, 流场中沙粒水平速度分布规律以及气体速度场在沙粒反作用下的变化情 况等, 与实验结果相符合, 结果表明该方法不仅精度较高, 而且计算量较小, 适于求解风沙运动问题乃至其他 气粒两相流动问题. 关键词: 光滑粒子流体动力学, 有限体积, 耦合方法, 风沙运动 PACS: 02.70.–c, 47.11.–j, 47.61.Jd, 47.57.Gc DOI: 10.7498/aps.63.130202 运用该技术测量了风洞中跃移颗粒的速度及能量 1 引 言 分布, 得出了颗粒均值速度随高度呈幂指数函数变 化规律. 尽管高速摄像技术已经广泛用于单颗粒 风沙运动由于易造成沙漠化, 导致土壤风蚀、 粒子运动速度的测量, 但是这种方法存在的最大的 沙尘暴等一系列重大的环境问题, 逐渐引起人们的 问题为仅对尺寸较大的颗粒有效. 对处于蠕动中 高度关注. 目前对于风沙运动的研究主要有野外观 5 mm 以下高度以及稀疏流动中70 mm 以上高度的 测、实验测量、数学建模与数值模拟四种方式. 在 沙粒的速度较难获得. 随着多普勒粒子分析技术的 实验研究方面, 最早由Bagnold 1 借鉴空气动力学 发展, 该问题逐渐得到了解决. Dong 等5 运用粒 的分析手段, 开展了大量的野外观测与风洞实验 子动态分析仪(PDA)对风洞中颗粒掠过沙床表面 研究, 奠定了风沙物理学研究的基础. 高速摄像技 时速度的概率分布、均值速度随高度的变化关系以 术作为一种有效测量颗粒速度的方法, 已经用于 及颗粒的湍动能等进行了研究, 结果表明迎风速度 单颗粒轨迹、沉积体破坏、表层塌滑现象以及颗粒 呈现一种高斯函数分布形态, 垂直速度由于颗粒撞 与颗粒间碰撞等问题的理解. White 和Mounla 2 击壁面以及颗粒与颗粒之间的碰撞作用分布复杂. 研究了弗劳德数对于风洞中沙粒跃移过程的影响. Kang 等6 采用相位多普勒粒子分析仪(PDPA)对 3 Greely 等 通过野外实验, 运用高速摄像技术测量 风洞中沙床表面沙粒撞击和起跳速度概率分布以