%%----------- mass transfer and heat transfer ---------------------

clc
clear

% ------ physical properties --------

mio=1.305e-3;   % viscosity of water at 10C (pa.s)  
ro=999;        % density of water at 10C (kg/m^3) 
cp=4191;       % heat capacity of water at 10C (J/kg.K) 
k=0.5786;      % thernal conductivity of water at 10C (W/m.K)  
Dm=0.75e-9;   % Diffusion Coefficient of Carbamazepine in water at 10C (m^2/s)    

%--------- Geometrical Model ----------

a=600e-6;                               % pore size (m) 
porosity=0.94;                         % void fraction (-)
Ds=a*((4*(1-porosity))/(3*pi)).^0.5;   % strut diameter (m)
Sv=4/Ds*(1-porosity);                 % external surface area per unit volume of bed (m^2/m^3)   

% --------------- mass transfer and heat transfer coefficient --------------------
D=0.019;          % reactor (tube) diameter (m)
%F=[100 200 300 400 500];
F=250;
F=F*(1e-6/60);   % flow rate (m^3/s)
A=pi*D^2/4;       % cross section area (m^2)
u=F/A;            % velocity (m^2/s)


Re=(ro*u*Ds)/mio;    % Reynold number (-)
Sc=mio/(ro*Dm);      % Schmidt number (-)
Pr=(cp*mio)/k;       % Prandtl number (-)

Sh= 1.1*(Re.^0.43)*(Sc.^0.333)     % Sherwood number (-)
Nu= 0.9*(Re.^0.42)*(Pr.^0.333);    % Nusselt number (-)

K=Sh*Dm/Ds; K=K';                       % mass transfer coefficient (m/s)
h=Nu*k/Ds; h=h';                        % heat transfer coefficient (W/m^2.K)

Re=Re'; Sh=Sh';Nu=Nu';u=u';
AA=Sh/Sc.^0.333;