انتقال حرارت یک بعدی در فین

[فاطمه یاس]

با استفاده از نرم افزار متلب و چند خط برنامه نویسی می توان انتقال حرارت یک فین را که در سیال قرار گرفته محاسبه نمود.برای این فین از طریق جابجایی حرارت به سیال منتقل می شود.
کد های زیر را در یک M-File وارد نمایید:

[LEFT]۱٫   clear all
۲٫   clc
۳٫   %input data
۴٫   l=input(‘Enter length of fin: ‘);
۵٫   n=input(‘Enter number of nod: ‘);
۶٫   r=input(‘Enter radius of fin: ‘);
۷٫   T_b=input(‘Enter temperature of base: ‘);
۸٫   T_f=input(‘Enter temperature of fluid: ‘);
۹٫   k=input(‘Enter thermal conductivity of fin: ‘);
۱۰٫   h=input(‘Enter convection coeficient of fluid: ‘);
۱۱٫   %calculation parameter
۱۲٫   dx=l/(n-1);
۱۳٫   %calculate temperature distribution
۱۴٫   T=zeros(1,n);
۱۵٫   T(1)=T_b;
۱۶٫   iteration=0;
۱۷٫   R=1;
۱۸٫   while R~=0
۱۹٫       TT=T;
۲۰٫       for i=2:n-1
۲۱٫           T(i)=(T(i-1)+T(i+1)+(2*h*T_f*dx^2)/(k*r))/(2+(2*h*dx^2)/(k*r));
۲۲٫       end
۲۳٫       T(n)=(T(n-1)+((h*dx)/(k)+(h*dx^2)/(k*r))*T_f)/(1+(h*dx)/(k)+(h*dx^2)/(k*r));
۲۴٫       R=sum(abs(T-TT));
۲۵٫       iteration=iteration+1;
۲۶٫   end
۲۷٫   %output
۲۸٫   iteration
۲۹٫   T
۳۰٫   plot(T)[/LEFT]

تصویری از محیط اجرای کد در متلب:

​