gap_x=0.05; %
gap_t=600; %
goal=0.01; %
L=16; %
t=24600; %
Din=0.1; %
Dout=0.2; %
din=0.097; %
dout=0.103; %
Vin=9.6; %
Vout=2.4; %
Vo=4.8; %
B=0.7; %
Hc=0.005; %
S=1.35; %
Hr=0.0008; %
Hins=0.03; %
Ag=pi*((dout/2)^2-(din/2)^2); %
Cp=1.005; %
XIa=0.94; %
XIg=0.02; %
XIr=0.09; %
XIc=0.01; %
TOg=0.9; %
TOc=0.91; %
FAg=0.86; %
FAr=0.03; %
FAc=0.9; %
JUr=0.98; %
DAg=1.25; %
DAr=177; %
DAc=1.25; %
FAins=0.03; %
DAins=0.034; %
Cb=5.67*10^(-8); %
K=0.4; %
Tamb=xlsread('.xls','','c2:c43'); %
I=xlsread('.xls','','d2:d43'); %
vfeng=xlsread('.xls','','e2:e43'); %
Tin=xlsread('.xls','','f2:f43'); %
Qa1=I*B*gap_x*TOc*JUr*K*TOg*XIa; %(1)
Qg1=I*B*gap_x*TOc*JUr*K*XIg; %(1)
Qa2=I*B*gap_x*TOc*JUr*(1-K)*TOg*XIa; %(2)
Qg2=I*B*gap_x*TOc*JUr*(1-K)*XIg; %(2)
Qc=I*B*gap_x*XIc; %
Qr=I*B*gap_x*TOc*XIr; %
To=0.0552*((Tamb+273.15).^1.5); %
Nub=round(L/gap_x)+1; %
Nubt=round(t/gap_t)+1; %
tout=ones(Nubt,Nub); %
t1=ones(Nubt,Nub); %(1)
t2=ones(Nubt,Nub); %(2)
tg1=ones(Nubt,Nub); %1
tg2=ones(Nubt,Nub); %2
tr=ones(Nubt,Nub); %
tc=ones(Nubt,Nub); %
tins=ones(Nubt,Nub); %
Pouf1=9.55*10^(-6)*t1.^2-0.0044*t1+1.29; %(1)
Pr1=-0.000185*t1+0.7073; %(1)Pr
Pouf2=9.55*10^(-6)*t2.^2-0.0044*t2+1.29; %(2)
Pr2=-0.000185*t2+0.7073; %(2)Pr
ND1=(0.1026*t1+13.04)*10^(-6); %(1)
ND2=(0.1026*t2+13.04)*10^(-6); %(2)
LM1=(0.00754*t1+2.443) *10^(-2); %(1)
LM2=(0.00754*t2+2.443) *10^(-2); %(2)
Nu1=0.0214*((Vo*din./ND1).^0.8-100).*(Pr1.^0.4).*(((t1+273.15)./(tg1+273.15)).^0.45)*(1+(din/L)^(2/3));
Nu2=0.0214*((Vo*din./ND2).^0.8-100).*(Pr2.^0.4).*(((t2+273.15)./(tg2+273.15)).^0.45)*(1+(din/L)^(2/3));
hf1=Nu1.*LM1/din; %(1)
hf2=Nu2.*LM2/din; %(2)
hcgr1=3.25+0.0085*(tg1-tr)/(2*dout); %(1)
hcgr2=3.25+0.0085*(tg2-tr)/(2*dout); %(2)
hcgc1=3.25+0.0085*(tg1-tc)/(2*dout); %(1)
hcgc2=3.25+0.0085*(tg2-tc)/(2*dout); %(2)
hrgr1=Cb*(((tg1+273.15).^2+(tr+273.15).^2).*((tg1+273.15)+(tr+273.15)))/(1/FAr-1+(1/FAg-1)*S/pi/dout); %(1)
hrgr2=Cb*(((tg2+273.15).^2+(tr+273.15).^2).*((tg2+273.15)+(tr+273.15)))/(1/FAr-1+(1/FAg-1)*S/pi/dout); %(2)
hrgc1=Cb*(((tg1+273.15).^2+(tc+273.15).^2).*((tg1+273.15)+(tc+273.15)))/(1/FAc+(1/FAg-1)*B/pi/dout); %(1)
hrgc2=Cb*(((tg2+273.15).^2+(tc+273.15).^2).*((tg2+273.15)+(tc+273.15)))/(1/FAc+(1/FAg-1)*B/pi/dout); %(2)
hcco=5.7+3.8*vfeng; %
hcro=5.7+3.8*vfeng; %
hrrc=Cb*((tc+273.15).^2+(tr+273.15).^2).*((tc+273.15)+(tr+273.15))/(1/FAc-1+(1/FAr-1)*B/S); %
tout(:,1)=Tin; %(1)
t1(:,1)=Tin(:,1); %(1)
t2(:,1)=Tin(:,1); %(2)
tg1(:,1)=(Qg1(:,1)+pi*din*gap_x*hf1(:,1).*t1(:,1)+S*gap_x*hcgr1(:,1).*tr(:,1)+S*gap_x*hrgr1(:,1).*tr(:,1)+B*gap_x*hcgc1(:,1).*tc(:,1)+B*gap_x*hrgc1(:,1).*tc(:,1)+Ag*DAg/gap_x*tg1(:,2))./(pi*din*gap_x*hf1(:,1)+S*gap_x*hcgr1(:,1)+S*gap_x*hrgr1(:,1)+B*gap_x*hcgc1(:,1)+B*gap_x*hrgc1(:,1)+Ag*DAg/gap_x); %1
tg2(:,1)=(Qg2(:,1)+pi*din*gap_x*hf2(:,1).*t2(:,1)+S*gap_x*hcgr2(:,1).*tr(:,1)+S*gap_x*hrgr2(:,1).*tr(:,1)+B*gap_x*hcgc2(:,1).*tc(:,1)+B*gap_x*hrgc2(:,1).*tc(:,1)+Ag*DAg/gap_x*tg2(:,2))./(pi*din*gap_x*hf2(:,1)+S*gap_x*hcgr2(:,1)+S*gap_x*hrgr2(:,1)+B*gap_x*hcgc2(:,1)+B*gap_x*hrgc2(:,1)+Ag*DAg/gap_x); %2
tc(:,1)=(Qc(:,1)+B*gap_x*hcco(:,1).*Tamb(:,1)+B*gap_x*FAc*Cb*(To(:,1)-273.15).*((tc(:,1)+273.15).^2+To(:,1).^2).*((tc(:,1)+273.15)+To(:,1))+5*B*gap_x*tr(:,1)+B*gap_x*hrrc(:,1).*tr(:,1)+B*gap_x*hcgc1(:,1).*tg1(:,1)+B*gap_x*hrgc1(:,1).*tg1(:,1)+B*gap_x*hcgc2(:,1).*tg2(:,1)+B*gap_x*hrgc2(:,1).*tg2(:,1)+B*Hc*DAc/gap_x*tc(:,2))./(B*gap_x*hcco(:,1)+B*gap_x*FAc*Cb*((tc(:,1)+273.15).^2+To(:,1).^2).*((tc(:,1)+273.15)+To(:,1))+5*B*gap_x+B*gap_x*hrrc(:,1)+B*gap_x*hcgc1(:,1)+B*gap_x*hrgc1(:,1)+B*gap_x*hcgc2(:,1)+B*gap_x*hrgc2(:,1)+B*Hc*DAc/gap_x); %
tr(:,1)=(Qr(:,1)+S*gap_x*hcro(:,1).*Tamb(:,1)+S*gap_x*FAr*Cb*(To(:,1)-273.15).*((tr(:,1)+273.15).^2+To(:,1).^2).*((tr(:,1)+273.15)+To(:,1))+5*B*gap_x*tc(:,1)+B*gap_x*hrrc(:,1).*tc(:,1)+S*gap_x*hcgr1(:,1).*tg1(:,1)+S*gap_x*hrgr1(:,1).*tg1(:,1)+S*gap_x*hcgr2(:,1).*tg2(:,1)+S*gap_x*hrgr2(:,1).*tg2(:,1)+S*Hr*DAr/gap_x*tr(:,2))./(S*gap_x*hcro(:,1)+S*gap_x*FAr*Cb*((tr(:,1)+273.15).^2+To(:,1).^2).*((tr(:,1)+273.15)+To(:,1))+5*B*gap_x+B*gap_x*hrrc(:,1)+S*gap_x*hcgr1(:,1)+S*gap_x*hrgr1(:,1)+S*gap_x*hcgr2(:,1)+S*gap_x*hrgr2(:,1)+S*Hr*DAr/gap_x); %
times = 0;
J=1;
while J==1
if times == 100
break;
end
times = times + 1;
Tout=tout; %
T1=t1; %1
T2=t2; %2
Tg1=tg1; %1
Tg2=tg2; %2
Tr=tr; %
Tc=tc; %
Tins=tins; %
for j=1:1:Nubt %
for i=2:1:Nub-1 %
T1(j,i)=(Qa1(j,1)-0.5*Cp*((din/2)^2)*pi*Vo*Pouf1(j,i+1).*T1(j,i+1)+0.5*Cp*((din/2)^2)*pi*Vo*Pouf1(j,i-1).*T1(j,i-1))./(pi*din*hf1(j,i)*gap_x)+Tg1(j,i); %1
end
end
for j=1:1:Nubt %
for i=2:1:Nub-1 %
T2(j,i)=(Qa2(j,1)-0.5*Cp*((din/2)^2)*pi*Vo*Pouf2(j,i+1).*T2(j,i+1)+0.5*Cp*((din/2)^2)*pi*Vo*Pouf2(j,i-1).*T2(j,i-1))./(pi*din*hf2(j,i)*gap_x)+Tg2(j,i); %2
end
end
for j=1:1:Nubt %
for i=2:1:Nub-1 %
Tg1(j,i)=(Qg1(j,1)+pi*din*gap_x*hf1(j,i).*T1(j,i)+S*gap_x*hcgr1(j,i).*Tr(j,i)+S*gap_x*hrgr1(j,i).*Tr(j,i)+B*gap_x*hcgc1(j,i).*Tc(j,i)+B*gap_x*hrgc1(j,i).*Tc(j,i)+Ag*DAg/gap_x*(Tg1(j,i-1)+Tg1(j,i+1)))./(pi*din*gap_x*hf1(j,i)+S*gap_x*hcgr1(j,i)+S*gap_x*hrgr1(j,i)+B*gap_x*hcgc1(j,i)+B*gap_x*hrgc1(j,i)+2*Ag*DAg/gap_x); %1
end
end
for j=1:1:Nubt %
for i=2:1:Nub-1 %
Tg2(j,i)=(Qg2(j,1)+pi*din*gap_x*hf2(j,i).*T2(j,i)+S*gap_x*hcgr2(j,i).*Tr(j,i)+S*gap_x*hrgr2(j,i).*Tr(j,i)+B*gap_x*hcgc2(j,i).*Tc(j,i)+B*gap_x*hrgc2(j,i).*Tc(j,i)+Ag*DAg/gap_x*(Tg2(j,i-1)+Tg2(j,i+1)))./(pi*din*gap_x*hf2(j,i)+S*gap_x*hcgr2(j,i)+S*gap_x*hrgr2(j,i)+B*gap_x*hcgc2(j,i)+B*gap_x*hrgc2(j,i)+2*Ag*DAg/gap_x); %1
end
end
for j=1:1:Nubt %
for i=2:1:Nub-1 %
Tc(j,i)=(Qc(j,1)+B*gap_x*hcco(j,1).*Tamb(j,1)+B*gap_x*FAc*Cb*(To(j,1)-273.15).*((Tc(j,i)+273.15).^2+To(j,1).^2).*((Tc(j,i)+273.15)+To(j,1))+5*B*gap_x*Tr(j,i)+B*gap_x*hrrc(j,i)*Tr(j,i)+B*gap_x*hcgc1(j,i).*Tg1(j,i)+B*gap_x*hrgc1(j,i).*Tg1(j,i)+B*gap_x*hcgc2(j,i).*Tg2(j,i)+B*gap_x*hrgc2(j,i).*Tg2(j,i)+B*Hc*DAc/gap_x*(Tc(j,i-1)+Tc(j,i+1)))/(B*gap_x*hcco(j,1)+B*gap_x*FAc*Cb*((Tc(j,i)+273.15).^2+To(j,1).^2).*((Tc(j,i)+273.15)+To(j,1))+5*B*gap_x+B*gap_x*hrrc(j,i)+B*gap_x*hcgc1(j,i)+B*gap_x*hrgc1(j,i)+B*gap_x*hcgc2(j,i)+B*gap_x*hrgc2(j,i)+2*B*Hc*DAc/gap_x); %
end
end
for j=1:1:Nubt %
for i=2:1:Nub-1 %
Tr(j,i)=(Qr(j,1)+S*gap_x*hcro(j,1).*Tamb(j,1)+S*gap_x*FAr*Cb*(To(j,1)-273.15).*((Tr(j,i)+273.15).^2+To(j,1).^2).*((Tr(j,i)+273.15)+To(j,1))+5*B*gap_x*Tc(j,i)+B*gap_x*hrrc(j,i)*Tc(j,i)+S*gap_x*hcgr1(j,i).*Tg1(j,i)+S*gap_x*hrgr1(j,i).*Tg1(j,i)+S*gap_x*hcgr2(j,i).*Tg2(j,i)+S*gap_x*hrgr2(j,i).*Tg2(j,i)+S*Hr*DAr/gap_x*(Tr(j,i-1)+Tr(j,i+1)))/(S*gap_x*hcro(j,1)+S*gap_x*FAr*Cb*((Tr(j,i)+273.15).^2+To(j,1).^2).*((Tr(j,i)+273.15)+To(j,1))+5*B*gap_x+B*gap_x*hrrc(j,i)+S*gap_x*hcgr1(j,i)+S*gap_x*hrgr1(j,i)+S*gap_x*hcgr2(j,i)+S*gap_x*hrgr2(j,i)+2*S*Hr*DAr/gap_x); %
end
end
if norm(abs(T1-t1),1)>goal norm(abs(T1-t1),1)>goal %
tout=Tout; %
t1=T1; %1
t2=T2; %2
tg1=Tg1; %1
tg2=Tg2; %2
tr=Tr; %
tc=Tc; %
tins=Tins; %
J=1;
else
J=0;
end
end
Tout=(T1+T2)/2;
xlswrite('.xls',Tout(:,Nub-1),'sheet2','c2');