《高浓度固液两相流中泥沙分布的修正公式》由会员分享，可在线阅读，更多相关《高浓度固液两相流中泥沙分布的修正公式（5页珍藏版）》请在金锄头文库上搜索。
1、1高浓度固液两相流中泥沙分布的修正公 式摘要：本文在过去低浓度固液两相流动理论研究的基础上, 提出了对高浓度固液两相流中泥沙浓度分布的一个修正模型，模型的预测结果与有关实测数据较为吻合。 关键词：高浓度固液两相流 修正 浓度分布 1 前言迄今为止，关于固液两相流中颗粒浓度分布规律的研究已有很多，人们提出了许多理论和公式 (Rouse, 1937; Itakura Kishi, 1980; Hunt, 1954)。其中，倪晋仁和王光谦(1989)提出的统一公式，能较好地概括各种公式。实际上，实测结果多数表现为两种浓度分布，其中一种是人们较为熟知的底部浓度最高而顶部浓度最低的分布模式; 另一种则是
2、容易被人们忽视的，即最大浓度出现在底部以上的分布模式。采用常规的理论将很难解释上述两种浓度分布模式的存在。迄今为止，几乎所有的理论都无法定量来描述它们，这意味着我们必须采用新的方法来描述它们。为此，Wang Ni (1991)采用动理论方法求解了低浓度两相流中的泥沙浓度分布。由于问题的复杂性，过去对于高浓度情形进行的理论研究相对较少。在已有的少数研究中，Bagnold (1954) 的研究和 Michalik (1973)的工作是很有参考价值的，这也为以后的研究奠定了基础。2 修正模型通过类比气体分子运动理论中的相关处理方法，固体颗粒的运动可用Boltzmann 方程描述。这样，单位体积的颗粒
3、数沿垂向的分布公式可表示为(Wang Ni, 1991)2（1）其中 A 为未知参数； 是与 Lagrange 算子有关的参数；,Z*为与水沙运动特征相关的参数。注意到（2）（3）（4）其中 3（5）（6）（7）在上述各公式中，a 代表垂线上参考点的相对坐标;=y/D，相对高度; D, 特征长度尺度，对于封闭圆管 D 为管道内径，对于矩形断面的管道 D=h (h 为水深的一半)。对(1)式求导并令其为零，可以得到垂线最大颗粒数对应的相对坐标m=(-1)/Z*(8)这意味着最大浓度的分布位置取决于 和 Z*。当 =1 时，垂线最大颗粒数的位置对应于床面, 这时浓度分布为 II 型; 对于较大的
4、来说，垂线最大颗粒数的位置在床面以上, 这时浓度分布为 I 型。在悬移质运动研究中，我们知道 Z 与 /U*成线性关系。 是 f/p 和 d/D 的函数。这里 Z*与通常所说的悬浮指标 Z 类似，尽管它们的确切含义不尽相同。对于低浓度固液两相流，Z*和 可表示为4Z*=5(/U*(9)(10)其中 是颗粒的沉降速度；d 为颗粒直径，f,p 分别为液相和固相的密度；U*为摩阻流速。对于低浓度水沙流来说，固体颗粒间的作用不太明显，泥沙的悬浮主要依靠水流的紊动作用；对于高浓度水沙流来说，固体颗粒间的作用变得越来越显著，因此必须对反映颗粒所受上升作用的参数 进行修正。根据 Bagnold (1954)
5、 的理论，颗粒间的作用可通过以下参数 0 来反映(11)其中 是颗粒的平均体积浓度；Cvm 是固液两相流的最大可能浓度，一般取为0.74。对于高浓度流动情况，可将参数 进行相应修正, 如根据实测资料将之修正为 5(12)便可结合以上各式得出高浓度固液两相流垂向浓度分布的修正模型。 3 模型的验证为对现行模型进行验证，我们采用了 Michalik(1973)的高浓度固液两相流的实测数据。Michalik 采用的固液混合物由沙 (p=2.65g/cm3)和水(f=1.00g/cm3)组成，平均流速=4.00m/s。沙的特征直径(中值粒径)d50=0.45mm。Michalik 的试验在一个内径为
6、D=200mm 的工业用管道中进行，唯一的可变参数是混合物的平均体积浓度 ，取值范围为 0.15 至 0.54。设水温为 1520，动力粘滞系数为=0.01cm2/s，雷诺数为 Re = D/=8105。根据沙的直径和水温，可以得出对应泥沙沉降速度 =6.15cm/s，由 Chezy 公式(13)结合 Nikuradse 公式6(14)通过迭代计算，得出 b=0.0121 或 U*=0.0389, =0.1556 m/s。其中, b 为阻力系数，C 为 Chezy 系数。 通过给定的参数，可以给出计算的浓度分布，计算值与 Michalik 实测数据的比较如图 1。由于目前很难十分精确地从理论上
7、来描述高浓度固液两相流，并且由于测量上的困难，所以图中所示的分析解和实测值存在的差异应该是可以容许的。4 结论固液两相流中颗粒运动的动理论为低浓度和高浓度的两相流研究提供了很好的理论基础。对于高浓度的情况，必须对其依据的低浓度模型进行修正。修正后的模式与实测结果较为吻合。图 1 公式与实测值的对比验证参 考 文 献 1 BAGNOLD, R.A. (1954), Experiments on A gravity-free dispersion oflarge solid spheres in Newtonian fluid under shear, Proc., Royal Soc. Lond
8、on, Ser. A, 225,49-63. 2 CHAPMAN, S. COWLING, T.G. (1970), The Mathematical Theory of Non-Uniform Cases, 3rd ed. Cambridge University Press. 3 Hunt, J.N. (1954), The turbulent transport of suspended sediment in openchannels, Proc. R. Soc. London, 224(1158), 322-335. 4 Itakura, T. Kishi (1980), Open
9、channel flow with suspended sediments, Jour. of Hydraulic. Div., Am. Soc. Civ. Eng., 8,1325-1343. 5 MICHALIK, A. (1973), Density patterns of the inhomogeneous liquids in the industrial pipelines measured by means of radiometric scanning. La Houille Balanche, No.1, 53-59. 6 NI, J.R. (1989), A general
10、ized formula on the vertical distributions of suspended sediment concentration, the Fourth International Symposium on River Sedimentation, November 1-5, China, 603-610.7 倪晋仁, 王光谦。论悬移质泥沙浓度垂线分布存在的两种类型及其产生原因。水利学报, 1987，(7)： 60-68. 78 OWEN, P.R. (1969), Pneumatic transport, Jour. of Fluid Mechanics, London39, 407-432. 9 WANG, G.Q. NI, J.R. (1991), Kinetic theory for particle concentration distribution in two phase flow, Jour. Engrg. Mechanics, 116(12), 2738-2748.
