Example CD13-1: Maximum Mixedness Model with Asymmetric RTD (Multiple Reactions) |POLVER05_0
d(ca)/d(z) =   -(-ra+(ca-cao)*EF)   #
d(cb)/d(z) =   -(-rb+(cb-cbo)*EF)   #
d(cc)/d(z) =   -(-rc+(cc-cco)*EF)   #
d(cd)/d(z) =  -(-rd+(cd-cdo)*EF)  #
d(ce)/d(z) =  -(-re+(ce-ceo)*EF)  #
d(F)/d(z) =   -E   #
cao =   1   #
cbo =  1  #
cco =   0   #
cdo =  0  #
ceo =  0  #
k1 =    1   #
k2 =    1   #
k3 =  1  #
ra =  -k1*ca*cb-k2*ca  #
rc =  k1*ca*cb  #
rb =  -k1*ca*cb-k3*cb*cd  #
rd =  k2*ca-k3*cb*cd  #
re =  k3*cb*cd  #
lam =   6-z   #
E1 =   0.47219*lam^4-1.30733*lam^3+0.31723*lam^2+0.85688*lam-0.20909   #
E2 =   3.83999*lam^6-58.16185*lam^5+366.2097*lam^4-1224.66963*lam^3+2289.84857*lam^2+2265.62125*lam-925.4646   #
EF =   E/(1-F)   #
E =   if(lam<=1.26)then(E1)else(if(lam<=2.8)then(E2)else(E3))   #
E3 =   0.00410*lam^4-0.07593*lam^3+0.52276*lam^2-1.59457*lam+1.84445   #
z(0)=0
ca(0)=1
cb(0)=0
cc(0)=0
cd(0)=0
ce(0)=0
F(0)=.99
z(f)=6