PFR Problem 13, Peters et al. (2003)

r = 1.987;  
kS = 
      Exp[14.76 - 55000/r/temp];   
kS1=Exp[
        14.76 - 55000/r/(800 + 460)]  
eqConst = 0.10*Exp[1200/r/temp];  
eqConst1 = 0.10*Exp[1200/r/(800 + 460)]

aS = 
      1.25;   
cS0 = 
      0.036*45*0.35/32;   
pHi = 
      0.85*2500;  
pHo = 0.95  pHi;   
cS = 
      cS0*(1 - xS);
pH = (pHi + pHo)/2;
pHSo = 
      pHi - pHo;   
pHS = (0 + pHSo)/
        2;
rate = 
        kS1*aS*cS^2*pH/
              cS0/(1 + eqConst1*pHS);
soln = 
      Solve[ratio == cS0*Integrate[1/rate, {xS, 0, 0.74}], 
        ratio]

soln1 = soln[[1, 1, 2]]

slurryRho = 45*0.35 + 65*0.65;
q = 50000/24/4.2*100/35*2000/
          slurryRho  
amtCatal = soln1*q/
          2000  
reactorVol = amtCatal*2000/42
diam = 
      15;   
reactorLength = reactorVol/(Pi*
            diam^2/4)

lrr3@psu.edu (01/21/2008)