AP Physics C: Mechanics Master Class
Work-Energy Theorem with Velocity-Dependent Drag Forces
1. Introduction & AP Exam Weight
The AP Physics C: Mechanics exam regularly tests top-tier analytical ability through non-conservative force problems involving dynamic variables. While basic Newtonian mechanics deals with constant forces, real-world systems—such as high-speed vehicles, projectiles, and fluid dynamics applications—involve velocity-dependent drag forces ($F_d \propto v$ or $F_d \propto v^2$).
Understanding how velocity-dependent drag forces interact with the Work-Energy Theorem represents the divide between a Score 4 and a Score 5. These problems account for approximately 10–15% of the AP Physics C: Mechanics exam, frequently anchoring the calculus-heavy Free Response Question (FRQ #1 or #2).
To achieve a 5, you cannot simply rely on $W = Fd \cos\theta$. You must apply differential calculus, execute spatial chain-rule transformations ($a = v \frac{dv}{dx}$), setup differential equations, and evaluate path integrals for non-conservative work:
$$W_{\text{net}} = \Delta K = \int_{x_i}^{x_f} F_{\text{net}}(x, v) \, dx$$
2. Deep Concept Breakdown
Theoretical Foundations & The Spatial Chain Rule
The generalized Work-Energy Theorem states that the net work done on a particle equals its change in kinetic energy:
$$W_{\text{net}} = W_{\text{conservative}} + W_{\text{non-conservative}} = \Delta K$$
When a resistive force depends on velocity (e.g., $F_d(v) = -bv$ for laminar flow or $F_d(v) = -cv^2$ for turbulent flow), time-based integration ($a = \frac{dv}{dt}$) yields $v(t)$, which then requires another integration to yield $x(t)$.
To evaluate work directly as a function of position without explicitly solving for time, you must convert the acceleration differential equation into position space using the Spatial Chain Rule:
$$a = \frac{dv}{dt} = \frac{dv}{dx} \frac{dx}{dt} = v \frac{dv}{dx}$$
Case Study: Quadratic Drag Force Analysis
Consider a particle of mass $m$ launched horizontally on a frictionless surface with initial velocity $v_0$, subject to a quadratic drag force $F_d(v) = -cv^2$.
+-------------------+
| Mass (m) | ---> v_0
+-------------------+
<--- F_d = -c*v^2
=================================================== (Frictionless Surface)
Step 1: Differential Equation Setup in Space
Applying Newton's Second Law along the axis of motion:
$$\Sigma F_x = m a_x \implies -c v^2 = m v \frac{dv}{dx}$$
Assuming $v \neq 0$, divide both sides by $v$:
$$-c v = m \frac{dv}{dx}$$
Step 2: Separation of Variables & Integration
Separate variables to integrate $v$ over position $x$ from $x = 0$ to $x$:
$$\int_{v_0}^{v(x)} \frac{1}{v} \, dv = -\frac{c}{m} \int_{0}^{x} \, dx'$$
$$\ln \left( \frac{v(x)}{v_0} \right) = -\frac{c}{m} x \implies v(x) = v_0 e^{-\frac{c}{m}x}$$
Step 3: Direct Work Integration vs. Change in Kinetic Energy
Method A: Direct Work Integration $$W_{\text{drag}} = \int_{0}^{x} F_d(x') \, dx' = \int_{0}^{x} \left( -c [v(x')]^2 \right) dx'$$
Substitute $v(x') = v_0 e^{-\frac{c}{m}x'}$:
$$W_{\text{drag}} = -c v_0^2 \int_{0}^{x} e^{-\frac{2c}{m}x'} \, dx' = -c v_0^2 \left[ -\frac{m}{2c} e^{-\frac{2c}{m}x'} \right]_{0}^{x}$$
$$W_{\text{drag}} = \frac{1}{2} m v_0^2 \left( e^{-\frac{2c}{m}x} - 1 \right)$$
Method B: Work-Energy Theorem $$\Delta K = \frac{1}{2} m [v(x)]^2 - \frac{1}{2} m v_0^2 = \frac{1}{2} m \left( v_0 e^{-\frac{c}{m}x} \right)^2 - \frac{1}{2} m v_0^2 = \frac{1}{2} m v_0^2 \left( e^{-\frac{2c}{m}x} - 1 \right)$$
Notice that $W_{\text{drag}} = \Delta K$, rigorously confirming the Work-Energy Theorem for variable non-conservative forces.
Python Verification Simulation
Below is a Python simulation using scipy.integrate.solve_ivp that numerically solves the spatial differential equation, computes the work done via numerical integration, and verifies $W_{\text{net}} = \Delta K$.
import numpy as np
from scipy.integrate import solve_ivp, cumulative_trapezoid
import matplotlib.pyplot as plt
# Physical Parameters
m = 2.0 # Mass in kg
c = 0.5 # Drag coefficient in kg/m
v0 = 20.0 # Initial velocity in m/s
x_max = 10.0 # Maximum position in meters
# 1. Differential Equation: dv/dx = - (c/m) * v
def dv_dx(x, v):
return - (c / m) * v
# Solve ODE in spatial domain
x_eval = np.linspace(0, x_max, 500)
sol = solve_ivp(dv_dx, [0, x_max], [v0], t_eval=x_eval)
x = sol.t
v = sol.y[0]
# 2. Kinetic Energy Calculation
K = 0.5 * m * v**2
delta_K = K - K[0]
# 3. Work Done by Drag Force Calculation: Integral(F_drag * dx)
F_drag = -c * v**2
W_drag = cumulative_trapezoid(F_drag, x, initial=0)
# Verify Work-Energy Equivalence numerically
max_error = np.max(np.abs(W_drag - delta_K))
print(f"Maximum absolute error between W_drag and Delta K: {max_error:.6e} Joules")
# Plotting Results
plt.figure(figsize=(10, 5))
plt.plot(x, delta_K, 'r--', label=r'$\Delta K(x)$ (Kinetic Energy Change)', linewidth=2.5)
plt.plot(x, W_drag, 'b-', label=r'$W_{drag}(x) = \int F_{drag} dx$', linewidth=1.5)
plt.title('Validation of Work-Energy Theorem under Velocity-Dependent Drag Force', 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.show()
3. Common AP Exam Pitfalls & Score 5 Scoring Rubric Nuances
Pitfall Analysis
- Misapplying Constant Force Equations: Writing $W = F_d \cdot x = -c v^2 x$. Because $v$ varies continuously with $x$, this earns zero points on AP FRQ rubrics. Force must be inside the integral: $\int F(x) \, dx$.
- Incorrect Differential Substitution: Substituting $a = \frac{dv}{dt}$ and integrating with respect to $x$ without applying the spatial chain rule ($a = v \frac{dv}{dx}$ or $dx = v dt$).
- Improper Integration Limits: Swapping initial and final limits or forgetting the negative sign on drag forces, yielding positive work for dissipative resistive forces.
- Confusing Terminal Velocity Contexts: Assuming $F_{\text{net}} = 0$ everywhere instead of strictly when $a = 0$.
Score 4 vs. Score 5 Solution Nuances
| Solution Step | Score 4 Candidate | Score 5 Candidate |
|---|---|---|
| Setup of Acceleration | Uses $a = \frac{dv}{dt}$, tries to solve $v(t)$, then plugs $v(t)$ directly into $\int F dx$ without converting variables. | Immediately applies spatial chain rule: $a = v \frac{dv}{dx}$, setting up a single spatial differential equation. |
| Separation of Variables | Omits differentials (e.g., writes $\int \frac{1}{v} = \int -\frac{c}{m}$) or leaves mixed variables ($v$ and $x$) on the same side. | Rigorously isolates variables: $\frac{1}{v} dv = -\frac{c}{m} dx$, preserving proper mathematical notation. |
| Definite Integration | Integrates indefinitely with $+ C$, often forgetting to evaluate $C$ using initial conditions $(x=0, v=v_0)$. | Uses definite integrals with matching boundaries: $\int_{v_0}^{v(x)} \frac{1}{v'} dv' = -\frac{c}{m} \int_{0}^{x} dx'$. |
| Work-Energy Proof | Assumes $W = \Delta K$ without proving equivalence when explicitly asked to derive work via calculus. | Proves equivalence via both explicit force integration ($\int F dx$) and kinetic energy delta ($\Delta K$). |
4. Stanford University Placement Pathway
At Stanford University, high academic performance on the AP Physics C: Mechanics exam unlocks key academic advancement opportunities:
AP Physics C: Mechanics (Score 5)
│
▼
Exempts PHYSICS 41 (4 Units)
(Mechanics Prerequisite)
│
▼
Immediate Winter Placement:
PHYSICS 43 (Electricity & Magnetism)
│
▼
Accelerated Majors & Electives:
├── AA 100 (Aeronautics & Astronautics)
├── ME 70 / AA 103 (Fluid Mechanics)
└── CME 102 (Ordinary Differential Equations)
- Exemption Granted: Direct credit/exemption for PHYSICS 41: Mechanics (4 quarter units).
- Accelerated Track: Bypassing PHYSICS 41 allows immediate enrollment in PHYSICS 43: Electricity & Magnetism during the Winter Quarter of your freshman year.
- Engineering Prerequisite Advantage: Fulfilling the mechanics core requirement early frees up 4 units toward major requirements in School of Engineering tracks, including:
- Aeronautics & Astronautics (AA)
- Mechanical Engineering (ME)
- Bioengineering (BioE)
- Academic Rigor Alignment: Stanford's fluid mechanics (ME 70 / AA 103) and dynamics (ENGR 14) courses rely heavily on differential formulations of work, power dissipation, and drag dynamics. Mastering variable force calculus gives you a distinct advantage in advanced coursework like CME 102 (Ordinary Differential Equations for Engineers).
5. High-Yield Practice Problem & Step-by-Step Solution
Problem Statement
A projectile of mass $m$ is fired vertically upward into the air with an initial speed $v_0$. As the projectile travels upward, it experiences a quadratic air resistance force of magnitude $F_d = c v^2$, where $c$ is a positive constant and $v$ is the speed of the projectile. Acceleration due to gravity $g$ is constant and downward.
^ +y
| +-------+
| | m | v (upward)
| +-------+
| |
| v F_d = c*v^2 AND F_g = m*g
===================================================
- Derive an expression for the velocity $v(y)$ of the projectile as a function of height $y$ above its launch point during its upward ascent.
- Derive an expression for the maximum height $y_{\text{max}}$ reached by the projectile.
- Calculate the work done by the drag force $W_{\text{drag}}$ during the ascent from $y = 0$ to $y = y_{\text{max}}$ using the Work-Energy Theorem.
Step-by-Step Solution Checklist & AP Marking Scheme
Part 1: Deriving $v(y)$
- Step 1: Free-Body Diagram & Newton's Second Law During upward motion, both gravity and drag point downward ($-\hat{j}$ direction):
$$\Sigma F_y = -m g - c v^2 = m a_y$$
- Step 2: Apply the Spatial Chain Rule Substitute $a_y = v \frac{dv}{dy}$:
$$-m g - c v^2 = m v \frac{dv}{dy}$$
- Step 3: Separate Variables Rearrange to isolate $v$ terms on the left and $y$ terms on the right:
$$\frac{m v}{m g + c v^2} \, dv = -dy$$
- Step 4: Integrate Both Sides Set limits from $(y=0, v=v_0)$ to $(y, v)$:
$$\int_{v_0}^{v(y)} \frac{m v}{m g + c v^2} \, dv = -\int_{0}^{y} \, dy'$$
Use $u$-substitution for the left integral: Let $u = m g + c v^2 \implies du = 2 c v \, dv \implies v \, dv = \frac{du}{2c}$.
$$\frac{m}{2c} \int_{u_0}^{u} \frac{1}{u} \, du = -y \implies \frac{m}{2c} \ln\left( \frac{m g + c [v(y)]^2}{m g + c v_0^2} \right) = -y$$
- Step 5: Solve for $v(y)$
$$\ln\left( \frac{m g + c [v(y)]^2}{m g + c v_0^2} \right) = -\frac{2cy}{m}$$
$$\frac{m g + c [v(y)]^2}{m g + c v_0^2} = e^{-\frac{2cy}{m}}$$
$$m g + c [v(y)]^2 = (m g + c v_0^2) e^{-\frac{2cy}{m}}$$
$$v(y) = \sqrt{ \frac{1}{c} \left[ (m g + c v_0^2) e^{-\frac{2cy}{m}} - m g \right] }$$
Part 2: Deriving Maximum Height $y_{\text{max}}$
-
At maximum height, instantaneous vertical velocity is zero: $v(y_{\text{max}}) = 0$.
-
Set the expression inside the square root to zero:
$$(m g + c v_0^2) e^{-\frac{2c y_{\text{max}}}{m}} - m g = 0$$
- Isolate the exponential term:
$$e^{-\frac{2c y_{\text{max}}}{m}} = \frac{m g}{m g + c v_0^2}$$
- Take the natural logarithm of both sides:
$$-\frac{2c y_{\text{max}}}{m} = \ln \left( \frac{m g}{m g + c v_0^2} \right) = -\ln \left( 1 + \frac{c v_0^2}{m g} \right)$$
- Solve for $y_{\text{max}}$:
$$y_{\text{max}} = \frac{m}{2c} \ln \left( 1 + \frac{c v_0^2}{m g} \right)$$
Part 3: Computing $W_{\text{drag}}$ via Work-Energy Theorem
- Apply the Work-Energy Theorem between $y = 0$ and $y = y_{\text{max}}$:
$$W_{\text{net}} = W_g + W_{\text{drag}} = \Delta K$$
- Calculate individual energy terms:
$$\Delta K = K_f - K_i = 0 - \frac{1}{2} m v_0^2 = -\frac{1}{2} m v_0^2$$
$$W_g = \int_{0}^{y_{\text{max}}} (-m g) \, dy = -m g y_{\text{max}}$$
- Substitute $W_g$ and $\Delta K$ into the Work-Energy Theorem:
$$-m g y_{\text{max}} + W_{\text{drag}} = -\frac{1}{2} m v_0^2$$
$$W_{\text{drag}} = m g y_{\text{max}} - \frac{1}{2} m v_0^2$$
- Substitute the derived expression for $y_{\text{max}}$:
$$W_{\text{drag}} = m g \left[ \frac{m}{2c} \ln \left( 1 + \frac{c v_0^2}{m g} \right) \right] - \frac{1}{2} m v_0^2$$
$$W_{\text{drag}} = \frac{m^2 g}{2c} \ln \left( 1 + \frac{c v_0^2}{m g} \right) - \frac{1}{2} m v_0^2$$
AP FRQ Scoring Rubric Checklist
| Criteria | Point Allocation | Status |
|---|---|---|
| Newton's 2nd Law & Chain Rule | +1 Pt: Correctly states $-mg - cv^2 = m v \frac{dv}{dy}$. | $\checkmark$ |
| Separation of Variables | +1 Pt: Correctly isolates variables $\frac{m v}{mg + cv^2} dv = -dy$. | $\checkmark$ |
| Definite Integration | +1 Pt: Integrates using appropriate limits or solves for integration constant $C$. | $\checkmark$ |
| Correct $v(y)$ Expression | +1 Pt: Explicitly isolates $v(y)$ in a mathematically equivalent form. | $\checkmark$ |
| Setting $v=0$ for $y_{\text{max}}$ | +1 Pt: Recognizes condition for peak height and isolates $y_{\text{max}}$. | $\checkmark$ |
| Work-Energy Application | +1 Pt: Sets up $W_{\text{drag}} = \Delta K - W_g$ using $y_{\text{max}}$. | $\checkmark$ |
| Total Available Points | 6 Points (Full FRQ Section) | Score 5 |