Untitled
unknown
matlab
7 months ago
972 B
2
Indexable
Never
% Givens L=8; W=1; q=200; T0=20; V0=20*5/18; x=0:0.1:L; delta(1)=0; for i=2:length(x) Tsa=40; % surface temp assumption while true T_f=(Tsa+T0)/2; [rho]=air(T_f,1013.25,[],'rho'); [mu]=air(T_f,1013.25,[],'mu'); [k]=air(T_f,1013.25,[],'k'); [Pr]=air(T_f,1013.25,[],'Pr'); Re=rho*V0*x(i)/mu; if Re>500000 Nu=0.0308*Re^0.8*Pr^(1/3); d=(0.381*Re^(-0.2)-10256*Re^(-1))*x(i); else Nu=0.453*Re^0.5*Pr^(1/3); d=(5*x(i))/sqrt(Re); end H=k*Nu/x(i); T=(q/H)+T0; if T==Tsa break; else Tsa=T; end end Ts(i)=T; h(i)=H; delta(i)=d; end plot(x,delta) xlabel('Cartage length','fontsize',18) ylabel('Boundry layer','fontsize',18) title('Viscous boundary layer against the cartage length'); grid on