Kvant

Rough automatic draft, low quality: expect mistakes and awkward wording. It will be replaced by a better translation. Original →

20

This article talks about some of the properties of the ellipse and the different laws of planetary motion, the most important of which are called Kepler’s laws.

Gravity
The law of inertia now seems to be the most elementary and natural law - the study of mechanics will begin with it. However, in order to understand it, it was necessary to distract ourselves from what was happening directly before our eyes. Under normal conditions, the movement of the body in the horizontal plane is hampered by friction and air resistance; When moving in a vertical plane, the forces of air resistance and gravity act on the body. It was necessary to get used to their influence, at least mentally. It took many years.

Only after the law of inertia—the law of body motion under idealized conditions—was understood, could it be possible to consider the movement of bodies under the action of forces.

The law of falling bodies on Earth was established by Galileo (1564-1642). He was the first to set himself the task of understanding how the speed of a falling body changes in successive periods of time. Thus, he introduced acceleration into mechanics and established that the accelerations of all bodies falling to the ground are identical and constant. But the fact that the cause of acceleration is gravity was only later understood by Newton.

Almost this years after the publication of Galileo’s work, Newton discovered the law of universal gravitation (1687), according to which the force of gravity is not constant, but decreases with increasing distance between bodies as the square of this distance.

If physics had historically developed in the order in which it is presented in textbooks, then only after Newton’s discovery would it be possible to move on to the study of the laws of planetary motion. But the development of science does not follow the logic of the textbook. Much earlier than Newton, these laws were established by Galileo’s contemporary Johannes Kepler (1571-1630).

Kepler’s meanings in mechanics were extremely meager in comparison with those of a modern schoolchild. Not only did he not know Newton’s laws (they had not yet been discovered), but he also had a very vague idea of ​​what force was. All that Kepler had in his hands were very neat ancient observations of the movements of celestial bodies, which his teacher Tycho Brahe had conducted for many years *) and he himself continued to conduct. Only extraordinary hard work and confidence in the existence of simple laws governing the movement of the planets allowed him to discover the laws that have since bear his name.

Kepler’s laws are not so strong even now. This is mainly due to the fact that ellipses - the curves along which the planets move - are poorly understood by the average reader. In reality, obtaining Kepler’s laws is not so difficult. To do this, you need to use two laws of mechanics: the law of conservation of energy (which was not known to either Kepler or Newton) and the law of conservation of angular momentum. The latter, as we will see, is equivalent to Kepler’s second law: the areas slowed down by the planet’s radius vector over equal periods of time are the same.

*) Tycho Brahe was one of the last astronomers to observe the movement of celestial bodies with the naked eye.

22

Law of Conservation of Energy

On the surface of the Earth the force of gravity is mgmg. If we introduce the potential energy of a body into the gravitational field U=mghU = mgh (h is the distance from the body to the Earth), then we can say that when the body moves, the sum of its kinetic and potential energy E=mv22+mghE = \frac{mv^2}{2} + mgh does not change.

However, the expression for UU that we wrote is valid only at small distances from the body to the Earth, when hh is small compared to the radius of the Earth and gravity can be considered constant. If the height becomes greater, then the force of gravity decreases with distance. According to Newton’s law

∣F∣=γmMR2,|F| = \gamma \frac{mM}{R^2},

where γ\gamma is the gravitational constant, MM is the mass of the Earth and RR is the distance from the body to the center of the Earth.

If we use the exact expression for gravity, then the potential energy must be replaced by U=−γmMRU = -\gamma \frac{mM}{R}. Let’s check it out. Let us remember that the work done by a body during movement is equal to the decrease.

Fig. 2.

Fig. 3.

reducing its potential energy. We calculate the work done by the body feeding from the height R2R_2 to the height R1R_1 from the center of the Earth (Fig. 2). To do this, you need to calculate the difference in potential energies of the body:

работа=U2−U1=\text{работа} = U_2 - U_1 =

=γmM(−1R2+1R1)== \gamma m M \left( -\frac{1}{R_2} + \frac{1}{R_1} \right) =

