PDA

View Full Version : solving a partial DE


jimbo007
May26-04, 05:22 AM
hi all
i have been trying to solve to following problem,


\frac{\partial^2 u}{\partial x^2} - \frac{\partial^2 u}{\partial y^2}
+ 2\frac{\partial u}{\partial x} + u = 0



u=u(x,y)


after a bit of work using the change of variables


\zeta=\zeta(x,y)=y-x


and


\eta=\eta(x,y)=y+x


i obtain

- 4\frac{\partial u}{\partial \zeta \partial \eta} = 2\frac{\partial u}{\partial\zeta} - 2\frac{\partial u}{\partial\eta}-u

but i am unsure how to solve this, i used maple to solve this problem and it gave out a fairly harmless answer so i am pretty sure there would be any easy way to solve the above equation.

could someone kindly show me how to obtain a solution to this problem

arildno
May26-04, 05:33 AM
I would suggest a separation of variables solution,
u(x,y)=f(x)g(y)

Max0526
May28-04, 04:58 PM
Hi;
The first substitution to use is "running wave" substitution:

z=x-ct

and you will have:

\frac{\partial^2 u}{\partial x^2} - \frac{\partial^2 u}{\partial y^2}
+ 2\frac{\partial u}{\partial x} + u = 0


u_{xx}-u_{yy}+2u_x+u=0


z=x-ct


u_{zz}-c^2u_{zz}+2u_z+u=0


(1-c^2)u_{zz}+2u_z+u=0

The last equation is a linear ODE of the 2nd order while c<>+1 and c<>-1. Otherwise we have the first order ODE

2u_z+u=0

Anyways both equations are simple and can be solved easily.
Still, this is not the only way of reduction of the given PDE (there could be another substitutions v=v(u,x,y), z=z(u,x,y) that will give us ordinary DEs). If you need some other forms for v and z, let me know, I can find some for you.
Best of luck,
Max.

arildno
May28-04, 05:51 PM
Here's a rather general solution using separation of variables:

1.
u(x,y)=f(x)g(y)\to\frac{f''+2f'+f}{f}=\frac{g''}{g }

2. Since LHS is a function of x, wheras RHS is a function of y, equality is only possible if they are equal to some constant k.
I set a subscript k on the functions at this moment:
\frac{f_{k}''+2f_{k}'+f_{k}}{f_{k}}=k\to{f}_{k}(x) =A_{k}e^{-x-\sqrt{k}x}+B_{k}e^{-x+\sqrt{k}x}
\frac{g_{k}''}{g_{k}}=k\to{g}_{k}(x)=C_{k}e^{-\sqrt{k}y}+D_{k}e^{\sqrt{k}y}

3. Since the original differential equation is linear, a sum of two solutions (for example, for different k-values) is also a solution, and we may write:

u(x,y)=
e^{-x}\int_{-\infty}^{\infty}\mathcal{C}_{(k,+)}Cosh(\sqrt{k}(x +y))+\mathcal{C}_{(k,-)}Cosh(\sqrt{k}(x-y))dk+
e^{-x}\int_{-\infty}^{\infty}\mathcal{S}_{(k,+)}Sinh(\sqrt{k}(x +y))+\mathcal{S}_{(k,-)}Sinh(\sqrt{k}(x-y))dk

4. Negative values of k will yield trigonometric functions; the unknown functions of k,
\mathcal{C}_{(k,+)},\mathcal{C}_{(k,-)},\mathcal{S}_{(k,+)},\mathcal{S}_{(k,-)}
may be suited to satisfy boundary conditions.

jimbo007
May28-04, 11:28 PM
thanks for your help guys
i have the problem solved