Orbital motions and the Centre of Mass

Orbital motions of astronomical bodies and their Centre of Mass from different reference frames: a conceptual step between the Geocentric and Heliocentric models

André G C Guerra111Corresponding author. and Paulo Simeão Carvalho Departamento de Física e Astronomia, Centro de Física do Porto, Faculdade de Ciências, Universidade do Porto, Rua do Campo Alegre, s/n, 4169-007 Porto, Portugal Departamento de Física e Astronomia, UEC, Faculdade de Ciências, Universidade do Porto, Rua do Campo Alegre, s/n, 4169-007 Porto, Portugal IFIMUP-IN, Rua do Campo Alegre, s/n, 4169-007 Porto, Portugal aguerra@fc.up.pt psimeao@fc.up.pt
Abstract

The motion of astronomical bodies and the centre of mass of the system is not always well perceived by students. One of the struggles is the conceptual change of reference frame, which is the same that held back the acceptance of the Heliocentric model over the Geocentric one. To address the question, the notion of centre of mass, motion equations (and their numerical solution for a system of multiple bodies), and change of referential is introduced. The discussion is done based on conceptual and real world examples, using the solar system. Consequently, through the use of simple “do it yourself” methods and basic equations, students can debate complex motions, and have a wider and potentially effective understanding of physics.

: \PED

,

  • May 2016

Keywords: Reference Frame, Centre of Mass, Motion Equations, Numerical Solution, Geocentric Model, Heliocentric Model

1 Introduction

The notion of referential is one of the main causes why students have difficulties interpreting motions in astronomy. This is not completely surprising, since the change of referential involves abstract reasoning, which is usually not completely developed in all population, especially young students [1]. Indeed, several examples of misconceptions were identified among students in many contents of physics and astronomy that require formal thinking [2, 3].

Looking back to the history of physics, the meticulous observations of nature by Aristotle (384 – 322 BC) lead him to a philosophical description of the universe. This was based on a simple, logic and very appealing common sense language. For Aristotle, the Earth was the centre of the universe; the Earth was the only referential he knew. Such geocentric description was lately improved by Ptolemy (c. 100 – c. 170), in his work Almagest [4, 5, 6, 7], where the motions of the planets, Moon and Sun around the Earth, were presented as a physical model supported by mathematical descriptions of the orbits.

The geocentric model had several advantages at that time: i. it did not contradict Aristotle’s descriptions, ii. it was supported by the religious ideals of the middle age in Europe, and iii. it was based on a local frame of reference (the Earth) and followed a concrete reasoning. That is why it prevailed for about 14 centuries.

When Copernicus (1473 – 1543) took the (previously) rejected ideas of Aristarchus of Samos (c. 310 – c. 230 BC) and proposed the Heliocentric model, he and his followers quickly realized they had two major issues to overcome: a huge and visible religious dispute, and a not-so-visible but still important cognitive barrier. The cognitive barrier arose because it was necessary to put oneself outside the Earth to fully understand the observations taken from it. It was this combination of factors that maintained the debate of Geocentric and Heliocentric-based models up to the moment when humankind sent satellites out of the Earth in the 20th century.

Nowadays, students still remain confused whenever changes of referential are needed to understand astronomical phenomena, such as moon phases, seasons, eclipses, or the motion of planets around a common centre of mass (CM).

Some simple but potentially effective ideas concerning the topic of referential were developed during the 11th Summer School of Physics of University of Porto for high secondary level students, and are presented here. The concepts addressed are: CM of a system, trajectory of bodies in the CM referential, and absolute trajectory of the same bodies. In this approach, some computing is used to simulate the trajectories of the bodies, and to discuss the implications of the initial conditions (mass and speed) on the description of the phenomena.

2 Theory and Computational Method

Before diving into the subject of referential changes, it is important to introduce the theory needed for solving the problems numerically [8].

Consider a system composed by a number of distinct bodies, randomly distributed, orbiting each other. For this setup, it is useful to define the location of the CM of the system. The bodies’ motion due to gravity forces can be computed with respect to a system centered on the CM, or to an inertial reference system.

