流体力学多相流自学作业.docx
《流体力学多相流自学作业.docx》由会员分享,可在线阅读,更多相关《流体力学多相流自学作业.docx(11页珍藏版)》请在冰豆网上搜索。
流体力学多相流自学作业
多相流及其应用
1.
欧阳光明(2021.03.07)
2.两相与多相流的定义与分类
在物理学中物质有固、液、气和等离子四态或四相。
单相物质的流动称为单相流,两种混合均匀的气体或液体的流动也属于单相流。
同时存在两种及两种以上相态的物质混合体的流动就是两相或多相流。
在多相流动力学中,所谓的相不仅按物质的状态,而且按化学组成、尺寸和形状等来区分,即不同的化学组成、不同尺寸和不同形状的物质都可能归属不同的相。
在两相流研究中,把物质分为连续介质和离散介质。
气体和液体属于连续介质,也称连续相或流体相。
固体颗粒、液滴和气泡属于离散介质,也称分散相或颗粒相。
流体相和颗粒相组成的流动叫做两相流动。
自然界和工业过程中常见的两相及多相流主要有如下几种,其中以两相流最为普遍。
(1)气液两相流
气体和液体物质混合在一起共同流动称为气液两相流。
它又可以分单组分工质如水—水蒸气的汽液两相流和双组分工质如空气—水气液两相流两类,前者汽、液两相都具有相同的化学成分,后者则是两相各具有不同的化学成分。
单组分的汽液两相流在流动时根据压力和温度的变化会发生相变。
双组分气液两相流则一般在流动中不会发生相变。
自然界中如下雨时的风雨交加,湖面和海面上带雾的上升气流、山区大气中的云遮雾罩。
生活中沸腾的水壶中的循环,啤酒及汽水等夹带着气泡从瓶中注人杯子的流动等都属于气液两相流。
现代工业设备中广泛应用着气液两相流与传热的原理和技术,如锅炉、核反应堆蒸汽发生器等汽化装置,石油、天然气的管道输送,大量传热传质与化学反应工程设备中的各种蒸发器、冷凝器、反应器、蒸馏塔、汽提塔,各式气液混合器、气液分离器和热交换器等,都广泛存在气液两相流与传热现象。
(2)气固两相流
气体和固体颗粒混合在一起共同流动称为气固两相流。
空气中夹带灰粒与尘土、沙漠风沙、飞雪、冰雹,在动力、能源、冶金、建材、粮食加工和化工工业中广泛应用的气力输送、气流千燥、煤粉燃烧、石油的催化裂化、矿物的流态化焙烧、气力浮选、流态化等过程或技术,都是气固两相流的具体实例。
严格地说,固体颗粒没有流动性,不能作流体处理。
但当流体中存在大量固体小粒子流时,如果流体的流动速度足够大,这些固体粒子的特性与普通流体相类似,即可以认为这些固体颗粒为拟流体,在适当的条件下当作流体流动来处理。
在流体力学中,尽管流体分子间有间隙,但人们总是把流体看着是充满整个空间没有间隙的连续介质。
由于两相流动研究的不是单个颗粒的运动特性,而是大量颗粒的统计平均特性。
虽然颗粒的数密度(单位混合物体积中的颗粒数)比单位体积中流体分子数少得多,但当悬浮颗粒较多时,人们仍可设想离散分布于流体中颗粒是充满整个空间而没有间隙的流体。
这就是常用的拟流体假设。
(3)液固两相流
液体和固体颗粒混合在一些共同流动称液固两相流。
自然界和工业中的典刑实例有夹带泥沙奔流的江河海水,动力、化工,采矿、建筑等工业工程中广泛使用的水力输送,矿浆、纸浆、泥浆、胶桨等浆液流动等。
其它像火电厂锅炉的水力除渣管道中的水渣混合物流动,污水处理与排放中的污水管道流动等。
(4)液液两相流
两种互不相溶的液体混合在一起的流动称液液两相流。
油田开采与地面集输、分离、排污中的油水两相流,化工过程中的乳浊液流动、物质提纯和萃取过程中大量的液液混合物流动均是液液两相流的工程实例。
(5)气液液、气液固和液液固多相流
气体、液体和固休颗粒混合在一起的流动称气液固三相流。
气体与两种不能均匀混合、互不相溶的液体混合物在一起的共同流动称为气液液三相流;两种不能均匀混合、互不相溶的液体与固体颗粒混合在一起的共同流动称为液液固三相流。
图1两相流的范围
3.两相流的研究和处理方法
与普通流体动力学类似,研究两相流问题的方法可以分为理论研究和实验研究两方面。
从理论分析方法来看,仍然存在微观和宏观两种观点。
微观分析法就是从分子运动论出发,利用Boltzman(波尔兹曼)方程和统计平均概念及其理论,建立两相流中各相的基本守恒方程。
宏观分析法,就是以连续介质假设为基础,将两相流中各相都视为连续介质流体,根据每一相的质量、动量和能量宏观守恒方程以及相间相互作用,建立两相流的基本方程组,再利用这些两相流基本方程组去研究分析各种具体的两相流问题。
在许多实际问题中,我们所关心的不是单个物质粒子的运动而是大量粒子运动所产生的总结果,也就是所谓的宏观量,如压强、密度、温度、平均流速等。
从宏观观点分析两相流的方法又可以分为3类。
(1)模型法
即假定相互扩散作用是连续进行,其基本观点是:
(1)两相流混合物体中的每一点都同时被两相所占据;
(2)混合物的热力学和输运特性取决于各相的特性和浓度;
(3)各相以自己的质量速度中心移动,相间相互扩散作用反映在模型内。
(2)容积法
假定过程处于平衡状态,可用平衡方程式进行描述,基本方法是:
(1)认为流动是一维的;
(2)对一个有限容积写出质量、动量和能量守衡方程;
(3)守衡方程即可按混合物写出,也可按单独相列出。
(3)平均法
假定过程处于平衡状态,用平均的守衡方程进行描述,类似低通滤波的方法。
上述3种方法的共同点就是不考虑局部的和瞬时的特性,仅考虑相界面上流体微粒集中的相互作用,即宏观动力学。
4.两相与多相流的专用术语与基本持性参数
一般用下角标1和2分别表示两相流中的两种相或组分,对液—气两相系统用下角l和g区分,对流体—固体两相系统用下角标f和s区分。
通常选定组分或相2为分散相或为分层流动中的轻相。
3.1质量流量、质量流速和质量相含率(相分数)
质量流量是指单位时间内流过通道总流通截面积的流体质量,用W表示。
两相流总的质量流量是各相质量流量之和,各相的质量流量分别用W1和W2表示,所以有
W=W1+W2(kg/s)
质量流速是单位流通截面积上的质量流量,用G表示,如流道的总流通截面积为A,各相所占的流通截面积分别为A1和A2,则有
G=W/AG1=W1/AG2=W2/A
各相质量流量与总质量流里之比称为质量相含率或质量相分数。
在气液两相流系统中,气体的质量相含率俗称为质量含气率或干度,用x表示;液体的质量相含率俗称为质量含液率。
气液两相流的质量含气率与质量含液率之和是1,所以有
x=Gg/G=Wg/W
1-x=Gl/G=Wl/W
3.2容积流量、容积流速和容积相含率
容积流量是指单位时间流过通道总流通截面积的流体容积,用Q表示。
两相流总的容积流量同样是各相容积流量Q1和Q2之和,即
Q=Q1+Q2=(W1/ρ1+W2/ρ2)
容积流速是单位流通截面积上的容积流量,又称折算速度,是容积流量除以通道总流通面积A,用J表示。
J=Q/A=(Q1+Q2)/A=J1+J2
J1=Q1/A=W1/(ρ1A)
J2=Q2/A=W2/(ρ2A)
各相的折算速度在两相流中是十分重要并常用的一个术语和参数,它表示两相混合物中任何一相单独流过整个通道截面积时的速度,称为该相的折算速度。
容积相含率是指各相容积流量与总容积流量之比。
在气液两相流系统气相的容积相含率又称为容积含气率,用β表示,液相的容积相含率称为容积含液率。
3.3各相真实流速
各相容积流量除以流动中各相各自所占流通截面积即为各相的真实流速。
用vi表示(i=1,2,或i=l,g,为气液两相流;i=f,s为流体—固相两相流)。
vi=Qi/Ai
3.4真实相含率或截面相含率
某相的流动在任意流通截面上所占通道截面积与总的流通截面积之比称作该相的真实相含率或截面相含率。
对气液两相流,气相的真实相含率又称为截面含气率、真实含气率或空隙率,用α表示。
而液相所占截面积与总流通截面积之比称为截面含液率,用(1-α}表示。
即;
α=Ag/A1-α=Al/A
3.5滑动比、滑移速度、漂移速度和漂移流率
两相流中各相真实速度的比值称为滑动比。
气液两相流的滑动比用S表示,是气相真实速度与液相真实速度之比
S=vg/vl
滑移速度是指两相流各相真实速度的差,用vs表示
vs=vg-vl=Jg/α-Jl/(1-α)
漂移速度是指轻相(如气相)速度与两相混合物平均速度vH之差,用vD表示
vD=vg-vH
两相混合物平均速度指当滑动比S=1时两相混合物的速度。
漂移流率是指滑移速度vs的式两边乘以通分后的分母项,消去分母后的等式,用jD表示,有
jD=(vg-vl)α(1-α)=Jg(1-α)-Jlα
3.6两相混合物的密度与比容
两相流体的密度有两种表示法:
(1)流动密度。
指单位时间内流过截面的两相混合物的质量与容积之比,即
ρ0=W/Q
(2)真实密度。
指流动瞬间任一流动截面上两相流混合物的密度,用ρm表示,定义如下:
ρm=αρg+(1-α)ρm
由截面含气率α与容积含气率β及滑移比S各自的定义可推得:
显然当S=1,即vg=vl,两相间无相对速度时,α=β流动密度才等于真实密度。
比容是密度的倒数,因此,两相混合物的比容为:
流动比容
;
真实比容
5.气液两相流的流型和流型图
4.1垂直上升管中的气液两相流流型及其流型图
气液两相流在垂直管中上升流动时的几种常见流型。
(1)细泡状流型
细泡状流动是最常见的流型之一。
其特征为在液相中带有散布在液体中的细小气泡。
直径在1mm以下的气泡是球形的。
直径在1mm以上的气泡外形是多种多样的。
(2)气弹状流型
气弹状流型由一系列气弹组成。
气弹端部呈球形而尾部是平的。
在两气弹之间夹有小气泡而气弹与管壁之间存在液膜。
(3)块块流型
当管内气速增大时,气弹发生分裂形成块状流型。
此时大小不一的块状气体在液流中以混乱状态流动。
(4)带纤维的环状流型
在这种流型中,管壁上液膜较厚且含有小气泡。
管子核心部分主要是气体,但在气流中含有由被气体从液膜带走的细小液滴形成的长条纤维。
(5)环状流型
在这种流型中,管壁上有一层液膜,管子核心部分为带有自液膜卷入的细小液滴的气体。
环状流型都发生在较高气体流速时。
在受热管道中,当管壁温度高到足以使管壁液膜汽化时,流动结构就会发展到壁上无液膜,只有气相中还含有细小液滴的雾状流型。
图3表示有单组分气液两相流体在垂直上升受热管中的流型和管壁热流密度的关系。
在图中,温度低于饱和温度的液体以固定流量进入各受热管。
各受热管的热流密度依次自左向右逐渐增加。
由图3可见,随着热流密度的增大,各管中的沸腾点逐渐移向管子进口,各管中的流型也逐渐由单相液体、细泡状流型、气弹状流型、块状流型、环状流型、雾状流型一直发展到干饱和蒸汽和过热蒸汽流动。
在气液两相流中,在两相流量、流体的物性值(密度、粘度、表面张力等)、管道的几何形状,管道尺寸以及热流密度确定的条件下,要判断管内气液两相流的流型可应用流型图。
流型图主要是根据试验资料总结而成的,因而应用流型图时不应超出获得该流型图的试验范围。
在判断垂直上升的流型图中,图4所示的流型图得到较为广泛的应用。
此图适用于空气—水和汽—水两相流,是在管子内径为31.2mm的管子中用压力为0.14-0.54MPa的空气—水混合物为工质得出的,此图和应用压力为3.45-6.9MPa的汽水混合物在管径为12.7mm管子中得出的试验数据相符良好,所以也可适用于上述参数的汽水混合物。
图中横坐标为脚ρLJ2L纵坐标为ρGJ2G可分别按下列两式计算:
5.2垂直下降管中的气液两相流流型及其流型图
在垂直管中气液两相一起往下流动时的流型示于图5。
这些流型是由空气—水混合物的试验结果得出的。
在气液两相做下降流动时的细泡状流型和上升流动时的细泡状流型不同。
前者的细泡集中在管子核心部分,而后者则散布在整个管子截面上。
如液相流量不变而使气相流量增大,则细泡将聚集成气弹。
下降流动时的气弹状流型比上升流动时稳定。
下降流动时的环状流动有几种流型,在气相及液相流量小时,有一层液膜沿管壁下流,核心部分为气相,这称为下降液膜流型。
当液相流量增大,气相将进人液膜,这称为带气泡的下降液膜流型,当气液两相流量都增大时会出现块状流型。
在气相流量较高时能发展为核心部分为雾状流动,壁面有液膜的雾式环状流型。
图6示出的是作下降流动的气液两相流流型图。
该图是以空气和多种液体混合物作试验得出的,试验管径为25.4mm,试验压力为0.17MPa。
图6选用
作横坐标,用
作纵坐标。
,
Y为液相物性系数,
—液相动力粘度Pa·s;
—20℃,0.1MPa时水的动力粘度Pa·s;
—液相密度kg/m3;
—20℃,0.1MPa时水的密度kg/m3;
—液相表面张力N/m;
—20℃,0.1MPa时水的表面张力N/m
5.3水平管中的气液两相流流型及其流型图
气液两相流体在水平管中流动时的流型种类比垂直管中的多。
这主要是由于重力的影响使气液两相有分开流动的倾向造成的。
气液两相流体在水平管中流动时流型大致可分为6种,即:
细泡状流型、气塞状流型、分层流型、波状分层流型、气弹状流型及环状流型。
图7中有这6种流型的示意图。
在细泡状流型中,由于重力的影响,细泡大都位于管子上部。
当气体流量增加时,小气泡合并成气塞形成气塞状流型。
分层流型发生于气液两相的流量均小时甲此时,气液两相分开流动,两相之间存在一平滑的分界面。
当气相流量较高时,两相分界面上出现流动波,形成波状分层流型。
气相流量再增大会形成气弹状流型,但此时气弹偏向管子上部。
当气相流量很高而液相流量较低时就出现环状流型。
水平管中的气液两相流流型也可按相应的流型图确定。
在现有的各种水平管流型图中,Baker流型图建立最早且得到广泛应用,特别在石油工业和冷凝工程设计中。
图8所示为坐标经过改进后的Baker流型图。
横坐标为JLФ,纵坐标为JG/λ。
其中JG、JL分别为气相及液相的折算质量流速;系数λ和Ф值与流体的物性有关。
6.界面现象
5.1概述
两相与多相流中另一个最基本最显著的特征是混合物系统中存在明确的相界面,并且伴随着混合物的流动,这些相界面的状态、形状也处于不断的运动或变化发展中。
相界面是指将两种不同相的物质分隔开的区域,在此区域内的物质的特征和性质不同于相邻的区域。
界面的邻近区域内物质的密度和浓度均有显著的差异,在接近两相分隔的临界点,密度在与界面相垂直的方向上是位置的连续函数。
也就是说界面是一个三维区域,其厚度可能有几个分子直径大小或者更厚一些,这已经为统计力学的计算结果所证实。
所谓界面现象,是指与相界面上的动量、能量和质量传递相关的所有效应。
许多物理、化学作用都发生在表面和界面上,例如沸腾和凝结的状态会因界面上存在即便是极少量的杂质而发生明显的改变,泡沫的稳定性、湖面的波度以及镜体表面的汽水雾罩等现象都是发生在界面上的物理、化学变化。
5.2速度边界条件
如果不考虑相变或传质、速度相容性在单组分和双组分的两相流动中都是相同的。
但在有相变时,就可能在界面上存在一有限的速度跳跃。
流体流经界面必须满足连续性条件,如果密度有变化则速度也会有变化。
囚此,相对于图9所示的相界面,在平行和垂直相界面的两个方向上都必须满足如下的相容性条件:
(1)连续切向速度要求
v1t=v2t
(2)穿透相界面的连续性要求
ρ1v1n=ρ2v2n=m
式中:
m表示垂直通过相界面的流体质量流量。
通过相变过程中的平均连续速度被假设为垂直于相界面的。
如果界面还正在移动、就将界面的移动速度简单地叠加到上述速度上去。
5.3应力边界条件
表面张力是界面上每单位面积的自由能,即形成一单位表面积所需外界做的功,因为将物质分子从体相内移动到表面区是要克服分子间的吸力而做功的。
也可以把表面张力想象为表面作用在每单位长度截线上的力。
界面稳定存在的必要条件是表面张力大于0,如果小于0,那么偶然的涨落就可导致表面区不断扩大,最后使一种物质完全分散到另一种物质中,两种气体之间混合的情况就是这样。
通过界面的应力场的连续性因表面张力效应而改变。
如果表面张力是均匀的,则通过界面的剪切应力也是连续的。
但杂质、灰尘或表面活性剂的存在、或则沿界面的温度梯度等都会引起表面张力梯度的形成。
这种情况下,通过界面的剪切应力存在一等于表面张力梯度的突增量。
7.两相流基本数理模型
两相流中各相在空间和时间上随机扩散,同时存在动态的相互作用。
对于这种复杂的二维两相瞬态问题,完整的解析解目前还无法导出。
人们先后提出了多种数理模型,从最简单的均相流模型一直到最复杂的双流体模型。
6.1均相流动模型(HomogeneousFlowModel)
均相流动模型把气液两相混合物看作是一种均匀介质,相间没有相对速度,流动参数取两相相应参数的平均值。
在此基础土,可将两相流视为具有平均流体特性的单相流对待。
在均相模型中采用了两个基本假设:
(1)两相间处于热力学平衡状态,即两相具有相同的温度并且都处于饱和状态;
(2)气液两相的流速相等,即为均匀流。
其数学表达式及守恒方程如下:
uG=uL=uH
S=uG/uL=1
在通道微元段
上的动量守恒方程为:
动量增加率=动量流出率-动量流入率+动量存积率=作用于控制体的力之和
Pp为通道圆周
对于等截面通道的稳定流动,这个方程简化为
能量方程并不直接用于压降的计算,主要用来计算绝热流动中的局部干度,能量守衡方程为:
能量流人率=能量流出率+能量存积率
化简后
e—单位质量流体的对流传热热能
实际上气液两相的流速并不相等,只有在高含气量或很小含气量时两相流速才近似相等,因此这一模型实际上只适用于泡状流和雾状流。
从前人的研究中可以看出,均相模型的计算结果同实验值有较大出人,偏差随质量流速的减小而增大。
这种偏差的存在是因为均相模型假设两相之间没有速度差异。
当质量流速较小时,浮力效应显著,引起两相速度间相当大的差异;而质量流速较大时,液相湍动的结果使得两相的混合更加均匀,因此质量流速
增大时偏差减小。
6.2分相流动模型(SeparatedFlowM0del)
分相流动模型是将气、液两相都当作连续流体分别来处理,并考虑了两相之间的相互作用。
其基本假设是:
(1)两相间保持热力学平衡;
(2)气液两相的速度为常量,但不一定相等。
假定气液两相都以一定的平均速度在流道中流动。
分相流动模型在一定程度上考虑了两相间的相互作用,计算结果比均相模型理想。
当两相平均流动速度相等时,分相流动即可转化为均相模型。
因此可将均相模型视为分相模型的一个特殊情况。
分相流动模型适用于两相间存在微弱耦合的场合,如分层流和环状流。
分相流动的基本方程推导如下:
先写出每一相的各个守恒方程,然后将每一组的守恒方程进行合并,成为混合物的总平衡方程这也是在确定压降时被人们广泛采用的分析方法。
连续性方程,也就是质量守恒方程的数学式为
Ge是单位长度上液体转化为气体的转化率。
分相流动模型的动量守衡方程同样可以从各相的动量方程导出,混合物动量方程
分相流动模型的两相能量守恒方程式为:
6.3漂移模型(Drift-FluxModel)
这一模型主要由Zuber等人提出。
它是在热力学平衡的假设下,建立在两相平均速度场基础上的一种模型。
漂移模型提出了一个漂移速度的概念,当两相流以某一混合速度流动时,气相相对于这个混合速度有一个漂移速度,液体则有一个反向的漂移速度以保持流动的连续性。
在守恒方程组中将相间相对速度以漂移速度来考虑,通过附加的气相连续方程来描写气液两相流动。
6.4两流体模型(TwoFluidModel)
由于各相的动力学性质不完全相同或浓度分布不均匀,气、液两相的运动存在相当大的差异。
上述几种模型中,均相模型完全没有考虑两相的差异。
分相模型和漂移模型在一定程度上引入了气液两相的相互作用,但仍过于简单而无法精确描述两相的运动与空间分布。
目前一致公认最为完善可靠的模型是两流体模型。
它可以用欧拉方法、拉格朗日方法来描述,每种方法都有其固有的优缺点。
工程实际中采用的方法一般基于积分法或时均紊流模型,往往不能提供瞬态的流动结构。
而这对于离散相的输运恰恰又是非常关键的。
只有拉格朗日方法能够考虑离散相和瞬态流动结构之间复杂的相互作用,但其表达式尚存在某种程度的不确定性,需要很多理想化的假设。
两流体模型将每一种流体都看作是充满整个流场的连续介质,针对两相分别写出质量、动量和能量守恒方程,通过相界面间的相互作用(动量、能量和质量的交换)将两组方程祸合在一起。
这种方法只需假设每相在局部范围内都是连续介质,不必引人其它人为假设,而且对两相流的种类和流型没有任何限制,适用于可当做连续介质研究的任何二元混合物,所建立的两相流方程是目前最全面完整的,求得的解中包含的信息丰富完全。
但两流体模型包含的变量多、方程复杂,因此求解比较困难。
8.多相流应用
7.1高效旋风式分离器
其办法是把圆筒壁面做成粗糙面,使湍流增强,碰撞壁面机会增多。
旋风分离器可用于浮渣分离、豆奶脱臭、粉尘分离、喷漆分离等有强刺激气体的场合。
7.2高浓度多相流泵
常用的泵多按单相流理论设计。
对于高浓度多相流,由于流体的密度、粘度都发生了很大的变化,且随着各相浓度的不同,固相或气相均有不同变化即使旋转速度不变化。
即由于粘度变化,雷诺数相似律不存在了因此。
高粘流体泵的设计应按流体的流变参数考虑。
我国蔡保元教授设汁的两相流泵只考虑固、液两相的相对运动,把固相作为动的边界条件,山于两者的惯性力不同,便会发生“相对阻塞”和“相对抽吸”作用。
在抽浑水时.两相流泵的效率比清水泵高但该两相流泵没考虑粘度的影响,而在高粘度下这种影响是不可忽视的。
7.3三相流垂直提升一泵
为了开采海底矿产,国外已应用三相流垂直提升技术采集锰结核,可把5000米深的海底的锰块提升上来。
这一技术也可用于深井开采、井底清淤、河道疏浚、污染底泥处理等方面。
参考文献
[1]郭烈锦.两相与多相流动力学[M].西安:
西安交通大学出版社.