=γmM(R2−R1R1R2).= \gamma m M \left( \frac{R_2 - R_1}{R_1 R_2} \right).

If R2−R1≪R1R_2 - R_1 \ll R_1 (R2≈R1R_2 \approx R_1), then the product R1R2R_1 R_2 can be replaced by R2R^2 - the square of the geometric mean R1R_1 and R2R_2 (for R2≈R1R_2 \approx R_1: R≈R2≈R1R \approx R_2 \approx R_1), and R2−R1R_2 - R_1 can be replaced by ll - the path traversed by the body. Then the last formula takes the form

работа=сила×путь=γmMR2l.\text{работа} = \text{сила} \times \text{путь} = \gamma \frac{mM}{R^2} l.

force path of gravity

This means that the expression for UU, which we wrote above, is the real potential energy of the body.

Now we need to write the formula for the total energy of the body differently:

E=mv22−γmMR=m(v22−γMR).E = \frac{mv^2}{2} - \gamma \frac{mM}{R} = m \left( \frac{v^2}{2} - \gamma \frac{M}{R} \right).

Ellipse

Let’s start with the general definition of an ellipse (Fig. 3). An ellipse is the geometric locus of lobes, the sum of the distances of which to two lobes (the foci of the ellipse) is constructed.

23

Fig. 4.

From this definition we can deduce an interesting “optical” property of the ellipse: if a light source is placed at one of the focuses of the mirror ellipse ), then all the rays coming out of it will be collected at another focus of the ellipse*). This means that the angles between the straight lines going to the point of the ellipse from its foci and the tangent to the ellipse are equal to each other (Fig. 4).

Let us prove that the ellipse has one more property, which is even convenient to take as the definition of an ellipse. To do this, we draw a tangent tt to the ellipse at point AA and lower the redistributions from both foci of the ellipse to this tangent (Fig. 5). We extend one of the redistributions beyond the straight line tt to a distance F2′N2F'_2 N_2 equal to N2F2N_2 F_2. Let us now connect the points F2′F'_2 and F1F_1 with the point AA. Let us prove that F2′AF1F'_2 AF_1 is a straight line.

This is immediately clear from the fact that α1=α2\alpha_1 = \alpha_2 (by construction) and α1=α3\alpha_1 = \alpha_3 (according to the optical property of the ellipse). Therefore, α2=α3\alpha_2 = \alpha_3, and therefore F2′AF1F'_2 AF_1, is a straight line.

*) An ellipside is the geometric locus of a point, the sum of the distances from which to two fixed points is constructed. Such a surface is obtained by growing an ellipse around its axis.

**) The easiest way to prove this property of an ellipse is to use Fermat’s principle (see «Kvant» No. 11, 1970).

Fig. 5.

Further, since AF2′=AF2AF'_2 = AF_2, then F1F2′=F1A+AF2′=F1A+AF2F_1 F'_2 = F_1 A + AF'_2 = F_1 A + AF_2, that is, the segment F1F2′F_1 F'_2 is equal to the sum of the distances from the ellipse point to the foci. Let’s denote this sum 2a2a. Let us now lower the redistribution F1BF_1 B onto the line F2′F2′F'_2 F'_2, calculate the length of the leg F1BF_1 B from two right triangles - brown and yellow - and equate both expressions for F1BF_1 B to each other:

(2a)2−(ρ1+ρ2)2=(2c)2−(ρ2−ρ1)2(2a)^2 - (\rho_1 + \rho_2)^2 = (2c)^2 - (\rho_2 - \rho_1)^2.

Opening the brackets, we get

ρ1ρ2=a2−c2\rho_1 \rho_2 = a^2 - c^2.

So, the product ρ1ρ2\rho_1 \rho_2 does not depend on the point at which the tangent is drawn. Let’s draw it at the end of the semi-minor axis (Fig. 6). In this case a2−b2=c2a^2 - b^2 = c^2. Therefore ρ1ρ2=b2\rho_1 \rho_2 = b^2.

Thus, we got the main theorem:

The product of the distances from two foci to any tangent of the ellipse is a constant value equal to the square of the semi-minor axis.

