index 10. Mechanical similarity

10. Mechanical similarity


Multiplication of the Lagrangian by any constant clearly does not affect the equations of motion. This fact (already mentioned in §2. The principle of least action) makes possible, in a number of important cases, some useful inferences concerning the properties of the motion, without the necessity of actually integrating the equations.

Such cases include those where the potential energy is a homogeneous function of the co-ordinates, i.e. satisfies the condition

\[\label{equations/10.1}\tag*{(10.1)}U(\alpha\mathbf{r}_1,\alpha\mathbf{r}_2,\ldots,\alpha\mathbf{r}_n)=\alpha^kU(\mathbf{r}_1,\mathbf{r}_2,\ldots,\mathbf{r}_n)\]

where \(\alpha\) is any constant and \(k\) the degree of homogeneity of the function.

Let us carry out a transformation in which the co-ordinates are changed by a factor \(\alpha\) and the time by a factor \(\beta: \boldsymbol{r}_a\rightarrow \alpha\boldsymbol{r}_a, t\rightarrow \beta t\). Then all the velocities \(\boldsymbol{v}_a = \mathrm{d}^{}\boldsymbol{r}_a/\mathrm{d}^{}t\) are changed by a factor \(\alpha/\beta\), and the kinetic energy by a factor \(\alpha^2/\beta^2\). The potential energy is multiplied by \(\alpha^k\). If \(\alpha\) and \(\beta\) are such that \(\alpha^2/\beta^2 = \alpha^k\), i.e. \(\beta = \alpha^{1-k/2}\), then the result of the transformation is to multiply the Lagrangian by the constant factor \(\alpha^k\), i.e. to leave the equations of motion unaltered.

A change of all the co-ordinates of the particles by the same factor corresponds to the replacement of the paths of the particles by other paths, geometrically similar but differing in size. Thus we conclude that, if the potential energy of the system is a homogeneous function of degree k in the (Cartesian) co-ordinates, the equations of motion permit a series of geometrically similar paths, and the times of the motion between corresponding points are in the ratio

\[\label{equations/10.2}\tag*{(10.2)}t'/t=(l'/l)^{1-k/2}\]

where \(l'/l\) is the ratio of linear dimensions of the two paths. Not only the times but also any mechanical quantities at corresponding points at corresponding times are in a ratio which is a power of \(l'/l\). For example, the velocities, energies and angular momenta are such that

\[\label{equations/10.3}\tag*{(10.3)}v'/v=(l'/l)^{k/2},\quad E'/E=(l'/l)^{k},\quad M'/M=(l'/l)^{1+k/2}\quad\]

The following are some examples of the foregoing.

As we shall see later, in small oscillations the potential energy is a quadratic function of the co-ordinates (\(k=2\)). From (10.2) \(t'/t=(l'/l)^{1-k/2}\) we find that the period of such oscillations is independent of their amplitude.

In a uniform field of force, the potential energy is a linear function of the co-ordinates (see (5.8) \(U=-\mathbf{F}\cdot\mathbf{r}\) ), i.e. \(k = 1\). From (10.2) \(t'/t=(l'/l)^{1-k/2}\) we have \(t'/t = \sqrt{l'/l}\). Hence, for example, it follows that, in fall under gravity, the time of fall is as 2the square root of the initial altitude.

In the Newtonian attraction of two masses or the Coulomb interaction of two charges, the potential energy is inversely proportional to the distance apart, i.e. it is a homogeneous function of degree \(k = 1\). Then \(t'/t = (l'/l)^{3/2}\), and we can state, for instance, that the square of the time of revolution in the orbit is as the cube of the size of the orbit (Kepler’s third law).

If the potential energy is a homogeneous function of the co-ordinates and the motion takes place in a finite region of space, there is a very simple relation between the time average values of the kinetic and potential energies, known as the virial theorem.

Since the kinetic energy \(T\) is a quadratic function of the velocities, we have by Euler’s theorem on homogeneous functions \(\sum_a \boldsymbol{v}a\cdot\partial T/\partial\boldsymbol{v}_a = 2T\), or, putting \(\partial T/\partial \boldsymbol{v}_a = \boldsymbol{p}_a\), the momentum,

\[\label{equations/10.4}\tag*{(10.4)}2T=\sum_a\mathbf{p}_a\mathbf{v}_a=\frac{\mathrm{d}^{}}{\mathrm{d}^{}t}(\sum_a\mathbf{p}_a\mathbf{r}_a)-\sum_a\mathbf{r}_a\dot{\mathbf{p}}_a\]

Let us average this equation with respect to time. The average value of any function of time \(f(t)\) is defined as

\[ \bar{f}=\lim_{\tau\rightarrow\inf}\frac{1}{\tau}\int_0^\tau f(t)\mathrm{d}^{}t \]

It is easy to see that, if \(f(t)\) is the time derivative \(\mathrm{d}^{}F(t)/\mathrm{d}^{}t\) of a bounded function \(F(t)\), its mean value is zero. For

\[ \bar{f} =\lim_{\tau\rightarrow\inf}\frac{1}{\tau}\int_0^\tau \frac{\mathrm{d}^{}F}{\mathrm{d}^{}t}\mathrm{d}^{}t =\lim_{\tau\rightarrow\inf}\frac{F(\tau)-F(0)}{\tau} =0 \]

Let us assume that the system executes a motion in a finite region of space and with finite velocities. Then \(\sum_a\boldsymbol{p}_a\cdot\boldsymbol{r}_a\) is bounded, and the mean value of the first term on the right-hand side of (10.4) \(2T=\sum_a\mathbf{p}_a\mathbf{v}_a=\frac{\mathrm{d}^{}}{\mathrm{d}^{}t}(\sum_a\mathbf{p}_a\mathbf{r}_a)-\sum_a\mathbf{r}_a\dot{\mathbf{p}}_a\) is zero. In the second term we replace \(\boldsymbol{\dot{p}}_a\) by \(-\partial U/\partial \boldsymbol{r}_a\) in accordance with Newton’s equations (5.3) \(m_a\mathrm{d}^{}\mathbf{v}_a/\mathrm{d}^{}t=-\partial U/\partial \mathbf{r}_a\) , obtaining1

\[\label{equations/10.5}\tag*{(10.5)}2\overline{T}=\overline{\sum_a\mathbf{r}_a\cdot\partial U/\partial \mathbf{r}_a}\]

If the potential energy is a homogeneous function of degree \(k\) in the radius vectors \(\boldsymbol{r}_a\), then by Euler’s theorem equation (10.5) \(2\overline{T}=\overline{\sum_a\mathbf{r}_a\cdot\partial U/\partial \mathbf{r}_a}\) becomes the required relation:

\[\label{equations/10.6}\tag*{(10.6)}2\overline{T}=k\overline{U}\]

Since \(\bar{T}+\bar{U}=\bar{E}=E\), the relation (10.6) \(2\overline{T}=k\overline{U}\) can also be expressed as

\[\label{equations/10.7}\tag*{(10.7)}\overline{U}=2E/(k+2),\qquad \overline{T}=kE/(k+2)\]

which express \(\bar{U}\) and \(\bar{T}\) in terms of the total energy of the system.

In particular, for small oscillations (\(k=2\)) we have \(\bar{T} = \bar{U}\), i.e. the mean values of the kinetic and potential energies are equal. For a Newtonian interaction \((k = - 1)\), \(2\bar{T} = - \bar{U}\), and \(E=-\bar{T}\), in accordance with the fact that, in such an interaction, the motion takes place in a finite region of space only if the total energy is negative (see §15 Kepler's problem).