Differential Equations in Oscillations & Circuits - College Physics

SECTION 1

Historical Context & Motivation

The story of differential equations in oscillatory systems stretches back to the dawn of modern mechanics, when natural philosophers first attempted to describe repetitive motion in precise mathematical language. Robert Hooke's empirical observation that restoring forces are proportional to displacement, combined with Newton's second law, produced one of the earliest physically motivated differential equations. Over the next two centuries, mathematicians such as Euler, Lagrange, and Fourier developed the analytical tools needed to solve these equations in progressively more complex settings. The true power of the framework became evident when James Clerk Maxwell and later physicists realized that the same second-order ODE governing a mass on a spring also governs charge oscillations in an electrical circuit. This mathematical unification remains one of the most elegant results in classical physics and continues to underpin modern engineering disciplines from acoustics to telecommunications.

1676 • Hooke's Law

Robert Hooke publishes ut tensio, sic vis ("as the extension, so the force"), establishing the linear restoring-force model that leads directly to the simple harmonic oscillator equation.

1687 • Newton's Principia

Isaac Newton's second law, F = ma, provides the dynamical framework that converts Hooke's force law into the second-order differential equation m d²x/dt² = −kx.

1855 • Lord Kelvin's LC Circuit Analogy

William Thomson (Lord Kelvin) predicts that a capacitor discharging through an inductor will produce oscillating current, drawing an explicit analogy to a mechanical oscillator and sparking the study of electrical oscillations.

1891 • Hertz & Tesla: AC Circuits

The practical exploitation of oscillatory circuits accelerates with Hertz's demonstration of electromagnetic waves and Tesla's development of alternating-current systems, both relying on solutions to driven RLC differential equations.

1920s • Quantum Harmonic Oscillator

Schrödinger and Heisenberg extend the classical harmonic oscillator to quantum mechanics, showing that quantized energy levels emerge from the same differential equation with modified boundary conditions.

The central question this lesson addresses is deceptively simple: why do so many physically distinct systems — springs, pendulums, and RLC circuits — obey the same differential equation, and how do we solve it? Understanding this unifying structure will equip you with a single mathematical toolkit that applies across mechanics, electromagnetism, and beyond.

SECTION 2

Core Principles & Definitions

Before diving into solutions, it is essential to establish the foundational ideas that connect oscillatory mechanics to circuit theory. A second-order linear ordinary differential equation (ODE) is the mathematical object at the heart of both domains. In mechanics, Newton's second law combined with a restoring force produces such an equation; in circuits, Kirchhoff's voltage law applied around a loop with inductors, resistors, and capacitors yields an identical form. The key physical ingredients — inertia, dissipation, and a restoring mechanism — map one-to-one between the two systems, making the analogy far more than a pedagogical convenience: it is a structural isomorphism.

Simple Harmonic Motion (SHM)

Oscillation about equilibrium with no energy loss. The displacement satisfies d²x/dt² + ω₀²x = 0, producing sinusoidal solutions of constant amplitude at the natural frequency ω₀.

Damped Oscillation

A dissipative term (proportional to velocity in mechanics, or resistance in circuits) introduces exponential decay. Three regimes arise depending on the damping ratio ζ: underdamped, critically damped, and overdamped.

Driven (Forced) Oscillation

An external periodic force or voltage source adds a non-homogeneous term. The steady-state response exhibits resonance when the driving frequency matches the natural frequency, dramatically amplifying the response.

Mechanical–Electrical Analogy

Mass ↔ inductance, damping coefficient ↔ resistance, spring constant ↔ 1/capacitance, displacement ↔ charge. This mapping allows any circuit problem to be rephrased as a mechanics problem and vice versa.

✦

KEY TAKEAWAY

Think of a second-order ODE as a recipe with three ingredients: something that resists changes in motion (inertia/inductance), something that bleeds energy (friction/resistance), and something that pulls back toward equilibrium (spring/capacitor). Whether the oscillating quantity is a position coordinate or an electrical charge, the recipe — and therefore the solution — is identical. This is analogous to how a single blueprint can describe buildings made of wood or steel: the materials differ, but the structural logic does not.

