MATLAB: Unsteady difussion in spherical coordinates boundary conditions pdepe

boundary conditionsMATLABpdepesphericalunsteady difussion

Hi, I'm having trouble defining the boundary conditions of my unsteady difussion problem in spherical coordinates using the pdepe function of MATLAB. the boundary conditions are as follows: for r=0 and t>0, dC/dt=0 for r=R and t>0, C=0
Right now Im trying with pl=0, ql=0, pr=0, qr=1 but I just get a flat line. Any help would be greatly appreciated.

Best Answer

To define dC/dt = 0 at r=0 does not make sense. Maybe you mean dC/dr = 0 ?
If this is the case, try
pl=0;
ql=1;
pr=u;
qr=0;
Best wishes
Torsten.
Related Question