坡面流格子boltzmann方法与preissmann隐式差分法-农业机械学报.pdfVIP

２０１４年 １０月 农 业 机 械 学 报 第 ４５卷 第 １０期 ｄｏｉ：１０．６０４１／ｊ．ｉｓｓｎ．１０００１２９８．２０１４．１０．０２１ 坡面流格子Ｂｏｌｔｚｍａｎｎ方法与Ｐｒｅｉｓｓｍａｎｎ隐式差分法模拟 １ ２ ３ ４ 张小娜 　冯　杰 　张东辉 　刘宁宁 （１．南京信息工程大学水文气象学院，南京 ２１００４４；２．中国水利水电科学研究院水资源研究所，北京 １０００４４； ３．江苏科技大学船舶与海洋工程学院，镇江 ２１２００３；４．山东省水利工程管理局，济南 ２５００１３） 摘要：以Ｄ１Ｑ５速度模型为例，时间多尺度分析为手段通过待定系数法来确定平衡态分布函数，将格子 Ｂｏｌｔｚｍａｎｎ 方法应用于坡面流运动方程；通过理想算例，以解析解为标准，比较了格子 Ｂｏｌｔｚｍａｎｎ方法与应用广泛的Ｐｒｅｉｓｓｍａｎｎ ４点隐式差分法的计算精度。研究表明，对于模拟坡脚断面水深和单宽流量过程，格子 Ｂｏｌｔｚｍａｎｎ方法的 Ｄ１Ｑ５模 型的计算精度高于 Ｐｒｅｉｓｓｍａｎｎ４点隐式差分法，尤其在退水阶段。对于达到平衡时间之前的坡面水深和坡面单宽 流量，格子 Ｂｏｌｔｚｍａｎｎ方法的Ｄ１Ｑ５模型的计算精度也高于 Ｐｒｅｉｓｓｍａｎｎ４点隐式差分法。但对于达到平衡时间之后 的坡面水深和坡面单宽流量，格子 Ｂｏｌｔｚｍａｎｎ方法的Ｄ１Ｑ５模型在ｘ＝１ｍ点上出现了较大的相对误差，计算精度逊 于 Ｐｒｅｉｓｓｍａｎｎ４点隐式差分法。将格子 Ｂｏｌｔｚｍａｎｎ方法应用于坡面流运动时，弛豫时间选择范围以［１，１２］ｓ为宜。 关键词：坡面流　格子 Ｂｏｌｔｚｍａｎｎ方法　Ｐｒｅｉｓｓｍａｎｎ隐式差分法　Ｄ１Ｑ５速度模型　弛豫时间 中图分类号：Ｐ６４２１１；ＴＶ１２２ 文献标识码：Ａ 文章编号：１０００１２９８（２０１４）１００１３２０８ 的关注［５－１２］。已有研究人员成功地运用格子 　　引言 Ｂｏｌｔｚｍａｎｎ方法，对一些复杂的流动与传热问题进行 土壤水蚀是一个世界性的环境问题，危害十分 了卓有成效的研究，但迄今有关格子 Ｂｏｌｔｚｍａｎｎ方 严重，容易引起土壤肥力减弱、农作物产量下降、河 法应用于坡面水流运动的研究寥寥无几，格子 床泥沙淤积、水质污染等。坡面水流是土壤水蚀过 Ｂｏｌｔｚｍａｎｎ方法和 Ｐｒｅｉｓｓｍａｎｎ隐式差分法在坡面水 程的主要动力，准确模拟坡面水流运动过程，对土壤 流运动过程的应用对比研究更是寥寥无几。本文拟 侵蚀发生、发展机理研究，进而对其进行预报具有重 运用格子 Ｂｏｌｔｚｍａｎｎ方法和 Ｐｒｅｉｓｓｍａｎｎ隐式差分法 ［１］ 要的意义 。 求解坡面流运动方程，通过算例进行对比分析并深 坡面水流运动属于流体运动的范畴，水流运动 入探讨。 是大量粒子微观运动的宏观表现。传统的计算流体 １　坡面流运动方程 力学方法通常建立在宏观连续介质假设的基础上， 即将流体假设为连续介质并用 Ｎａｖｉｅｒ Ｓｔｏｋｅｓ方程 坡面流运动过程可用运动波方程进行描述，即 组及其简化形式作为控制方程，然后采用有限差分、 ｑ ｈ   ＋ ＝ｒ ｘ ｔ 有限体积、有限元等离散方法把描述流体运动的连