Frozen orbits for satellites close to an Earth-like planet

We say that a planet is Earth-like if the coefficient of the second order zonal harmonic dominates all other coefficients in the gravity field. This paper concerns the zonal problem for satellites around an Earth-like planet, all other perturbations excluded. The potential contains all zonal coefficientsJ ₂ throughJ ₉. The model problem is averaged over the mean anomaly by a Lie transformation to the second order; we produce the resulting Hamiltonian as a Fourier series in the argument of perigee whose coefficients are algebraic functions of the eccentricity — not truncated power series. We then proceed to a global exploration of the equilibria in the averaged problem. These singularities which aerospace engineers know by the name of frozen orbits are located by solving the equilibria equations in two ways, (1) analytically in the neighborhood of either the zero eccentricity or the critical inclination, and (2) numerically by a Newton-Raphson iteration applied to an approximate position read from the color map of the phase flow. The analytical solutions we supply in full to assist space engineers in designing survey missions. We pay special attention to the manner in which additional zonal coefficients affect the evolution of bifurcations we had traced earlier in the main problem (J ₂ only). In particular, we examine the manner in which the odd zonalJ ₃ breaks the discrete symmetry inherent to the even zonal problem. In the even case, we find that Vinti's problem (J ₄+J ₂ ² =0) presents a degeneracy in the form of non-isolated equilibria; we surmise that the degeneracy is a reflection of the fact that Vinti's problem is separable. By numerical continuation we have discovered three families of frozen orbits in the full zonal problem under consideration; (1) a family of stable equilibria starting from the equatorial plane and tending to the critical inclination; (2) an unstable family arising from the bifurcation at the critical inclination; (3) a stable family also arising from that bifurcation and terminating with a polar orbit. Except in the neighborhood of the critical inclination, orbits in the stable families have very small eccentricities, and are thus well suited for survey missions.     

A New Pumping Mechanism for A Series of Class II Methanol Masers in the J 0 - J -1 E Transitions

A new pumping mechanism – methanol masers without population inversion is presented in this paper. It can be used to explain the formation of a series of J ₀ → J _-1 E methanol masers, while the 2₁ → 3₀ A ⁺ methanol masers are regarded as a driving coherent micrwave field. In the new mechanism, the intensities of J ₀ - J _-1 E methanol masers are increased with the decreasing transition frequencies (or with rotational number J, approximately). These results agree with Slysh et al. (1995) and Slysh et al. (1999) J ≤ 5 observations for G3345.01+1.79 and W48, in which both J ₀ → J _-1 E and 2₁ → 3₀ A ⁺ methanol masers are detected coincidentally. Other astronomical conditions, such as magnetic field, 2₁ → 3₀ A ⁺ coherent radition, incoherent pumping rate by thermal radition and so on are also discussed. The new mechanism can operate as a complement to other ordinary maser pumping mechanisms for some class II methonal maser sources. This revised version was published online in July 2006 with corrections to the Cover Date.     

Analytical Approach to the Motion of a Lunar Artificial Satellite

The canonical equations of motion of an artificial lunar satellite are formulated including the effects of the asphericity of the Moon comprising the harmonics J ₂, J ₂₂, J ₃, J ₃₁, J ₄ andJ ₅, the oblateness of the Earth up to the second zonal harmonic, as well as the disturbing function due to the attractions of the Earth and of the Sun (terms are retained up to order 10^-6 for the higher orbits and 10^-8 for the lower orbits). This revised version was published online in July 2006 with corrections to the Cover Date.     

Methanol 72–81A+ emission and some maser series in Orion KL

