Lecture The Energy principle
1.10 The Energy principle
1.10.1 Conservation of energy
In many cases, the ODE describing a physical system can be simplified.
Consider an object with mass $m$ subject to a force, $F(t).$ From Newton's Second Law of Motion, the object's trajectory is described as
$ \ds m\frac{d^2 y}{dt^2}=F(y,t). $
Often, the force $F$ depends only on position $y(t).$
1.10.1 Conservation of energy
Often, the force $F$ depends only on position $y(t).$
We write
$ \ds m\frac{d^2 y}{dt^2}=F(y(t)).~~ (\star) $
By the Fundamental Theorem of Calculus, if $F$ is a continuous function, it has an anti-derivative $\phi(y)$ . Hence, $\displaystyle \frac{d\phi}{dy}=F(y).$
1.10.1 Conservation of energy
We multiply both sides by $\displaystyle \frac{dy}{dt}$
$ \ds m\frac{d^2 y}{dt^2}\frac{dy}{dt} =\frac{d\phi}{dy}\frac{dy}{dt}. $
Applying the chain rule, the LHS becomes
$\ds m\frac{dy}{dt}\frac{d^2 y}{dt^2}$ $\ds =\frac{d}{dt}\left[\frac{1}{2}m\left(\frac{dy}{dt}\right)^2\right]$
1.10.1 Conservation of energy
while the RHS is
$\ds\frac{d\phi}{dy}\frac{dy}{dt}$ $\ds =\frac{d\phi(y(t))}{dt}. $
Equation $m\dfrac{d^2 y}{dt^2}=F(y,t)$ may then be written as
$\ds\frac{d}{dt}\left[\frac{1}{2}m\left(\frac{dy}{dt}\right)^2\right] = \frac{d}{dt}\phi(y(t))$
1.10.1 Conservation of energy
Equation $m\dfrac{d^2 y}{dt^2}=F(y,t)$ may then be written as
$\ds\frac{d}{dt}\left[\frac{1}{2}m\left(\frac{dy}{dt}\right)^2\right] = \frac{d}{dt}\phi(y(t))$
$\ds\Ra \frac{d}{dt}\left(\frac{1}{2}m\left(\frac{dy}{dt}\right)^2-\phi(y(t)) \right)= 0$
$\ds\Ra \frac{1}{2}m\left(\frac{dy}{dt}\right)^2-\phi(y(t)) = E$
1.10.1 Conservation of energy
$\ds\Ra \frac{1}{2}m\left(\frac{dy}{dt}\right)^2-\phi(y(t)) = E$
$E$ is the total energy, or simply the energy of the system. It is constant in time, or in physics parlance, conserved.
To simplify the equation, let's introduce a function $V(y)=-\phi(y).$ Energy can now be expressed as the sum of two terms, rather than the difference.
$\ds \frac{1}{2}m\left(\frac{dy}{dt}\right)^2+V(y(t)) =E.$
1.10.1 Conservation of energy
$\ds \frac{1}{2}m\left(\frac{dy}{dt}\right)^2+V(y(t)) =E.$
The first term, $\displaystyle \frac{1}{2}m\left(\frac{dy}{dt}\right)^2$ is a function of velocity (mass is a constant). It is the contribution to total energy due to the object's motion, and is called kinetic energy (K.E). Note that $\left(\frac{dy}{dt}\right)^2 \ge 0,$ so kinetic energy is non-negative.
The second term, $\displaystyle V(y(t))$ is a function of position. It's called potential energy (P.E.), because it has the potential to be converted into kinetic energy.
1.10.1 Conservation of energy
There are several advantages to adopting this formalism:
1.10.1 Conservation of energy
We can re-state the result in another way.
If we write
$y_f =y(t_f)$ (Final position) | $\ds v_f =\frac{dy}{dt}(t_f)$ (Final velocity) |
$y_0 =y(t_0)$ (Initial position) | $\ds v_0=\frac{dy}{dt}(t_0)$ (Initial velocity) |
1.10.1 Conservation of energy
Then we have
$\ds\frac{1}{2}mv_f^2+V(y_f)$ $\ds =E $ $\ds=\frac{1}{2}mv_0^2+V(y_0)$
$\ds\Ra \left(\frac{1}{2}mv_f^2 - \frac{1}{2}mv_0^2 \right)$ $\ds+\,\left(V(y_f)-V(y_0)\right)$ $=0$
$\ds\Ra \Delta \text{K.E.}+\Delta \text{P.E.}=0$
This calculation shows energy cannot be created or destroyed, but can be changed from one form to another. 🤯
1.10.2 Potential Energy diagrams
Since kinetic energy is non-negative, it follows that total energy can never be less than potential energy. We can visualise this for a physical system by graphing $V$ against $y$ using a potential energy diagram:
1.10.2 Potential Energy diagrams
Since kinetic energy is non-negative, it follows that total energy can never be less than potential energy. We can visualise this for a physical system by graphing $V$ against $y$ using a potential energy diagram:
1.10.2 Potential Energy diagrams
Potential energy diagrams provide a visual description of a system's dynamics. Picture a ball rolling on the hilly terrain imposed by the shape of $V(y).$
If the ball is placed where
$V$ has an upward slope ($V^\prime(y)>0$), your
intuition tells you it should be pushed
to the left.
1.10.2 Potential Energy diagrams
Indeed, since the force is $\displaystyle F(y)=-V^\prime(y)\lt 0,$ your intuition is correct. For the same reason, if $V$ has a downward slope, it is pushed to the right.
1.10.2 Potential Energy diagrams
For $a\lt y\lt b,$ $V(y)\lt E.$ This means kinetic energy is positive, and hence motion can take place here. The same is true when $y>c$.
On the other hand, the system does not have enough energy to occupy positions $y\lt a$, or $b\lt y\lt c.$
1.10.2 Potential Energy diagrams
Positions $a,$ $b,$ and $c$ are called turning points, because an object in this system reverses its direction here. Since $V(y)=E$ at these points, kinetic energy must be zero, and so an object must be temporarily at rest.
1.10.2 Potential Energy diagrams
A system with the same potential energy function $V(y)$ but a different total energy $E^\star$ will behave differently. In the graph below, motion is only possible for \( y>e.\)
1.10.2 Potential Energy diagrams
Example: Show that equilibrium solutions to equation $m\dfrac{d^2y}{dt}= F\left(y(t)\right)$ can only occur where kinetic energy is $0,$ and only at critical points of $V$.
Let $y=\alpha $ be an equilibrium solution.
$\Ra y'=0,$ and so $\ds K.E.=\frac{1}{2}m\cdot 0^2$ $=0$
Also $y'' = 0 $ and so $\ds F = m\cdot 0$ $=-V^\prime(\alpha)$ $\Ra V^\prime(\alpha)=0$
Thus, $V^\prime(\alpha) =0,$ so $\alpha $ must be a critical point of $V.$
An object at point $d$ in the graph above, with total energy $E^\star,$ represents an equilibrium solution.
1.10.3 Conservation of energy under constant gravity
Consider the acceleration of a particle under a constant gravitational force. (For motion near the Earth's surface, gravitational force is approximately constant). By Newton's second law ($ma=F$), the ODE governing this sytem is
$\ds m \frac{d^2 y}{dt^2}=-mg.~~(\#)$
There is a negative sign on the right of this equation because gravity pulls down towards the earth, but $y$ is measured as distance up from the the earth. Now also recall that $V,$ the potential energy has been defined as the negative of the integral of the force with respect to distance.
1.10.3 Conservation of energy under constant gravity
There is a negative sign on the right of this equation because gravity pulls down towards the earth, but $y$ is measured as distance up from the the earth. Now also recall that $V,$ the potential energy has been defined as the negative of the integral of the force with respect to distance.
$V(y)\ds =-\int F(y)\, dy\qquad$
$\ds =-\int \left(-mg\right)\, dy$
$\ds =mgy+c.\quad \;$
1.10.3 Conservation of energy under constant gravity
Setting $c$ equal to $0$, the ODE ($\#$) may be written;
$\ds \frac{1}{2}m \left(\frac{dy}{dt}\right)^2+mgy = E$
As $y$ decreases, the P.E. of the particle decreases and the K.E. increases so that the total energy is conserved.
1.10.3 Conservation of energy under constant gravity
The energy diagram looks like this:
The potential energy function $V(y)$ is a straight line, with slope $mg$. Given a fixed total energy $E$, the system cannot occupy positions $y>y_{\text{max}}.$ The turning point (and thus the maximum height an object can reach) is $y_{\text{max}}$.
1.10.3 Conservation of energy under constant gravity
Example: A projectile is launched vertically upward from ground level, with initial speed $y=98$m/s. What is the maximum height the projectile will reach?
$\ds \frac{dy}{dt}(0)=98, \quad y(0)=0.$
Hence, using the energy equation, with $g=9.8,$ the initial conditions are
$\ds \frac{1}{2}m \left(98\right)^2+m\times g\times 0 $ $\ds=\frac{m}{2}98^2=E.$
1.10.3 Conservation of energy under constant gravity
$\ds \Ra \frac{m}{2}98^2=E$
At maximum height, $\ds \frac{dy}{dt}=0$ and so
$\ds \frac{1}{2}m \left(0\right)^2+mgy_{\text{max}}$ $ \ds =m\times9.8\times y_{\text{max}}$ $ \ds =E$ $ \ds=\frac{m}{2}98^2$
$\ds \Ra y_{\text{max}}=\frac{\frac{m}{2}98^2}{m\times 9.8}$ $\ds =\frac{98^2}{2\times9.8}$ $\ds =490\text{ metres.}$
1.10.4 Elastic potential energy
Recall the undamped mass-spring system with equation;
$\ds m\frac{d^2 x}{dt^2} =-kx.$
From page (76), the solution of this ODE is
$\ds x(t)=A\cos(\omega t-\phi).$
where $\omega =\sqrt{k/m}$.
1.10.4 Elastic potential energy
The potential energy corresponding to the spring force $F=-kx$ can be calculated as
$\ds V(x) =-\int (-kx)\, dx$
$\ds\;\; =\frac{1}{2}kx^2+c$
$\ds\quad\;\; =\frac{1}{2}m\omega^2 x^2+c$
1.10.4 Elastic potential energy
Setting $c=0,$ the system is characterised by the energy equation
$\ds \frac{1}{2}m\left(\frac{dx}{dt}\right)^2+\frac{1}{2}m\omega^2x^2 =E.~\clubsuit$
Since $\displaystyle x(t)=A\cos(\omega t-\phi),$ $\ds \frac{dx}{dt}=-A\omega\sin(\omega t-\phi).$ Equation $\clubsuit$ becomes
$\ds E = \frac{1}{2}m\left(A^2 \omega^2 \sin^2(\omega t-\phi) \right)$ $\ds +\frac{1}{2} m\omega^2 \left(A^2\cos^2(\omega t-\phi)\right)$
1.10.4 Elastic potential energy
Setting $c=0,$ the system is characterised by the energy equation
$\ds \frac{1}{2}m\left(\frac{dx}{dt}\right)^2+\frac{1}{2}m\omega^2x^2 =E.~\clubsuit$
Since $\displaystyle x(t)=A\cos(\omega t-\phi),$ $\ds \frac{dx}{dt}=-A\omega\sin(\omega t-\phi).$ Equation $\clubsuit$ becomes
$\ds E = \frac{1}{2}m\left(A^2 \omega^2 \sin^2(\omega t-\phi) \right)$ $\ds +\frac{1}{2} m\omega^2 \left(A^2\cos^2(\omega t-\phi)\right)$
$ \ds\quad =\frac{1}{2} m\omega^2 A^2 \big(\big. $ $ \ds \sin^2(\omega t-\phi)+\cos^2(\omega t-\phi) $ $ \ds \big.\big) $
$\ds\quad =\frac{m}{2} \omega^2 A^2 .$
1.10.4 Elastic potential energy
The energy of this system is proportional both to square of the angular frequency, and to the square of the amplitude.
The amplitude $A$ is the maximum extent of the object's oscillation.
1.10.4 Elastic potential energy
The amplitude $A$ is the maximum extent of the object's oscillation.
Hence, it represents a turning point, where $V(y)=E,$ and kinetic energy is $0$.
1.10.4 Elastic potential energy
Hence we can calculate:
$\ds V(A)$ $\ds=\frac{m}{2}\omega^2 A^2$ $\ds = E$
1.10.4 Elastic potential energy
Example: Recall the undamped, oscillating spring system from page 78. It had a mass of $9$kg, and a spring constant of $4$N/m. It was pulled down $1$m and given an initial upward kick of $-0.5$m/s.
How fast is the mass moving when it passes through its equilibrium position?
From page 78 (workbook) we determined that
$\ds \omega =2/3,$ and $\ds A=\frac{5}{4}.$ Recall also $\ds E=\frac{m}{2}\omega^2 A^2.$
1.10.4 Elastic potential energy
From page 78 (workbook) we determined that
$\ds \omega =2/3,\,$ and $\,\ds A=\frac{5}{4}.$ Recall also $\,\ds E=\frac{m}{2}\omega^2 A^2.$
At the equilibrium position
$\ds E =\frac{m}{2}\left(\frac{dx}{dt}\right)^2+\frac{m}{2}\omega^2\times 0^2$ $=\ds\frac{m}{2}\left(\frac{dx}{dt}\right)^2$
Hence $\ds \; \frac{m}{2}\omega^2A^2$ $\ds=\frac{m}{2}\left(\frac{dx}{dt}\right)^2$
$\ds \Ra \frac{dx}{dt}=\pm\sqrt{\omega^2A^2}$ $=\pm\omega A$ $\ds=\pm\left(\frac{2}{3}\right)\left(\frac{5}{4}\right)$ $\ds=\pm\frac{5}{6}.$
1.10.5 Main points