Mechanical Waves: A Mathematical Analysis

1. Theoretical Analysis and Derivations

1.1 Wave Kinematics and the Wave Function

A mechanical wave is a disturbance that propagates through a medium, transferring energy and momentum without transferring matter. The displacement of a point in the medium from its equilibrium position is a function of both position $x$ and time $t$, denoted by $y(x, t)$.

1.1.1 General Form of the Wave Function

For a wave propagating with a constant speed $v$ in the positive $x$-direction, the function must maintain its shape over time. This leads to the general form:

$$ y(x, t) = f(x - v t) \tag{1} $$

where $f$ is any function that describes the shape of the wave profile. A wave propagating in the negative $x$-direction is given by $y(x, t) = f(x + v t)$.

1.1.2 Harmonic Wave Function (Sinusoidal)

The most common type of wave for analysis is the harmonic (sinusoidal) wave. For a wave propagating in the positive $x$-direction, the displacement is:

$$ y(x, t) = A \cos(k x - \omega t + \phi) \tag{2} $$

Here, $A$ is the amplitude, $k$ is the wave number, $\omega$ is the angular frequency, and $\phi$ is the phase constant. The fundamental relationships connecting these parameters are:

$$ k = \frac{2\pi}{\lambda} \quad ; \quad \omega = 2\pi f = \frac{2\pi}{T} \quad ; \quad v = \frac{\omega}{k} = \lambda f \tag{3} $$

where $\lambda$ is the wavelength, $f$ is the frequency, and $T$ is the period.

1.2 The Classical Wave Equation (CWE)

The classical wave equation is a linear, second-order partial differential equation that governs the propagation of many types of waves, including mechanical waves.

1.2.1 Proof of the CWE from the General Solution

Consider the general wave function $y(x, t) = f(x - v t)$. Let $u = x - v t$. The wave function is $y(u)$. We apply the chain rule to find the partial derivatives with respect to $x$ and $t$.

First-order partial derivatives:

$$ \frac{\partial y}{\partial x} = \frac{dy}{du} \frac{\partial u}{\partial x} = f'(u) (1) = f'(u) \tag{4} $$
$$ \frac{\partial y}{\partial t} = \frac{dy}{du} \frac{\partial u}{\partial t} = f'(u) (-v) = -v f'(u) \tag{5} $$

Second-order partial derivatives:

