<?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">IJAA</journal-id><journal-title-group><journal-title>International Journal of Astronomy and Astrophysics</journal-title></journal-title-group><issn pub-type="epub">2161-4717</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ijaa.2015.53019</article-id><article-id pub-id-type="publisher-id">IJAA-59334</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>
 
 
  Use of an Energy-Like Integral to Study the Motion of an Axi-Symmetric Satellite under Drag and Radiation Pressure
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>hmed</surname><given-names>Mostafa</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Department of Mathematics, Faculty of Science, Ain Shams University, Cairo, Egypt</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>ahmed.mahmoud@guc.edu.eg</email></corresp></author-notes><pub-date pub-type="epub"><day>27</day><month>08</month><year>2015</year></pub-date><volume>05</volume><issue>03</issue><fpage>148</fpage><lpage>154</lpage><history><date date-type="received"><day>13</day>	<month>May</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>29</month>	<year>August</year>	</date><date date-type="accepted"><day>1</day>	<month>September</month>	<year>2015</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>
 
 
  The axi-symmetric satellite problem including radiation pressure and drag is treated. The equations of motion of the satellite are derived. An energy-like is given for a general drag force function of the polar angle 
  <em>θ</em>, and then it is used to find a relation for the orbit equation of the satellite with initial conditions satisfying the vanishing of arbitrarily choosing higher derivatives of the velocity.
 
