首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
Julio A. Fernández 《Icarus》1980,42(3):406-421
The orbital evolution of 500 hypothetical comets during 109 years is studied numerically. It is assumed that the birthplace of such comets was the region of Uranus and Neptune from where they were deflected into very elongated orbits by perturbations of these planets. Then, we adopted the following initial orbital elements: perihelion distances between 20 and 30 AU, inclinations to the ecliptic plane smaller than 20°, and semimajor axes from 5 × 103 to 5 × 104 AU. Gravitational perturbations by the four giant planets and by hypothetical stars passing at distances from the Sun smaller than 5 × 105 AU are considered. During the simulation, somewhat more than 50% of the comets were lost from the solar system due to planetary or stellar perturbations. The survivors were removed from the planetary region and left as members of what is generally known as the cometary cloud. At the end of the studied period, the semimajor axes of the surviving comets tend to be concentrated in the interval 2 × 104 < a < 3 × 104 AU. The orbital planes of the comets with initial a ≧ 3 × 104AU acquired a complete randomization while the others still maintain a slight predominance of direct orbits. In addition, comet orbits with final a < 6 × 104AU preserve high eccentricities with an average value greater than 0.8 Most “new” comets from the sample entering the region interior to Jupiter's orbit had already registered earlier passages through the planetary region. By scaling up the rate of paritions of hypothetical new comets with the observed one, the number of members of the cometary cloud is estimated to be about 7 × 1010 and the conclusion is drawn that Uranus and Neptune had to remove a number of comets ten times greater.  相似文献   

2.
The influence of aerodynamic drag and the geopotential on the motion of the satellite 1964-52B is considered. A model of the atmosphere is adopted that allows for oblateness, and in which the density behaviour approximates to the observed diurnal variation. A differential equation governing the variation of the eccentricity, e, combining the effects of air drag with those of the Earth's gravitational field is given. This is solved numerically using as initial conditions 310 computed orbits of 1964-52B.The observed values of eccentricity are modified by the removal of perturbations due to luni-solar attraction, solid Earth and ocean tides, solar radiation pressure and low-order long-periodic tesseral harmonic perturbations. The method of removal of these effects is given in some detail. The behaviour of the orbital eccentricity predicted by the numerical solution is compared with the modified observed eccentricity to obtain values of atmospheric parameters at heights between 310 and 430 km. The daytime maximum of air density is found to be at 14.5 hours local time. Analysis of the eccentricity near 15th order resonance with the geopotential yielded values of four lumped geopotential harmonics of order 15, namely: 109C1,015 = ?78.8 ± 7.0, 109S1,015 = ?69.4 ± 5.3, 109C?1,215 = ?41.6 ± 3.5109S?1,215 = ?26.1 ± 8.9, at inclination 98.68°.  相似文献   

3.
From published ground-base, spacecraft, and rocket photometry and polarimetry of the zodiacal light, a number of optical and physical parameters have been derived. It was assumed that the number density, mean particle size, and albedo vary with heliocentric distance, and shown that average individual interplanetary particles have a small but definite opposition effect, a mean single-scattering albedo in the V band at 1-AU heliocentric distance of 0.09 ± 0.01, and a zero-phase geometric albedo of 0.04. Modeled by a power law, both albedos decrease with increasing heliocentric distance as r?0.54. The corresponding exponents for changes in mean particle size and number density are related in a simple way. The median orbital inclination of zodiacal light particles with respect to the ecliptic is 12°, close to the observed median value for faint asteroids and short-period comets. Furthermore, the color of dust particles and its variation with solar phase angle closely resemble those of C asteroids. These findings are, at least, consistent with the zodiacal cloud originating primarily from collisions among asteroids. Finally, a value of ?1018?ErmE g was derived for the mass of the zodiacal cloud, where ?E is the mean particle radius (in micrometers) at 1-AU-heliocentric distance. For extinction in the ecliptic, Δm = 10?5??12mag was obtained, where ? is the solar elongation in degrees.  相似文献   

