首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
As part of the NEAR Radio Science investigation, a global solution that includes both spherical and ellipsoidal harmonic gravity fields of Eros, Eros pole and rotation rate, Eros ephemeris, and landmark positions from the optical data was generated. This solution uses the entire one-year in orbit collection of X-band radiometric tracking (Doppler and range) from the Deep Space Network and landmark tracking observations generated from the NEAR spacecraft images of Eros. When compared to a constant density shape model, the gravity field shows a nearly homogeneous Eros. The Eros landmark solutions are in good agreement with the Eros shape model, and they reduce the center-of-mass and center-of-figure offset in the z direction to 13 m. Most of the NEAR spacecraft orbits are determined in all directions to an accuracy of several meters. The solution for the ephemeris of Eros constrains the mass of Vesta to 18.2±0.4 km3/s2 and reduces the uncertainty in the Earth-Moon mass ratio.  相似文献   

2.
Small body surface gravity fields via spherical harmonic expansions   总被引:1,自引:0,他引:1  
Conventional gravity field expressions are derived from Laplace’s equation, the result being the spherical harmonic gravity field. This gravity field is said to be the exterior spherical harmonic gravity field, as its convergence region is outside the Brillouin (i.e., circumscribing) sphere of the body. In contrast, there exists its counterpart called the interior spherical harmonic gravity field for which the convergence region lies within the interior Brillouin sphere that is not the same as the exterior Brillouin sphere. Thus, the exterior spherical harmonic gravity field cannot model the gravitation within the exterior Brillouin sphere except in some special cases, and the interior spherical harmonic gravity field cannot model the gravitation outside the interior Brillouin sphere. In this paper, we will discuss two types of other spherical harmonic gravity fields that bridge the null space of the exterior/interior gravity field expressions by solving Poisson’s equation. These two gravity fields are obtained by assuming the form of Helmholtz’s equation to Poisson’s equation. This method renders the gravitational potentials as functions of spherical Bessel functions and spherical harmonic coefficients. We refer to these gravity fields as the interior/exterior spherical Bessel gravity fields and study their characteristics. The interior spherical Bessel gravity field is investigated in detail for proximity operation purposes around small primitive bodies. Particularly, we apply the theory to asteroids Bennu (formerly 1999 RQ36) and Castalia to quantify its performance around both nearly spheroidal and contact-binary asteroids, respectively. Furthermore, comparisons between the exterior gravity field, interior gravity field, interior spherical Bessel gravity field, and polyhedral gravity field are made and recommendations are given in order to aid planning of proximity operations for future small body missions.  相似文献   

3.
We study the effect of systematic variations in stellar parallaxes over the celestial sphere on the results of a kinematic analysis of stellar proper motions. Our approach is based on the representation of stellar parallaxes by scalar spherical harmonics and on the decomposition of stellar proper motions into a system of vector spherical harmonics. We derive theoretical relations that relate the coefficients of the decomposition of stellar proper motions into toroidal and spheroidal harmonics to the coefficients of the decomposition of stellar parallaxes into scalar spherical harmonics. We have established that the systematic variations of parallaxes over the celestial sphere distort all parameters of the linear Ogorodnikov-Milne model and can be responsible for the appearance of beyond-the-model harmonics. We have performed a kinematic analysis of the proper motions of blue-white and red giants based on Hipparcos data. The parallaxes of blue-white giants show a strong dependence on Galactic latitude (with predominant contraction along the Galactic equator). In contrast, the deviations of the parallaxes from the mean for red giants are localized only in two regions of the celestial sphere. For these samples, the effect of parallax variations over the celestial sphere on kinematic parameters has turned out to be comparable to their rms errors. The global solutions performed using both samples have revealed strong beyond-the-model kinematic effects described by second-order toroidal harmonics and third-order spheroidal harmonics. Using the solutions performed separately in the northern and southern Galactic hemispheres, we have established that not the systematic variations of parallaxes over the celestial sphere but the retardation of Galactic rotation with increasing distance of stars from the principal Galactic plane is mainly responsible for the appearance of these harmonics. Based on these samples of stars, we have estimated the magnitude of the vertical Galactic rotation velocity gradient to be 18.0±2.9 and 22.7±2.2 km s?1 kpc?1, respectively.  相似文献   

