Difference between revisions of "Isothermal reactor design  2013"
Kevin Dunn (talk  contribs) 
Kevin Dunn (talk  contribs) m 

Line 61:  Line 61:  
* [http://learnche.mcmaster.ca/media/3K42013Class06C.mp3 Audio] and [http://learnche.mcmaster.ca/media/3K42013Class06C.mp4 video] recording of the class  * [http://learnche.mcmaster.ca/media/3K42013Class06C.mp3 Audio] and [http://learnche.mcmaster.ca/media/3K42013Class06C.mp4 video] recording of the class  
=== 25 February 2013 (07A) ===  
* [http://learnche.mcmaster.ca/media/3K42013Class07A.mp3 Audio] and [http://learnche.mcmaster.ca/media/3K42013Class07A.mp4 video] recording of the class  
The example covered in class is based on example 48 in F2006 and example 62 in F2011.  
<rst>  
The 3 ODE's are:  
.. math::  
\dfrac{dF_A}{dV} &= r_A\\ \dfrac{dF_B}{dV} &= r_B  R_B \\ \dfrac{dF_C}{dV} &= r_C  
where :math:`r_A = r_B = r_C` and :math:`r_A = k\left(C_A  \dfrac{C_B C_C}{K_C} \right)`, and :math:`R_B = k_\text{diff}C_B`.  
* :math:`k = 0.01\,\text{s}^{1}`  
* :math:`k_\text{diff} = 0.005\,\text{s}^{1}`  
* :math:`K_C = 50\,\text{mol.m}^{3}`  
We derived earlier in the course that  
.. math::  
C_A = C_\text{TO}\left(\dfrac{F_A}{F_T}\right)\left(\dfrac{P}{P_0}\right)\left(\dfrac{T_0}{T}\right)  
Assuming isothermal and isobaric conditions in the membrane:  
.. math::  
C_A = C_\text{T0}\left(\dfrac{F_A}{F_T}\right)  
where :math:`F_T = F_A + F_B + F_C` and :math:`C_\text{T0} = \dfrac{P_0}{RT_0}`  
Using all of the above derivations, we can set up our numerical integration.  
</rst>  
{ class="wikitable"  
  
! MATLAB  
! Python  
  
 width="50%" valign="top"   
In a file called '''<tt>membrane.m</tt>''':  
<syntaxhighlight lang="matlab">  
sdfsdf  
</syntaxhighlight>  
[[Image:PlotsMATLAB.png  550px]]  
 width="50%" valign="top"   
<syntaxhighlight lang="python">  
asdasd  
</syntaxhighlight>  
[[Image:PlotsPython.png550px]]  
} 
Revision as of 00:11, 26 February 2013
Class date(s):  04 February to 14 February  
 
 
 
 
 
 F2011: Chapter 5 and 6
 F2006: Chapter 4
04 February 2013 (05A)
 General problem solving strategy for reactor engineering
 Audio and video recording of the class
06 February 2013 (05B)
 The Ergun equation derivation
 Audio and video recording of the class
07 February 2013 (05C)
 Notes used during the class
 The spreadsheet with the Ergun equation example. Use it to try
 different lengths of reactor
 different catalyst particle sizes
 different pipe diameters
 gas properties (e.g. density)
 to see the effect on pressure drop in the packed bed.
11 February 2013 (06A)
 Audio and video recording of the class
 Codes to solve the example in class are available on the page software for integrating ODEs.
14 February 2013 (06C): midterm review
25 February 2013 (07A)
The example covered in class is based on example 48 in F2006 and example 62 in F2011. <rst> The 3 ODE's are:
.. math::
\dfrac{dF_A}{dV} &= r_A\\ \dfrac{dF_B}{dV} &= r_B  R_B \\ \dfrac{dF_C}{dV} &= r_C
where :math:`r_A = r_B = r_C` and :math:`r_A = k\left(C_A  \dfrac{C_B C_C}{K_C} \right)`, and :math:`R_B = k_\text{diff}C_B`.
 :math:`k = 0.01\,\text{s}^{1}`
 :math:`k_\text{diff} = 0.005\,\text{s}^{1}`
 :math:`K_C = 50\,\text{mol.m}^{3}`
We derived earlier in the course that
.. math:: C_A = C_\text{TO}\left(\dfrac{F_A}{F_T}\right)\left(\dfrac{P}{P_0}\right)\left(\dfrac{T_0}{T}\right)
Assuming isothermal and isobaric conditions in the membrane:
.. math:: C_A = C_\text{T0}\left(\dfrac{F_A}{F_T}\right)
where :math:`F_T = F_A + F_B + F_C` and :math:`C_\text{T0} = \dfrac{P_0}{RT_0}`
Using all of the above derivations, we can set up our numerical integration. </rst>
MATLAB  Python 

In a file called membrane.m: sdfsdf

asdasd
