Solving Two PDEs To Derive Traveling Wave Solutions

  • Thread starter Thread starter sigmund
  • Start date Start date
  • Tags Tags
    Derive Pdes Wave
sigmund
Messages
23
Reaction score
0
I have a system of two PDEs:

y_t+(h_0v)_x=0 \quad (1a)
v_t+y_x=0 \quad (1b),

where h_0 is a constant.

Then I want to show that (1) has traveling wave solutions of the form

y(x,t)=f(x-ut) \quad (2a)
v(x,t)=g(x-ut) \quad (2b),

where u is the propagation velocity.

Differentiating (1a) w.r.t. x and (1b) w.r.t. t I conclude that
\frac{\partial^2v}{\partial t^2}=h_0\frac{\partial^2v}{\partial x^2},
but I am not sure that this is useful. I have tried a linear change of variables in order to arrive at a solution. Introducing the variables \alpha=ax+bt and \beta=cx+dt I arrive at a solution for v, namely either v=C(x+t) or v=C(x-t), depending on how the constants a,b,c,d are chosen. Here C is an arbitrary function. Indeed this looks like (2b), but I would like to put forward a more rigorous argument.

I would like some thoughts on this problem, and any hint on how to solve it is welcome. And is it possible to say anything about the relationship between f and g in (2)?
 
Last edited:
Physics news on Phys.org
sigmund said:
I have a system of two PDEs:

y_t+(h_0v)_x=0 \quad (1a)
v_t+y_x=0 \quad (1b),

where h_0 is a constant.

Then I want to show that (1) has traveling wave solutions of the form

y(x,t)=f(x-ut) \quad (2a)
v(x,t)=g(x-ut) \quad (2b),

where u is the propagation velocity.

Differentiating (1a) w.r.t. x and (1b) w.r.t. t I conclude that
\frac{\partial^2v}{\partial t^2}=h_0\frac{\partial^2v}{\partial x^2},
but I am not sure that this is useful. I have tried a linear change of variables in order to arrive at a solution. Introducing the variables \alpha=ax+bt and \beta=cx+dt I arrive at a solution for v, namely either v=C(x+t) or v=C(x-t), depending on how the constants a,b,c,d are chosen. Here C is an arbitrary function. Indeed this looks like (2b), but I would like to put forward a more rigorous argument.

I would like some thoughts on this problem, and any hint on how to solve it is welcome. And is it possible to say anything about the relationship between f and g in (2)?

I'm not sure what more you want. Yes, differentiating as you did eliminates y from the equation leaving
\frac{\partial^2v}{\partial t^2}=h_0\frac{\partial^2v}{\partial x^2}
which is the wave equation. It is well known that
f(x-\sqrt{h_0}t)
and
f(x+\sqrt{h_0}t)
are solutions for any twice differentiable function f, and you can show that by substituting them into your equation.
Of course, if you differentiate the first equation with respect to t and the second with respect to x you can eliminate v, getting the same equation for y.

If all you want to do is show that
(1) has traveling wave solutions of the form

y(x,t)=f(x-ut) \quad (2a)
v(x,t)=g(x-ut) \quad (2b)

you don't really need to do all that. Just substitute those functions: yx= f ' and yt= -u f '; vx= g' and vt= -ug'. Substitute those into the equations and you can see that they satisfy both equations as long as u2= h0.
 
HallsofIvy, thank you for the reply.

I guess it suffices to substitute the assumed solutions into the system, and show that they satisfy the differential equations. This is of course the straightforward way to do it. As Einstein said: "Everything should be done as simple as possible, but no simpler."

What can we say about the relationship between f and g? Obviously, we can read it from the equation, but can we say something in general?
 
Not without initial conditions. For example, if we have some given initial shape, F(x), and are told that the wire or string or whatever is held in that shape and then released (the partial derivative with respect to t is 0 when t=0, a very common situation) then we know that f(x)+ g(x)= F(x) and that cf '(x)- cg'(x)= 0. From the second equation, f '(x)= g'(x) so f(x)= g(x)+ C. From the first then g(x)+ C+ g(x)= F(x) or 2g(x)= F(x)- C so g(x)= (1/2)F(x)- (1/2)C. Then f(x)= g(x)+ C= (1/2)F(x)+ (1/2)C. The solution is (1/2)F(x+ ct)+ (1/2)F(x-ct) (The constants cancel). That is, the initial shape divides into two 1/2 size waves going in opposite directions.

For any functions, f and g, it is possible to come up with initial conditions that will yield them so nothing can be said in general.
 
There are two things I don't understand about this problem. First, when finding the nth root of a number, there should in theory be n solutions. However, the formula produces n+1 roots. Here is how. The first root is simply ##\left(r\right)^{\left(\frac{1}{n}\right)}##. Then you multiply this first root by n additional expressions given by the formula, as you go through k=0,1,...n-1. So you end up with n+1 roots, which cannot be correct. Let me illustrate what I mean. For this...
Back
Top