4.
A set of spherical harmonics is the most widely used representation of the Earth’s gravity potential. This series converges outside and on the surface of a reference sphere enveloping the Earth. However, the Earth’s surface is better approximated by the reference ellipsoid—a compressed ellipsoid of revolution that covers the entire Earth. The gravity potential can be expanded in a series over ellipsoidal harmonics on the surface of the reference ellipsoid and on the surface of other external confocal ellipsoids of revolution. In contrast to spherical harmonics, depending on the associated Legendre functions of the first kind, ellipsoidal harmonics depend also on the associated Legendre functions of the second kind. The latter contain the very slowly converging hypergeometric Gauss series. The number of series increases with increasing the order of their derivatives. In this work, we derived new series for the gravitational potential of the Earth and its derivatives over ellipsoidal harmonics. Starting from the first order derivative, all the series corresponding to higher order derivatives depend on the same two hypergeometric Gauss series. The latter converges considerably faster than that for the hypergeometric series previously used when computing the gravity potential and its derivatives.  相似文献   

5.
Determination of Shape, Gravity, and Rotational State of Asteroid 433 Eros   总被引:5,自引:0,他引:5  
Prior to the Near Earth Asteroid Rendezvous (NEAR) mission, little was known about Eros except for its orbit, spin rate, and pole orientation, which could be determined from ground-based telescope observations. Radar bounce data provided a rough estimate of the shape of Eros. On December 23, 1998, after an engine misfire, the NEAR-Shoemaker spacecraft flew by Eros on a high-velocity trajectory that provided a brief glimpse of Eros and allowed for an estimate of the asteroid's pole, prime meridian, and mass. This new information, when combined with the ground-based observations, provided good a priori estimates for processing data in the orbit phase.After a one-year delay, NEAR orbit operations began when the spacecraft was successfully inserted into a 320×360 km orbit about Eros on February 14, 2000. Since that time, the NEAR spacecraft was in many different types of orbits where radiometric tracking data, optical images, and NEAR laser rangefinder (NLR) data allowed a determination of the shape, gravity, and rotational state of Eros. The NLR data, collected predominantly from the 50-km orbit, together with landmark tracking from the optical data, have been processed to determine a 24th degree and order shape model. Radiometric tracking data and optical landmark data were used in a separate orbit determination process. As part of this latter process, the spherical harmonic gravity field of Eros was primarily determined from the 10 days in the 35-km orbit. Estimates for the gravity field of Eros were made as high as degree and order 15, but the coefficients are determined relative to their uncertainty only up to degree and order 10. The differences between the measured gravity field and one determined from a constant density shape model are detected relative to their uncertainty only to degree and order 6. The offset between the center of figure and the center of mass is only about 30 m, indicating that Eros has a very uniform density (1% variation) on a large scale (35 km). Variations to degree and order 6 (about 6 km) may be partly explained by the existence of a 100-m, regolith or by small internal density variations. The best estimates for the J2000 right ascension and declination of the pole of Eros are α=11.3692±0.003° and δ=17.2273±0.006°. The rotation rate of Eros is 1639.38922±0.00015°/day, which gives a rotation period of 5.27025547 h. No wobble greater than 0.02° has been detected. Solar gravity gradient torques would introduce a wobble of at most 0.001°.  相似文献   

6.
Tsvetkov  A. S.  Amosov  F. A. 《Astronomy Letters》2020,46(8):509-517
Astronomy Letters - The technique of spherical harmonics, both scalar and vector ones, has long been applied to analyze the astronomical data on a sphere, for example, in the representation of...  相似文献   

7.
This paper presents an analytic solution of the equations of motion of an artificial satellite, obtained using non singular elements for eccentricity. The satellite is under the influence of the gravity field of a central body, expanded in spherical harmonics up to an arbitrary degree and order. We discuss in details the solution we give for the components of the eccentricity vector. For each element, we have divided the Lagrange equations into two parts: the first part is integrated exactly, and the second part is integrated with a perturbation method. The complete solution is the sum of the so-called “main” solution and of the so-called “complementary” solution. To test the accuracy of our method, we compare it to numerical integration and to the method developed in Kaula (Theory of Satellite Geodesy, Blaisdell publ. Co., New York. 1966), expressed in classical orbital elements. For eccentricities which are not very small, the two analytical methods are almost equivalent. For low eccentricities, our method is much more accurate.  相似文献   

