Figure 7.4:

Dimensionless concentration versus dimensionless radial position for different values of the Thiele modulus.

Code for Figure 7.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
# Converted from cvsr1.m - concentration vs radius (sphere, various phi)
import numpy as np
from misc import save_ascii

phi  = [0.1, 0.5, 1., 2.]
small = 1e-6
npts  = 100
r     = np.linspace(small, 3., npts).reshape(-1, 1)
Csa   = 1.
rs    = 1.
tmp   = r.copy()

for k in range(len(phi)):
    Ca = Csa * 3 * rs / r * np.sinh(phi[k] * (r/rs)) / np.sinh(3*phi[k])
    tmp = np.column_stack([tmp, Ca])

save_ascii('cvsr1.dat', tmp)