SECTION 3

Visual Explanation — The Mechanical–Electrical Parallel

Left: a mass–spring–damper system with mass m, spring constant k, and damping coefficient b. Right: a series RLC circuit with inductor L, resistor R, and capacitor C. The mapping table at bottom shows the one-to-one correspondence between mechanical and electrical quantities.

The diagram above illustrates the deep structural parallel between a mass–spring–damper system and a series RLC circuit. In the mechanical system, the mass provides inertia (it resists acceleration), the dashpot dissipates energy through viscous friction, and the spring supplies a restoring force proportional to displacement. In the electrical system, the inductor resists changes in current (analogous to inertia), the resistor dissipates energy as heat, and the capacitor stores energy in its electric field while providing a restoring "force" that opposes charge accumulation. When you write Newton's second law for the mechanical system and Kirchhoff's voltage law for the circuit, you obtain differential equations of identical form — a fact captured in the mapping table at the bottom of the figure.

SECTION 4

Mathematical Framework

We now formalize the differential equations that govern both mechanical oscillators and electrical circuits. The derivation proceeds from first principles — Newton's second law for mechanics and Kirchhoff's voltage law for circuits — and arrives at a single canonical form. The solution technique involves the characteristic equation, a quadratic whose roots dictate the qualitative behavior of the system.

Mechanical Oscillator

DAMPED HARMONIC OSCILLATOR (MECHANICS)

m d²x/dt² + b dx/dt + kx = F(t)

where m = mass (kg), b = damping coefficient (kg/s), k = spring constant (N/m), F(t) = external driving force.

Series RLC Circuit

DAMPED OSCILLATOR (RLC CIRCUIT)

L d²q/dt² + R dq/dt + q/C = V(t)

where L = inductance (H), R = resistance (Ω), C = capacitance (F), q = charge on the capacitor, V(t) = applied EMF.

Canonical Form & Characteristic Equation

Dividing the mechanical equation by m and the circuit equation by L, both can be cast into the standard form shown below, where ω₀ is the natural angular frequency and γ is the damping parameter.

STANDARD FORM

d²x/dt² + 2γ dx/dt + ω₀² x = f(t)

Mechanical: ω₀ = √(k/m), γ = b/(2m). Circuit: ω₀ = 1/√(LC), γ = R/(2L). The damping ratio is ζ = γ/ω₀.

CHARACTERISTIC EQUATION

r² + 2γr + ω₀² = 0 → r = −γ ± √(γ² − ω₀²)

The discriminant Δ = γ² − ω₀² determines the regime: Δ < 0 → underdamped (oscillatory decay), Δ = 0 → critically damped (fastest non-oscillatory return), Δ > 0 → overdamped (sluggish exponential return).

💡

Why the Characteristic Equation Works

We assume a trial solution of the form x(t) = Ae^rt. Substituting into the homogeneous ODE turns every derivative into a power of r, reducing the differential equation to an algebraic (quadratic) equation. The roots r₁ and r₂ then dictate whether the solution is oscillatory (complex roots), exponentially decaying (equal real roots), or a sum of two distinct exponentials (distinct real roots). This is the ansatz method, one of the most powerful techniques in the theory of linear ODEs.

SECTION 5

Damping Regimes in Detail

The character of the solution to the damped oscillator equation changes qualitatively depending on the value of the damping ratio ζ = γ/ω₀. This section presents all three regimes — underdamped, critically damped, and overdamped — alongside their explicit solutions and a visual comparison of the resulting waveforms. Physically, these regimes correspond to very different behaviors: a bell that rings (underdamped), a door closer that shuts smoothly (critically damped), and honey on a spoon that oozes without oscillating (overdamped). Recognizing which regime applies is crucial in both physics and engineering design.

Displacement x(t) versus time for the three damping regimes, all starting from the same initial displacement x₀ with zero initial velocity. The underdamped curve (cyan) oscillates within an exponentially decaying envelope. The critically damped curve (amber) returns to equilibrium in the shortest time without overshooting. The overdamped curve (red) decays more slowly, approaching zero asymptotically.

