Physics C: Mechanics • Score 5 Strategy

Work-Energy Theorem with Velocity-Dependent Drag Forces Guide: AP Physics C: Mechanics Score 5 for Harvard University

AP Physics C: Mechanics Mastery Guide

Work-Energy Theorem with Velocity-Dependent Drag Forces


1. Introduction & AP Exam Weight

The interaction between conservative field forces and velocity-dependent non-conservative drag forces represents one of the most sophisticated, high-yielding conceptual frameworks on the AP Physics C: Mechanics Exam. While introductory physics relies on constant acceleration models, the College Board heavily favors non-linear differential equations where force vector fields depend dynamically on system velocity $\vec{v}(t)$.

This topic typically constitutes 15–23% of the Free-Response Question (FRQ) section, frequently appearing in FRQ 1 or FRQ 2 as a hybrid dynamic-energy calculus question.

Core Conceptual Scope

Mastery of this domain demands transitioning from algebraic energy equations to dynamic variable-separable differential integrals—a hallmark skill distinguishing top-tier applicants to elite physical science and engineering programs.


2. Deep Concept Breakdown

Theoretical Framework & Differential Mechanics

The Work-Energy Theorem asserts that the net work performed on a particle equals the change in its kinetic energy:

$$W_{\text{net}} = \Delta K = K_f - K_i$$

When non-conservative forces like velocity-dependent drag act on an object, the total work splits into conservative ($W_c$) and non-conservative ($W_{\text{nc}}$) components:

$$W_c + W_{\text{nc}} = \Delta K \implies W_{\text{nc}} = \Delta K - W_c = \Delta E_{\text{mech}}$$

Consider a particle of mass $m$ subject to a constant gravitational force and a linear velocity-dependent drag force $\vec{F}_d = -b\vec{v}$. If the particle is launched vertically downward with initial velocity $v_0$ from $y = 0$:

$$\sum F_y = mg - bv = m a = m v \frac{dv}{dy}$$

To find velocity as a function of position $y$, we establish a variable-separable differential equation:

$$m v \frac{dv}{dy} = mg - bv \implies \frac{m v}{mg - bv} \, dv = dy$$

Derivation of Spatial Velocity $v(y)$ via Dynamic Integration

$$\int_{v_0}^{v(y)} \frac{m v}{mg - bv} \, dv = \int_0^y dy$$

Using the substitution $u = mg - bv \implies du = -b \, dv$ and $v = \frac{mg - u}{b}$:

$$\int_{u_0}^{u} \frac{m \left(\frac{mg - u}{b}\right)}{u} \left(-\frac{du}{b}\right) = -\frac{m}{b^2} \int_{u_0}^{u} \left( \frac{mg}{u} - 1 \right) du = y$$

$$-\frac{m}{b^2} \left[ mg \ln\left(\frac{u}{u_0}\right) - (u - u_0) \right] = y$$

Substituting back $u = mg - bv$ and $u_0 = mg - bv_0$:

$$-\frac{m}{b^2} \left[ mg \ln\left(\frac{mg - bv}{mg - bv_0}\right) + b(v - v_0) \right] = y$$

Exact Non-Conservative Work Integral

The work done by the drag force as a function of position $y$ can be expressed either spatially or temporally:

$$W_{\text{drag}} = \int_0^y \vec{F}_d \cdot d\vec{y} = \int_0^y (-bv) \, dy$$

Since $dy = v \, dt$:

$$W_{\text{drag}} = -b \int_0^t v(t)^2 \, dt$$

Alternatively, invoking the Work-Energy Theorem directly:

$$W_{\text{drag}} = \Delta K - W_g = \frac{1}{2}m\left(v(y)^2 - v_0^2\right) - mgy$$


Computational Simulation (Python)

To validate non-conservative work dissipation and energy conservation numerically, high-achieving students utilize numerical integration schemes such as RK2 (Midpoint Method) or RK4. Below is a production-grade Python script modeling velocity-dependent drag work dissipation.

import numpy as np

def simulate_drag_energy(m: float, g: float, b: float, v0: float, y_target: float, dt: float = 1e-5):
    """
    Numerically integrates dynamic equations of motion with linear drag v(y) 
    and verifies the Work-Energy Theorem: W_nc = Delta K - W_g.
    """
    # State variables
    y = 0.0
    v = v0
    t = 0.0

    # Energy terms
    K_initial = 0.5 * m * v0**2
    W_drag = 0.0

    while y < y_target:
        # Dynamic acceleration
        a = g - (b / m) * v

        # Incremental work done by drag over displacement dy = v * dt
        dy = v * dt
        F_drag = -b * v
        dW_drag = F_drag * dy

        # State update (Euler-Cromer Integration for stability)
        v += a * dt
        y += dy
        t += dt
        W_drag += dW_drag

    K_final = 0.5 * m * v**2
    Delta_K = K_final - K_initial
    W_gravity = m * g * y

    # Verification of Work-Energy balance
    E_error = abs(Delta_K - (W_gravity + W_drag))

    print(f"--- Numerical Results at y = {y:.3f} m ---")
    print(f"Final Velocity (v): {v:.4f} m/s")
    print(f"Delta K:            {Delta_K:.4f} J")
    print(f"Gravity Work (Wg):  {W_gravity:.4f} J")
    print(f"Drag Work (W_drag): {W_drag:.4f} J")
    print(f"Absolute Residual:  {E_error:.8e} J")