2.1 Centre of mass of a system

The CM of a multiple bodies system can be defined as the unique point, within the system, which can be used to describe the system’s response to external forces.

The mass of the CM () is the sum () of the mass of all bodies, each with mass ,

(1)

The position () and velocity () of the CM are related to each body location () and velocity (), and can be computed by

(2)
(3)

2.2 Motion Equations

To calculate the motion of the bodies in the system, Newton’s second law has to be used,

(4)

that relates the forces acting on a body () and its acceleration ().

Admitting two bodies (of mass and ) in free space, the gravitational force acting on body 1 is

(5)

where is the unit vector that points from body 1 to body 2, and is the gravitational constant. Using the notion of unit vector in equation 5, the result is

(6)

The acceleration of body 1 can thus be obtain from both equations 4 and 6,

(7)

The same can be applied to body 2, just changing to and inverting the direction of in equation 7.

2.2.1 Numerical Solution

There are many methods to solve equation 7. Here a numerical one is used because a general solution for a system of multiple bodies is searched, and there is no analytical solution.

The simplest numerical method is the Euler one (and can even be applied with a spreadsheet editor [9]). It states that for a general function , it is possible to compute its value at an instant , with an initial value (at instant ), using the relation

(8)

where is the function derivative at instant , and .

Applying this method to the multiple body system, and knowing that , the velocity of a body can be estimated (according to equation 8) by

(9)

where from equation 7, the acceleration for body , subject to the gravity of the others, is

(10)

where is the vector pointing from body to body . Applying the same scheme to the position (), it yields

(11)

The Euler method is considered a first-order numerical method, where the global error is proportional to the step size (). As the examples latter discussed involve astronomical bodies (where the steps used can be of the order of days), an adjustment to the above formulae can be made, to have a greater accuracy. One simple adjustment is to compute the position with a second order derivative, using the following equation,

(12)

The above sequence of equations (910 and 12) is repeated successively, until a pre-determined final instant is reached. A way to verify how much error is introduced by the numerical solution, is to calculate the variation in the total mechanical energy [10].

2.3 Frames of Reference

All previous equations describe the motion in an inertial reference system. This frame does not have translational or rotational acceleration, relative to the “fixed stars” [11]. The Heliocentric frame is an example of an inertial referential for the Earth-Moon system, where the Sun is at the origin of the reference frame.

Sometimes, the motion of the bodies can be simplified when described at the CM reference system (where the origin is at the CM). In this case, it is necessary to impose mathematically that is zero, as well as . This implies that the velocities and position of the bodies in the CM referential, respectively and , are related to the absolute ones by,

(13)
(14)

This change of reference frame is a simple translation of the inertial reference frame, and can be applied to any general point by replacing and by the respective velocity and position of that point. No rotation is introduced here.

3 Examples

With the equations presented in the previous section, several multiple bodies systems can be solved. Some examples, explored with students at the 11th Summer School of Physics, are presented here.

The first one is a generic two body system (which can be two stars or two planets). The next ones correspond to real examples in the solar system, which include the Sun, and the planets Earth, Moon, Mars and Jupiter. To deal with real examples, the Jet Propulsion Laboratory HORIZON Ephemerides, corresponding to 11th April 2016, were used to extract the bodies’ real position and velocity [12].

3.1 Example 1 – Two Bodies General System

The CM of a two bodies system is always lying between them, closer to the more massive body, or in the middle if they have the same mass. Admitting one of the bodies has an initial velocity (table 1), students easily understand that both will start moving around each other, due to the gravitational interaction between them. The trajectory of the bodies can be computed with a simple spreadsheet editor, as can be consulted in reference [13].

Body Mass [] Position (, ) [] Velocity (, ) []
1 59.72 (0, 0) (0, 0)
2 7.346 (357.0, 0) (0, 0.2500)
Table 1: Two bodies system initial conditions

It is therefore not surprising that students conclude the both bodies describe, in the CM reference system, closed elliptical trajectories around a common point, the CM (figure 1).

