Damped harmonic oscillator revisited: a new approach to energy decay in the case of Coulomb, Stokes, and Newton damping

This paper presents a novel, simplified analytical framework for deriving accurate approximate formulas for energy decay in damped harmonic oscillators under Coulomb, Stokes, and Newton damping, while also offering an exact solution for Stokes damping that bypasses standard differential equation methods, making the approach suitable for undergraduate and advanced high school physics education.

The Gist

Imagine a child on a swing. In a perfect, frictionless world, once you give them a push, they would swing back and forth forever, reaching the exact same height every time. This is the "Harmonic Oscillator" in its purest form.

But in the real world, things slow down. Air resistance, friction on the chains, or a parent's hand gently stopping the swing all act as damping forces. They steal energy from the swing, causing the child to eventually come to a halt.

This paper by Robert Pezera and Karlo Lelas is like a new, simpler guidebook for understanding how that swing loses its energy, depending on what is slowing it down. The authors looked at three specific ways a swing can be slowed and found a clever, unified way to predict exactly how much energy is left at any given moment, without needing to solve incredibly difficult math problems.

Here is a breakdown of their three "villains" of motion and their new approach:

1. The Three Types of "Drag"

The authors compare three different ways energy is stolen from the system:

• Coulomb Damping (The "Sticky Floor"):
  Imagine the swing is moving on a rough carpet. The friction here is constant. It doesn't matter if you are moving fast or slow; the carpet drags on you with the same stubborn force.

  • The Result: The swing loses a fixed amount of height every time it swings. It's like a staircase going down. Eventually, the swing stops abruptly because the friction becomes stronger than the swing's ability to move. It doesn't fade away slowly; it just stops.

• Stokes Damping (The "Thick Honey"):
  Imagine the swing is moving through thick honey or water. The faster you go, the harder the fluid pushes back. If you double your speed, the resistance doubles.

  • The Result: This is the classic "exponential decay" taught in most physics classes. The swing loses a percentage of its energy each cycle. It gets smaller and smaller, theoretically never quite reaching zero, but getting so tiny it looks like it stopped.

• Newton Damping (The "Wind Resistance"):
  Imagine the swing is moving very fast through the air. At high speeds, air resistance doesn't just double; it quadruples (it's proportional to the square of the speed).

  • The Result: This is the most violent thief. If you push the swing hard, it loses a massive chunk of energy immediately. But once it slows down, the air resistance drops off quickly. It's a "front-loaded" energy loss: big loss at the start, then a slow fade.

2. The New "Shortcut" Method

Traditionally, figuring out exactly how a swing behaves with these forces requires solving very complex, second-order differential equations (a type of math that often scares students).

The authors say: "Why solve the whole puzzle when we can just look at the energy?"

They used a clever trick:

1. The Assumption: They assumed that even while the swing is slowing down, the ratio of its kinetic energy (speed) to its total energy still looks a lot like a perfect, undamped swing. It's like assuming a tired runner still takes steps that look mostly like a sprinter's steps, just slower.
2. The Shortcut: Instead of tracking the position of the swing second-by-second, they tracked the energy directly. They realized that the rate at which energy is lost depends on how fast the swing is moving.
3. The Magic: By using this "energy-first" perspective, they derived simple formulas that predict the swing's behavior with amazing accuracy.

3. Why This Matters (The "Aha!" Moment)

The paper isn't just about math; it's about teaching and understanding the world better.

• Simplicity for Students: You don't need to be a math wizard to understand how a swing stops. The authors show that by focusing on energy (a concept everyone understands intuitively), you can predict complex behaviors without getting bogged down in heavy calculus.
• Real-World Accuracy: Their formulas work surprisingly well, even when the damping is strong.
  • For the Sticky Floor (Coulomb), they confirmed the swing stops after a specific number of swings and calculated exactly how much energy is "frozen" in the system when it stops.
  • For the Wind (Newton), they showed that even though the math is usually a nightmare, their energy-based shortcut gives a clear picture of how the swing slows down rapidly at first and then drifts to a stop.

The Big Picture Analogy

Think of the three damping types as three different ways a bank account can be drained:

1. Coulomb (Sticky Floor): A flat fee is charged every time you make a transaction, regardless of the amount. Eventually, you run out of money and the account freezes.
2. Stokes (Honey): A percentage tax is taken every time you move money. The more you have, the more you lose, but it slows down gradually.
3. Newton (Wind): A massive penalty is charged if you move money quickly, but almost nothing if you move slowly. You lose a fortune in the first few seconds, then the losses taper off.

In summary: This paper gives us a new, simpler lens to look at how things slow down in the real world. It proves that you don't need to be a master of complex equations to understand the physics of energy loss; sometimes, just looking at the "energy budget" is enough to see the whole picture clearly.

Technical Summary

1. Problem Statement

The paper addresses the pedagogical and analytical challenges associated with modeling the Damped Harmonic Oscillator (DHO) under three distinct damping mechanisms:

1. Coulomb Damping (Constant): Friction force independent of velocity ($F_d \propto -\text{sgn}(v)$).
2. Stokes Damping (Linear): Friction force proportional to velocity ($F_d \propto -v$).
3. Newton Damping (Quadratic): Friction force proportional to the square of velocity ($F_d \propto -v|v|$).

