MathGroup Archive 2010

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

Search the Archive

Re: Help me : Solve a simple PDE in Mathematica

  • To: mathgroup at smc.vnet.net
  • Subject: [mg110507] Re: Help me : Solve a simple PDE in Mathematica
  • From: Oliver Ruebenkoenig <ruebenko at wolfram.com>
  • Date: Tue, 22 Jun 2010 06:59:28 -0400 (EDT)

Hi,

On Mon, 21 Jun 2010, schochet123 wrote:

> 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

http://library.wolfram.com/infocenter/Conferences/7549/

Perhaps this helps,

Oliver

> 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).
>>
>
>


  • Prev by Date: Re: Named patterns forbidden in Except?
  • Next by Date: Re: Ballistic launch trajectory code
  • Previous by thread: Re: Help me : Solve a simple PDE in Mathematica
  • Next by thread: Re: Help me : Solve a simple PDE in Mathematica