Solve substrate concentration profile in a spherical gel (solution via "odesolve")
Instructor: Nam Sun Wang
Substrate diffusion in a spherical bead in the dimensionless form:
B.C.:
Dimensionless model parameters and rate expression:
Transform the above equation into two 1st-order ODEs:
B.C.:
Plot of substrate profile
Compute the effectiveness factor, which is (observed rate / max rate without mass transfer limitation):
Another approach without decomposing into two 1st-order ODEs.
B.C.:
Double check B.C.
The following does not work because of 1/r at r=0
B.C.:
However, one can cheat without affecting numerical accuracy by starting at a value slightly above r=0.
B.C.:
Double check B.C.