4.
I consider the range of Hill stability in the restricted circular problem of three bodies when the larger one of the two principal bodies has a finite oblateness. I show that the range r satisfies the equation
r = 1? μCcr ? 3μ + μr2 ? v (1? μ)r?3 ±(2 + 3v)(1 ? μ 1 + 3vr2r,
where μ is the mass parameter and v is an oblateness parameter. This result is applied to the solar system, the Earth-Moon system and binary star systems. It is then shown that, all the inner planets of the solar system, the great majority of asteroids and some short-period comets are Hill stable, that direct artificial satellites of the Earth are more stable than retrograde ones, and that contact binaries possess cores between which no mass exchange takes place.  相似文献   

5.
A theory is presented for charged-particle collection by a cylindrical conducting object, such as a spacecraft or an electrostatic probe, which is moving transversely through a collisionless plasma, such as those in the upper atmosphere and space. The calculation is approximate, using symmetric potential profiles which are exact for the infinite-cylinder stationary case. Theoretical current predictions are presented for ratios of collector potential to electron thermal energy c/kTe from 0 to ?25, for ion-to-electron temperature ratios Ti/Tc = 1 and 0.5, ratio of collector radius to electron Debye length rc/λD from 0 to 100, and ratio of flow speed to ion thermal speed Si = U/(2kTi/mi12) from 0 to 10. Comparisons with existing exact calculations by other authors show that none of these fulfil all of the requirements for nontrivial comparison. Appropriate parameter ranges for future exact calculations are thereby suggested. These are as follows: (a) rc/λD should be large enough that the collector not be in or near orbit-limited conditions; (b) the ratio Si2/¦χc, i¦ of ion directed energy to potential energy change in the sheath, should be close to unity or if
Si2/¦χc,i¦? 1, then Si ? 1
.  相似文献   

6.
7.
We study numerically and analytically the conditions leading to the truncation, at the 2:1 resonance, of a disk of infrequently colliding particles surrounding the primary of a binary system. We focus on systems with small mass ratios, q, such as the Sun-Jupiter system with q = 10?3. Previous studies showed that if collisions are frequent with respect to the orbital period, truncation 3nly occurs if the Reynolds number is greater than q?2. This corresponds to particle eccentricity, e, less than of order q for a particle disk of optical depth unity. In thepresent case collisions are less frequent than q?23 orbital periods (the period of the forced eccentricity at the 2:1 resonance), and truncation occurs and (Kirkwood) gaps are produced only if e is less than some critical value which we estimate to be of order q59, or ~0.02 for the Sun-Jupiter case. We mention several means whereby the eccentricities may have been subsequently increased to their observed values.  相似文献   

8.
Editorial     
The Galilean satellites Io, Europa, and Ganymede interact through several stable orbital resonances where λ1 ? 2λ2 + ω1 = 0, λ1 ? 2λ2 + ω2 = 180°, λ2 ? 2λ3 + ω2 = 0 and λ1 ? 3λ2 + 2λ3 = 180°, with λi being the mean longitude of the ith satellite and ωi the longitude of the pericenter. The last relation involving all three bodies is known as the Laplace relation. A theory of origin and subsequent evolution of these resonances outlined earlier (C. F. Yoder, 1979b, Nature279, 747–770) is described in detail. From an initially quasi-random distribution of the orbits the resonances are assembled through differential tidal expansion of the orbits. Io is driven out most rapidly and the first two resonance variables above are captured into libration about 0 and 180° respectively with unit probability. The orbits of Io and Europa expand together maintaining the 2:1 orbital commensurability and Europa's mean angular velocity approaches a value which is twice that of Ganymede. The third resonance variable and simultaneously the Laplace angle are captured into libration with probability ~0.9. The tidal dissipation in Io is vital for the rapid damping of the libration amplitudes and for the establishment of a quasi-stationary orbital configuration. Here the eccentricity of Io's orbit is determined by a balance between the effects of tidal dissipation in Io and that in Jupiter, and its measured value leads to the relation k1?1/Q1 ≈ 900kJ/QJ with the k's being Love numbers, the Q's dissipation factors, and f a factor to account for a molten core in Io. This relation and an upper bound on Q1 deduced from Io's observed thermal activity establishes the bounds 6 × 104 < QJ < 2 × 106, where the lower bound follows from the limited expansion of the satellite orbits. The damping time for the Laplace libration and therefore a minimum lifetime of the resonance is 1600 QJ years. Passage of the system through nearby three-body resonances excites free eccentricities. The remnant free eccentricity of Europa leads to the relation Q2/?2 ? 2 × 10?4 QJ for rigidity μ2 = 5 × 1011 dynes/cm2. Probable capture into any of several stable 3:1 two-body resonances implies that the ratio of the orbital mean motions of any adjacent pair of satellites was never this large.A generalized Hamiltonian theory of the resonances in which third-order terms in eccentricity are retained is developed to evaluate the hypothesis that the resonances were of primordial origin. The Laplace relation is unstable for values of Io's eccentricity e1 > 0.012 showing that the theory which retains only the linear terms in e1 is not valid for values of e1 larger than about twice the current value. Processes by which the resonances can be established at the time of satellite formation are undefined, but even if primordial formation is conjectured, the bounds established above for QJ cannot be relaxed. Electromagnetic torques on Io are also not sufficient to relax the bounds on QJ. Some ideas on processes for the dissipation of ideal energy in Jupiter yield values of QJ within the dynamical bounds, but no theory has produced a QJ small enough to be compatible with the measurements of heat flow from Io given the above relation between Q1 and QJ. Tentative observational bounds on the secular acceleration of Io's mean motion are also shown not to be consistent with such low values of QJ. Io's heat flow may therefore be episodic. QJ may actually be determined from improved analysis of 300 years of eclipse data.  相似文献   

9.
The beam cyclotron instability and electron acoustic instability, driven by cross-tail current and inhomogeneity in density and magnetic field, are found to be unstable in the earth's magnetic tail region. The anomalous resistivities due to these instabilities are found to be of the order of (10?1?10?3e?1e being the electron gyro frequency). It is also suggested that the non-linear saturation of the beam cyclotron instability may lead to conditions favourable for exciting ion acoustic instability.  相似文献   

10.
Gravitational perturbations in semimajor axis, eccentricity, and inclination resulting from close planetesimal encounters (near 1 AU) out to 10 Tisserand sphere of influence radii were calculated by two- and three-dimensional numerical integration. These are compared with the results of treating the encounter as a two-body problem, as is customary in Monte Carlo calculations of orbital evolution and in numerical and analytical studies of planetary accumulation. It is found that for values of (VVe) ? 0.35 (V = relative velocity, Ve = escape velocity of largest body), the two-body body approximation fails to describe the outcome of individual encounters. In this low-velocity region, the two-body “gravitational focusing” cross section is no longer valid; “anomalous gravitational focusing” often leads to bodies on distant unperturbed trajectories becoming close encounters and vice versa. In spite of these differences, average perturbations given by the two-body approximation are valid within a factor of 2 when VVe > 0.07. In this same velocity range the “Arnold extrapolation,” whereby a few very close encounters are used to estimate the effect of many more distant encounters, is found to be a useful approximation.  相似文献   

11.
The influence of aerodynamic drag and the geopotential on the motion of the satellite 1964-52B is considered. A model of the atmosphere is adopted that allows for oblateness, and in which the density behaviour approximates to the observed diurnal variation. A differential equation governing the variation of the orbital inclination combining the effects of air drag with those of the Earth's gravitational field is given.The 310 observed values of inclination are modified by the removal of perturbations due to luni-solar attraction, solid Earth and ocean tides, solar radiation pressure, low-order long-periodic tesseral harmonic perturbations and changes due to precession. The method of removal of these effects is given in some detail.The variations in inclination due to drag are analysed to give four values of the average atmospheric rotation rate at heights of 296–476 km at latitude 0–54°. These values are as expected from previous analyses.The analysis of the change in inclination due to solar radiation pressure shows that this rapidly tumbling cylindrical satellite may be considered as equivalent to a spherical satellite of a given area-to-mass ratio.Analysis of the inclination near 15:1 resonance with the geopotential yields values of lumped geopotential harmonics of order 15 and 30, namely, 109C?0.115 = ?31.2 ± 2.3 109S?0.115 = ?4.4 ± 3.2 109C?0.230 = 39.0 ± 10.7 109S?0.230 = 51.8 ± 10.0  相似文献   

12.
Ariel 1, the first international satellite, was launched on 26 April 1962, into an orbit inclined at 53.85° to the equator, with an initial perigee height near 390 km. On 8 May 1973 the orbit passed through 15th-order resonance and has been determined, with the RAE orbit refinement program PROP, at eight epochs between February and August 1973 using 500 observations.The orbital inclinations during the time of 15th-order resonance, as given by these eight orbits and 31 U.S. Navy orbits, were fitted with a theoretical curve using the THROE computer program, the best fit giving 109C?15 = ?370 ± 14 and 109S15 = ?114 ± 31.The values of eccentricity were also successfully fitted using THROE, and the results are discussed.  相似文献   

13.
Oort cloud comets occasionally obtain orbits which take them through the planetary region. The perturbations by the planets are likely to change the orbit of the comet. We model this process by using a Monte Carlo method and cross sections for orbital changes, i.e. changes in energy, inclination and perihelion distance, in a single planet-comet encounter. The influence of all major planets is considered. We study the distributions of orbital parameters of observable comets, i.e. those which have perihelion distance smaller than a given value. We find that enough comets are captured from the Oort cloud in order to explain the present populations of short period comets. The median value of cos i for the Jupiter family is 0.985 while it is 0.27 for the Halley types. The results may explain the orbital features of short period comets, assuming that the active lifetime of a comet is not much greater than 400 orbital revolutions.  相似文献   

14.
Oort cloud comets occasionally obtain orbits which take them through the planetary region. The perturbations by the planets are likely to change the orbit of the comet. We model this process by using a Monte Carlo method and cross sections for orbital changes, i.e. changes in energy, inclination and perihelion distance, in a single planet-comet encounter. The influence of all major planets is considered. We study the distributions of orbital parameters of observable comets, i.e. those which have perihelion distance smaller than a given value. We find that enough comets are captured from the Oort cloud in order to explain the present populations of short period comets. The median value of cos i for the Jupiter family is 0.985 while it is 0.27 for the Halley types. The results may explain the orbital features of short period comets, assuming that the active lifetime of a comet is not much greater than 400 orbital revolutions.  相似文献   

15.
Secular perturbations of asteroids are derived for mean motion resonance cases under the assumptions that the disturbing planets are moving along circular orbits on the same plane and that critical arguments are fixed at stable equilibrium points. Under these assumptions the equations of motion are reduced to those of one degree of freedom with the energy integral. Then equi-energycurves on (2g-X) plane (g and X being, respectively, the argument of perihelion and (1–e2)1/2) are derived for given values of the two constant parameters, the semi-major axis and =(1–e2)1/2 cos i, and the variations of the eccentricity and the inclination as functions of the argument of perihelion are graphically estimated. In fact this method is applied to numbered asteroids with commensurable mean motions to estimate the ranges of the variations of orbital elements.The same method is also applied to the Pluto-Neptune system and the results are found to agree with those of numerical integrations and show that the argument of perihelion of Pluto librates around 90°.  相似文献   

16.
The orbit of Cosmos 837 rocket (1976-62E) has been determined at 36 epochs between January and September 1978, using the RAE orbit refinement program PROP 6 with about 3000 observations. The inclination was 62.7° and the eccentricity 0.039. The orbital accuracy achieved was between 30m and 150m, both radial and crosstrack. The orbit was near 29:2 resonance in 1978 (exact resonance occurred on 14 May) and the values of orbital inclination obtained have been analysed to derive lumped 29th-order geopotential harmonic coefficients, namely:
109C0,229 = ? 10 ± 15
and
109S0,229 = ?76 ± 12
. These will be used in future, when enough results at different inclinations have accumulated, to determine individual coefficients of order 29. The values of lumped harmonics obtained from analysis of the values of eccentricity were not well defined, because of the high correlations between them and the errors in removing the very large perturbation (31 km) due to odd zonal harmonics.  相似文献   

17.
T.E. Cravens  A.E.S. Green 《Icarus》1978,33(3):612-623
The intensities of radiation from the inner comas of comets which are composed primarily of water and carbon monoxide have been calculated. Only “airglow” emissions initiated by the absorption of extreme ultraviolet radiation have been considered. The photoionizations of H2O, CO, CO2, and N2 are the most important emission sources, although photoelectron excitation is also considered. Among the emission features for which intensities were calculated are H2O+ (A?2A1?X?2B1), CO+ (first negative), CO (fourth positive), CO (Cameron), CO2+ (B?2?u?X?2IIg), N2 (Vegard-Kaplan), N2+ (first negative), and OI (1304 Å). In the inner coma (collision region) these airglow mechanisms are shown to be possible competitors with the usually assumed resonance scattering and flourescence excitation mechanisms which are appropriate for the outer coma and tail.  相似文献   

18.
The magnetopause, the boundary layer, or current sheath, which separates the magnetosphere from the solar wind, is the particular interaction considered in this paper.The collision free electron skin depth, ξe = cωpe, where c is the velocity of light and ωpe, is the plasma frequency, gives a classical measure of the penetration depth of a collisionless plasma by an electromagnetic field. This penetration depth is small compared with the dimensions of the magnetosphere and hence the boundary layer may be conveniently considered in one dimension.In General all one dimensional solutions lie within an order of magnitude of the value of ξe, the only exception being the important one, in which the electric field perpendicular to the current sheath plane is not present, either due to a particular trapped particle distribution or due to a short circuiting end effect. For this exception the thickness is increased by the factor (mii/me)12.The current sheath solutions discussed are equilibrium solutions but not necessarily stable equilibrium solutions.The extension of the models to three dimensions has a larger effect than might at first be expected. The effect may be intuitively understood as a consequence of flux conservation in the sheath. The one dimensional solutions then correspond to the current sheath profiles at the thinnest point of the three dimensional sheath.  相似文献   

19.
Six times of maxima of the ultrashort-period cepheid variable EH Librae were measured in 1980 May to June and in 1981 January, with a three-channel photocounting high-speed photoelectric photometer. These, together with all the photoelectric times of maxima over the past 30 years, are used to re-examine the nature of the change of the period. We found that we can fix the times of maxima by the following formula
Tmax = T0+P0E+12βE2+AsinEP0E0
where T0 = HJD 2433438.6088 and P0 = 0.0884132445 d represent the initial maximum epoch and the pulsation period, β = ?2.8 × 10?8/yr; A = 0.0015 d, P0 = 6251 d = 17.1 yr are the semi-amplitude and the period of the sine curve, and E is the number of periods elapsed since T0, and (E0 = 70700).If we interpret this 17.1 year periodicity as a modulation of the phase of maximum by binary motion, then the semi-amplitude of the orbital radial velocity variation is K = 2πasini/E0 = 0.45 km/s and the mass function is
f(m)=m32sin3i(m1m2)2=(asini)3E20=6 x 10?5M
  相似文献   

20.
Previous work has parameterized the pitch angle dependence of the charge-exchange lifetime τ of ring current ions in terms of γ, the power of the cosine of the mirror latitude λm of the particles, such that τ(λm)τ(0) ≌cosγ λm at given L. Using the atomic hydrogen density model of Johnson and Fish, previous authors have suggested values of γ = 5 or 6. We here evaluate γ as a function of λm and L using the more recent Chamberlain density models, and show that γ = 3?4 is more appropriate over most of the pitch angle and L range. Consequently, ion distributions in the ring current decay phase are expected to become rather less anisotropic in pitch angle due to chargeexchange than previously believed. We have also investigated the use of several other simple approximate analytic forms for τ(λm)τ(0), one of which gives far better agreement with the numerical results than the cosγ λm, variation, and should hence be used in future studies.  相似文献   

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

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