<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">AM</journal-id><journal-title-group><journal-title>Applied Mathematics</journal-title></journal-title-group><issn pub-type="epub">2152-7385</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/am.2012.39142</article-id><article-id pub-id-type="publisher-id">AM-22991</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Restricted Three-Body Problem in Different Coordinate Systems
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>ohammed</surname><given-names>A. Sharaf</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Aisha</surname><given-names>A. Alshaery</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Astronomy, Faculty of Science, King Abdulaziz University, Jeddah, KSA</addr-line></aff><aff id="aff2"><addr-line>Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah, KSA</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>sharaf_adel@hotmail.com(OAS)</email>;<email>alshaary@hotmail.com(AAA)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>27</day><month>09</month><year>2012</year></pub-date><volume>03</volume><issue>09</issue><fpage>949</fpage><lpage>953</lpage><history><date date-type="received"><day>July</day>	<month>10,</month>	<year>2012</year></date><date date-type="rev-recd"><day>August</day>	<month>10,</month>	<year>2012</year>	</date><date date-type="accepted"><day>August</day>	<month>17,</month>	<year>2012</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  In this paper of the series, the equations of motion for the spatial circular restricted three-body problem in sidereal spherical coordinates system were established. Initial value procedure that can be used to compute both the spherical and Cartesian sidereal coordinates and velocities was also developed. The application of the procedure was illustrated by numerical example and graphical representations of the variations of the two sidereal coordinate systems.
 
