소개글
[열전달] Unsteady state, One dimension, Double layer에 대한 자료입니다.
목차
1. System
2. Condition & Solve Equations
3. Result (Graph & Eigenvalues)
4. Coding
5. References
본문내용
beta=zeros; %베타값은 n의 변화에 따라 10개 값이 들어가는 매트릭스로 구성
nu=zeros(20,10); %베타(t,n)에 따라 변하므로 값의 개수가 같다.
gamma=zeros(20,10); %상동
N=zeros(20,10);
p=zeros(20,10);
%syms beta;
k1=0.6; %thermal conductivity of water layer
k2=10; %thermal conductivity of stainless layer
h3=10; %film coefficient
k2star=10; %seperation constant
al1=0.14*10^(-6); %thermal diffusivity of water layer
al2=3.5*10^(-5); %thermal diffusivity of stainless layer
a=0.03; %length of ice cream layer
b=0.035; %length of wheat layer %여기까지 상수 입력
u1=0; %n에 따른 합을 구하기 위해 초기값0 설정
u2=0;
%k1,k2,h3,k2star,al1,al2,a,b는 상수이므로 매트릭스 연산을 할 필요 없음.
H=(b*h3)/k2star; %상수
K=(k1/k2)*((al2/al1)^0.5); %상수
for i=1:1:10;
syms beta1
betax=beta1-((i*pi)+atan(k2star*beta1/(h3*al2^0.5)))/(((a-b)/al2^0.5)-(a/al1^0.5)); %Eigen value 식-미지수로 이루어져있다.
beta1=double(solve(betax));
for t=1:1:20; %time 5,10,15,20,25,30,35,40,45,50
beta(t,i)=beta1;
end
%betax=beta-((2*pi)+atan(k2star*beta/(h3*al2^0.5)))/(((a-b)/al2^0.5)-(a/al1^0.5)); %Eigen value 식-미지수로 이루어져있다.
%beta=double(solve(betax));
%beta
end
beta
for t=1:1:20;
참고문헌
5. References
- Yunus A. Cengel, Heat transfer, McGraw-Hill c2003.
- M.Necati Ozisik, Heat conduction, New York:Wiley c1980.
- Erwin Kreyszig, Kreyszig 공업수학, 범한서적주식회사 2007.
- Lindon C. Thomas, Heat transfer, Prentice-Hall International 1992.