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

2.
We describe a collection of results obtained by numerical integration of orbits in the main problem of artificial satellite theory (theJ 2 problem). The periodic orbits have been classified according to their stability and the Poincaré surfaces of section computed for different values ofJ 2 andH (whereH is thez-component of angular momentum). The problem was scaled down to a fixed value (–1/2) of the energy constant. It is found that the pseudo-circular periodic solution plays a fundamental role. They are the equivalent of the Poincaré first-kind solutions in the three-body problem. The integration of the variational equations shows that these pseudo-circular solutions are stable, except in a very narrow band near the critical inclincation. This results in a sequence of bifurcations near the critical inclination, refining therefore some known results on the critical inclination, for instance by Izsak (1963), Jupp (1975, 1980) and Cushman (1983). We also verify that the double pitchfork bifurcation around the critical inclination exists for large values ofJ 2, as large as |J 2|=0.2. Other secondary (higher-order) bifurcations are also described. The equations of motion were integrated in rotating meridian coordinates.  相似文献   

3.
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.  相似文献   

4.
In a previous paper, we developed a technique for estimating the inner eccentricity in coplanar hierarchical triple systems on initially circular orbits, with comparable masses and with well-separated components, based on an expansion of the rate of change of the Runge-Lenz vector. Now, the same technique is extended to non-coplanar orbits. However, it can only be applied to systems with I 0 < 39.23° or I 0 > 140.77°, where I is the inclination of the two orbits, because of complications arising from the so-called ‘Kozai effect’. The theoretical model is tested against results from numerical integrations of the full equations of motion. This revised version was published online in July 2006 with corrections to the Cover Date.  相似文献   

5.
We study the problem of critical inclination orbits for artificial lunar satellites, when in the lunar potential we include, besides the Keplerian term, the J 2 and C 22 terms and lunar rotation. We show that, at the fixed points of the 1-D averaged Hamiltonian, the inclination and the argument of pericenter do not remain both constant at the same time, as is the case when only the J 2 term is taken into account. Instead, there exist quasi-critical solutions, for which the argument of pericenter librates around a constant value. These solutions are represented by smooth curves in phase space, which determine the dependence of the quasi-critical inclination on the initial nodal phase. The amplitude of libration of both argument of pericenter and inclination would be quite large for a non-rotating Moon, but is reduced to <0°.1 for both quantities, when a uniform rotation of the Moon is taken into account. The values of J 2, C 22 and the rotation rate strongly affect the quasi-critical inclination and the libration amplitude of the argument of pericenter. Examples for other celestial bodies are given, showing the dependence of the results on J 2, C 22 and rotation rate.  相似文献   

6.
In this work, we have simulated orbits of a particle moving in gravitational field of the Sun-Jupiter system. The effect of solar radiation pressure, including Poynting Robertson drag, on the evolution of particle orbits in phase space have been studied for different values of the parameter β 1 (the ratio of radiation to gravitational force) and initial conditions. Characteristics of various computed trajectories have been studied using wavelet transform (WT), Fourier transform (FT) and Poincare surface of section method. We use wavelet analysis to identify transitions of a trajectory in time-frequency plane and further apply it to classify it as regular or chaotic in phase space. Unlike the Fourier transform method (FT), we observe that the wavelet transform (WT) also provides a basis to identify ‘sticky’ trajectories in the present dynamical system.  相似文献   

