Multivariate Calculus
&
Ordinary Differential Equations

Lecture The Energy principle



1 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

1.10.1 Conservation of energy

There are several advantages to adopting this formalism:

  1. Newton's second law of motion is a second order ODE. The above technique allows any system with a continuous, position-dependent force to be described as a separable, first order ODE. Any such system is solvable in principle.
  2. Energy $E$ is seen as an integration constant, calculated from initial position and velocity.
  3. Given knowledge of the energy, position can be calculated from velocity, and vice-versa.


1 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

1.10.4 Elastic potential energy

Hence we can calculate:

$\ds V(A)$ $\ds=\frac{m}{2}\omega^2 A^2$ $\ds = E$



1 Ordinary Differential Equations

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 Ordinary Differential Equations

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 Ordinary Differential Equations

1.10.5 Main points

  • You should understand how to find the energy equation for any position-dependent force.
  • You should be able to determine position and velocity from such an equation.
  • You should also be able to draw a potential energy diagram, for a position-dependent force.
  • You should be able to interpret system behaviour and determine equilibrium solutions from an energy diagram.


Credits