首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
In the method of variation of parameters we express the Cartesian coordinates or the Euler angles as functions of the time and six constants. If, under disturbance, we endow the “constants” with time dependence, the perturbed orbital or angular velocity will consist of a partial time derivative and a convective term that includes time derivatives of the “constants”. The Lagrange constraint, often imposed for convenience, nullifies the convective term and thereby guarantees that the functional dependence of the velocity on the time and “constants” stays unaltered under disturbance. “Constants” satisfying this constraint are called osculating elements. Otherwise, they are simply termed orbital or rotational elements. When the equations for the elements are required to be canonical, it is normally the Delaunay variables that are chosen to be the orbital elements, and it is the Andoyer variables that are typically chosen to play the role of rotational elements. (Since some of the Andoyer elements are time-dependent even in the unperturbed setting, the role of “constants” is actually played by their initial values.) The Delaunay and Andoyer sets of variables share a subtle peculiarity: under certain circumstances the standard equations render the elements nonosculating. In the theory of orbits, the planetary equations yield nonosculating elements when perturbations depend on velocities. To keep the elements osculating, the equations must be amended with extra terms that are not parts of the disturbing function [Efroimsky, M., Goldreich, P.: J. Math. Phys. 44, 5958–5977 (2003); Astron. Astrophys. 415, 1187–1199 (2004); Efroimsky, M.: Celest. Mech. Dyn. Astron. 91, 75–108 (2005); Ann. New York Acad. Sci. 1065, 346–374 (2006)]. It complicates both the Lagrange- and Delaunay-type planetary equations and makes the Delaunay equations noncanonical. In attitude dynamics, whenever a perturbation depends upon the angular velocity (like a switch to a noninertial frame), a mere amendment of the Hamiltonian makes the equations yield nonosculating Andoyer elements. To make them osculating, extra terms should be added to the equations (but then the equations will no longer be canonical). Calculations in nonosculating variables are mathematically valid, but their physical interpretation is not easy. Nonosculating orbital elements parameterise instantaneous conics not tangent to the orbit. (A nonosculating i may differ much from the real inclination of the orbit, given by the osculating i.) Nonosculating Andoyer elements correctly describe perturbed attitude, but their interconnection with the angular velocity is a nontrivial issue. The Kinoshita–Souchay theory tacitly employs nonosculating Andoyer elements. For this reason, even though the elements are introduced in a precessing frame, they nevertheless return the inertial velocity, not the velocity relative to the precessing frame. To amend the Kinoshita–Souchay theory, we derive the precessing-frame-related directional angles of the angular velocity relative to the precessing frame. The loss of osculation should not necessarily be considered a flaw of the Kinoshita–Souchay theory, because in some situations it is the inertial, not the relative, angular velocity that is measurable [Schreiber, K. U. et al.: J. Geophys. Res. 109, B06405 (2004); Petrov, L.: Astron. Astrophys. 467, 359–369 (2007)]. Under these circumstances, the Kinoshita–Souchay formulae for the angular velocity should be employed (as long as they are rightly identified as the formulae for the inertial angular velocity).  相似文献   

2.
We continue the study undertaken in Efroimsky [Celest. Mech. Dyn. Astron. 91, 75–108 (2005a)] where we explored the influence of spin-axis variations of an oblate planet on satellite orbits. Near-equatorial satellites had long been believed to keep up with the oblate primary’s equator in the cause of its spin-axis variations. As demonstrated by Efroimsky and Goldreich [Astron. Astrophys. 415, 1187–1199 (2004)], this opinion had stemmed from an inexact interpretation of a correct result by Goldreich [Astron. J. 70, 5–9 (1965)]. Although Goldreich [Astron. J. 70, 5–9 (1965)] mentioned that his result (preservation of the initial inclination, up to small oscillations about the moving equatorial plane) was obtained for non-osculating inclination, his admonition had been persistently ignored for forty years. It was explained in Efroimsky and Goldreich [Astron. Astrophys. 415, 1187–1199 (2004)] that the equator precession influences the osculating inclination of a satellite orbit already in the first order over the perturbation caused by a transition from an inertial to an equatorial coordinate system. It was later shown in Efroimsky [Celest. Mech. Dyn. Astron. 91, 75–108 (2005a)] that the secular part of the inclination is affected only in the second order. This fact, anticipated by Goldreich [Astron. J. 70, 5–9 (1965)], remains valid for a constant rate of the precession. It turns out that non-uniform variations of the planetary spin state generate changes in the osculating elements, that are linear in , where is the planetary equator’s total precession rate that includes the equinoctial precession, nutation, the Chandler wobble, and the polar wander. We work out a formalism which will help us to determine if these factors cause a drift of a satellite orbit away from the evolving planetary equator.By “precession,” in its most general sense, we mean any change of the direction of the spin axis of the planet—from its long-term variations down to nutations down to the Chandler wobble and polar wander.  相似文献   

