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?

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