共查询到20条相似文献，搜索用时 31 毫秒
1.
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. 相似文献
2.
A regularization method for integrating the equations of motion of small N-body systems is discussed. We select a chain of interparticle vectors in such a way that the critical interactions requiring regularization are included in the chain. The equations of motion for the chain vectors are subsequently regularized using the KS-variables and a time transformation. The method has been formulated for any number of bodies, but the most important application appears to be in the four-body problem which is therefore discussed in detail. 相似文献
3.
Sergei M. Poleshchikov 《Celestial Mechanics and Dynamical Astronomy》2003,85(4):341-393
The sets of L-matrices of the second, fourth and eighth orders are constructed axiomatically. The defining relations are taken from the regularization of motion equations for Keplerian problem. In particular, the Levi-Civita matrix and KS-matrix are L-matrices of second and fourth order, respectively. A theorem on the ranks of L-transformations of different orders is proved. The notion of L-similarity transformation is introduced, certain sets of L-matrices are constructed, and their classification is given. An application of fourth order L-matrices for N-body problem regularization is given. A method of correction for regular coordinates in the Runge–Kutta–Fehlberg integration method for regular motion equations of a perturbed two-body problem is suggested. Comparison is given for the results of numerical integration in the problem of defining the orbit of a satellite, with and without the above correction method. The comparison is carried out with respect to the number of calls to the subroutine evaluating the perturbational accelerations vector. The results of integration using the correction turn out to be in a favorable position. 相似文献
4.
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. 相似文献
5.
Time transformations involving momenta in addition to the coordinates are studied from the points of view of stabilization and regularization of the equations of motion. The generalization of Sundman's transformation by using the potential function to transform the time is further generalized by using the Lagrangian function for the same purpose. The possibility of the stabilization of the equations of motion is investigated similarly to Stiefel's and Baumgarte's recent results but instead of a factorial, an additive control function is introduced in all equations of motion. The relation between the original and new independent variables is integrated by a modification of Ebert's theorem and it is shown that the new independent variable is Hamilton's principal function. Numerical examples illustrate the method and seem to indicate that the computation of close approach trajectories benefit especially by the transformations discussed. The Appendix offers an analytic treatment regarding the stabilization of the constant of energy. 相似文献
6.
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. 相似文献
7.
The chain regularization method (Mikkola and Aarseth 1990) for high accuracy computation of particle motions in smallN-body systems has been reformulated. We discuss the transformation formulae, equations of motion and selection of a chain of interparticle vectors such that the critical interactions requiring regularization are included in the chain. The Kustaanheimo-Stiefel (KS) coordinate transformation and a time transformation is used to regularize the dominant terms of the equations of motion. The method has been implemented for an arbitrary number of bodies, with the option of external perturbations. This formulation has been succesfully tested in a generalN-body program for strongly interacting subsystems. An easy to use computer program, written inFortran, is available on request. 相似文献
8.
Szebehely's renowned equation given in 1974, allowing for potential determination from a given orbit or family of orbits, is proved to be equivalent with an equation deduced in 1963 by Drǎmbǎ. This basic equation in the inverse problem of dynamics, for which the denomination of Drǎmbǎ –Szebehely equation is proposed, is generalized for the motion in the n-dimensional Euclidean space. A method for the determination of the potential function from motion equations is extended to this space. 相似文献
9.
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. 相似文献
10.
Jörg Waldvogel 《Celestial Mechanics and Dynamical Astronomy》2008,102(1-3):149-162
Quaternions have been found to be the ideal tool for describing and developing the theory of spatial regularization in Celestial Mechanics. This article corroborates the above statement. Beginning with a summary of quaternion algebra, we will describe the regularization procedure and its consequences in an elegant way. Also, an alternative derivation of the theory of Kepler motion based on regularization will be given. Furthermore, we will consider the regularization of the spatial restricted three-body problem, i.e. the spatial generalization of the Birkhoff transformation. Finally, the perturbed Kepler motion will be described in terms of regularized variables. 相似文献
11.
A. G. Mavraganis 《Celestial Mechanics and Dynamical Astronomy》1987,42(1-4):169-174
The author's aim is to achieve global regularization in the Magnetic-Binary problem by suitably transforming the state-time space of the system. The functions which perform the change of the physical time and the geometrical figures of the system, are connected by a special relation leaving the form of the equations of motion invariant. Additionally, a proposition for generalization of the process is discussed in an aspect as well, of how much such a regularization is profitable. 相似文献
12.
Using the rectangular equations of motion for the restricted three-body problem a comparison is made of the integration of these equations by the Encke method and by a set of perturbational equations. Each set of differential equations is integrated using Taylor series expansions where the coefficients of the powers of time are determined by recurrence relations. It is shown that for very small perturbations the use of the perturbational equations is more efficient than the use of the Encke method. A discussion is also given of when Cowell's method is more efficient than either of these techniques. 相似文献
13.
Arnold Stokes 《Celestial Mechanics and Dynamical Astronomy》1977,16(1):27-34
A regularization of Kepler's problem due to Moser is used to ‘stabilize’ the equations of motion, that is, imbed a particular solution of Kepler's problem in a Lyapounov stable system. 相似文献
14.
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. 相似文献
15.
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. 相似文献
16.
Joachim W. Baumgarte 《Celestial Mechanics and Dynamical Astronomy》1976,13(2):247-251
In order to reduce the error growth during a numerical integration, a method of stabilization, of the differential equations of the Keplerian motion is offered. It is characterized by the use of the eccentric anomaly as independent variable in such a way that the time transformation is given by a generalized Lagrange formalism. The control terms in the equations of motion obtained by this modified Lagrangian give immediately a completely Lyapunov-stable set of differential equations. In contrast to other publications, here the equation of time integration is modified by a control term which leads to an integral which defined the time element for the perturbed Keplerian motion.This paper was supported by the National Research Council and the National Aeronautics and Space Administration and also by the Deutsche Forschungsgemeinschaft. It was presented at the Flight Mechanics/Estimation Theory Symposium, Goddard Space Flight Center, Greenbelt, Md., April 15–16, 1975. 相似文献
17.
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. 相似文献
18.
Taeyoung Lee Melvin Leok N. Harris McClamroch 《Celestial Mechanics and Dynamical Astronomy》2007,98(2):121-144
Equations of motion, referred to as full body models, are developed to describe the dynamics of rigid bodies acting under
their mutual gravitational potential. Continuous equations of motion and discrete equations of motion are derived using Hamilton’s
principle. These equations are expressed in an inertial frame and in relative coordinates. The discrete equations of motion,
referred to as a Lie group variational integrator, provide a geometrically exact and numerically efficient computational method
for simulating full body dynamics in orbital mechanics; they are symplectic and momentum preserving, and they exhibit good
energy behavior for exponentially long time periods. They are also efficient in only requiring a single evaluation of the
gravity forces and moments per time step. The Lie group variational integrator also preserves the group structure without
the use of local charts, reprojection, or constraints. Computational results are given for the dynamics of two rigid dumbbell
bodies acting under their mutual gravity; these computational results demonstrate the superiority of the Lie group variational
integrator compared with integrators that are not symplectic or do not preserve the Lie group structure. 相似文献
19.
Murray Silver 《Celestial Mechanics and Dynamical Astronomy》1975,11(1):39-41
A short derivation is given of the regularized equations of motion for the perturbed two-body problem. This method is then applied to the slightly modified time transformation dt/ds=r/ω. 相似文献
20.
N. Borderies 《Celestial Mechanics and Dynamical Astronomy》1977,16(3):291-308
This paper deals with the Adams-Moulton-Cowell multistep integrator, as described by Oestwinter and Cohen (1972). In order to evaluate the accuracy of the method, we started to test it in the case of the unperturbed two-body motion; numerical instability may arise by integrating first order systems. The accuracy is improved by applying a Sundmann transformation of the independent variable. The algorithm is then modified such that the equations of pure keplerian motion are integrated with respect to the new independent variable without truncation error; numerical experiments show the considerable improvement of accuracy and the reduction of computing time for Keplerian motion.If terms of the disturbing function of the Earth are added to the central potential, the time-transformation is less effective. With a modification of this time-transformation as given by Moynot in 1971, it is possible to reduce the propagation of the truncation error in the J2 problem. 相似文献