I am trying to model mass transfer in a circular object. I only have one mass balance and 3 boundary and initial conditions.
Program gives an error of "infinite expression encountered" and I couldn't find where the problem is.
Mass balance is:
ep*Deff*(2)*D[cp[r, t], r] + r*ep*Deff*D[cp[r, t], r, r] == r*ep*D[cp[r, t], t] + r*qp*AA*D[cp[r, t], t]
Boundary and initial conditions:
cp[r, 0] == 0,
Derivative[1, 0][cp][0, t] == 0,
-Deff*Derivative[1, 0][cp][rp, t] == kg*(3*ec/rp)* 0.0024 - cp[rp, t])
Thank you for your help!