24

Kepler’s first law

Now we can consider the movement of the planet. For simplicity, we will not prove that the planet’s orbit is an ellipse, but, assuming this and using the laws of conservation of the moment of motion and energy, we will find it and make sure that the movement along ellipses in a gravitational field does not contradict the laws of conservation.

Let us note two positions of the planet AA and A′A', symmetrical relative to the center of the ellipse (Fig. 7). We can say it another way: point AA is the true position of the planet at some point in time, and A′A' is obtained from the first by transferring the Sun to another focus, while the planet remains in the same place.

Let’s write down the law of conservation of energy. Marking the quantities with index 1 for point AA and index 2 for point A′A', we obtain (E1=E2=E)∗(E_1 = E_2 = E)^*.

−E=mv122−γmMr1=mv222−γmMr2.-E = \frac{mv_1^2}{2} - \gamma \frac{mM}{r_1} = \frac{mv_2^2}{2} - \gamma \frac{mM}{r_2}.

The angular momentum of the planet is equal to the product of the amount of motion mv1mv_1 (or mv2mv_2) and the “play” - the impact parameter ρ1\rho_1 (or ρ2\rho_2), while the total energy is T+U=−E<0T + U = -E < 0.

How is the angular momentum for both positions of the planet the same:
L=mv1ρ1=mv2ρ2.L = mv_1\rho_1 = mv_2\rho_2.
Let us calculate the product of the kinetic energies of the planet at points AA and A′A':
T1T2=mv122⋅mv222=(−E+αr1)×(−E+αr2).T_1T_2 = \frac{mv_1^2}{2} \cdot \frac{mv_2^2}{2} = \left(-E + \frac{\alpha}{r_1}\right) \times \left(-E + \frac{\alpha}{r_2}\right).
(Here α=γMm\alpha = \gamma Mm.)
Replacing the velocities in this equation with their expressions through the moment of quantity of motion, after simple transformations we obtain:
14L4m2ρ12ρ22−E2=−α(r1+r2)E+α2r1r2.\frac{1}{4} \frac{L^4}{m^2\rho_1^2\rho_2^2} - E^2 = \frac{-\alpha(r_1 + r_2)E + \alpha^2}{r_1r_2}.

It is easy to see that this equality holds if the planet’s orbit is an ellipse. Indeed, for this you need to put
r1+r2=αE=2ar_1 + r_2 = \frac{\alpha}{E} = 2a
And
ρ1ρ2=L22mE=b2.\rho_1\rho_2 = \frac{L^2}{2mE} = b^2.

Both of these equalities are nothing more than the two definitions of an ellipse given above.

Thus, our trajectory is an ellipse with semi-axes: large
a=γmM2E и малой b=L2mE.a = \gamma \frac{mM}{2E} \text{ и малой } b = \frac{L}{\sqrt{2mE}}.

26

Fig. 8.

Then the speed of the planet will be dimensionless. Speed ​​is a number that shows how many times v0v_0 is contained in the speed of the planet:

β1=v1v0иβ2=v2v0.\beta_1 = \frac{v_1}{v_0} \quad \text{и} \quad \beta_2 = \frac{v_2}{v_0}.

According to the previous speeds β1\beta_1 and β2\beta_2 are related by the relation

β1=1β2.\beta_1 = \frac{1}{\beta_2}.

What does it mean? If a particle moves at some point in the orbit with a speed of β\beta, then at a point symmetric relative to the center of the orbit, its speed will be equal to 1β\frac{1}{\beta}

We described the trajectory of the planet, specifying its position (coordinates) at each moment in time. You can do it differently - set the speed of the planet. The curve that we get in this case is called the velocity graph (we can say that this is a trajectory in velocity space). Since the product of the segments β1\beta_1 and β2\beta_2 has been constructed, the speed graph of the planet is a circle (Fig. 8) (remember the theorem: “A point is given inside a circle. Then for all chords passing through this point, the product of the lengths of the segments on which this point divides the chords are the same”).

