열전달 예제 4.3(Heat transfer incropera_EXAMPLE 4.3) 매틀랩 코드

소개글

ex 4,3) A major objective in advancing gas turbine engine technologies is to increase the temperature
limit associated with operation of the gas turbine blades. ~ At what location is the temperature a maximum?

예제 풀이에 나와있는 식으로 만들었습니다. 매틀랩에 그대로 붙여넣기하면 나옵니다.
파일 두번째 줄에 dx값만 1, 1/2, 1/4 …로 바꿔주면 됩니다.

목차

없음

본문내용

%행렬크기(m)
dx=1/8;
n=1/dx;
m=(5*n+1)*(2*n+1)+n*(2*n+1);
A=99*eye(m);

%%%%% 1. 첫번째줄(대류)-완성!
for i=1:5*n+1
A(i,i)=-2*(2+0.04*dx); %%첫번째 박스 대각성분
A(i,i+1)=1; %%대각성분의 오른쪽
A(i+1,i)=1; %%대각성분의 왼쪽
A(i,i+(5*n+1))=2; %%대각성분 오른쪽 멀리
end

A(1,1)=-(2+0.04*dx); %%(1,1)성분
A(5*n+1,5*n+1)=-(2+0.04*dx); %%(6,6)성분
A(1,5*n+2)=1; %%(1,7)성분
A(5*n+1,2*(5*n+1))=1; %%(6,12)성분
A(5*n+1,5*n+2)=0; %%(6,7)성분

<중 략>

%%%%% 3. 세번째줄,점13~점15사이에 생기는 점-완성!
w=2*n-1;
for i=1:2*n
A((w+1)*(5*n+1)+i,(w+1)*(5*n+1)+i)=-4; %-4 %세번째박스 대각성분
A((w+1)*(5*n+1)+i,(w+1)*(5*n+1)-1+i)=1; %1 %대각성분 왼쪽
A((w+1)*(5*n+1)+i,(w+1)*(5*n+1)+1+i)=1; %1 %대각성분 오른쪽
A((w+1)*(5*n+1)+i,w*(5*n+1)+i)=1; % %대각성분 왼쪽멀리
A((w+1)*(5*n+1)+i,(w+2)*(5*n+1)+i)=1; %1 %대각성분 오른쪽 멀리
end

A((w+1)*(5*n+1)+1,(w+1)*(5*n+1))=0; %%(13,12)
A((w+1)*(5*n+1)+1,(w+1)*(5*n+1)+2)=2; %%(13,14)

%세번째줄, 15점부터 끝까지-완
w=2*n-1;
for i=2*n+1:5*n+1
A((w+1)*(5*n+1)+i,(w+1)*(5*n+1)+i)=-2*(2+0.008*dx); %네번째 박스 대각성분
A((w+1)*(5*n+1)+i,(w+1)*(5*n+1)+1+i)=1; %%대각성분 오른쪽

참고 자료

없음