The Error in Rayleigh's Approximative Period

This paper establishes rigorous a priori bounds for the exact period of Rayleigh's stretched string equation to prove that Rayleigh's approximation overestimates the true period, providing explicit inequalities that demonstrate the relative error is directly proportional to the initial displacement and inversely proportional to the initial stretch.

The Gist

Imagine you have a guitar string stretched tight between two points. Now, imagine tying a heavy weight right in the middle of that string. If you pull the weight down and let it go, it will bob up and down like a pendulum.

For over a century, scientists have used a famous, simple formula by Lord Rayleigh to guess how long it takes for that weight to complete one full bounce (its "period"). Rayleigh's logic was: "If the bounce is small, the string doesn't stretch much more than it already is, so let's just pretend the tension is constant."

It's a handy shortcut, like using a flat map to navigate a small neighborhood. But as Mark Villarino points out in this paper, flat maps lie when you try to navigate a mountain.

Here is what this paper actually does, translated into everyday language:

1. The Problem: The "Flat Map" vs. The "Mountain"

Rayleigh's formula assumes the string stays perfectly straight and tight, ignoring the fact that as the weight moves, the string actually stretches a tiny bit more, changing the tension.

• Rayleigh's View: The string is a rigid rod.
• Reality: The string is a rubber band that gets tighter the more you pull it.

Because Rayleigh ignored this extra stretching, his formula always predicts the bounce will be slower than it actually is. It's like guessing a car trip will take an hour because you ignored the traffic; in reality, you might get there faster if the road was clear.

2. The Solution: Drawing a Safety Net

Villarino didn't just say, "Rayleigh is wrong." He did something much more useful: He built a safety net.

He proved mathematically that Rayleigh's answer is always too high, but he also calculated exactly how much too high it could be. He created two invisible walls:

