Figure 7.17:

Effectiveness factor versus Thiele modulus for different values of the Biot number; first-order reaction in a spherical pellet.

Code for Figure 7.17

Text of the GNU GPL.

main.py

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
# Converted from bioteta.m - effectiveness factor vs Thiele modulus, varying Biot
import numpy as np
from misc import save_ascii

Binv = [0., 0.5, 2., 10.]
nphi = 100
phi  = np.logspace(-2, 2, nphi).reshape(-1, 1)
tmp  = phi.copy()

for k in range(len(Binv)):
    eta = (1./phi * (1./np.tanh(3*phi) - 1./(3*phi)) /
           (1 + Binv[k] * (phi/np.tanh(3*phi) - 1./3)))
    tmp = np.column_stack([tmp, eta])

save_ascii('bioteta.dat', tmp)