共查询到20条相似文献，搜索用时 15 毫秒
1.
Kourosh Parand Mehran Nikarya Jamal Amani Rad 《Celestial Mechanics and Dynamical Astronomy》2013,116(1):97-107
The Lane–Emden type equations are employed in the modeling of several phenomena in the areas of mathematical physics and astrophysics. These equations are categorized as non-linear singular ordinary differential equations on the semi-infinite domain $[0,\infty )$ . In this research we introduce the Bessel orthogonal functions as new basis for spectral methods and also, present an efficient numerical algorithm based on them and collocation method for solving these well-known equations. We compare the obtained results with other results to verify the accuracy and efficiency of the presented scheme. To obtain the orthogonal Bessel functions we need their roots. We use the algorithm presented by Glaser et al. (SIAM J Sci Comput 29:1420–1438, 2007) to obtain the $N$ roots of Bessel functions. 相似文献
2.
H. B. Papo 《Earth, Moon, and Planets》1973,8(4):539-545
The differential equations of rotational motion of the Moon are solved by numerical integration methods. Euler's dynamical equations transformed to a convenient form are treated by techniques analogous to ordinary orbit determination procedures. The proposed method is fully consistent with the ephemeris of the Moon and can utilize a variety of observational material for the solution of the selected parameters. The parameters are grouped into three distinct groups, namely: --The physical libration angles of the Moon and their time rates at an arbitrary initial epoch. --Physical constants featuring the principal moments of ineria of the Moon. --Parameters associated with the particular observational material being used. Examples are given of comparison between the proposed method and Eckhardt's 1970 model of the physical librations of the Moon. The merits of the new method are discussed in the light of conventional data sources like Earth-based or satellite-based photography as well as newly available data types like Laser ranging to retroreflectors on the Moon. 相似文献
3.
Satellite orbits around a central body with arbitrary zonal harmonics are considered in a relativistic framework. Our starting point is the relativistic Celestial Mechanics based upon the first post-Newtonian approximation to Einstein’s theory of gravity as it has been formulated by Damour et al. (Phys Rev D 43:3273–3307, 1991; 45:1017–1044, 1992; 47:3124–3135, 1993; 49:618–635, 1994). Since effects of order \((\mathrm{GM}/c^2R) \times J_k\) with \(k \ge 2\) for the Earth are very small (of order \( 7 \times 10^{-10}\,\times \,J_k\)) we consider an axially symmetric body with arbitrary zonal harmonics and a static external gravitational field. In such a field the explicit \(J_k/c^2\)-terms (direct terms) in the equations of motion for the coordinate acceleration of a satellite are treated first with first-order perturbation theory. The derived perturbation theoretical results of first order have been checked by purely numerical integrations of the equations of motion. Additional terms of the same order result from the interaction of the Newtonian \(J_k\)-terms with the post-Newtonian Schwarzschild terms (relativistic terms related to the mass of the central body). These ‘mixed terms’ are treated by means of second-order perturbation theory based on the Lie-series method (Hori–Deprit method). Here we concentrate on the secular drifts of the ascending node \(<\!{\dot{\Omega }}\!>\) and argument of the pericenter \(<\!{\dot{\omega }}\!>\). Finally orders of magnitude are given and discussed. 相似文献
4.
Previously, we have considered the equations of motion of the three-body problem in a Lagrange form (which means a consideration of relative motions of 3-bodies in regard to each other). Analysing such a system of equations, we considered the case of small-body motion of negligible mass m 3 around the second of two giant-bodies m 1, m 2 (which are rotating around their common centre of masses on Kepler’s trajectories), the mass of which is assumed to be less than the mass of central body. In the current development, we have derived a key parameter η that determines the character of quasi-circular motion of the small third body m 3 relative to the second body m 2 (planet). Namely, by making several approximations in the equations of motion of the three-body problem, such the system could be reduced to the key governing Riccati-type ordinary differential equations. Under assumptions of R3BP (restricted three-body problem), we additionally note that Riccati-type ODEs above should have the invariant form if the key governing (dimensionless) parameter η remains in the range 10?2 Open image in new window 10?3. Such an amazing fact let us evaluate the forbidden zones for Moon’s orbits in the inner solar system or the zones of distances (between Moon and Planet) for which the motion of small body could be predicted to be unstable according to basic features of the solutions of Riccati-type. 相似文献
5.
The Tunguska event on 30 June 1908 has been subjected to much speculation within different fields of research. Publication of the results of the 1961 expedition to the Tunguska area (Florensky, 1963) supports that a cometary impact caused the event. Based on this interpretation, calculations of the impactor energy release and explosion height have been reported by Ben-Menahem (1975), and velocity, mass, and density of the impactor by Petrov and Stulov (1975). Park (1978) and Turco et al., 1981, Turco et al., 1982, used these numbers to calculate a production of ca. 30 × 106 tons of NO during atmospheric transit. This paper presents a high-resolution study of nitrate concentration in the Greenland ice sheet in ca. 10 years covering the Tunguska event. No signs of excess nitrate are found in three ice cores from two different sites in Greenland in the years following the Tunguska event. By comparing these results with results for other aerosols generally found in the ice, the lack of excess NO3? following the Tunguska event can be interpreted as indicating that the impactor nitrate production calculated by Park (1978) and Turco et al., 1981, Turco et al., 1982 are 1–2 orders of magnitude too high. To explain this it is suggested, from other lines of reasoning, that the impactor density determined by Petrov and Stulov (1975) probably is too low. 相似文献
6.
In this paper, we have investigated the plane symmetric space-time with wet dark fluid (WDF), which is a candidate for dark energy, in the framework of f (R,T) gravity Harko et al. 2011, Phys. Rev. D, 84, 024020), where R and T denote the Ricci scalar and the trace of the energy–momentum tensor respectively. We have used the equation of state in the form of WDF for the dark energy component of the Universe. It is modeled on the equation of state p = ω(ρ ? ρ ?). The exact solutions to the corresponding field equations are obtained for power-law and exponential volumetric expansion. The geometrical and physical parameters for both the models are studied. Also, we have discussed the well-known astrophysical phenomena, namely the look-back time, proper distance, the luminosity distance and angular diameter distance with red shift. 相似文献
7.
We study the physical behavior of a five dimensional non-static spherically symmetric cosmological models in the presence of massive strings in the framework of \(f(R,T)\) gravity proposed by Harko et al. (Phys. Rev. D 84:024020, 2011). Here \(R\) is the Ricci scalar and \(T\) is the trace of the stress energy tensor and the fifth dimension is not observed because it is compact. We solve the field equations (i) using a relation between the scale factors given by Samantha and Dhal (Int. J. Theor. Phys. 52:1334, 2013) and (ii) equations of state for string models. The models obtained correspond to \(p\)-string, geometric string and massive string models in this modified theory in five dimensions. Cosmological parameters of the models are determined and their dynamical properties are discussed. 相似文献
8.
Edgar Everhart 《Celestial Mechanics and Dynamical Astronomy》1974,10(1):35-55
The solutions of \(\ddot x = F(x,t)\) , and also \(\dot x = F(x,t)\) , are developed in truncated series in timet whose coefficients are found empirically. The series ending in thet 6 term yields a position at a final prechosen time that is accurate through 9th order in the sequence size. This is achieved by using Gauss-Radau and Gauss-Lobatto spacings for the several substeps within each sequence. This timeseries method is the same in principle as implicit Runge-Kutta forms, including some not described previously. In some orders these methods are unconditionally stable (A-stable). In the time-series formulation the implicit system converges rapidly. For integrating a test orbit the method is found to be about twice as fast as high-order explicit Runge-Kutta-Nyström-Fehlberg methods at the same accuracies. Both the Cowell and the Encke equations are solved for the test orbit, the latter being 35% faster. It is shown that the Encke equations are particularly well-adapted to treating close encounters when used with a single-sequence integrator (such as this one) provided that the reference orbit is re-initialized at the start of each sequence. This use of Encke equations is compared with the use of regularized Cowell equations. 相似文献
9.
Evolutionary tracks from the zero age main sequence to the asymptotic giant branch were computed for stars with initial masses 2 M ⊙ ≤ M ZAMS ≤ 5 M ⊙ and metallicity Z = 0.02. Some models of evolutionary sequences were used as initial conditions for equations of radiation hydrodynamics and turbulent convection describing radial stellar pulsations. The early asymptotic giant branch stars are shown to pulsate in the fundamental mode with periods 30 day ? Π ? 400day. The rate of period change gradually increases as the star evolves but is too small to be detected (Π?/Π < 10?5 yr?1). Pulsation properties of thermally pulsing AGB stars are investigated on time intervals comprising 17 thermal pulses for evolutionary sequences with initial masses M ZAMS = 2 M ⊙ and 3 M ⊙ and 6 thermal pulses for M ZAMS = 4 M ⊙ and 5 M ⊙. Stars with initial masses M ZAMS ≤ 3 M ⊙ pulsate either in the fundamental mode or in the first overtone, whereas more massive red giants (M ZAMS ≥ 4 M ⊙) pulsate in the fundamental mode with periods Π ? 103 day. Most rapid pulsation period change with rate ?0.02 yr?1 ? Π?/Π ? ?0.01 yr?1 occurs during decrease of the surface luminosity after the maximum of the luminosity in the helium shell source. The rate of subsequent increase of the period is Π?/Π ? 5 × 10?3 yr?1. 相似文献
10.
11.
12.
A. Elipe J. I. Montijano L. Rández M. Calvo 《Celestial Mechanics and Dynamical Astronomy》2017,129(4):415-432
In this note a study of the convergence properties of some starters \( E_0 = E_0(e,M)\) in the eccentricity–mean anomaly variables for solving the elliptic Kepler’s equation (KE) by Newton’s method is presented. By using a Wang Xinghua’s theorem (Xinghua in Math Comput 68(225):169–186, 1999) on best possible error bounds in the solution of nonlinear equations by Newton’s method, we obtain for each starter \( E_0(e,M)\) a set of values \( (e,M) \in [0, 1) \times [0, \pi ]\) that lead to the q-convergence in the sense that Newton’s sequence \( (E_n)_{n \ge 0}\) generated from \( E_0 = E_0(e,M)\) is well defined, converges to the exact solution \(E^* = E^*(e,M)\) of KE and further \( \vert E_n - E^* \vert \le q^{2^n -1}\; \vert E_0 - E^* \vert \) holds for all \( n \ge 0\). This study completes in some sense the results derived by Avendaño et al. (Celest Mech Dyn Astron 119:27–44, 2014) by using Smale’s \(\alpha \)-test with \(q=1/2\). Also since in KE the convergence rate of Newton’s method tends to zero as \( e \rightarrow 0\), we show that the error estimates given in the Wang Xinghua’s theorem for KE can also be used to determine sets of q-convergence with \( q = e^k \; \widetilde{q} \) for all \( e \in [0,1)\) and a fixed \( \widetilde{q} \le 1\). Some remarks on the use of this theorem to derive a priori estimates of the error \( \vert E_n - E^* \vert \) after n Kepler’s iterations are given. Finally, a posteriori bounds of this error that can be used to a dynamical estimation of the error are also obtained. 相似文献
13.
S. G. Zhuravlev 《Celestial Mechanics and Dynamical Astronomy》1981,25(3):297-315
In our article (Zhuravlev, 1979) a formal method of constructing conditionally periodic solutions of canonical systems of differential equations with a quick-rotating phase in the case of sharp commensurability was presented. The existence of stationary (or periodic) solutions of an averaged system of differential equations corresponding to the initial system of differential equations is necessary for an effective application of the method for different problems.Evidently, the stationary solutions do not always exist but in numerous papers on stationary solutions (oscillations or motions), the conditions of existence of such solutions are very often not considered at all. Usually a simple assumption is used that the stationary solutions do exist.Otherwise it is well known that Poincaré's theory of periodic solutions (Poincaré, 1892) let one set up conditions of existence of periodic solutions in different systems of differential equations. Particularly, in papers,Mah (1949, 1956), see alsoexmah (1971), the necessary and sufficient conditions of the existence of periodic solutions of (non-canonical) systems of differential equations which are close to arbitrary non-linear systems are given. For canonical autonomous systems of differential equations the conditions of existence of periodic solutions and a method of calculation are presented in the paperMepmah (1952).In our paper another approach is given and the conditions of existence of stationary solutions of canonical systems of differential equations with a quick-rotating phase are proved. For this purpose Delaunay-Zeipel's transformation and Poincaré's small parameter method are used. 相似文献
14.
E. A. Roth 《Celestial Mechanics and Dynamical Astronomy》1980,21(2):155-155
It is assumed that the dynamical system can be represented by equations of the form $$\begin{gathered} \dot x = \varepsilon _i f_i (x,y) \hfill \\ \dot y = u(x,y) + \varepsilon _i g_i (x,y) \hfill \\ \end{gathered} $$ as this is the case for the Lagrange equations in celestial mechanics. The perturbation functionsf i andg i may also depend on the timet. The fast angular variabley is now taken as independent variable. Using perturbation theory and expanding in Taylor series the differential equations for the zeroth, first, second, ... order approximations are obtained. In the stroboscopic method in particular the integration is performed analytically over one revolution, say from perigee to perigee. By the rectification step applied tox andt, the initial values for the next revolution are obtained. It is shown how the second order terms can be determined for the various perturbations occurring in satellite theory. The solution constructed in this way remains valid for thousands of revolutions. An important feature of the method is the small amount of computing time needed compared with numerical integration. 相似文献
15.
An improved version of the 3D stellar reddening map in a space with a radius of 1200 pc around the Sun and within 600 pc of the Galactic midplane is presented. As in the previous 2010 and 2012 versions of the map, photometry with an accuracy better than 0.05 m in the J and Ks bands for more than 70 million stars from the 2MASS catalogue is used in the new version. However, the data reduction technique is considerably more complicated. As before, an analysis of the distribution of stars near the main-sequence turnoff on the (J ? Ks)?Ks diagram, where they form a distribution maximum, provides a basis for the method. The shift of this maximum, i.e., the mode (J ? Ks), along (J ? Ks) and Ks, given the spatial variations of the mean dereddened color (J ? Ks)0 of these stars, is interpreted as a growth of the reddening with increasing distance. The main distinction of the new method is that instead of the fixed mean absolute magnitude, dereddened color, distance, and reddening for each cell, the individual values of these quantities are calculated for each star by iterations when solving the system of equations relating them. This has allowed one to increase the random accuracy of the map to 0.01 m and its spatial resolution to 20 pc in coordinates and distance and to 1° in longitude and latitude. Comparison with other reddening estimates for the same spatial cells and Gaia DR1 TGAS stars shows that the constructed map is one of the best maps for the space under consideration. Its systematic errors have been estimated to be σ(E(J ? Ks)) = 0.025 m , or σ(E(B ? V)) = 0.04 m . The main purpose of the map is to analyze the characteristics of Galactic structures, clouds, and cloud complexes. For this purpose, the reddening map within each spatial cell has also been computed by analyzing the reddening along each line of sight. 相似文献
16.
The analysis of relative motion of two spacecraft in Earth-bound orbits is usually carried out on the basis of simplifying assumptions. In particular, the reference spacecraft is assumed to follow a circular orbit, in which case the equations of relative motion are governed by the well-known Hill–Clohessy–Wiltshire equations. Circular motion is not, however, a solution when the Earth’s flattening is accounted for, except for equatorial orbits, where in any case the acceleration term is not Newtonian. Several attempts have been made to account for the \(J_2\) effects, either by ingeniously taking advantage of their differential effects, or by cleverly introducing ad-hoc terms in the equations of motion on the basis of geometrical analysis of the \(J_2\) perturbing effects. Analysis of relative motion about an unperturbed elliptical orbit is the next step in complexity. Relative motion about a \(J_2\)-perturbed elliptic reference trajectory is clearly a challenging problem, which has received little attention. All these problems are based on either the Hill–Clohessy–Wiltshire equations for circular reference motion, or the de Vries/Tschauner–Hempel equations for elliptical reference motion, which are both approximate versions of the exact equations of relative motion. The main difference between the exact and approximate forms of these equations consists in the expression for the angular velocity and the angular acceleration of the rotating reference frame with respect to an inertial reference frame. The rotating reference frame is invariably taken as the local orbital frame, i.e., the RTN frame generated by the radial, the transverse, and the normal directions along the primary spacecraft orbit. Some authors have tried to account for the non-constant nature of the angular velocity vector, but have limited their correction to a mean motion value consistent with the \(J_2\) perturbation terms. However, the angular velocity vector is also affected in direction, which causes precession of the node and the argument of perigee, i.e., of the entire orbital plane. Here we provide a derivation of the exact equations of relative motion by expressing the angular velocity of the RTN frame in terms of the state vector of the reference spacecraft. As such, these equations are completely general, in the sense that the orbit of the reference spacecraft need only be known through its ephemeris, and therefore subject to any force field whatever. It is also shown that these equations reduce to either the Hill–Clohessy–Wiltshire, or the Tschauner–Hempel equations, depending on the level of approximation. The explicit form of the equations of relative motion with respect to a \(J_2\)-perturbed reference orbit is also introduced. 相似文献
17.
Yu. A. Fadeyev 《Astronomy Letters》2012,38(4):260-270
Excitation of radial oscillations in population I (X = 0.7, Z = 0.02) red supergiants is investigated using the solution of the equations of radiation hydrodynamics and turbulent convection. The core helium burning stars with masses 8M ⊙ ≤ M ≤ 20M ⊙ and effective temperatures T eff < 4000 K are shown to be unstable against radial pulsations in the fundamental mode. The oscillation periods range between 45 and 1180 days. The pulsational instability is due to the κ-mechanism in the hydrogen and heliumionization zones. Radial pulsations of stars with mass M < 15M ⊙ are strictly periodic with the light amplitude ΔM bol ≤ 0?5. The pulsation amplitude increases with increasing stellar mass and for M > 15M ⊙ the maximum expansion velocity of outer layers is as high as one third of the escape velocity. The mean radii of outer Lagrangean mass zones increase due to nonlinear oscillations by ≤30% in comparison with the initial equilibrium. The approximate method (with uncertainty of a factor of 1.5) to evaluate the mass of the pulsating red supergiant with the known period of radial oscillations is proposed. The approximation of the pulsation constant Q as a function of the mass-to-radius ratio is given. Masses of seven galactic red supergiants are evaluated using the period-mean density relation. 相似文献
18.
In this paper, the initial value problem of space dynamics in universal Stumpff anomaly \(\psi\) is set up and developed in analytical and computational approach.For the analytical expansions, the linear independence of the functions \(\mathrm{U}_{{ j}} { (\psi;\varsigma)}\); \({j=0,1,2,3}\) are proved. The differential and recurrence equations satisfied by them and their relations with the elementary functions are given. The universal Kepler equation and its validations for different conic orbits are established together with the Lagrangian coefficients. Efficient representations of these functions are developed in terms of the continued fractions.
For the computational developments we consider the following items:
Finally we give summary on the computational design for the initial value problem of space dynamics in universal Stumpff anomaly. This design based on the solution of the universal Kepler’s equation by an iterative schemes of quadratic up to any desired order \(\ell\). 相似文献
- 1.Top-down algorithm for continued fraction evaluation.
- 2.One-point iteration formulae.
- 3.Determination of the coefficients of Kepler’s equation.
- 4.Derivatives of Kepler’s equation of any integer order.
- 5.Determination of the initial guess for the solution of the universal Kepler equation.
19.
Attitude stability of spacecraft subjected to the gravity gradient torque in a central gravity field has been one of the most fundamental problems in space engineering since the beginning of the space age. Over the last two decades, the interest in asteroid missions for scientific exploration and near-Earth object hazard mitigation is increasing. In this paper, the problem of attitude stability is generalized to a rigid spacecraft on a stationary orbit around a uniformly-rotating asteroid. This generalized problem is studied via the linearized equations of motion, in which the harmonic coefficients $C_{20}$ and $C_{22}$ of the gravity field of the asteroid are considered. The necessary conditions of stability of this conservative system are investigated in detail with respect to three important parameters of the asteroid, which include the harmonic coefficients $C_{20}$ and $C_{22}$ , as well as the ratio of the mean radius to the radius of the stationary orbit. We find that, due to the significantly non-spherical shape and the rapid rotation of the asteroid, the attitude stability domain is modified significantly in comparison with the classical stability domain predicted by the Beletskii–DeBra–Delp method on a circular orbit in a central gravity field. Especially, when the spacecraft is located on the intermediate-moment principal axis of the asteroid, the stability domain can be totally different from the classical stability domain. Our results are useful for the design of attitude control system in the future asteroid missions. 相似文献
20.
Based on the kinetic theory, Landau damping of dust acoustic waves (DAWs) propagating in a dusty plasma composed of hybrid nonthermal nonextensive distributed electrons, Maxwellian distributed ions and negatively charged dust grains is investigated using Vlasov-Poisson’s equations. The characteristics of the DAWs Landau damping are discussed. It is found that the wave frequency increases by decreasing (increasing) the value of nonextensive (nonthermal) parameter, \(q\) (\(\alpha \)). It is recognized that \(\alpha \) plays a significant role in observing damping or growing DAW oscillations. For small values of \(\alpha \), damping modes have been observed until reaching a certain value of \(\alpha \) at which \(\omega _{i}\) vanishes, then a growing mode appears in the case of superextensive electrons. However, only damping DAW modes are observed in case of subextensive electrons. The present study is useful in the space situations where such distribution exists. 相似文献