7.
The spectra of ‘stretching numbers’ (or ‘local Lyapunov characteristic numbers’) are different in the ordered and in the chaotic domain. We follow the variation of the spectrum as we move from the centre of an island outwards until we reach the chaotic domain. As we move outwards the number of abrupt maxima in the spectrum increases. These maxima correspond to maxima or minima in the curve a(θ), where a is the stretching number, and θ the azimuthal angle. We explain the appearance of new maxima in the spectra of ordered orbits. The orbits just outside the last KAM curve are confined close to this curve for a long time (stickiness time) because of the existence of cantori surrounding the island, but eventually escape to the large chaotic domain further outside. The spectra of sticky orbits resemble those of the ordered orbits just inside the last KAM curve, but later these spectra tend to the invariant spectrum of the chaotic domain. The sticky spectra are invariant during the stickiness time. The stickiness time increases exponentially as we approach an island of stability, but very close to an island the increase is super exponential. The stickiness time varies substantially for nearby orbits; thus we define a probability of escape Pn(x) at time n for every point x. Only the average escape time in a not very small interval Δx around each x is reliable. Then we study the convergence of the spectra to the final, invariant spectrum. We define the number of iterations, N, needed to approach the final spectrum within a given accuracy. In the regular domain N is small, while in the chaotic domain it is large. In some ordered cases the convergence is anomalously slow. In these cases the maximum value of ak in the continued fraction expansion of the rotation number a = [a0,a1,... ak,...] is large. The ordered domain contains small higher order chaotic domains and higher order islands. These can be located by calculating orbits starting at various points along a line parallel to the q-axis. A monotonic variation of the sup {q}as a function of the initial condition q0 indicates ordered motions, a jump indicates the crossing of a localized chaotic domain, and a V-shaped structure indicates the crossing of an island. But sometimes the V-shaped structure disappears if the orbit is calculated over longer times. This is due to a near resonance of the rotation number, that is not followed by stable islands. This revised version was published online in July 2006 with corrections to the Cover Date.  相似文献   

8.
Linked twist maps of ergodic theory are modified to provide a model for chaotic orbits in ‘double bowl’ potentials in conservative systems, e.g. direct ‘exchange’ orbits in the planar restricted problem of three bodies. The maps are typically ergodic (and mixing), and provide a clue to the nature of relaxation or diffusion in chaotic Hamiltonian systems.  相似文献   

9.
We investigate on the relative inclination of the planets B and C orbiting the pulsar PSR B1257+12. First, we show that the third Kepler’s law does represent an adequate model for the orbital periods P of the planets, because other Newtonian and Einsteinian corrections are orders of magnitude smaller than the accuracy in measuring P B/C. Then, on the basis of available timing data, we determine the ratio sin i C/ sin i B = 0.92±0.05 of the orbital inclinations i B and i C independently of the pulsar’s mass M. It turns out that coplanarity of the orbits of B and C would imply a violation of the equivalence principle. Adopting a pulsar mass range 1 ≲ M ≲ 3, in solar masses (supported by present-day theoretical and observational bounds for pulsar’s masses), both face-on and edge-on orbital configurations for the orbits of the two planets are ruled out; the acceptable inclinations for B span the range 36 deg ≲ i B ≲ 66 deg, with a corresponding relative inclination range 6 deg ≲ (i C − i B) ≲ 13 deg.  相似文献   

10.
The critical inclination is of special interest in artificial satellite theory. The critical inclination can maintain minimal deviations of eccentricity and argument of pericentre from the initial values, and orbits at this inclination have been applied to some space missions. Most previous researches about the critical inclination were made under the assumption that the oblateness term J 2 is dominant among the harmonic coefficients. This paper investigates the extension of the critical inclination where the concept of the critical inclination is different from that of the traditional sense. First, the study takes the case of Venus for instance, and provides some preliminary results. Then for general cases, given the values of argument of pericentre and eccentricity, the relationship between the multiplicity of the solutions for the critical inclination and the values of J 2 and J 4 is analyzed. Besides, when given certain values of J 2 and J 4, the relationship between the multiplicity of the solutions for the critical inclination and the values of semimajor axis and eccentricity is studied. The results show that for some cases, the value of the critical inclination is far away from that of the traditional sense or even has multiple solutions. The analysis in this paper could be used as starters of correction methods in the full gravity field of celestial bodies.  相似文献   