$$ \frac{\partial^2 y}{\partial x^2} = \frac{\partial}{\partial x} \left( \frac{\partial y}{\partial x} \right) = \frac{\partial}{\partial x} (f'(u)) = f''(u) \frac{\partial u}{\partial x} = f''(u) \tag{6} $$
$$ \frac{\partial^2 y}{\partial t^2} = \frac{\partial}{\partial t} \left( \frac{\partial y}{\partial t} \right) = \frac{\partial}{\partial t} (-v f'(u)) = -v \left( f''(u) \frac{\partial u}{\partial t} \right) = -v (f''(u) (-v)) = v^2 f''(u) \tag{7} $$

By comparing equations (6) and (7), we can eliminate $f''(u)$ to obtain the one-dimensional Classical Wave Equation:

$$ \frac{\partial^2 y}{\partial x^2} = \frac{1}{v^2} \frac{\partial^2 y}{\partial t^2} \tag{8} $$

1.3 Speed of Transverse Waves on a String

We now derive the wave speed $v$ by applying Newton's Second Law to an infinitesimal segment of a string under tension $T$.

1.3.1 Physical Derivation of Wave Speed $v$

Consider a small segment of length $dx$ of a string with linear mass density $\mu$. The tension $T$ acts tangentially to the string at both ends of the segment. We assume the vertical displacement $y$ is small, so the tension $T$ remains approximately constant.

The net vertical force $dF_y$ on the segment is the difference between the vertical components of the tension at $x+dx$ and $x$.

$$ dF_y = T \sin(\theta_{x+dx}) - T \sin(\theta_x) \tag{9} $$

For small angles (small transverse displacement), $\sin\theta \approx \tan\theta$. The slope of the string at any point is $\tan\theta = \frac{\partial y}{\partial x}$.

$$ dF_y \approx T \left( \frac{\partial y}{\partial x} \Big|_{x+dx} - \frac{\partial y}{\partial x} \Big|_x \right) \tag{10} $$

According to Newton's Second Law, $dF_y = (dm) a_y$. Since $dm = \mu \, dx$ and $a_y = \frac{\partial^2 y}{\partial t^2}$:

$$ T \left( \frac{\partial y}{\partial x} \Big|_{x+dx} - \frac{\partial y}{\partial x} \Big|_x \right) = \mu \, dx \, \frac{\partial^2 y}{\partial t^2} \tag{11} $$

Dividing by $dx$ and taking the limit as $dx \to 0$, the term on the left side becomes the definition of the second derivative of the slope with respect to $x$:

$$ \lim_{dx \to 0} \frac{1}{dx} \left[ \frac{\partial y}{\partial x} \Big|_{x+dx} - \frac{\partial y}{\partial x} \Big|_x \right] = \frac{\partial^2 y}{\partial x^2} \tag{12} $$

Substituting this back into equation (11), we get the CWE specifically for a string:

$$ T \frac{\partial^2 y}{\partial x^2} = \mu \frac{\partial^2 y}{\partial t^2} \quad \implies \quad \frac{\partial^2 y}{\partial x^2} = \frac{\mu}{T} \frac{\partial^2 y}{\partial t^2} \tag{13} $$

Comparing this result (13) with the general CWE (8), we find that $\frac{1}{v^2} = \frac{\mu}{T}$. Thus, the wave speed is:

$$ v = \sqrt{\frac{T}{\mu}} \tag{14} $$

The speed of a mechanical wave is determined solely by the properties of the medium: the restoring force (Tension $T$) and the inertia (linear mass density $\mu$).

1.4 Energy and Power Transmission

Waves carry energy. We calculate the rate at which energy is transferred across a point, known as the wave power $P$.

1.4.1 Kinetic Energy Density

The infinitesimal kinetic energy $dK$ of a segment $dx$ is $dK = \frac{1}{2} (dm) v_y^2$, where $v_y = \frac{\partial y}{\partial t}$ is the transverse velocity.

$$ dK = \frac{1}{2} (\mu \, dx) \left(\frac{\partial y}{\partial t}\right)^2 \quad \implies \quad \frac{dK}{dx} = \frac{1}{2} \mu \left(\frac{\partial y}{\partial t}\right)^2 \tag{15} $$

For a harmonic wave $y(x, t) = A \cos(k x - \omega t)$, the velocity is $\frac{\partial y}{\partial t} = \omega A \sin(k x - \omega t)$.

$$ \frac{dK}{dx} = \frac{1}{2} \mu \omega^2 A^2 \sin^2(k x - \omega t) \tag{16} $$

1.4.2 Potential Energy Density

The potential energy $dU$ stored in the segment is due to the stretching of the string. The length change $\Delta L = ds - dx$ is approximated using the binomial expansion of $ds = dx \sqrt{1 + \left(\frac{\partial y}{\partial x}\right)^2}$ as $\Delta L \approx \frac{1}{2} \left(\frac{\partial y}{\partial x}\right)^2 dx$. The potential energy $dU = T \Delta L$:

$$ \frac{dU}{dx} = \frac{1}{2} T \left(\frac{\partial y}{\partial x}\right)^2 \tag{17} $$

Using $T = \mu v^2$ and $v=\omega/k$ (which gives $T = \mu \omega^2/k^2$), we find that the potential energy density is equal to the kinetic energy density:

$$ \frac{dU}{dx} = \frac{1}{2} \mu \omega^2 A^2 \sin^2(k x - \omega t) \tag{18} $$

Crucially, $\frac{dK}{dx} = \frac{dU}{dx}$ for a travelling harmonic wave.

1.4.3 Average Power Transmitted

The instantaneous power $P$ is $P(x, t) = -T \frac{\partial y}{\partial x} \cdot \frac{\partial y}{\partial t}$. After substitution and averaging $\langle \sin^2(\theta) \rangle = \frac{1}{2}$, the average power $\bar{P}$ transmitted is found to be:

$$ \bar{P} = \frac{1}{2} \mu v \omega^2 A^2 \tag{19} $$

The average power is proportional to the medium's inertia ($\mu$), wave speed ($v$), the square of the angular frequency ($\omega^2$), and the square of the amplitude ($A^2$).

1.5 Principle of Superposition and Standing Waves

The principle of superposition states that the resultant displacement $y_R(x, t)$ is the vector sum of individual wave displacements: $y_R(x, t) = y_1(x, t) + y_2(x, t)$.

1.5.1 Standing Waves Derivation

A standing wave results from the superposition of two identical waves traveling in opposite directions: $y_1(x, t) = A \sin(k x - \omega t)$ and $y_2(x, t) = A \sin(k x + \omega t)$.

$$ y_S(x, t) = A [\sin(k x - \omega t) + \sin(k x + \omega t)] \tag{20} $$

Using the trigonometric identity $\sin \alpha + \sin \beta = 2 \sin\left(\frac{\alpha+\beta}{2}\right) \cos\left(\frac{\alpha-\beta}{2}\right)$:

$$ y_S(x, t) = \left[ 2 A \sin(k x) \right] \cos(\omega t) \tag{21} $$

This form shows the wave is stationary: the spatial dependence is separated from the time dependence. The amplitude of oscillation is dependent only on position $x$: $A_S(x) = 2 A \sin(k x)$.

2. Problems (10)

These problems range from intermediate calculations to advanced conceptual derivations. Solutions follow in the next section.

  1. P1. (Intermediate) A wave on a string is described by $y(x, t) = 0.05 \cos(4\pi x - 20\pi t + \pi/4)$, where $x$ and $y$ are in meters and $t$ is in seconds. Determine the phase velocity $v$, wavelength $\lambda$, and the maximum transverse velocity $v_{y, \text{max}}$ of any element of the string.
  2. P2. (Intermediate) A $10 \text{ m}$ long string with a mass of $100 \text{ g}$ is under tension of $400 \text{ N}$. If a wave pulse is generated at one end, calculate the time $t$ required for the pulse to reach the other end.
  3. P3. (Intermediate) Two harmonic waves are described by $y_1 = 3.0 \sin(2x - 5t)$ and $y_2 = 3.0 \sin(2x - 5t + \pi/3)$. Find the amplitude $A_R$ of the resultant wave formed by their superposition.
  4. P4. (Intermediate) The average power $\bar{P}$ transmitted by a harmonic wave on a string is $10 \text{ W}$. If the tension $T$ is doubled and the amplitude $A$ is halved, what is the new average power $\bar{P}'$, assuming the linear density $\mu$ and angular frequency $\omega$ remain constant?
  5. P5. (Intermediate) A standing wave on a string is described by $y(x, t) = 0.1 \sin(1.5\pi x) \cos(30\pi t)$, where $x$ and $y$ are in meters. Find the distance $\Delta x$ between a node and the nearest antinode.
  6. P6. (Advanced) Given the 3D wave equation: $$\nabla^2 \Psi = \frac{1}{v^2} \frac{\partial^2 \Psi}{\partial t^2}$$ Assume a spherical wave solution $\Psi(r, t) = \frac{f(r-vt)}{r}$. Show that the product $\Phi(r, t) = r\Psi(r, t)$ satisfies the one-dimensional wave equation: $$\frac{\partial^2 \Phi}{\partial r^2} = \frac{1}{v^2} \frac{\partial^2 \Phi}{\partial t^2}$$
  7. P7. (Advanced) A damped harmonic wave is described by $y(x, t) = A_0 e^{-\alpha x} \cos(k x - \omega t)$, where $\alpha$ is the damping coefficient. Calculate the instantaneous power $P(x, t)$ transferred across the point $x$. Express the average power $\bar{P}(x)$ in terms of the initial average power $\bar{P}_0 = \frac{1}{2} \mu v \omega^2 A_0^2$. (Assume $\alpha$ is small such that $\alpha/k \approx 0$).
  8. P8. (Advanced) A semi-infinite string ($x \ge 0$) has two linear mass densities: $\mu_1$ for $x < 0$ and $\mu_2$ for $x \ge 0$. A wave of amplitude $A_I$ is incident from the left ($x < 0$). Determine the amplitude reflection coefficient $R = A_R / A_I$ and the transmission coefficient $T = A_T / A_I$. Use the boundary conditions that the displacement $y(0, t)$ and the slope $\frac{\partial y}{\partial x}\big|_{x=0}$ must be continuous.
  9. P9. (Irodov-like) A uniform flexible string of total length $L$ and mass $M$ is suspended vertically from a rigid support. A transverse wave pulse is generated at the bottom end of the string. Determine the total time $t$ required for the pulse to travel from the bottom end to the support. (Hint: The tension $T(y)$ varies with height $y$ from the bottom).
  10. P10. (Irodov-like) An elastic bar with density $\rho$ and Young's modulus $E$ is rigidly fixed at one end ($x=0$) and free at the other ($x=L$). The longitudinal wave speed in the bar is $v = \sqrt{E/\rho}$. A standing wave is set up in the bar. Find the two boundary conditions for the displacement $\xi(x, t)$ and determine the possible resonance frequencies $f_n$ for this bar.

3. Solutions

3.1 Solutions to Intermediate Problems (P1 - P5)

P1. Wave function is $y(x, t) = A \cos(k x - \omega t + \phi)$. Comparing with the given equation, $k = 4\pi \text{ rad/m}$ and $\omega = 20\pi \text{ rad/s}$. $$\text{Phase Velocity: } v = \frac{\omega}{k} = \frac{20\pi}{4\pi} = 5.0 \text{ m/s}$$ $$\text{Wavelength: } \lambda = \frac{2\pi}{k} = \frac{2\pi}{4\pi} = 0.5 \text{ m}$$ $$\text{Max Transverse Velocity: } v_{y, \text{max}} = A\omega = (0.05) (20\pi) = \pi \approx 3.14 \text{ m/s}$$
P2. The time required is $t = L/v$, where $v = \sqrt{T/\mu}$. $$\text{Linear density: } \mu = M/L = 0.100 \text{ kg} / 10 \text{ m} = 0.01 \text{ kg/m}$$ $$\text{Wave speed: } v = \sqrt{400 \text{ N} / 0.01 \text{ kg/m}} = 200 \text{ m/s}$$ $$\text{Time: } t = 10 \text{ m} / 200 \text{ m/s} = 0.05 \text{ s}$$
P3. The resultant amplitude $A_R$ for two waves with phase difference $\phi$ is $A_R = 2 A \left| \cos\left(\frac{\phi}{2}\right) \right|$. Here $A=3.0$ and $\phi = \pi/3$. $$A_R = 2 (3.0) \cos\left(\frac{\pi/3}{2}\right) = 6.0 \cos\left(\frac{\pi}{6}\right)$$ $$A_R = 6.0 \left( \frac{\sqrt{3}}{2} \right) = 3\sqrt{3} \approx 5.20 \text{ m}$$
P4. Average power $\bar{P}$ is proportional to $\sqrt{T}$ and $A^2$. Specifically, $\bar{P} \propto \sqrt{\mu T} \omega^2 A^2$. $$\frac{\bar{P}'}{\bar{P}} = \frac{\sqrt{T'} (A')^2}{\sqrt{T} A^2} = \frac{\sqrt{2T} (A/2)^2}{\sqrt{T} A^2}$$ $$\frac{\bar{P}'}{\bar{P}} = \frac{\sqrt{2}}{1} \cdot \frac{1}{4} = \frac{\sqrt{2}}{4}$$ $$\bar{P}' = 10 \text{ W} \cdot \frac{\sqrt{2}}{4} = \frac{5\sqrt{2}}{2} \approx 3.54 \text{ W}$$
P5. The distance between a node and the nearest antinode is $\lambda/4$. The wave number $k$ is $1.5\pi$. $$\text{Wavelength: } \lambda = \frac{2\pi}{k} = \frac{2\pi}{1.5\pi} = \frac{4}{3} \text{ m}$$ $$\text{Distance between node and antinode: } \Delta x = \frac{\lambda}{4} = \frac{1}{4} \left( \frac{4}{3} \text{ m} \right) = \frac{1}{3} \text{ m}$$

3.2 Solutions to Advanced Problems (P6 - P8)

P6. We are given $\Psi(r, t) = \frac{f(u)}{r}$, where $u=r-vt$. We define $\Phi = r\Psi = f(u)$. $$\text{Partial derivatives of } \Phi \text{ with respect to } r \text{ and } t:$$ $$\frac{\partial \Phi}{\partial r} = f'(u) \frac{\partial u}{\partial r} = f'(u)$$ $$\frac{\partial^2 \Phi}{\partial r^2} = f''(u) \frac{\partial u}{\partial r} = f''(u)$$ $$\frac{\partial \Phi}{\partial t} = f'(u) \frac{\partial u}{\partial t} = -v f'(u)$$ $$\frac{\partial^2 \Phi}{\partial t^2} = -v f''(u) \frac{\partial u}{\partial t} = (-v) f''(u) (-v) = v^2 f''(u)$$ $$\text{Substituting into the 1D wave equation: } \frac{\partial^2 \Phi}{\partial r^2} = f''(u) \quad ; \quad \frac{1}{v^2} \frac{\partial^2 \Phi}{\partial t^2} = \frac{1}{v^2} v^2 f''(u) = f''(u)$$ $$\text{Thus: } \frac{\partial^2 \Phi}{\partial r^2} = \frac{1}{v^2} \frac{\partial^2 \Phi}{\partial t^2}$$
P7. Damped wave: $y(x, t) = A_0 e^{-\alpha x} \cos(k x - \omega t)$. $$\frac{\partial y}{\partial x} = -\alpha A_0 e^{-\alpha x} \cos(\dots) + k A_0 e^{-\alpha x} \sin(\dots)$$ $$\frac{\partial y}{\partial t} = \omega A_0 e^{-\alpha x} \sin(\dots)$$ $$\text{Using the approximation } \alpha/k \approx 0 \text{ (i.e., } \frac{\partial y}{\partial x} \approx k A_0 e^{-\alpha x} \sin(\dots) \text{):}$$ $$\text{Instantaneous Power: } P(x, t) = -T \frac{\partial y}{\partial x} \frac{\partial y}{\partial t} \approx T (k A_0 e^{-\alpha x} \sin) (\omega A_0 e^{-\alpha x} \sin)$$ $$P(x, t) = T k \omega A_0^2 e^{-2\alpha x} \sin^2(k x - \omega t)$$ $$\text{Average Power: } \bar{P}(x) = \langle P(x, t) \rangle = \frac{1}{2} T k \omega A_0^2 e^{-2\alpha x}$$ $$\text{Since } \bar{P}_0 = \frac{1}{2} T k \omega A_0^2 \text{ (from } x=0 \text{):}$$ $$\bar{P}(x) = \bar{P}_0 e^{-2\alpha x}$$
P8. Boundary conditions at $x=0$ require continuity of displacement ($y$) and slope ($\frac{\partial y}{\partial x}$). $$y_I + y_R = y_T \implies A_I + A_R = A_T \tag{i}$$ $$\frac{\partial y_I}{\partial x} + \frac{\partial y_R}{\partial x} = \frac{\partial y_T}{\partial x} \implies k_1 (A_I - A_R) = k_2 A_T \tag{ii}$$ $$\text{From (i) and (ii), the coefficients are:}$$ $$\text{Reflection Coefficient: } R = \frac{A_R}{A_I} = \frac{k_1 - k_2}{k_1 + k_2}$$ $$\text{Transmission Coefficient: } T = \frac{A_T}{A_I} = \frac{2 k_1}{k_1 + k_2}$$ $$\text{Where } k_1 = \omega/v_1 \text{ and } k_2 = \omega/v_2 \text{ (since frequency } \omega \text{ is conserved).}$$

3.3 Solutions to Irodov-like Problems (P9 - P10)

P9. Let $\mu = M/L$ be the linear density. The tension $T(y)$ at height $y$ (measured from the bottom) is the weight of the string segment below it. $$T(y) = (\mu y) g = \frac{M g}{L} y$$ $$\text{Wave speed: } v(y) = \sqrt{\frac{T(y)}{\mu}} = \sqrt{\frac{(M g / L) y}{M/L}} = \sqrt{g y}$$ $$\text{Time element: } dt = \frac{dy}{v(y)} = \frac{dy}{\sqrt{g y}}$$ $$\text{Total time: } t = \int_{0}^{L} \frac{dy}{\sqrt{g y}} = \frac{1}{\sqrt{g}} \int_{0}^{L} y^{-1/2} dy$$ $$t = \frac{1}{\sqrt{g}} \left[ 2 y^{1/2} \right]_{0}^{L} = 2 \sqrt{\frac{L}{g}}$$
P10. Longitudinal wave displacement $\xi(x, t)$. $$\text{1. Fixed end at } x=0 \text{ (Displacement is zero): } \xi(0, t) = 0$$ $$\text{2. Free end at } x=L \text{ (Strain } \partial \xi / \partial x \text{ is zero): } \frac{\partial \xi}{\partial x}\Big|_{x=L} = 0$$ $$\text{Applying } \xi(0, t) = 0 \text{ to the general solution } \xi(x, t) = (C \sin(k x) + D \cos(k x)) \cos(\omega t) \text{ yields } D = 0.$$ $$\text{Applying } \frac{\partial \xi}{\partial x}\Big|_{x=L} = 0 \text{ to } \xi(x, t) = C \sin(k x) \cos(\omega t) \text{ yields:}$$ $$k C \cos(k L) \cos(\omega t) = 0 \implies \cos(k L) = 0$$ $$\text{Quantization condition: } k L = (n + 1/2)\pi \quad \text{for } n = 0, 1, 2, \dots$$ $$\text{Since } k = 2\pi f/v \text{ and } v = \sqrt{E/\rho} \text{:}$$ $$\frac{2\pi f_n}{v} L = (n + 1/2)\pi$$ $$\text{Resonance Frequencies: } f_n = \frac{(2n+1) v}{4L} = \frac{(2n+1)}{4L} \sqrt{\frac{E}{\rho}}$$