if __name__ == "__main__":
    # Parameters: m=1.5kg, g=9.8m/s^2, b=0.45kg/s, v0=2.0m/s, y_target=10.0m
    simulate_drag_energy(m=1.5, g=9.8, b=0.45, v0=2.0, y_target=10.0)

3. Common AP Exam Pitfalls & Score 5 Scoring Rubric Nuances

Score 4 vs. Score 5 Performance Contrast

Topic Score 4 Student Approach Score 5 Elite Student Approach
Work Integrals with Drag Assumes force is static and writes $W_d = F_d \cdot d = (-bv)d$, attempting to plug in average velocity $\bar{v}$. Treats drag work strictly as a dynamic integral $W_d = -\int b v^2 dt$ or isolates it via $W_{\text{nc}} = \Delta K - W_c$.
Differential Setup Writes $a = \frac{dv}{dt}$ when asked for position-dependent velocity, getting stuck in time-domain variable separation. Immediately substitutes $a = v\frac{dv}{dx}$, setting up separable integrals with spatial boundary conditions ($x_0 \to x_f$).
Terminal Limits Forgets that as $t \to \infty$, $a \to 0$, leading to incorrect asymptotic integration boundaries. Explicitly solves for $v_T = \frac{mg}{b}$ prior to integration to define explicit integral convergent bounds.
Vector Calculus Alignment Neglects dot product signs; explicitly double-counts negative signs, yielding positive work for dissipative drag forces. Formally evaluates $\vec{F}_d \cdot d\vec{r} = F_d \, dr \cos(180^\circ) = -F_d \, dr$, ensuring work is strictly negative.

Scoring Rubric Nuances (AP FRQ Analysis)

Consider an official AP Physics C FRQ prompt involving non-constant forces:

AP Physics C AP Rubric Breakdown (3-Point Sub-part)

  1. [1 Point] Differential Equation Setup:
  2. Criterion: Writes a valid Newton's 2nd Law equation incorporating dynamic derivative definitions.
  3. Score 5 Execution: Explicitly writes $m v \frac{dv}{dx} = -b v^2 - \mu m g$. Earns point instantly.
  4. [1 Point] Separation of Variables & Limits:
  5. Criterion: Correctly separates variables with proper integration limits matching state variables.
  6. Score 5 Execution: Writes $\int_{v_0}^{v_f} \frac{m v}{b v^2 + \mu m g} dv = -\int_0^D dx$.
  7. [1 Point] Correct Final Algebraic Expression:
  8. Criterion: Correctly integrates natural logarithms and evaluates boundaries without sign errors.

Pitfall Alert: Never drop explicit integration limits or leave the constant of integration $+C$ un-evaluated. AP Readers immediately deduct the final point if $+C$ remains without initial condition substitution ($v(0) = v_0$).


4. Harvard University Placement Pathway

At Harvard University, achieving a 5 on AP Physics C: Mechanics along with AP Physics C: Electricity & Magnetism provides substantial placement leverage within the School of Engineering and Applied Sciences (SEAS) and the Department of Physics.

                           [AP Physics C: Mechanics (Score 5)]
                                            │
                                            ▼
                       [Exemption: SEAS Introductory Mechanics]
                                            │
               ┌────────────────────────────┴────────────────────────────┐
               ▼                                                         ▼
   [Physics 15b: E&M & Waves]                                 [Physics 16: Honors Mechanics]
  (Accelerated Physics Sequence)                             (Special Relativity Emphasis)
               │                                                         │
               └────────────────────────────┬────────────────────────────┘
                                            ▼
                             [Advanced Computational Physics]

Course Exemption & Acceleration Details

  1. Exempted Requirement:
  2. Fulfills the introductory mechanics requirement equivalent to Physics 15a (Introductory Mechanics and Special Relativity) or AP-level SEAS foundational prerequisites.
  3. Subsequent Placement:
  4. Physics 15b (Electromagnetism and Waves): High-achieving STEM concentrators bypass standard sequence bottlenecks and directly enter Physics 15b in their first year.
  5. Physics 16 (Honors Mechanics and Relativity): Students wanting deep theoretical exposure use their strong calculus-based mechanics foundation from AP Physics C to enroll directly in Harvard's famously rigorous honors track.

Bridge to Advanced Coursework