11.
The vertical stability character of the families of short and long period solutions around the triangular equilibrium points of the restricted three-body problem is examined. For three values of the mass parameter less than equal to the critical value of Routh (μ R ) i.e. for μ = 0.000953875 (Sun-Jupiter), μ = 0.01215 (Earth-Moon) and μ = μ R = 0.038521, it is found that all such solutions are vertically stable. For μ > (μ R ) vertical stability is studied for a number of ‘limiting’ orbits extended to μ = 0.45. The last limiting orbit computed by Deprit for μ = 0.044 is continued to a family of periodic orbits into which the well known families of long and short period solutions merge. The stability characteristics of this family are also studied.  相似文献   

12.
We study the existence, linear stability and bifurcations of what we call the Sitnikov family of straight line periodic orbits in the case of the restricted four-body problem, where the three equal mass primary bodies are rotating on a circle and the fourth (small body) is moving in the direction vertical to the center mass of the other three. In contrast to the restricted three-body Sitnikov problem, where the Sitnikov family has infinitely many stability intervals (hence infinitely many Sitnikov critical orbits), as the “family parameter” ż0 varies within a finite interval (while z 0 tends to infinity), in the four-body problem this family has only one stability interval and only twelve 3-dimensional (3D) families of symmetric periodic orbits exist which bifurcate from twelve corresponding critical Sitnikov periodic orbits. We also calculate the evolution of the characteristic curves of these 3D branch-families and determine their stability. More importantly, we study the phase space dynamics in the vicinity of these orbits in two ways: First, we use the SALI index to investigate the extent of bounded motion of the small particle off the z-axis along its interval of stable Sitnikov orbits, and secondly, through suitably chosen Poincaré maps, we chart the motion near one of the 3D families of plane-symmetric periodic orbits. Our study reveals in both cases a fascinating structure of ordered motion surrounded by “sticky” and chaotic orbits as well as orbits which rapidly escape to infinity.  相似文献   

13.
Finding relative satellite orbits that guarantee long-term bounded relative motion is important for cluster flight, wherein a group of satellites remain within bounded distances while applying very few formationkeeping maneuvers. However, most existing astrodynamical approaches utilize mean orbital elements for detecting bounded relative orbits, and therefore cannot guarantee long-term boundedness under realistic gravitational models. The main purpose of the present paper is to develop analytical methods for designing long-term bounded relative orbits under high-order gravitational perturbations. The key underlying observation is that in the presence of arbitrarily high-order even zonal harmonics perturbations, the dynamics are superintegrable for equatorial orbits. When only J 2 is considered, the current paper offers a closed-form solution for the relative motion in the equatorial plane using elliptic integrals. Moreover, necessary and sufficient periodicity conditions for the relative motion are determined. The proposed methodology for the J 2-perturbed relative motion is then extended to non-equatorial orbits and to the case of any high-order even zonal harmonics (J 2n , n ≥ 1). Numerical simulations show how the suggested methodology can be implemented for designing bounded relative quasiperiodic orbits in the presence of the complete zonal part of the gravitational potential.  相似文献   

14.
Frozen orbits are always important foci of orbit design because of their valuable characteristics that their eccentricity and argument of pericentre remain constant on average. This study investigates quasi-circular frozen orbits and examines their basic nature analytically using two different methods. First, an analytical method based on Lagrangian formulations is applied to obtain constraint conditions for Martian frozen orbits. Second, Lie transforms are employed to locate these orbits accurately, and draw the contours of the Hamiltonian to show evolutions of the equilibria. Both methods are verified by numerical integrations in an 80 × 80 Mars gravity field. The simulations demonstrate that these two analytical methods can provide accurate enough results. By comparison, the two methods are found well consistent with each other, and both discover four families of Martian frozen orbits: three families with small eccentricities and one family near the critical inclination. The results also show some valuable conclusions: for the majority of Martian frozen orbits, argument of pericentre is kept at 270° because J 3 has the same sign as J 2; while for a minority of ones with low altitude and low inclination, argument of pericentre can be kept at 90° because of the effect of the higher degree odd zonals; for the critical inclination cases, argument of pericentre can also be kept at 90°. It is worthwhile to note that there exist some special frozen orbits with extremely small eccentricity, which could provide much convenience for reconnaissance. Finally, the stability of Martian frozen orbits is estimated based on the trace of the monodromy matrix. The analytical investigations can provide good initial conditions for numerical correction methods in the more complex models.  相似文献   

