5. Solutions to the Wave Equation¶
The wave equation interestingly is a rare example of a partial differential equation (PDE) which exhibits fairly simple solutions, despite its complexity (in general PDEs are quite difficult to solve exactly).
5.1. Analysis of the Wave Equation PDE¶
One way to find solutions is to rewrite the wave equation in a factorised form (by analogy to the difference of two squares)
which we can rewrite in terms of an intermediate variable \(v\):
The form of Equation (5.2) suggests that we can write \(v = v(y)\), where \(y = y(x,\,t)\), meaning that
which can be satisfied by \(y = x + ct\) and so \(v = v(x + ct)\). Solving Equation (5.1), which is now a sourced PDE:
The function \(u\) here will have a homogeneous and an inhomogeneous part:
We find \(u_h\) satisfies:
which we can write as \(u_h(z)\), where \(z = z(x,\,t)\), meaning that
which can be satisfied by \(z = x - ct\) and so \(u_h = u_h(x - ct)\). For \(u_{PI}\) we can see that:
and through a similar process as Equation (5.3), we can show that that \(u_{PI} = u_{PI}(x+ct)\).
The general solution to the wave equation is therefore:
where \(A,\,B\) are functions that we still have not yet found.
We call these travelling wave solutions and we can interpret these two functions as left and right moving wave solutions here, to see how consider the functions first at \(t = 0\):
if we now move to a time \(t \rightarrow 1\), then the function has the form:
For \(u(t,x)\) to stay invariant under this change:
For the \(A\) function, \(x \rightarrow x - c \), suggesting it is a left moving wave.
For the \(B\) function, \(x \rightarrow x + c\) suggesting it is a right moving wave.
5.2. Separable Ansatz for Wave Equation PDE¶
Although this analysis gives us the functional dependence of the functions \(A,\,B\), it does not narrow down their form. We can however use the fact that the PDE is also separable into two ODEs to solve here, our ansatz for the form of \(u\):
Plugging this into the different sides of the wave equation, we find two ODEs result:
this allows us to separate out the variables of the wave equation, which we will write in the form:
These two equations are also equal to a constant. How? Well the LHS here does not vary when we change variable \(t\), so the right hand side cannot either! Equally the RHS does not vary when changing variable \(x\), so the LHS cannot either - meaning there is no dependence on \(x,\,t\) on either side, making them both equal to the same constant.
We are free to pick whichever constant we like here, but lets choose a negative value \(-a^2,\, a\in \mathbb{R}\).
Picking the constant is mathematically arbitrary, but for physical oscillating solutions, we need a negative constant.
We can solve for \(f(x)\) equation to find:
We can do the same with \(g(t)\):
and putting these together means we have:
where we have grouped together the first two terms as the left hand moving wave solutions and last term terms as the right hand moving wave solutions, which matches the functional form of Equation (5.5).
One thing we should try to take care about when writing functions like this is to endure that the arguments of the functions are dimensionless - consider the function \(f(x) = e^x\) where \(x\) is a length, if we look at the Maclaurin expansion:
adding a dimensionless number to a length to an area to a volume really is meaningless! Therefore we should really write functions in the form:
which in each case has a dimensionless quantity being manipulated and \(a\) carries the same units as \(x\).
So lets rewrite \(u(x,\,t)\) in the form:
where \(k\) is a quantity with units \([\textrm{length}]^{-1}\) and \(\omega\) is a quantity with units \([\textrm{time}]^{-1}\). We can make this change, with no less of generality, as we have not fixed a value of constant \(a\) previously and also because it is a more useful way to understand wave behaviour.
We can reduce the form of these wave solutions further, consider the addition of two complex exponentials:
since we know through trigonometric angle identites that:
and therefore we can make the association:
and therefore we can write:
where \(A'\) is some constant to be found. So if we use the condition that we will eventually just take the real parts of our complex solutions, it dramatically simplifies the form of the equations. We call writing our complex solutions of oscillating functions in this way the phasor representation.
Therefore we rewrite Equation (5.6):
where \(A',\, B'\) are constants to be found.
5.3. Features of Wave Solutions¶
Lets examine these quantities \(k,\, \omega\) a little more closely too, start with a right travelling wave \(u\):
where \(u(x=0,\, t=0) = u_0\).
If we shift \(t \rightarrow t + 2\pi/\omega\), then we see that \(u(x,\,t)\) remains invariant, due to the periodic nature of \(e^{i\omega t}\). This suggests that \(2\pi/\omega\) is the time period of the waves oscillation \(T = 2\pi/\omega\) and in fact it is clear that \(\omega = 2\pi/T = 2\pi f\) is just the angular frequency of the wave.
Likewise if we shift \(x \rightarrow x + 2\pi/k\), then \(u(x,\,t)\) remains invariant, due to the periodic nature of \(e^{i k x}\) and this suggests that \(2\pi/k\) is the wavelength of the wave, hence \(k = 2\pi/\lambda\). This is known as the wave vector of the wave.
Although it seems strange that a scalar quantity is being called a vector, this is only true in the one dimensional case, for a right travelling wave propagating in three dimensions:
where \({\bf k},\, \,{\bf u},\, {\bf u_0}\) are now vector quantities, with \({\bf k} = (k_x,\, k_y,\, k_z),\ {\bf x} = (x,\, y,\, z)\), and \(k = |\bf{k}| = 2\pi/\lambda\). \({\bf u_0}\) is known as the waves Polarisation Vector (more on this later).
The quantity \(\Phi = kx-\omega t\) is often known as the phase of the wave and comparing these between two travelling waves \(|\Phi_1 - \Phi_2|\) allows us to look at the phase difference between waves, which is important in phenomena such as interference. All of these considerations also apply to left hand moving waves also.