budynek aiut AIUT sp. z o. o.
★  HenrykDot.com   ★
is the online companion to a series of books published by AIUT under the common main title
"Physics of My Imaginary Space-Time" by Henryk Dot.
POLSKI

last update
01/02/2023
Physics 3 - Return
Return to Physics 3

Physics 4 -
Home
Introduction

Physics 4 - Chapter 1
Motion
Orbital motion

Contact email: henryk.dot(at)aiut.com
"Subject" should begin with
a digit corresponding to the day
of the week., e.g. Sunday=7




Books published by AIUT
are found in libraries according to the list of compulsory copies.

Second Edition of "Fizyka 3"
ISBN 978-83-926856-1-6
Fizyka
can be bought in Warsaw
in the Academic Bookstore
PW Publishing House
Noakowskiego street 18/20

and in Katowice
in the bookstore "Liber"
Bankowa street 11.
(area of Silesian University)

English edition of "Physics"
ISBN 978-83-926856-2-3
Fizyka
is also in libraries
and the distribution method should be asked wydawca@aiut.com.

2.  Orbital Motion

     Considerations on orbital motion will come useful for further sections of the book.
The motion of a single object will be considered. In a real world, there are multiple objects interacting with each other. Discussion of a single object is sufficient to present the essence of the issue. Perhaps many graduate students will find this presented approach as preferable way to deal with this topic.
We will consider motion in which the acceleration is described by the relation: 

1x18 (1.18)
where:
           r    - position vector,
           r - position vector length (modulus),
       1x18a - functional relation: acceleration dependence on distance.
In our considerations, b will take constant values, as for example in the case of gravity  b = kg M  or for case of atom  b = Z⋅ e2 / (4⋅π⋅me⋅ ε0)
where:
           kg - gravitational constant,
           M - mass of the central object,
           Z - number of protons in a nucleus,
           e - electron charge,
           me - electron mass,
           ε0 - vaccum permittivity.
Many of consideration/conclusions are still correct when b is not a constant but
a function of  r  i.e.,  b = b(r).

For the sake of clarity, we begin our considerations with an assumption that the central object is stationary and its interaction is immediate. This has no significant impact on
the validity of the conclusions drawn at this stage.

In some problems, a consideration of the motion of the central object and the delay of its interaction will be very important, but we will deal with these in another section.

Orbital motion is conveniently analyzed in polar coordinates    r  and  α  such as:
    x = r ⋅ cos(α). , (1.19)
    y = r ⋅ sin(α) . (1.20)
We obtain then the following relations
1x21, (1.21)
1x22, (1.22)
1x23, (1.23)
1x24, (1.24)
Inserting (1.23) and (1.24) respectively into (1.18) we obtain:
1x25, (1.25)
and
1x26. (1.26)
Now, performing multiplications (1.25) and (1.26) by cos(α) and sin(α) respectively,
then adding, and then subtracting both equations by sides and finally replacing
the expression 1x26a with new variable 1x26b we obtain very clear and simple formulas (1.27) and (1.28).
2w2 (1.27)
2w2 (1.28)
We can also rewrite formula (1.27) by expressing ω by 1x27a , where vs - circumferential velocity or in other words tangential velocity to a circle with radius r
1x29. (1.29)
It should be noted that formula (1.28) does not use any variables besides r and ω and which is very important for farther considerations it does not contain the value b !
We just obtained a complete set of orbital formulas (1.27) and (1.28). Very often instead of these two, only one formula is presented applicable to circular orbit:
1x30 (1.30)
It represents the specific case of formula (1.29), when 1x30b .
However, formula (1.30) still can be useful for us.
When we rewrite it in the form as below,
1x31 (1.31)
it can be noticed that (using interesting relation from formula (1.17)) an object in orbital motion on circular orbit contains kinetic energy equal to a half of the energy required to leave a field.
The velocity momentum is constant not only for orbital motion. If an object moves
in a straight line with a constant speed, the constant velocity momentum is conserved
in respect to any point lying outside the straight line of motion.
The principle of conservation of angular momentum is directly related to the principle of conservation of velocity momentum.

Formula (1.28) can be written as:
1x32 (1.32)
which results in
1x33 (1.33)
and after substitution of ω by 1x27a
1x34. (1.34)
It is the same as
1x35 (1.35)
thus
1x36 (1.36)
and this allows us to express our desired rule without using a mass m.
r ⋅ vs = k = const (1.37)
Obviously, since  r ⋅ vs = k = const thus m ⋅ vs ⋅ r  must be a constant value as well.
Formula (1.37) is equivalent to its three-dimensional vector representation:
r    × v    = k    (1.38)
Relation shown in (1.37) can be used as constant in the formula (1.29), which
then be expressed as:
1x39 (1.39)
If k and b are constant, the expression 1x39a is also constant and it determines certain radius, which we can denote as 1x39b and then express formula (1.39) as:
1x40 (1.40)