15.
For the n-centre problem of one particle moving in the potential of attracting centres of small mass fixed in an arbitrary smooth potential and magnetic field, we prove the existence of periodic and chaotic trajectories shadowing sequences of collision orbits. In particular, we obtain large subshifts of solutions of this type for the circular restricted 3-body problem of celestial mechanics. Poincaré had conjectured existence of the periodic ones and given them the name ‘second species solutions’. This revised version was published online in July 2006 with corrections to the Cover Date.  相似文献   

16.
In this concise review of the recent developments in relativistic shock theory in the Universe we restrict ourselves to shocks that do not exhibit quantum effects. On the other hand, emphasis is given to the formation of shocks under both non-magnetised and magnetised conditions. We only briefly discuss particle acceleration in relativistic shocks where much of the results are still preliminary. Analytical theory is rather limited in predicting the real shock structure. Kinetic instability theory is briefed including its predictions and limitations. A recent self-similar relativistic shock theory is described which predicts the average long-term shock behaviour to be magnetised and to cause reasonable power-law distributions for energetic particles. The main focus in this review is on numerical experiments on highly relativistic shocks in (i) pair and (ii) electron-nucleon plasmas and their limitations. These simulations do not validate all predictions of analytic and self-similar theory and so far they do not solve the injection problem and the self-modification by self-generated cosmic rays. The main results of the numerical experiments discussed in this review are: (i) a confirmation of shock evolution in non-magnetised relativistic plasma in 3D due to either the lepton-Weibel instability (in pair plasmas) or to the ion-Weibel instability; (ii) the sensitive dependence of shock formation on upstream magnetisation which causes suppression of Weibel modes for large upstream magnetisation ratios σ>10−3; (iii) the sensitive dependence of particle dynamics on the upstream magnetic inclination angle θ Bn , where particles of θ Bn >34° cannot escape upstream, leading to the distinction between ‘subluminal’ and ‘superluminal’ shocks; (iv) particles in ultra-relativistic shocks can hardly overturn the shock and escape to upstream; they may oscillate around the shock ramp for a long time, so to speak ‘surfing it’ and thereby becoming accelerated by a kind of SDA; (v) these particles form a power-law tail on the downstream distribution; their limitations are pointed out; (vi) recently developed methods permit the calculation of the radiation spectra emitted by the downstream high-energy particles; (vii) the Weibel-generated downstream magnetic fields form large-amplitude vortices which could be advected by the downstream flow to large distances from the shock and possibly contribute to an extended strong field region; (viii) if cosmic rays are included, Bell-like modes can generate upstream magnetic turbulence at short and, by diffusive re-coupling, also long wavelengths in nearly parallel magnetic field shocks; (ix) advection of such large-amplitude waves should cause periodic reformation of the quasi-parallel shock and eject large-amplitude magnetic field vortices downstream where they contribute to turbulence and to maintaining an extended region of large magnetic fields.  相似文献   

