《高浓度固液两相流地运动特性研究》由会员分享，可在线阅读，更多相关《高浓度固液两相流地运动特性研究（12页珍藏版）》请在金锄头文库上搜索。
1、word高浓度固液两相流的运动特性研究倪晋仁1,2，黄湘江1,2(1.大学环境科学中心；2.水沙科学教育部重点实验室)摘要：利用固体颗粒运动的动理论，通过改变颗粒浓度可以考察非粘性颗粒在水流中运动的典型微观和宏观运动特性。本文分别对微观的颗粒速度分布函数变化和由此衍生的诸如颗粒平均速度、颗粒脉动速度和单位体积颗粒数垂线分布等宏观变量的变化进展了系统比拟。研究结果明确：动理论能够比传统理论获得更详细的微观和宏观信息，也更适合研究高浓度固液两相流运动特性，颗粒运动微观和宏观特性在颗粒浓度超过一定阈值后会发生本质的变化，但临界颗粒浓度值(阈值)在不同的计算和实验条件下会有一定的差异。关键词：高浓度挟
2、沙水流，微观，宏观，特性，运动学理论基金项目：国家自然科学基金资助项目(49625101)作者简介：倪晋仁(1963-)，男，某某山阴人，教授，主要从事环境科学与泥沙方面的研究。高浓度固液两相流在生产实践中经常遇到。河流中的泥沙含量高，可能导致河道淤积、河床抬高和洪水频率增加1。高浓度固液两相流的流动和输运特性与低浓度固液两相流有着很大的不同。高浓度挟沙水流经常表现出非牛顿流体的特性2，不同于低浓度时的牛顿流体。以往对于高浓度固液两相流的描述多基于宾汉塑性体模型或拜格诺的膨胀体模型3，4。就含有粘性颗粒的高浓度固液两相流而言，中国学者提出了许多关于屈服应力和宾汉粘性系数的经验表达式，这些表达式
3、中大都采用颗粒浓度和反映颗粒大小组分的变量。Chen5曾对这方面的研究工作进展了全面的评述。就含有非粘性颗粒的高浓度固液两相流而言，以往的研究6多从Bagnold3的颗粒离散应力概念出发。Chen7的粘塑体模型包含了以上两种情况。最近，新的流变模型研究又有进展，并用于描述高浓度挟沙水流的复杂特性，参见Chen8和Brufau9等。通常描述固液两相流的连续介质理论10能够合理地描述流体和颗粒的宏观运动特性，但不能充分解释颗粒与颗粒的相互作用，更不能描述颗粒运动的微观特性。采用基于Boltzmann方程的动理论能够很好地描述个体颗粒运动和颗粒之间相互作用的微观特性。这个方法类比自气体分子运动论，一
4、旦微观的颗粒速度分布函数，固液两相流的微观和宏观特性都可得到很好的认识。尽管动理论过去多被用于描述低浓度固液两相流，近年来该方法已经被用于高浓度固液两相流研究中。例如，倪晋仁和王光谦1114曾应用动理论研究高浓度固液两相流中悬浮颗粒垂向分布。本文如此将动理论的应用扩展到研究高浓度固液两相流的主要微观运动特性(如颗粒运动速度分布函数变化)和宏观运动特性(如颗粒平均速度、颗粒脉动速度和单位体积颗粒数等)。为此，从颗粒速度分布函数的微观信息入手，探讨颗粒浓度由低向高变化时固体颗粒特性的响应变化。1 颗粒微观和宏观特性的主要变量王光谦和倪晋仁15，16曾在低浓度固液两相流研究中引入了动理论，并建议固相
5、颗粒类比气体分子运动用Boltzmann方程(1)来描述。方程右边是反映颗粒碰撞影响的积分项。在低浓度固液两相流中，颗粒碰撞影响较小，积分碰撞项通常被忽略。颗粒速度分布函数f=f(vi,xi,t)是在空间坐标xi和时间t颗粒速度介于vi和vi+dvi的颗粒数目，其中dvi=dv1dv2dv3，vi是颗粒的随机速度，Fi是作用在颗粒上的单位质量力，它包括重力和液相的作用力。颗粒速度分布函数f能够很好地反映颗粒运动的特性。颗粒速度分布函数的任何变化都将引起一系列颗粒运动宏观特性的变化12，13。例如，单位体积颗粒数目为n=fdvi(2)相密度为=mfdvi(3)其中m是单个颗粒的质量。颗粒平均速度
6、为i=1/mvifdvi(4)颗粒的脉动速度为vi=vi-i(5)颗粒脉动速度的均方值为2=1/m(vi-i)(vi-i)fdvi(6)2 颗粒运动的微观特性颗粒速度分布函数是颗粒运动最重要的微观特性。王光谦和倪晋仁16在平衡条件下针对低浓度固液两相流，通过忽略式(1)中的复杂碰撞项得到了一个与Champan & Cowling所推导形式类似的方程。当积分碰撞项很小并可被忽略时，得到了相应于低浓度固液两相流情形下(颗粒体积浓度小于0.05)的颗粒速度分布函数f0=n0exp-/2u2L/u2*3/2L/pH/DCd(1-Cl/Cd)(1-u)2-(u-u0)2-/2u2L/u2*3/2L/PH
7、/DCd(1-Cl/Cd)(1-u)2-(u-u0)2du(7)其中，u=v/uL和u0=v0/uL./u*)/1+，(8)u0/u*(9)n0=naexp(-A*(-a)，(10)A*=15(1-v)n1/u*(11)式中：v0为颗粒参考速度；为颗粒沉降速度；u*为剪切速度；D为颗粒粒径；H为水深；uL为流体的平均流速；Cd为综合阻力系数；Cl为综合升力系数；P为颗粒密度；L为流体密度；为颗粒垂线平均的体积浓度；a为参考点位的相对高度；g为重力加速度;n1为指数。在高浓度固液两相流中碰撞项不能再被忽略。在这种情况下，可以采用与分子气体运动论中类似的方法进展简化=-f-f0/(12)其中，f0
8、是低浓度条件下的颗粒速度分布函数，与高浓度下的颗粒速度分布函数f不同；为松弛因子，可以表示为=l/vz(13)其中，l与固液两相流中固体颗粒在垂向运动的平均自由程成正比，可近似地视为一常数(实际可能是颗粒特征的函数)；vz是等向速度场中的颗粒速度。考虑恒定流动，如此Boltzmann方程简化为f=f0-1/vz(vi-Fi)(14)作为一阶近似，再次采用气体分子运动理论中的处理方法12,14，可以得到高浓度和低浓度固液两相流条件下的颗粒速度分布函数之间的关系f=f0-l-mglf0(15)或者f=f0-L-ALf0(16)其中，L=l/H，=z/H，A=kmgH，k为反映颗粒浓度影响的系数；H
9、为水深；为Lagrange系数15。为与作者以前的工作具有可比性和一致性，本文仍采用Michalik17的实验作为计算条件。Michalik的实验是在内径为200mm的方管中进展的，挟沙水流由密度为p3的固体颗粒和密度为L3的水组成的固液混合流。颗粒的代表粒径为d50=0.45mm。颗粒的垂线平均体积浓度(2/s,对应的颗粒雷诺数Re=VD=8105。在给定的颗粒粒径和温度条件下，颗粒沉降速度为=6.15cm/s。液相流体的速度分布可以根据修正的窦国仁公式(1987)计算，即uL/u*=2.5ln1+u*z/5+7.05(u*z/1+u*z/)2+2.5(u*z/1+u*z/)+0.51-co
10、s(z/H)(17)这一处理方法可能在高含沙量条件下带来一定的误差，但根据作者采用其它流速分布公式进展比拟，对所得结论不会有明显影响。在具体计算时，这方面值得继续改良。当颗粒浓度很高时，流体相关参数的测量存在很大困难。在求解方程过程中，边界条件假设符合如下公式fb=Nb1/exp(-(u-ub)2/22b)(18)其中，fb为床面位置处的颗粒速度分布函数；Nb为床面处的单位体积颗粒数；b为床面处颗粒速度分布函数的标准偏差；ub为床面处颗粒分布函数的平均速度。在过去研究成果L=3的根底上，进一步考虑高浓度条件下颗粒浓度的影响，并采用经验系数k=30(-0.5)2进展修正。L的大小间接反映与平均自
11、由程的比例关系，与颗粒连续两次碰撞所用时间有关。取床面作为参考点，即a=0,便可通过求解上述诸方程得到相应的数值解(见图1和图2)。假如定义颗粒速度概率密度分布函数为p(u)=(1n)f，如此会相应用到两个参数，一个是对应颗粒速度概率密度分布函数取最大值时的特征速度uc，另一个反映概率密度分布函数形态的标准偏差。当颗粒浓度给定时，图1和图2仅给出了沿垂直方向上相对水深为z/H=0.1，0.5和0.9的局部计算结果；颗粒平均体积浓度的变化如此给出了=0.15，0.27，0.31，0.42和0.54的代表情形。图1 不同颗粒浓度条件下的概率密度分布函数p(u)=(1/n)f变化分析明确，尽管对应于
12、最大颗粒速度概率密度分布函数值的特征速度uc从床面向上呈不断增加的趋势，但在垂直方向的任一位置上根本不受颗粒平均浓度变化的影响，见图3。然而，颗粒平均浓度的变化却明显地影响着标准偏差的变化。由图4可见，沿垂线任意高度上，随的变化可以大致分为三个阶段：即对应于0.42的第三阶段。在第一阶段，随的增加而缓慢减小；在第二阶段，随的增加而根本不变；在第三阶段，随的增加而急剧减小。 应该指出，不同阶段划分的阈值=0.31和0.42本身仅仅是在本文的实验条件下获得的，在其它实验条件下可能会有所不同。但是，这里所揭示的高浓度固液两相流中颗粒浓度超过一定阈值后发生的微观运动特性变化却是非常具有启发性的。这一颗
13、粒运动的微观特性很难从一般的连续介质理论得到。从动理论获得的微观运动特性不仅为理解颗粒运动的细节增加了新的内容，而且也为研究如下颗粒运动的宏观特性奠定了理论根底。3 颗粒运动的宏观特性3.1 颗粒平均速度垂直分布 颗粒运动的平均速度是一个十分重要的宏观特性，可根据式(4)的定义和式(16)给出的颗粒速度分布函数表达式求解，计算结果见图5。图2 不同垂直高度下的概率密度分布函数p(u)=(1/n)f变化可以看出，颗粒平均速度随颗粒浓度的变化根本表现为：当颗粒浓度没有超过一定阈值前，颗粒平均速度的垂直分布与常规的固液两相流垂直分布变化相似；但是，随着颗粒浓度的继续增加，颗粒平均速度因颗粒相互作用形
14、式的变化而快速增加，突出表现为颗粒向上聚集趋势的增加和最大颗粒平均速度出现位置的上移。二者的界限大致在颗粒浓度=0.42左右，这与颗粒运动微观特性分析中讨论的第二个阈值相近。3.2 颗粒脉动速度垂直分布 颗粒脉动速度的垂直分布可根据式(6)的定义和式(16)给出的颗粒速度分布函数求解，计算结果见图6。由图可见，在颗粒浓度0.31的X围内，颗粒脉动速度的垂直分布随颗粒浓度的增加而缓慢变化；当0.310.42时，颗粒脉动速度随颗粒浓度的增加急剧增加，特别是对于靠近水流外表的区域。高浓度固液两相流中颗粒脉动速度发生本质性变化的阈值与前面讨论的颗粒速度分布函数变化对应的阈值具有较好的对应关系。3.3
15、颗粒浓度垂直分布 颗粒浓度(或颗粒数)的垂直分布可以很好地反映颗粒运动的状态，是颗粒运动分析时不可缺少且相对容易测量的重要宏观特性。颗粒浓度(或颗粒数)可根据式(2)的定义和式(16)给出的颗粒速度分布函数表达式求解，也可以对式(16)在速度空间上直接积分得到。仍然采用Michalik17的测量数据进展计算，可以得到如图7所示的颗粒浓度分布。由图可见，沿垂线上对应于最大颗粒浓度的位置随颗粒平均浓度的增加而显著上升，这有别于一般的低浓度固液两相流。发生显著变化的对应颗粒平均浓度大约为0.42。图3 特征速度uc随颗粒浓度和位置的变化图4 标准偏差和平均颗粒浓度的关系图5 颗粒平均速度垂向分布随平均颗粒浓度的变化图6 颗粒脉动速度随颗粒平均浓度的变化由于对高浓度固液两相流中微观特性的直接测量在技术上存在很大的困难，因此要对颗粒速度分布函数进展直接验证是不易做到的。鉴于此，我们可以用测量到的颗粒运动宏观特征之一颗粒浓度垂直分布进展间接验证。由M
