[SciPy-user] first order pde
Donald Fredkin
drfredkin at ucsd.edu
Fri Oct 5 14:46:21 EDT 2007
Issa Karambal wrote:
> Hi all,
>
> I am having problem to solve numerically a hyperbolic equation like
> u_x=-pi*u_t, u(x,0)=exp(sin(x)), u(0,t)=u(2pi,t)=1. I used the method
> of lines; so I discretized the spatial domain using the centered
> finite difference method and thereafter I used the implicit Euler
> method, but I am still having problem after long time integration.
> If anyone has an idea on how I can solve numerically my problem?
>
> THANKS,
>
> issa
You can use the method of characteristics. The example you give is
easily solved numerically (and there is trouble with the boundary
condition at x = 2pi, as you might expect). If you must solve the
equation numerically, I'd suggest an upwind differencing method.
Don
--
More information about the SciPy-User
mailing list