高浓度固液两相流的运动特性研究

高浓度固液两相流的运动特性研究倪晋仁1,2，黄湘江1,2

(1.北京大学环境科学中心；2.水沙科学教育部重点实验室)摘要：利用固体颗粒运动的动理论，通过改变颗粒浓度可以考察非粘性颗粒在水流中运动的典型微观和宏观运动特性。本文分别对微观的颗粒速度分布函数变化和由此衍生的诸如颗粒平均速度、颗粒脉动速度和单位体积颗粒数垂线分布等宏观变量的变化进行了系统比较。研究结果表明：动理论能够比传统理论获得更详细的微观和宏观信息，也更适合研究高浓度固液两相流运动特性，颗粒运动微观和宏观特性在颗粒浓度超过一定阈值后会发生本质的变化，但临界颗粒浓度值(阈值)在不同的计算和实验条件下会有一定的差别。关键词：高浓度挟沙水流，微观，宏观，特性，运动学理论基金项目：国家自然科学基金资助项目(49625101)

作者简介：倪晋仁(1963-)，男，山西山阴人，教授，主要从事环境科学及泥沙方面的研究。高浓度固液两相流在生产实践中经常遇到。河流中的泥沙含量高，可能导致河道淤积、河床抬高和洪水频率增加[1]。高浓度固液两相流的流动和输运特性与低浓度固液两相流有着很大的不同。高浓度挟沙水流经常表现出非牛顿流体的特性[2]，不同于低浓度时的牛顿流体。以往对于高浓度固液两相流的描述多基于宾汉塑性体模型或拜格诺的膨胀体模型[3，4]。就含有粘性颗粒的高浓度固液两相流而言，中国学者提出了许多关于屈服应力和宾汉粘性系数的经验表达式，这些表达式中大都采用颗粒浓度和反映颗粒大小组分的变量。Chen[5]曾对这方面的研究工作进行了全面的评述。就含有非粘性颗粒的高浓度固液两相流而言，以往的研究[6]多从Bagnold[3]的颗粒离散应力概念出发。Chen[7]的粘塑体模型包含了以上两种情况。最近，新的流变模型研究又有进展，并用于描述高浓度挟沙水流的复杂特性，参见Chen[8]和Brufau[9]等。通常描述固液两相流的连续介质理论[10]能够合理地描述流体和颗粒的宏观运动特性，但不能充分解释颗粒与颗粒的相互作用，更不能描述颗粒运动的微观特性。采用基于Boltzmann方程的动理论能够很好地描述个体颗粒运动和颗粒之间相互作用的微观特性。这个方法类比自气体分子运动论，一旦微观的颗粒速度分布函数已知，固液两相流的微观和宏观特性都可得到很好的认识。尽管动理论过去多被用于描述低浓度固液两相流，近年来该方法已经被用于高浓度固液两相流研究中。例如，倪晋仁和王光谦[11～14]曾应用动理论研究高浓度固液两相流中悬浮颗粒垂向分布。本文则将动理论的应用扩展到研究高浓度固液两相流的主要微观运动特性(如颗粒运动速度分布函数变化)和宏观运动特性(如颗粒平均速度、颗粒脉动速度和单位体积颗粒数等)。为此，从颗粒速度分布函数的微观信息入手，探讨颗粒浓度由低向高变化时固体颗粒特性的响应变化。1颗粒微观和宏观特性的主要变量王光谦和倪晋仁[15，16]曾在低浓度固液两相流研究中引入了动理论，并建议固相颗粒类比气体分子运动用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是单个颗粒的质量。颗粒平均速度为i=1/ρ∫mvifdvi(4)颗粒的脉动速度为v′i=vi-i(5)颗粒脉动速度的均方值为′2=1/ρ∫m(vi-i)(vi-i)fdvi(6)2颗粒运动的微观特性颗粒速度分布函数是颗粒运动最重要的微观特性。王光谦和倪晋仁[16]在平衡条件下针对低浓度固液两相流，通过忽略式(1)中的复杂碰撞项得到了一个与Champan&Cowling所推导形式类似的方程。当积分碰撞项很小并可被忽略时，得到了相应于低浓度固液两相流情形下(颗粒体积浓度小于0.05)的颗粒速度分布函数图1不同颗粒浓度条件下的概率密度分布函数p(u)=(1/n)f变化分析表明，尽管对应于最大颗粒速度概率密度分布函数值的特征速度uc从床面向上呈不断增加的趋势，但在垂直方向的任一位置上基本不受颗粒平均浓度变化的影响，见图3。然而，颗粒平均浓度的变化却明显地影响着标准偏差σ的变化。由图4可见，沿垂线任意高度上，σ随的变化可以大致分为三个阶段：即对应于<0.31的第一阶段，对应于=0.31～0.42的第二阶段，和对应于>0.42的第三阶段。在第一阶段，σ随的增加而缓慢减小；在第二阶段，σ随的增加而基本不变；在第三阶段，σ随的增加而急剧减小。应该指出，不同阶段划分的阈值=0.31和0.42本身仅仅是在本文的实验条件下获得的，在其它实验条件下可能会有所不同。但是，这里所揭示的高浓度固液两相流中颗粒浓度超过一定阈值后发生的微观运动特性变化却是非常具有启发性的。这一颗粒运动的微观特性很难从一般的连续介质理论得到。从动理论获得的微观运动特性不仅为理解颗粒运动的细节增加了新的内容，而且也为研究下列颗粒运动的宏观特性奠定了理论基础。3颗粒运动的宏观特性3.1颗粒平均速度垂直分布颗粒运动的平均速度是一个十分重要的宏观特性，可根据式(4)的定义和式(16)给出的颗粒速度分布函数表达式求解，计算结果见图5。图2不同垂直高度下的概率密度分布函数p(u)=(1/n)f变化可以看出，颗粒平均速度随颗粒浓度的变化基本表现为：当颗粒浓度没有超过一定阈值前，颗粒平均速度的垂直分布与常规的固液两相流垂直分布变化相似；但是，随着颗粒浓度的继续增加，颗粒平均速度因颗粒相互作用形式的变化而快速增加，突出表现为颗粒向上聚集趋势的增加和最大颗粒平均速度出现位置的上移。二者的界限大致在颗粒浓度=0.42左右，这与颗粒运动微观特性分析中讨论的第二个阈值相近。3.2颗粒脉动速度垂直分布颗粒脉动速度的垂直分布可根据式(6)的定义和式(16)给出的颗粒速度分布函数求解，计算结果见图6。由图可见，在颗粒浓度<0.31的范围内，颗粒脉动速度的垂直分布随颗粒浓度的增加而缓慢变化；当0.31<<0.42时，颗粒脉动速度垂直分布与颗粒浓度变化关系不大；当>0.42时，颗粒脉动速度随颗粒浓度的增加急剧增加，特别是对于靠近水流表面的区域。高浓度固液两相流中颗粒脉动速度发生本质性变化的阈值与前面讨论的颗粒速度分布函数变化对应的阈值具有较好的对应关系。3.3颗粒浓度垂直分布颗粒浓度(或颗粒数)的垂直分布可以很好地反映颗粒运动的状态，是颗粒运动分析时不可缺少且相对容易测量的重要宏观特性。颗粒浓度(或颗粒数)可根据式(2)的定义和式(16)给出的颗粒速度分布函数表达式求解，也可以对式(16)在速度空间上直接积分得到。仍然采用Michalik[17]的测量数据进行计算，可以得到如图7所示的颗粒浓度分布。由图可见，沿垂线上对应于最大颗粒浓度的位置随颗粒平均浓度的增加而显著上升，这有别于一般的低浓度固液两相流。发生显著变化的对应颗粒平均浓度大约为0.42。图3特征速度uc随颗粒浓度和位置的变化图4标准偏差σ和平均颗粒浓度的关系图5颗粒平均速度垂向分布随平均颗粒浓度的变化图6颗粒脉动速度随颗粒平均浓度的变化由于对高浓度固液两相流中微观特性的直接测量在技术上存在很大的困难，因此要对颗粒速度分布函数进行直接验证是不易做到的。鉴于此，我们可以用测量到的颗粒运动宏观特征之一——颗粒浓度垂直分布进行间接验证。由Michalik的测量数据与本文计算结果进行的比较(见图7)说明，测量和计算结果相当吻合。4结论颗粒运动的动理论为研究高浓度固液两相流的微观特性和宏观特性提供了理论基础。一旦获得微观特性的信息，就可相应得知所有传统的颗粒运动宏观特性。本文重点探讨了微观颗粒速度分布函数及其由它推演出的颗粒运动宏观变量(如颗粒平均速度、颗粒脉动速度和颗粒浓度分布等)随颗粒浓度不断增加发生的变化规律。尽管本文计算时采用的实验参数可以变化，但是由此得到的有关高浓度固液两相流的特性很难由传统的理论获得，而且有关规律对研究者们具有启发性。(1)颗粒速度概率密度分布函数的峰值随颗粒平均浓度增加而减小，但是与颗粒速度概率密度分布函数峰值对应的特征速度uc与颗粒平均浓度的变化关系不大，尽管uc随着垂直坐标的增大向上增加。(2)颗粒速度概率密度分布函数中的标准偏差σ主要受颗粒平均浓度的控制。对水流中任意给定的垂向高度上，σ随颗粒平均体积浓度的变化被两个阈值(即下临界阈值和上临界阈值)分为三个阶段。当颗粒平均浓度小于下临界阈值(本文条件下为=0.31)时，σ随颗粒平均浓度的增加而缓慢减少；当颗粒平均浓度大于上临界阈值(本文条件下为=0.42)时，σ随颗粒平均浓度的增加而急剧减少；当颗粒平均浓度介于两个阈值之间时，σ变化不明显。(3)颗粒运动的宏观特性基本上由颗粒运动的微观运动特性决定。图7测量和计算的颗粒浓度分布在固液两相流进入高浓度固液两相流范围且当颗粒平均大于上临界阈值(本文条件下为=0.42)时，颗粒运动的宏观特性变量，如颗粒平均速度、颗粒脉动速度和颗粒浓度垂向分布等变化剧烈，甚至出现反常变化趋势。(4)应用颗粒运动的动理论能够获得高浓度固液两相流中有关颗粒运动的详尽的微观和宏观特性信息，这些是采用传统连续介质理论和依靠修正低浓度固液两相流研究结果无法做到的。参考文献：[1]钱宁，万兆惠，钱意颖.黄河高含沙水流问题[J].清华大学学报，1979，19(2)：27-34.[2]KriegerIM,DoughertyTJ.AmechanismfornonNewtonianflowsinsuspensionsofrigidspheres[J].Trans.theSocietyofRheology,1959，3：137-152.[3]BagnoldRA.ExperimentsonagravityfreedispersionoflargesolidspheresinNewtonianfluidundershear.[J].Proc.,RoyalSoc.London,Ser.A,1954，225：49-63.[4]TakahashiT.Debrisflowsonprismaticopenchannel[J].J.Hydr.Div.,ASCE,1980,106(3):381-396.[5]ChenCL.Generalsolutionforviscoplasticdebrisflow[J].J.Hydr.Engrg.,ASCE,1988,114(3):259-282.[6]TakahashiT.Mechanicalcharacteristicsofdebrisflows[J].J.Hydr.Div.ASCE,1978,104(8):1153-1169.[7]ChenCL.Generalizedviscoplasticmodelingofdebrisflow[J].J.Hydr.Engrg.,ASCE,1988,114(3):237-258.[8]ChenCL.Rheologicalequationinasymptoticregimesofgranularflow[J].J.Engrg.Mech.,ASCE,1998,124(3):301-310.[9]BrufauP,Garcia-NavarroP,ChilardiPetal.1Dmathematicalmodelingofdebrisflow[J].J.Hydr.Res.,IAHR,2000,38(6):435-446.[10]DrewDA.Mathematicalmodelingoftwo-phaseflows[J].A.Rev.FluidMech.1983,15：261-291.[11]倪晋仁，王光谦，张红武.固液两相流基本理论及其最新应用[M].北京：科学出版社，1991.[12]倪晋仁，王光谦.高浓度恒定固液两相流运动机理探析：Ⅰ.理论[J].水利学报，2000，(5)：22-26.[13]倪晋仁，王光谦.高浓度恒定固液两相流运动机理探析：Ⅰ.应用[J].水利学报，2000，(5)：27-32.[14]NiJR,WangGQ，BorthwickAGL.Thekinetictheoryforverticalprofileofparticleconcentrationindiluteanddensesolid-liquidflows[J].J.Hydr.Engrg.,ASCE,2000，126(12):893-903.[15]WangGQ,NIJR.Kinetictheoryforparticleconcentrationdistributionintwo-phaseflows[J].J.Engrg.Mech.,ASCE,1990,116(1):2738-2748.