Figure 1: Two body system motions on the CM reference frame. Data was computed for a time interval of 90 days (only 43 days are shown here).

However, as the bodies’ motion influence the velocity of the CM, the result in the inertial referential is a helical-like motion of the lower massive body, while the more massive one experiences a composition of translation in space and rotation around the CM (figure 2).

Figure 2: Two body system motions on an inertial reference frame. Data was computed for a time interval of 90 days.

Such weird trajectories are counter-intuitive and provide a good first discussion about why we observe unlike trajectories in different reference systems.

3.2 Example 2 – Earth-Moon System

The Earth-Moon system has some similarities to the previous example, but now the Earth is about ten times more massive (the Moon’s mass is 1.23% of the Earth). Therefore, the CM of this system is much closer to the Earth.

In the CM referential, the Moon describes an elliptical motion, whereas the Earth does not seem to move (figure (a)a). However, a closer observation of the Earth reveals that in fact it moves, but the CM of the system is always located inside the planet (figure (b)b).

(a) Earth and Moon seen from the CM
(b) The CM seen from Earth
Figure 3: The Earth-Moon system motion. In (a) the Earth graphically overlaps the CM. Data was computed for a time interval of 90 days (only 27 days are shown here).

Because the CM lies inside the Earth, the first astronomers had the perception that the Moon orbited the Earth and not a common point (the CM)! That only happened because their astronomical observations were taken from the Earth local referential!

3.3 Example 3 – Sun-Earth-Moon System

Moving from the Earth-Moon system to the more general Sun-Earth-Moon system, we have to considerer that the Earth, in fact, orbits around the Sun and therefore has already a velocity in space.

At this stage, students are already aware that the most massive body (the Sun) is affected by the gravitational forces of the other two (Earth and Moon). Thus, taking into account the previous example, they can predict that the centre of the Sun experiments some motion on the inertial reference system, which means that the solar system as a whole moves in space! What they really may find surprising is that this motion is very small when compared to the Sun dimension (which has about 99.8% of the total solar system mass), and this is why we state that the planets orbit only around the Sun in the solar system.

The relative motions of the Earth and Moon can also be observed from the CM of this system of bodies. Figure 4 shows a zoom of the trajectories of these bodies while they orbit the Sun. Students can see that the Moon’s trajectory is a composition of its orbit around the Earth, and the Earth around the Sun. Therefore, the Moon intersects the trajectory of the Earth each time it completes an orbital revolution, resulting in an oscillating motion around the Earth trajectory.

The final stage of this topic is to let students understand why the (apparent) retrograde motion is seen from the Earth. For this purpose, a more complex system must be projected.

Figure 4: A close-up of the Earth-Moon system, orbiting the Sun, seen from an inertial reference frame. Data was computed for a time interval of 400 days (only a few days are shown here).

3.4 Example 4 – More Complex Systems

The principles for the mechanics of a more complex system do not differ from those in the previous examples.

The example considered here is the motions of the Earth, Mars and Jupiter around the Sun. The resulting trajectories of these bodies on an inertial reference frame can be seen in figure 5.

Figure 5: The motion of the Sun-Earth-Mars-Jupiter system, as seen from an inertial reference frame. In this image the Sun graphically overlaps the CM. Data was computed for a time interval of 4333 days.

Even though they influence each other, through gravitational forces, the planets follow almost circular orbits around the Sun, as was expected. However, this is only true when seen from an inertial reference frame, or the Sun since his motion is negligible when compared to the planets’ ones.

To evaluate the (apparent) motion of those planets as seen from a local referential (the Earth), students only need to apply equations (13) and (14), replacing and by and , respectively, to obtain the velocity and position in the Earth reference frame.

Plots for the case of Sun-Earth-Mars and Sun-Earth-Jupiter systems are given in figures 6 and 7, respectively.