Methanol 7₂–8₁ A ⁺ is mapped for the first time in Orion KL. Analysing the observed data and solving the statistical equilibrium and radiative transfer equations, it is concluded that line series ofJ ₂–(J+1)₁ A ⁺ (J=7,8,9) is in quasi-thermal emission rather than the masers in Orion KL. The maser spots of methanolJ ₂–J ₁ E (J=6,7) and 8₀–7₁ A ⁺ are distributed in the northeast part of the contour plot of 7₂–8₁ A ⁺. The physical conditions of the regions of maser seriesJ ₀–(J–1)₁ A ⁺ (J=7,8,9) are discussed. Also from the calculation results another maser seriesJ ₁–(J–1)₂ A ^– (J=10,11,12) that might coexist with maser seriesJ ₂–J ₁ E, is found. The sizes of the 2-dimension Gaussian fit plots of methanol 7₂–8₁ A ⁺ and HCOOCH₃ 10(0,10)–9(0,9)A are almost the same, and the main parts overlap each other.     

Jupiter’s moment of inertia: A possible determination by Juno

The moment of inertia of a giant planet reveals important information about the planet’s internal density structure and this information is not identical to that contained in the gravitational moments. The forthcoming Juno mission to Jupiter might determine Jupiter’s normalized moment of inertia NMoI = C/MR² by measuring Jupiter’s pole precession and the Lense–Thirring acceleration of the spacecraft (C is the axial moment of inertia, and M and R are Jupiter’s mass and mean radius, respectively). We investigate the possible range of NMoI values for Jupiter based on its measured gravitational field using a simple core/envelope model of the planet assuming that J₂ and J₄ are perfectly known and are equal to their measured values. The model suggests that for fixed values of J₂ and J₄ a range of NMoI values between 0.2629 and 0.2645 can be found. The Radau–Darwin relation gives a NMoI value that is larger than the model values by less than 1%. A low NMoI of ∼0.236, inferred from a dynamical model (Ward, W.R., Canup, R.M. [2006]. Astrophys. J. 640, L91–L94) is inconsistent with this range, but the range is model dependent. Although we conclude that the NMoI is tightly constrained by the gravity coefficients, a measurement of Jupiter’s NMoI to a few tenths of percent by Juno could provide an important constraint on Jupiter’s internal structure. We carry out a simplified assessment of the error involved in Juno’s possible determination of Jupiter’s NMoI.     

Construction d'orbites periodiques de satellites dans un systeme d'axes tournants,I

Résumé On développe une méthode de construction d'orbites périoldiques dans un système d'axes tournants, pour un satellite gravitant autour d'un sphéroide. Les orbites sont quasi circulaires,i est l'inclinaison sur le plan équatorial de la planète. Pour les petites inclinaisons, la solution est donnée jusqu'aux termes enJ ₂ ² etJ ₄.Ce modèle peut être appliqué aux satellites de Saturne. Des valeurs observées des longitudes des noeuds ascendants de Mimas et Téthys, on donne une estimation des valeurs deJ ₂ etJ ₄ du potentiel de Saturne. La valeur deJ ₂ est très sensible aux valeurs adoptées pour le rayon équatorial de la planète.

---

Construction of periodic orbits of satellites in a moving system of axes, I

We give an algorithm for the construction of periodic orbits in a rotating frame for the cases of satellites moving around an oblate planet.The orbits are near to the circular case; the asymptotic developments of the periodic solutions are completely calculated for the termsJ ₂ andJ ₄ of the potential. The solutions for small inclinations are given up toJ ₂ ² .The families of solutions depend on three parameters: the semi-major axis, the inclination of the generating orbit and the initial position on this orbit.These solutions can be applied to the motion of the Saturnian satellites. From the observed longitudes of the ascending nodes of Mimas and Tethys, we estimate the valuesJ ₂ andJ ₄ of the Saturnian potential, the value ofJ ₂ very strongly depends on the adopted value of the planet's equatorial diameter.

     

On the theory of the oblateness of the Sun

