MathGroup Archive 2011

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

Search the Archive

Re: Odd behaviour of solution of PDE

  • To: mathgroup at smc.vnet.net
  • Subject: [mg116606] Re: Odd behaviour of solution of PDE
  • From: Alan Ford <fabio.sattin at igi.cnr.it>
  • Date: Mon, 21 Feb 2011 19:28:18 -0500 (EST)
  • References: <201102211034.FAA22078@smc.vnet.net> <ijthmr$mn1$1@smc.vnet.net>

>
> this is difficult to say without seeing the actually code. Could you send
> it?


Hi, thanks for your help.

For instance, even this simple piece of code shows the problem:


Clear[T]


DT[x_,t_] = D[T[x,t],x];

DecayTimeScale = 0.015;
DTBordo = -80 (1- Exp[-t/DecayTimeScale]);

eps = 0.01;
x0 = 1;

TTD = T[x,t] /. Flatten[NDSolve[{ D[ x D[T[x,t],x] ,{x,1}] == x
D[T[x,t],t],
                                 DT[eps,t] == 0, T[1,t] == DTBordo,
T[x,0] == 0}, T[x,t] , {x,eps,x0}, {t,0,0.5}];


Plot[ TTD /. x -> 0.1 , {t , 0, 0.5}, PlotRange -> All]



Try using first x0 = 1 , and then x0 = 0.999


Thanks


Fabio




>
> Oliver- Nascondi testo citato
>
> - Mostra testo citato -



  • Prev by Date: Rational[a,b] vs Rational[1,2]
  • Next by Date: Re: Misprint in the Documentation?
  • Previous by thread: Re: Odd behaviour of solution of PDE
  • Next by thread: Re: Odd behaviour of solution of PDE