油滴沉降速度计算.docx

油气集输油滴沉降速度计算梁静2014010875储运14-2班油滴沉降速度计算油滴在密度为ρg气相中所受到的重力为…………4-2气体对油滴的阻力…………………………………………………….4-3油滴作匀速沉降时F=R联立得………………………………………………………4-4雷诺数……………………………………………………………………………….4-5按相关式计算………………………………………………4-6迭代法设CD=0.34，由式（4-4）计算的该液滴的沉降速度vd；由求得的vd按式（4-5）求Re；由Re按式（4-6）求CD;由CD按式（4-4）求，与上一个vd进行比较，若在控制误差范围内，计算所得的即vd为欲求的沉降速度；否则，返回步骤（2）直至前后两次求得的vd（或CD）在控制误差范围内。程序：#includestdio.h #includemath.hint main() {double p1,p2,u,d,c,a1,a2,v1,v,c1,re,m;printf(液体密度);scanf(%lf,p1);printf(气体密度);scanf(%lf,p2);printf(气体粘度 );scanf(%lf,u);printf(油滴直径 );scanf(%lf,d);c1=0.34;do{c=c1;a1=4*9.8*d/1000000*(p1-p2)/3/c/r2;v=pow(a1,0.5);re=v*d*p2/u/1000000;m=pow(re,0.5);c1=24/re+3/m+0.34;a2=4*9.8*d/1000000*(p1-p2)/3/c/p2;v1=pow(a2,0.5);c1=c1+0.01;}while(fabs(v1-v)=0.001);printf(油滴沉降速度%8.4f\n,v);}按流态分区计算流态Re范围CD计算式层流过渡流湍流0.440.1层流区斯托克斯公式……………………………………….4-7过渡区阿伦公式……………………….4-8湍流区 CD=0.44 牛顿公式…………………..4-9判别层流和其他流态临界油滴粒径……………..4-10判别过渡区流和湍流临界油滴粒径……………..4-11根据给定分离条件和油气性质算出d1和d2后，可根据任一粒径的油滴处于何种流态下降，并从式（4-7）~（4-9）中选择相应的计算式，计算该油滴的沉降速度。程序：#includestdio.h #includemath.hint main() {double p1,p2,u,d,d1,d2,v,m1,m2,m4,a1,a2,a3,a4,a5,a6;printf(液体密度);scanf(%lf,p1);printf(气体密度);scanf(%lf,p2);printf(气体粘度 );scanf(%lf,u);printf(油滴直径 );scanf(%lf,d);m1=1000000000000000000*u*u/(p2*9.8*(p1-p2));a1=pow(m1,0.333333);d1=3.3*a1;d2=43.5*a1;if(d=d1){v=d*d*9.8*(p1-p2)/18/u;} else{if(d=d2){m2=9.8*d*(p1-p2)/p2;a2=pow(a2,0.5);v=1.74*a2;}else{a3=pow(d/1000000,1.143);m4=p1-p2;a4=pow(m4,0.714);a5=pow(u,0.428);a6=pow(p2,0.286);v=0.153*5.1019*a3*a4/a5/a6;}}printf(油滴沉降速度%8.4f\n,v); }图解法联解式（4-4）（4-5）得……………………………….4-12将图4-6改画成4-7形式，由上式求得后，由图4-7查得，代入式（4-4）求得油滴的沉降终速vd。令,经回归可用以下相关式表示式中系数见下表：a13.96699b11859.95045a222.56013b277.47627a396.46966b310.02429程序：#includestdio.h #includemath.hint main() { double p1,p2,u,d,cr,c,m1,m2,m3,a1,a2,a3,b,v;printf(3.图解法\n请输入所需数据\n);printf(液体密度 );scanf(%lf,p1);printf(气体密度 );scanf(%lf,p2);printf(气体粘度 );scanf(%lf,u);printf(油滴直径 );scanf(%lf,d);cr=4*9.8*d*d*d*p2*(p1-p2)/(3*u*u*100000000000000