8.
We consider the dissipative evolution of a spherical magnetic vortex with a force-free internal structure, located in a resistive medium and held in equilibrium by the potential external field. The magnetic field inside the sphere is force-free (the model of Chandrasekhar in Proc. Natl. Acad. Sci. 42, 1, 1956). Topologically, it is a set of magnetic toroids enclosed in spherical layers. A new exact MHD solution has been derived, describing a slow, uniform, radial compression of a magnetic spheroid under the pressure of an ambient field, when the plasma density and pressure are growing inside it. There is no dissipation in the potential field outside the sphere, but inside the sphere, where the current density can be high enough, the magnetic energy is continuously converted into heat. Joule dissipation lowers the magnetic pressure inside the sphere, which balances the pressure of the ambient field. This results in radial contraction of the magnetic sphere with a speed defined by the conductivity of the plasma and the characteristic spatial scale of the magnetic field inside the sphere. Formally, the sphere shrinks to zero within a finite time interval (magnetic collapse). The time of compression can be relatively small, within a day, even for a sphere with a radius of about 1 Mm, if the magnetic helicity trapped initially in the sphere (which is proportional to the number of magnetic toroids in the sphere) is quite large. The magnetic system is open along its axis of symmetry. On this axis, the magnetic and electric fields are strictly radial and sign-variable along the radius, so the plasma will be ejected along the axis of magnetic sphere outwards in both directions (as jets) at a rate much higher than the diffusive one, and the charged particles will be accelerated unevenly, in spurts, creating quasi-regular X-ray spikes. The applications of the solution to solar flares are discussed.  相似文献   

9.
Spherical harmonics are the natural parameters for the Earth's gravity field as sensed by orbiting satellites, but problems of resolution arise because the spectrum of effects is narrow and unique to each orbit. Comprehensive gravity models now contain many hundreds of thousands of observations from more than thirty different near-Earth artificial satellites. With refinements in tracking systems, newer data is capable of sensing the spherical harmonics of the field experienced by these satellites to very high degree and order. For example, altimeter, laser and satellite-tracking-satellite systems contain gravitational information well above present levels of satellite gravity field recovery (l = 20), but significant aliasing results because the orbital parameters are too restricted compared to the large number of spherical harmonics.It is shown however that the unique spectrum of information for each satellite contained within a comprehensive spherical harmonic model can be represented by simple gravitational constraint equations (lumped harmonics). All such constraints are harmonic in the argument of perigee (ω) with constants determinable directly from tracking data or reconstituted from the comprehensive solution:
(C1, S1) = (Co, So) + Σi = 1 (CCi, SCi) cos i ω + (CSi, SSi) sin i ω
. The constants are simple linear combinations of the geopotential harmonics. Through these lumped harmonics any satellite gravity field can be decomposed and then uniformly extended to any degree or tailored to a given orbit without reintegration of the trajectory and variational equations. They also make possible the inclusion of information into the field from special deep resonance passages, long arc zonal analyses, and satellites unique to other models. Numerous examples of the derivation, combination, extension and tailoring of the harmonics are presented. The importance of using data spanning an apsidal period is emphasized.  相似文献   

10.
Using Hill's variables, an analytical solution of a canonical system of six differential equations describing the motion of a satellite in the gravitational field of the earth is derived. The gravity field, expanded into spherical harmonics, has to be expressed as a function of the Hill variables. The intermediary is chosen to include the main secular terms. The first order solution retains the highly practical formal structure of Kaula's linear solution, but is valid for circular orbits and provides of course a spectral decomposition of radius vector and radial velocity. The resulting eccentricity functions are much simpler than the Hansen functions, since a series evaluation of the Kepler equation is avoided. The present solution may be extended to higher order solutions by Hori's perturbation method.  相似文献   

