MathGroup Archive 2009

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

Search the Archive

Re: NDSolve initial value problem

  • To: mathgroup at smc.vnet.net
  • Subject: [mg99203] Re: NDSolve initial value problem
  • From: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
  • Date: Wed, 29 Apr 2009 06:41:32 -0400 (EDT)
  • Organization: Uni Leipzig
  • References: <gt90l4$l1t$1@smc.vnet.net>
  • Reply-to: kuska at informatik.uni-leipzig.de

Hi,

the parts of Piecewise[] outside of the interval y in[-1,1]
are total useless and will not help to make consistent
boundary conditions.

Regards
   Jens

Murat Havzal wrote:
> Dear Mathematica users;
> 
> 
> 
> I am trying to solve liquid heat/mass transfer equation with mixed boundary
> conditions.
> 
> My code looks like this:
> 
> 
> 
> pe=0.5;
> 
> u[y]=1-y^2;
> 
> 
> 
> sol=NDSolve[
> 
> {
> 
> pe*u[y]*D[T[x,y],x]==D[T[x,y],y,y],
> 
> 
> 
> T[0,y]==Piecewise[{{-y,y<=-1},{0,-1< y<1},{1,1<=y}}],
> 
> 
> 
> T[x,1]==1,
> 
> (D[T[x,y],y]/.y->-1)==-1
> 
> 
> 
> },
> 
> T,{x,0,1},{y,-1,1},SolveDelayed->True];
> 
> 
> 
> This returns inconsistent initial boundary conditions error.
> 
> I also tried to make the piecewise initial condition, a numerical function
> namely:
> 
> 
> 
> initial[y_?NumericQ]:=Piecewise[{{-y,y<=-1},{0,-1< y<1},{1,1<=y}}]
> 
> 
> 
> 
> 
> However this did not work, too. I understand that there is a similar
> question already asked about this subject,
> 
> and I tried my best to convert it to my problem but couldn't. I couldn't
> find the code but I also tried to give the
> 
> initial piecewise function as a interpolating function, it returned the same
> response.
> 
> Any help would be appreciated.
> 
> 
> 
> Thanks.
> 
> 
> 
> Murat Havzal
> 
> 


  • Prev by Date: Re: Select and Cases Give Different Answers
  • Next by Date: Re: Select and Cases Give Different Answers
  • Previous by thread: NDSolve initial value problem
  • Next by thread: Re: NDSolve initial value problem