In this paper we first emphasize why it is important to know the successive zonal harmonics of the Sun's figure with high accuracy: mainly fundamental astrometry, helioseismology, planetary motions and relativistic effects. Then we briefly comment why the Sun appears oblate, going back to primitive definitions in order to underline some discrepancies in theories and to emphasize again the relevant hypotheses. We propose a new theoretical approach entirely based on an expansion in terms of Legendre's functions, including the differential rotation of the Sun at the surface. This permits linking the two first spherical harmonic coefficients (J ₂ and J ₄) with the geometric parameters that can be measured on the Sun (equatorial and polar radii). We emphasize the difficulties in inferring gravitational oblateness from visual measurements of the geometric oblateness, and more generally a dynamical flattening. Results are given for different observed rotational laws. It is shown that the surface oblateness is surely upper bounded by 11 milliarcsecond. As a consequence of the observed surface and sub-surface differential rotation laws, we deduce a measure of the two first gravitational harmonics, the quadrupole and the octopole moment of the Sun: J ₂=−(6.13±2.52)×10⁻⁷ if all observed data are taken into account, and respectively, J ₂=−(6.84±3.75)×10⁻⁷ if only sunspot data are considered, and J ₂=−(3.49±1.86)×10⁻⁷ in the case of helioseismic data alone. The value deduced from all available data for the octopole is: J ₄=(2.8±2.1)×10⁻¹². These values are compared to some others found in the literature. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1005238718479     

A New Reference Equipotential Surface, and Reference Ellipsoid for the Planet Mars

Using the shape model of Mars GTM090AA in terms of spherical harmonics complete to degree and order 90 and gravitational field model of Mars GGM2BC80 in terms of spherical harmonics complete to degree and order 80, both from Mars Global Surveyor (MGS) mission, the geometry (shape) and gravity potential value of reference equipotential surface of Mars (Areoid) are computed based on a constrained optimization problem. In this paper, the Areoid is defined as a reference equipotential surface, which best fits to the shape of Mars in least squares sense. The estimated gravity potential value of the Areoid from this study, i.e. W ₀ = (12,654,875 ± 69) (m²/s²), is used as one of the four fundamental gravity parameters of Mars namely, {W ₀, GM, ω, J ₂₀}, i.e. {Areoid’s gravity potential, gravitational constant of Mars, angular velocity of Mars, second zonal spherical harmonic of gravitational field expansion of Mars}, to compute a bi-axial reference ellipsoid of Somigliana-Pizzetti type as the hydrostatic approximate figure of Mars. The estimated values of semi-major and semi-minor axis of the computed reference ellipsoid of Mars are (3,395,428 ± 19) (m), and (3,377,678 ± 19) (m), respectively. Finally the computed Areoid is presented with respect to the computed reference ellipsoid.     

Cs studies of dense cores in regions of high mass star formation: a survey of southern molecular masers

We searched for the CSJ = 2 – 1 emission towards 29 southern H₂O and H₂O/OH masers and 1 OH maser with the SEST radio telescope. We detected and mapped 24 CS emitting regions probably associated with 27 H₂O masers. The C³⁴SJ = 2 – 1 and COJ = 1 – 0 lines were also observed at the grid positions closest to the CS peaks. Four cores were mapped in the CSJ = 5 – 4 and C³⁴SJ = 2 – 1 lines.     