11.
The gravity field of Venus has been modeled by a spherical harmonic expansion of the potential to degree and order seven. The estimates of these coefficients were obtained by combining information from 43 short arcs (4 hr) of line-of-sight Doppler data centered at periapsis. The data arcs were distributed in longitude and time over more than two circulations of Venus by the Pioneer Venus Orbiter subperiapsis point which was confined to the band of latitudes from 14°N to 17°N. Convergence of the solution has been assured by iterating upon the initial estimate. All estimates were performed with zero a priori information on the gravity coefficients. Since the altitude of periapsis for most of the orbits was within the sensible Venusian atmosphere, drag effects on the estimated harmonics have been removed using an exponential atmosphere density model. Estimates of the mass parameter (GM) of Venus using this dataset are also evaluated.  相似文献   

12.
A method for a kinematic analysis of stellar radial velocities using spherical harmonics is proposed. This approach does not depend on the specific kinematic model and allows both low-frequency and high-frequency kinematic radial velocity components to be analyzed. The possible systematic variations of distances with coordinates on the celestial sphere that, in turn, are modeled by a linear combination of spherical harmonics are taken into account. Theoretical relations showing how the coefficients of the decomposition of distances affect the coefficients of the decomposition of the radial velocities themselves have been derived. It is shown that the larger the mean distance to the sample of stars being analyzed, the greater the shift in the solar apex coordinates, while the shifts in the Oort parameter A are determined mainly by the ratio of the second zonal harmonic coefficient to the mean distance to the stars, i.e., by the degree of flattening of the spatial distribution of stars toward the Galactic plane. The distances to the stars for which radial velocity estimates are available in the CRVAD-2 catalog have been decomposed into spherical harmonics, and the existing variations of distances with coordinates are shown to exert no noticeable influence on both the solar motion components and the estimates of the Oort parameter A, because the stars from this catalog are comparatively close to the Sun (no farther than 500 pc). In addition, a kinematic component that has no explanation in terms of the three-dimensional Ogorodnikov-Milne model is shown to be detected in the stellar radial velocities, as in the case of stellar proper motions.  相似文献   

13.
We have developed a theory of the rotation of the Moon, for the purpose of obtaining libration series explicitly dependent upon lunar gravitational field model parameters. A nonlinear model is used in which the rigid Moon, whose motion in space is that of the main problem of lunar theory, and whose gravity potential is represented through its third degree harmonics, is torqued by the Earth and Sun. The analytical series are then obtained as Poisson series. Numerical comparisons with Eckhardt's solution are briefly exposed.  相似文献   

14.
Theory of the motion of an artificial Earth satellite   总被引:1,自引:0,他引:1  
An improved analytical solution is obtained for the motion of an artificial Earth satellite under the combined influences of gravity and atmospheric drag. The gravitational model includes zonal harmonics throughJ 4, and the atmospheric model assumes a nonrotating spherical power density function. The differential equations are developed through second order under the assumption that the second zonal harmonic and the drag coefficient are both first-order terms, while the remaining zonal harmonics are of second order.Canonical transformations and the method of averaging are used to obtain transformations of variables which significantly simplify the transformed differential equations. A solution for these transformed equations is found; and this solution, in conjunction with the transformations cited above, gives equations for computing the six osculating orbital elements which describe the orbital motion of the satellite. The solution is valid for all eccentricities greater than 0 and less than 0.1 and all inclinations not near 0o or the critical inclination. Approximately ninety percent of the satellites currently in orbit satisfy all these restrictions.  相似文献   

15.
A method for determining the velocity field parameters free from the distortions due to the systematic variations of stellar parallaxes over the celestial sphere is proposed. The method is based on the approximation of parallaxes as a function of coordinates on the sphere using spherical harmonics and can be applied in those cases where the solar motion cannot be eliminated from the stellar proper motions. Numerical experiments have shown that our method is able to obtain accurate coordinates of the solar apex and to calculate the kinematic parameters of the Ogorodnikov-Milne model to within three coefficients of the decomposition of parallaxes into first-order spherical harmonics. Examples of applying the method to the stellar proper motions of the Hipparcos catalogue, which admits checking the results using trigonometric parallaxes, are provided. Such a check has been found to yield a positive result only for nearby stars at heliocentric distances that do not exceed 400 pc and for which the parallaxes were determined with a relative error of at least 30%. An interesting feature of this method is the possibility to construct the shape of the figure which is formed by the deviations of the parallaxes from the sphere corresponding to the average parallaxes of the stars under consideration. It should be specially emphasized that all of this is done in the complete absence of information about the stellar parallaxes. The “solar terms” of the stellar proper motions that are formed by the products of the parallaxes by the solar motion components relative to the centroid of stars are the main source of information about the parallaxes here.  相似文献   