</p></abstract><kwd-group><kwd>Artificial Satellite</kwd><kwd> Drag Effect</kwd><kwd> Radiation Pressure</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>The classical two body problem is one of the most important topics in the field of celestial mechanics, specially the applications of the theory of artificial satellites. Since Brouwer and Hori [<xref ref-type="bibr" rid="scirp.59334-ref1">1</xref>] , so many works have been made to study the problem with different factors considered e.g. Mittleman and Jezewski [<xref ref-type="bibr" rid="scirp.59334-ref2">2</xref>] and Jezewski and Mittleman [<xref ref-type="bibr" rid="scirp.59334-ref3">3</xref>] , Danby [<xref ref-type="bibr" rid="scirp.59334-ref4">4</xref>] , Leach [<xref ref-type="bibr" rid="scirp.59334-ref5">5</xref>] , Gorringe and Leach [<xref ref-type="bibr" rid="scirp.59334-ref6">6</xref>] , McMahon and Scheeres [<xref ref-type="bibr" rid="scirp.59334-ref7">7</xref>] etc.</p><p>Marvaganis [<xref ref-type="bibr" rid="scirp.59334-ref8">8</xref>] studied the motion of an almost constant-speed two body problem under the effect of air resistance. The drag force was taken in the form of Danby’s drag, while Marvaganis and Michalakis [<xref ref-type="bibr" rid="scirp.59334-ref9">9</xref>] studied the two body problem in the existence of Danby’s drag and where the bigger body was radiating. They used a Laplace-like integral to derive the orbit equation. El-Shaboury and Mostafa [<xref ref-type="bibr" rid="scirp.59334-ref10">10</xref>] studied the problem of an axi-symmetric satellite under drag and radiation pressure by first neglecting the effect of axi-symmetry of the satellite, and then adding it as a perturbation to the problem.</p><p>In this work, an attempt is made to get a solution for the problem of an axi-symmetric satellite under drag and radiation pressures, which all the effects are included in the equation of motion from the beginning by using energy like integral. A relation for the orbit equation is derived first for a general air drag function and then for the case of Danby’s drag. Finally, the solution of an almost constant speed satellite has been given.</p></sec><sec id="s2"><title>2. The Equation of Motion and the Integral of Angular Momentum</title><p>The equation of motion of an axi-symmetric satellite under the gravitational force of a spherical body with an additional force due to the resistance force and radiation pressure can be modeled such as Mavraganis and Michalakis (1994), and El-Shaboury and Mostafa (2014).</p><disp-formula id="scirp.59334-formula441"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x5.png"  xlink:type="simple"/></disp-formula><p>The air resistance is taken as a general function R of the polar coordinates<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x6.png" xlink:type="simple"/></inline-formula>, and the definitions of the involved parameters are as following:</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x7.png" xlink:type="simple"/></inline-formula>is the gravitational constant.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x8.png" xlink:type="simple"/></inline-formula>is the radiation constant, where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x9.png" xlink:type="simple"/></inline-formula></p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x10.png" xlink:type="simple"/></inline-formula>are the principal moments of inertia of the satellite (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x11.png" xlink:type="simple"/></inline-formula>about the symmetric axis), and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x12.png" xlink:type="simple"/></inline-formula> is the direction cosine of the radius vector with respect to the axis of the satellite. For simplification, we will rename</p><disp-formula id="scirp.59334-formula442"><graphic  xlink:href="http://html.scirp.org/file/3-4500458x13.png"  xlink:type="simple"/></disp-formula><p>Now, vector product Equation (1) with r, and remembering that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x14.png" xlink:type="simple"/></inline-formula> we get,</p><disp-formula id="scirp.59334-formula443"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x15.png"  xlink:type="simple"/></disp-formula><p>which gives immediately,</p><disp-formula id="scirp.59334-formula444"><graphic  xlink:href="http://html.scirp.org/file/3-4500458x16.png"  xlink:type="simple"/></disp-formula><p>This expression admits a first vector which is the constant direction <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x17.png" xlink:type="simple"/></inline-formula> of the angular momentum<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x18.png" xlink:type="simple"/></inline-formula>. Therefore the motion is planar. This enables us to simplify Equation (2) by writing</p><disp-formula id="scirp.59334-formula445"><label>. (3a)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x19.png"  xlink:type="simple"/></disp-formula><p>Let the resistance be a general function of the polar angel <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x20.png" xlink:type="simple"/></inline-formula> divided by the square of the radial distance r, i.e.</p><disp-formula id="scirp.59334-formula446"><label>(3b)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x21.png"  xlink:type="simple"/></disp-formula><p>We get,</p><disp-formula id="scirp.59334-formula447"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x22.png"  xlink:type="simple"/></disp-formula><p>where, h is the constant angular momentum in the absence of the drag force, and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x23.png" xlink:type="simple"/></inline-formula>.</p></sec><sec id="s3"><title>3. The Energy-Like Invariant of Motion</title><p>Now, let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x24.png" xlink:type="simple"/></inline-formula> be the angle between the radial distance <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x25.png" xlink:type="simple"/></inline-formula> and the distance p of the origin O from the tangent (<xref ref-type="fig" rid="fig1">Figure 1</xref>).</p><p>The vector equation of motion is thus resolved to</p><disp-formula id="scirp.59334-formula448"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x26.png"  xlink:type="simple"/></disp-formula><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Geometrical meaning of the angle φ</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-4500458x27.png"/></fig><disp-formula id="scirp.59334-formula449"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x28.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x29.png" xlink:type="simple"/></inline-formula> is the radius of curvature. From the definition of the angle<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x30.png" xlink:type="simple"/></inline-formula>, we have</p><disp-formula id="scirp.59334-formula450"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x31.png"  xlink:type="simple"/></disp-formula><p>and using the substitution</p><disp-formula id="scirp.59334-formula451"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x32.png"  xlink:type="simple"/></disp-formula><p>Equation (5) becomes,</p><disp-formula id="scirp.59334-formula452"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x33.png"  xlink:type="simple"/></disp-formula><p>This gives,</p><disp-formula id="scirp.59334-formula453"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x34.png"  xlink:type="simple"/></disp-formula><p>where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x35.png" xlink:type="simple"/></inline-formula>. It is clear that the above equation gives the energy integral in the absence of resistance, radiation pressure and oblateness. However in the absence of resistance only, we will still have the invariant of motion</p><disp-formula id="scirp.59334-formula454"><label>(10a)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x36.png"  xlink:type="simple"/></disp-formula></sec><sec id="s4"><title>4. A Relation for the Orbit Equation</title><p>Taking R in the form<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x37.png" xlink:type="simple"/></inline-formula>, where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x38.png" xlink:type="simple"/></inline-formula> is an arbitrary function of the angle<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x39.png" xlink:type="simple"/></inline-formula>. Equation (10) gives,</p><disp-formula id="scirp.59334-formula455"><label>, (11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x40.png"  xlink:type="simple"/></disp-formula><p>which implies by using Equation (4)</p><disp-formula id="scirp.59334-formula456"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x41.png"  xlink:type="simple"/></disp-formula><p>In order to integrate the required integration, we expand V in Taylor series of the polar angle<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x42.png" xlink:type="simple"/></inline-formula>,</p><disp-formula id="scirp.59334-formula457"><label>(13a)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x43.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.59334-formula458"><label>(13b)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x44.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x45.png" xlink:type="simple"/></inline-formula> is the initial value of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x46.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x47.png" xlink:type="simple"/></inline-formula> Then we have</p><disp-formula id="scirp.59334-formula459"><graphic  xlink:href="http://html.scirp.org/file/3-4500458x48.png"  xlink:type="simple"/></disp-formula><p>We write,</p><disp-formula id="scirp.59334-formula460"><graphic  xlink:href="http://html.scirp.org/file/3-4500458x49.png"  xlink:type="simple"/></disp-formula><p>where,</p><disp-formula id="scirp.59334-formula461"><graphic  xlink:href="http://html.scirp.org/file/3-4500458x50.png"  xlink:type="simple"/></disp-formula><p>Thus,</p><disp-formula id="scirp.59334-formula462"><label>, (14a)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x51.png"  xlink:type="simple"/></disp-formula><p>where,</p><disp-formula id="scirp.59334-formula463"><label>(14b)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x52.png"  xlink:type="simple"/></disp-formula><p>Therefore, we have the integration,</p><disp-formula id="scirp.59334-formula464"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x53.png"  xlink:type="simple"/></disp-formula><p>Thus, using Equations (12), (14a), and (15) we get a relation for the orbit equation in the form,</p><disp-formula id="scirp.59334-formula465"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x54.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x55.png" xlink:type="simple"/></inline-formula> is given by Equation (14b), <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x56.png" xlink:type="simple"/></inline-formula>is given by Equation (13b) and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x56.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x57.png" xlink:type="simple"/></inline-formula></p><p>Equation (16) describes a relation for the orbit equation of an axi-symmetric satellite with oblateness coeffi-</p><p>cient k under radiation pressure of coefficient <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x58.png" xlink:type="simple"/></inline-formula> and air drag whose function is given by<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x59.png" xlink:type="simple"/></inline-formula>, where</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x60.png" xlink:type="simple"/></inline-formula>can be chosen arbitrary.</p><p>The convergence of the involved series is guaranteed for initial velocity satisfying the vanishing of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x61.png" xlink:type="simple"/></inline-formula></p><p>for all n &gt; N, where N can be chosen arbitrary.</p><sec id="s4_1"><title>4.1. The Case of Danby Drag</title><p>In the special case of Danby’s drag (Dabny, 1962), <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x62.png" xlink:type="simple"/></inline-formula>, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x63.png" xlink:type="simple"/></inline-formula> where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x64.png" xlink:type="simple"/></inline-formula> is a constant, the required integrations reduce to:</p><disp-formula id="scirp.59334-formula466"><label>(17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x65.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.59334-formula467"><label>(18)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x66.png"  xlink:type="simple"/></disp-formula><p>Substituting,</p><disp-formula id="scirp.59334-formula468"><label>(19)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x67.png"  xlink:type="simple"/></disp-formula><p>we get the integration (18) in the form<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x68.png" xlink:type="simple"/></inline-formula>. To evaluate this integral, we distinguish between two cases</p><p>for n when it is even or odd,</p><p>When n is even, we write</p><disp-formula id="scirp.59334-formula469"><graphic  xlink:href="http://html.scirp.org/file/3-4500458x69.png"  xlink:type="simple"/></disp-formula><p>and then we use the expansion</p><disp-formula id="scirp.59334-formula470"><graphic  xlink:href="http://html.scirp.org/file/3-4500458x70.png"  xlink:type="simple"/></disp-formula><p>thus we get</p><disp-formula id="scirp.59334-formula471"><graphic  xlink:href="http://html.scirp.org/file/3-4500458x71.png"  xlink:type="simple"/></disp-formula><p>And when n is odd, we write</p><disp-formula id="scirp.59334-formula472"><graphic  xlink:href="http://html.scirp.org/file/3-4500458x72.png"  xlink:type="simple"/></disp-formula><p>then we use the expansion</p><disp-formula id="scirp.59334-formula473"><graphic  xlink:href="http://html.scirp.org/file/3-4500458x73.png"  xlink:type="simple"/></disp-formula><p>thus we get</p><disp-formula id="scirp.59334-formula474"><graphic  xlink:href="http://html.scirp.org/file/3-4500458x74.png"  xlink:type="simple"/></disp-formula><p>We can collect the two cases together in one case to get,</p><disp-formula id="scirp.59334-formula475"><label>(20)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x75.png"  xlink:type="simple"/></disp-formula><p>Substituting from Equations (17 - 20) into Equation (16), we get after simplification</p><disp-formula id="scirp.59334-formula476"><label>(21)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x76.png"  xlink:type="simple"/></disp-formula><p>Equation (21) gives a relation for the orbit equation of an axi-symmetric satellite under the gravitational effect of a radiating body and air resistance described by Danby’s drag.</p></sec><sec id="s4_2"><title>4.2. The Case of an Almost Constant Speed Satellite</title><p>If the satellite is of almost constant speed, then we assume that the first derivative is of small value, and all the higher derivatives to be zero (e.g. Mavraganis, 1991), we get the solution</p><disp-formula id="scirp.59334-formula477"><label>(22)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-4500458x77.png"  xlink:type="simple"/></disp-formula><p>Equation (22) is a special case of Equation (21) when the satellite is of almost constant speed.</p></sec></sec><sec id="s5"><title>5. Conclusions</title><p>In this paper, the motion of an axi-symmetric satellite under the effect of a radiating body in the presence of air drag is studied. An energy-like integral for the problem has been evaluated using a Taylor expansion for the velocity around the initial value of the polar angel. The convergence of the integral is guaranteed by the assump-</p><p>tion that the derivative <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-4500458x78.png" xlink:type="simple"/></inline-formula> for all n &gt; N, for an arbitrary N.</p><p>The energy-like integral has been used to get a relation for the orbit equation of the satellite. The relation is derived first for a general air drag function and then for the case of Danby’s drag. Finally, the solution of an almost constant speed satellite has been given.</p></sec><sec id="s6"><title>Cite this paper</title><p>AhmedMostafa, (2015) Use of an Energy-Like Integral to Study the Motion of an Axi-Symmetric Satellite under Drag and Radiation Pressure. International Journal of Astronomy and Astrophysics,05,148-154. doi: 10.4236/ijaa.2015.53019</p></sec></body><back><ref-list><title>References</title><ref id="scirp.59334-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Brouwer, D. and Hori, G. (1961) Theoretical Evaluation of Atmospheric Drag Effects in the Motion of an Artificial Satellite. The Astronomical Journal, 66, 193-225. http://dx.doi.org/10.1086/108399</mixed-citation></ref><ref id="scirp.59334-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Mittleman, D. and Jezwski, D. (1982) An Analytic Solution to the Classical Two-Body Problem with Drag. Celestial Mechanics and Dynamical Astronomy, 28, 401-413. http://dx.doi.org/10.1007/BF01372122</mixed-citation></ref><ref id="scirp.59334-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Jezwski, D. and Mittleman, D. (1983) Integrals of Motion for the Classical Two-Body Problem with Drag. International Journal of Non-Linear Mechanics, 18, 119-124. http://dx.doi.org/10.1016/0020-7462(83)90039-2</mixed-citation></ref><ref id="scirp.59334-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Danby, G.M.A. (1962) Fundamentals of Celestial Mechanics. MacMillan, New York.</mixed-citation></ref><ref id="scirp.59334-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Leach, P.G.L. (1987) The First Integrals and Orbit Equation for the Kepler Problem with Drag. Journal of Physics A: Mathematical and General, 20, 1997-2004. http://dx.doi.org/10.1088/0305-4470/20/8/019</mixed-citation></ref><ref id="scirp.59334-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Gorringe, V.M. and Leach, P.G.L. (1988) Hamiltonlike Vectors for a Class of Kepler Problems with Drag. Celestial Mechanics and Dynamical Astronomy, 41, 125-130. http://dx.doi.org/10.1007/BF01238757</mixed-citation></ref><ref id="scirp.59334-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">McMahon, J. and Scheeres, D. (2010) Secular Orbit Variation due to solar Radiation Effects: A Detailed Model for BYORP. Celestial Mechanics and Dynamical Astronomy, 106, 261-300. http://dx.doi.org/10.1007/s10569-009-9247-9</mixed-citation></ref><ref id="scirp.59334-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Mavraganis, A.G. (1991) The Almost Constant-Speed Two-Body Problem with Resistance. Celestial Mechanics and Dynamical Astronomy, 51, 395-405. http://dx.doi.org/10.1007/BF00052930</mixed-citation></ref><ref id="scirp.59334-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Mavraganis, A.G. and Michalakis, D.G. (1994) The Two-Body Problem with Drag and Radiation Pressure. Celestial Mechanics and Dynamical Astronomy, 58, 393-403. http://dx.doi.org/10.1007/BF00692013</mixed-citation></ref><ref id="scirp.59334-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">El-Shaboury, S.M. and Mostafa, A. (2014) The Motion of Axisymmetric Satellite with Drag and Radiation Pressure. Astrophysics and Space Science, 352, 515-519. http://dx.doi.org/10.1007/s10509-014-1975-y</mixed-citation></ref></ref-list></back></article>