MathGroup Archive 2010

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Help me : Solve a simple PDE in Mathematica

  • To: mathgroup at
  • Subject: [mg110505] Re: Help me : Solve a simple PDE in Mathematica
  • From: schochet123 <schochet123 at>
  • Date: Mon, 21 Jun 2010 04:55:19 -0400 (EDT)
  • References: <hvmvoo$pp1$>

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.


On Jun 21, 9:10 am, thaihang le <thaihang... at> 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).

  • Prev by Date: Re: Workbench fails to build documentation
  • Next by Date: Re: Difficulty with saving Package .m files
  • Previous by thread: Help me : Solve a simple PDE in Mathematica
  • Next by thread: Re: Help me : Solve a simple PDE in Mathematica