- +1-530-264-8006
- info@tutorsglobe.com

18,76,764

Questions

Asked

21,311

Experts

9,67,568

Questions

Answered

Start Excelling in your courses, Ask an Expert and get answers for your homework and assignments!!

Submit Assignment2015 © Tutors Globe. All rights reserved.

## Satellite Orbits in the Two-Body System, Physics tutorial

Concept of Two-Body Problem:In practice, there are several investigations which involve more than two bodies. For example, let a hypothetical satellite placed in orbit between the earth and mars. The earth, the sun and mars apply forces on satellite. Therefore, four bodies (Earth, Sun, Mars and satellite itself) are engaged. This is no longer a two-body problem. It is a four-body problem. Luckily, it is frequently allowable to split many-body problem in numerous two-body problems. This technique makes it feasible for approximate values of parameters involved to be computed with relative ease. For example, four-body problem given above can be split in 3 different two-body problems comprising of (a) two body problem involving the earth and the satellite only (b) two-body problem involving sun and the satellite only (c) a two-body problem involving mars and the satellite only.

When satellite is under main attraction of the earth, the two-body problem involving the earth and the satellite becomes about valid. Similar statement holds for each of the remaining two-body problems. Though, in carry-over region where earth and sun attract the satellite with the similar order of magnitude of force, three-body treatment becomes essential. Similar assertion also holds in carry-over region where mars and sun attract the satellite with about the same force. Luckily, these carry-over regions constitute only small fraction of total journey of satellite in its orbit. As a result, much calculation time is saved using this technique.

Solution to Two-Body Problem:Let two bodies of masses m1 and m2 separated by the linear displacement r

^{→}. Newton's second law of motion defines thatF

^{→}= d(mv^{→})/dt or F^{→}= md^{2}r^{→}/dt^{2 }or md^{2}r^{→}/dt^{2}= ∑i=1kF^{→}iForce of attraction, F

_{1}^{→}, on m_{1}is directed along vector r^{→}. Force of attraction, F_{2}, on m_{2}is in opposite direction.By Newton's third law,

F

_{1}^{→ }= -F_{2}^{→}Also, by Newton's law of gravitation,

F

_{1}^{→ }= Gm_{1}m_{2}/r_{2}(r^{→}/r)Again

F

_{2}^{→}= -G(m_{1}m_{2}/r_{2})(r^{→}/r)Let vectors r

_{1}^{→}and r_{2}^{→ }be directed from some fixed reference point O to particles of masses m_{1}and m_{2}, respectively.Now,

F

^{→}= m(d^{2}r1^{→}/dt^{2})so that

F

_{1}^{→}= m_{1}(d^{2}r^{→}1/dt^{2}) and F_{2}^{→}= m_{2}(d^{2}r^{→}2/dt^{2})Equation (1) can then be rewritten as

m

_{1}(d^{2}r_{1}^{→}/dt^{2}) = (Gm_{1}m_{2}/r_{2})(r^{→}/r)Equation (2) can also be rewritten as

m

_{2}(d^{2}r_{2}^{→}/dt^{2}) = (Gm_{1}m_{2}/r_{2})(r^{→}/r)By Adding, we get

m

_{1}(d^{2}r_{2}^{→}/dt^{2}) + m_{2}((d^{2}r_{2}^{→}/dt^{2})) = 0by definition, the centre-of-mass C of the system of masses m

_{1}and m_{2}is given by MR^{→}= m_{1}r_{1}^{→}+ m_{2}r_{2}^{→}where M = m_{1}+ m_{2 }and R^{→}is the position vector of C Hence, by Eq.(8), you getM(dR

^{→}/dt) = a^{→}Also by equation we getMR

^{→ }= a^{→}t + b^{→}These two relations show that centre of mass of the system moves with constant velocity (i.e., with constant speed in straight line).

Angular Momentum Integral:Equations may be written as

d

^{2}r^{→}/dt^{2}= Gm_{2}(r^{→}/r_{3}) and d^{2}r_{2}^{→}/dt^{2}= -Gm_{1}(r^{→}/r_{3}) respectivelyAfter solving equation gives

d

^{2}(r^{→}_{1}- r_{2}^{→})/dt^{2}= G(m_{1}+ m_{2})(r^{→}/r_{3}) But from figurer

_{1}^{→}- r_{2}= -r^{→}Hence d^{2}r^{→}/dt^{2}+ μr^{→}/r_{3}= 0where μ = G(m

_{1}+ m_{2})Taking the vector product of r→ with Equation you obtain

r

^{→}xd^{2}r^{→}/dt^{2}+ μr^{→}/r_{3}x r^{→}= 0i.e.

Integrating, we have r

^{→ }x dr^{→}/dt = h^{→ }where h^{→}is constant vector.Equation is angular momentum integral. Now, as h

