油滴沉降速度计算

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); }
油气集输
油滴沉降速度计算
梁静 2014010875 储运 14-2 班
油滴沉降速度计算
油滴在密度为 ρg 气相中所受到的重力为F =
������������ 2 ������ 2
3 ������������ ������
6
������0 − ������������ ������…………4-2
气体对油滴的阻力R = ������������ ������ ������ ������������ …………………………………………………….4-3 4 2 油滴作匀速沉降时 F=R 联立得������������ = 雷诺数������������
查表得������������
4������������������ ������0 − ������������ ������������ = 3������������ ������������
0.5
end
程序： #include<stdio.h> #include<math.h> int 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*1000000000000000000); m1=-cr/1859.95045; m2=-cr/77.47627; m3=-cr/10.02429; a1=pow(2.71828,m1); a2=pow(2.71828,m2); a3=pow(2.71828,m3); c=3.96699*a1+22.56013*a2+96.46966*a3; b=4*9.8*d/1000000*(p1-p2)/3/c/p2; v=pow(b,0.5); printf("油滴沉降速度 %6.4f\n",v); "); "); "); ");
4������������ ������ ������ 0 −������ ������ 3������������ ������ ������ 0.5
………………………………………………………4-4
=
������������ ������ ������ ������ ������ ������ ������
0.5
程序： #include<stdio.h> #include<math.h> int 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("油滴直径 "); ");
2
4.得 Re 与 Ar 数的关系 Re 范围 Re ≤ 2 2 < Re ≤ 500 Re > 500 Ar 范围 Ar ≤ 36 36 < ������������ ≤ 83 × 10 Ar > 83 × 103
3
Re 与 Ar 关系 Re = 0.056Ar Re = 0.153Ar0.714 Re = 1.74Ar0.5
2． 按流态分区计算 流态 层流 过渡流 湍流
2 ��e 范围 Re ≤ 2 2 < Re ≤ 500 500 < Re ≤ 2 × 105 Re > 2 × 105
������ 0 −������ ������
CD 计算式 24Re−1 18.5Re−0.6 0.44 0.1
3
������ ������ ������ ������ 0 −������ ������
……………..4-10
判别过渡区流和湍流临界油滴粒径������2 = 43.5
2 ������ ������
1
3
������ ������ ������ ������ 0 −������ ������
……………..4-11
1.143 0.153������0.714 ������������ ������0 − ������������ ������������ = 0.428 0.286 ������������ ������������
0.714
Y ������ < ������2
N
end
������������������ ������0 − ������������ ������������ = 1.74 ������������
3
……………………………………………………………………………….4-5
1． 按相关式计算
������������ = ������������ + ������������ 0.5 + 0.34………………………………………………4-6
迭代法 (1) 设 CD=0.34，由式（4-4）计算的该液滴的沉降速度 vd； (2) 由求得的 vd 按式（4-5）求 Re； (3) 由 Re 按式（4-6）求 CD; (4) 由 CD 按式（4-4）求，与上一个 vd 进行比较，若在控制误差范围内，计算所得的即 vd 为欲求的沉降速度； (5) 否则，返回步骤（2）直至前后两次求得的 vd（或 CD）在控制误差范围内。 dd,ρ 0,ρ g,g 已知 CD=0.34,vd=0
24
4������������������ ������0 − ������������ ������������ = 3������������ ������������
0.5
Y Vd=v end ������������ − ������ < ������ N
������������ =
24 3 + 0.5 + 0.34 ������������ ������������
������������ =
������������ ������������ ������������ ������������
程序：
#include<stdio.h> #include<math.h> int 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); } "); ");
用阿基米德准数求油滴沉降速度时，先根据欲求沉降速度的油滴直径 dd 和分离条件下的油 气物性参数，求出 Ar，后在表中选择相应的公式算出 Re，再按 Re 定义式（4-5）计算出油 滴沉降速度 vd，这样避免了繁琐的迭代。
根据给定分离条件和油气性质算出 d1 和 d2 后，可根据任一粒径的油滴处于何种流态下降， 并从式（4-7）~（4-9）中选择相应的计算式，计算该油滴的沉降速度。
2 ������������ 1 3
������1 = 3.3
μ g,ρ 0,ρ g,g 已知
������������ ������ ������0 − ������������
}
4． 阿基米德准数法
3 令：Ar = ������������ ������0 − ������������ ������������������ ������2 ������ ,Ar 为无因次数，称阿基米德数。