16.
Starting from a complex operator of derivation, we give expressions for derivatives of arbitrary order of the gravity potential with respect to rectangular coordinates. These expressions have a form similar to the original potential expanded in spherical harmonics and are free of singularity at the poles. Computing sets of numerical coefficients once for all, we can compute the derivatives with a very limited work: the same functions are used to compute all derivatives by means of a unique parametrized formula. This is very comfortable for further algebraic manipulations. Numerical tests prove the accuracy and the efficiency of the algorithm derived from our formula to compute the gravity acceleration vector and the gravity gradient tensor.  相似文献   

17.
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 0 = (12,654,875 ± 69) (m2/s2), is used as one of the four fundamental gravity parameters of Mars namely, {W 0, GM, ω, J 20}, 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.  相似文献   

18.
The gravity field dedicated satellite missions like CHAMP, GRACE, and GOCE are supposed to map the Earth's global gravity field with unprecedented accuracy and resolution. New models of the Earth's static and time-variable gravity fields will be available every month as one of the science products from GRACE. A method for the efficient gravity field recovery is presented using in situ satellite-to-satellite observations at altitude and results on static as well as temporal gravity field recovery are shown. Considering the energy relationship between the kinetic energy of the satellite and the gravitational potential, the disturbing potential observations can be computed from the orbital state vector, using high-low GPS tracking data, low–low satellite-to-satellite GRACE measurements, and data from 3-axis accelerometers. The solution method is based on the conjugate gradient iterative approach to efficiently recover the gravity field coefficients and approximate error covariance up to degree and order 120 every month. Based on the monthly GRACE noise-only simulation, the geoid was obtained with an accuracy of a few cm and with a resolution (half wavelength) of 160 km. However, the geoid accuracy can become worse by a factor of 6–7 because of spatial aliasing. The approximate error covariance was found to be a very good accuracy measure of the estimated coefficients, geoid, and gravity anomaly. The temporal gravity field, representing the monthly mean continental water mass redistribution, was recovered in the presence of measurement noise and high frequency temporal variation. The resulting recovered temporal gravity fields have about 0.3 mm errors in terms of geoid height with a resolution of 670 km.  相似文献   

19.
The theoreticl treatment of several geophysical problems presupposes the solution of field equations of the magnetic field in the Earth's mantle. The field equations are given in a scalar form for a spherical model of the Earth. It will be shown that analytical solutions are possible in all cases. The boundary conditions are discussed with regard to the dynamo process in the Earth's core and the existence of a field representation, which is investigated on the Earth's surface. The question is discussed, to what extend the mantle's field is given by this field representation, when some special assumptions about the Earth's model are made.  相似文献   

20.
We solve the problem on a kinematic analysis of the three-dimensional velocity field of stars from zonal catalogues, i.e., catalogues in which the stars are presented at all right ascensions in some declination zones. We have constructed a system of vector spherical harmonics with the properties of completeness and orthogonality for a chosen declination zone. We suggest a method that allows the Ogorodnikov-Milne model parameters in the Galactic coordinate system to be estimated by analyzing the proper motions and radial velocities of stars in the equatorial coordinate system. The vector spherical harmonics are shown to have the following advantages over the standard approach based on a direct leastsquares estimation of the parameters for a specific model. First, in contrast to the standard approach, the new method can reveal all systematic components of the velocity field irrespective of a particular model. Second, it allows one to get rid of the correlation between the sought-for parameters, which presents a serious problem for the conventional method in the case of zonal catalogues. Third, the method of vector spherical harmonics allows the kinematic parameters to be estimated at least by two techniques. Comparison of these two solutions makes it possible to test the standard kinematic model for compatibility with the observational data. The developed method has been tested on the basis of numerical experiments and applied for a kinematic analysis of the proper motions of Tycho-2 stars in the southern hemisphere for which the parallaxes can be estimated using data from the Tycho-2 Spectral Type Catalogue.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号