Use of an Energy-Like Integral to Study the Motion of an Axi-Symmetric Satellite under Drag and Radiation Pressure

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 θ, 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.


Introduction
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 [1], so many works have been made to study the problem with different factors considered e.g.Mittleman and Jezewski [2] and Jezewski and Mittleman [3], Danby [4], Leach [5], Gorringe and Leach [6], McMahon and Scheeres [7] etc.
Marvaganis [8] 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 [9] 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 [10] 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.
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.

The Equation of Motion and the Integral of Angular Momentum
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). ( ) The air resistance is taken as a general function R of the polar coordinates , r θ , and the definitions of the in- volved parameters are as following: µ is the gravitational constant.β is the radiation constant, where 0 1.

β < < ,
C A are the principal moments of inertia of the satellite ( C about the symmetric axis), and γ is the direc- tion cosine of the radius vector with respect to the axis of the satellite.For simplification, we will rename ( )( ) Now, vector product Equation (1) with r, and remembering that , = × H r r  we get, ( ) which gives immediately, 0.

× = H H 
This expression admits a first vector which is the constant direction H = H e H of the angular momentum H . Therefore the motion is planar.This enables us to simplify Equation (2) by writing ( ) .

The Energy-Like Invariant of Motion
Now, let φ be the angle between the radial distance r and the distance p of the origin O from the tangent (Figure 1).The vector equation of motion is thus resolved to ( ) ( ) where ρ is the radius of curvature.From the definition of the angle φ , we have and using the substitution Equation ( 5) becomes, ( ) This ( ) where ( ) 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

A Relation for the Orbit Equation
Taking R in the form ( ) , where ( ) f θ is an arbitrary function of the angle θ .Equation (10) gives, which implies by using Equation ( 4) A. Mostafa 151 In order to integrate the required integration, we expand V in Taylor series of the polar angle θ , ( ) where 0 θ is the initial value of θ and ( ) where, where, Therefore, we have the integration, ( Thus, using Equations ( 12), (14a), and (15) we get a relation for the orbit equation in the form, where n c is given by Equation (14b), n a is given by Equation (13b) and

( ) ( )
where α is a constant, the required integrations reduce to: we get the integration (18) in the form d n y y y a − + ∫ .To evaluate this integral, we distinguish between two cases for n when it is even or odd, When n is even, we write And when n is odd, we write

∑ ∫
We can collect the two cases together in one case to get, ( ) ( ) ( ) Substituting from Equations (17 -20) into Equation ( 16), we get after simplification 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.

The Case of an Almost Constant Speed Satellite
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 ( ) ( ) ( ) Equation ( 22) is a special case of Equation (21) when the satellite is of almost constant speed.

Conclusions
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 assumption that the derivative 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.

=
for all n > N, for an arbitrary N.