[화공생명공학 요소실험] Residence Time Distribution (R T D) RTD 메인

소개글

[화공생명공학 요소실험] Residence Time Distribution (R T D) RTD 메인에 대한 보고서 자료입니다.

목차

Ⅰ.Data & Result
　
　1. 실험 Data
　　ⓐ 4월 14일
　　ⓑ 4월 27일　
　2. CSTR
　　ⓐ 200 RPM
　　ⓑ 700 RPM
　　ⓒ 1500 RPM
　
　3. PBR
　　ⓐ 100%
　　ⓑ 70%
　　ⓒ 40%
　　ⓓ 값

Ⅱ. Discussion

Ⅲ. Reference

Ⅳ. Matlab source
　예시- 4월 14일 데이터 기준

본문내용

);
end
end
end
Tprime(n,n-1)=beta3prime(n-1);
Tprime(n,n)=beta3prime(n-1);
% 행렬 S prime2
for i=1:n
for j=1:n-2
if i==j+1
Sprime2(i,j)=beta1prime2(i);
end
end
end
for j=1:n-2
Sprime2(n,j)=0;
end
% 행렬 prime, prime2
A=Sprime*inv(N)*R+Tprime;
B=Sprime2*inv(N)*R;
% Boundary Condition
BC=[zeros(1,n-3) (1/6)*( h(n-1) - (3*(1-k)^2) / h(n-1) ) ]*inv(N)*R + [zeros(1,n-2) -1/h(n-1) 1/h(n-1)];
d1=BC(1);
dr=BC(2:n-1);
dn=BC(n);
m1=-d1/dn;
m2=-dr/dn;
% coefficient of ODE
M=-u*A+D*B;
M1=M(2:n-1,1);
Mr=M(2:n-1,2:n-1);
Ml=M(2:n-1,n);
Nr=Mr+Ml*m2;
N1=M1+Ml*m1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% discretized model
Co = [0.00042;0.00042;0.00042;0.00042;0.00042];
Tspan = [0 : 1 : 128];
[t,yr] = ode45('PFR_plot2',Tspan,Co);
concn_model = yr(:,end) ;
%% simulation을 통해 구한 농도값과 실험값을 동시에 그리기
clc
clear all
global Nr N1
% 구간 나눔
n = 7 ;
%공정상수
u = 5.17 ;
D = 1 ;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% far-side BC
z = [0:1/(n-2):1 1+1/(n-2)];
k = 10 ;
in = linspace(0,1,n-1);
% 노드 간격 h
for i=1:n-1
h(i)=z(i+1)-z(i);
end
% 행렬 N
for i=1:n-2
for j=1:n-2
if i==j
N(i,j)=2*(h(i)+h(i+1));
elseif i==j+1
N(i,j)=h(i);
elseif i==j-1
N(i,j)=h(j);
end
end
end
% 행렬 R
for i=1:n-2
for j=1:n
if i==j
R(i,j)=6*(1/h(i));
elseif i==j-1
R(i,j)=-6*((h(i)+h(i+1))/(h(i)*h(i+1)));
elseif i==j-2
R(i,j)=6*(1/h(i+1));
end
end
end
% 계수 beta prime1
for i=1:n-1
beta1prime(i)=-(1/3)*h(i);
beta2prime(i)=-(1/6)*h(i);
beta3prime(i)=-(1/h(i));
beta4prime(i)=(1/h(i));
end
% 계수 beta prime2
for i=1:n-1
beta1prime2(i)=1;
beta2prime2(i)=0;
beta3prime2(i)=0;
beta4prime2(i)=0;
end
% 행렬 S prime
for i=1:n
for j=1:n-2
if i==j
Sprime(i,j)=beta2prime(i);
elseif i==j+1
Sprime(i,j)=beta1prime(i);
end
end
end
Sprime(n,n-2)=h(n-1)/6;
% 행렬 T prime
for i=1:n-1
for j=1:n
if i==j
Tprime(i,j)=beta3prime(i);
elseif i==j-1
Tprime(i,j)=beta4prime(i);
end
end
end
Tprime(n,n-1)=beta3prime(n-1);
Tprime(n,n)=beta3prime(n-1);
% 행렬 S prime2
for i=1:n
for j=1:n-2
if i==j+1
Sprime2(i,j)=beta1prime2(i);
end
end
end
for j=1:n-2
Sprime2(n,j)=0;
end
% 행렬 prime, prime2
A=Sprime*inv(N)*R+Tprime;
B=Sprime2*inv(N)*R;
% Boundary Condition
BC=[zeros(1,n-3) (1/6)*( h(n-1) - (3*(1-k)^2) / h(n-1) ) ]*inv(N)*R + [zeros(1,n-2) -1/h(n-1) 1/h(n-1)];
d1=BC(1);
dr=BC(2:n-1);
dn=BC(n);
m1=-d1/dn;
m2=-dr/dn;
% coefficient of ODE
M=-u*A+D*B;
M1=M(2:n-1,1);
Mr=M(2:n-1,2:n-1);
Ml=M(2:n-1,n);
%% Nr 과 N1
Nr=Mr+Ml*m2;
N1=M1+Ml*m1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% discretized model
% ODE model
Co = [0.00042;0.00042;0.00042;0.00042;0.00042];
Tspan = [0 : 1 : 128];
[t,yr] = ode45('PFR_plot2',Tspan,Co);
yr_exp=csvread('C:\Users\김정욱\Desktop\알티디\매트랩연구\PBR-1\PBR-100.csv'); %% experimental data
figure,plot(1:length(yr),yr(:,end),1:length(yr_exp),yr_exp)
xlabel('Time');ylabel('Concentration');legend('Estimated by model','Experimental data');
% axis([0 50 0 20]);
function dyr = PFR_plot2(t,yr);
global Nr N1
if t>=65 & t<=73
dyr = Nr*yr + N1*(0.00042+0.00124875);
else dyr = Nr*yr + N1*0.00042;
end

키워드

CSTR,   RTD,   RTD 메인,   화공생명공학 요소실험 ,   Residence Time Distribution