^{→}is constant, it should point in one direction or all values of t. Therefore, motion of one body about other should lie in plane defined by th direction of h.Mathematical form of Kepler's Second Law:If polar coordinates r and θ are taken in this plane velocity component along radius vector joining m

_{1}to m_{2}is r^. Also, perpendicular component of radius vector is rθ^, where dot replaces d/dt.By defining I

^{→}and J^{→}as the unit vectors along and perpendicular to the radius vector, the resultant velocity r→^ becomesr

^{→}^ = I^{→}r+ J^{.}^{→}rθ^{.}Therefore, using Equations we can easily verify each of the given steps:

r

^{→ }x dr^{→}/dt = h^{→}r

^{→ }x r^{→}^ = h^{→}r

^{→ }x(I^{→}r+ J^{.}^{→}rθ) = h^{.}^{→}I

^{→}rx(I^{→}r+ J^{.}^{→}rθ) = h^{.}^{→}(I

^{→}rxI^{→})r^{.}+(I^{→}r x J^{→})rθ= h^{.}^{→}0 + (K

^{→}r)rθ= h^{.}^{→}Thus K

^{→}r^{2}θ= h^{.}^{→}We may thus write:

r

^{2}θ= h^{.}Here constant h is twice the rate at which radius vector defines the area. This is mathematical form of Kepler's second law. For the circular orbit, rate of change of area may, of course, be evaluated as follows:

A = πr

^{2}dA/dt = dA/dr x dr/dt

2πrdr/dt => 2πrr

^{.}Energy Conservation Equation of System:If scalar product of r

^{→.}with equation is now taken, we getr+ μ/r

^{→.}.d^{2}r^{→}/dt^{2}+ μr^{→.}.r^{→}/r^{3}= 0 i.e. r^{→.}.r^{→..}^{3}r^{→.}.r^{→}= 0This may be integrated to provide

1/2r

^{→.}.r^{→.}-μ/r = cHere C is a constant v.

I.e., 1/2v

^{2 }- μ/r = cHere v is velocity.

Equation is energy conservation equation of the system. Quantity C isn't the total energy of system. Though, quantity 1/2v

^{2}is associated to kinetic energy whereas quantity -μ/r is associated to potential energy of the system.Satellite Obits in Two-Body System:The components of acceleration along and perpendicular to radius vector are

r

- rθ^{..}^{.}^{2}and 1/2(d/dt)(r^{2}θ) respectively. Equation can be written as:^{.}I

^{→}(r- rθ^{..}^{.}^{2}) + J^{→}[(1/r)(d/dt)(r^{2}θ)] + (μ/r^{.}^{3})(I^{→}r) = 0Equating coefficients of vectors, we get:

r

- rθ^{..}^{.}^{2}= -μ/r^{2}and (1/r)(d/dt)(rθ^{.}^{2}) = 0Integration of equation provides angular momentum integral

r

^{2}θ= h^{.}Substituting u = 1/r and eliminating time t between equations we get:

d

^{2}u/dt^{2}+ u = μ/h^{2}General solution of equation is as follows:

u = μ/h

^{2}+ Acos(θ-ω)Here A and ω are two constants of integration. Reintroducing r in equation it becomes:

1/r = μ/h

^{2}+ Acos(θ-ω) orr = [1/(μ/h

^{2}+ Acos(θ-ω))] orr = [(h

^{2}/μ)/1-(Ah^{2}/μ)cos(θ-ω)]But, generally, polar equation of the conic section may be written as:

r = P/(1 + ecos(θ-ω)) so that p = h

^{2}/μ and e = Ah^{2}/μTherefore, solution of two-body problem is a conic section. It comprises Kepler's first law as special case. Particularly, orbit of a satellite about earth is categorized by value of eccentricity e. Four cases to be remembered are:

i) For 0 < e < 1, orbit is ellipse.

ii) For e = 1, orbit is parabola.

iii) For e >1, orbit is hyperbola.

iv) For e = 0, orbit is circle

Tutorsglobe: A way to secure high grade in your curriculum (Online Tutoring)Expand your confidence, grow study skills and improve your grades.

Since 2009, Tutorsglobe has proactively helped millions of students to get better grades in school, college or university and score well in competitive tests with live, one-on-one online tutoring.

Using an advanced developed tutoring system providing little or no wait time, the students are connected on-demand with a tutor at www.tutorsglobe.com. Students work one-on-one, in real-time with a tutor, communicating and studying using a virtual whiteboard technology. Scientific and mathematical notation, symbols, geometric figures, graphing and freehand drawing can be rendered quickly and easily in the advanced whiteboard.

Free to know our price and packages for online physics tutoring. Chat with us or submit request at info@tutorsglobe.com