Regime	Condition	General Solution (homogeneous)
Underdamped	ζ < 1 (γ < ω₀)	x(t) = e^−γt [A cos(ω_dt) + B sin(ω_dt)], ω_d = √(ω₀² − γ²)
Critically Damped	ζ = 1 (γ = ω₀)	x(t) = (A + Bt) e^−γt
Overdamped	ζ > 1 (γ > ω₀)	x(t) = A e^r₁t + B e^r₂t, r₁,₂ = −γ ± √(γ² − ω₀²)

Solutions for each damping regime of the homogeneous second-order ODE.

🔔

PHYSICAL INTUITION

Imagine striking a bell inside a room filled with air (underdamped — you hear a tone that fades), then filling the room with thick honey (overdamped — the bell barely moves). Critical damping is the boundary: just enough viscosity to suppress oscillation while allowing the fastest possible return to rest. Engineers design car shock absorbers and galvanometer needles to operate at or near critical damping for exactly this reason.

SECTION 6

Worked Example — Underdamped RLC Circuit

Consider a series RLC circuit with L = 0.50 H, R = 20 Ω, and C = 50 μF. The capacitor is initially charged to q₀ = 100 μC and there is no initial current. We wish to determine: (a) the damping regime, (b) the oscillation frequency, and (c) the charge on the capacitor as a function of time.

Underdamped RLC Discharge

1
Step 1 — Identify Given Values and Relevant Equation

We have L = 0.50 H, R = 20 Ω, C = 50 × 10⁻⁶ F, q(0) = 100 × 10⁻⁶ C, and i(0) = dq/dt|_t=0 = 0. The governing equation is L d²q/dt² + R dq/dt + q/C = 0 (homogeneous, no external source).
2
Step 2 — Compute Natural Frequency ω₀ and Damping Parameter γ

ω₀ = 1/√(LC) = 1/√(0.50 × 50 × 10⁻⁶) = 1/√(2.5 × 10⁻⁵) = 1/(5.0 × 10⁻³) = 200 rad/s. γ = R/(2L) = 20/(2 × 0.50) = 20 s⁻¹.

ω₀ = 200 rad/s, γ = 20 s⁻¹
3
Step 3 — Determine the Damping Regime

Since γ = 20 s⁻¹ < ω₀ = 200 rad/s, we have ζ = γ/ω₀ = 0.10 < 1. The circuit is underdamped, so the solution will be oscillatory with exponential decay.

Underdamped (ζ = 0.10)
4
Step 4 — Compute the Damped Frequency ω_d

ω_d = √(ω₀² − γ²) = √(200² − 20²) = √(40000 − 400) = √39600 ≈ 199 rad/s. Note that because ζ is small, the damped frequency is very close to the natural frequency.

ω_d ≈ 199 rad/s (f_d ≈ 31.7 Hz)
5
Step 5 — Write the General Solution and Apply Initial Conditions

The general underdamped solution is q(t) = e^−γt [A cos(ω_dt) + B sin(ω_dt)]. IC 1: q(0) = A = 100 μC. IC 2: dq/dt|_t=0 = −γA + ω_dB = 0, so B = γA/ω_d = (20)(100 μC)/199 ≈ 10.1 μC.

q(t) = $e^{−20t}$ [100 cos(199t) + 10.1 sin(199t)] μC

SECTION 7

Strengths & Limitations of the Linear Model

The linear second-order ODE model is remarkably versatile, but it rests on assumptions that may break down in real-world situations. Understanding where the model excels and where it fails is essential for applying it responsibly in both physics and engineering contexts.

Strengths	Limitations
Exact closed-form solutions exist for all three damping regimes, making analysis fast and transparent.	Linearity assumes small displacements or small signals; large-amplitude oscillations (e.g., a pendulum at 90°) require nonlinear corrections.
The mechanical–electrical analogy lets engineers prototype circuits with mechanical models and vice versa.	Real components are frequency-dependent — inductor resistance, capacitor leakage, and nonlinear spring behavior are not captured.
Superposition applies: the response to a complex input can be decomposed into simpler parts via Fourier analysis.	Superposition fails for nonlinear oscillators (e.g., Duffing oscillator, Van der Pol oscillator), which exhibit chaos and limit cycles.
Resonance phenomena are predicted quantitatively — critical for radio tuning, MRI, and structural engineering.	Near resonance, the model may predict infinite amplitudes (undamped driven case), requiring damping or nonlinear saturation for realism.

