소개글
[기계항공공학실험] 온도실험에 대한 자료입니다.
목차
□ 수치해석 리포트
∘ Plot the 1-D temperature profile with analytical solution
(temperature vs fin length)
∘ Derive a discrete equation with the energy balance method and plot the temperature profile with Finite Differential Method(FDM)
(temp. vs (x,y))
∘ Compare '1D temperature profile' with 'temperature profile by FDM' and comment about the results.
□ 온도실험 결과리포트
∘ 실험목적 / 장치 / 방법
∘ DATA & RESULT
- Calibration, Fin의 2차원 온도분포
- analytic solution 의Temperature profile
- Thermocouple로 측정된 Temperature profile
- TLC의 HUE값으로 Calibration한 Temperature profile
- FDM으로 얻어진 data의 Temperature profile
∘ DISCUSSION
- Fin Effectiveness와 Fin Efficiency
□ Matlab 소스코드
본문내용
(4) FDM으로 얻은 data의 Temperature profile
clear all
close all
T_fin = 28.8;
T_b = 34.8;
T_inf = 19;
k = 401;
h = 3.4752;
dx = 0.002;
d = 0.002;
A = eye(3750);
B = zeros(3750,1);
C = zeros(3750,1);
T = zeros(150,50);
for i=1:1:25
B(i,1) = T_fin;
end
for i=3726:1:3750
B(i,1) = T_b;
end
for i=26:1:3725
A(i,i-25) = d;
A(i,i-1) = d;
A(i,i) = -4*d-2*(h*dx*dx/k);
A(i,i+1) = d;
A(i,i+25) = d;
B(i,1) = -2*(h*dx*dx/k)*T_inf;
end
for i=2:1:149
A(25*i-24,25*i-25) = 0;
A(25*i-24,25*i-24) = -2*(h*dx*(d+dx)/k+2*d);
A(25*i-24,25*i-23) = 2*d;
B(25*i-24,1) = -2*(h*dx*(d+dx)/k)*T_inf;
A(25*i,25*i) = -3*d-2*(h*dx*dx/k);
A(25*i,25*i+1) = 0;
end
C = inv(A)*B
for i=1:1:150
for j=1:1:25
T(i,j) = C(25*i+j-25,1);
end
for j=26:1:50
T(i,j) = T(i,51-j);
end
end
surf(T)
title('2D Fin Analysis(T_fin is known)');
xlabel('wide'),ylabel('length'),zlabel('temperature');