Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2009

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

Search the Archive

NDSolve initial value problem

  • To: mathgroup at smc.vnet.net
  • Subject: [mg99154] NDSolve initial value problem
  • From: Murat Havzalı <gezginorman at gmail.com>
  • Date: Wed, 29 Apr 2009 03:46:54 -0400 (EDT)

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: Nested numerical integration
  • Previous by thread: Re: Inverse function warnings
  • Next by thread: Re: NDSolve initial value problem