A new outlook on the `Differential Theory' of the solar Quadrupole Moment and Oblateness

The modeling of the quadrupole moment J ₂ and of the oblateness , two key solar parameters, derives from the development in successive spherical harmonics of the gravitational potential. These harmonics are representative of the shape of the Sun, by studying the local distortion of the internal layers, under their distribution of mass and velocity. The first aim of this paper is to study, over the radius r and the colatitude , the structure of the internal layers of the Sun through a geometrical approach, considering J ₂ and under a differential form. The second aim is to determine their theoretical values, after integration over r and , taking the best available models of density and rotation into account constrained by helioseismic data. The novelty of our approach lies in the analysis of the profiles of the two above-mentioned solar parameters, under differential form, from the core to the surface. This analysis allows us to comply with the physical processes located in the transition regions, namely the tachocline and maybe a new subsurface layer which could be called the leptocline. The profiles of tildeJ ₂ show two sharp decreases, which are directly connected to the shear layers located at 0.7 R and beneath the surface. The profiles of tilde show five changes of curvature, which seem to be connected to solar processes, such as the matter circulation flows, seismic events or the storage of the magnetic field, phenomena taking place in the transition regions. These sets of profiles allow us to propose a configuration scenario composed of a double layer. Moreover, as a result of the theoretical determination of tildeJ ₂ and tilde, the values at the surface of the quadrupole moment and of the oblateness can be deduced, which are 1.60×10^–7 and 8.77×10^–6, respectively. As a result of an analysis of available data, we may admit J ₂=(2.0±0.4)×10^–7. The theoretical computations of J ₂ and at the surface will be confronted in the near future with the values measured in space by means of the Picard microsatellite. This mission should permit one to measure at the same time both the total solar irradiance and the latitudinal diameters in any position angle (after removing the passing spots or faculae at the limb).     

The Earth-figure perturbations in the Lunar theory

We compute the perturbations on the motion of the Moon due to the shape of the Earth. The zonal terms inJ ₂,J ₃, andJ ₄ are considered. The accuracy is estimated at 3×10^–5 and the results compared with previous theories.     

Long-term motion of resonant satellites with arbitrary eccentricity and inclination

A first-order, semi-analytical method for the long-term motion of resonant satellites is introduced. The method provides long-term solutions, valid for nearly all eccentricities and inclinations, and for all commensurability ratios. The method allows the inclusion of all zonal and tesseral harmonics of a nonspherical planet.We present here an application of the method to a synchronous satellite includingonly theJ ₂ andJ ₂₂ harmonics. Global, long-term solutions for this problem are given for arbitrary values of eccentricity, argument of perigee and inclination.     

Analytical theory of earth's artificial satellites (A.T.E.A.S.)

The problem of A.T.E.A.S. is treated, for the zonal perturbations, in its Hamiltonian form. The method consists in eliminating angular variables from the Hamiltonian function. Nearly identity canonical transformations are used, first to remove short periodic terms, second to remove long periodic terms. The general solution, up toJ ₂ ³ , is represented by the generators of the transformations and by the mean motions of averaged variables, known up toJ ₂ ⁴ . Open expressions in the eccentricity are avoided as far as possible. It permits to obtain a closed second order theory with closed third order mean motions.Proceedings of the Sixth Conference on Mathematical Methods in Celestial Mechanics held at Oberwolfach (West Germany) from 14 to 19 August, 1978.     

A noncanonical analytic solution to theJ 2 perturbed two-body problem

The motion of a satellite subject to an inverse-square gravitational force of attraction and a perturbation due to the Earth's oblateness as theJ ₂ term is analyzed, and a uniform, analytic solution correct to first-order inJ ₂, is obtained using a noncanonical approach. The basis for the solution is the transformation and uncoupling of the differential equations for the model. The resulting solution is expressed in terms of elementary functions of the independent variable (the ‘true anomaly’), and is of a compact and simple form. Numerical results are comparable to existing solutions.     

GPS-Based On-Board Orbit Determination of a Satellite Using Extended H ∞ Filtering Algorithms

In this paper the extended H filtering algorithms for the design of the GPS-based on-board autonomous navigation system for a low earth orbit (LEO) satellite are introduced. The dynamic process models for the estimation of position, velocity and acceleration from the GPS measurements are established. The nominal orbit of the small LEO satellite is determined by using the 7th–8th order Runge—Kutta algorithms. Three filtering approaches are applied to smooth the orbit solutions, respectively, based upon the simulated GPS pseudo range observables using the Satellite Navigation Tool Box. The simulation shows that the observed orbit errors obtained by using the extended H filtering algorithms can be reduced to a lower level than the observed orbit errors in the sense of RMS within 12 h of tracking time by using the H filtering algorithms and the extended Kalman filtering algorithms under the appropriately designed parameters. Based upon the position errors predicted by the three filtering algorithms after the last observation, we find that the extended H filtering algorithm provides the least position errors of the user satellite.This revised version was published online in October 2005 with corrections to the Cover Date.     

A note on the algorithm of symplectic integrators

For a Hamiltonian that can be separated into N+1(N\geq 2) integrable parts, four algorithms can be built for a symplectic integrator. This research compares these algorithms for the first and second order integrators. We found that they have similar local truncation errors represented by error Hamiltonian but rather different numerical stability. When the computation of the main part of the Hamiltonian, H ₀, is not expensive, we recommend to use S ^* type algorithm, which cuts the calculation of the H ₀ system into several small time steps as Malhotra(1991) did. As to the order of the N+1 parts in one step calculation, we found that from the large to small would get a slower error accumulation. This revised version was published online in July 2006 with corrections to the Cover Date.     

A new estimate of the quadrupole moment of the sun

Recent solar observations at Pic du Midi are reported that yield a value of J ₂=(2.57 ± 2.36) x 10^–6 for the quadrupole moment of the Sun. These observations were conducted from July 1993 to July 1994 after several improvements of the scanning heliometer. This instrument operates by fast photoelectric scans of opposite limbs of the Sun quasi-simultaneously, which provides the distance between both inflection points of the limb profiles. Any number of solar diameters in any position angle can be measured within a time interval short enough to minimize the scattering of the observational parameters. Errors due to atmospheric deterioration are discussed. From our results, compared to previous values obtained by other authors, it can be concluded than an upper limit for J ₂ is probably 1.0 × 10^-5.     

Unsteady MHD flow past a porous plate with step function change in suction at slip flow regime

Unsteady two-dimensional hydromagnetic flow of an electrically conducting viscous incompressible fluid past a semi-infinite porous flat plate with step function change in suction velocity is studied allowing a first order velocity slip at the boundary condition. The solution of the problem is obtained in closed form and the results are discussed with the aid of graphs for various parameters entering in the problem.Notations B intensity of magnetic field - H magnetic field parameter,H=(M+1/4)^1/2–1/2 - h rarefaction parameter - L ₁ slip coefficient; ;I, mean free path of gas molecules;f, Maxwell's reflection coefficient - M magnetic field parameter - r suction parameter - t time - t dimensionless time - u velocity of the fluid - u dimensionless velocity of the fluid - U velocity of the fluid at infinity - v suction velocity - v ₁ suction velocity att<=0 - v ₂ suction velocity att>0 - x distance parallel to the plate - y distance normal to the plate - y nondimensional distance normal to the plate - v kinematic viscosity - electric conductivity of the fluid - density of the fluid - shear stress at the wall - nondimensional shear stress at the wall - erf error function - erfc complementary error function     

Regularization and linearization of the equations of motion in central force-fields

Within the framework of linear and regular celestial mechanics, a wide class of central force field problems are considered. We take as potential function a polynomial whose variable is the reciprocal of the distance from the origin, and, as regularizing function the square root of a certain reciprocal polynomial, related to the potential function. The relations between the coefficients of both functions, in order to arrive to linear equations, are given. In particular, the case of a quintic polynomial is analyzed, and an application is made to artificial earth satellites considering harmonicsJ ₂,J ₃,J ₄.     

Candidate red clump giants in the PPMX catalogue

This paper presents the RCGP catalogue of more than 0.5 million candidate red clump stars with the limiting magnitude K _s = 9.5 ^m . These stars are selected from the PPMX catalogue as the most probable red clump members by analyzing the color-reduced proper motion diagrams built from the proper motions given in PPMX and J, K _s -photometry given in the 2MASS catalogue. Reddening of the selected stars is used to find extinction in the K _s -band and to consider it in the further analysis. The two-dimensional galactic rotation model generalized by Ogorodnikov is used to investigate the tangential velocity field of the selected red clump members, most of which are thin disk stars located within 1.5 kpc from the sun. The values of kinematic parameters and solar components are determined as a function of stellar heights above the galactic equatorial plane and their heliocentric distances.