//shell and tube condenser simulation
//Wang Bingming compile
//1-2型管殼式換熱器,管子排列方式采用正三角形錯排,折流板厚度為5mm,換熱管的材料采用銅


#include <stdio.h>
#include <math.h>
main()
{
//冷凝器的基本尺寸	
printf("1 冷凝器的基本尺寸\n\n");
int Ls,d0,n0;//Ls殼體長度,d0換熱管外徑,n0換熱管子總數(shù)
float Ds,a,fai,derta0,L,t,b,f1,Mco2,rouco2,uco2;//Ds殼體內(nèi)徑,a換熱管排列修正系數(shù),fai管板充滿系數(shù),derta0換熱管厚度,L換熱管長,t換熱管中心距,b折流板間距
                                    //f1管側(cè)每程的通道面積,Mco2 CO2質(zhì)量流量,rouco2 CO2密度,uco2 CO2換熱管內(nèi)流速
printf("換熱管排列修正系數(shù)a為 ");
scanf("%f",&a);
printf("管板充滿系數(shù)fai為 ");
scanf("%f",&fai);
printf("換熱管厚度derta0為 ");
scanf("%f",&derta0);
printf("換熱管長L為 ");
scanf("%f",&L);
printf("換熱管外徑d0為 ");
scanf("%d",&d0);
printf("CO2質(zhì)量流量Mco2為 ");
scanf("%f",&Mco2);
printf("CO2的密度rouco2為 ");
scanf("%f",&rouco2);
printf("CO2管內(nèi)流速uco2為 ");
scanf("%f",&uco2);
f1=Mco2/(rouco2*uco2);
n0=2*4*f1/(3.14*(d0-2*derta0)*(d0-2*derta0))+1;
printf("換熱管子總數(shù)為%d\n",n0);
uco2=2*4*Mco2/(3.14*(d0-2*derta0)*(d0-2*derta0)*n0);
printf("CO2的實際管內(nèi)流速uco2為 %f\n",uco2);
t=int(1.25*d0);
Ds=int(a*t*sqrt(n0/fai))+1;
Ls=5*Ds;
b=0.5*Ds;
printf("殼體內(nèi)徑Ds為%g\n殼體長度Ls為%d\n換熱管中心距t為%g\n折流板間距b為%g\n\n\n",Ds,Ls,t,b);


//冷卻介質(zhì)的進口溫度以及其溫升的選擇
//此處冷卻介質(zhì)為氨
printf("2 冷卻介質(zhì)的進口溫度以及其溫升的選擇\n\n");
int drtat;//drtat為冷卻介質(zhì)的溫升
float t11,t111;//t11為冷卻介質(zhì)的進口溫度,t111為冷卻介質(zhì)出口溫度
printf("冷卻介質(zhì)的溫升drtat為 ");
scanf("%d",&drtat);
printf("冷卻介質(zhì)的進口溫度t11為 ");
scanf("%f",&t11);
t111=t11+drtat;
printf("冷卻介質(zhì)出口溫度t111為%g\n\n\n",t111);


//冷凝器中的污垢熱阻
printf("3 冷凝器中污垢的熱阻\n\n");
float gama0,gamai;
printf("管外熱阻gama0為 ");
scanf("%f",&gama0);
printf("管內(nèi)熱阻gamai為 ");
scanf("%f",&gamai);
printf("\n\n");


//冷凝器的設(shè)計計算
printf("4 冷凝器的設(shè)計計算\n\n");

}