展 1999 年第 1 期.pdf

　 辑第 14 卷第 1 期　　　　　 水 动 力 学 研 究 与 进 展 　　　 . , . 14, . 1 A Ser A V o l N o 　　1999 年 3 月　　　　　　　 　　　　　 . , 1999 JOU RNAL O F H YDROD YNAM ICS M ar 任意流场中稀疏颗粒运动方程的 数值解法及其应用 黄社华　　李　炜 (武汉水利电力大学, 武汉 430072) 程良骏 (华中理工大学, 武汉 430074) 　　摘　要　任意流场中稀疏颗粒运动方程是一非线性的微分积分方程, 极难获得解析解。本文 首先对方程中Basset 力项构造了适当的计算格式, 然后详细研究了颗粒运动方程的数值计算方法 并提出称之为P (E C ) k 多步法的差分格式, 简要分析了该格式的收敛性和代数精度, 给出了时域离 散步长必须满足的条件。最后应用 P (E C ) k 格式实际计算均匀流场中和二维无界平板垂直射流驻 点区稀疏颗粒的运动, 得到了一些有意义的结论, 将均匀流的计算结果与其解析解的比较证明了 本文数值计算方法的正确性和良好精度。 关键词　固液两相流, 颗粒运动方程, P (E C ) k 数值方法 分类号　 359 O 1　引言 在自然界和现代工程技术领域中有相当多的一类流动现象属于稀疏固液两相流范畴, 如 水利水电工程和各种流体动力设备中固体颗粒的沉积、输运和对设备的磨损问题, 环境工程中 污染物的扩散问题, 单相和多相湍流激光测速技术中粒子的跟随性问题等, 无不与其中稀疏颗 粒复杂的运动特性有关。因此研究任意流场中稀疏颗粒的运动规律具有重要的理论和工程实 际价值。 在文献[1, 2 ] 中, 作者通过分析各种流动条件下稀疏刚性颗粒运动的相间作用力并进行 量级比较, 得到了任意流场中稀疏颗粒运动方程的一般形式, 此方程考虑的因素比较全面, 因 而具有较为广泛的适用性。为便于本文的研究, 现给出该方程的矢量形式如下: d u 3 p f f d u f ( ) ( ) d t = 1 - f + d t + CD 4D u - up u - up s s s d (u - u ) 9C t d (u - u ) d f p B f ef f p + Cm + ∫ d 2 d t D 0 s s t - ( ) 本文