Figure 7.6:

Effectiveness factor versus Thiele modulus for a first-order reaction in a sphere (log-log scale).

Code for Figure 7.6

Text of the GNU GPL.

main.py

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
# Converted from etasphere.m - eta vs phi for sphere, with asymptotes
import numpy as np
from misc import save_ascii

nphi   = 50
phi    = np.logspace(-2, 2, nphi).reshape(-1, 1)
eta    = (1./np.tanh(3*phi) - 1./(3*phi)) / phi
line1  = 1./phi
phiasy = np.logspace(-2, np.log10(2.), nphi).reshape(-1, 1)
line2  = np.ones((nphi, 1))
tmp    = np.column_stack([phi, eta, line1, phiasy, line2])

save_ascii('etasphere.dat', tmp)