共查询到20条相似文献，搜索用时 15 毫秒
1.
The idea of using various L-matrices in numerical integration of the regular equations, which describe the motion of small bodies of the Solar System, is developed. The problem of the optimal position of the radius vector and velocity at numerical integration in the KS-coordinate system is posed. The solution of this problem, which reduces the number of calculations of the vector of perturbing accelerations, is given. The transformation providing this optimal solution is suggested, and the results of numerical integration are given. 相似文献
2.
The parameters of L matrices are applied to the numerical integration of regular equations describing the motion of minor bodies in the Solar System. The problem of the optimal choice of the regularizing change of variables is formulated in the context of the numerical integration of the equations of motion using the Runge–Kutta–Fehlberg method. Arbitrary perturbations are taken into account. This problem is completely solved in the case of planar motion. The solution of the optimization problem reduces the amount of computations needed to determine the vector of perturbing accelerations. Results of numerical integrations are given. 相似文献
3.
Guy Janin 《Celestial Mechanics and Dynamical Astronomy》1974,10(4):451-467
For computing highly eccentric (e0.9) Earth satellite orbits with special perturbation methods, a comparison is made between different schemes, namely the direct integration of the equations of motion in Cartesian coordinates, changes of the independent variable, use of a time element, stabilization and use of regular elements. A one-step and a multi-step integration are also compared.It is shown that stabilization and regularization procedures are very helpful for non or smoothly perturbed orbits. In practical cases for space research where all perturbations are considered, these procedures are no longer so efficient. The recommended method in these cases is a multi-step integration of the Cartesian coordinates with a change of the independent variable defining an analytical step size regulation. However, the use of a time element and a stabilization procedure for the equations of motion improves the accuracy, except when a small step size is chosen. 相似文献
4.
One of the main difficulties encountered in the numerical integration of the gravitationaln-body problem is associated with close approaches. The singularities of the differential equations of motion result in losses of accuracy and in considerable increase in computer time when any of the distances between the participating bodies decreases below a certain value. This value is larger than the distance when tidal effects become important, consequently,numerical problems are encounteredbefore the physical picture is changed. Elimination of these singularities by transformations is known as the process of regularization. This paper discusses such transformations and describes in considerable detail the numerical approaches to more accurate and faster integration. The basic ideas of smoothing and regularization are explained and applications are given. 相似文献
5.
An appropriate generalization of the Jacobi equation of motion for the polar moment of inertia I is considered in order to study the N-body problem with variable masses. Two coupled ordinary differential equations governing the evolution of I and the total energy E are obtained. A regularization scheme for this system of differential equations is provided. We compute some illustrative
numerical examples, and discuss an average method for obtaining approximate analytical solutions to this pair of equations.
For a particular law of mass loss we also obtain exact analytical solutions. The application of these ideas to other kind
of perturbed gravitational N-body systems involving drag forces or a different type of mass variation is also considered.
This revised version was published online in July 2006 with corrections to the Cover Date. 相似文献
6.
Johannes Hagel 《Celestial Mechanics and Dynamical Astronomy》1995,63(2):205-225
The problem of stability of the Lagrangian pointL
4 in the circular restricted problem of three bodies is investigated close to the 1 : 2 commensurability of the long and short period libration. By stability we define boundedness of the solution for a given initial finite displacement from the equilibrium point as function of the mass parameter close to the commensurability. A rigorous treatment close to the resonance condition is possible using a transformation that diagonalizes the matrix related to the linear part of the equations of motion. The so obtained equations are further transformed to action angle type variables. Then using an isolated resonance approach, only the slowly varying terms are kept in the equations and two independent isolating first integrals can be found. These integrals finally enable us to solve the stability problem in an exact way. The so obtained results are compared to numeric integration of the equations of motion and are found to be in perfect agreement. 相似文献
7.
This paper discusses the formulation and the numerical integration of large systems of differential equations occurring in the gravitational problem ofn-bodies.Different forms of the pertinent differential equations of motion are presented, and various regularizing and smoothing transformations are compared. Details regarding the effectiveness and the efficiency of the Kustaanheimo-Stiefel and of other methods are discussed. In particular, a method is described in which some of the phase variables are treated in the regularized system and others in the ordinary system. This mixed method of numerical regularization offers some advantages.Several numerical integration techniques are compared. A high order Runge-Kutta method, Steffensen's method, and a finite difference method are investigated, especially with regard to their adaptability to regularization.The role of integrals and integral invariants is displayed in controlling the accuracy of the numerical integration.Numerical results are described with 5, 25 and 500 bodies participating. These examples compare the various integration techniques, several regularization methods and different logics in treating binaries. 相似文献
8.
Trojan type orbits in the system of two gravitational centers with variable separation are studied within the framework of the restricted problem of three bodies. The backward numerical integration of the equations of motion of the bodies starting in the triangular libration pointsL
4 andL
5 (reverse problem) finds the breakdown of librations as the separation decreases because of the mass gain of the smaller component and an approach of the body of negligible, mass to the latter up to the distance below its sphere of action with a relative velocity approximately equal to the escape one on this sphere. The breakdown of librations aboutL
5 occurs under the mass gain of the smaller component considerably larger than in the case ofL
4 and implications are made for the asymmetry of the number of librators aboutL
4 andL
5 in the solar system (Greeks and Trojans). Other parameters of the libration motion near 1/1 commensurability are obtained, namely, the asymmetry of the libration amplitudes about the triangular points as well as the values of periods and amplitudes within the limits of those for real Trojan asteroids. Trojans could be supposedly, formed inside the Proto-jupiter and escape during its intensive mass loss. 相似文献
9.
T. A. Heppenheimer 《Celestial Mechanics and Dynamical Astronomy》1973,7(2):177-194
Out-of-plane motion about libration points is studied within the framework of the elliptic restricted three-body problem. Nonlinear motion in the circular restricted problem is given to third order in the out-of-plane amplitudeA
z by Jacobi elliptic functions. Linear motion in the elliptic problem is studied using Mathieu's and Hill's equations. Additional terms needed for a complete third-order theory are found using Lindsted's method. This theory is constructed for the case of collinear libration points; for the case of triangular points, a third-order nonlinear solution is given separately in terms of Jacobi elliptic functions. 相似文献
10.
R. K. Das A. K. Shrivastava B. Ishwar 《Celestial Mechanics and Dynamical Astronomy》1988,45(4):387-393
The differential equations of motion of the elliptic restricted problem of three bodies with decreasing mass are derived. The mass of the infinitesimal body varies with time. We have applied Jeans' law and the space-time transformation of Meshcherskii. In this problem the space-time transformation is applicable only in the special case whenn=1,k=0,q=1/2. We have applied Nechvile's transformation for the elliptic problem. We find that the equations of motion of our problem differ from that of constant mass only by a small perturbing force. 相似文献
11.
We consider the non-canonical Hamiltonian dynamics of a gyrostat in Newtonian interaction with n spherical rigid bodies. Using the symmetries of the system we carry out two reductions. Then, working in the reduced problem,
we obtain the equations of motion, a Casimir function of the system and the equations that determine the relative equilibria.
Global conditions for existence of relative equilibria are given. Besides, we give the variational characterization of these
equilibria and three invariant manifolds of the problem; being calculated the equations of motion in these manifolds, which
are described by means of a canonical Hamiltonian system. We give some Eulerian and Lagrangian equilibria for the four body
problem with a gyrostat. Finally, certain classical problems of Celestial Mechanics are generalized. 相似文献
12.
Byron D. Tapley 《Celestial Mechanics and Dynamical Astronomy》1970,2(3):319-333
A completely regular form for the differential equations governing the three-dimensional motion of a continuously thrusting space vehicle is obtained by using the Kustaanheimo-Stiefel regularization. The differential equations for the thrusting rocket are transformed using the K-S transformation and an optimal trajectory problem is posed in the transformed space. The canonical equations for the optimal motion in the transformed space are regularized by a suitable change of the independent variable. The transformed equations are regular in the sense that the differential equations do not possess terms with zero divisors when the motion encounters a gravitational force center. The resulting equations possess symmetry in form and the coefficients of the dependent variables are slowly varying quantities for a low-thrust space vehicle.Presented at the Conference on Celestial Mechanics, Oberwolfach, Germany, August 17–23, 1969. 相似文献
13.
The regularization of a new problem, namely the three-body problem, using ‘similar’ coordinate system is proposed. For this
purpose we use the relation of ‘similarity’, which has been introduced as an equivalence relation in a previous paper (see
Roman in Astrophys. Space Sci. doi:, 2011). First we write the Hamiltonian function, the equations of motion in canonical form, and then using a generating function,
we obtain the transformed equations of motion. After the coordinates transformations, we introduce the fictitious time, to
regularize the equations of motion. Explicit formulas are given for the regularization in the coordinate systems centered
in the more massive and the less massive star of the binary system. The ‘similar’ polar angle’s definition is introduced,
in order to analyze the regularization’s geometrical transformation. The effect of Levi-Civita’s transformation is described
in a geometrical manner. Using the resulted regularized equations, we analyze and compare these canonical equations numerically,
for the Earth-Moon binary system. 相似文献
14.
Numerical integration of unstable differential equations should be avoided since a numerical error during thenth step produces erroneous initial values for the next step and thus deteriorates the subsequent integration in an unstable manner. A method is offered to stabilize the equations of motion corresponding to a given HamiltonianH by transformingH into a new HamiltonianH
* which is equivalent to the Hamiltonian of a harmonic oscillator. In contrast to other methods of stabilization the realm of canonical mechanics is thus not abandoned. Perturbations are discussed and as examples the Keplerian motion and the motion of a gyroscope are presented. 相似文献
15.
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. 相似文献
16.
Victor R. Bond 《Celestial Mechanics and Dynamical Astronomy》1974,10(3):303-318
The Newtonian differential equations of motion for the two-body problem can be transformed into four, linear, harmonic oscillator equations by simultaneously applying the regularizing time transformation dt/ds=r and the Kustaanheimo-Stiefel (KS) coordinate transformation. The time transformation changes the independent variable from time to a new variables, and the KS transformation transforms the position and velocity vectors from Cartesian space into a four-dimensional space. This paper presents the derivation of uniform, regular equations for the perturbed twobody problem in the four-dimensional space. The variation of parameters technique is used to develop expressions for the derivatives of ten elements (which are constants in the unperturbed motion) for the general case that includes both perturbations which can arise from a potential and perturbations which cannot be derived from a potential. These element differential equations are slightly modified by introducing two additional elements for the time to further improve long term stability of numerical integration.Originally presented at the AAS/AIAA Astrodynamics Specialists Conference, Vail, Colorado, July 1973 相似文献
17.
D. V. Mikryukov 《Astronomy Letters》2018,44(5):337-350
A system of averaged equations of planetary motion around a central star is constructed. An astrocentric coordinate system is used. The two-planet problem is considered, but all constructions are easily generalized to an arbitrary number N of planets. The motion is investigated in modified (complex) Poincarécanonical elements. The averaging is performed by the Hori–Deprit method over the fast mean longitudes to the second order relative to the planetary masses. An expansion of the disturbing function is constructed using the Laplace coefficients. Some terms of the expansion of the disturbing function and the first terms of the expansion of the averaged Hamiltonian are given. The results of this paper can be used to investigate the evolution of orbits with moderate eccentricities and inclinations in various planetary systems. 相似文献
18.
This paper investigates the triangular libration points in the photogravitational restricted three-body problem of variable
mass, in which both the attracting bodies are radiating as well and the infinitesimal body vary its mass with time according
to Jeans’ law. Firstly, applying the space-time transformation of Meshcherskii in the special case when q=1/2, k=0, n=1, the differential equations of motion of the problem are given. Secondly, in analogy to corresponding problem with constant
mass, the positions of analogous triangular libration points are obtained, and the fact that these triangular libration points
cease to be classical ones when α≠0, but turn to classical L
4 and L
5 naturally when α=0 is pointed out. Lastly, introducing the space-time inverse transformation of Meshcherskii, the linear stability of triangular
libration points is tested when α>0. It is seen that the motion around the triangular libration points become unstable in general when the problem with constant
mass evolves into the problem with decreasing mass. 相似文献
19.
R. A. Krikorian 《Astrophysics》2003,46(4):496-501
It is shown that the equation of motion Du
j/Ds = (e/mc
2)F
ji
u
i , a natural generalization to the curved spacetime of the Heaviside-Lorentz law of ponderomotive force, is equivalent to the metric independent and covariant Van Dantzig's equations of motion dx
j [jpi] = 0 or L
v
p
i = 0, where p
i is the conjugate momentum 4-vector and v a vector determined by the condition p
i
v
i, only with respect to holonomic coordinates. With respect to an anholonomic system, the Heaviside-Lorentz equation is a particular case of the VD equations valid for a privileged class of anholonomic frames, those consisting of orthogonal unit vectors. 相似文献
20.
Victor A. Avdyushev 《Celestial Mechanics and Dynamical Astronomy》2003,87(4):383-409
Mainly, the author focuses on Baumgarte's method and its applications in satellite, asteroid, stellar and planetary problems. In the paper arguments are given for the use of energy relations for stabilization in the elliptical two-body problem. Stabilizing properties of Baumgarte's equations and others are discussed. A simple approach is proposed for stabilizing the equations of almost circular motion. By using Baumgarte's technique, the author derives stabilized equations of perturbed restricted three-body problem. It is shown experimentally that stabilization in the problems mentioned above can raise the accuracy of numerical integration by several orders. 相似文献