با استفاده از نرم افزار متلب و چند خط برنامه نویسی می توان انتقال حرارت یک فین را که در سیال قرار گرفته محاسبه نمود.برای این فین از طریق جابجایی حرارت به سیال منتقل می شود. کد های زیر را در یک M-File وارد نمایید: ۱٫ 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) تصویری از محیط اجرای کد در متلب: فایل آماده را می توانید از لینک زیر دانلود نمایید: رمز فایل فشرده: www.irchemeng.ir حجم فایل: ۳۴٫۷۶Kb منبع:انجمن مهندسی شیمی