Strengths and limitations of the linear second-order ODE model for oscillations and circuits.

⚖️

KEEP IN PERSPECTIVE

The linear ODE model is physics' equivalent of a map at 1:100,000 scale — extraordinarily useful for navigation but unable to show every pothole. For most undergraduate-level problems and a wide range of engineering applications, the linear approximation is accurate and sufficient. When it isn't — in high-power circuits, large-angle pendulums, or chaotic systems — you will reach for perturbation theory, numerical solvers, or entirely different mathematical frameworks.

SECTION 8

Connection to Advanced Theory

The second-order ODE you have mastered in this lesson is a gateway to several advanced topics. In upper-division courses and graduate work, the same equation reappears in new guises — sometimes with additional spatial dimensions, sometimes in the complex-frequency domain, and sometimes quantized. The table below maps key extensions from your current knowledge to their advanced counterparts.

This Lesson (Classical, Linear)	Advanced Extension
Damped harmonic oscillator (one mass, one spring)	Coupled oscillators & normal modes — N masses on springs lead to N coupled ODEs, diagonalized via matrix methods.
Lumped-element RLC circuit	Transmission lines & wave equations — distributed L and C per unit length produce a partial differential equation (the telegraph equation).
Time-domain ODE: d²x/dt² + 2γ dx/dt + ω₀²x = f(t)	Laplace & Fourier transforms — algebraize the ODE in the s-domain or ω-domain, turning differentiation into multiplication.
Classical harmonic oscillator (continuous energy)	Quantum harmonic oscillator — same potential V = ½kx², but Schrödinger's equation yields quantized energy levels Eₙ = (n + ½)ℏω₀.
Resonance at ω = ω₀	Parametric resonance & Floquet theory — the spring constant itself varies periodically, leading to exponential growth under specific conditions.

How the topics of this lesson connect to advanced coursework.

One particularly elegant bridge is the impedance formalism taught in AC circuit analysis. By writing voltages and currents as complex exponentials V(t) = V₀ e^jωt, the differential equation reduces to a simple algebraic equation involving complex impedances Z_L = jωL, Z_C = 1/(jωC), and Z_R = R. This phasor approach is essentially a special case of the Fourier transform method and allows you to solve driven RLC problems without ever explicitly solving the ODE — a tremendous computational advantage.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

A series RLC circuit and a mass–spring–damper system both have ζ = 0.5. Without solving any equations, describe the qualitative behavior of both systems after they are displaced from equilibrium and released. Why is the behavior the same despite the physical systems being completely different?

Show answer

Both systems are underdamped (ζ < 1), so they will oscillate about equilibrium with an amplitude that decreases exponentially over time. The oscillation frequency is slightly lower than the natural frequency due to the damping. The behavior is identical because both systems obey the same canonical second-order ODE: d²x/dt² + 2γ dx/dt + ω₀²x = 0. The mathematical structure — not the physical medium — determines the qualitative dynamics, which is the essence of the mechanical–electrical analogy.

Question 2: BASIC CALCULATION

A mass of 0.25 kg is attached to a spring with spring constant k = 100 N/m. There is no damping. Calculate the natural frequency ω₀ (in rad/s) and the period T of oscillation.

Show answer

The natural angular frequency is ω₀ = √(k/m) = √(100/0.25) = √400 = 20 rad/s. The period is T = 2π/ω₀ = 2π/20 ≈ 0.314 s. This is a simple application of the undamped SHM equation d²x/dt² + ω₀²x = 0, where the solution is x(t) = A cos(ω₀t + φ) with angular frequency 20 rad/s and period approximately 0.31 seconds.

Question 3: INTERMEDIATE

A series RLC circuit has L = 10 mH, C = 1 μF, and R = 10 Ω. (a) Find ω₀, γ, and ζ. (b) Determine the damping regime. (c) If R is increased to 400 Ω, what is the new damping regime?

