压降计算中摩阻系数的确定.docxVIP

压降计算中摩阻系数的确定 多相流技术的计算主要包括流型分类、固液体积和压力降。流型划分和持液率计算是基础, 压降计算是最终目的。工艺计算中, 压降计算是确定多相管流管径、泵和终端分离器大小的基础, 因此准确计算沿线压降具有重要的意义。沿程摩阻系数是压降计算的重要内容之一, 其计算的准确性直接决定了压降计算的准确性。因此有必要对多相管流的沿程摩阻系数进行分析。 多相管流的沿程摩阻系数的确定比单相流体复杂得多。主要是因为多相管流中不仅存在气体或液体与管壁之间的相互作用, 气液界面之间还存在质量交换和能量交换。目前多相管流的沿程摩阻系数还只能采用以实验为基础的经验或半经验关系式来计算。流型不同, 沿程摩阻系数计算方法就不同。因此, 在计算摩阻系数之前, 需确定多相管流的流型。1990年Xiao等人的流型判别法为目前准确性较好、适用性广的流型判别法。Xiao等人将流型分为分层流 (分层光滑流和分层波状流) 、段塞流、环状流和分散泡状流4种流型。 1. 气体或液体与管壁之间的相互作用 在多相管流中, 流型不同, 沿程摩阻系数计算的方法就不同。有的流型仅需用到气体或液体与管壁之间的相互作用, 比如分散泡状流流型;也有的流型不仅用到气体或液体与管壁之间的相互作用, 还要用到气液界面之间的相互作用, 如分层流、环状流和段塞流。气体或液体与管壁之间的相互作用表现为流体与管壁之间的剪切作用, 同样气液之间的相互作用也表现为气液界面之间的剪切作用。目前, 大多数多相管流研究者认为:在多相管流中, 气体或液体与管壁之间的沿程摩阻系数的计算可以采用单相流体的沿程摩阻系数的计算方法, 而气液界面之间的沿程摩阻系数却只能采用多相管流中的经验或半经验关系式。 1.1 计算模式 在分层流、环状流和段塞流中, 气体或液体与管壁之间的沿程摩阻系数可完全采用单相流体的计算方法。在文献中, 将单相流体的流动区域分为层流区Ⅰ、层流与紊流之间的过渡区Ⅱ、紊流光滑区Ⅲ和紊流混合摩擦区Ⅳ及紊流粗糙区Ⅴ。 在Ⅰ区的沿程摩阻系数λ计算如下: 在Ⅱ区范围小且流动复杂, 目前尚无计算方法, 在精度要求不高的情况下, 可采用Ⅰ区计算方法, 也可以取Ⅰ和Ⅲ区计算的平均值。 在Ⅲ区, 可采用Miller公式来计算: 在Ⅳ区, 可采用Chen氏表达式来计算: 在Ⅴ区, 可将隐式ColebrookWhite方程化为显式: 在Ⅲ区和Ⅳ区的分界判断条件为 在Ⅳ区和Ⅴ区的分界判断条件为 上述各式中,ε为相对粗糙度,λ为沿程摩阻系数。 由于分层流、环状流和段塞流中, 液体和气体没有完全占据管道的整个截面, 因此在分别计算液体和气体的沿程摩阻系数λl和λg时, 相应的雷诺数Rel和Reg计算如下: 式中:λl\,λg——液体和气体对应的沿程摩阻系数; wsl\,wsg——液体和气体对应的折算速度, m/s; Dl\,Dg——液体和气体对应的水力直径, m; ρl\,ρg——液体和气体对应的密度, kg/m3; μl\,μg——液体和气体对应的动力粘度, Pa·s。 Dl和Dg的计算如下: 由1976年TaitelDukler模型 其中 式中:X2——Lockhart-Martilli参数, 无因次; Y——Chisholm参数, 无因次; Sg\,Sl——气体湿周、液体湿周, m; Si——管路横截面上气液分界线长度, m; β——管线倾角; hl——液位高度, m。 由于式 (8) 为隐式方程, 可通过 二分法进行迭代, 求解液位高度hˉhˉl。求解时在某些情况下会出现多重根, 通常取最小值为实际值, 然后利用式 (16) 计算Dl和Dg;计算出Rel\,Reg和Dl\,Dg后, 代入式 (1) ～ (6) , 就可求出λl和λg。 1.2 流体速度和雷诺数 当多相管流流型为分散泡状流时, 由于气液两相之间没有滑动, 与管壁接触的流体只有液体, 所以用气液混合物的平均雷诺数Rem和D代入式 (1) ～式 (6) 中, 求出λl。Rem的计算如下: 式中Hl为持液率, H1=wsl/ (wsl+wsg) (18) 2. 气液界面的摩阻系数 2.1 沿程摩阻系数 Xiao等人将Andritsos和Baker方法结合起来, 提出一种计算λi的新方法。如果D≤0.127m时, 采用Andritsos和Baker等关系式计算λi, 即 wsgt为分层光滑流转变为波状流的气相折算速度; 式中:P——工作压力, Pa; hl/D通过式 (3) 迭代求解。 如果D0.127m时, 先采用Baker等提出的方法计算εi: 其中Nwe=ρgw2slsl2εi/σ;Nu=μ2ll2(ρgεiσ) (24) 式中:εi——气液界面的相对粗糙度, 其值在ε～0.25 (hl/D) 之间; Nwe\,Nu—