No Title |POLVER05_0 |1
# pressure drop in packed bed using Ergun equation
# units for calculations are in m, kg, s, mol, Pa
d(P) / d(z) = -((1-phi)*v/(Dp*phi^3*Ac))*(B + C) # Ergun equation for packed bed
P(0) = 12E5 # feed pressure (Pa)
B = 150*(1-phi)*mu/Dp # laminar flow term
C = 1.75*massflow/Ac # turbulent flow term
T = 273 +150 #temperature (K)
phi = 0.4 # void fraction
mu = 2.3E-5 # viscosity [kg/(m s)]
Dp = 2E-3 # diameter of catalyst particles (m)
D = 5E-2 #inside diameter of tube (m)
massflow = 0.0139 # total mass flow through bed (kg/s)
z(0) = 0
z(f) = 4.1 # Length of packed bed (m)
FTo = 0.257 # total molar flow rate through reactor (mol/s)
R = 8.314 # gas constant [Pa m3/(mol K)]
Ac = 3.14159*D^2/4 # cross section area of inside of tube (m2)
v = FTo*R*T/P # volumetric flow rate (m3/s)
Pbar = P*1E-5 # Pbar is pressure (bar)
vperh = 3600*v # volumetric flow rate (m3/h)
# Reaction A --> B first order
# Isothermal
# Material balance
d(FA) / d(z) = -k*CA*Ac
FA(0) = 0.257
FAo = 0.257
k = 0.05
CA = FA*P/(FTo*R*T)
X = (FAo-FA)/FAo