While the Stokes (linear) case is standard in undergraduate curricula, the Coulomb and Newton cases are often treated with complex numerical methods or omitted due to the difficulty of solving their non-linear differential equations. The authors aim to derive approximate analytical formulas for energy decay in all three cases that are mathematically accessible to undergraduate students while maintaining high accuracy compared to exact or numerical solutions.

2. Methodology

The authors employ a novel approach that bypasses the standard method of solving second-order differential equations for position $x(t)$. Instead, they utilize an energy-based framework derived from fundamental calculus and physical intuition.

• Energy Dissipation Rate: They start with the rate of mechanical energy loss, $\frac{dE_m}{dt} = v \cdot F_{drag}$, where $E_m$ is total mechanical energy and $F_{drag}$ is the damping force.
• Dimensionless Formulation: The equations are converted into dimensionless forms using $\tau = \omega_0 t$ and normalized energy variables.
• The Core Approximation: The central hypothesis is that for weak damping, the ratio of kinetic energy to total mechanical energy oscillates similarly to the undamped case:
  $$\frac{E_k(\tau)}{E_m(\tau)} \approx \sin^2(\tau)$$
  This allows the authors to express the velocity-dependent damping terms solely in terms of total energy $E_m$ and the oscillatory function $\sin(\tau)$, transforming the problem into a first-order differential equation for energy.
• Separation of Variables: By substituting the approximation into the energy dissipation rate, the authors derive separable differential equations for $E_m(\tau)$ for each damping type, solving them piecewise over half-cycles to ensure continuity.

3. Key Contributions

• Unified Analytical Framework: The paper provides a single, consistent methodology to derive energy decay laws for linear, constant, and quadratic damping without solving the full second-order equation of motion.
• Exact Solution for Stokes via Ansatz: In Appendix A, the authors demonstrate a simplified method to derive the exact solution for Stokes damping using a trial function (ansatz) $x(t) = f(t)\cos(\omega_d t + \phi)$. This avoids the characteristic polynomial method typically taught in differential equations courses, making the derivation accessible to students with only basic calculus knowledge.
• Closed-Form Approximations: The authors derive explicit, closed-form expressions for the total mechanical energy envelope $E_m(\tau)$ for Coulomb and Newton damping, which were previously difficult to obtain analytically.

4. Key Results

A. Coulomb Damping (Constant Friction)

• Energy Decay: The total mechanical energy decays quadratically with time. The envelope is given by:
  $$\tilde{E}_m(\tau) = \left( 1 - \frac{2}{\pi}\gamma_0 \tau \right)^2$$
  where $\gamma_0$ is the dimensionless damping coefficient.
• Amplitude Decay: The amplitude decreases linearly with time (by a fixed amount every half-cycle).
• Stopping Time: Unlike linear damping, the system comes to a complete stop in finite time. The authors derive a precise condition for the stopping time based on the ratio of static to kinetic friction coefficients.

B. Stokes Damping (Linear Friction)

• Energy Decay: The energy decays exponentially, consistent with standard theory:
  $$E_m(\tau) \approx e^{-\gamma_1 \tau}$$
  The authors show that the exact solution includes a small oscillatory modulation term, which their approximation captures effectively in the weak damping limit.

C. Newton Damping (Quadratic Friction)

• Energy Decay: The total mechanical energy decays according to an inverse square law:
  $$\tilde{E}_m(\tau) = \left( 1 + \frac{4}{3\pi}\gamma_2 \tau \right)^{-2}$$
• Behavior: The dissipation is strongest at high velocities (large amplitudes). As the amplitude decreases, the system effectively transitions into a weak damping regime.
• Accuracy: The approximation remains valid even for moderate damping strengths ($\gamma_2 \approx 0.6$), provided the system undergoes several oscillations.

5. Significance and Pedagogical Impact

• Undergraduate Accessibility: The methods used rely on elementary derivatives, integration, and physical reasoning rather than advanced differential equation techniques. This makes the complex physics of non-linear damping accessible to introductory physics and engineering students.
• Conceptual Clarity: By focusing on energy rather than displacement, the paper highlights the physical mechanism of dissipation (work done by friction) more intuitively. It demonstrates that exponential decay is not universal; different damping laws lead to polynomial or inverse-square decay.
• Experimental Validation: The derived formulas are highly accurate when compared to exact numerical simulations. The authors suggest these formulas are ideal for designing laboratory exercises where students can use motion tracking software to verify energy decay laws for different friction types (e.g., sliding blocks vs. air resistance).
• Theoretical Insight: The paper illustrates the interplay between theory and experiment, showing how a theoretical model derived for weak damping can predict behavior in regimes (like critical or overdamped Stokes) that are difficult to observe experimentally without prior theoretical guidance.

In conclusion, Pezer and Lelas provide a robust, simplified, and highly accurate analytical toolkit for understanding damped harmonic motion across different physical regimes, significantly lowering the mathematical barrier to entry for studying non-linear damping in physics education.