
285Homogeneous Reactors
function pe = cal_pe(theta_bar,sigma_sqr)
r = sigma_sqr/(theta_bar
∧
2) ;
pe_o = 2/r ;
pe_n = 2/(r + (2/pe_o
∧
2)*(1-exp(-pe_o))) ;
while abs(pe_n - pe_o) > = 0.0001*pe_o
pe_o = pe_n ;
pe_n = 2/(r + (2/pe_o
∧
2)*(1-exp(-pe_o))) ;
end
pe = pe_n ;
r_cal = (2/pe) - (2/pe
∧
2)*(1-exp(-pe)) ;
PROGRAM: non_id_dead_bypass.m
% Program to detect dead space and bypass in nonideal cstr and estimate
% conversion of I order reaction
clear all ;
% impluse tracer test data
imp_tr_data = [ 0 0.05 0.1 0.2 0.3 0.4 0.5 0.8 1.0 2.0 3.0 4.0 5.0 6.0
7.0 8.0 9.0 10.0 12.0 14.0 16.0 18.0 20.0 22.0 24.0 26.0 ;
0 3.38 3.34 3.3 3.2 3.13 3.05 3.00 2.74 2.61 2.0 1.54 ...