• The Lower Wall: The absolute fastest the weight could possibly bounce.
• The Upper Wall: The absolute slowest it could possibly bounce (which is Rayleigh's answer).

The true answer is guaranteed to be somewhere between these two walls.

3. The "Surprise" Ingredients

You might think the error depends on how heavy the weight is or what the string is made of (steel vs. copper). Villarino found something surprising: It doesn't.

The error depends entirely on geometry:

1. How much was the string stretched initially? (Is it loose or tight?)
2. How far did you pull the weight down?

Think of it like a trampoline. If you pull the mat down a little bit on a trampoline that is already pulled tight, the error is tiny. But if the trampoline is barely stretched (loose) and you pull it down a lot, the physics get weird, and the simple formula fails spectacularly.

4. The "Aha!" Moment: When Rayleigh Fails

The paper includes a scary example to prove its point.

• Scenario: A steel wire is pulled tight with a huge force, but the extra stretch caused by the weight is microscopic (almost zero).
• Rayleigh's Prediction: The bounce takes 0.22 seconds.
• The Reality: The bounce takes 0.18 seconds.
• The Result: Rayleigh was off by 25%.

Why? Because when the initial stretch is tiny, the "extra stretch" caused by the movement becomes the dominant factor. Rayleigh's assumption that "stretch is negligible" collapses like a house of cards.

The Takeaway

This paper is like a mechanic telling you: "That old shortcut you use to fix your car? It works fine for a smooth highway, but if you hit a pothole, it breaks. Here is a new rule that tells you exactly how bad the pothole will be, so you know when to stop using the shortcut."

Villarino gives us a precise way to measure the "lie" in Rayleigh's famous formula, ensuring that engineers and scientists know exactly when they can trust the old math and when they need to do the hard work of the real calculation.

Technical Summary

Here is a detailed technical summary of the paper "The Error in Rayleigh's Approximative Period" by Mark B. Villarino.

1. Problem Statement

The paper addresses a classic problem in mathematical physics involving the vertical oscillation of a mass $m$ attached to the center of a horizontal elastic wire.

• Physical Setup: A wire of natural length $2L_0$and spring constant$\sigma$is stretched to a length$2L$. A mass$m$(where$m$is much greater than the wire's mass) is attached to the center and displaced vertically by$y(t)$. Gravity and damping are neglected.
• The Governing Equation: Using Hooke's Law and Newton's Second Law, the exact equation of motion is derived as a nonlinear differential equation:
  $$\ddot{y} + \frac{2\sigma}{m} \left( \frac{\sqrt{L^2 + y^2} - L_0}{\sqrt{L^2 + y^2}} \right) y = 0$$
  This equation is complex and lacks an exact analytical solution in terms of elementary functions.
• Rayleigh's Approximation: Lord Rayleigh proposed a simplification by assuming that for small vibrations, the additional stretching is negligible, allowing the tension $T$ to be treated as a constant ($T = \sigma(L - L_0)$). This reduces the system to a simple harmonic oscillator:
  $$\ddot{y} + \frac{2T}{mL}y = 0$$
  yielding an approximate period $P = 2\pi \sqrt{\frac{mL}{2T}}$.
• The Gap: While Rayleigh's approximation is widely used, the literature lacks a rigorous error analysis quantifying how close the approximate period $P$ is to the true period $\mathcal{P}$.

2. Methodology

Villarino employs a rigorous mathematical approach to bound the error without relying on numerical simulation alone.

• Derivation of the Exact Period: By transforming the second-order nonlinear ODE into a separable first-order equation for velocity ($v = \dot{y}$), the author derives an integral formula for the true period $\mathcal{P}$. This integral involves a radical function and is identified as a complicated elliptic integral that cannot be computed exactly using elementary functions.
• Bounding the Integrand: Instead of solving the integral, the author establishes upper and lower bounds for the integrand's radical term. By bounding the term $\frac{1}{\sqrt{L_0} - \frac{2}{\sqrt{L^2+y^2} + \sqrt{L^2+y_0^2}}}$ within the integral, the author derives strict inequalities for the true period $\mathcal{P}$.
• Relative Error Analysis: The author defines the relative error as $R = \frac{P - \mathcal{P}}{\mathcal{P}}$. By comparing the bounds of the true period $\mathcal{P}$ with Rayleigh's approximate period $P$, the paper derives explicit inequalities for $R$.
• Series Expansion: To simplify the complex radical bounds into a usable formula, the author utilizes Taylor series expansions (specifically for $\sqrt{1+u}$) to approximate the behavior of the error term for small displacements.

3. Key Contributions and Results

The paper provides the first rigorous a priori bounds for Rayleigh's approximation, yielding three main theorems:

• Theorem 1 (Direction and Bounds of Error):

  • Rayleigh's approximation always overestimates the true period ($P > \mathcal{P}$, meaning $R < 0$).
  • The relative error is bounded by:
    $$\sqrt{\frac{\frac{1}{L_0} - \frac{2}{L+\sqrt{L^2+y_0^2}}}{\frac{1}{L_0} - \frac{1}{L}}} < \frac{P}{\mathcal{P}} < \sqrt{\frac{\frac{1}{L_0} - \frac{1}{\sqrt{L^2+y_0^2}}}{\frac{1}{L_0} - \frac{1}{L}}}$$
  • Key Insight: The error depends only on the geometric ratios $L_0/L$ and $y_0/L$. It is independent of the mass $m$ and the material's Young's modulus (provided Hooke's law holds).

• Theorem 2 (Explicit Error Formula):
  The absolute relative error $|R|$ is shown to be directly proportional to the square of the initial fractional displacement and inversely proportional to the initial stretch:
  $$|R| = \Theta \cdot \frac{L_0}{L - L_0} \cdot \left( \frac{y_0}{2L} \right)^2$$
  where the proportionality factor $\Theta$ is bounded by:
  $$\frac{1}{2 - (y_0/2L)^2} \cdot \frac{1}{1 + \frac{1}{32}\frac{L_0}{L-L_0}} \leq \Theta \leq 1$$

• Theorem 3 (Exact Period Formula):
  The paper provides the exact integral representation for the period, confirming it is a non-elementary elliptic integral:
  $$\mathcal{P} = 4 \sqrt{\frac{m}{2\sigma L_0}} \int_0^{y_0} \frac{1}{\sqrt{y_0^2 - y^2}} \sqrt{\frac{1}{L_0} - \frac{2}{\sqrt{L^2+y^2} + \sqrt{L^2+y_0^2}}} \, dy$$

4. Significance and Examples

The paper demonstrates that Rayleigh's approximation can be dangerously inaccurate under specific conditions, specifically when the initial stretch ($L - L_0$) is very small relative to the displacement.

• Example 1 (Valid Approximation): For a wire stretched from 1152mm to 1250mm with a moderate displacement, the relative error was found to be 3.4%, which fell comfortably within the theoretical bounds (2.5% to 5%).
• Example 2 (Failure of Approximation): For a steel wire with a very high spring constant (resulting in a tiny stretch of 0.00064m) and a 2cm displacement, Rayleigh's approximation yielded a 25% error.
  • Explanation: The error term contains the factor $\frac{L_0}{L-L_0}$. When the stretch ($L-L_0$) is tiny, this factor becomes huge, causing the approximation to fail. The theoretical lower bound correctly predicted an error of at least 15.7%, explaining the large deviation.

Conclusion

Villarino's work fills a significant gap in classical mechanics literature by providing rigorous error bounds for a famous approximation. The results show that while Rayleigh's formula is useful for large stretches, it is unreliable for stiff wires with small initial tensions. The derived formulas allow engineers and physicists to predict the accuracy of the approximation based solely on geometric ratios before performing calculations.