4.10b - Modelling with differential equations

4.10b - Modelling with differential equations

A differential equation turns a statement about rates or forces into a prediction. You will construct and interpret models in kinematics and other contexts, solve suitable first order models, and distinguish a model you can formulate from one you can solve analytically.

Variable force and motion

Can a constant driving force produce a speed limit? It can when resistance grows with velocity. To understand that result, first translate a simpler changing force into an equation for motion.

You already know

For a particle of constant mass, Newton's second law is resultant force equals mass times acceleration. Choose a positive direction before assigning force and velocity signs. Integrating acceleration gives velocity; integrating velocity gives displacement.

Let xx be the signed displacement in the chosen direction, measured at time tt. A dot means differentiation with respect to time; it does not mean multiplication.

Kinematics and Newton's second law

For constant mass m>0m>0, in a fixed direction,

v=x˙=dxdt,a=v˙=x¨=dvdt=d2xdt2v=\dot{x}=\frac{dx}{dt},\qquad a=\dot{v}=\ddot{x}=\frac{dv}{dt}=\frac{d^2x}{dt^2}

and

Fresultant=mv˙=mx¨F_{\mathrm{resultant}}=m\dot{v}=m\ddot{x}

The force and velocity are signed components along that direction.
Not in the formula booklet: learn it.

When the force changes, the constant-acceleration formulae do not apply. Construct the derivative equation and integrate it instead.

Worked example 1: a force proportional to time

Question. A particle of mass 2 kg2\,\mathrm{kg} moves along a straight line. Its resultant force at time tt seconds is 6t N6t\,\mathrm{N} in the positive direction. Initially v=1 m s−1v=1\,\mathrm{m\,s^{-1}} and x=0x=0. Find vv and xx in terms of tt.

Plan. Use Newton's second law, integrate for velocity and displacement, and use each initial condition.

