MathGroup Archive 2001

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

Search the Archive

Boundary value problem with NDSolve

  • To: mathgroup at smc.vnet.net
  • Subject: [mg27414] Boundary value problem with NDSolve
  • From: "ronehwa" <ronehwa.eo85g at nctu.edu.tw>
  • Date: Sun, 25 Feb 2001 00:53:37 -0500 (EST)
  • Organization: NCTU-CC
  • Sender: owner-wri-mathgroup at wolfram.com

Hi Dear All,
Consider the Following question:

eq1=D[f[z,t],t]+v D[f[z,t],z]
eq2=D[g[z,t],t]-v D[g[z,t],z]
twowave=NDSolve[ { eq1 == 0,eq2 == 0,
                             f[z,0] == z+0.32*132/(1-0.32^2), g[z,0]
== -z+132/(1-0.32^2),
                             f[0,t] == 0.32*g[0,t], 0.32*f[100,t] ==
g[100,t] },
                             {f,g},{z,0,100},{t,0,80}] ;

eq1==0 represents f[z,t] is a right-going propagation wave,
 and eq2==0 represents g[z,t] is a left propagation wave.
v is the propagation velocity.
At the boundary z=100, the right-going wave refects 32% its amount and turn
left to propagate,
At the boundary z=0, the left-going wave also refects 32% its amount and
turn right to propagate.

After executing twowave, the mathematica 4.0 return the following message:
"NDSolve::"bcart": "Warning: An insufficient number of boundary conditions
have been specified. Artificial boundary effects may be present in the
solution."
I think the key to the question should be the boundary condition. How can I
set the the boundary condition properly?
Please help me solve this problem, I will deeply appreciate your help.


ronehwa.




  • Prev by Date: Re: Re: page breaks ?
  • Next by Date: Re: How to define a polyfunction?
  • Previous by thread: RE: Re: converting hex strings
  • Next by thread: Re: Boundary value problem with NDSolve