Very difficult partial differential equation

Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
33 replies · 31K views
Thanks for your reply. The equation you got for X(x) is the same I tried to solve. Let's do it again:

[tex]\lambda^2+2VX'\lambda + (V^2-c^2)X''=0[/tex]

has solution

[tex]X(x)=\exp\left(r_{\pm} x\right)[/tex]

With [tex]r=\frac{V\lambda \pm c\sqrt{\lambda^2}}{V^2-c^2}[/tex]

Now you can distinguish three cases:
  • [tex]\lambda[/tex] purely real
  • [tex]\lambda[/tex] purely imaginary
  • [tex]\lambda[/tex] = a + bi

[tex]\lambda[/tex] purely real
Then
[tex]r=\frac{-V\lambda \pm c\lambda}{V^2-c^2}=\frac{-V\lambda \pm c\lambda}{(V+c)(V-c)}=\frac{-\lambda}{V \pm c}[/tex]
So
[tex]X(x)=A\exp\left(\frac{-\lambda}{V+c}x\right)+B\exp\left(\frac{-\lambda}{V-c}x\right)[/tex]
With boundary conditions: [tex]X(0)=X(L)=0[/tex] it follows that A + B = 0 and
[tex]X(L)=A\exp\left(\frac{-\lambda}{V+c}L\right)-A\exp\left(\frac{-\lambda}{V-c}L\right)=0\Leftrightarrow A=0[/tex]

It doesn't matter which values lambda takes ([tex]\lambda>0,\lambda<0, \lambda=0[/tex]). [tex]\lambda[/tex] cannot be purely real.

[tex]\lambda[/tex] purely imaginary
Say [tex]\lambda = b i[/tex] Then [tex]r=\frac{-b V i \pm c b i}{V^2-c^2}[/tex].
And
[tex]X(x)=A\exp\left(\frac{-b V i + c b i}{V^2-c^2}x \right)+B\exp\left(\frac{-b V i - c b i}{V^2-c^2}x \right)[/tex]
[tex]X(0)=0[/tex] gives: B=-A, so
[tex]X(x)=A\exp\left(\frac{-bVi}{V^2-c^2}x \right) \left(\exp\left(\frac{c b i}{V^2-c^2}x \right)-\exp\left(\frac{- c b i}{V^2-c^2}x \right)\right)[/tex]
[tex]X(x)=2iA\exp\left(\frac{-bVi}{V^2-c^2}x \right) \left(\sin\left( \frac{bcx}{V^2-c^2}\right))\right)[/tex]
[tex]X(L)=0[/tex] gives [tex]b=\frac{n\pi \left( V^2 -c^2 \right)}{cL}[/tex]
so [tex]\lambda = \frac{n\pi \left( V^2 -c^2 \right)}{cL} i[/tex]

In the end then:

[tex]u(x,t)= \sum_{n=1}^{\infty} a_n \sin\left(\frac{n \pi x}{L} \right)\exp\left(\frac{n\pi V x}{cL}i \right) \exp\left(\frac{n\pi (V^2-c^2) t}{cL}i \right)[/tex]

Did you get this as a solution?

I haven't tried the third case yet.
 
Physics news on Phys.org
Mrx said:
Thanks for your reply. The equation you got for X(x) is the same I
[tex]u(x,t)= \sum_{n=1}^{\infty} a_n \sin\left(\frac{n \pi x}{L} \right)\exp\left(\frac{n\pi V x}{cL}i \right) \exp\left(\frac{n\pi (V^2-c^2) t}{cL}i \right)[/tex]

Did you get this as a solution?

Yes, that's the solution I got! It's tough but we've made it! Note that you can write this in sine and cosine and you can determine the (in this case 2) coefficients via the initial conditions.
 
Okay great! thank you for your help.

Do you know what this equation ectually represents? I can't find anything about the physical meaning, or a derivation of this equation. I tried searching on conveyor belt equation but got nothing.
 
spinoza1989 said:
Could anyone explain a little bit of the physical meaning of this PDE?

If V=0 than you have the wave equation. So c is the speed of the wave (in the belt) and u is the displacement in the y-direction. So what does it means that V=0, I think it means the belt is not moving. What is very strange to me because V is a constant and not time depending.

In the solution given by
Mrx said:
[tex]u(x,t)= \sum_{n=1}^{\infty} a_n \sin\left(\frac{n \pi x}{L} \right)\exp\left(\frac{n\pi V x}{cL}i \right) \exp\left(\frac{n\pi (V^2-c^2) t}{cL}i \right)[/tex]

you can bring the exp-parts together and you get a x-v't where v' is a the speed of the wave relative to the frame. (c is speed relative to the belt and V the speed of the belt relative to the frame).

[tex]2 V u_x_t + V^2 u_x_x[/tex] is the "source therm" I think, first of all since you have just a wave equations and the other therms are source therms. Second reason is that the units of [tex]u_x_t[/tex] are Hertz and mabye it has to do with the rotations of one of the wheels.

The 2V has something to do with the speed between upper belt and lower (returning) belt.

This is all what I could thought about.
Maybe somebody can make a 3dplot in maple?

Greetz Bendavid2