If we replace the speed β\beta with 1β\frac{1}{\beta} at each point of the godgraph, then the godgraph will not change, but simply the surface speed at 180∘180^\circ. This property is called relative surface motion invariance.

If the planet’s trajectory is circular, then β=1\beta = 1 (prove it!). In this case, the godgraph is a circle with its center at the origin.

Circulation period

In order to calculate the orbital period - the time of a complete revolution of a planet around Coln or an artificial satellite around the Earth, it is necessary to calculate the time of passage of each section of the trajectory (separately, since the speed changes all the time) and then sum it over all sections. Of course, this cannot be done simply, but all difficulties can be circumvented.

From Kepler’s second law we know that in time Δt\Delta t the radius vector of the planet “sweeps” an area equal to L2mΔt\frac{L}{2m} \Delta t (recall that L=constL = \text{const}). This means that the period of revolution can be calculated by dividing the area of ​​the ellipse by the constant speed of “sweeping”

T=SэллипсаL/2m.T = \frac{S_{\text{эллипса}}}{L/2m}.

Since the ellipse is obtained from a circle sloping in one direction ba\frac{b}{a} times∗^*), then the area of ​​the ellipse is equal to

Sэллипса=Sкрута⋅ba=πa2ba=πab.S_{\text{эллипса}} = S_{\text{крута}} \cdot \frac{b}{a} = \pi a^2 \frac{b}{a} = \pi ab.

From these formulas we get

T=2πabL.T = \frac{2\pi ab}{L}.

But b=L2Emb = \frac{L}{\sqrt{2Em}}, therefore

T=πa2mET = \pi a \sqrt{\frac{2m}{E}}

∗^*) See the article by I. N. Bronstein “Ellipse” («Kvant» No. 9, 1970).

27

or, substituting here the expressions for a(a=γMm2E)a \left( a = \gamma \frac{Mm}{2E} \right) and replacing (2Em)3/2\left( \frac{2E}{m} \right)^{3/2} with v03v_0^3, we get
T=2γMm2Em2E=2πγMv03.T = 2\gamma M \frac{m}{2E} \sqrt{\frac{m}{2E}} = 2 \frac{\pi \gamma M}{v_0^3}.
From this formula it is clear that the period of revolution depends only on the value v0v_0. If several identical satellites are released from the same point above the Earth’s surface in different directions, but with the same speed v0v_0, then all of them will simultaneously (after a time TT) converge at the same point.

Kepler’s third law
If now in the formula (*) for the period of revolution of the planet together the energy EE substitute her expression with a larger one (E=γmM2a)\left( E = \gamma \frac{mM}{2a} \right), then we get
T=2πa3/21γM.T = 2\pi a^{3/2} \frac{1}{\sqrt{\gamma M}}.
Therefore, T2T^2 is proportional to a3a^3, that is, the square of the planet’s orbital period is proportional to the cube of the large half-ellipse. And this is Kepler’s third law.

Fig. 9.
An ellipse is the geometric locus of the centers of a circle passing through the focus and touching a given circle (guide circle).

In this article we talked a little about geometry, optics and mechanics. Now solve the problems.

Optical problem
In one of the focuses of the ellipson, the inner surface of which is mirrored, there is an exact light source. Construct an image of this source for an observer located at the second focus of the ellipson.

Geometric problem
Show that the ellipse is the locus of the centers of a circle passing through the focus and touching a given circle (guide circle) (Fig. 9).

At figure it is clear that the bases of the perpendiculars ρ1\rho_1 and ρ2\rho_2 lie on the “principal circle”. Are you listening to this?

Space challenge
The orbit of the Earth’s satellite has the following parameters: greatest distance RR, shortest distance rr. Calculate the total energy and angular momentum of the satellite. Find the speed of the satellite at the ends of the major and minor axes.

Beautiful formation
The drawing on the cover of the magazine is constructed like this: a point is selected inside a circle, and a circle smaller than it is drawn from it, as from the center. Rays are drawn from the center of this second circle until they meet the larger circle. Through a point on each ray lying in the middle between two circles, a chord is drawn perpendicular to this ray. The spot inside has the shape of an ellipse.