Hello!
I am trying to use the pore diffusion model which is given by:
(equation 1)where q is the amount adsorbed, Ds is the diffusion coefficient of adsorbate in the solid phase, r is the particle radius coordinate. This model is usually applied on the an adsorbent particle.
The initial condition is:
and the boundary conditions are:
where, rp is the particle radius, kf is the external mass transfer coefficient, rho_p is the particle density, Ct is the concentration in the bulk fluid at time t, same as "C" in the code in the forum linked above and Cs is the concentration of adsorbate on the surface of the particle.
The dqdt is replaces the dqdt in the following equation:
(equation 2)and has the same boundary and initial conditions as given in the link.
I modified the code but I am confused. Equation 1 and 2 have to be solved simultaneosuly, but this introduces an additional spatial parameter r. So one of the ways this maybe possible is by averaging the values of q at all r at a particular time instant to calculate dqdt and then replacing it in the equations 2 to calculate dCdt. However, this requires equation 1 to be fully solved before starting the calculations for dCdt.
I need help implementing this code. I could really use some.