Figure 6: The motion of the Sun-Earth-Mars system, as seen from the Earth. In this image the Sun graphically overlaps the CM. Data was computed for a time interval of 2748 days (only 1776 days are shown here).
Figure 7: The motion of the Sun-Earth-Jupiter system, as seen from the Earth. In this image the Sun graphically overlaps the CM. Data was computed for a time interval of 4333 days.

The loops observed on Mars and Jupiter trajectories correspond to apparent retrograde motions. In ancient times these loops lead Ptolemy to formulate the epicycle theory for the Geocentric model, as an attempt to explain the phenomena observed for the five known planets at that time: Mercury, Venus, Mars, Jupiter and Saturn. Instead of making the Geocentric model more accepted, the complexity of such theory weakened it.

Students can thus realise that the Heliocentric model proposed by Copernicus offers a simpler explanation of the planets’ motion, clearly revealing why the retrograde motions are observed from the Earth.

To obtain all these plots an educational program [14] was developed (for a x64 Windows operating system), based on MATLAB software. This program allows to see step by step the bodies’ motions, which makes it easier and more engaging for students to understand, especially when they observe Mars and Jupiter making retrograde motions, as seen from the Earth.

4 Conclusion

Students usually give very little importance to the meaning of Centre of Mass of a system, because they tend to confine it to the particle model problems. The approach described in the context of astronomy gives a new and wider perspective of understanding of the role of CM in the dynamics of a system of bodies. In particular, it promotes students’ reasoning about how the trajectories of bodies can be seen in different reference frames, which is essential to understand natural phenomena such as the Coriolis effect or the Earth tides.

This approach also allows students to realise that apparently complex motions can be explained with basic and simple laws. In fact, the most interesting about this approach is that students, with very simple mathematics but significant physics reasoning, can obtain these results by themselves, only considering gravitational interactions and no pre-established formulas.

The use of a simple spreadsheet where a convenient, but accessible, manipulation of equations can be made to solve the problems numerically, is a very engaging strategy that involves actively students in team-based learning and results in an effective understanding of physics.

Acknowledgements

The work of A. G., and P. C., is supported by the Fundação para a Ciência e a Tecnologia (Portuguese Agency for Research) fellowship PD/BD/113536/2015, and Project UID/NAN/50024/2013, respectively.

References

  • [1] Shayer M and Adey P 1983 Towards a Science of Science Teaching: Cognitive development and curriculum demand (London: Heinemann Educational Books) ISBN 0435578251
  • [2] Driver R, Guesne E and Tiberghien A 2000 Childen’s Ideas in Science (Philadelphia: Open University Press) ISBN 0335150403
  • [3] Arons A B 1997 Teaching Introductory Physics (New York: John Wiley & Sons, Inc.) ISBN 9780471137078
  • [4] Serzedello M 2012 Física!? Isso para mim é Grego (São Paulo: Hn Ed. e Publieditorial) ISBN 9788560521388 coleção Produção Académica
  • [5] de la Cotardière P 2010 História das Ciências (Lisboa: Edições Texto & Grafia, Lda.) ISBN 9789898285188 portuguese version of Histoire des Sciences, Tallandier, 2004
  • [6] 2016 Ptolemy URL https://en.wikipedia.org/wiki/Ptolemy
  • [7] 2016 Almagest and references within. URL https://en.wikipedia.org/wiki/Almagest
  • [8] Burge E J 1971 Physics Education 6 003 ISSN 00319120
  • [9] Slegr J and Kraus I 2012 Physics Education 47 309–312 ISSN 0031-9120
  • [10] King W H 1989 Physics Education 24 309 ISSN 00319120
  • [11] Curtis H D 2010 Orbital Mechanics For Engineering Students 2nd ed (Amsterdam: Butterworth-Heinemann) ISBN 9780123747785
  • [12] Jet Propulsion Laboratory 2016 HORIZONS System URL http://ssd.jpl.nasa.gov/?horizons
  • [13] Guerra A G C 2016 Excel Numeric Solution (address to be provided by the journal)
  • [14] Guerra A G C 2016 Orbital Motions (address to be provided by the journal)
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
164913
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description