共查询到20条相似文献，搜索用时 46 毫秒
1.
G. H. Born E. J. Christensen L. K. Seversike 《Celestial Mechanics and Dynamical Astronomy》1974,9(1):41-53
The concept of employing osculating reference position and velocity vectors in the numerical integration of the equations of motion of a satellite is examined. The choice of the reference point is shown to have a significant effect upon numerical efficiency and the class of trajectories described by the differential equations of motion. For example, when the position and velocity vectors on the osculating orbit at a fixed reference time are chosen, a universal formulation is yielded. For elliptical orbits, however, this formulation is unattractive for numerical integration purposes due to Poisson terms (mixed secular) appearing in the equations of motion. Other choices for the reference point eliminate this problem but usually at the expense of universality. A number of these formulations, including a universal one, are considered here. Comparisons of the numerical characteristics of these techniques with those of the Encke method are presented. 相似文献
2.
The methods for analytical determination of partial derivatives of the current parameters of motion with respect to their initial values are described. The methods take into account principal perturbations and are based on the use of the osculating and superosculating intermediate orbits constructed earlier by the author. These orbits ensure the first-, second-, and third-order contact to the real trajectory at the initial time. The solution for parameters of the intermediate motion and partial derivatives of these parameters is given in a universal closed form. The partial derivatives on long time intervals are computed using a step-by-step procedure combined with the Encke method of special perturbations, in which the intermediate orbits are used as the reference. The numerical results show that the new approach can be efficiently used for solving the problem of differential correction of orbits of asteroids and comets on the basis of observational data. 相似文献
3.
Numerical methods are usually used for the computation of ephemerides with perturbations for the precise orbital determination of an artificial satellite. But their numerical stability will be encountered in a long arc. In this case the use the improved Encke special perturbation methods has been suggested. The results of this paper show that Encke's method does indeed have a certain effectiveness, but cannot yet completely resolve the numerical stability, and the more efficient method is to use the energy integral or its variational relation to control the growth of the along-track error in general numerical calculations so that the aim of stabilization can be achieved. 相似文献
4.
A small generalization of the equation of motion for the Poynting-Robertson effect is tested in order to find the significance of new terms. The test is made for dust particles ejected at perihelia of the orbit of the comet Encke. The particles are released at the speed of 40 m s?1. Gravitational perturbations of planets, Poynting-Robertson effect and solar corpuscular radiation (solar wind) are considered. Other nongravitational effects may be represented by new terms in the suggested form of the nongravitational force. Various values of normal and transversal components of the perturbing nongravitational force are used. The final results of numerical integrations are compared with those obtained on the basis of the Poynting-Robertson effect. 相似文献
5.
A small generalization of the equation of motion for the Poynting-Robertson effect is tested in order to find the significance of new terms. The test is made for dust particles ejected at perihelia of the orbit of the comet Encke. The particles are released at the speed of 40 m s–1. Gravitational perturbations of planets, Poynting-Robertson effect and solar corpuscular radiation (solar wind) are considered. Other nongravitational effects may be represented by new terms in the suggested form of the nongravitational force. Various values of normal and transversal components of the perturbing nongravitational force are used. The final results of numerical integrations are compared with those obtained on the basis of the Poynting-Robertson effect. 相似文献
6.
The stability of collinear and triangular libration points is investigated in the photogravitational elliptic restricted three-body problem, in which two primary bodies emit light energy simultaneously. The conditions of stability of the collinear and triangular libration points are obtained based on a linearized set of equations of perturbed motion for various values of the eccentricity of the Keplerian orbits and the mass ratio of the primary bodies. The maximal numerical value is defined for the eccentricity at which a stable libration point can still exist. It is demonstrated how the parametric resonance causes an instability of collinear and triangular libration points; the evolution of the origination of the instability zones is traced. The minimal eccentricity value is found at which zones of instability of triangular libration points arise. 相似文献
7.
Using the rectangular equations of motion for the restricted three-body problem a comparison is made of the integration of these equations by the Encke method and by a set of perturbational equations. Each set of differential equations is integrated using Taylor series expansions where the coefficients of the powers of time are determined by recurrence relations. It is shown that for very small perturbations the use of the perturbational equations is more efficient than the use of the Encke method. A discussion is also given of when Cowell's method is more efficient than either of these techniques. 相似文献
8.
Edgar Everhart 《Celestial Mechanics and Dynamical Astronomy》1974,10(1):35-55
The solutions of \(\ddot x = F(x,t)\) , and also \(\dot x = F(x,t)\) , are developed in truncated series in timet whose coefficients are found empirically. The series ending in thet 6 term yields a position at a final prechosen time that is accurate through 9th order in the sequence size. This is achieved by using Gauss-Radau and Gauss-Lobatto spacings for the several substeps within each sequence. This timeseries method is the same in principle as implicit Runge-Kutta forms, including some not described previously. In some orders these methods are unconditionally stable (A-stable). In the time-series formulation the implicit system converges rapidly. For integrating a test orbit the method is found to be about twice as fast as high-order explicit Runge-Kutta-Nyström-Fehlberg methods at the same accuracies. Both the Cowell and the Encke equations are solved for the test orbit, the latter being 35% faster. It is shown that the Encke equations are particularly well-adapted to treating close encounters when used with a single-sequence integrator (such as this one) provided that the reference orbit is re-initialized at the start of each sequence. This use of Encke equations is compared with the use of regularized Cowell equations. 相似文献
9.
A.H. Delsemme 《Icarus》1975,24(1):95-110
In order to prepare a flyby mission to Comet Encke, six different sources of information on the possible chemical composition of the cometary nucleus are compared. These are: the neutral and charged radicals and molecules observed in cometary spectra; the chemical composition of type I carbonaceous chondrites; the meteor spectra; the metallic ions collected in the upper atmosphere and correlated with the meteor shower associated with Comet Encke; and finally the volatile molecules observed in a volatile-rich sample of lunar soil, that were interpreted as a possible cometary impact. Possible molecular abundances for the volatile fraction of Comet Encke are tentatively proposed. 相似文献
10.
P. E. Moran 《Celestial Mechanics and Dynamical Astronomy》1973,7(1):122-135
Using the rectangular equations of motion for the restricted three-body problem a comparison is made of the Encke and Cowell methods of integration. Each set of differential equations is integrated using Taylor series expansions where the coefficients of the powers of time are determined by recurrence relationships. It is shown that for fairly highly eccentric orbits in which the perturbing force is less than one thousandth of the two-body force the Encke method achieves a considerable saving in machine time. This is also true for almost circular orbits when low or moderate accuracy is required. When very high accuracy is required, however, the Cowell method is faster unless the perturbing force is less than 10–6 of the two-body force. There is little difference in the accuracy of the two methods, the Cowell method being slightly more accurate when a low or moderate accuracy criterion is imposed. 相似文献
11.
Reduction the secular solution to periodic solution in the generalized restricted three-body problem
Elbaz I. Abouelmagd M. E. Awad E. M. A. Elzayat Ibrahim A. Abbas 《Astrophysics and Space Science》2014,350(2):495-505
The aim of the present work is to find the secular solution around the triangular equilibrium points and reduce it to the periodic solution in the frame work of the generalized restricted thee-body problem. This model is generalized in sense that both the primaries are oblate and radiating as well as the gravitational potential from a belt. We show that the linearized equation of motion of the infinitesimal body around the triangular equilibrium points has a secular solution when the value of mass ratio equals the critical mass value. Moreover, we reduce this solution to periodic solution, as well as some numerical and graphical investigations for the effects of the perturbed forces are introduced. This model can be used to examine the existence of a dust particle near the triangular points of an oblate and radiating binary stars system surrounded by a belt. 相似文献
12.
Masanori Iye Tatsuo Ueda Masafumi Noguchi Shinko Aoki 《Astrophysics and Space Science》1983,89(2):363-376
A new scheme has been devised to calculate discrete unstable global shearing spiral modes of gaseous disk models of galaxies. The scheme makes use of the Legendre expansion of eigen functions and the problem of stability analysis is reduced to an eigenvalue problem of an infinite matrix. The spiral patterns of these shearing wave solutions of linearized equations change their form in the course of time due to the differential rotation of the equilibrium disk. These shearing wave solutions are presumed to have intermediate characteristics between so-called density-waves and material arms. Comparison between these shearing modes and the non-shearing normal modes for a series of disk models is presented. 相似文献
13.
On the origin of the unusual orbit of Comet 2P/Encke 总被引：1，自引：0，他引：1
The orbit of Comet 2P/Encke is difficult to understand because it is decoupled from Jupiter—its aphelion distance is only 4.1 AU. We present a series of orbital integrations designed to determine whether the orbit of Comet 2P/Encke can simply be the result of gravitational interactions between Jupiter-family comets and the terrestrial planets. To accomplish this, we integrated the orbits of a large number of objects from the trans-neptunian region, through the realm of the giant planets, and into the inner Solar System. We find that at any one time, our model predicts that there should be roughly 12 objects in Encke-like orbits. However, it takes roughly 200 times longer to evolve onto an orbit like this than the typical cometary physical lifetime. Thus, we suggest that (i) 2P/Encke became dormant soon after it was kicked inward by Jupiter, (ii) it spent a significant amount of time inactive while rattling around the inner Solar System, and (iii) it only became active again as the ν6 secular resonance drove down its perihelion distance. 相似文献
14.
This paper describes N-body simulations of two regions of the saturnian ring system and examines what we might expect the Cassini orbiter to see in those areas. The first region is the edge of the Encke gap in the A ring that is perturbed by the satellite, Pan. Our previous simulations of this region neglected particle self-gravity [Lewis and Stewart, 2000a, Bull. Am. Astron. Soc. 34, 883]. Here we examine the interactions of the wakes caused by Pan with the wakes that form from local gravitational instabilities. We find that the two phenomena do not normally coexist and predict that measurements of particle sizes between the moon wakes should reflect the true particle size distribution of the region and not what is caused by gravitational aggregation. The region between the Encke gap edge and the first wake peak is an exception to this rule because our simulations exhibit the formation of exceptionally large gravity-induced wakes in this region. We also describe simulations of the F ring and explain the nature of braid-like structures that form naturally when the ring is perturbed by a single moon on an eccentric orbit. Finally, we discuss the very dynamic nature of the F ring system and how this should be taken into account when interpreting observations and even when planning future observations of this system. 相似文献
15.
Stellar occultations by Saturn’s rings observed with the Visual and Infrared Mapping Spectrometer (VIMS) onboard the Cassini spacecraft reveal that dusty features such as the F ring and the ringlets in the Encke and the Laplace Gaps have distinctive infrared transmission spectra. These spectra show a narrow optical depth minimum at wavelengths around 2.87 μm. This minimum is likely due to the Christiansen Effect, a reduction in the extinction of small particles when their (complex) refractive index is close to that of the surrounding medium. Simple Mie-scattering models demonstrate that the strength of this opacity dip is sensitive to the size distribution of particles between 1 and 100 μm across. Furthermore, the spatial resolution of the occultation data is sufficient to reveal variations in the transmission spectra within and among these rings. In both the Encke Gap ringlets and F ring, the opacity dip weakens with increasing local optical depth, which is consistent with the larger particles being concentrated near the cores of these rings. The Encke Gap ringlets also show systematically weaker opacity dips than the F ring and Laplace Gap ringlet, implying that the former has a smaller fraction of grains less than ∼30 μm across. However, the strength of the opacity dip varies most dramatically within the F ring; certain compact regions of enhanced optical depth lack an opacity dip and therefore appear to have a greatly reduced fraction of grains in the few-micron size range. Such spectrally-identifiable structures probably represent a subset of the compact optically-thick clumps observed by other Cassini instruments. These variations in the ring’s particle size distribution can provide new insights into the processes of grain aggregation, disruption and transport within dusty rings. For example, the unusual spectral properties of the F-ring clumps could perhaps be ascribed to small grains adhering onto the surface of larger particles in regions of anomalously low velocity dispersion. 相似文献
16.
The main limit to the time span of a numerical integration of the planetary orbits is no longer set by the availability of computer resources, but rather by the accumulation of the integration error. By the latter we mean the difference between the computed orbit and the dynamical behaviour of the real physical system, whatever the causes. The analysis of these causes requires an interdisciplinary effort: there are physical model and parameters errors, algorithm and discretisation errors, rounding off errors and reliability problems in the computer hardware and system software, as well as instabilities in the dynamical system. We list all the sources of integration error we are aware of and discuss their relevance in determining the present limit to the time span of a meaningful integration of the orbit of the planets. At present this limit is of the order of 108 years for the outer planets. We discuss in more detail the truncation error of multistep algorithms (when applied to eccentric orbits), the coefficient error, the method of Encke and the associated coordinate change error, the procedures used to test the numerical integration software and their limitations. Many problems remain open, including the one of a realistic statistical model of the rounding off error; at present, the latter can only be described by a semiempirical model based upon the simpleN
2 formula (N=number of steps, =machine accuracy), with an unknown numerical coefficient which is determined only a posteriori. 相似文献
17.
George W. Wetherill 《Icarus》1979,37(1):96-112
The steady-state distribution of orbits of Apollo-Amor objects is calculated for a variety of possible sources. These include asteroids near the inner edge of the belt, cometary orbits similar to Encke, and hypothetical extinct cometary orbits with perihelia larger than that of Encke. In all but one case, the steady-state distributions are similar for all these sources, and predict Amor/Apollo ratios of 1.5 to 3. These ratios are lower than those predicted by work in which the effects of the ν6 secular resonance were not considered. These results are in general agreement with observation, although the higher (～3) Amor/Apollo ratios found for many of the sources may turn out to be unacceptably high. The absolute number of Apollo-Amors observed is found to require an injection rate of ～15 objects/(106 years). This rate is easily achieved if the present existence of Encke is assumed to be a reasonably probable event, and if Encke becomes a ～1-km-diameter Apollo object following exhaustion of its volatile material; best estimates of the injection rate from the asteroid belt [～1.5/(106 years)] are too low. Hence a dominant cometary component is suggested. The predicted number of Apollo objects in small (q < 1.0 AU, a < 1.4 AU orbits is in agreement with observation. Predicted lunar and terrestrial cratering rates agree approximately with observation. An unexplained difference between the lunar and terrestrial results is probably caused by uncertainties in the scaling laws or crater counts used. This discrepancy precludes an exact test of these calculations using cratering data. 相似文献
18.
In early January 1888, James E. Keeler was one of the first astronomers to work with the very new Lick Observatory 36-in. refractor. On January 7 while observing Saturn visually on a night of very fine seeing, he discovered a narrow, dark “division” in the outer part of the A ring. Despite repeated attempts, neither Keeler nor any of the other Lick observers saw this gap again until over a year later, on March 2, 1889, another night of extremely good seeing. On that occasion not only Keeler, but also E. S. Holden, J. M. Schaeberle, and E. E. Barnard all observed “Mr. Keeler's division,” as Barnard called it. It could only be seen using very high magnification with this large telescope, at a site known to be excellent, on the nights of very best definition. This gap is not the same as the feature which J. F. Encke had earlier discovered and described as a low-contrast division nearly in the middle of the A ring, and had drawn as nearly the same width as Cassini's division. Later visual observations by B. Lyot and A. Dollfus, again on nights of fine seeing with large telescopes, showed that the Encke division is complex. To them, with the best resolution, it appeared as three wide minima of light, fuzzy, and of low contrast, with a narrow, well-marked minimum of light at its outer edge. The outer edge is just where Keeler placed his gap, although he did not see the low-contrast structure in the Encke division. The images, with much superior resolution obtained from the Pioneer and Voyager space probes, show that the Encke division is even more complex than Lyot and Dollfus realized, but confirm the narrow Keeler feature as a true gap in the outer part of the A ring. 相似文献
19.
Jozef Klačka 《Earth, Moon, and Planets》1994,65(2):191-196
The action of the solar electromagnetic radiation on the moving interplanetary dust particles in its more complete form than the special case known as the Poynting-Robertson effect is theoretically discussed in application to meteoroid stream of comet Encke.Normal and transversal components of the perturbing nongravitational force are used due to the action of the solar electromagnetic radiation. It is shown that the normal component of the force is negligible. However, transversal component is very important: it can probably completely explain all the observed meteoroid streams situated along the orbit of comet Encke (and, possibly, some asteroids) as the product of the comet Encke alone. Much shorter time is required for producing such a meteoroid stream than is a general conception.If the idea about the significance of the transversal component of the nongravitational force (may be, not produced by electromagnetic radiation) is correct, it may have important consequences for our understanding of ageing of comets, global evolution of the cometary (and, partially, asteroidal) system, and, of course, for a long-term evolution of small interplanetary particles. 相似文献
20.
A comparison of the ultraviolet spectrum of periodic Comet Encke (1980 XI), recorded by the IUE between 24 October and 5 November 1980 with similar spectra of short- and long-period comets shows the gaseous composition of P/Encke to be nearly identical to that of the other comets observed by the IUE. If P/Encke is indeed the remains of a once giant comet, this similarity implies a homogeneous radial structure for the cometary ice nucleus. The OH brightness distribution shows a spatial variation similar to the visible fan-shaped image of the comet, suggestive of a nonuniform distribution of volatile ices on the surface of the nucleus. The total derived water production rate appears to be a factor of 5 higher than that derived from HI Lyman-α observations made during the 1970 apparition and shows a variation with heliocentric distance (r) as r?3.3 over the range 0.81 to 1.02 AU. 相似文献