气液两相流雾环状流动的压降计算
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
scanf("%lf,%lf",&M1,&M11);
printf("请分别输入两相流的液相和气相的密度(kg/m^3):rou1,rou11\n");
scanf("%lf,%lf",&rou1,&rou11);
printf("请分别输入两相流的液相和气相的动力粘度（Pa*s）:u1,u11\n");
double lmdl,lmdg,Qg2,Qg21,QlE2,RelE,lmdlE,dpfdzlE;
double WoE1,a,b,D,A,dpfdz,dpfdz1,Y,E;
printf("请输入垂直上升管的管径D(m),xgm\n");
scanf("%lf,%lf",&D,&xgm);
printf("请分别输入两相流的液相和气相的质量流量(kg/s):M1,M11\n");
lmdlE=0.3164*pow(RelE,-0.25);
WoE1=M1*(1-E)/(rou1*A);
dpfdzlE=lmdlE*rou1*WoE1*WoE1/(2*D);
QlE2=dpfdz1/dpfdzlE;
a=1-pow(QlE2,-0.5);
Qg2=(1+75*(1-a))*pow(a,-2.5)*(M11+E*M1)/M11*pow(1-2*a/(1-a)*rou11/rou1*M1*(1-E)/M11,2);
dpfdz=Qg2*dpfdzg;
}
b=(1-pow(a,0.5))*D/2;
printf("摩擦压降dpf/dz：%lf pa/m\n",dpfdz);
printf("截面含气率a：%lf\n",a);
printf("平均液膜厚度b：%lf m\n",b);
}
五、计算算例及结果
气-水混合物向wk.baidu.com流过内径 的垂直管，物性参数如下：
dpfdz=Qg2*dpfdzg;
while(fabs(dpfdz-dpfdz1)>=0.01)
{
dpfdz1=dpfdz;
Qg2=dpfdz/dpfdzlE;
a=1-pow(Qg2,-0.5);
Qg2=(1+75*(1-a))*pow(a,-2.5)*(M11+E*M1)/M11*pow(1-2*a/(1-a)*rou11/rou1*M1*(1-E)/M11,2);
1.预先估计摩阻梯度值
先计算气相和液相的折算速度 和
据此求出雷诺数
根据雷诺数判断流形，再按布拉修斯公式求出摩阻系数 和 ，即
然后用公式计算出分液相及分气相的摩阻梯度 和
再按下式求出马蒂内里参数
按奇斯霍姆式求出分气相折算系数
最后估算出两相压降梯度
2.求窜流比值
计算出参数
然后根据窜流比 的关系计算出 ，即
{
doublelmd;
if(Re<=1000)
{
lmd=64/Re;
}
else
{
lmd=0.3164*pow(Re,-0.25);
}
returnlmd;
}
main()
{
double M1,M11,rou1,rou11,u1,u11,xgm;
double jf,jg,Rel,Reg,dpfdzl,dpfdzg,X2;
3.计算气流夹带水滴工况下的分液相折算系数
先算出计入窜流量的雷诺数
据此算出摩阻系数
算出气流夹带水滴时的分液相摩阻梯度
式中， 为扣除被气流带走液滴后的液相折算速度，即
然后算出计入窜流量的分液相折算系数
4.按下式算出雾环流状的截面含气率，即
5.计算两相摩阻梯度
先按公式算出 ，即
再据此求得
将所得结果与步骤1的预估值进行比较，若两者相等或其差值在容许误差范围之内，则可认为计算结束。否则，将所得的 值带入步骤3中的公式
# include <math.h>
# define pi 3.1415926
double CL1(double Y)
{
double E;
if(Y<=4)
{
E=0.005515*pow(Y,2.858);
}
else
{
E=0.400038*pow(Y-4,0.2875);
}
return E;
}
double CL2(double Re)
求：摩擦压降 、截面含气率
解：计算过程和结果如下：
请输入垂直上升管的管径D(m),xgm
0.03,0.072
请分别输入两相流的液相和气相的质量流量(kg/s):M1,M11
0.2,0.1
请分别输入两相流的液相和气相的密度(kg/m^3):rou1,rou11
1000,1.64
请分别输入两相流的液相和气相的动力粘度（Pa*s）:u1,u11
Harbin Institute of Technology
课程设计说明书（论文）
课程名称：气液两相流
设计题目：雾环状流动的压降计算
院系：能源学院
班级：0902301
设计者：王小亮
学号：**********
*******
设计时间：2012年6月10日
哈尔滨工业大学
气液两相流雾环状流动的压降计算
一、计算步骤
scanf("%lf,%lf",&u1,&u11);
A=pi*D*D/4;
jf=M1/(rou1*A);
jg=M11/(rou11*A);
Rel=rou1*jf*D/u1;
Reg=rou11*jg*D/u11;
lmdl=CL2(Rel);
lmdg=CL2(Reg);
dpfdzl=lmdl*rou1*jf*jf/(2*D);
重新计算 ，再依次经过步骤4、5算出新的 。如此重复迭代，直到最后的压降梯度值与其前一次算出压降梯度值的差值的绝对值满足所需精度（0.01）为止。然后用最后所得的截面含气率 ，通过公式
计算出平均夜膜厚度 。
二、程序流程图
三、变量名说明
分液相折算系数（ ）；
四、计算程序
# include <stdio.h>
0.001,0.000018
摩擦压降dpf/dz：15127.046754 pa/m
截面含气率a：0.977168
平均液膜厚度b：0.000172 m
dpfdzg=lmdg*rou11*jg*jg/(2*D);
X2=dpfdzl/dpfdzg;
Qg21=1+20*pow(X2,0.5)+X2;
dpfdz1=Qg21*dpfdzg;
Y=jg*u11/xgm*pow(rou11/rou1,0.5)*10000;
E=CL1(Y);
RelE=Rel*(1-E);
printf("请分别输入两相流的液相和气相的密度(kg/m^3):rou1,rou11\n");
scanf("%lf,%lf",&rou1,&rou11);
printf("请分别输入两相流的液相和气相的动力粘度（Pa*s）:u1,u11\n");
double lmdl,lmdg,Qg2,Qg21,QlE2,RelE,lmdlE,dpfdzlE;
double WoE1,a,b,D,A,dpfdz,dpfdz1,Y,E;
printf("请输入垂直上升管的管径D(m),xgm\n");
scanf("%lf,%lf",&D,&xgm);
printf("请分别输入两相流的液相和气相的质量流量(kg/s):M1,M11\n");
lmdlE=0.3164*pow(RelE,-0.25);
WoE1=M1*(1-E)/(rou1*A);
dpfdzlE=lmdlE*rou1*WoE1*WoE1/(2*D);
QlE2=dpfdz1/dpfdzlE;
a=1-pow(QlE2,-0.5);
Qg2=(1+75*(1-a))*pow(a,-2.5)*(M11+E*M1)/M11*pow(1-2*a/(1-a)*rou11/rou1*M1*(1-E)/M11,2);
dpfdz=Qg2*dpfdzg;
}
b=(1-pow(a,0.5))*D/2;
printf("摩擦压降dpf/dz：%lf pa/m\n",dpfdz);
printf("截面含气率a：%lf\n",a);
printf("平均液膜厚度b：%lf m\n",b);
}
五、计算算例及结果
气-水混合物向wk.baidu.com流过内径 的垂直管，物性参数如下：
dpfdz=Qg2*dpfdzg;
while(fabs(dpfdz-dpfdz1)>=0.01)
{
dpfdz1=dpfdz;
Qg2=dpfdz/dpfdzlE;
a=1-pow(Qg2,-0.5);
Qg2=(1+75*(1-a))*pow(a,-2.5)*(M11+E*M1)/M11*pow(1-2*a/(1-a)*rou11/rou1*M1*(1-E)/M11,2);
1.预先估计摩阻梯度值
先计算气相和液相的折算速度 和
据此求出雷诺数
根据雷诺数判断流形，再按布拉修斯公式求出摩阻系数 和 ，即
然后用公式计算出分液相及分气相的摩阻梯度 和
再按下式求出马蒂内里参数
按奇斯霍姆式求出分气相折算系数
最后估算出两相压降梯度
2.求窜流比值
计算出参数
然后根据窜流比 的关系计算出 ，即
{
doublelmd;
if(Re<=1000)
{
lmd=64/Re;
}
else
{
lmd=0.3164*pow(Re,-0.25);
}
returnlmd;
}
main()
{
double M1,M11,rou1,rou11,u1,u11,xgm;
double jf,jg,Rel,Reg,dpfdzl,dpfdzg,X2;
3.计算气流夹带水滴工况下的分液相折算系数
先算出计入窜流量的雷诺数
据此算出摩阻系数
算出气流夹带水滴时的分液相摩阻梯度
式中， 为扣除被气流带走液滴后的液相折算速度，即
然后算出计入窜流量的分液相折算系数
4.按下式算出雾环流状的截面含气率，即
5.计算两相摩阻梯度
先按公式算出 ，即
再据此求得
将所得结果与步骤1的预估值进行比较，若两者相等或其差值在容许误差范围之内，则可认为计算结束。否则，将所得的 值带入步骤3中的公式
# include <math.h>
# define pi 3.1415926
double CL1(double Y)
{
double E;
if(Y<=4)
{
E=0.005515*pow(Y,2.858);
}
else
{
E=0.400038*pow(Y-4,0.2875);
}
return E;
}
double CL2(double Re)
求：摩擦压降 、截面含气率
解：计算过程和结果如下：
请输入垂直上升管的管径D(m),xgm
0.03,0.072
请分别输入两相流的液相和气相的质量流量(kg/s):M1,M11
0.2,0.1
请分别输入两相流的液相和气相的密度(kg/m^3):rou1,rou11
1000,1.64
请分别输入两相流的液相和气相的动力粘度（Pa*s）:u1,u11
Harbin Institute of Technology
课程设计说明书（论文）
课程名称：气液两相流
设计题目：雾环状流动的压降计算
院系：能源学院
班级：0902301
设计者：王小亮
学号：**********
*******
设计时间：2012年6月10日
哈尔滨工业大学
气液两相流雾环状流动的压降计算
一、计算步骤
scanf("%lf,%lf",&u1,&u11);
A=pi*D*D/4;
jf=M1/(rou1*A);
jg=M11/(rou11*A);
Rel=rou1*jf*D/u1;
Reg=rou11*jg*D/u11;
lmdl=CL2(Rel);
lmdg=CL2(Reg);
dpfdzl=lmdl*rou1*jf*jf/(2*D);
重新计算 ，再依次经过步骤4、5算出新的 。如此重复迭代，直到最后的压降梯度值与其前一次算出压降梯度值的差值的绝对值满足所需精度（0.01）为止。然后用最后所得的截面含气率 ，通过公式
计算出平均夜膜厚度 。
二、程序流程图
三、变量名说明
分液相折算系数（ ）；
四、计算程序
# include <stdio.h>
0.001,0.000018
摩擦压降dpf/dz：15127.046754 pa/m
截面含气率a：0.977168
平均液膜厚度b：0.000172 m
dpfdzg=lmdg*rou11*jg*jg/(2*D);
X2=dpfdzl/dpfdzg;
Qg21=1+20*pow(X2,0.5)+X2;
dpfdz1=Qg21*dpfdzg;
Y=jg*u11/xgm*pow(rou11/rou1,0.5)*10000;
E=CL1(Y);
RelE=Rel*(1-E);