Group Abstract Group Abstract

Message Boards Message Boards

Partial differential equation (Diffusion equation in spherical coordinates)

Dodji Cohovi, Student
Posted 10 years ago

Can I please get a help with partial differential equation in mathematica? I tried but got no where. It's a diffusion in a spherical coordinate. Here is the how I tried to solve it and IO added the attachment.

rin = 0.003; rout = 0.00345; cmeq = 0.0245;

sol1 = DSolve[{D[cm[t, r], t] == (1/r)^2* D[difcoef*r^2*D[cm (t, r), r], r], cm (0, r) == cmo, 
cm (t, rout) == cmeq, D[cm (t, rin), r] == 0}, cm (t, r), {t, r}]
Attachments:
POSTED BY: Dodji Cohovi
1 Reply
Sort By:

I don't get anything from DSolve, but NDSolve gives a solution with proper syntax:

With[{difcoefr = 1, cmo = 1, rin = 1, rout = 2, cmeq = 1}, 
 NDSolve[{D[cm[t, r], t] == (1/r)^2*D[difcoefr^2 D[cm[t, r], r], r], 
   cm[0, r] == cmo, 
   cm[t, rout] == cmeq, (D[cm[t, r], r] /. r -> rin) == 0}, 
  cm, {t, 0, 1}, {r, rin, rout}]]
POSTED BY: Gianluca Gorni
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Tag limit exceeded
Note: Only the first five people you tag will receive an email notification; the other tagged names will appear as links to their profiles.
Reply Preview
Attachments
Remove
or Discard