Figure 6.4:

Conversion of A versus reactor temperature.

Figure 6.4

Code for Figure 6.4

Text of the GNU GPL.

main.py



 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
# Converted from reversible_cstr.m
import numpy as np
from misc import save_ascii

k_10 = 3e7
E_1  = 5838.
K_10 = 1.9e-11
G_1  = -9059.
c_Af = 4.
Q_f  = 250.
V_R  = 2000.

T       = np.linspace(250, 400, 50)
k_1     = k_10 * np.exp(-E_1 / T)
K_1     = K_10 * np.exp(-G_1 / T)

c_A     = (Q_f/V_R + k_1/K_1) * c_Af / (Q_f/V_R + k_1*(1 + 1./K_1))
convvec = (c_Af - c_A) / c_Af
table   = np.column_stack([T, c_A, convvec])

save_ascii('reversible_cstr.dat', table)