Show answer

(a) ω₀ = 1/√(LC) = 1/√(10 × 10⁻³ × 1 × 10⁻⁶) = 1/√(10⁻⁸) = 10⁴ rad/s = 10,000 rad/s. γ = R/(2L) = 10/(2 × 0.01) = 500 s⁻¹. ζ = γ/ω₀ = 500/10,000 = 0.05. (b) Since ζ = 0.05 < 1, the circuit is underdamped and will exhibit oscillatory discharge with very slow decay. (c) With R = 400 Ω: γ = 400/(0.02) = 20,000 s⁻¹, ζ = 20,000/10,000 = 2 > 1, so the circuit becomes overdamped. The crossover to critical damping occurs at R = 2√(L/C) = 2√(10⁻² / 10⁻⁶) = 2 × 100 = 200 Ω.

Question 4: APPLIED

An automotive shock absorber must be critically damped to prevent passenger discomfort. The car's suspension acts as a spring with k = 40,000 N/m, and the effective mass per wheel is 400 kg. What damping coefficient b (in kg/s) is required for critical damping?

Show answer

For critical damping, ζ = 1, which requires γ = ω₀. We compute ω₀ = √(k/m) = √(40,000/400) = √100 = 10 rad/s. The damping parameter is γ = b/(2m), so b = 2mγ = 2mω₀ = 2(400)(10) = 8,000 kg/s. This ensures the suspension returns to equilibrium as quickly as possible without oscillation — exactly the compromise between comfort and stability that automotive engineers target. In practice, slight underdamping (ζ ≈ 0.7–0.9) is sometimes preferred for ride comfort.

Question 5: CRITICAL THINKING

Consider an undamped LC circuit (R = 0) driven by an external voltage V(t) = V₀ cos(ωt), where ω is exactly equal to the natural frequency ω₀ = 1/√(LC). The particular solution to the driven ODE grows without bound. (a) Show mathematically why this happens. (b) Explain why infinite charge is physically unrealizable. (c) What role does even a small resistance play in resolving this paradox?

Show answer

(a) The ODE is L d²q/dt² + q/C = V₀ cos(ω₀t). Dividing by L: d²q/dt² + ω₀²q = (V₀/L) cos(ω₀t). Since the driving frequency matches the natural frequency, the standard particular solution A cos(ω₀t) fails (it solves the homogeneous equation). We must use the modified ansatz

q_p

(t) = (V₀t)/(2Lω₀) sin(ω₀t), which grows linearly in time — this is unbounded resonance. (b) In reality, infinite charge is impossible because: (i) the capacitor will undergo dielectric breakdown at high voltage, (ii) component nonlinearities increase with amplitude, and (iii) radiation losses occur at high oscillation amplitudes. (c) Even a tiny resistance R adds a 2γ dq/dt term, converting the equation into a damped driven oscillator. The steady-state amplitude becomes V₀/(R ω₀) — finite and inversely proportional to R. This shows that damping limits the resonance peak; the quality factor Q = ω₀/(2γ) quantifies how sharp and tall that peak is.

SUMMARY

Lesson Summary

Summary

This lesson demonstrated that mechanical oscillators (mass–spring–damper systems) and series RLC circuits are governed by the same second-order linear ODE: d²x/dt² + 2γ dx/dt + ω₀²x = f(t). The mechanical–electrical analogy maps mass to inductance, damping coefficient to resistance, and spring constant to inverse capacitance, while displacement corresponds to charge. Solving the characteristic equation r² + 2γr + ω₀² = 0 reveals three distinct damping regimes — underdamped (oscillatory decay), critically damped (fastest non-oscillatory return), and overdamped (sluggish exponential decay) — classified by the damping ratio ζ.

The driven oscillator introduces the phenomenon of resonance, where the response amplitude peaks when the driving frequency matches ω₀. This framework extends naturally to coupled oscillators, transmission lines, Laplace transforms, and the quantum harmonic oscillator. Mastering the second-order ODE equips you with a universal tool: once you can solve it for a spring, you can solve it for any linear system in physics and engineering that exhibits inertia, dissipation, and a restoring mechanism.