3.
Using a consistent perturbation theory for collisionless disk-like and spherical star clusters, we construct a theory of slow modes for systems having an extended central region with a nearly harmonic potential due to the presence of a fairly homogeneous (on the scales of the stellar system) heavy, dynamically passive halo. In such systems, the stellar orbits are slowly precessing, centrally symmetric ellipses (2: 1 orbits). Depending on the density distribution in the system and the degree of halo inhomogeneity, the orbit precession can be both prograde and retrograde, in contrast to systems with 1: 1 elliptical orbits where the precession is unequivocally retrograde. In the first paper, we show that in the case where at least some of the orbits have a prograde precession and the stellar distribution function is a decreasing function of angular momentum, an instability that turns into the well-known radial orbit instability in the limit of low angular momenta can develop in the system. We also explore the question of whether the so-called spoke approximation, a simplified version of the slow mode approximation, is applicable for investigating the instability of stellar systems with highly elongated orbits. Highly elongated orbits in clusters with nonsingular gravitational potentials are known to be also slowly precessing 2: 1 ellipses. This explains the attempts to use the spoke approximation in finding the spectrum of slow modes with frequencies of the order of the orbit precession rate. We show that, in contrast to the previously accepted view, the dependence of the precession rate on angular momentum can differ significantly from a linear one even in a narrow range of variation of the distribution function in angular momentum. Nevertheless, using a proper precession curve in the spoke approximation allows us to partially “rehabilitate” the spoke approach, i.e., to correctly determine the instability growth rate, at least in the principal (O(α T−1/2) order of the perturbation theory in dimensionless small parameter α T, which characterizes the width of the distribution function in angular momentum near radial orbits.  相似文献   

4.
The present paper addresses the existence of J 2 invariant relative orbits with arbitrary relative magnitude over the infinite time using the Routh reduction and Poincaré techniques in the J 2 Hamiltonian problem. The current research also proposes a novel numerical searching approach for J 2 invariant relative orbits from the dynamical system point of view. A new type of Poincaré mapping is defined from different central manifolds of the pseudo-circular orbits (parameterized by the Jacobi energy E, the polar component of momentum H z and the measure of distance Δr between the fixed point and its central manifolds) to the nodal periods T d and the drifts of longitude of the ascending node during one period (ΔΩ), which differs from Koon et al.’s (AIAA 2001) definition on central manifolds parameterized by the same fixed point. The Poincaré mapping is surjective because it compresses the three-dimensional variables into two-dimensional images, and the mapping degenerates into a bijective mapping in consideration of the fixed points. An iteration algorithm to the degenerated bijective mapping is proposed from the continuation procedure to perform the ergodic representation of E- and H z -contour maps on the space of T d –ΔΩ. For the surjective mapping with Δr ≠ 0, different pseudo-circular or elliptical orbits may share the same images. Hence, the inverse surjective mapping may achieve non-unique variables from a single image, which makes the generation of J 2 invariant relative orbits possible. The pseudo-circular or elliptical orbits generated from the surjective mapping will be defined in different meridian planes. Hence, the critical contribution of the present paper is the assignment of J 2 invariant relative orbits to different invariant parameters E and H z depending on the E- and H z -contour map, which will hold J 2 invariant relative orbits for extended durations. To investigate the high-order nonlinearity neglected by previous studies, a formation configuration with a large magnitude of 500 km is successfully generated from the theory developed in the present work, which is beyond the scope of the linear conditions of J 2 invariant relative orbits. Therefore, the existence of J 2 invariant relative orbit with an arbitrary relative magnitude over the infinite time is achieved from the dynamical system point of view.  相似文献   

5.
The role of the angular momentum in the regular or chaotic character of motion in an axially symmetric quasar model is examined. It is found that, for a given value of the critical angular momentumL zc , there are two values of the mass of the nucleusM n for which transition from regular to chaotic motion occurs. The [L zc – M n ] relationship shows a linear dependence for the time independent model and an exponential dependence for the evolving model. Both cases are explained using theoretical arguments together with some numerical evidence. The evolution of the orbits is studied, as mass is transported from the disk to the nucleus. The results are compared with the outcomes derived for galactic models with massive nuclei.  相似文献   

6.
Techniques from dynamical systems theory have been applied to the construction of transfers between unstable periodic orbits that have different energies. Invariant manifolds, trajectories that asymptotically depart or approach unstable periodic orbits, are used to connect the initial and final orbits. The transfer asymptotically departs the initial orbit on a trajectory contained within the initial orbit’s unstable manifold and later asymptotically approaches the final orbit on a trajectory contained within the stable manifold of the final orbit. The manifold trajectories are connected by the execution of impulsive maneuvers. Two-body parameters dictate the selection of the individual manifold trajectories used to construct efficient transfers. A bounding sphere centered on the secondary, with a radius less than the sphere of influence of the secondary, is used to study the manifold trajectories. A two-body parameter, κ, is computed within the bounding sphere, where the gravitational effects of the secondary dominate. The parameter κ is defined as the sum of two quantities: the difference in the normalized angular momentum vectors and eccentricity vectors between a point on the unstable manifold and a point on the stable manifold. It is numerically demonstrated that as the κ parameter decreases, the total cost to complete the transfer decreases. Preliminary results indicate that this method of constructing transfers produces a significant cost savings over methods that do not employ the use of invariant manifolds.  相似文献   

7.
We study numerically the asymmetric periodic orbits which emanate from the triangular equilibrium points of the restricted three-body problem under the assumption that the angular velocity ω varies and for the Sun–Jupiter mass distribution. The symmetric periodic orbits emanating from the collinear Lagrangian point L 3, which are related to them, are also examined. The analytic determination of the initial conditions of the long- and short-period Trojan families around the equilibrium points, is given. The corresponding families were examined, for a combination of the mass ratio and the angular velocity (case of equal eigenfrequencies), and also for the critical value ω = 2
, at which the triangular equilibria disappear by coalescing with the inner collinear equilibrium point L 1. We also compute the horizontal and the vertical stability of these families for the angular velocity parameter ω under consideration. Series of horizontal–critical periodic orbits of the short-Trojan families with the angular velocity ω and the mass ratio μ as parameters, are given.  相似文献   

8.
The values of the initial velocity of the meteoroids ejected from the parent bodies are small and as a result, the most of the young meteoroid streams have similar orbits to their parent bodies. Assuming that the members of the observed meteor stream evolved under the influence of gravitational perturbations mostly, Pittich [1991, Proceedings of the Conference on Dynamic of Small Bodies of the Solar System, Polish-Slovak Conference, Warsaw, October 25–28, 1988, pp. 55-61], Williams [1996, Earth, Moon, Planets 72, 321–326; 2001, Proceedings of the Meteoroids 2001 conference, Kiruna, Sweden, August 6–10, 2001, pp. 33–42] estimated the ejection velocities of the stream meteoroids. Equation relating the ejection velocity Δυ and the change Δa of the semi-major axis, Williams (2001), was applied with two slightly different variations. In the first one (M1) as Δa the difference between the mean orbit of the stream and the orbit of the parent body was substituted, in the second one (M2), as Δa the dispersion of semi-major axes around the mean orbit of the stream was used. The results obtained by these two methods are not free from discrepancies, partly explained by the particular orbital structure of the stream. Kresak [1992, Contrib. Astron. Obs. Skalnate Pleso 22, 123–130] strongly criticized the attempts to determine the initial velocities of the stream using the statistics of the meteor orbits. He argued that this is essentially impossible, because the dispersion of the initial velocities are masked by much larger measuring errors and by the accumulated effects of planetary perturbations. In our paper, we study the reliability of M1 and M2 methods. We made a numerical experiment consisting of formation of several meteor streams and their dynamical evolution over 5000 years. We ejected meteoroids particles from the comets: 1P/Halley, 2P/Encke, 55P/Tempel-Tuttle, 109P/Swift-Tuttle and from minor planets (3200) Phaethon and 2002 SY50. During the integration, the ejection velocities were estimated using both M1 and M2 methods. The results show that the velocities obtained by M1 method are unstable: too high or too low, when compared with the known ejection velocities at the time of the stream formation. On the other hand, the velocities obtained using M2 method are too small, mostly. In principle, M2 estimates the dispersion of the distribution of the ejection velocities around the mean value, not the mean value itself. Applying more accurate Equation relating Δυ and Δa we decreased the bias of the results, but not their variation observed during the evolution of the streams and the parent bodies. We have found that the variability of the estimated ejection velocities was caused mainly by the gravitational changes of the semi-major axis and eccentricity of the parent body. In brief, we have found that the reliability of the results obtained by M1 or M2 method are low, and have to be used with great care.  相似文献   

9.
Field equations in the presence of a perfect fluid distribution are obtained in a scalar-tensor theory of gravitation proposed by Saez and Ballester (Phys. Lett. 113, 1985, 467) with the aid of Einstein–Rosen cylindrically symmetric metric. A static vacuum model and a non-static stiff fluid model are presented. The physical and geometrical properties of the stiff fluid model are studied.  相似文献   

10.
A new concept of three dimensional non-Keplerian trajectories with double angular momentum reversal is investigated with high performance solar sails. The main discussion of this paper is about such 3D solar inverse orbits with inner constraints. The problem is addressed in a time optimal control framework solved by an indirect method. Two typical solar inverse orbits have been achieved and presented in a 3D non-dimensional dynamic model in the Heliocentric Inertial Frame. Starting from the Earth orbit ecliptic plane, a sailcraft in the inverse orbit exhibits a butterfly shape trajectory. As such, the new orbits are symmetrical with respect to a plane which contains the Sun-perihelion line. The relation of the sail attitude angles between the two symmetrical parts of the orbits are used to reduce the simulation effort. The quasi-heliostationary property at its aphelia is demonstrated with variation of the orbital radius. Evolutions of the orbital velocity and optimal sail orientations are also outlined and discussed to benefit future design work. As is suited for space observation guaranteed by its butterfly shape, the inverse orbits are thoroughly studied in terms of the concerned parameters. The discussion of the parametric influence is ranked in order as perihelion distance r E , required maximum position z max, perihelion position z f and the sail lightness number β. Suitable ranges of each parameter are adopted to illustrate the orbital variation trend. Through numerical simulations the features of such inverse orbits are further emphasized to provide an initial reference for future researchers.  相似文献   

11.
We used a multipolar code to create, through dissipationless collapses of systems of 106 particles, two cuspy self-consistent triaxial stellar systems with γ ≈ 1. One of the systems has an axial ratio similar to that of an E4 galaxy and it is only mildly triaxial (T = 0.914), while the other one is strongly triaxial (T = 0.593) and its axial ratio lies in between those of Hubble types E5 and E6. Both models rotate although their total angular momenta are zero, i.e., they exhibit figure rotation. The angular velocity is very small for the less triaxial model and, while it is larger for the more triaxial one, it is still comparable to that found by Muzzio (Celest Mech Dynam Astron 96(2):85–97, 2006) to affect only slightly the dynamics of a similar model. Except for minor evolution, probably caused by unavoidable relaxation effects of the N-body code, the systems are highly stable. The potential of each system was subsequently approximated with interpolating formulae yielding smooth potentials, stationary in frames that rotate with the models. The Lyapunov exponents could then be computed for randomly selected samples of the bodies that make up the two systems, allowing the recognition of regular and of partially and fully chaotic orbits. Finally, the regular orbits were Fourier analyzed and classified using their locations on the frequency map. Most of the orbits are chaotic, and by a wide margin: less than 30% of the orbits are regular in our most triaxial model. Regular orbits are dominated by tubes, long axis ones in the less triaxial model and short axis tubes in the more triaxial one. Most of the boxes are resonant (i.e., they are boxlets), as could be expected from cuspy systems.  相似文献   

12.
On closed but non-geometrically similar orbits   总被引：1，自引：1，他引：0  
As is well known, the existence of geometrically similar orbits for a particle moving under a central conservative force is a consequence of the fact that the corresponding potential energy is a homogeneous function of the co-ordinates. In this paper, we consider a particular non-homogeneous potential of the form V = U + W, where U and W are homogenous functions of degrees −1 and −2, respectively, because, for this potential, the search of closed orbits, for discrete values of the angular momentum, is straightforward. We focus our attention on these daisy-like charming orbits and graphically show the consequences of the impossibility of geometrical similarity.  相似文献   

13.
We study the change of phase space structure of the rectilinear three-body problem when the mass combination is changed. Generally, periodic orbits bifurcate from the stable Schubart periodic orbit and move radially outward. Among these periodic orbits there are dominant periodic orbits having rotation number (n − 2)/n with n ≥ 3. We find that the number of dominant periodic orbits is two when n is odd and four when n is even. Dominant periodic orbits have large stable regions in and out of the stability region of the Schubart orbit (Schubart region), and so they determine the size of the Schubart region and influence the structure of the Poincaré section out of the Schubart region. Indeed, with the movement of the dominant periodic orbits, part of complicated structure of the Poincaré section follows these orbits. We find stable periodic orbits which do not bifurcate from the Schubart orbit.  相似文献   

14.
In the previous paper (Li et al. in Phys. Lett. B 666:125–130, 2008), we show the solutions of Einstein equations with static spherically-symmetric quintessence-like matter surrounding a global monopole. Furthermore, this monopole become a black hole with quintessence-like matter and a deficit solid angle when it is swallowed by an ordinary black hole. We study its quasinormal modes by WKB method in this paper. The numerical results show that both the real part of the quasinormal frequencies and the imaginary part decrease as the state parameter w, for scalar and gravitational perturbations. And we also show variations of quasinormal frequencies of scalar and gravitational fields via different ε (deficit solid angel parameter) and different ρ 0 (density of static spherically-symmetric quintessence-like matter at r=1), respectively.  相似文献   

15.
Summary Informations on 736 pairs of visual binaries are given in the form of the Tables. The General Catalogue gives ephemerides (t,θ°,ρ″) in 20 years (1984–2003) for each pair with the figures of its apparent elliptical orbit where the positions of secondary component relative to the primary one at different epochs are indicated.The General Catalogue contains four parts: Part one — Source and Grading of Orbit; Part two — Ephemérides and Atlas of Apparent Orbits; Part three — Classifications of Visual Binary Stars; Part four — Indexes of Visual Binary Stars. The principle of calculation and the statistical data are presented in this paper. There are seven statistical tables, giving the elements distribution of true and apparent orbits, grading distribution of orbits, number distribution with different physical property of component of binary star. The number of binary stars in anyone constellation, the number of binary stars brighter than 6m.5. The number of binary stars nearer than 25 parsec.  相似文献   

16.
We have studied periodic orbits generated by Lagrangian solutions of the restricted three body problem when one of the primaries is an oblate body. We have determined the periodic orbits for different values of μ, h and A (h is energy constant, μ is mass ratio of the two primaries and A is an oblateness factor). These orbits have been determined by giving displacements along the tangent and normal to the mobile coordinates as defined by Karimov and Sokolsky (Celest. Mech. 46:335, 1989). These orbits have been drawn by using the predictor-corrector method. We have also studied the effect of oblateness by taking some fixed values of μ, A and h. As starters for our method, we use some known periodic orbits in the classical restricted three body problem.  相似文献   

17.
H. Kiliç 《Solar physics》2009,255(1):155-162
The short-term periodicities in sunspot numbers, sunspot areas, and flare index data are investigated in detail using the Date Compensated Discrete Fourier Transform (DCDFT) for the full disk of the Sun separately over the rising, the maximum, and the declining portions of solar cycle 23 (1996 – 2006). While sunspot numbers and areas show several significant periodicities in a wide range between 23.1 and 36.4 days, the flare index data do not exhibit any significant periodicity. The earlier conclusion of Pap, Tobiska, and Bouwer (1990, Solar Phys. 129, 165) and Kane (2003, J. Atmos. Solar-Terr. Phys. 65, 1169), that the 27-day periodicity is more pronounced in the declining portion of a solar cycle than in the rising and maximum ones, seems to be true for sunspot numbers and sunspot area data analyzed here during solar cycle 23.  相似文献   

18.
A time-dependent model for the energy of a flaring solar active region is presented based on an existing stochastic jump-transition model (Wheatland and Glukhov in Astrophys. J. 494, 858, 1998; Wheatland in Astrophys. J. 679, 1621, 2008 and Solar Phys. 255, 211, 2009). The magnetic free energy of an active region is assumed to vary in time due to a prescribed (deterministic) rate of energy input and prescribed (random) jumps downwards in energy due to flares. The existing model reproduces observed flare statistics, in particular flare frequency – size and waiting-time distributions, but modeling presented to date has considered only the time-independent choices of constant energy input and constant flare-transition rates with a power-law distribution in energy. These choices may be appropriate for a solar active region producing a constant mean rate of flares. However, many solar active regions exhibit time variation in their flare productivity, as exemplified by NOAA active region (AR) 11029, observed during October – November 2009 (Wheatland in Astrophys. J. 710, 1324, 2010). Time variation is incorporated into the jump-transition model for two cases: (1) a step change in the rates of flare transitions, and (2) a step change in the rate of energy supply to the system. Analytic arguments are presented describing the qualitative behavior of the system in the two cases. In each case the system adjusts by shifting to a new stationary state over a relaxation time which is estimated analytically. The model exhibits flare-like event statistics. In each case the frequency – energy distribution is a power law for flare energies less than a time-dependent rollover set by the largest energy the system is likely to attain at a given time. The rollover is not observed if the mean free energy of the system is sufficiently large. For Case 1, the model exhibits a double exponential waiting-time distribution, corresponding to flaring at a constant mean rate during two intervals (before and after the step change), if the average energy of the system is large. For Case 2 the waiting-time distribution is a simple exponential, again provided the average energy of the system is large. Monte Carlo simulations of Case 1 are presented which confirm the estimate for the relaxation time and the expected forms of the frequency – energy and waiting-time distributions. The simulation results provide a qualitative model for observed flare statistics in AR 11029.  相似文献   

19.
In publications presenting analytical results on the non-coplanar motion of a circumbinary planet it was shown that the unperturbed elliptical orbit of the planet undergoes simultaneously two kinds of the precession: the precession of the orbital plane and the precession of the orbit in its own plane. It is also well-known that there is also the relativistic precession of the planetary orbit in its own plane. In the present paper we study a combined effect of the all of the above precessions. For the general case, where the planetary orbit is not coplanar with the stars orbits, we analyzed the dependence of the critical inclination angle ic, at which the precession of the planetary orbit in its own plane vanishes, on the angular momentum L of the planet. We showed that the larger the angular momentum, the smaller the critical inclination angle becomes. We presented the analytical result for ic(L) and calculated the value of L, for which the critical inclination value becomes zero. For the particular case, where the planetary orbit is not coplanar with the stars orbits, we demonstrated analytically that at a certain value of the angular momentum of the planet, the elliptical orbit of the planet would become stationary: no precession. In other words, at this value of the angular momentum, the relativistic precession of the planetary orbit and its precession, caused by the fact that the planet revolves around a binary (rather than single) star, cancel each other out. This is a counterintuitive result.  相似文献   

20.
In the second part of this papersee [V. Yu. Terebizh, Astrofizika,40, 139 (1997)] we give results of an auxiliary nature for a first-order autoregression process. A formal statement of the problem of estimating the spectral density of a time series is given as an inverse problem of mathematical physics. The numbering of sections in this and the next three papers in this series is a continuation of the numbering in the first paper, printed in the preceding issue of the journal [Astrofizika,40, 139–148 (1997)]. Translated from Astrofizika, Vol. 40, No. 2, pp. 273–280, April–June, 1 997.  相似文献   

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

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