2dvdt=6tNewton’s second lawv=32t2+Cintegrate 3tC=1v(0)=1x=t+12t3+Dintegrate velocityD=0x(0)=0\begin{aligned} 2\frac{dv}{dt}&=6t &&\text{Newton's second law}\\ v&=\frac32t^2+C &&\text{integrate }3t\\ C&=1 && v(0)=1\\ x&=t+\frac12t^3+D &&\text{integrate velocity}\\ D&=0 && x(0)=0 \end{aligned}

Why this step? The given force is the resultant, so no extra force is added to it; divide by the mass to obtain acceleration.

Answer. v=1+32t2 m s−1v=1+\dfrac32t^2\,\mathrm{m\,s^{-1}} and x=t+12t3 mx=t+\dfrac12t^3\,\mathrm{m}, for t≥0t\ge0.

Check. Differentiating gives x˙=1+3t2/2=v\dot{x}=1+3t^2/2=v and 2v˙=6t2\dot{v}=6t, and both initial conditions hold.

Add only a constant opposing contribution of 4 N4\,\mathrm{N} to the resultant force. Its sign must survive both integrations.

Worked example 2: a signed changing resultant

Question. The same particle now has resultant force (6t−4) N(6t-4)\,\mathrm{N}, measured positively in the same direction, with the same initial conditions. Find vv and xx.

Plan. Retain the negative force term when forming acceleration, then integrate twice.

2v˙=6t−4signed resultantv=32t2−2t+Cintegrate accelerationC=1v(0)=1x=12t3−t2+t+Dintegrate velocityD=0x(0)=0\begin{aligned} 2\dot{v}&=6t-4 &&\text{signed resultant}\\ v&=\frac32t^2-2t+C &&\text{integrate acceleration}\\ C&=1 && v(0)=1\\ x&=\frac12t^3-t^2+t+D &&\text{integrate velocity}\\ D&=0 && x(0)=0 \end{aligned}

Why this step? A negative resultant means negative acceleration, not automatically negative velocity; the initial motion still matters.

Answer. v=1−2t+32t2 m s−1v=1-2t+\dfrac32t^2\,\mathrm{m\,s^{-1}} and x=t−t2+12t3 mx=t-t^2+\dfrac12t^3\,\mathrm{m}, for t≥0t\ge0.

Check. At t=0t=0, acceleration is −2 m s−2-2\,\mathrm{m\,s^{-2}} while velocity is +1 m s−1+1\,\mathrm{m\,s^{-1}}, consistent with initially slowing down.

Signs encode the direction of the force; taking its magnitude too early changes the model.

Watch out

Wrong: 2v˙=∣6t−4∣2\dot{v}=|6t-4|, replacing the signed resultant by its size. Right: 2v˙=6t−42\dot{v}=6t-4 throughout the model's stated interval.

Change the constant opposing term once more, and keep velocity signed.

Velocity dependent resistance

A force depending on velocity makes the unknown function appear on both sides of Newton's law. With resistance proportional to velocity, the resulting equation is linear and the integrating factor method applies.

Suppose a constant drive F>0F>0 acts positively and a resistive force has signed component −kv-kv, with k>0k>0. The minus sign makes resistance oppose the motion. Newton's law is mv˙=F−kvm\dot{v}=F-kv. At a steady velocity, v˙=0\dot{v}=0, so the drive and resistance balance.

Terminal velocity

A terminal velocity is a limiting steady velocity approached as the resultant force tends to zero. In the model mv˙=F−kvm\dot{v}=F-kv, it is F/kF/k. For vertical falling motion with downward positive, the corresponding model is mv˙=mg−kvm\dot{v}=mg-kv, giving terminal velocity mg/kmg/k.

Terminal velocity is a prediction of the force law. It is not an extra force to add to the equation.

These terminal-velocity results are not in the formula booklet; obtain them by setting the resultant force to zero.

[DIAGRAM: asset_name: Drive and resistance approach a balance; asset_slug: h245_4_10_terminal_velocity; caption: For a 3 kg particle with 12 N drive and resistance 3v N, velocity rises towards 4 metres per second while acceleration tends to zero; description: Left panel shows a particle with a rightward arrow labelled drive 12 N, a leftward arrow labelled resistance 3v N, and a rightward positive direction. Right panel plots v = 4(1 - exp(-t)) for t from 0 to 5, starting at the origin and approaching a dashed horizontal terminal line v = 4. The half-terminal point (ln 2, 2) is marked and labelled.]
For a 3 kg particle with 12 N drive and resistance 3v N, velocity rises towards 4 metres per second while acceleration tends to zero

Worked example 3: approaching terminal velocity

Question. A 3 kg3\,\mathrm{kg} particle starts from rest on a horizontal line. A constant drive of 12 N12\,\mathrm{N} acts positively and the resistance has signed component −3v N-3v\,\mathrm{N}, where vv is in m s−1\mathrm{m\,s^{-1}}. Find v(t)v(t) and its limiting value.

Plan. Resolve forces, divide by the mass and solve the linear equation with its initial condition.

3v˙=12−3vresolve positivelyv˙+v=4normaliseddt(etv)=4etfactor etetv=4et+Cintegrate0=4+Cv(0)=0v=4(1−e−t)divide by et\begin{aligned} 3\dot{v}&=12-3v &&\text{resolve positively}\\ \dot{v}+v&=4 &&\text{normalise}\\ \frac{d}{dt}(e^tv)&=4e^t &&\text{factor }e^t\\ e^tv&=4e^t+C &&\text{integrate}\\ 0&=4+C && v(0)=0\\ v&=4(1-e^{-t}) &&\text{divide by }e^t \end{aligned}

Why this step? Resistance is part of the resultant, so the equation is 3v˙=12−3v3\dot{v}=12-3v, rather than 3v˙=123\dot{v}=12.

Answer. v=4(1−e−t) m s−1v=4(1-e^{-t})\,\mathrm{m\,s^{-1}}, tending to 4 m s−14\,\mathrm{m\,s^{-1}} as t→∞t\to\infty.

Check. 3v˙=12e−t=12−3[4(1−e−t)]3\dot{v}=12e^{-t}=12-3[4(1-e^{-t})]. At the limiting velocity, resistance is 3(4)=12 N3(4)=12\,\mathrm{N}, balancing the drive.

Change only the driving force from 12 N12\,\mathrm{N} to 18 N18\,\mathrm{N}. The rate of approach stays the same, but the limiting velocity changes.

Worked example 4: changing the drive

Question. The same particle, starting from rest, has drive 18 N18\,\mathrm{N} and the same resistance law. Find its velocity and the time to reach half its terminal velocity.

Plan. Solve the new linear equation, then set the velocity equal to half its limiting value.

v˙+v=6Newton’s law divided by 3etv=6et+Cfactor et and integrateC=−6v(0)=0v=6(1−e−t)divide by et3=6(1−e−t)⇒e−t=12half the limitt=ln⁡2take logarithms\begin{aligned} \dot{v}+v&=6 &&\text{Newton's law divided by }3\\ e^tv&=6e^t+C &&\text{factor }e^t\text{ and integrate}\\ C&=-6 && v(0)=0\\ v&=6(1-e^{-t}) &&\text{divide by }e^t\\ 3=6(1-e^{-t})&\Rightarrow e^{-t}=\frac12 &&\text{half the limit}\\ t&=\ln2 &&\text{take logarithms} \end{aligned}

Why this step? Half-terminal velocity is half the new limit 66, rather than half the previous limit 44.

Answer. v=6(1−e−t) m s−1v=6(1-e^{-t})\,\mathrm{m\,s^{-1}}; the required time is ln⁡2 s\ln2\,\mathrm{s}.

Check. At t=ln⁡2t=\ln2, e−t=1/2e^{-t}=1/2, so v=3 m s−1v=3\,\mathrm{m\,s^{-1}}.

In fact, the first particle also reaches half its terminal velocity at ln⁡2 s\ln2\,\mathrm{s}. Increasing the drive raises both velocities in the same proportion; it does not change the time needed to complete the same fraction of the approach in this linear model. This conclusion depends on the assumed resistance law.

Continue this lesson.

You’re previewing NovaLearn without an account. Create a free account to read the rest.

Create free account