气液分离器
一、概述
管柱式气液旋流分离器是一种带有倾斜切向入口及气体、液体出口的垂直管。它依靠旋流离心力实现气、液两相分离,与传统的重力式分离器相比,具有结构紧凑、重量轻、投资节省成本等优点,是代替传统容积式分离器的新型分离装置。在气液两相旋流分析的基础上,建立了预测分离性能的机理模型,该模型包括了入口分离模型、旋涡模型、气泡及液滴轨迹模型;依据机理模型,提出了管柱式旋流分离器工艺设计技术指标和工艺步骤.设计根据管柱式旋流分离器的机理模型以及设计工况,完成了管柱式旋流分离器的结构设计、强度分析、理论校核、焊接工艺设计以及分离器内气液两相流的数值模拟,为工程设计和理论设计提供一定的理论依据。
2、旋流式分离器的结构及工作原理
旋流分离器,是一种利用离心沉降原理将非均相混合物中具有不同密度的相分离的机械分离设备。旋流分离器的基本构造为一个分离腔、一到两个入口和两个出口。分离腔主要有圆柱形、圆锥形、柱-锥形三种基本形式。入口有单入口和多入口几种,但在实践中,一般只有单入口和双入口两种。就入口与分离腔的连接形式来分,入口又有切向入口和渐开线入口两种。出口一般为两个,而且多为轴向出口,分布在旋流分离器的两端。靠近进料端的为溢流口,远离进料端的为底流口。在具有密度差的混合物以一定的方式及速度从入口进入旋流分离器后,在离心力场的作用下,密度大的相被甩向四周,并顺着壁面向下运动,作为底流排出;密度小的相向中间迁移,并向上运动,最后作为溢流排出,分离示意图如图1。这样就达到了分离的目的。旋流分离技术可用于液液分离、气液分离、固液分离、气固分离等。本文设计的旋流分离器用于石油钻井中钻井液的气液分离。
3、旋流式分离器的优缺点
在石油化工中装置中,有各种各样的分离器,其中以立式重力气液分离器最为常见,这种气液分离器具有结构简单、操作可靠等持点。立式重力式分离器的主体为一立式圆筒体,多相流一般从该筒体中段进入,顶部为气流出口,底部为液体出口,其结构简图见图2-2。
虽然旋流式气液分离技术在石油化工方面的应用要晚得多,但与常规的重力式分离相比较,它具有很多优点:
① 分离效率高,由于分离原理的不同使得旋流式分离器具有很高的分离效率;
② 成本低,占用空间较小、维护费用少、能耗低、不需要任何帮助分离的介质;
③ 安装灵活方便,旋流器可以任何角度安装;
④ 工作连续、可靠,操作维护方便,一旦设计、调试好, 就可自动、稳定地工作。
旋流式气液分离器有以上优点,但也有如下缺点:
① 由于旋流器内流体的流动产生一定的剪切作用,如果参数设计不当,容易将液滴(油滴或水滴) 打碎乳化而恶化分离过程;
② 通用性较差。不同的分离要求、不同的处理物料的性质往往需要不同结构尺寸或操作条件的旋流器,因此旋流器往往不能互换使用。在欠平衡钻井中,使用旋流式气液分离器分离钻井液中的气体,能充分发挥该离器优点,同时又能有效的避免它的缺点。因此,旋流式气液分离用于分离钻中的气体具有广阔前景。
第三章 管柱式气液旋流分离器
3.1简介
管柱式旋流分离器(Gas-Liquid Cylindrical Cyclone,简写GLCC)是一种新型分离装置(见图3-1)。它既没有可移动部件,也无需内部装置。气液混合物由切向入口进入旋流分离器后形成的旋流产生了比重力高出许多倍的离心力,由于气液相密度不同,所受离心力差别很大,重力、离心力和浮力联合作用将气体和液体分离。液体沿径向被推向外侧,并向下由液体出口排出;而气体则运动到中心,并向上由气体出口排出。这一成本低、重量轻的新型分离器在替代传统容器式分离器方面具有很大的吸引力。在油和气的流量分别为16000m3d和1980Mm3d,分离压力为680kPa的分离工况下,若分别采用管柱式旋流分离器、传统容器型立式和卧式分离器,模拟计算表明管柱式旋流分离器结构尺寸为1.5m×6m,相当于同等规模的传统立式分离器(2.7m×10.5m)的一半左右,相当于传统卧式分离器(5.8m×29m)的四分之一左`右。
图3-1 GLCC结构图
3.2气液分离机理
管柱式旋流气液分离器由入口区(段)、入口分流区、漩涡区、气泡区、液滴区、气相和液相出口配管等部分组成(图3-2)。
3.2.1 入口区
由于管柱式旋流分离器主要依靠旋流产生的离心力实现气液的高效分离,而入口结构决定了分离器的入口气液相分布及气液相切向入口速度的大小。因此入口结构是影响管柱式旋流分离器分离特性的关键因素之一。
气液相流速的不同,油、气两相或油、气、水多相流在入口管和喷嘴内可能呈现分层流、段塞流、分散气泡流或环状流等多种流型。Kouba的实验研究表明,采用向下倾斜的入口管,保证入口管流型呈现分层流将在很大程度上改善气液分离效果、扩展管柱式旋流分离器的适用范围,最佳倾斜角为-27°。而传统分离器采用的入口结构通常为垂直于筒体的结构(目前很多分离器采用的分气包亦为类似结构),实验证明采用垂直结构的管柱式旋流分离器,气液分离效果差,工作范围大约是倾斜向下的入口结构分离器的一半。入口管向下倾斜,在重力作用下有利于形成分层流,实现气液两相的初步分离。
图3-2 GLCC结构图
3.2.2 入口分流区
入口分流区即与入口槽连接的筒体部分(如图3-2)。气液相经入口槽进入入口分流区实现气液的初步分离,上部的气相、下部的液相分别沿筒壁旋转形成旋流场。
3.2.3 漩涡区
经初步分离的液相以切向速度沿筒壁旋转形成强制旋流。除壁面附近,该旋流可看作刚体转动,如图3-3所示,微元控制体以恒定角速度
旋转。流动是轴对称的、各分速度沿z轴变化不大,因此假设微元控制体在轴向、切向上没有加速度,径向加速度
,应用欧拉方程式得径向和轴向运动微分方程:
图3-3 旋流控制体受力分析
(1)
(2)
式中:
为含气泡液相混合物密度;R为筒体半径,
R=d
2。
控制体压力微分方程:
(3)
把方程(1)、(2)代入方程(3),并沿等压漩涡表面积分(
)得:
(4)
整理上式即可得到漩涡方程:
(5)
(6)
利用漩涡方程积分可得入口分流区下方液相总体积,应用液相体积守恒:
(7)
式中:
指平衡液位高度(见图3-2)。
把方程(5)代入方程(7),积分得到常数C,漩涡形状方程如下:
Z=
(8)
从漩涡方程(8)求解漩涡高度
:
=
EMBED Equation.3 (9)
3.2.4 气泡区和液滴区
分离器上部,气相为连续相,液滴分散其中,称液滴区;下部,液相为连续相,气泡分散其中,称气泡区。连续相做涡旋运动,由于气液相密度差,分散相粒子(气泡或液滴)与连续相间存在滑脱。研究分散相粒子的运动轨迹,可以分析分离特性。
1 .气泡轨迹分析
在漩涡区,较大直径的气泡容易被掳获分出,因此气泡轨迹的研究区域是从漩涡底部开始的涡流区。假设液相涡流区小气泡均匀分布,气泡的径向滑脱速度为
EMBED Equation.3 ,轴向滑脱速度为
,忽略连续相-液相的径向速度,轴向速度即液相折算速度
,这样,气泡的径向绝对速度为
=-
,轴向(向下)的绝对速度为
,显然,实现气泡分离的必要条件是壁面处气泡运动至中心气核的时间应少于气泡随连续相—液相流出分离器的时间。
2 .液滴轨迹分析
经过入口分流区初步分离后的旋转气流携带液滴进入分离器上部的液滴区,与气泡区不同的是,由于分散粒子液滴的密度大于连续相气相密度,液滴被甩向器壁。假设径向液滴滑脱速度为
,轴向滑脱速度为
,忽略气相的径向流速,气相的轴向速度即折算速度为
,这样,液滴径向绝对速度为
,轴向(向上)的绝对速度为
,液滴能够分出的必要条件是分离器中心处的液滴到达器壁的时间应短于液滴随旋转气流流出分离器的时间。到达器壁的液滴,在旋转气流的作用下,将在器壁上形成螺旋状薄层液流沿器壁向下流动分出,完成气液分离。
3 .连续相旋流特性
分离器内连续相旋流场分布对气液分离至关重要。连续相切向速度
、径向速度
和轴向速度
构成了旋流场。由于粘性耗散和壁面摩擦阻力的作用,切向速度沿轴向逐渐减小,采用旋流强度表征旋流场的衰减特性。旋流强度
定义为在某一高度截面上连续相切向动量通量与总轴向动量通量的比。
(10)
式中:
是轴向平均流速,下标c表示连续相。
旋流强度沿轴向变化的经验关系式:
(11 )
式中动量通量之比:
(12)
式中:
是连续相(气相或液相)入口质量流量;
是器内对应连续相(气相或液相)的质量流量;
指分离器横截面积;
指对应连续相入口实际流通面积。
根据实验与数值模拟研究,径向流速通常较切向流速和轴向流速低二个数量级,对气液分离效率的影响非常小,因此忽略径向速度。沿旋转半径方向上,切向速度的变化非常显著,按照其不同变化规律,可以将旋流分成两个旋转区域:近壁面处的自由涡旋区和中心强制涡旋区。忽略壁面附近的自由涡旋区,将旋流场视为强制涡旋运动,切向流速近似呈线性分布;将轴向流动视为匀速运动,轴向流速等于平均流速。对于工程应用,以上的假设是趋于保守的。切向速度分布:
(13)
将切向速度表达式代入旋流强度定义式得:
(14)
根据定义Ω是分离器截面位置(z)的函数,确定了某截面的旋流强度,即可求得壁面处(忽略壁面附近的自由涡旋区)的最大切向速度。
(15)
3.3 气泡、液滴轨迹模型
假设:①气泡或液滴分散粒子为球形,在运动中不变形;②粒子间及粒子与器壁间没有相互作用力;③气泡或液滴的运动为定常流动;④作用在气泡或液滴上的各种力的合力为零,分散粒子做匀速运动;⑤分离系统与外界没有热交换,视为等温系统;没有漩涡脱落现象。
径向上,分散粒子受离心力和阻力作用,根据受力平衡,气泡径向滑脱速度:
(16)
同理,液滴径向滑脱速度:
(17)
式中:
、
为气泡和液滴直径;
、
为气液连续相切向速度;
、
为气泡和液滴阻力系数。
阻力系数采用文献提出的关系式求解:
(18)
(19)
式中:
,
轴向上气泡受重力和阻力的作用,但由于连续相为油相,柔占度大,滑脱速度小,其流态一般总处于层流区,因Stokes公式适于求解气泡的滑脱速度:
(20)
同样,液滴在轴向上受重力和阻力作用,根据受力平衡,液滴轴向滑脱速度:
EMBED Equation.3 (21)
气泡或液滴所受阻力与气泡或液滴的合成速度有关,其滑脱速度分别为:
EMBED Equation.3 (23)
(24)
根据上面的气泡和液滴轨迹分析,
t时间间隔中,气泡或液滴在径向和轴的
位移(见图3-4)
图3-4 气泡和液滴轨迹示意图
消去
,
,即为气泡和液滴轨迹控制方程,积分得气泡或液滴的轴向位移,
(25)
(26)
根据气泡或液滴实现分离的必要性,Δ
或Δ
即为气泡或液滴区的最小高度。
第4章 结构设计分析
4.1 入口设计分析
由于管柱式旋流分离器主要依靠旋流产生的离心力实现气液的高效分离,而入口结构决定了分离器的气液分布及其初始切向入口速度的大小,因此入口结构和尺寸是影响管柱式旋流分离器实现气液分离的关键因素。管柱式旋流分离器入口主要由入口管、喷嘴和入口槽3部分组成。
4.1.1入口管分析
气液相流速的不同,油气两相或油气水多相流在入口管内可能呈现分层流、段塞流、分散气泡流或环状流等多种流型。实验研究表明,采用向下倾斜的入口管,保证入口管流型呈现分层流将在很大程度上改善气液分离效果,扩展管柱式旋流分离器的适用范围。而传统分离器采用的入口结构通常为垂直于筒体的结构(目前很多分离器采用的分气包亦为类似结构),采用垂直结构的管柱式旋流分离器实验证明,与倾斜向下的入口结构相比,气液分离效果明显变差,工作范围大约减小一半。入口管向下倾斜,在重力作用下有利于形成分层流,实现气液两相的初步分离,同时,向下的倾斜结构使经过初步分离的液相在入口下方旋转一圈后形成旋流场,避免了对气相向分离器上方运动的阻塞。故入口管采用倾斜入口.
入口管倾角以-27°为宜,管长取1.0~1.5m,入口管直径的选取应保证流型为分层流,由Taitel和Dukler预测模型确定,分层流转变为间歇流或环状流的判别准则为:
(27)
式中:
EMBED Equation.3
式中D为入口管直径,A是它的横截面积,h
是无量纲液位高。迭代求解准则方程,D和h
作为迭代参量,直至准则方程左端小于1,din即为满足分层流条件最小入口管直径。
4.1.2入口喷嘴分析
喷嘴是入口段最后一个影响进入分离器气液相流速分布和入口切向流速大小的因素。通过对3种不同入口开槽结构(矩形、同心圆形及新月形)的初步实验发现,同心圆形喷嘴(缩口管)结构的分离特性最差,而矩形结构喷嘴的分离效果最佳,新月形结构喷嘴的分离效果与矩形喷嘴接近,由于矩形槽结构加工困难,推荐采用新月形结构。
入口喷嘴面积的选取应保证入口液相流速在4.5-6ms之间。液相流速过小将难以发挥旋流离心分离的作用,但液相流速过大将形成过高的漩涡区,在筒体中过早出现气相夹带液滴和液相夹带气泡现象,影响分离效果。
4.1.3双入口分析
双倾斜入口将入口流预分为两股流动:低入口的富液流和高入口的富气流。双入口的试验表明中等大小的气体流量(在入口处段塞流转为分层流)下,气体带液率有明显降低,当气体流量较高时(在入口处为环空流),无多大变化。
4.2主体结构设计分析
4.2.1入口位置
对于没有液位控制的GLCC,将入口段定位于靠近液面的上方是至关重要的。最新的许多试验都表明,单入口GLCC的最佳液面大约在距离入口下方1~3Ld处。过低的液面,如距离入口处远大于3Ld,会导致切向入口速度的过度衰减,影响GLCC的性能。如果液面高于入口,气体会通过液体而溢出,造成更多液体的携带。
4.2.2最佳外形比
外形比是指GLCC的长径比。GLCC的尺寸影响其性能及造价。对于一个给定的直径,GLCC中入口上方的长度提供了液流扰动的容量,而入口下方的长度则决定了用于从液体中分离气泡的存留时间。另外,离心力和浮力的大小与直径成反比,切向速度衰减与长度成正比。由于这一现象的复杂性,最近才刚刚提出了一套决定最佳外形比的基本标准。
4.2.3旋流体锥度
针对反锥型、正锥型和圆柱型旋流体进行的研究表明,对于气液分离,圆柱型旋流体要稍优于反锥型和正锥型结构。
4.3 出口管设计分析
气液相出口管线的配置可根据气液相流量、配置仪表的要求确定。建议气相出口流速取3~30ms,液相出口流速取1.2~12ms。容器内平衡液位应低于入口0 .3m,分离气液相的汇合点一般低于入口点0.3~0.5m,以保持正的静水压头。若分离器配置控制系统,汇合点位置可以高于入口。
第5章 分离器工艺设计
5.1 计算分离器直径
考虑分离器上部的气相分离部分,分离器直径的选取应避免气相中夹带液滴,以气相折算速度表示,即气相折算速度不能大于气流中出现液滴时的临界速度。气相临界速度是:
(28)
式中We是无因次Weber数,它决定于液滴的尺寸,这里取值7。
对于分离器入口以下的液相分离部分,应充分发挥离心分离特性,避免液相中夹带气泡。研究表明,保持液相入口切向流速
和液相流速
的比为40时,旋流分离效果最佳。通常切向流速一般取6ms,显然临界液相流速
=0.15ms。对于高气油比的分离工况,分离器直径:
(29)
对于低气油比的分离工况,分离器直径:
(30)
式中
和
分别是分离工况下气体流量和液体流量。
根据入口喷嘴分析,取液相入口切向速度为:
=6ms
由上分析知当
EMBED Equation.3 =40时, 旋流分离效果最佳,故液相流速:
=0.15ms
已知:日处理液量:900 m3d=0.0104m3s, 日处理气量:90m3d=0.00104 m3s,又知本设计工况是低气油比工况,故由公式(30)得
分离器的直径:
=
2(m)
经圆整后取 d
=300mm
5.2分离器的高度计算
5.2.1确定上半部分的高度
以入口为界分离器分为上下两个部分。经过旋流离心分离,部分液体可能以旋流液膜的形式向上爬升或以液滴的形式随气流向上运动产生气相夹带液滴的现象。上半部分的高度应足够高,甚至在极端流型—段塞流工况下能够吸收液相流量的波动,避免气相夹带液滴的现象。参考传统立式分离器的沉降分离段的处理方法,结合现场应用经验,管柱式旋流分离器的上半部分高度一般取1~1.5m,根据要求分离器的日处理量,设计半部分的高度为:()
{ float v1,v2,v3,r,R,R1,p1,p2,p,s,m,D,n1,n2,n3,n4,n5,cd;
R=0.15,p1=1.185,p2=1500.00,s=0.05, D=0.0005;
scanf("%f,%f",&v1,&r);
printf("\nv1=%f,r=%f",v1,r);
n1=pow(rR,0.9);
p=p1+(p2-p1)*n1;
v3=6*n1;
R1=p*v1*Ds;
n2=pow(R1,0.657);
n3=pow(R1,-1.09);
cd=24*(1+0.173*n2)R1+0.413(1+16300*n3);
v2=sqrt(4*(p-p1)*v3*v3*D(3*p2*r*cd));
m=v1-v2;
if(m==0) printf("\nv1=%f,r=%f",v1,r);
else{
do
{v1=v1-2*m3;
R1=p*v1*Ds;
n4=pow(R1,0.657);
n5=pow(R1,-1.09);
cd=24*(1+0.173*n4)R1+0.413(1+16300*n5);
v2=sqrt(4*(p-p1)*v3*v3*D(3*p2*r*cd));
}
while((fabs(v1-v2))
1e-6);
printf("\nv1=%f,r=%f",v1,r);
}
}
运行结果如下:
r=0.020000, v1=0.000498
r=0.300000,v1=0.001994
r=0.045000,v1=0.005379
r=0.060000,v1=0.009649
r=0.075000,v1=0.016942
r=0.090000,v1=0.026239
r=0.105000,v1=0.036193
r=0.120000,v1=0.049866
r=0.135000,v1=0.063890
r=0.150000,v1=0.084456
应用excel绘出10点的关系曲线如下图5-1所示:
图5-1
根据图示v1与r的曲线关系,v1和r之间近似指数关系,故设
(35)
对上式两边同时取对数得:
(36)
记
,b=lgk 则
r
v
r*r
lgv
rlgv
0.02
0.000498
0.0004
-3.30277
-0.06606
0.03
0.001994
0.0009
-2.7003
-0.08101
0.045
0.005379
0.002025
-2.2693
-0.10212
0.06
0.009649
0.0036
-2.01552
-0.12093
0.075
0.016942
0.005625
-1.77104
-0.13283
0.09
0.026239
0.0081
-1.58105
-0.14229
0.105
0.036193
0.011025
-1.44138
-0.15134
0.12
0.049866
0.0144
-1.3022
-0.15626
0.135
0.06389
0.018225
-1.19457
-0.16127
0.15
0.084456
0.0225
-1.07337
-0.16101
0.83
0.295106
0.0868
-18.6515
-1.27512
0.78
0.292614
0.0855
-12.6484
-1.12805
关于
曲线分析图如下所示:
图5-2
由
的曲线关系知
近似直线关系,故采用最小二乘法确定常数a,b.
由方程组:
(37)
(38)
即
(39)
(40)
把已知代入方程组得:
0.086625a+0.825b=0.284776 (41)
0.825a+10b=12.165652 (42)
解方程组得:a=15.24,b=-3.13
k=0.0007
(43)
气泡在GLCC内穿过的整个轴向距离:
EMBED Equation.3
GLCC各部分计算高度总结如下表:
液滴区(m )
入口分流区(m)
旋涡区(m)
气泡区(m)
总高(m)
1.00
0.12
0.46
1.51
3.09
5.3分离能力计算
① 径向方向上
旋流中气泡受力如图5-3示,由受力可知,
图5-3
相对运动微方程:
(47)
式中:
为气泡质量,
为相对滑移速度,
为气泡直径,
为半径r处的旋流速度。假设旋流为强迫旋涡,即
。则上式可写为:
(48)
解得:
(
为常数) (49)
初始条件:
时,
(50)
设
为时间常数,
,式中第一项
,故
趋于终端沉降速度:
又
,故
(t=0时,
)
若不考虑旋转时的能量损失,则
(
为切向入口速度)
, 解得:
一般地,当
时,近似认为气泡已迁移到中心,对应的时间为最小驻留时间
①
一般气泡从边壁到中心的平均移动速度
,即
②
分离器的处理量为Q,则分离器内液体占据的最小体积
③
式中V为旋流器的容积。
② 竖直方向上
计算在最小驻留时间
内混合体运动的距离时,可忽略气泡与液体的相对滑移,认为气泡随液体一起在重力的作用下向下运动,可得在
内,气泡下降的距离:
④
即旋流器的最小长度为
。
设计参数:
由以上参数和式①
= 2 \* GB3 ②
= 3 \* GB3 ③
= 4 \* GB3 ④计算可得:
旋流器的半径
。当分离器的设计直径
时,相应地
,
5.4平衡计算及液封高度的确定
图5-4 平衡液位示意图
气液分离器如图5-3所示,来自井口的气液混合物进入分离器后产生旋流,在离心力的作用下实现分离,气体占据液滴区并经排气管从气体出口排出,而液体占据气泡区并经液体管从泥浆出口排出。在分离器工作的过程中,可能出现两种极端情况:(1)气量特别大,以至气体压力p1足以克服排液管线的能量的消耗,使气体以连续相从泥浆出口排出;(2)液量特别大,为了克服排液管线的能量消耗,液面上升,有可能出现液相从排气管线排出。在设计设计气液分离器时,根据设计参数的要求,要求避免这两种情况的发生。这就需要研究气液相在输送过程的能量平衡,要求正确的确定液封高度z2和气液混合物的进口高度。
5.4.1 气室压力p1的确定
在管中流动的气体由于温度和压力沿管长变化,其流速和密度也会有显著的变化.根据流体力学理论稳态流动的能量方程有:
+
+
EMBED Equation.3
EMBED Equation.3 d
=0 (51)
通用气体定律给出:
pV=
RT或
=
(52)
设计参数所对应的压力为p
、温度T
、体积V
,那么:
=
将上式带入(2),得
pV=
T
故管内流速
等于
=
EMBED Equation.3
EMBED Equation.3 (53)
将式(3)和式(2)代入式(1),有
pdp – (
)
EMBED Equation.3 dp +
EMBED Equation.3 (
)
Td
=0
假设流动为等温的,积分上式得:
p1
- 2
(
)
EMBED Equation.3 ln(
)=p
EMBED Equation.3 +
EMBED Equation.3 (
)
EMBED Equation.3 (54)
根据Weymouth提出的公式计算,即
=
(55)
上列各式中:
---------------------------------泥浆密度
---------------------------------气室内绝对压力
----------------------------------排气管出口绝对压力
----------------------------------排气管内径
-----------------------------------排气管长度
M--------------------------摩尔质量
m--------------------------气体质量
R--------------------------气体常数
T--------------------------华氏温度
------------------------------------出气管路沿程阻力系数
设排气管的直径DG=20mm,此时
,则
=
EMBED Equation.3
已知日处理气量为:
90m
d=0.00104 m
s,
,
,M=29,
,R=8.315代入(55)式得
(56)
编程计算程序如下:
#include
0.选择第一断面为分离器内液面,第二断面为泥浆出口处,建立两断面的能量平衡关系如下:
+Z1+
=
+Z2+
(1+
+
) (57)
式中:
p
-----2断面压力 Z1-----1断面处液柱高度
Z2-----2断面处液柱高度
EMBED Equation.3 -----1断面处泥浆流速
EMBED Equation.3 -----2断面处泥浆流速 D
-----排液管内径
L
-----排液管长度
-----出液管路沿程阻力系数
由(57)式变形得:Z1-Z2=
+
(1+
+
) -
(58)
(1) 当DL=40mm时,
R
代入(58)式得
,显然不满足条件。
(2) 当DL=80mm时,
R
代入(58)式得
,为了保证液体的顺利排除必须保持一定的压头,同时排液管排出的液体相当于自由流,因此要尽可能的使Z1-Z2
,故所取DL过大。
(3) 当DL=65时,
R
代入(58)式得
,
(4) 当DL=50mm时,
R
代入(58)式得
,由因实际的工况,如果选择次要求的高度太高。
故综上所述,当DL=65mm时比较合理的。
故取Z1=1 .9m,Z2=1.00m
5.5入口管的设计计算
1.入口管倾角取-27°,管长取1m,入口直径的选取是为了保证入口管内呈现分层流,入口管直径根据Taitel&Duker模型确定。
由分流层转变为断续流或环状流的判别式:
(59)
式中:
EMBED Equation.3
式中D为入口管直径,
为管内横截面上液面离管底的高度,A
为管内气相所占的横截面积。D和
作为迭代量,直至准则方程左端小于1,D即为满足分层流条件最小入口直径。
对判别式变形得:
(60)
已知
代入上式得:
(61)
用C语言编程求解最小进口管直径,计算程序如下
#include0.04&&d<0.3&&n<1)
printf("\nn=%f,d=%f",n,d);
}
}
运行结果:
n=0.710000 ,d=0.042261
n=0.720000,d=0.044783
n=0.730000,d=0.047534
n=0.740000,d=0.050547
n=0.750000,d= 0.053856
n=0.760000,d=0.057856
n=0.770000,d=0.061541
n=0.780000,d=0.066028
n=0.790000,d=0.071038
n=0.800000,d=0.076659
n=0.810000,d=0.083000
n=0.820000,d=0.090196
n=0.830000,d=0.098418
n=0.839999,d=0.107880
n=0.849999,d=0.118860
n=0.859999,d=0.131719
n=0.869999,d=0.146938
n=0.879999,d=0.165166
n=0.889999,d=0.187308
n=0.899999,d=0.214625
n=0.909999,d=0.248998
n=0.919999,d=0.293245
故计算求得最小直径: D=40mm
1. 入口槽选型,根据分离器工艺设计分析入口槽选用新月形,同时入口喷嘴的截面积的确定保证入口液相流速为6 ms.
2. 综合考虑本设计分离器的结构,最后确定为单进口。
本设计中GLCC工艺尺寸综合参数如下表:
分离工况
GLCC尺寸
入口管
GLCC回路
(m3s)
(m3s)
气液比
工作压力(MPa)
直径
(mm)
上部高度(mm)
下部高度(mm)
直径
(mm)
排液管直径(mm)
排气管直径(mm)
0.0104
0.00104
1:10
0.8
300
1000
2.3
40
65
20
第六章 强度设计
6.1筒体设计
1. 设计参数如下表:
设计压力(MPa)
最大工作压力
(MPa)
工作温度
容器尺寸(mm)
工作介质
0.8
1.2
常温
(立式)
有毒、易燃
2.筒体厚度计算
首先根据单层圆筒厚度计算公式:
(62)
式中
计算厚度,mm
设计直径
焊接接头系数
计算压力,MPa
设计温度下许用应力
由设计参数,筒体的材料是:16MnR.
根据文献[26]查表D1知,在设计温度下16MnR.的许用应力,当厚度为6-16mm时,[
]
=170MPa,当厚度为16-36时,
EMBED Equation.3 =163MPa
因为容器所盛介质易燃,所以压力容器所有的焊缝采用全焊透结构,根据文献[26]查表4-3知
设计压力:
假设计算厚度为6-16mm,许用应力[
]
=170MPa,则
=
设计厚度:
对于16Mn钢板偏差
因而可以取名义厚度
但是对低合金钢制容器,规定不包括腐蚀余量的最小厚度
不小于3mm,若加上2mm的腐蚀余量名义厚度至少取5mm.又根据钢材的规格取名义厚度
应力判别式:
(63)
式中
计算压力,MPa
设计直径
设计温度下许用应力
--------有效厚度
已知
则代入上式得
故满足强度要求
6.2封头设计
封头采用标准椭圆形封头.材料:16MnR,当厚度为6~16mm时,
,当厚度为16~36mm时,
1. 椭圆形状系数
:
(64)
式中:
------------------------------------椭圆长轴
------------------------------------直边短轴
而对于标准椭圆形封头
,故
=
EMBED Equation.3
2. 计算椭圆的厚度
:
(65)
式中: k----------------------------------椭圆形状系数
pc---------------------------------设计压力
Di---------------------------------------------- ---圆筒直径
------------------------------设计温度下的许用应力
--------------------------------焊接接头系数
假设计算壁厚为6~16mm,则
=170MPa,
于是有
EMBED Equation.3
显然壁厚在假设范围内,则设计壁厚
对于16MnR,钢板负偏差
,故取名义厚度
.但对于低合金钢制压力容器规定不包括腐蚀余量的最小壁厚不小于3mm,加上2mm的腐蚀余量,名义厚度至少应取5mm.又由钢材标准规格,名义厚度应取为 6mm.
3.强度校核
根据应力强度判别式
:
(66)
式中: pc---------------------------------设计压力
---------------------------------------------- ---圆筒直径
------------------------------ 有效壁厚
------------------------------设计温度下的许用应力
把已知条件代入得:
故强度满足要求
4.计算椭圆封头的最大许用工作应力
:
式中:
------------------------------设计温度下的许用应力
--------------------------------焊接接头系数,对封头要求100%的全焊缝
=1.00
------------------------------有效壁厚
Di---------------------------------------------- ---圆筒直径
k----------------------------------椭圆形状系数
代入已知量得:
=
故封头满足最大工作压力的要求
6.3法兰的校核计算
6.3.1排气管法兰校核计算
1.法兰型式及密封面选择
考虑到设计压力较高,介质不允许泄漏,拟选用平焊法兰,凸凹密封面。
法兰材料:16Mn(锻件),
垫片材料:石棉橡胶板,s=3mm,y=11MPa,m=2.00
初拟法兰的尺寸如下所示
D
D=100 D
=70 D
=50
H=10 R=9
垫片尺寸:
垫片宽度:
N=
EMBED Equation.3
垫片基本密封宽度:
垫片有效密封宽度:
故垫片的计算宽度:
2.螺栓的计算
螺栓选用的材料是35号钢,
预紧时螺栓载荷:
操作时螺栓载荷:
EMBED Equation.3
拟采用螺栓4个 ,则
,
根据螺纹标准以及管道法兰螺栓的最小直径,取螺纹根径
,相当于M10的螺栓。
实际螺栓总面积:
螺栓平均间距:
s=
根据文献[28]查表5-10,对于M12-的螺栓,最小螺栓间距为
允许最大螺栓间距:
故
,即所选螺柱直径符合安装和密封要求
3.法兰强度计算
预紧时螺柱载荷:W=
M
=
.m
操作时
垫片载荷:、
压力载荷:
介质静压轴向载荷:
作用在法兰端面上的总力矩:
EMBED Equation.3
法兰应力及度校核
法兰形状系数:
根据文献[28]查表(5-28),(5-29),(5-30),(5-31)得:
f=1.00 F=0.852 V=0.348 T=1.00 U=1.9 Y=1.8 Z=1.32
于是
EMBED Equation.3
轴向应力:
=0.053MPa
径向应力:
周向应力:
MPa
拟定法兰尺寸及选材合适,可安全使用
因为管道法兰的内径符合工称直径系列,应选用标准管道法兰。根据参考文献[24]管道法兰标准选用公称压力为1.6MPa带颈平焊法兰,材料:16Mn,常温下允许的工作压力为1.6MPa.管道法兰标记:管法兰MFM1.6-22 SH3406-92.
6.3.2排液管管法兰校核计算
考虑到设计压力较高,介质不允许泄漏,拟选用平焊法兰,凸凹密封面。
法兰材料:16Mn(锻件),
垫片材料:石棉橡胶板,s=3mm,y=11MPa,m=2.00
初拟法兰的尺寸如下所示
D
D=150 D
=150 D
=80
H=10 R=17.5
垫片尺寸:
垫片宽度:
N=
垫片基本密封宽度:
垫片有效密封宽度:
故垫片的计算宽度:
2.螺栓的计算
螺栓选用的材料是35号钢,
预紧时螺栓载荷:
操作时螺栓载荷:
EMBED Equation.3
拟采用螺栓4个 ,则
,
根据螺纹标准以及管道法兰螺栓的最小直径,取螺纹根径
,相当于M16的螺栓。
实际螺栓总面积:
螺栓平均间距:
s=
根据文献[28]查表5-10,对于M16的螺栓,最小螺栓间距为
允许最大螺栓间距:
故
,即所选螺柱直径符合安装和密封要求
3.管道法兰强度计算
预紧时螺柱载荷:W=
M
=
.
m
操作时
垫片载荷:、
压力载荷:
介质静压轴向载荷:
作用在法兰端面上的总力矩:
EMBED Equation.3
法兰应力及度校核
法兰形状系数:
根据文献[28]查表(5-28),(5-29),(5-30),(5-31)得:
f=1.00 F=0.84 V=0.32 T=1.32 U=2.5 Y=2.32 Z=1.40
于是
EMBED Equation.3
轴向应力:
径向应力:
周向应力:
MPa
拟定法兰尺寸及选材合适,可安全使用
因为管道法兰的内径符合工称直径系列,应选用标准管道法兰。根据参考文献[24]管道法兰标准选用公称压力为1.6MPa带颈平焊法兰,材料:16Mn,常温下允许的工作压力为1.6MPa.管道法兰标记:管法兰MFM1.6-49.5 SH3406-92.
6.4计算设备重量载荷
(1) 设备操作重量:
(67)
----------------------容器壳体和支座的重量
----------------------容器内部构件的重量
----------------------容器保温材料的重量
----------------------平台、扶梯的重量
-----------------------操作时容器内物料的重量
-----------------------人孔、接管、法兰等附件的重量
=243.02+0+0+111+243.7+73.06=670.78kgf
(2)计算分离器各段载荷,所取截面如图1所示
图6-1
各段载荷计算如下表所示
各段载荷
0-1
1-2
2-3
(kgf)
69.6
379.8
221.38
6.5计算风载荷和风弯矩
(1) 计算风载荷
:
(68)
以1-2段为例计算风载荷:
----------------------设备所在地的基本风载荷,取乌鲁木齐的基本风载荷进行计算
=60
---------------------高度变化系数,由[2]表16-6得
=0.78
----------------------1-2段分离器高度,
=2m
-----------------------空气动力系数,对圆筒形设备
=0.7
-----------------------风振系数,
-----------------------计算段顶截面距地面高度系数, 由[2]表16-4得
=0.35
-------------------------设备基本振型自振周期变化系数
因为本设计气液分离器是等直径厚壁设备,
故
由[2]表16-5得
=1,
----------------直立设备有效直径,
------------------容器各段外直径,mm
------------------容器各段保温层厚度,mm
-------------------笼式扶梯当量宽度,mm
--------------------操作平台当量宽度,mm
=
m
=
同理可求出分离器各段风载荷于下表所示
计算各段载荷
0-1
1-2
2-3
0.7
1
0.35
0.35
0.35
1.35
1.35
1.35
(kgf)
60
0.78
0.78
0.78
0.8
2
1.2
平台数
0
1
1
0
0.37
0.617
0.312
1.05
1.546
11.04
92.87
82.05
(2)计算风弯矩
容器任意计算截面I-I的风弯矩计算公式
:
(69)
分离器各截面的风弯矩计算结果如下表:
截面
弯矩
0-0
1-1
2-2
6.6地震弯矩计算
任意截面的地震弯矩
:
(70)
等直径、等壁厚设备任意截面计算截面I-I地震弯矩
:
(71)
底部截面的地震弯矩:
(72)
(当
时,视设备为柔性结构,则需考虑高振型的影响,在进行稳定或其他验算时,所取的地震弯矩值应为上列计算值的1.25倍)
式中:c--------------------------结构影响系数,对直立圆筒形容器取c=0.5
----------------------地震影响系数,这里取最大地震时影响系数
=0.945
Q0 ----------------------自计算截面以上的操作载荷
----------------------塔顶至第 I段(含第i段)的高度,m
;根据GB12470-70选用焊剂的牌号HJ431型号HJ401-H08A.
3.焊接采用埋弧焊,对焊缝进行100%的射线探伤检测,要求符合GB3325-87中的Ⅱ级为合格.
7.2接管与筒体焊接结构设计
1. 筒体与接管的焊接接头型式和尺寸选用G2 GB20583-1998,其示结构如图2所示.
2. 根据GBT518-95《低合金钢电焊条》选用焊条牌号J502,型号E5003。
3. 焊接采用手工电弧焊,对焊缝进行100%的射线探伤检测,要求符合GB3325-87中的Ⅱ级为合格.
7.3带颈平焊法兰与接管焊接结构设计
1. 带颈平焊法兰与接管焊接接头尺寸选用F6(JB4700~4703-92),其结构如图4所示。
2.根据GBT518-95《低合金钢电焊条》选用焊条牌号J507,型号E5015。
3.焊后对焊缝进行100%的射线探伤检测,要求符合GB 3325-87中的Ⅱ级为合格.
第八章 气液两相流场的数值模拟
8.1 数值计算方法简介
计算流体力学作为流体力学研究中的一门新兴分支,正在工业和科研领域内发挥越来越重要的作用。将CFD工具运用到分离机械的研究中,也成为工程技术人员改进设计、提高效率的有效手段,是CFD应用的前沿。一些成熟的算法,模型也以商业软件的形式出现在工程及科研领域。相比研究单位自行开发的计算程序,商业计算软件一般具有以下特点:
· 通用性广。由于商业软件面向的用户对象广泛,处理的实际问多种多样,因此其覆盖的应用范围要尽可能广。
· 计算稳定性好。多数软件经过不同研究领域内的算例测试,对不同类型的问题具有较好的适应能力。
· 使用方便,商业软件都提供了比较友好的用户界面,方便用户的使用。
一般商业软件也存在一些明显的不足,例如:算法相对陈旧,不能紧跟CFD研究领域内的最新成果;与不同行业内的实际要求存在一定的距离,难以将各研究单位已有的研究成果结合到商业软件中。这在一定程度上限制了商业软件在工程实际中的应用。
本文运用Fluent软件对离心是分离器的内流场进行分析计算,Fluent公司是享誉世界的最大计算流体力学(CFD)软件供应商,Fluent软件能够精确地模拟无粘流、层流、湍流、化学反应、多相流等复杂的饿流动现象。应用领域包括:航空航天、汽车设计、生物医药、化学处理、石油天然气,发电系统、电子半导体、蜗轮设计、HVAC、玻璃加工等。FLUENT具有精度高,收敛快,稳定性好等特点,同时可通过添加拥护自定义的函数(UDF)解决实际具体问题。
Gambit是前置处理器,能针对极其复杂的几何外行生成三维四面体,六面体的非结构化网络及混合网络。该模块还具有方便的网络检查功能,对网络单元体积、扭曲率、长细比等影响收敛和稳定的参数进行统计并生成报告。
8.1.1 控制方程
对于所有的流动问题,FLUENT需要求解质量和动量守恒方程。对于热传导或可压缩流动,需要解能量守恒的附加方程。对于包括组分混合和反应的流动,需要解组分守恒方程或者使用PDF模型来解混合分数的守恒方程。当流动是湍流是,还要解附加的输运方程。
FLUENT可以在惯性坐标系(无加速坐标系)和具有加速度的参考坐标系(旋转坐标系)中建立流动模型。在旋转坐标系中,通过建立一个与旋转设备一起运动的相对坐标系来建摸,近似认为流体旋转角速度为常数,旋转边界相对与参考系静止。
在旋转坐标系中,绝对速度v或相对速度v的关系如下:
(80)
其中
─角速度向量(即旋转坐标系的角速度);
─旋转坐标系中的位置向量
旋转坐标系中的质量守恒方程(连续性方程):
(81)
其中
源项
上假如到连续性相的第二相质量(比方说由于液滴的蒸发,质量发生变化),源项也可以是自定义源项。
惯性坐标系中的动量方程:
(82)
其中p是静压,τ是应力张量,ρ
和
分别是中立体积力和外部体积力(如离散相相互作用的产生的升力)。
包含了其它的模型相关源项,如多孔介质和自定义源项。
应力张量
(83)
其中:μ─分子粘性;
I─单位张量。
旋转坐标系中的动量方程为:
(84)
由于在FLUENT中忽略了
项,因此不能用动量方程的相对速度表达式准确的计算随时间变化的角速度。
8.1.2湍流模型
湍流是由大小不同尺度的旋涡组成,对时间和空间都是非线性的随机运动。它最本质的特征是“湍动”,既随机的脉动。湍流流场是无数不同尺寸涡旋相互掺混的流动场。湍流的出现影响着整个流场的速度、压力、温度和物质浓度的分布。
目前,关于紊流的数值计算可分为细观模拟和统观模拟。完全模拟、大涡流模拟属于细观模拟。完全模拟(Direct Numberical Simulation DNS)在湍流尺度的网格尺寸内求解N ─S方程而不使用任何湍流模型。该方法必须采用很少的时间和空间步长。如文献估算,长。如文献估算,对一个涡旋进行数值计算,至少要设置十个节点,这样对于一个小尺寸范围内的紊流运动要在1cm 的流场中布置10个节点。显然,完全模拟在短期尚无法用于求解实际工程中的复杂湍流问题。大涡流模拟(Large Eddy Simulation LES)是在大涡尺度的网格尺寸内求结N─S方程,由于计算量仍很大,只能模拟一些简单流动,如弯道等,目前也不能直接用于工程。而基于求解雷诺时均方程的模拟,即统观模拟,利用某些假设将雷诺时均方程中高阶未知关联项或者时均量来表达,从而使雷诺时均方程封闭。因此统观模拟当前成为解决工程实际问题的有效手段。
所谓湍流模型理论就是以雷诺平均运动方程为基础,依靠理论与经验引进一些模拟假设,建立一组描述湍流平均值的封闭方程组的计算方法。该模型的平均行为,应与实际的湍流统计平均行为基本一致。
1925年Prandtl提出的动量传递理论以及后来提出的自由剪切层模型、泰勒的涡量传递理论的冯-卡门的相似性理论等一系列半经验理论,其基本思想都是建立关于雷诺应力的模型假设,使雷诺平均运动方程得以封闭。由于只考虑了一阶湍流统计量的动力学微分方程,即平均运动方程,仅引人附加的代数关系而没有引进任何高阶统计量的微分方程,因此属于零方程模型。这种模型对运动的预测性很差,不适用于有回流的复杂运动。
Kolmogorov和Prandtl提出的单方程模型虽较半经验理论有所改进,且Cd和Cu值也较容易确定,但L值的确定并不比混合长度的确定容易。因而单方程模型同样只适用于简单流动,不适用于带回流的复杂流动。
在所有的双方程模型中,K-ε双方程模型的应用最为普遍,先后由周培源(1945)、Harlow-Nakayama(1968) 、Jones-Launder(1972)提出。在进一步简化的模型中,人们放弃给雷诺应力建立方程的想法,将它们直接用推广的Bossinesq涡粘性模式来表示。大量的预报及实验结果对照表明,K-ε模型可以成功或基本成功用于以下几种情形:无浮力平面射流、平壁边界层、管流、通道流或喷管内流动、无旋涡及弱旋的二维和三维流动。由于它采用了同向性湍流疏运的假设,故它不适应具有非同向性湍流输运的强旋流。
在离心分离器内流场中,因切速度远远大于径向和轴向分量,因而,通常采用雷诺应力模型(RSM)、代数应力模型(ASM)或RNG模型来代替K-ε模型。
雷诺应力模型RSM对雷诺应力及通量项采用微分方程直接求解,具有很大通用性,这一模型的优点在于可准确地考虑各向异性效应。
代数应力模型(ASM)是雷诺应力模型在一定条件下的简化表达式,表达式形式随简化条件而异,由应力代数表达式加上κ及ε方程构成,因此又称为扩展的K-ε模型或者2个半方程模型。应用该模型可避开求解雷诺应力方程所面临的十分复杂的计算工作,又较好地预报了流动的各向异性的特点。
RNG K-ε模型是一种修正的K-ε模型,在文献[33][34]中有较详细的讨论。它从原始的基本方程推导而来,其中使用了所谓的“Kolmogorov”数学技巧。Kolmogorov定律:E(k)=k 。RNG模型是一个更一般,更基本的模型,尤其对强旋流流场及高曲率流线的离心分离器有着很好的改进效果。
在FLUENT中有多种湍流模型可选择,包括常用的Spalart-Allmaras单方程模型、标准K-ε模型、和K-o模型、雷诺应力模型以及大旋涡模型。选择合适的湍流模型是数值计算中及其重要的环节,不同湍流模型的使用会产生不同的流场细节。本文应用RNGK-ε模型预测离心式分离器的内流场。
8.1.3多相流模型
多相流的数值计算方法有两种:欧拉—拉格朗日(Euler-Lagrange)方法和欧拉—欧拉(Euler- Euler)方法。
1.拉格朗日离散项模型
在Fluentzhong中的拉格朗日离散项模型遵循欧拉—拉格朗日方法。流体相为连续相,直接求解时均N-S方程,而离散相通过计算流场中大量的粒子,气泡或是液滴的运动得到。离散相和流体相之间有动量、质量和能量的交换。该模型假设离散相(第二相)的体积比率很低。粒子或液滴运行轨迹的计算是独立的,被安排在流相计算的指定间隙完成。这样处理能较好的符合喷雾干燥,煤和液体燃料燃烧,以及一些粒子负载的流动情况,但是不适合用于液-液混合物,流化床和其它第二相体积率不容忽略的情形。
2.欧拉—欧拉模型
在欧拉—欧拉方法中,不同的相被处理成互相贯穿的连续介质。由于一种相所占的体积不能被其他相占有,故引入相体积分数(phasic volume fraction)。体积分数是时间和空间的连续函数,各相的体积分数之和等于1。从各相的守恒方程可以推导出一组方程,这些方程对于所有的相都具有类似的形式。从实验得到的数据可以建立一些特定的关系,从而使上述方程封闭,另外,对于小颗粒流(granular flows),可以通过应用分子运动论的理论使方程封闭。在FLUENT中,有三种欧拉—欧拉多相流模型:流体体积模型(VOF),混合(Mixture)模型以及欧拉(Euler)模型。混合模型和欧拉模型主要用于模拟相间的混合和分离。
(1)流体体积模型(VOF)
VOF模型通过求解单独的动量方程和处理穿过区域的每一相的体积分数来模拟两种或三种不能混合的流体。典型的应用包括预测流体中大气泡的运动和气液界面的稳态和瞬态变化。
(2)混合模型(Mixture)
混合模型求解混合物的动量方程,并通过相对速度描述离散相。混合模型是欧拉模型在几种情况下的很好替代。当颗粒相广泛分布或界面规律未知时,完善的多相流模型是不切实际的。在这种情况下,混合模型能取得较好的结果。
混合模型的应用包括低负载的粒子负载流,气泡流,沉降,以及旋风分离等多相流。混合模型也可用于没有相对速度的均匀多相流。
(3)欧拉模型(Eulerian)
欧拉模型是Fluent中最复杂的多相流模型。它建立含有n个动量方程和连续方程的方程组来求解各相。压力项和各界面交换系数耦合在一起,耦合方式依赖于所含相的情况。颗粒流(液—固)与非颗粒流(液—液)的处理是不同的。对于颗粒流,可应用分子运动理论求得流动特性。不同相之间的动量交换也依赖于混合物的类别。该模型的应用包括气泡柱,上浮,颗粒悬浮,以及流化床。
3.多相流求解策略
解决多相流问题,先选择最能符合实际的流体模式。然后根据不同的模式,选择恰当的多相流模型。如在气—液或液—液多相流中,大致有泡状流,弹状流,断塞流及自由液面流四种模式。对于离心分离器,主要是分离在来连续液相中分散的小气泡(含气率为10%左右),气泡的运动一般认为是泡状流,可选用混合模型或欧拉模型。以下为混合模型和欧拉模型的求解策略。
启动混合模型求解,采用0.2或更小的滑流速度欠松弛因子开始计算。如果解显示出好的收敛趋势,可逐渐增加欠松弛因子。初始计算是可以不求解题解分数和滑流速度方程,当计算收敛后,在求解这些方程。
启动欧拉模型求解,为了提高收敛性,在求解多相流模型前可以先获得初始解再继续计算。初始解可通过混合模型得到。由于欧拉模型对内存及收敛要求较高,在下面的数值计算中,多相流模型选用混合模型。
8.1.4 数值计算方法
目前关于数值计算的方法有:有限差分法、有限元法、边界元法等。在CFD中,有限差分法占主导地位,且最为成熟,目前已经发展了多种收敛性好、精度高的离散格式,较常用的有:Taylor(泰勒)展开法、有限容积法。与有限差分法相比,有限单元法在流体力学应用还远远不够成熟,易出现数值溢出、数值发散等问题。Flent采用有限容积离散动量方程,速度和压力耦合采用SIMPLE、SIMPLEC及PISO算法。对流项差分格式纳入了一阶迎风、中心差分及QUICK等格式。在多种情况下,中心格式的计算稳定性好于迎风格式。代数方程可以采用多重网格及最小残差法。
对于旋转坐标系中的流动问题,求解过程经常变得不稳定。这是因为当旋转项的影响很大时,动量方程高度耦合。高速的旋转会产生一个很大的径向压力梯度,从而推动流体向轴向和径向流动,并在流场中形成旋涡或旋度的分布。为了改善求解的稳定性,可以采用以下方法:
●采用PRESTO方法离散压力项,这种方法非常适合于旋流中的大压力梯度。
●减小速度的欠松弛因子,将设为0.3~0.5或更低。
●用一个较低的旋转速度开始计算,然后逐渐增加,最后获得合理的结果。
8.1.5 边界条件的处理
边界条件包括流动变量在在边界处的值。FLUENT中可处理的边界条件包括:进口边界条件(压力、速度、质量进口、进风口、进气扇、压力出口、压力远场边界条件,质量进口、通风口、固体壁面、对称面、周期边界、轴,内部区域:流体、固体;内部表面边界:风扇、换热器、多孔跳跃、壁面、内部类型。
8.2 计算前处理
计算前处理主要包括离心式分离器三维流道模型的建立及网格的划分。在本文采用FLLUENT的前处理器Gambit来实现。
8.2.1 建立三维模型
图8-1
根据理论计算尺寸建立分离器的流道模型。除了直接在Gambit工作环境下建立模型外,还可以利用其他CAD软件如UG、PROE、CAXA、SOLID WORKS等生成后缀为.igs的图形文件。通过软件之间的接口,调入到Gambit内进行网格划分。
本章中用于计算的几何模型具体结构参数如下:
离心式气液分离器入口管直径:
=40mm
离心式气液分离器入口管当量长度:
=1000mm
离心式气液分离器气体出口管直径:
=20mm
离心式气液分离器气体出口管当量长度:
=200mm
离心式气液分离器液体出口管直径:
=65mm
离心式气液分离器液体出口管当量长度:
=1614mm
离心式气液分离器壳体内径:d=300mm
离心式气液分离器总体高度:H=3300mm
8.2.2网格划分
FLUENT作为一种求解器,在二维问题中可以使用三角形、四边形或混合单元组成的网格,在三维问题中可以使用四面体、六面体、金字塔形、楔形单元,或者两种单元的混合。Gambit具有强大的自动生成网格的功能,能快速生成结构化、非结构化网格。由于计算区域形状比较复杂,本文采用C结构化网格(混合网格)。
在分离器进口管附近的流动区域,流场的速度场、压力场梯度比较大,网格要求加密。同时,液液出口管径和气体出口管径相对较小,网格尺寸也要相应的变小.因为技术计算机的内存有限,流动区域的其他地方的网格相应的较疏,以保证网格数不会太大。在Gambit工作环境中,可以利用size function控制网格尺寸由以上加密的位置向空间增大。图8-2为用此法生成的计算模型。
图8-2 CFD模型
8.3求解方法和边界条件
数值计算使用三维粘性雷诺时均N-S方程求解软件FLUENT,运用分离的隐式求解方法、RNGk-e模型,混合的多相流模型,同时考虑了离心式气液分离器内流体的强选旋流动,对压力离散采用PRESTO法,对动量方程采用二阶迎风差分格式进行离散,对速度和压力耦合采用SIMPLE算法。
按不同的求解方法,FLUENT的耦合式和分离式。在分离式求解法中,独立有序地求解速度和压力各变量。
边界条件的确定也是数值计算的关键。在本文中需要涉及到壁面边界和进出口边界。进口的边界取混合液的入口速度,同时确定第二相的体积分数。气体出口和液体出口设为压力出口,表压为0。固体壁面定义为无滑移壁面,筒体的内表面为旋转壁面。对于紊流需要说明性质。假设流场为全湍流场,速度入口的流动充分的发展,入口的湍流强度为5%,湍流粘性比为10%。当残差下降到
时,并且进、出口的流量误差小于5%时认为计算收敛。
8.4数值计算结果
气液分离器内主要包括速度场、压力场及各相体积分数的分布。下面以本设计工况下分离器内流场分布为例分析其特点。
工作介质为一定含气量的钻井液,每天的处理量:液体900
d,气体90
d,入口流速为6ms,其中A为入口管截面面积。钻井液密度1500kg
,粘度0.05 Pa.s来液中气体体积分数为9.09%,气泡直径0.8mm,气体的密度1.185kg
以下为静压力场、速度场以及体积分数的分布情况.
图8-3 纵向剖面泥浆压力场分布图
图8-4混合液的速度场
图8-5 混合液的速度云图
图8-6 混合液的壁面速度分布图
图8-7进口附近壁面速度场
图8-8 Z=1.0m横向截面流场矢量图
图8-9 Z=1.7m横向截面流场矢量图
图8-10 排液口流场矢量图
图8-11 排气口附近流场流速矢量分布图
由图8-3可以看出分离器内整体压力分布情况,在分离器的底部压力最大,在液相出口压力最小。由图8-4、8-5、8-6可以很明显的看出钻井液进入分离气后,钻井液的速度场在进口处和液相出口的速度最大。由图8-7、8-8、8-9可以清楚的看出钻井液通过进口管进入分离器后确实产生产生了强烈旋流,沿径向速度矢量逐渐增大。图8-10、8-11逼真的反映了液相和气相的速度场的矢量分布,气相、液相流向各自的出口,在离心式分离器中实现了对气液的分离。
混合液发生气液分离后,各相按一定的比例浓度分布。液相的体积分数如图8-12所示。由图8-12可以看出,气液混合物经过分离后,在分离器的上半部分泥浆的体积分数几乎是零,泥浆在分离器的下半部分和排液管中的体积分数明显高于进口管附近的体积分数,特别是在分离器的底部和排液管线中泥浆的体积分数几乎是100%,十分有利的证明了钻井液在分离器内实现了气液的分离。由图8-13可以看出在气体出口管内充满气体,携带液量可以忽略。在垂直于旋转轴的截面内,含液量随着半径的增大而增大,这是因为离心力与半径的平方成正比关系,在强旋流场中,高密度的液相沿径向向外运动抛向壁面,较小密度的气相向中心聚集转而向上运动。
图8-12 纵向剖面泥浆体积分数分布图
图8-13纵向剖面气体体积分数分布
8.5计算分离效率
分离效率计算公式
:
式中:
--------------------------------分离器的分离效率,%
-------------------------------除气后钻井液密度,
-------------------------------不含气体的钻井液密度,
-------------------------------气侵后的钻井液密度,
已知
设计气液分离器的分离效率:
=
第九章 结束语
本文从结构上对用于钻井中泥浆处理的旋流式分离器进行了设计;并对分离器的各重要零部件进行了设计、强度计算以及流场分析。该管柱式旋流式分离器可分离气泡小,效率较高,分离器体积相对较小,维修量小,工作稳定。旋流式分离器具有这些特点,使得该分离器在钻井过程中用于泥浆的处理具有实用性和可靠性,并具有广阔的前景。
参考文献
[1]董大勤,袁凤隐主编,压力容器与化工设备实用手册,第一版。 北京:化学工业出版社,2003.2
[2]贺匡国主编,化工容器及设备简明设计手册,北京:化学工业出版社,1989.12
[3]卓震主编,化工容器及设备,第一版,北京:中国石化出版社,2000.2
[4]邹广华,刘强主编,过程装备制造及检测,第一版,北京:化学工业出版社,2003.8
[5]张德姜,赵勇主编,石油化工管道设计与安装,第一版, 北京:中国石化出版社,2002.2
[6]韩占中,王敬,兰小平主编,FLUENT-流体工程仿真计算实例应用,第一版, 北京:北京理工大学出版社,2004.6
[7]魏崇光,郑晓梅主编,化工工程制图,第一版, 北京:化学工业出版社,2000.8
[8]雷文平等主编,Solid Works 2001,第一版, 北京:清华大学出版社,2002.1
[9]刘朝儒,高政一等主编,机械制图,第四版,北京:高等教育出版社 ,2000.12
[10]成大先主编,机械设计手册,第四版,北京:机械工业出版社,2003.1
[11]曹学文, 黄庆宣 新型管柱式气液旋流分离器 天然气工业 2002 22(2)
[12]学成, 山本恭二 开孔圆筒式液体旋流分离器的研究 通风除尘 1998 17(1)
[13]朱浩东, 杨敏 国内外旋流分离器特点及发展方向 石油机械 1994 22(12)
[14]陈丽萍. 立式重力气液分离器工的工艺设计 天然气化工:C1化学与化工 1999, 24(2)
[15]K.E.Arnold, 王晓凤. 李太平结构紧凑的分离器系列 国外油田工程 2000,(8)
[16]贾建贞, 张宝彦, 刘景云. 欠平衡钻井液气分离器的研制与应用 石油机械 2001,(7)
[17]M.F.Schubert, 李太平. 水力旋流分离系统的发展 国外油田工程 1999,(7)
[18]袁惠新, 王跃进. 旋流分离技术在石油、石化工业中的应用 化工设备与防腐蚀 2002,(3)
[19]曹学文,黄庆宣. 新型管柱式气液旋流分离器设计与应用, 天然气工业 2002, 22(2)
[20]吴粤燊主编, 压力容器安全技术手册, 北京:机械工业出版社,1999.3.
[21]张康达,洪起超主编. 压力容器手册(上下), 北京:劳动人事出版社, 1988.2
[22]家祯主编, 压力容器材料实用手册(第一版), 北京:化学工业出版社工业装备和信息工程出版社 , 2000. 7
[23]吴粤荣主编 . 压力容器安全技术(第二版), 北京:化学工业出版社, 1993.11
[24]董振仁,魏新立主编. 过程装备成套技术,第一版。 北京: 化学工业出版社, 2001. 11
[25]郭年祥主编, 化工过程及设备, 北京:冶金工业出版社, 第一版,2003.3.
[26]郑津洋, 董其五, 桑芝富主编, 过程设备设计,第一版, 北京:化学工业出版社, 2001.3
[27] (美)耶帝什.特.夏编, 气液固反应器设计, 北京:烃加工出版社 , 1989. 9.
[28]陈国理主编, 压力容器及化工设备, 广州:华南理工大学出版社, 1988.11
[29]余国宗主编, 化工机械工程手册, 北京: 化学工业出版社, 2001. 7
[30]潭天恩等主编,化工原理,第二版. 北京:化学工业出版社, 2003. 1
[31]袁恩熙主编,工程流体力学,第一版. 北京:石油工业出版社,1986 .1
[32]同济大学主编,高等数学,第四版。北京:高等教育出版社。2002. 2
[33]龚伟安, 钻井液固相控制技术及设备,北京:石油大学出版社,1995 .4
毕业设计(论文)原创性声明和使用授权说明
原创性声明
本人郑重承诺:所呈交的毕业设计(论文),是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知,除文中特别加以标注和致谢的地方外,不包含其他人或组织已经发表或公布过的研究成果,也不包含我为获得 及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体,均已在文中作了明确的说明并表示了谢意。
作 者 签 名: 日 期:
指导教师签名: 日 期:
使用授权说明
本人完全了解 大学关于收集、保存、使用毕业设计(论文)的规定,即:按照学校要求提交毕业设计(论文)的印刷本和电子版本;学校有权保存毕业设计(论文)的印刷本和电子版,并提供目录检索与阅览服务;学校可以采用影印、缩印、数字化或其它复制手段保存论文;在不以赢利为目的前提下,学校可以公布论文的部分或全部内容。
作者签名: 日 期:
学位论文原创性声明
本人郑重声明:所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。除了文中特别加以标注引用的内容外,本论文不包含任何其他个人或集体已经发表或撰写的成果作品。对本文的研究做出重要贡献的个人和集体,均已在文中以明确方式标明。本人完全意识到本声明的法律后果由本人承担。
作者签名:
日期: 年 月 日
学位论文版权使用授权书
本学位论文作者完全了解学校有关保留、使用学位论文的规定,同意学校保留并向国家有关部门或机构送交论文的复印件和电子版,允许论文被查阅和借阅。本人授权 大学可以将本学位论文的全部或部分内容编入有关数据库进行检索,可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。
涉密论文按学校规定处理。
作者签名:
日期: 年 月 日
导师签名: 日期: 年 月 日
致 谢
时间飞逝,大学的学习生活很快就要过去,在这四年的学习生活中,收获了很多,而这些成绩的取得是和一直关心帮助我的人分不开的。
首先非常感谢学校开设这个课题,为本人日后从事计算机方面的工作提供了经验,奠定了基础。本次毕业设计大概持续了半年,现在终于到结尾了。本次毕业设计是对我大学四年学习下来最好的检验。经过这次毕业设计,我的能力有了很大的提高,比如操作能力、分析问题的能力、合作精神、严谨的工作作风等方方面面都有很大的进步。这期间凝聚了很多人的心血,在此我表示由衷的感谢。没有他们的帮助,我将无法顺利完成这次设计。
首先,我要特别感谢我的知道郭谦功老师对我的悉心指导,在我的论文书写及设计过程中给了我大量的帮助和指导,为我理清了设计思路和操作方法,并对我所做的课题提出了有效的改进方案。郭谦功老师渊博的知识、严谨的作风和诲人不倦的态度给我留下了深刻的印象。从他身上,我学到了许多能受益终生的东西。再次对周巍老师表示衷心的感谢。
其次,我要感谢大学四年中所有的任课老师和辅导员在学习期间对我的严格要求,感谢他们对我学习上和生活上的帮助,使我了解了许多专业知识和为人的道理,能够在今后的生活道路上有继续奋斗的力量。
另外,我还要感谢大学四年和我一起走过的同学朋友对我的关心与支持,与他们一起学习、生活,让我在大学期间生活的很充实,给我留下了很多难忘的回忆。
最后,我要感谢我的父母对我的关系和理解,如果没有他们在我的学习生涯中的无私奉献和默默支持,我将无法顺利完成今天的学业。
致 谢
四年的大学生活就快走入尾声,我们的校园生活就要划上句号,心中是无尽的难舍与眷恋。从这里走出,对我的人生来说,将是踏上一个新的征程,要把所学的知识应用到实际工作中去。
回首四年,取得了些许成绩,生活中有快乐也有艰辛。感谢老师四年来对我孜孜不倦的教诲,对我成长的关心和爱护。
学友情深,情同兄妹。四年的风风雨雨,我们一同走过,充满着关爱,给我留下了值得珍藏的最美好的记忆。
在我的十几年求学历程里,离不开父母的鼓励和支持,是他们辛勤的劳作,无私的付出,为我创造良好的学习条件,我才能顺利完成完成学业,感激他们一直以来对我的抚养与培育。
最后,我要特别感谢我的导师刘望蜀老师、和研究生助教吴子仪老师。是他们在我毕业的最后关头给了我们巨大的帮助与鼓励,给了我很多解决问题的思路,在此表示衷心的感激。老师们认真负责的工作态度,严谨的治学精神和深厚的理论水平都使我收益匪浅。他无论在理论上还是在实践中,都给与我很大的帮助,使我得到不少的提高这对于我以后的工作和学习都有一种巨大的帮助,感谢他耐心的辅导。在论文的撰写过程中老师们给予我很大的帮助,帮助解决了不少的难点,使得论文能够及时完成,这里一并表示真诚的感谢。
本科生毕业设计(论文)化要求
第一部分 学生应遵守以下规范要求
一、毕业设计论文说明
1. 毕业设计论文独立装订成册,内容包括:
(1) 封面(题目、学生姓名、指导教师姓名等)
(2) 中、外文内容摘要
(3) 正文目录(含页码)
(4) 正文(开始计算页码)
(5) 致谢
(6) 参考文献
(7) 附录
2. 中、外文内容摘要包括:课题来源,主要设计,实验方法,本人主要完成的成果。要求不少于400汉字,并译成外文。
3. 毕业设计论文页数为45页-50页。
4. 纸张要求:毕业设计说明书(论文报告)应用标准B5纸单面打字成文。
5. 文字要求:文字通顺,语言流畅,无错别字。
6. 图纸要求:毕业设计图纸应使用计算机绘制。图纸尺寸标注应符合国家标准。图纸应按“规范”叠好。
7. 曲线图表要求:所有曲线、图表、流程图、程序框图、示意图等不得徒手画,必须按国家规定标准或工程要求绘制。
8. 参考文献、资料要求:参考文献总数论文类不少于10篇、,应有外文参考文献。文献应列出序号、作者、文章题目、期刊名、年份、出版社、出版时间等。
二、外文翻译
1. 完成不少于2万印刷符的外文翻译。译文不少于5千汉字。
2. 译文内容必须与题目(或专业内容)有关,由指导教师在下达任务书时指定。
3. 译文应于毕业设计中期2月底前完成,交指导教师批改。
4. 将原文同译文统一印成B5纸规格装订成册,原文在前,译文在后。
三、形式审查
5月15日前,将毕业设计论文上交指导教师,审查不合格者,不能参加答辩。
四、准备答辩
答辩前三天,学生要将全部材料(包括光盘、论文)统一交指导教师。
关于毕业论文格式的要求
为方便统一、规范论文格式,现将学院的相关要求做如下强调、补充:
1. 基本要求
纸型: B5纸(或16开),单面打印;
页边距: 上2.54cm,下2.54cm,左2.5cm,右2.5cm;
页眉:1.5cm,页脚1.75cm,左侧装订
正文字体:汉字和标点符号用“宋体”,英文和数字用“Times New Roman”,字号小四;
图号1-1,指第1章第1个图
在图的前部要有文字说明(如图1-1所示)
表号3-5,指第3章第5个表
在表的前部要有文字说明(如表3-5所示)
图、表的标注字体大小是五号宋体
行距: 固定值20;
页码: 居中、小五、底部。
2. 封面格式
封皮: 大连理工大学城市学院(二号、黑体、居中)
本科生毕业设计(论文)(二号、黑体、居中)
学 院:(四号、黑体、居中、下划线:电子与自动化学院)
专 业:(四号、黑体、居中、下划线、专业名字之间无空格)
学 生:(四号、黑体、居中、下划线,名字是2个字的中间空1个字、3个或3个以上字的中间无空格)
指导教师:(四号、黑体、居中、下划线,名字是2个字的中间空1个字、3个或3个以上字的中间无空格,两位指导教师的中间用顿号“、”)
完成日期:(四号、黑体、居中、下划线,如:2009年5月25日)
(注意:5个下划线两端也是对齐的,单倍行距)
内 封:大连理工大学城市学院本科生毕业设计(论文)(四号、黑体)
题目 (二号、黑体、居中);
总计 毕业设计(论文) 页(五号、宋体)
表格 表(五号、宋体)
插图 幅 (五号、宋体)
(注意:页数正常不少于40页,优秀论文原则上不少于45页)
3. 中外文摘要
中文摘要:标题“摘 要” (三号、黑体、居中、中间空1个字)
正文(不少于400字)
关键词 (五号、黑体):3-5个主题词(五号),中间用分号“;”隔开。
外文摘要 (另起一页):标题“Abstract” (三号、黑体、居中)
正文 (必须用第三人称)
关键词: Key words(五号、黑体):3-5个主题词(五号)与中文关键词对应,中间用分号“;”隔开。
4. 目录
标题 “目录”(三号、黑体、居中);
章标题(四号、黑体、居左);
节标题(小四、宋体);
页码 (小四、宋体);
二、三级目录分别缩近1和2个字;
四级目录不在“目录”中体现,在正文中也不是单独一行,可以黑体(没有句号),然后空2个字接正文;
注意:正文中每章开头要另起一页;
“目录”下方中间的页码和摘要一样统一用罗马字,顺接摘要的。
摘要 目录加页眉
5. 论文正文
页眉: 论文题目(居中、小五、黑体);
章标题(三号、黑体、居中);
节标题(四号、黑体、居左);
正文
程序用“Times New Roman”,字号小四;
6. 参考文献
标题:“参考文献”(小四、黑体、居中)
参考文献的著录,按文稿中引用顺序排列,并注意在文内相应位置用上标标注,如:……的函数。
示例如下:(字体为五号、宋体)
期刊类:[序号]作者1,作者2,……作者n。文章名。期刊名(版本),出版年,卷次(期次)。页次
图书类:[序号]作者1,作者2,……作者n。书名。版本。出版地:出版者,出版年。页次
会议论文集:[序号]作者1,作者2,……作者n。论文集名。出版地:出版者,出版年。页次
网上资料:[序号]作者1,作者2,……作者n。文章名。网址。发表时间
7. 其它
量和单位的使用:必须符合国家标准规定,不得使用已废弃的单位(如高斯(G和Gg)、亩、克分子浓度(M)、当量能度(N)等)。量和单位不用中文名称,而用法定符号表示。
图表及公式:插图宽度一般不超过10cm,表名(小四)置上居中,图名(小四)置下居中。标目中物理量的符号用斜体,单位符号用正体,坐标标值线朝里。标值的数字尽量不超过3位数,或小数点以后不多于1个“0”。如用30Km代替30000m,用5µg代替0.005mg等,并与正文一致。图和表的编号从前至后顺序排列,图的编号及说明位于图的下方,居中;表的编号及说明位于表的上方,居中。公式编号加圆括号,居行尾。图表中的字体不应大于正文字体。注意:图表标题中的数字也是“Times New Roman”。
8.论文依次包括:封皮、内封、中文摘要、英文摘要、目录、正文、结论、致谢、参考文献、(附录),不要落项。
9.注意:上面没有说“加粗”的“黑体”,均为“黑体不加粗”。
补充:
1.答辩要求:自述15分钟,回答问题10分钟,自述要求使用PPT
答辩内容:
1).论文题目
2).设计内容
3).设计方案
4).如何完成设计
工作原理
软件或硬件设计
制作\调试\安装
5).存在不足,今后努力的方向
6).致谢
3.最后上交学生装订好的论文、光盘、表、成绩单
4.光盘里的文件夹命名为:学号_姓名_年级专业班级
文件夹里包括的文件有:论文、ppt、英文翻译
1) 论文的文件名格式:学号_姓名_年级专业班号_题目(论文)_完成日期doc
2) ppt的文件名格式:学号_姓名_年级专业班号_题目(ppt)_完成日期ppt
3) 英文翻译的文件名格式:学号_姓名_年级专业班号_题目(英文翻译)_完成日期doc
例如:
答辩问题5个, 侧重总体思路一个
软件或硬件一个
翻译一个
其他2个
� EMBED Equation.3 ���
�
2
62
_1179610761.unknown
_1179923386.unknown
_1181202859.unknown
_1181299557.unknown
_1235387458.unknown
_1235387634.unknown
_1235387716.unknown
_1235387797.unknown
_1235387833.unknown
_1235396813.unknown
_1235387804.unknown
_1235387826.unknown
_1235387779.unknown
_1235387791.unknown
_1235387763.unknown
_1235387669.unknown
_1235387685.unknown
_1235387644.unknown
_1235387589.unknown
_1235387632.unknown
_1235387488.unknown
_1235387576.unknown
_1235367480.unknown
_1235387318.unknown
_1235387363.unknown
_1235387071.unknown
_1235367147.unknown
_1235367382.unknown
_1181299780.unknown
_1181299827.unknown
_1181293574.unknown
_1181296829.unknown
_1181298491.unknown
_1181299537.unknown
_1181299552.unknown
_1181298200.unknown
_1181296767.unknown
_1181296817.unknown
_1181296555.unknown
_1181268169.unknown
_1181273237.unknown
_1181284028.unknown
_1181268185.unknown
_1181268393.unknown
_1181229300.unknown
_1181229588.unknown
_1181268074.unknown
_1181268119.unknown
_1181267830.unknown
_1181229932.unknown
_1181229507.unknown
_1181229565.unknown
_1181229355.unknown
_1181229132.unknown
_1181229018.unknown
_1181229067.unknown
_1181228889.unknown
_1181203091.unknown
_1181203246.unknown
_1181203308.unknown
_1181203330.unknown
_1181203344.unknown
_1181203261.unknown
_1181203153.unknown
_1181202997.unknown
_1181203036.unknown
_1181202879.unknown
_1180130820.unknown
_1180151388.unknown
_1180194140.unknown
_1181198611.unknown
_1181199170.unknown
_1181202841.unknown
_1181198642.unknown
_1181198708.unknown
_1180194939.unknown
_1181198489.unknown
_1181198539.unknown
_1181181921.unknown
_1181182576.unknown
_1180195012.unknown
_1180194406.unknown
_1180194692.unknown
_1180194233.unknown
_1180188322.unknown
_1180188426.unknown
_1180188504.unknown
_1180188536.unknown
_1180194021.unknown
_1180188520.unknown
_1180188445.unknown
_1180188367.unknown
_1180188402.unknown
_1180188349.unknown
_1180152548.unknown
_1180188297.unknown
_1180152011.unknown
_1180152201.unknown
_1180151756.unknown
_1180131658.unknown
_1180132363.unknown
_1180132705.unknown
_1180151090.unknown
_1180132410.unknown
_1180131805.unknown
_1180132232.unknown
_1180131770.unknown
_1180130882.unknown
_1180131532.unknown
_1180131327.unknown
_1180131375.unknown
_1180131137.unknown
_1180130837.unknown
_1180130846.unknown
_1179923831.unknown
_1179928678.unknown
_1180130642.unknown
_1180130660.unknown
_1180130790.unknown
_1179928865.unknown
_1179929075.unknown
_1179928772.unknown
_1179927683.unknown
_1179927828.unknown
_1179923901.unknown
_1179923755.unknown
_1179923787.unknown
_1179923815.unknown
_1179923770.unknown
_1179923724.unknown
_1179923739.unknown
_1179923479.unknown
_1179620790.unknown
_1179922504.unknown
_1179922979.unknown
_1179923112.unknown
_1179923147.unknown
_1179923159.unknown
_1179923133.unknown
_1179923092.unknown
_1179922581.unknown
_1179922873.unknown
_1179922942.unknown
_1179922601.unknown
_1179922562.unknown
_1179922577.unknown
_1179922525.unknown
_1179922176.unknown
_1179922398.unknown
_1179922435.unknown
_1179922499.unknown
_1179922406.unknown
_1179922389.unknown
_1179922375.unknown
_1179922376.unknown
_1179922139.unknown
_1179922155.unknown
_1179922164.unknown
_1179922146.unknown
_1179922046.unknown
_1179921967.unknown
_1179620822.unknown
_1179617315.unknown
_1179618959.unknown
_1179619530.unknown
_1179619588.unknown
_1179619828.unknown
_1179619859.unknown
_1179619803.unknown
_1179619572.unknown
_1179619184.unknown
_1179619390.unknown
_1179618988.unknown
_1179618250.unknown
_1179618616.unknown
_1179618939.unknown
_1179618504.unknown
_1179617557.unknown
_1179617803.unknown
_1179617447.unknown
_1179612966.xls
图表2
-3.30277
-2.7003
-2.2693
-2.01552
-1.77104
-1.58105
-1.44138
-1.3022
-1.19457
-1.07337
lgv与半径r的关系曲线
Sheet1
r v r*r lgv rlgv
0.02 0.000498 0.0004 -3.30277 -0.0660554
0.03 0.001994 0.0009 -2.7003 -0.081009
0.045 0.005379 0.002025 -2.2693 -0.1021185
0.06 0.009649 0.0036 -2.01552 -0.1209312
0.075 0.016942 0.005625 -1.77104 -0.132828
0.09 0.026239 0.0081 -1.58105 -0.1422945
0.105 0.036193 0.011025 -1.44138 -0.1513449
0.12 0.049866 0.0144 -1.3022 -0.156264
0.135 0.06389 0.018225 -1.19457 -0.16126695
0.15 0.084456 0.0225 -1.07337 -0.1610055
0.83 0.295106 0.0868 -18.6515 -1.27511795
0.78 0.292614 0.0855 -12.64843 -1.12805355
Sheet1
气泡径向速度分布v-r
Sheet2
lgv与半径r的关系曲线
Sheet3
_1179614069.unknown
_1179617067.unknown
_1179617229.unknown
_1179616776.unknown
_1179613210.unknown
_1179613811.unknown
_1179614061.unknown
_1179613594.unknown
_1179613170.unknown
_1179612579.unknown
_1179612812.unknown
_1179612899.unknown
_1179612721.unknown
_1179612312.xls
图表1
0.000498
0.001994
0.005379
0.009649
0.016942
0.026239
0.036193
0.049866
0.06389
0.084456
气泡径向速度分布v-r
Sheet1
r v r*r lgv rlgv
0.02 0.000498 0.0004 -3.30277 -0.0660554
0.03 0.001994 0.0009 -2.7003 -0.081009
0.045 0.005379 0.002025 -2.2693 -0.1021185
0.06 0.009649 0.0036 -2.01552 -0.1209312
0.075 0.016942 0.005625 -1.77104 -0.132828
0.09 0.026239 0.0081 -1.58105 -0.1422945
0.105 0.036193 0.011025 -1.44138 -0.1513449
0.12 0.049866 0.0144 -1.3022 -0.156264
0.135 0.06389 0.018225 -1.19457 -0.16126695
0.15 0.084456 0.0225 -1.07337 -0.1610055
0.83 0.295106 0.0868 -18.6515 -1.27511795
0.78 0.292614 0.0855 -12.64843 -1.12805355
Sheet1
气泡径向速度分布v-r
Sheet2
lgv与半径r的关系曲线
Sheet3
_1179612478.unknown
_1179611912.unknown
_1174459050.unknown
_1177962728.unknown
_1178799882.unknown
_1178801057.unknown
_1178802066.unknown
_1178802571.unknown
_1178803661.unknown
_1178905073.unknown
_1178905563.unknown
_1178905662.unknown
_1178905870.unknown
_1179607784.unknown
_1178908640.unknown
_1178905760.unknown
_1178905628.unknown
_1178905300.unknown
_1178905412.unknown
_1178905216.unknown
_1178804414.unknown
_1178904878.unknown
_1178803853.unknown
_1178804371.unknown
_1178803181.unknown
_1178803531.unknown
_1178803593.unknown
_1178803438.unknown
_1178802723.unknown
_1178802925.unknown
_1178802948.unknown
_1178802797.unknown
_1178802620.unknown
_1178802486.unknown
_1178802504.unknown
_1178802074.unknown
_1178801825.unknown
_1178801973.unknown
_1178802009.unknown
_1178801887.unknown
_1178801562.unknown
_1178801801.unknown
_1178801330.unknown
_1178800660.unknown
_1178800885.unknown
_1178800907.unknown
_1178800732.unknown
_1178800464.unknown
_1178800550.unknown
_1178799938.unknown
_1178799133.unknown
_1178799327.unknown
_1178799481.unknown
_1178799842.unknown
_1178799382.unknown
_1178799293.unknown
_1178799305.unknown
_1178799267.unknown
_1178798450.unknown
_1178799089.unknown
_1178799121.unknown
_1178798582.unknown
_1177963397.unknown
_1178798386.unknown
_1177963189.unknown
_1177963368.unknown
_1176762025.unknown
_1177788199.unknown
_1177788388.unknown
_1177958610.unknown
_1177958667.unknown
_1177962461.unknown
_1177958647.unknown
_1177956972.unknown
_1177958470.unknown
_1177790765.unknown
_1177790916.unknown
_1177788848.unknown
_1176817240.unknown
_1177187905.unknown
_1177213731.unknown
_1177213902.unknown
_1177214033.unknown
_1177214290.unknown
_1177217432.unknown
_1177217459.unknown
_1177214332.unknown
_1177214504.unknown
_1177214195.unknown
_1177214262.unknown
_1177214164.unknown
_1177213984.unknown
_1177214018.unknown
_1177213917.unknown
_1177213803.unknown
_1177213821.unknown
_1177213765.unknown
_1177213493.unknown
_1177213632.unknown
_1177213681.unknown
_1177213593.unknown
_1177188024.unknown
_1177188086.unknown
_1177187997.unknown
_1176818025.unknown
_1177187830.unknown
_1177187872.unknown
_1177187493.unknown
_1176817784.unknown
_1176817888.unknown
_1176817525.unknown
_1176813330.unknown
_1176813976.unknown
_1176815824.unknown
_1176816908.unknown
_1176817128.unknown
_1176817158.unknown
_1176816972.unknown
_1176816002.unknown
_1176816368.unknown
_1176816459.unknown
_1176816892.unknown
_1176816469.unknown
_1176816395.unknown
_1176816081.unknown
_1176816157.unknown
_1176816054.unknown
_1176815941.unknown
_1176815966.unknown
_1176815891.unknown
_1176815699.unknown
_1176815741.unknown
_1176815786.unknown
_1176815718.unknown
_1176815572.unknown
_1176815671.unknown
_1176814290.unknown
_1176813618.unknown
_1176813867.unknown
_1176813932.unknown
_1176813770.unknown
_1176813471.unknown
_1176813592.unknown
_1176813449.unknown
_1176812423.unknown
_1176813090.unknown
_1176813152.unknown
_1176813287.unknown
_1176813117.unknown
_1176812744.unknown
_1176812989.unknown
_1176812709.unknown
_1176812215.unknown
_1176812325.unknown
_1176812411.unknown
_1176812295.unknown
_1176812140.unknown
_1176812159.unknown
_1176812108.unknown
_1174500271.unknown
_1176381755.unknown
_1176757626.unknown
_1176757825.unknown
_1176757865.unknown
_1176758087.unknown
_1176758177.unknown
_1176758011.unknown
_1176757844.unknown
_1176757715.unknown
_1176757776.unknown
_1176757687.unknown
_1176381923.unknown
_1176727848.unknown
_1176729107.unknown
_1176731841.unknown
_1176754903.unknown
_1176755903.unknown
_1176729770.unknown
_1176728651.unknown
_1176726426.unknown
_1176726460.unknown
_1176726694.unknown
_1176726408.unknown
_1176381821.unknown
_1176381855.unknown
_1176381795.unknown
_1176380165.unknown
_1176380848.unknown
_1176381134.unknown
_1176381459.unknown
_1176380953.unknown
_1176380795.unknown
_1176378855.unknown
_1176379031.unknown
_1176379509.unknown
_1176379652.unknown
_1176379783.unknown
_1176379909.unknown
_1176379670.unknown
_1176379550.unknown
_1176378980.unknown
_1176378373.unknown
_1176378781.unknown
_1176319904.unknown
_1174495357.unknown
_1174495691.unknown
_1174495891.unknown
_1174500135.unknown
_1174495877.unknown
_1174495449.unknown
_1174495683.unknown
_1174495438.unknown
_1174494503.unknown
_1174495033.unknown
_1174495314.unknown
_1174494756.unknown
_1174493888.unknown
_1174493948.unknown
_1174459249.unknown
_1174479968.unknown
_1174449079.unknown
_1174452201.unknown
_1174458007.unknown
_1174458228.unknown
_1174458722.unknown
_1174458845.unknown
_1174458495.unknown
_1174458084.unknown
_1174458124.unknown
_1174458039.unknown
_1174452413.unknown
_1174457963.unknown
_1174457987.unknown
_1174452457.unknown
_1174452377.unknown
_1174452392.unknown
_1174452320.unknown
_1174449649.unknown
_1174451509.unknown
_1174451758.unknown
_1174451772.unknown
_1174451524.unknown
_1174451462.unknown
_1174451496.unknown
_1174449748.unknown
_1174449360.unknown
_1174449605.unknown
_1174449621.unknown
_1174449404.unknown
_1174449307.unknown
_1174449338.unknown
_1174449238.unknown
_1174449287.unknown
_1174448073.unknown
_1174448682.unknown
_1174448871.unknown
_1174448888.unknown
_1174448715.unknown
_1174448288.unknown
_1174448390.unknown
_1174448121.unknown
_1174192277.unknown
_1174447760.unknown
_1174447926.unknown
_1174447960.unknown
_1174447875.unknown
_1174447202.unknown
_1174447599.unknown
_1174192296.unknown
_1174447121.unknown
_1146978764.unknown
_1146980429.unknown
_1174192051.unknown
_1174192229.unknown
_1174192256.unknown
_1174191631.unknown
_1174191989.unknown
_1174192020.unknown
_1174191948.unknown
_1146981165.unknown
_1146985044.unknown
_1174191610.unknown
_1146981568.unknown
_1146983377.unknown
_1146983413.unknown
_1146981567.unknown
_1146980778.unknown
_1146980922.unknown
_1146980718.unknown
_1146979553.unknown
_1146980168.unknown
_1146980239.unknown
_1146980314.unknown
_1146980225.unknown
_1146979661.unknown
_1146980123.unknown
_1146979594.unknown
_1146979356.unknown
_1146979418.unknown
_1146979473.unknown
_1146979397.unknown
_1146978929.unknown
_1146979238.unknown
_1146978781.unknown
_1146978901.unknown
_1146977934.unknown
_1146978243.unknown
_1146978731.unknown
_1146978076.unknown
_1146977735.unknown
_1146977762.unknown
_1146977719.unknown