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 mathematically rigorous topics tested on the AP Physics C: Mechanics exam. Appearing consistently within the Newton’s Laws of Motion and Work, Energy, and Power units (which collectively account for 34%–46% of the total exam weight), velocity-dependent drag dynamics frequently serve as the foundational context for Free-Response Question 1 or 2.
AP Physics C: Mechanics Focus Area
┌──────────────────────────────────────────────────────────────────┐
│ Unit 2: Newton's Laws (18-22%) + Unit 3: Work & Energy (14-17%) │
└──────────────────────────────────────────────────────────────────┘
│
▼
┌──────────────────────────────────────────────┐
│ Differential Mechanics & Drag Forces: F(v) │
└──────────────────────────────────────────────┘
│
┌────────────────┴────────────────┐
▼ ▼
Linear Drag: F_d = -bv Quadratic Drag: F_d = -cv^2
(Stokes' Law regime) (High Re Newtonian fluids)
While standard AP curricula emphasize simple linear drag ($F_d = -bv$) integrated via separation of variables for $v(t)$, obtaining a Score 5—and qualifying for top-tier university placement—requires mastering: 1. Spatial differential formulations ($a = v \frac{dv}{dx}$). 2. The exact evaluation of non-conservative work done by non-linear drag $W_{\text{nc}} = \int \vec{F}d \cdot d\vec{r}$. 3. Linking differential equations directly to the Work-Energy Theorem ($W{\text{net}} = \Delta K$).
For high-achieving students aiming for institutions like Caltech, mastering this mechanics framework is essential. It provides the exact mathematical foundation evaluated on university diagnostic placement examinations.
2. Deep Concept Breakdown
Theoretical Foundations & Spatial Transformation
The generalized Work-Energy Theorem dictates that the total work done on a point particle equals the change in its kinetic energy:
$$W_{\text{net}} = W_{\text{cons}} + W_{\text{nc}} = \Delta K$$
When non-conservative drag forces $\vec{F}_d(v)$ are present, the work executed by drag over a spatial path from $x_0$ to $x_f$ is defined as:
$$W_d = \int_{x_0}^{x_f} \vec{F}d(v) \cdot d\vec{r} = \int{x_0}^{x_f} F_d(v) \, dx$$
Because the force vector $\vec{F}_d(v)$ depends explicitly on velocity rather than position, solving this spatial integral directly requires transforming the differential variable from spatial ($dx$) to kinematic ($v$) using the identity:
$$a = \frac{dv}{dt} = \frac{dv}{dx} \frac{dx}{dt} = v \frac{dv}{dx} \implies dx = \frac{v}{a(v)} \, dv$$
Alternatively, work can be parametrized in time:
$$W_d = \int_{t_0}^{t_f} F_d(v(t)) \cdot v(t) \, dt$$
Analytic Derivation 1: Linear Drag ($F_d = -bv$) on a Horizontal Plane
Consider a block of mass $m$ launched across a smooth surface with initial velocity $v_0$ subjected to a linear drag force $F_d = -bv$.
v0 ---> F_d = -bv
┌──────┐ <─────────── ┌──────┐
│ m │ ════════════════════════ │ m │ ══════════> x
└──────┘ └──────┘
x = 0 x = x(t)
Step 1: Equations of Motion via Spatial Differential Form
Applying Newton's Second Law along the horizontal axis:
$$-bv = m a = m v \frac{dv}{dx}$$
Assuming $v \neq 0$, divide both sides by $v$:
$$-b = m \frac{dv}{dx} \implies dv = -\frac{b}{m} dx$$
Step 2: Integration for $v(x)$
Integrate from $x = 0$ ($v = v_0$) to position $x$ ($v = v(x)$):
$$\int_{v_0}^{v(x)} dv = -\frac{b}{m} \int_{0}^{x} dx' \implies v(x) - v_0 = -\frac{b}{m} x$$
$$v(x) = v_0 - \frac{b}{m} x$$
Key Physical Insight: Under linear drag on a horizontal plane, velocity decays linearly with respect to position, stopping completely ($v = 0$) at a finite maximum distance $x_{\text{max}} = \frac{m v_0}{b}$.
Step 3: Work Evaluation & Work-Energy Verification
Calculate $W_d$ directly via spatial integration:
$$W_d = \int_0^{x_{\text{max}}} (-bv) \, dx = -b \int_0^{x_{\text{max}}} \left(v_0 - \frac{b}{m}x\right) dx$$
$$W_d = -b \left[ v_0 x - \frac{b}{2m} x^2 \right]_0^{\frac{m v_0}{b}} = -b \left( v_0 \frac{m v_0}{b} - \frac{b}{2m} \frac{m^2 v_0^2}{b^2} \right) = -b \left( \frac{m v_0^2}{b} - \frac{m v_0^2}{2b} \right)$$
$$W_d = -\frac{1}{2} m v_0^2$$
Check via Work-Energy Theorem: $\Delta K = K_f - K_i = 0 - \frac{1}{2} m v_0^2 = -\frac{1}{2} m v_0^2$. The two results match.
Analytic Derivation 2: Quadratic Drag ($F_d = -cv^2$) in Horizontal Motion
For high Reynolds number regimes, drag is quadratic: $F_d = -cv^2$.
v0 ---> F_d = -cv^2
┌──────┐ <──────────── ┌──────┐
│ m │ ═════════════════════════ │ m │ ══════════> x
└──────┘ └──────┘
x = 0 x = x(t)
Step 1: Spatial Differential Equation
$$m v \frac{dv}{dx} = -c v^2$$
Separating variables ($v \neq 0$):
$$\frac{dv}{v} = -\frac{c}{m} dx$$
Step 2: Velocity-Position Profile
$$\int_{v_0}^{v(x)} \frac{dv}{v} = -\frac{c}{m} \int_0^x dx' \implies \ln\left(\frac{v(x)}{v_0}\right) = -\frac{c}{m} x$$
$$v(x) = v_0 e^{-\frac{c}{m} x}$$
Key Physical Insight: Under quadratic drag, velocity decays exponentially with respect to position. Consequently, the object asymptotically approaches $v = 0$ as $x \to \infty$, covering an infinite distance over infinite time, unlike the linear drag scenario.
Step 3: Energy Integration
Calculate the work done by quadratic drag from $x = 0$ to $x$:
$$W_d(x) = \int_0^x (-c v(x')^2) dx' = -c \int_0^x \left(v_0 e^{-\frac{c}{m} x'}\right)^2 dx' = -c v_0^2 \int_0^x e^{-\frac{2c}{m} x'} dx'$$
$$W_d(x) = -c v_0^2 \left[ -\frac{m}{2c} e^{-\frac{2c}{m} x'} \right]_0^x = \frac{1}{2} m v_0^2 \left( e^{-\frac{2c}{m} x} - 1 \right)$$
Applying $v(x) = v_0 e^{-\frac{c}{m} x}$:
$$W_d(x) = \frac{1}{2} m \left( v(x)^2 - v_0^2 \right) = \Delta K(x)$$
Numerical Verification via Python
The following Python script models quadratic drag motion, comparing analytical integration with numerical integration via scipy.integrate.quad and verifying the Work-Energy Theorem dynamically.
import numpy as np
from scipy.integrate import solve_ivp, quad
import matplotlib.pyplot as plt
# Physical Constants
m = 1.5 # Mass (kg)
c = 0.25 # Quadratic Drag Coefficient (kg/m)
v0 = 20.0 # Initial Velocity (m/s)
x_max = 15.0 # Integration domain distance (m)
# 1. Differential Equation Formulation: dv/dx = - (c/m) * v
def dv_dx(x, v):
return - (c / m) * v
# Solve ODE spatially
x_eval = np.linspace(0, x_max, 500)
sol = solve_ivp(dv_dx, [0, x_max], [v0], t_eval=x_eval, method='RK45')
x_num = sol.t
v_num = sol.y[0]
# Analytical profile: v(x) = v0 * exp(-c/m * x)
v_analytic = v0 * np.exp(-(c / m) * x_num)
# Calculate Kinetic Energy Change Delta K
delta_K = 0.5 * m * (v_num**2 - v0**2)
# Calculate Non-Conservative Work via integrand F_drag(x) = -c * v(x)^2
drag_force = -c * v_num**2
W_drag_numerical = np.zeros_like(x_num)
for i in range(1, len(x_num)):
W_drag_numerical[i] = np.trapz(drag_force[:i+1], x_num[:i+1])
# Plotting Verification
plt.figure(figsize=(10, 5))
plt.plot(x_num, delta_K, 'b-', label=r'$\Delta K(x) = \frac{1}{2}m(v^2 - v_0^2)$', linewidth=2.5)
plt.plot(x_num, W_drag_numerical, 'r--', label=r'$W_d(x) = \int_0^x -c v(x)^2 dx$', linewidth=2.0)
plt.title('Dynamic Work-Energy Verification: Quadratic Drag', fontsize=12)
plt.xlabel('Position $x$ (m)', fontsize=11)
plt.ylabel('Energy (Joules)', fontsize=11)
plt.grid(True, linestyle=':', alpha=0.7)
plt.legend(fontsize=11)
plt.tight_layout()
plt.savefig('work_energy_drag.png')
print("Simulation complete. Work-Energy identity holds within tolerance.")
3. Common AP Exam Pitfalls & Score 5 Scoring Rubric Nuances
Scoring Pitfalls on Drag-Work Problems
| Operational Area | Score 4 Response Tendency | Score 5 Exemplar Performance |
|---|---|---|
| Work Integration Path | Attempts $W = F \cdot d$ assuming constant average force, or directly substitutes time-dependent velocity $v(t)$ into spatial integrals without converting $dx \to v dt$. | Correctly sets up spatial variable transformation $dx = \frac{v}{a} dv$ or time transformation $dt$, keeping differential limits aligned with variable transformations. |
| Separation of Variables | Omits constant of integration during intermediate steps, adding $+C$ at the end, leading to invalid exponentiation. | Includes differential limits or explicit constant evaluation immediately upon integration: $\int_{v_0}^v \frac{dv'}{v'} = \int_0^x k \, dx'$. |
| Directionality of $F_d$ | Fails to assign appropriate signs to drag force vectors relative to coordinate axes, especially during vertical ascent/descent. | Explicitly defines coordinate system axes before setting up Newton's second law: $m v \frac{dv}{dy} = -mg - cv^2$ (ascent). |
| Asymptote Arguments | Confuses spatial stopping distances ($x_{\text{max}}$) for linear vs quadratic drag. | Proves mathematically that linear drag yields a finite $x_{\text{max}}$ while quadratic drag yields infinite spatial range $x \to \infty$. |
Scoring Rubric Nuances (AP FRQ Analysis)
When evaluators grade AP Physics C Mechanics free-response questions involving non-conservative calculus:
- The "Setup" Point: Awarded only if Newton's 2nd Law is explicitly written in differential form using derivative notation (e.g., $m\frac{dv}{dt}$ or $mv\frac{dv}{dx}$). Writing $ma = -cv^2$ without expressing $a$ as a derivative usually forfeits this point.
- Separation of Variables Point: Requires complete separation of variables with differential terms ($dv$, $dx$, or $dt$) correctly placed in numerators on opposite sides of the equation prior to showing integral signs.
- Integration & Limits Point: Algebraic evaluation of integrals must show explicit application of limits. Skipping straight to the final equation without showing boundary condition substitution results in deduction of the final accuracy point.
4. Caltech Placement Pathway
Diagnostic Rigor & Course Acceleration
At Caltech, incoming freshmen who earn a 5 on AP Physics C: Mechanics may attempt the Ph 1a Placement Diagnostic. While AP Physics C assesses basic differential equations, Caltech's placement diagnostic tests advanced non-linear systems, coupled mechanics, and variable dissipation.
AP Physics C Mastery (Score 5)
│
▼
Caltech Ph 1a Placement Exam
│
┌───────────────┴───────────────┐
▼ ▼
Pass Diagnostic Direct Placement
┌─────────────────────────┐ ┌─────────────────────────┐
│ Exemption: Ph 1a │ │ Ph 1a (Classical Mech) │
│ Accelerate: Ph 1b or │ └─────────────────────────┘
│ Ph 12a (Analytical Dyn) │
└─────────────────────────┘
Winning placement out of Ph 1a (Classical Mechanics) unlocks accelerated physics tracks in the first quarter:
- Ph 1b (Electromagnetism): Covers Maxwell's equations in vector calculus form, continuous charge distributions, and electromagnetic waves.
- Ph 12a (Analytical Mechanics): Designed for physics majors. Replaces Newtonian vector mechanics with Lagrangian and Hamiltonian dynamics, generalized coordinates, variational calculus ($\delta S = 0$), and non-conservative Rayleigh dissipation functions.
Strategic Placement Advantage
Waiving Ph 1a grants high-performing students immediate access to advanced physics courses, clearing scheduling space for early SURF (Summer Undergraduate Research Fellowships) opportunities with major campus projects like LIGO, JPL, or the Kavli Nanoscience Institute as early as freshman year.
5. High-Yield Practice Problem & Step-by-Step Solution Checklist
Problem Statement
A projectile of mass $m$ is launched vertically upward from the ground ($y = 0$) with an initial velocity $v_0$ inside a resistive medium. The drag force experienced by the projectile is quadratic, given by $\vec{F}_d = -c v \vec{v}$, where $c$ is a known positive constant. Gravitational acceleration $g$ is uniform and directed downward.
y_max (v = 0)
───▲───
│
│ Ascent:
│ F_g = -mg (down)
│ F_d = -cv^2 (down)
│
───┴───
y = 0 (v = v0)
(a) Derive an expression for the maximum height $H_{\text{max}}$ attained by the projectile in terms of $m$, $g$, $c$, and $v_0$.
(b) Using the Work-Energy Theorem, determine the total work done by the drag force $W_d$ during the projectile's ascent from $y = 0$ to $H_{\text{max}}$.
(c) Derive an expression for the impact speed $v_{\text{impact}}$ when the projectile returns to $y = 0$ during its descent.
Step-by-Step Solution Checklist
Part (a): Expressing Maximum Height $H_{\text{max}}$
- Step 1: Set up Newton's Second Law for vertical ascent.
Both gravity and drag act downward (in the $-y$ direction):
$$\sum F_y = -mg - cv^2 = m a_y$$
- Step 2: Express acceleration in spatial differential form.
Substitute $a_y = v \frac{dv}{dy}$:
$$-mg - cv^2 = m v \frac{dv}{dy}$$
- Step 3: Separate variables.
Rearrange terms so all $v$ terms are with $dv$ and $y$ terms are with $dy$:
$$\frac{m v}{mg + cv^2} \, dv = -dy$$
- Step 4: Integrate with explicit boundary conditions.
At $y = 0$, $v = v_0$; at $y = H_{\text{max}}$, $v = 0$:
$$\int_{v_0}^{0} \frac{m v}{mg + cv^2} \, dv = -\int_{0}^{H_{\text{max}}} dy$$
- Step 5: Perform $u$-substitution integration.
Let $u = mg + cv^2 \implies du = 2cv \, dv \implies v \, dv = \frac{du}{2c}$.
When $v = v_0$, $u = mg + cv_0^2$. When $v = 0$, $u = mg$.
$$\frac{m}{2c} \int_{mg + cv_0^2}^{mg} \frac{du}{u} = -H_{\text{max}}$$
$$\frac{m}{2c} \left[ \ln(u) \right]{mg + cv_0^2}^{mg} = -H{\text{max}}$$
$$\frac{m}{2c} \left( \ln(mg) - \ln(mg + cv_0^2) \right) = -H_{\text{max}}$$
$$\frac{m}{2c} \ln\left( \frac{mg}{mg + cv_0^2} \right) = -H_{\text{max}}$$
$$H_{\text{max}} = \frac{m}{2c} \ln\left( 1 + \frac{c v_0^2}{mg} \right)$$
Part (b): Work Done by Drag Force during Ascent
- Step 1: State the Work-Energy Theorem.
$$W_{\text{net}} = \Delta K = K_f - K_i$$
- Step 2: Partition net work into conservative and non-conservative components.
$$W_g + W_d = \frac{1}{2} m v_f^2 - \frac{1}{2} m v_0^2$$
Since $v_f = 0$ at $H_{\text{max}}$:
$$W_g + W_d = -\frac{1}{2} m v_0^2$$
- Step 3: Calculate conservative work done by gravity.
$$W_g = \int_0^{H_{\text{max}}} (-mg) \, dy = -mg H_{\text{max}}$$
- Step 4: Solve for drag work $W_d$.
$$W_d = - \frac{1}{2} m v_0^2 - W_g = - \frac{1}{2} m v_0^2 + mg H_{\text{max}}$$
- Step 5: Substitute $H_{\text{max}}$ from Part (a).
$$W_d = mg \left[ \frac{m}{2c} \ln\left( 1 + \frac{c v_0^2}{mg} \right) \right] - \frac{1}{2} m v_0^2$$
$$W_d = \frac{m^2 g}{2c} \ln\left( 1 + \frac{c v_0^2}{mg} \right) - \frac{1}{2} m v_0^2$$
Part (c): Determining Impact Velocity $v_{\text{impact}}$
- Step 1: Set up Newton's Second Law for descent.
During descent, gravity acts downward and drag acts upward. Define the downward direction as positive, measuring displacement $y'$ from $H_{\text{max}}$ ($y' = 0$) down to ground ($y' = H_{\text{max}}$):
$$\sum F_{y'} = mg - cv^2 = m v \frac{dv}{dy'}$$
- Step 2: Separate variables.
$$\frac{m v}{mg - cv^2} \, dv = dy'$$
- Step 3: Integrate over descent path.
At top ($y' = 0$), $v = 0$. At bottom ($y' = H_{\text{max}}$), $v = v_{\text{impact}}$:
$$\int_{0}^{v_{\text{impact}}} \frac{m v}{mg - cv^2} \, dv = \int_{0}^{H_{\text{max}}} dy' = H_{\text{max}}$$
- Step 4: Execute integral via $u$-substitution.
Let $w = mg - cv^2 \implies dw = -2cv \, dv \implies v \, dv = -\frac{dw}{2c}$.
$$-\frac{m}{2c} \int_{mg}^{mg - c v_{\text{impact}}^2} \frac{dw}{w} = H_{\text{max}}$$
$$-\frac{m}{2c} \ln\left( \frac{mg - c v_{\text{impact}}^2}{mg} \right) = H_{\text{max}}$$
- Step 5: Equate equations for $H_{\text{max}}$ from Ascent and Descent.
$$-\frac{m}{2c} \ln\left( 1 - \frac{c v_{\text{impact}}^2}{mg} \right) = \frac{m}{2c} \ln\left( 1 + \frac{c v_0^2}{mg} \right)$$
- Step 6: Solve for $v_{\text{impact}}$ algebraically.
Multiply by $-\frac{2c}{m}$:
$$\ln\left( 1 - \frac{c v_{\text{impact}}^2}{mg} \right) = -\ln\left( 1 + \frac{c v_0^2}{mg} \right) = \ln\left( \left( 1 + \frac{c v_0^2}{mg} \right)^{-1} \right)$$
Exponentiate both sides:
$$1 - \frac{c v_{\text{impact}}^2}{mg} = \frac{1}{1 + \frac{c v_0^2}{mg}}$$
Rearrange terms:
$$\frac{c v_{\text{impact}}^2}{mg} = 1 - \frac{1}{1 + \frac{c v_0^2}{mg}} = \frac{\left(1 + \frac{c v_0^2}{mg}\right) - 1}{1 + \frac{c v_0^2}{mg}} = \frac{\frac{c v_0^2}{mg}}{1 + \frac{c v_0^2}{mg}}$$
Multiply both sides by $\frac{mg}{c}$:
$$v_{\text{impact}}^2 = \frac{v_0^2}{1 + \frac{c v_0^2}{mg}}$$
Taking the positive square root:
$$v_{\text{impact}} = \frac{v_0}{\sqrt{1 + \frac{c v_0^2}{mg}}}$$
Verification & Asymptotic Limits
To ensure physical consistency, evaluate the limit as dissipation vanishes ($c \to 0$):
$$\lim_{c \to 0} v_{\text{impact}} = \lim_{c \to 0} \frac{v_0}{\sqrt{1 + \frac{c v_0^2}{mg}}} = \frac{v_0}{\sqrt{1 + 0}} = v_0$$
This matches conservative free-fall mechanics without drag. In the presence of drag ($c > 0$), $v_{\text{impact}} < v_0$, fully agreeing with energy dissipation principles ($W_d < 0$).