| Author | Message | |||
     Anonymous (anonymous) New member Username: anonymous Post Number: 2 Registered: 12-2003 |
I have implemented a spherical pellet in FlexPDE by mapping R on XY plane. The pellet has been modeled for the first quarter. Lets use Landa instead of R here: Landa^2=x^2+y^2 dLanda(U)=(dx(U)+dy(U))*sqrt(x^2+y^2)/(x+y) The code has beed attached Problems: 1- The result is far from the exact solution. 2- Using FlexPDE R coordinate for calculating exact solution shows another result which is not consistent with values calculed by Landa.
| |||
| Robert G. Nelson (rgnelson) Moderator Username: rgnelson Post Number: 82 Registered: 06-2003 |
The results don't match because you have algebraic errors in your derivation and typographical errors in your descriptor. But the easy way to do this is to declare cylindrical geometry, grid a quarter circle, and use the div(grad()) operator. |
Reaction and Diffusion in Spherical Pellet