374 Green Chemical Engineering
tau0 = 4 ;% time taken for complete conversion in HRS for
r0 = 4 ;% particle of size diameter mm
% CALCULATION
%_________________________________________________________________
vec_size = size(size_dist) ;
n_data = vec_size(1,1) ;
xb_bar = 0 ;
ws = 0 ;
for i = 1:n_data
r = size_dist(i,1) ; % particle size
w = size_dist(i,2) ; % weight fraction
if ((mechanism == 1) || (mechanism == 3))
tau = tau0*(r/r0);
end ;
if mechanism == 2
tau = tau0*(r/r0)^2;
end ;
theta_by_tau = theta_bar/tau ;
xb = cal_xb(theta_by_tau,mechanism) ;
xb_bar = xb_bar + xb*w ;
ws = ws + w ;
end ;
xb_bar = xb_bar/ws ;
% DISPLAY RESULTS
%_________________________________________________________________
fprintf('------------------------------------------------------------\n') ...