Re: Help me : Solve a simple PDE in Mathematica
- To: mathgroup at smc.vnet.net
- Subject: [mg110505] Re: Help me : Solve a simple PDE in Mathematica
- From: schochet123 <schochet123 at gmail.com>
- Date: Mon, 21 Jun 2010 04:55:19 -0400 (EDT)
- References: <hvmvoo$pp1$1@smc.vnet.net>
The reason for the message about a condition "not specified on single edge" is that the boundary condition b2 is not given on an edge of the domain. The fact that condition b3 is mentioned in the error message instead of b2 seems to be a bug in the error handling code. However, if you change the domain to {r,10^{-5),1} you will discover the sad fact that NDSolve does not handle pure multidimensional boundary value problems, only initial boundary value problems. You will therefore have to code it yourself or look for a person or a Package that has already done so. Note that the demonstration you cite uses an exact formula from some book, although they claim, without providing code that NDSolve obtains their solution, which is hard to tell since their boundary conditions seem either garbled or incomplete. Steve On Jun 21, 9:10 am, thaihang le <thaihang... at gmail.com> wrote: > > eqn = D[u[r,z],{z,2}]+D[u[r,z],{r,2}+D[u[r,z],{r,1}]*1/r == 0 > > b1 = ( D[u[r,z],{z,1}]/.z->0 ) ==0 > b2 = ( D[u[r,z],{r,1}]/.r->10^-5 ) ==0 > b3= u[r,2]==1 > b4 =u[2,z] ==1 > > NDSolve[{eqn,b1,b2,b3,b4},u,{r,0,2},{z,0,2}] ==> Error : u[2,z]=== 1 is > not specified on single edge > > and i dont use b4 : > NDSolve[{eqn,b1,b2,b3},u,{r,0,2},{z,0,2}] ===> Error : Number of > constraint (1) is not equal total diff (2). >