Example 6-10 Calculating Concentrations as a Function of Position for NH3 Oxidation in a PFR |POLVER05_0
d(Fa)/d(V) =   r1a+r2a+2*r4c/3   #
d(Fb)/d(V) =   1.25*r1a+.75*r2a+r3b   #
d(Fc)/d(V) =   -r1a+2*r3b+r4c   #
d(Fd)/d(V) =   -1.5*r1a-1.5*r2a-r4c   #
d(Fe)/d(V) =  -.5*r2a-5*r4c/6  #
d(Ff)/d(V) =  -2*r3b  #
Ft =  Fa+Fb+Fc+Fd+Fe+Ff  #
r1a =  -5*8*(Fa/Ft)*(Fb/Ft)^2  #
r2a =  -2*4*(Fa/Ft)*(Fb/Ft)  #
r4c =  -5*3.175*(Fc/Ft)*(Fa/Ft)^(2/3)  #
r3b =  -10*8*(Fc/Ft)^2*(Fb/Ft)  #
Ca =  2*Fa/Ft  #
V(0)=0
Fa(0)=10
Fb(0)=10
Fc(0)=0
Fd(0)=0
Fe(0)=0
Ff(0)=0
V(f)=10