This formula is very interesting, and together with (1.37) encompasses a different approach in presenting the property known as the Laplace-Runge-Lenz vector.


Let us try to find some interesting features of equation (1.40).
Let vr be denoted as a radial component of the velocity
1x41 (1.41)
We can then use the following relation
1x42 (1.42)
and express equation (1.40) as:
1x43 (1.43)
In this equation, the variables are conveniently separated, so we can integrate on both sides of the equation, yielding:
1x44 (1.44)
where C is an integration constant.
We can also eliminate the ½ coefficient on the left side of the equation:
1x45 (1.45)
To determine integration constant C we will take advantage of the fact that at certain extreme points r, radial speed equals zero. We can write the equation for these points:
1x46 (1.46)
which leads to a quadratic equation:
0 = - b ⋅ rk + 2⋅b ⋅ r + C ⋅ r2 (1.47)
The equation has two solutions:
1x48a (1.48a)
and
1x48b (1.48b)
Denoting:
1x49 (1.49)
formula(1.48) can be written as:
1x50 (1.50)
and
1x50a
We easily notice that
1x51 (1.51)
and
1x52. (1.52)
In addition (1.51) represents the sum of the distances from any point on the perimeter
to both foci of the ellipse which are used to plot the ellipse.
Thus equation (1.45) can be expressed as:
1x53 (1.53)
and then including (1.51) and (1.52)
1x54 (1.54)
The same result can be yielded when using energy formulas.
1x55 (1.55)
where:
     E – total energy,
     v – total speed.
For orbital motion, the total energy of such system is negative.
Calculations will be easier to conduct when we relate the energy to a unit of mass
and, we denote it with the index m i.e.
1x56 (1.56)
By substituting total speed by its components, we obtain:
1x57 (1.57)
which is equivalent to formula (1.45) best visible after equation expressions are rearranged
1x58 (1.58)
and it is the same equation since:  k2 = b ⋅ rk   and   2 ⋅ Em = C .
This prompted a suggestion to associate radius  r  with energy and assign the index E to it.

When deriving the formula for vr relation (1.55) was not used to avoid the use of the terms work and energy, which we had not defined before. They are referenced only to demonstrate full compatibility of mathematical derivations.


The new values introduced as radii  rk  and  rE   seem to better reflect the parameters
of the ellipse than the popularly used term– eccentricity.

Then the equation of an ellipse using   rE   and    rk   values:
1x59 (1.59)
and in parametric form, which is easy to plot:
1x60a (1.60a)
1x60b (1.60b)
It is useful to remember the characteristic values:
longer semi-axis of an ellipse     rE   , (1.61)
shorter semi-axis                     1x61a, (1.61a)
amplitude of a radius changes 1x61b . (1.61b)
Amplitude of a radius changes corresponds to the distance of the foci of the ellipse from the center of the ellipse and is directly related to commonly used concept - eccentricity.

The characteristic quantities can be supplemented with the time of the ellipse circulation, which in our case represents the period of radius change i.e. the transition time from
the maximum radius to the minimum radius and back to the maximum radius.
To obtain this value, it is sufficient to calculate half of this period that is, the transition time from the minimum to the maximum radius.
We will use the definition of radial velocity
1x62. (1.62)
Thus,
1x63a and 1x63b (1.63)
Substituting vr with formula (1.54) we obtain (1.54)
1x64 (1.64)
then, after solving the integral
1x64a
this equals to:
1x84b
and results with the answer
1x65. (1.65)
From formula (1.65) we derive Kepler’s Third Law
1x66. (1.66)

The description of orbital motion proposed here using two parameters,  rE  i  rk
characterizes the properties of orbital motion very well.
Keeping either parameter  rE  or  rk   constant, we obtain two families of orbits shown in the figures below, one with constant energy and the other with constant velocity momentum (Fig. 2).
The  rE  radius corresponds to the energy and also determines the radius against which
the value varies with amplitude 1x61b.

The radius  rk determines the radius with respect to which the direction of the force changes along the radius.. For  r > rk  here is an attraction, and for  r < rk   a repulsion.

For an ellipse family of constant energy, i.e., constant  rE , two extreme cases are interesting. First when rE = rk  and the orbit is circular and the second when rk  approaches zero and the orbit approaches a straight-line segment (shown as green line in the picture). The latter case is practically impossible, but still interesting.

To illustrate the effect of fixed  rE  and  rk   in both figures, the origin of the coordinate system is placed at the first focus.
Fig. 1 shows a family of orbits with equal energy – fixed rE .
and rk varies from zero to 0.9 * rE   in  0.1 * rE intervals .
x y Sorry, your browser does not support inline SVG.
Fig.1.    

Fig. 2 shows a family of orbits with equal velocity momentum – fixed rk .
rE varies from zero to 4.0 * rk  in  0.5 * rk intervals .
x y Sorry, your browser does not support inline SVG.
Fig.2.    






  © 2020 Henryk Dot -