Figure 7.5:

Effectiveness factor versus Thiele modulus for a first-order reaction in a sphere.

Code for Figure 7.5

Text of the GNU GPL.

main.py

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
# Converted from etanolog.m - eta vs phi (sphere, first order)
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
tmp  = np.column_stack([phi, eta])

save_ascii('etanolog.dat', tmp)