《两相流动的基本方程》由会员分享，可在线阅读，更多相关《两相流动的基本方程（67页珍藏版）》请在金锄头文库上搜索。
1、研究生课程：研究生课程：气液两相流理论气液两相流理论两相流动的基本方程两相流动的基本方程中国石油大学（北京）石油天然气工程学院中国石油大学（北京）石油天然气工程学院 檀朝东檀朝东 138013312551.1均相模型均相模型1.2分相模型分相模型2.解析模型基本方程解析模型基本方程1.两相流动的基本方程两相流动的基本方程2.1二流体模型二流体模型2.2混合模型混合模型2.3滑移模型、扩散模型、漂移模型滑移模型、扩散模型、漂移模型1.11.1均相模型均相模型均相模型把两相介质看成均匀介质，介 质的参数取两相平均参数，然后根据单相均 匀介质建立两相流基本方程。下面只讨论两相流均相流模型的动量方程
2、。l单相流的动量方程：摩擦剪切应力是线性的，)(11 rurrrzpgzuuzz z +=0=r R 0zRu r= 1.11.1均相模型均相模型单相流的动量方程可以写成：上式可见，流道压降由三部分组成： 摩擦压降加速压降重力压降021()zurr rrR=02z zupugzRz=+0024Fp zRD=azz zpuuuGzzz=gpgz=1.11.1均相模型均相模型用两相流的混合流速U和混合密度代替单 相方程中的和，这样就可得到两相流 均相模型的动量方程形式为：也需要按两相参数处理，其中，Re与 都需取描述两相流的平均数值。 均相模型属于经验模型范畴。cmzu002cmpUGgzRz=+
3、1.21.2分相模型分相模型分相模型是把两相分成两种单相流动（气相 和液相）。介质参数分别取各自的介质参 数，分相模型也是一种经验模型。 1.连续方程 液相入流面： 液相出流面：气相转化为液相的量：lMl ldMMdz+dzlidW C1.21.2分相模型分相模型 则液相连续方程为：为液相单位面积质量流量。 液相液相连续方程可写为：式中：为液相界面效应；为其周界；为空 泡率。(0l lllidMMMdWCdz+=dz)+(1)llllMuAAG=lG(1) +=0llliduA dz dWCdzldWiC1.21.2分相模型分相模型同理，气相连续方程为：其中，为气相界面效应，而且 ； 为潜热（
4、KJ/Kg）; q为热流密度（）+=0gggiduA dz dW Cdz/gldWdWq =2/kWmgdW1.21.2分相模型分相模型气液相连续方程相加，界面效应抵消，可得：分相模型中各相本身的方程存在界面效应，当分相方程合并成两相方程后，界面效应将不再存在，这是分相模型的缺陷。(1) 0llggdduAuAdzdz+=1.21.2分相模型分相模型2.动量方程 液相入流面： 动量：压力： 液相出流面： 动量：压力： 管壁剪切应力： 相间切应力： 重力分量： 忽略了界面质量变化引起的动量项。llM u(1)pA(+llllMdMudu+)(）(1)( (1) )dpApA dzdz+00llC
5、 dziiC dz(1)coslAgdz1.21.2分相模型分相模型根据动量定理动量定理，可以得出：整理上式：00(1) (1)( (1) )(1)cos(+llliilllllldpApApA dzdz AgdzC dzCdzMdMu duM u+=+)(）-00(1)(1) cos1(+llii lllllllCCdpgdzAAMdMu duMuAdz +=+)(）-1.21.2分相模型分相模型同理：气相动量方程为：式中：分别为气相、液相所接触的流 道壁面周界，且，为流道壁面 总周界，为界面周界。 气液相动量方程相加，且，可 得：00cos1(+ggii gggggggCCdpgdzAAM
6、dMuduM uAdz =+)(）-00,glCC000glCCC+=0CiC000lg=1.21.2分相模型分相模型等截面流道，00(1)cos1()glllggCdpgdzA dM uM uA dz +=+(1)lllMAu=gggMAu =2200cos(1)cmllggCdpdguudzAdz =+1.21.2分相模型分相模型则：X-质量气流率/()gguGx =(1)/ (1)lluGx=22 200(1)cos()(1)cm lgCdpdxxgGdzAdz =+1.21.2分相模型分相模型摩擦压降：加速压降：重力压降：可以看出，两相流和单相流是相似的。0000 222FCRdp d
7、zARR =22 2(1)()(1)algdpdxxGdzdz =+cosg tpdpgdz=1.21.2分相模型分相模型3.能量方程： 分相模型的能量方程，不考虑界面效应，只考虑 每一相的作用。 热力学第一定律：由热力学第一定律可得：dUQW=+UQW=+2(/2)sinq Wdhd ugdz=+其中，h为混合物比焓；W为对外做的功；q为外界对体系的传热；1.21.2分相模型分相模型e为比热学能，F为体系耗散功。体系不对外做功 ，w=0, 而可得：对z求导得：+dhde dpvdF=+qdepdv=+2(/2)sin0vdp dFd ugdz+=21(/2)sin0dpdFdugdzdzdz
8、+=2(/2)sin0dp dFdugdzdzdz+=1.21.2分相模型分相模型利用质量含气率x对动能加权：代入分相模型两相流能量方程得：2222233 2 2222(1)(1)()(1)(1)(1) (1)glglgluxux uGxGxxxxxG =+=+=+23322221(1)sin2(1)gldpdFGdxxgdzdzdz =+1.21.2分相模型分相模型能量方程包含三个部分，等号右侧第一项为能量耗散部分（并不完全等于前面动量方程 中摩擦压降所造成的能量损失）；第二项为介质动能部分；第三项为介质位能部分。整个过程中，密度不变（x不变或等温流动），则 能量方程变为：2332222(1
9、)sin2(1)gldpdFGdxxgdzdzdz =+1.31.3总结总结由以上推导可见，不管是均相模型还是混合后的分相模型，都没有反映出界面效应。它们不能用来研究流场中的局部特性，而只能研究流道的整体特性，或流道的对外效果。在工程中，有许多程序采用这种经验模型。使用分相模型时，对两相流摩擦阻力系数和空泡份额还需建立结构式。2. 2.解析模型基本方程解析模型基本方程2.12.1二流体模型二流体模型单相流体连续性微分方程：设k代表相的角码（k=g,l）,相的数目再多 时，k可取1，2，3，可以写出每相的连续 性微分方程为：()()()+0yxzuuu txyz+=()0ut+=()0k kku
10、t+=2.12.1二流体模型二流体模型l纳维-斯托克斯方程式（动量方程） 三维坐标下，单相流体x,y,z方向的动量微分方 程为：其中，X，Y，Z,分别为x,y,z轴所受的质量力。1()yxxxzxxpduXxyzdt +=1()xyyyzyypduYxyzdt+=1()yzxzzzzpduZxyzdt +=2.12.1二流体模型二流体模型质量力只有重力时，N-S方程的矢量式为()dupITgdt=+zdzduudyduudxduutuydzduudyduudxduutuxdzduudyduudxduutu dtduz zz yz xzy zy yy xyx zx yx xxvvv)()()(+
11、=duuuudtt=+（ ）2.12.1二流体模型二流体模型N-S动量方程可变为：则每一相的动量方程为：式中：为相速度向量；为相压力标量 I 为单位张量 ；为剪应力张量； 为重力加速度向量。=()uuupITgt+（ ）=()k kkkkkkkkuuup ITgt+（ ）kukpkgkT2.12.1二流体模型二流体模型3.能量方程式中：为比热力学能标量。22 (+)(+)22 ()kkk kkkkkkkkkkkuueeuqt p ITugu+= + +ke2.12.1二流体模型二流体模型在流场中，除了k相以外，还存在相间界面区， 如果认为每一相区中，介质特性和单相流一样是连续 的，那么界面区则
12、对每一相都不连续。为了描述整个流场特性，建立界面特性基本方 程，与相特性相同，界面特性是以界面区域均一特性 为基础的。为界面区内两相的界面；为的单位向量（法向量）；12aa，12,nn12aa，2.12.1二流体模型二流体模型为界面总厚度之和，通常 ；N为封闭端面的单位向量；为封闭端面；为 与的交线。界面控制体为与所包围的体积；如果 认为很小，则，对控制体的积 分特性为：12=+12=iiiiaiaiiV12nn= iV1221()()ikiiiikkikkk kVaiiiii VddVnuuT dadtNuuT d ddV =+ 2.12.1二流体模型二流体模型式中与为流场特性参数和流场体积
13、源。界面特性 方程也是三个（质量方程，动量方程，能量方程）。l对质量方程l对动量方程l对能量方程按不同情况决定。 上始终T为针对流场特性的流出率（ef lux）; 为相速度向量。把上式经过体积分和面积分的变换 （雷诺转换和格林变换）以及必要的简化，最后可得 出三个基本方程。1,0= ,ug=2,e u= +ku2.12.1二流体模型二流体模型假设则质量守恒方程中， 则质量方程对V的积分均为零。 则：进行一维处理，可以得出界面质量方程为：0 =1,0,0iikT=21()0kkkikkk kanuuT da =+=()()ggigllilglnuunuunn= =，()()gigliluuuu=
14、2.12.1二流体模型二流体模型动量守恒方程式同样进行处理，其中， ,可以得出界面动量方程为:能量方程式中， 得出界面能量方程式为：,0iiikug T=()()gigggglillluu upuup+=+22,ikkikkkkkkkuugu Tqp uu=+=+22()(2)()(2)ggigggggggllillllllluuUuuuquuUuuuq+= +2.12.1二流体模型二流体模型这样，就可得出三个界面方程，前面气液两相得出六个基本方程（连续，动量，能量方程），九个方程可以耦合成三组描述两相流场的基本方程组。这三组方程是单独描述气相、液相和界面的，它们是不连续的，不能直接求解，要达到求解的目的，必须把它们连续化连续化，使其成为类似单相流体的连续方程组。2.12.1二流体模型二流体模型所谓连续化，就是用平均的方法，把方程改造成一定的平均尺度内既保留原流场（不连续）的起伏变化特性，又表现出连续性。如右图，某个特性参数（或数组方程，包括代数方程或微分方程）在原流场中是不连续的，如折线A所示。经过平均后，便成B。这里的平均，可以按时间间隔进行，称为时间平均。也可以按空间（