17.
It is surprising that we hardly know only 4% of the universe. Rest of the universe is made up of 73% of dark-energy and 23% of dark-matter. Dark-energy is responsible for acceleration of the expanding universe; whereas dark-matter is said to be necessary as extra-mass of bizarre-properties to explain the anomalous rotational-velocity of galaxy. Though the existence of dark-energy has gradually been accepted in scientific community, but the candidates for dark-matter have not been found as yet and are too crazy to be accepted. Thus, it is obvious to look for an alternative theory in place of dark-matter. Milgrom (Astrophys. J. 270:365, 1983a; 270:371, 1983b) has suggested a ‘Modified Newtonian Dynamics (MOND)’ which appears to be highly successful for explaining the anomalous rotational-velocity. But unfortunately MOND lacks theoretical support. The MOND, in-fact, is (empirical) modification of Newtonian-Dynamics through modification in the kinematical acceleration term ‘a’ (which is normally taken as a=\fracv2ra=\frac{v^{2}}{r}) as effective kinematic acceleration aeffective = a m(\fracaa0)a_{\mathit{effective}} = a \mu(\frac{a}{a_{0}}), wherein the μ-function is 1 for usual-values of accelerations but equals to \fracaa0 ( << 1)\frac{a}{a_{0}} (\ll1) if the acceleration ‘a’ is extremely-low lower than a critical value a 0(10−10 m/s2). In the present paper, a novel variant of MOND is proposed with theoretical backing; wherein with the consideration of universe’s acceleration a d due to dark-energy, a new type of μ-function on theoretical-basis emerges out leading to aeffective = a(1 -K \fraca0a)a_{\mathit{effective}} = a(1 -K \frac{a_{0}}{a}). The proposed theoretical-MOND model too is able to fairly explain ‘qualitatively’ the more-or-less ‘flat’ velocity-curve of galaxy-rotation, and is also able to predict a dip (minimum) on the curve.  相似文献   

18.
Résumé On développe une méthode de construction d'orbites périoldiques dans un système d'axes tournants, pour un satellite gravitant autour d'un sphéroide. Les orbites sont quasi circulaires,i est l'inclinaison sur le plan équatorial de la planète. Pour les petites inclinaisons, la solution est donnée jusqu'aux termes enJ 2 2 etJ 4.Ce modèle peut être appliqué aux satellites de Saturne. Des valeurs observées des longitudes des noeuds ascendants de Mimas et Téthys, on donne une estimation des valeurs deJ 2 etJ 4 du potentiel de Saturne. La valeur deJ 2 est très sensible aux valeurs adoptées pour le rayon équatorial de la planète.
Construction of periodic orbits of satellites in a moving system of axes, I
We give an algorithm for the construction of periodic orbits in a rotating frame for the cases of satellites moving around an oblate planet.The orbits are near to the circular case; the asymptotic developments of the periodic solutions are completely calculated for the termsJ 2 andJ 4 of the potential. The solutions for small inclinations are given up toJ 2 2 .The families of solutions depend on three parameters: the semi-major axis, the inclination of the generating orbit and the initial position on this orbit.These solutions can be applied to the motion of the Saturnian satellites. From the observed longitudes of the ascending nodes of Mimas and Tethys, we estimate the valuesJ 2 andJ 4 of the Saturnian potential, the value ofJ 2 very strongly depends on the adopted value of the planet's equatorial diameter.
  相似文献   

19.
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.  相似文献   

20.
Sedimentation of particles in a fluid has long been used to characterize particle size distribution. Stokes’ law is used to determine an unknown distribution of spherical particle sizes by measuring the time required for the particles to settle a known distance in a fluid of known viscosity and density. In this paper, we study the effects of gravity on sedimentation by examining the resulting particle concentration distributed in an equilibrium profile of concentration C m,n above the bottom of a container. This is for an experiment on the surface of the Earth and therefore the acceleration of gravity had been corrected for the oblateness of the Earth and its rotation. Next, at the orbital altitude of the spacecraft in orbit around Earth the acceleration due to the central field is corrected for the oblateness of the Earth. Our results show that for experiments taking place in circular or elliptical orbits of various inclinations around the Earth the concentration ratio C m,n /C m,ave , the inclination seems to be the most ineffective in affecting the concentration among all the orbital elements. For orbital experiment that use particles of diameter d p =0.001 μm the concentration ratios for circular and slightly elliptical orbits in the range e=0–0.1 exhibit a 0.009 % difference. The concentration ratio increases with the increase of eccentricity, which increases more for particles of larger diameters. Finally, for particles of the same diameter concentration ratios between Earth and Mars surface experiments are related in the following way .  相似文献   

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

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