The dynamic mathematical techniques mastered in AP Physics C velocity-dependent drag problems—such as separable ordinary differential equations (ODEs), dynamic variable substitution, phase space trajectories, and non-conservative work line integrals—are directly applied in Physics 15b when calculating: * Velocity-dependent Lorentz Forces: $\vec{F} = q(\vec{E} + \vec{v} \times \vec{B})$ * Damped and driven harmonic oscillators: $m\frac{d^2 x}{dt^2} + b\frac{dx}{dt} + kx = F_0 \cos(\omega t)$ * Energy dissipation density in conductive media (Poynting Flux and Joule heating integrals)


5. High-Yield Practice Problem & Step-by-Step Solution Checklist

Problem Statement

A block of mass $m = 2.0 \text{ kg}$ is launched horizontally along a flat surface at initial speed $v_0 = 10.0 \text{ m/s}$ at position $x = 0$. The surface exerts both a constant kinetic friction force with coefficient $\mu = 0.20$ and a quadratic velocity-dependent air resistance force $F_d = -c v^2$, where $c = 0.10 \text{ kg/m}$. Take $g = 9.8 \text{ m/s}^2$.

                 Fd = -c*v^2      F_k = μ*m*g
                   <───          <───
               ┌──────────┐
  v0 = 10 m/s  │  Mass m  │ ───►  +x direction
  ────────────►│  2.0 kg  │
───────────────┴──────────┴───────────────────────────────
  x = 0                                           x = D

(a) Formulate a variable-separable differential equation for velocity as a function of position, $v(x)$.
(b) Solve the differential equation analytically to find an explicit expression for $v(x)$.
(c) Determine the exact distance $D$ traveled by the block before it comes to rest.
(d) Calculate the total work done by the quadratic drag force $W_{\text{drag}}$ during the interval from $x = 0$ to $x = D$. Prove that your result satisfies the Work-Energy Theorem.


Step-by-Step Solution Checklist

Part (a): Differential Equation Formulation

$$\sum F_x = -\mu m g - c v^2 = m a$$

$$m v \frac{dv}{dx} = -(\mu m g + c v^2)$$

$$\frac{m v}{\mu m g + c v^2} \, dv = -dx$$


Part (b): Analytical Integration for $v(x)$

$$\int_{v_0}^{v(x)} \frac{m v}{\mu m g + c v^2} \, dv = -\int_0^x dx'$$

$$\frac{m}{2c} \int_{u_0}^{u(x)} \frac{du}{u} = -x$$

$$\frac{m}{2c} \ln\left(\frac{\mu m g + c v(x)^2}{\mu m g + c v_0^2}\right) = -x$$

$$\frac{\mu m g + c v(x)^2}{\mu m g + c v_0^2} = e^{-\frac{2c}{m}x}$$

$$v(x)^2 = \left(v_0^2 + \frac{\mu m g}{c}\right) e^{-\frac{2c}{m}x} - \frac{\mu m g}{c}$$

$$v(x)^2 = (100 + 39.2) e^{-0.10 x} - 39.2 = 139.2 e^{-0.10 x} - 39.2$$

$$v(x) = \sqrt{139.2 e^{-0.10 x} - 39.2}$$


Part (c): Total Stopping Distance $D$

$$139.2 e^{-0.10 D} - 39.2 = 0$$

$$e^{-0.10 D} = \frac{39.2}{139.2} \approx 0.2816$$

$$-0.10 D = \ln(0.2816) \approx -1.2673 \implies D \approx 12.67 \text{ m}$$


Part (d): Non-Conservative Work Calculation & Theorem Verification

$$W_{\text{friction}} = -\mu m g D = -(0.20)(2.0)(9.8)(12.673) = -49.68 \text{ J}$$

$$W_{\text{net}} = W_{\text{friction}} + W_{\text{drag}} = \Delta K$$

$$-49.68 \text{ J} + W_{\text{drag}} = -100.0 \text{ J}$$

$$W_{\text{drag}} = -100.0 + 49.68 = -50.32 \text{ J}$$

$$W_{\text{drag}} = \int_0^D -c \left( 139.2 e^{-0.10 x} - 39.2 \right) dx$$

$$W_{\text{drag}} = -0.10 \left[ \frac{139.2}{-0.10} e^{-0.10 x} - 39.2 x \right]_0^{12.673}$$

$$W_{\text{drag}} = -0.10 \left[ \left(-1392 e^{-0.10(12.673)} - 39.2(12.673)\right) - \left(-1392 e^0 - 0\right) \right]$$

$$W_{\text{drag}} = -0.10 \left[ (-392 - 496.78) + 1392 \right] = -0.10 \left[ 503.22 \right] = -50.32 \text{ J}$$

Aiming for a Score 5 in Physics C: Mechanics?

Secure admission and advanced standing at top institutions like Harvard University with elite 1-on-1 AP STEM mentorship.

無料相談・学習プラン診断