</p></abstract><kwd-group><kwd>Spatial Restricted Circular Three Body Problem; Regularization; Coordinate Transformations</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>In a previous communication to this journal [<xref ref-type="bibr" rid="scirp.22991-ref1">1</xref>], hereafter referred to as Paper I we started our studies towards establishing new differential equations for the different forms of the three-body problem using some important coordinate systems. By this, we aims at obtaining differential equations (see [<xref ref-type="bibr" rid="scirp.22991-ref1">1</xref>] for details) which are: 1) Regular; 2) Suitable for the geometry to which they referred; 3) Producing slow variations in the coordinates during the orbital motion, a property which produces more stable numerical integration procedures. In Paper I, the equations of motion for spatial restricted circular three body problem in cylindrical coordinates system was established together with a computational algorithm that can be used to compute both the cylindrical and Cartesian coordinates and velocities. In the present paper, the equations of motion for the spatial circular restricted threebody problem in sidereal spherical coordinates system were established. Initial value procedure that can be used to compute both the spherical and Cartesian sidereal coordinates and velocities was also developed. The application of the procedure was illustrated by numerical example and graphical representations of the variations of the two sidereal coordinate systems.</p></sec><sec id="s2"><title>2. Circular Restricted Three-Body Problem in Sidereal System</title><p>If two of the bodies, say m<sub>1</sub> and m<sub>2</sub> in the three-body problem move in circular, coplanar orbits about their common center of mass and the mass say m<sub>3</sub> of the third body is too small to affect the motion of the other bodies, the problem of the motion of the third body is called the circular, restricted, three body problem. The two revolving bodies are called the primaries; their masses are arbitrary but have such internal mass distributions that they may be considered point masses.</p><p>The equations of motion of the third body in a dimensionless sidereal (inertial) coordinate <img src="1-7400967\5a2212fe-da2c-4867-b497-4761da64ceb7.jpg" /> system with the mean motion<img src="1-7400967\8f864d5b-8a6e-4d08-a731-bf337b278b22.jpg" />, are [<xref ref-type="bibr" rid="scirp.22991-ref2">2</xref>]</p><disp-formula id="scirp.22991-formula10122"><label>(1)</label><graphic position="anchor" xlink:href="1-7400967\1f2e8e01-a367-4fc4-876b-e4654d6344a4.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.22991-formula10123"><label>(2)</label><graphic position="anchor" xlink:href="1-7400967\2175c788-beab-401a-a3c8-40696c331b76.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.22991-formula10124"><label>(3)</label><graphic position="anchor" xlink:href="1-7400967\039d0e50-9c6a-420d-a163-54ae93a48336.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="1-7400967\7ca71f8b-b0aa-4f5a-9a8c-73d7f900044b.jpg" /> is given as</p><disp-formula id="scirp.22991-formula10125"><label>(4)</label><graphic position="anchor" xlink:href="1-7400967\551a5276-c305-4a15-9831-2f9a83e2ef40.jpg"  xlink:type="simple"/></disp-formula><p><img src="1-7400967\fcd14212-702a-467d-8d04-37f437d6cde0.jpg" />denotes the mass of the smaller primary when the total mass of the primaries has been normalized to unity.</p><disp-formula id="scirp.22991-formula10126"><label>(5)</label><graphic position="anchor" xlink:href="1-7400967\35ad997e-97d9-40e5-83be-9d7790e3b630.jpg"  xlink:type="simple"/></disp-formula><p>and <img src="1-7400967\98b4dec5-6fa3-4dfa-b3d6-372a3180ad79.jpg" /> are the distances of the third body from the primaries which are located at<img src="1-7400967\d787a7fe-2da9-46cb-8c90-860a788847f8.jpg" />;<img src="1-7400967\97ed7291-6a90-4052-b67f-7eb0ed53862d.jpg" />, these coordinates are functions of the time t and are given as</p><disp-formula id="scirp.22991-formula10127"><label>(6)</label><graphic position="anchor" xlink:href="1-7400967\2a2bce55-23d9-4b66-867b-b265e2e9ec2a.jpg"  xlink:type="simple"/></disp-formula></sec><sec id="s3"><title>3. The Equations of Motion in Sidereal Spherical Coordinate System</title><p>Corresponding to the Cartesian sidereal coordinate system<img src="1-7400967\5f6ecc73-0684-4d81-9380-3119a1062450.jpg" />, the coordinate system related to the system <img src="1-7400967\cbed7ac6-0605-4796-b6b2-f366eae2b727.jpg" /> by certain transformation, is also called sidereal coordinate system. In this respect the system <img src="1-7400967\56353c21-46b4-4ec0-9841-27350f7afd3d.jpg" /> of Equation (7) is called sidereal spherical coordinate system.</p><p>In what follows we shall establish, the differential equations for the spatial circular restricted three bodyproblem in sidereal spherical coordinate system.</p><sec id="s3_1"><title>3.1. Coordinate and Velocity Transformations</title><disp-formula id="scirp.22991-formula10128"><label>(7)</label><graphic position="anchor" xlink:href="1-7400967\e932a00c-3ceb-4cb4-ac17-a2cdd3f3b42d.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.22991-formula10129"><label>(8.1)</label><graphic position="anchor" xlink:href="1-7400967\670e9d10-d90a-4100-801d-bfc4d6ca02fa.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.22991-formula10130"><label>(8.2)</label><graphic position="anchor" xlink:href="1-7400967\cf33c41b-7b01-45f2-8542-538536805d7f.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.22991-formula10131"><label>(8.3)</label><graphic position="anchor" xlink:href="1-7400967\7a5ad54b-be7e-4930-abf9-b10dfb49299a.jpg"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.22991-formula10132"><label>(9)</label><graphic position="anchor" xlink:href="1-7400967\7b41b380-5059-4847-8a88-7d144ab06f49.jpg"  xlink:type="simple"/></disp-formula></sec><sec id="s3_2"><title>3.2. Inverse Transformations</title><p>From Equation (7) we have</p><disp-formula id="scirp.22991-formula10133"><label>(10)</label><graphic position="anchor" xlink:href="1-7400967\3aca0a63-24f0-4017-b187-a522f5de04cd.jpg"  xlink:type="simple"/></disp-formula><p>Differentiating the first and the third of Equation (10) and the third of Equation (7) with respect to the time t we get:</p><disp-formula id="scirp.22991-formula10134"><label>(11)</label><graphic position="anchor" xlink:href="1-7400967\9163d0f9-26fd-4767-af99-52f242fab7a4.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="1-7400967\f21465c4-1513-4035-a6a1-95482f208e80.jpg" /> and <img src="1-7400967\56778987-89f8-4a16-9ea9-c1fdf1ca8279.jpg" /> are given in terms of <img src="1-7400967\d0fc57ac-9392-4be0-b319-ab4ee2d7328a.jpg" /> and <img src="1-7400967\35efac9f-621e-4acf-9a49-537adc029871.jpg" /> from the previous equations.</p></sec><sec id="s3_3"><title>3.3. The Equations of Motion</title><p>The kinetic energy of a particle of unit mass in the spherical coordinate system is</p><disp-formula id="scirp.22991-formula10135"><label>(12)</label><graphic position="anchor" xlink:href="1-7400967\ecae0876-4a70-403e-aa9f-86b5fcda6143.jpg"  xlink:type="simple"/></disp-formula><p>By using the transformation equations (Equations (7)), the gravitational potential V could be expressed in term of<img src="1-7400967\0c31da57-eafe-4ae7-bb82-887a604b4488.jpg" />.</p><p><img src="1-7400967\9c4c247e-2e53-4eb2-81ab-61e31457aad4.jpg" /></p><p>Consequently, we deduce for the equations of motion in sidereal spherical coordinate system, the forms</p><disp-formula id="scirp.22991-formula10136"><label>(13.1)</label><graphic position="anchor" xlink:href="1-7400967\08683e3d-d9c7-4699-8487-2840c9fd5231.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.22991-formula10137"><label>(13.2)</label><graphic position="anchor" xlink:href="1-7400967\bcae9d67-bc8f-49b3-b035-66c29e921e8c.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.22991-formula10138"><label>(13.3)</label><graphic position="anchor" xlink:href="1-7400967\0882c2f0-7dd7-494d-96f7-6fc7e109a722.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="1-7400967\d6173ad8-a440-472d-aa88-5d55d24c58c6.jpg" /> are given as</p><disp-formula id="scirp.22991-formula10139"><label>(14)</label><graphic position="anchor" xlink:href="1-7400967\58207f37-0a29-4bd5-8c4d-115b032a907d.jpg"  xlink:type="simple"/></disp-formula><p><img src="1-7400967\78b1c2db-70cd-41d8-becd-078d889734d9.jpg" />, <img src="1-7400967\80ba201f-a665-4129-86a5-ea47565c5b47.jpg" />and<img src="1-7400967\ed39902c-9a37-49b8-abbc-ac2ec9dd41e4.jpg" />; <img src="1-7400967\2f36dc37-4766-4df2-bcd2-5fefbd796b31.jpg" />can be computed from Equation (7), while<img src="1-7400967\19fc8b96-9c6d-4b9b-bd0e-7f3622746a5b.jpg" />, <img src="1-7400967\9fae4467-e30c-4f1a-919b-95d72d678798.jpg" />and <img src="1-7400967\f96eba3f-008a-4181-b631-9bc436f306f1.jpg" /> can be computed from Equations (1)-(3), so we get</p><p><img src="1-7400967\96b2b9d4-32a6-4be6-b6ac-f2f1a2c2ae9d.jpg" /></p><p><img src="1-7400967\4043e05d-f210-482e-a03e-818413faef59.jpg" /></p><p><img src="1-7400967\b79d70ef-c361-4fc0-ab33-6443cb10b91c.jpg" />where</p><p><img src="1-7400967\effa4767-15d8-45a5-a926-622a15207e9b.jpg" /></p><p><img src="1-7400967\c9b91155-2eaa-469c-a1dc-9dfa08d4c379.jpg" /></p><p><img src="1-7400967\bc656b5b-6fb1-4826-804b-5f1c98fcc011.jpg" />,</p><p><img src="1-7400967\722fe881-4287-48c4-a294-e063d7fa730b.jpg" />,</p><p><img src="1-7400967\9842456a-73f6-4ff8-87c1-1c9c7c0d08c7.jpg" />,</p><p><img src="1-7400967\cb17b845-78c7-4a0c-a5c6-b47b50859ee0.jpg" /></p><p><img src="1-7400967\d34c03e1-4061-426d-98a7-58f8c8b19a7c.jpg" />.</p></sec></sec><sec id="s4"><title>4. Computational Development</title><sec id="s4_1"><title>4.1. Initial Value Procedure</title><p>In what follows, we shall establish a procedure that can be used to compute <img src="1-7400967\69621b9f-9a1f-49d5-a833-80693265ae28.jpg" /> (say) both:</p><p>1) The spherical sidereal coordinates and velocities<img src="1-7400967\500f899a-7333-49da-aaac-9201d596ab37.jpg" />, and 2) The Cartesian sidereal coordinates and velocities<img src="1-7400967\574465c3-bc60-46fa-a765-961ccffe5c04.jpg" />.</p><p>So, such procedure is a double usefulness computational algorithm, for which a differential solver can be used for the spherical sidereal six order system to obtain<img src="1-7400967\96c00656-3397-40f2-8976-6df6aa920a97.jpg" />. While the Cartesian sidereal coordinates and velocities <img src="1-7400967\d4d25817-8f6c-40b2-8d3c-a341ab0a0bba.jpg" /> are obtained by the substitutions in the direct transformation formulae (Equations (7) and (8)), rather than solving the six order system of Equations (1), (2) and (3). By this way, great time can be saved. &#160;</p><p>This initial value procedure using sidereal spherical coordinate system will be described through its basic points, input, output and computational steps.</p><p>Input: 1) <img src="1-7400967\9a9a05af-80c4-49a4-a3b3-a2b5aa2e1045.jpg" />at<img src="1-7400967\215c9858-981b-4715-a86a-ca1577bf04f7.jpg" />2) the final time <img src="1-7400967\79255ee4-c633-450e-ab54-6ca5395c363b.jpg" /></p><p>3) <img src="1-7400967\50760d79-4519-444a-81d5-aa0cd23fa7c2.jpg" /></p><p><img src="1-7400967\fb741df5-cd43-4210-95a9-621fd51ca606.jpg" /></p><p>Output: 1) <img src="1-7400967\f3fa378a-a9be-4c7f-b638-1438aa6668f8.jpg" /><img src="1-7400967\2897dc5d-3ff0-4cf1-9831-711e3e266d19.jpg" /></p><p>2) <img src="1-7400967\0675bcdd-999b-401e-b254-82a16e2bd2ef.jpg" /><img src="1-7400967\aecf03ad-0191-4507-94d3-edd9cd35b4f3.jpg" /></p><sec id="s4_1_1"><title>Computational steps</title><p>1) Using the given values <img src="1-7400967\bd7a3d02-7edc-4168-89ed-50154cb9f9a0.jpg" /> at <img src="1-7400967\89a02221-a1b2-42b6-8132-28e35bc0290d.jpg" /> and the inverse transformations to compute the initial values<img src="1-7400967\abf1e4df-f193-43d2-b8c8-64b3e7633061.jpg" />.</p><p>2) Using the partial derivatives <img src="1-7400967\abaf9e7a-1216-4d5c-be36-f48bccca92b0.jpg" /> (functions of<img src="1-7400967\e331b353-bddd-49e4-adfb-c10a4f3b946e.jpg" />) to construct the analytical forms of equations of motion as first order system.</p><table-wrap-group id="1"><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> The values of sidereal spherical coordinates and velocities</title></caption></table-wrap-group><table-wrap-group id="2"><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> The values of sidereal Cartesian coordinates and velocities</title></caption></table-wrap-group><p>3) Using the initial conditions<img src="1-7400967\3b7d3ad5-f47c-4a8e-947b-57432293fda5.jpg" /> from step 1 to solve numerically the above differential system of step 2 for <img src="1-7400967\3bcc9d9d-d582-40ce-aa1a-8088ac6c3f34.jpg" /> <img src="1-7400967\c174ce77-d818-4be4-b989-387fbe5cd8a9.jpg" />, (note that<img src="1-7400967\2d62633a-5da4-4c12-a63e-7b3f5630d0dc.jpg" />).</p><p>4) Using <img src="1-7400967\34b118a6-61f3-479b-adb5-68595c538209.jpg" /> from step 3 and the direct transformations of Equations (7) and (8) to compute numerically <img src="1-7400967\0a0cf51f-2dca-4dea-9d63-18bfa602db82.jpg" /> and <img src="1-7400967\ef0f0f6d-d857-4f54-864e-8811f3210888.jpg" /> <img src="1-7400967\84c16c7f-5339-498f-82a9-dabb4a348cef.jpg" />.</p><p>5) End.</p></sec></sec><sec id="s4_2"><title>4.2. Numerical Example</title><p>Consider the initial values</p><p><img src="1-7400967\f623debe-dda2-4f32-99de-156490454aee.jpg" /></p><p><img src="1-7400967\0f7776e2-25b2-480a-967d-26d19f772587.jpg" /></p><p><img src="1-7400967\8a2b8adf-0b12-426c-868a-38a66dc663ee.jpg" /></p><p><img src="1-7400967\6bf3ad7b-0479-4790-82d7-a1640716c5e7.jpg" /></p><p><img src="1-7400967\19d131ed-fade-40f2-80ff-c3cb68b66625.jpg" /></p><p><img src="1-7400967\55beb659-9942-4e99-9431-d45a3149e3fc.jpg" /></p><p><img src="1-7400967\35e33c02-e10d-4701-bfbc-966d14071d38.jpg" /></p><p><img src="1-7400967\e1f9b2c8-c9fa-4beb-9a46-e4c601a0f9b3.jpg" /></p><p>applying the above procedure we get the results as displayed in Tables 1 and 2.</p></sec><sec id="s4_3"><title>4.3. Graphical Representations</title><p><xref ref-type="fig" rid="fig1">Figure 1</xref> illustrates the time variations of the two sidereal coordinate systems <img src="1-7400967\63940872-ab60-4b4e-bd30-fc9651b4fc23.jpg" /> (left) and <img src="1-7400967\0af1a507-5443-41a2-87c2-6efa6c4efbc2.jpg" /> (right).</p></sec></sec><sec id="s5"><title>5. Conclusion</title><p>In this paper of the series, the equations of motion for the spatial circular restricted three-body problem in sidereal spherical coordinates system were established. Initial value procedure that can be used to compute both the spherical and Cartesian sidereal coordinates and velocities was also developed. The application of the procedure</p><p>was illustrated by numerical example and graphical representations of the variations of the two sidereal coordinate systems.</p></sec><sec id="s6"><title>REFERENCES</title></sec><sec id="s7"><title>NOTES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.22991-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">M. A. Sharaf and A. A. Alshaery, “Restricted Three- Body Problem in Cylindrical Coordinates System, (Paper I),” Journal of Applied Mathematics, 2012, in press.</mixed-citation></ref><ref id="scirp.22991-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">V. Szebehely, “Theory of Orbits,” Academic Press, New York, 1967.</mixed-citation></ref></ref-list></back></article>