these are the equation given:
d(Temp)/d(z)= B2 x RB + B3 x RC - B4 x (Temp - Tj); temperature in the reactor (K)(energy balance)
d(XB)/d(z)= B1 x RB; conversion to phthalic anhydride (mol balance)
d(XC)/d(z)= B1 x RC; conversion to CO and CO_2 (mol balance)
G= 4684; superficial mass velocity (kg/m^2-h)
MM= 0.02948; mean molecular weight (g/mol)
NA0= 9.27E-3; inlet mole fraction og o-xylene
N0= 0.208; mole fraction of oxygen
CP= 0.25; specific heat (kcal/kg K)
H1= -307; heat of reaction A --> B (kcal/mol)
H3= -1090; heat of reaction A --> C (kcal/mol)
RHOB= 1300; catalyst bulk density (kg/m^3)
DP= 3E-3; catalyst particle diameter (m)
U= 82.7; heat transfer coefficient (kcal/m^2 -h-K)
DT= 0.025; tube diameter (m)
R= 1.897; gas constant (cal/mol-K)
Tj=660; cooling jacket temp (K)
B1= RHOB x MM/(G x NA0); combined term in the mole balance eqs.
B2= RHOB x (-H1)/(G x CP); combined terms in the energy balance eqs.
B3= RHOB x (-H3)/(G x CP); combined terms in the energy balance eqs.
B4= 4 x U/(G x CP x DT); combined terms in the energy balance eqs.
XA= XB + XC; conversion of [I]o[I]-xylene
NA= NA0 x (1 - XA); mole fraction of [I]0[I]-xylene
NB= NA0 x XB; mole fraction of phthalic anhydride
K1= 100 x exp(-27000/(R x Temp) + 19.837); kinetic constant of the 1st reaction (kmol/kg(cat)-hK2= 100 x exp(-31400/(R x Temp) + 20.86); kinetic constant of the 2nd reaction (kmol/kg(cat)-hK3= 100 x exp(-28600/(R x Temp) + 18.97); kinetic constant of the 3rd reaction (kmol/kg(cat)-hRA= -(K1+K3) x NA x N0; rate of reaction of [I]o[I]-xylene (kmol/kg(cat)-h)
RB= K1 x NA x N0 - K2 x NB x N0; rate of generation of phthalic anhydride (kmol/kg(cat)-h)RC= K3 x NA x N0 + K2 x NB x N0; rate of generation of CO and CO_2 (kmol/kg(cat)-h)XB(0)= 0; initial conversion to phthalic anhydride
XC(0)= 0; initial conversion to Co and Co_2
Temp(0)= 650; initial temp. in the reactor (K)
z(0)= =
z(f)= 8; reactor length (m)
I wanted to graph temperature (K) at y-axis over reactor length (m) at x-axis being execute from the equation above using matlab ODE45 in maple. Can anyone help me? I need it to submit my assignment as appendix. I tried numerous method but I failed as I'm just a beginner in maple.
Thanks!
P/s: _ represent subscript & ^ represent superscript