iron - Ts1(i)=Tp(i)+273; Tf=(Ts1(i)+296)/2; Beta=1/Tf;

Info iconThis preview shows pages 1–2. Sign up to view the full content.

View Full Document Right Arrow Icon
p=2720; cp=875; qdot=120; q %initialize % Tp=296; T %evaluate h L=.200025; g=9.81; Tif=Tp; T Ts1=Tp+273; Tf=(Ts1+296)/2; Beta=1/Tf; Alpha=((6e-14)*Tf^3)-((6e-11)*Tf^2)+((2e-7)*Tf)-3e-5; Kv=((8e-11)*Tf^2)+((4e-8)*Tf)-3e-6; %Kinematic viscocity Tc=((2e-14)*Tf^4)-((1e-10)*Tf^3)+((1e-7)*Tf^2)+((4e-6)*Tf)+0.0137; %Thermal Conductivity (k) Pr=((3e-20)*Tf^6)-((4e-16)*Tf^5)+((1e-12)*Tf^4)-((3e-9)*Tf^3)+((3e-6)*Tf^2)- (0.0013*Tf)+0.8975; Ral=(g*Beta*(Ts1-Tif)*L^3)/(Alpha*Kv); if Ral<=10e9 Nul=0.68+((0.670*(Ral.^(1/4)))/((1+(0.492/Pr).^(9/16)).^(4/9))); else Nul=(0.825+((0.387*(Ral).^(1/6))/((1+(0.492/Pr).^(9/16)).^(8/27))))^2; end h=(Nul*Tc)/L; %heat Transfer coefficient Tp1=tp+2*((qdot-(.0176*h*(tp-296)+.87*5.67*10^(-8)*(tp.^4-296^4)*.0176))/ (p*.000105*cp)) i=2; Tp(2)=Tp1; t=2; while t<60 Tif=Tp(i);
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Background image of page 2
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: Ts1(i)=Tp(i)+273; Tf=(Ts1(i)+296)/2; Beta=1/Tf; Alpha=((6e-14)*Tf^3)-((6e-11)*Tf^2)+((2e-7)*Tf)-3e-5; Kv=((8e-11)*Tf^2)+((4e-8)*Tf)-3e-6; %Kinematic viscocity Tc=((2e-14)*Tf^4)-((1e-10)*Tf^3)+((1e-7)*Tf^2)+((4e-6)*Tf)+0.0137; %Thermal Conductivity (k) Pr=((3e-20)*Tf^6)-((4e-16)*Tf^5)+((1e-12)*Tf^4)-((3e-9)*Tf^3)+((3e-6)*Tf^2)-(0.0013*Tf)+0.8975; Ral=(g*Beta*(Ts1-Tif)*L^3)/(Alpha*Kv); if Ral&lt;=10e9 Nul=0.68+((0.670*(Ral.^(1/4)))/((1+(0.492/Pr).^(9/16)).^(4/9))); else Nul=(0.825+((0.387*(Ral).^(1/6))/((1+(0.492/Pr).^(9/16)).^(8/27))))^2; end h(i)=(Nul*Tc)/L; %heat Transfer coefficient Tp1(i)=tp(i)+2*((qdot-(.0176*h(i)*(tp(i)-296)+.87*5.67*10^(-8)*(tp(i).^4-296^4)*.0176))/(p*.000105*cp)) i=i+1; Tp(i)=Tp1(i-1); t=t+2; end e...
View Full Document

This note was uploaded on 12/14/2009 for the course MAE 4500 taught by Professor El-gizawy during the Spring '06 term at Missouri (Mizzou).

Page1 / 2

iron - Ts1(i)=Tp(i)+273; Tf=(Ts1(i)+296)/2; Beta=1/Tf;

This preview shows document pages 1 - 2. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online