I am gettin dizzy looking at your PDEs. Did you scan this in? The quality is horrible.
I would use MOL... Method of lines. Turn your PDEs into a system of ODEs. I did something like this with Comsol. Packed bed, three phase flow down the column and looking at diffusion into the spherical/cylindrical packed bed.
Good luck
