I've been trying to figure out this problem for a while and was wondering what you thought. I haveThe following pde is (approximately) the PDEleading order homogenized form of the local mass transport equation with a non-linear metabolism (in symmetrical spherical co-ordinates):
\begin{equation*} \frac{\partial c}{\partial t}+\frac{K}{r^2}\frac{\partial c}{\partial r}+\frac{Da~c}{c+n}=0 \end{equation*}
The initial condition is $c(r,0) = c_{in}$ and boundary conditionMy question is $\partial c(R_{in},t) \over \partial r $ $= 0 $, where $R_{in}$if this is the radius of an innerapplied over some 3D sphere (and $R_{out}$ the radius ofwith an outer and an inner sphere). $K$ and $Da$ are both constants. $r$ is the radius measured from then what would be the centre ofappropriate boundary/initial conditions? A no flux BC on the inner sphere outwards, and $t$ time.
I have attempted to solve this using the method of characteristics but the solution keeps reducingseems to only depend on time. And even then I don't end up usingimply that concentration should be zero throughout the second Neumann boundary condition.
Does anyone have any ideas as to where I'm going wrongdomain. Is it just the boundary conditions that(K,Da,n are causing the problem?constants).