首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Small tidal forces in the Earth–Moon system cause detectable changes in the orbit. Tidal energy dissipation causes secular rates in the lunar mean motion n, semimajor axis a, and eccentricity e. Terrestrial dissipation causes most of the tidal change in n and a, but lunar dissipation decreases eccentricity rate. Terrestrial tidal dissipation also slows the rotation of the Earth and increases obliquity. A tidal acceleration model is used for integration of the lunar orbit. Analysis of lunar laser ranging (LLR) data provides two or three terrestrial and two lunar dissipation parameters. Additional parameters come from geophysical knowledge of terrestrial tides. When those parameters are converted to secular rates for orbit elements, one obtains dn/dt = \(-25.97\pm 0.05 ''/\)cent\(^{2}\), da/dt = 38.30 ± 0.08 mm/year, and di/dt = ?0.5 ± 0.1 \(\upmu \)as/year. Solving for two terrestrial time delays and an extra de/dt from unspecified causes gives \(\sim \) \(3\times 10^{-12}\)/year for the latter; solving for three LLR tidal time delays without the extra de/dt gives a larger phase lag of the N2 tide so that total de/dt = \((1.50 \pm 0.10)\times 10^{-11}\)/year. For total dn/dt, there is \(\le \)1 % difference between geophysical models of average tidal dissipation in oceans and solid Earth and LLR results, and most of that difference comes from diurnal tides. The geophysical model predicts that tidal deceleration of Earth rotation is \(-1316 ''\)/cent\(^{2}\) or 87.5 s/cent\(^{2}\) for UT1-AT, a 2.395 ms/cent increase in the length of day, and an obliquity rate of 9 \(\upmu \)as/year. For evolution during past times of slow recession, the eccentricity rate can be negative.  相似文献   

2.
This paper is devoted to the study of the transfer problem from a libration point orbit of the Earth–Moon system to an orbit around the Moon. The transfer procedure analysed has two legs: the first one is an orbit of the unstable manifold of the libration orbit and the second one is a transfer orbit between a certain point on the manifold and the final lunar orbit. There are only two manoeuvres involved in the method and they are applied at the beginning and at the end of the second leg. Although the numerical results given in this paper correspond to transfers between halo orbits around the \(L_1\) point (of several amplitudes) and lunar polar orbits with altitudes varying between 100 and 500 km, the procedure we develop can be applied to any kind of lunar orbits, libration orbits around the \(L_1\) or \(L_2\) points of the Earth–Moon system, or to other similar cases with different values of the mass ratio.  相似文献   

3.
The aim of this work is to combine the model of orbital and rotational motion of the Moon developed for DE430 with up-to-date astronomical, geodynamical, and geo- and selenophysical models. The parameters of the orbit and physical libration are determined in this work from lunar laser ranging (LLR) observations made at different observatories in 1970–2013. Parameters of other models are taken from solutions that were obtained independently from LLR. A new implementation of the DE430 lunar model, including the liquid core equations, was done within the EPM ephemeris. The postfit residuals of LLR observations make evident that the terrestrial models and solutions recommended by the IERS Conventions are compatible with the lunar theory. That includes: EGM2008 gravitational potential with conventional corrections and variations from solid and ocean tides; displacement of stations due to solid and ocean loading tides; and precession-nutation model. Usage of these models in the solution for LLR observations has allowed us to reduce the number of parameters to be fit. The fixed model of tidal variations of the geopotential has resulted in a lesser value of Moon’s extra eccentricity rate, as compared to the original DE430 model with two fit parameters. A mixed model of lunar gravitational potential was used, with some coefficients determined from LLR observations, and other taken from the GL660b solution obtained from the GRAIL spacecraft mission. Solutions obtain accurate positions for the ranging stations and the five retroreflectors. Station motion is derived for sites with long data spans. Dissipation is detected at the lunar fluid core-solid mantle boundary demonstrating that a fluid core is present. Tidal dissipation is strong at both Earth and Moon. Consequently, the lunar semimajor axis is expanding by 38.20 mm/yr, the tidal acceleration in mean longitude is \(-25.90 {{}^{\prime \prime }}/\mathrm{cy}^2\), and the eccentricity is increasing by \(1.48\times 10^{-11}\) each year.  相似文献   

4.
Lunar frozen orbits, characterized by constant orbital elements on average, have been previously found using various dynamical models, incorporating the gravitational field of the Moon and the third-body perturbation exerted by the Earth. The resulting mean orbital elements must be converted to osculating elements to initialize the orbiter position and velocity in the lunar frame. Thus far, however, there has not been an explicit transformation from mean to osculating elements, which includes the zonal harmonic \(J_2\), the sectorial harmonic \(C_{22}\), and the Earth third-body effect. In the current paper, we derive the dynamics of a lunar orbiter under the mentioned perturbations, which are shown to be dominant for the evolution of circumlunar orbits, and use von Zeipel’s method to obtain a transformation between mean and osculating elements. Whereas the dynamics of the mean elements do not include \(C_{22}\), and hence does not affect the equilibria leading to frozen orbits, \(C_{22}\) is present in the mean-to-osculating transformation, hence affecting the initialization of the physical circumlunar orbit. Simulations show that by using the newly-derived transformation, frozen orbits exhibit better behavior in terms of long-term stability about the mean values of eccentricity and argument of periapsis, especially for high orbits.  相似文献   

5.
Vinti’s potential is revisited for analytical propagation of the main satellite problem, this time in the context of relative motion. A particular version of Vinti’s spheroidal method is chosen that is valid for arbitrary elliptical orbits, encapsulating \(J_2\), \(J_3\), and generally a partial \(J_4\) in an orbit propagation theory without recourse to perturbation methods. As a child of Vinti’s solution, the proposed relative motion model inherits these properties. Furthermore, the problem is solved in oblate spheroidal elements, leading to large regions of validity for the linearization approximation. After offering several enhancements to Vinti’s solution, including boosts in accuracy and removal of some singularities, the proposed model is derived and subsequently reformulated so that Vinti’s solution is piecewise differentiable. While the model is valid for the critical inclination and nonsingular in the element space, singularities remain in the linear transformation from Earth-centered inertial coordinates to spheroidal elements when the eccentricity is zero or for nearly equatorial orbits. The new state transition matrix is evaluated against numerical solutions including the \(J_2\) through \(J_5\) terms for a wide range of chief orbits and separation distances. The solution is also compared with side-by-side simulations of the original Gim–Alfriend state transition matrix, which considers the \(J_2\) perturbation. Code for computing the resulting state transition matrix and associated reference frame and coordinate transformations is provided online as supplementary material.  相似文献   

6.
Near-Earth asteroids have attracted attention for both scientific and commercial mission applications. Due to the fact that the Earth–Moon \(\hbox {L}_{1}\) and \(\hbox {L}_{2}\) points are candidates for gateway stations for lunar exploration, and an ideal location for space science, capturing asteroids and inserting them into periodic orbits around these points is of significant interest for the future. In this paper, we define a new type of lunar asteroid capture, termed direct capture. In this capture strategy, the candidate asteroid leaves its heliocentric orbit after an initial impulse, with its dynamics modeled using the Sun–Earth–Moon restricted four-body problem until its insertion, with a second impulse, onto the \(\hbox {L}_{2}\) stable manifold in the Earth–Moon circular restricted three-body problem. A Lambert arc in the Sun-asteroid two-body problem is used as an initial guess and a differential corrector used to generate the transfer trajectory from the asteroid’s initial obit to the stable manifold associated with Earth–Moon \(\hbox {L}_{2}\) point. Results show that the direct asteroid capture strategy needs a shorter flight time compared to an indirect asteroid capture, which couples capture in the Sun–Earth circular restricted three-body problem and subsequent transfer to the Earth–Moon circular restricted three-body problem. Finally, the direct and indirect asteroid capture strategies are also applied to consider capture of asteroids at the triangular libration points in the Earth–Moon system.  相似文献   

7.
We review the origin and evolution of the atmospheres of Earth, Venus and Mars from the time when their accreting bodies were released from the protoplanetary disk a few million years after the origin of the Sun. If the accreting planetary cores reached masses \(\ge 0.5 M_\mathrm{Earth}\) before the gas in the disk disappeared, primordial atmospheres consisting mainly of H\(_2\) form around the young planetary body, contrary to late-stage planet formation, where terrestrial planets accrete material after the nebula phase of the disk. The differences between these two scenarios are explored by investigating non-radiogenic atmospheric noble gas isotope anomalies observed on the three terrestrial planets. The role of the young Sun’s more efficient EUV radiation and of the plasma environment into the escape of early atmospheres is also addressed. We discuss the catastrophic outgassing of volatiles and the formation and cooling of steam atmospheres after the solidification of magma oceans and we describe the geochemical evidence for additional delivery of volatile-rich chondritic materials during the main stages of terrestrial planet formation. The evolution scenario of early Earth is then compared with the atmospheric evolution of planets where no active plate tectonics emerged like on Venus and Mars. We look at the diversity between early Earth, Venus and Mars, which is found to be related to their differing geochemical, geodynamical and geophysical conditions, including plate tectonics, crust and mantle oxidation processes and their involvement in degassing processes of secondary \(\hbox {N}_2\) atmospheres. The buildup of atmospheric \(\hbox {N}_2\), \(\hbox {O}_2\), and the role of greenhouse gases such as \(\hbox {CO}_2\) and \(\hbox {CH}_4\) to counter the Faint Young Sun Paradox (FYSP), when the earliest life forms on Earth originated until the Great Oxidation Event \(\approx \) 2.3 Gyr ago, are addressed. This review concludes with a discussion on the implications of understanding Earth’s geophysical and related atmospheric evolution in relation to the discovery of potential habitable terrestrial exoplanets.  相似文献   

8.
Gravitational capture is a useful phenomenon in the design of the low energy transfer (LET) orbit for a space mission. In this paper, gravitational lunar capture based on the Sun–Earth–Moon bicircular model (BCM) in the restricted four body problem is studied. By the mechanical analysis in the space near the Moon, we first propose a new parameter \(k\) , the corrected ratio of the radial force, to investigate the influence of the radial force on the capture eccentricity in the BCM. Then, a parametric analysis is performed to detect the influences on the corrected ratio \(k\) . Considering the restriction of time-of-flight and corrected ratio, we investigate, respectively, the minimum capture eccentricity and the corrected minimum capture eccentricity. Via numerical analysis, we discover two special regions on the sphere of capture, in which the capture point possesses the global minimum capture eccentricity and corrected capture eccentricity. They denote the optimal capture regions in terms of minimizing the fuel consumption of the maneuver. According to the results obtained, some suggestions on the design of the LET orbit are given.  相似文献   

9.
Jupiter’s satellites are subject to strong tidal forces which result in variations of the gravitational potential and deformations of the satellites’ surfaces on the diurnal tidal cycle. Such variations are described by the Love numbers \(k_2\) and \(h_2\) for the tide-induced potential variation due to internal mass redistribution and the radial surface displacement, respectively. The phase-lags \( \phi _{k_2}\) and \( \phi _{h_2}\) of these complex numbers contain information about the rheological and dissipative states of the satellites. Starting from interior structure models and assuming a Maxwell rheology to compute the tidal deformation, we calculate the phase-lags in application to Ganymede and Europa. For both satellites we assume a decoupling of the outer ice-shell from the deep interior by a liquid subsurface water ocean. We show that, in this case, the phase-lag difference \(\varDelta \phi = \phi _{k_2}- \phi _{h_2}\) can provide information on the rheological and thermal state of the deep interiors if the viscosities of the deeper layers are small. In case of Ganymede, phase-lag differences can reach values of a few degrees for high-pressure ice viscosities \({<}10^{14}\) Pa s and would indicate a highly dissipative state of the deep interior. In this case \(\varDelta \phi \) is dominated by dissipation in the high-pressure ice layer rather than dissipation within the ice-I shell. These phase lags would be detectable from spacecraft in orbit around the satellite. For Europa \(\varDelta \phi \) could reach values exceeding \(20^\circ \) and phase-lag measurements could help distinguish between (1) a hot dissipative silicate mantle which would in thermal equilibrium correspond to a very thin outer ice-I shell and (2) a cold deep interior implying that dissipation would mainly occur in a thick (several tens of km) outer ice-I shell. These measurements are highly relevant for ESA’s Jupiter Icy Moons Explorer (JUICE) and NASA’s Europa Multiple Flyby Mission, both targeted for the Jupiter system.  相似文献   

10.
Kan Liou  Chin-Chun Wu 《Solar physics》2016,291(12):3777-3792
Interplanetary magnetic field and solar wind plasma density observed at 1 AU during Solar Cycle 23?–?24 (SC-23/24) minimum were significantly smaller than those during its previous solar cycle (SC-22/23) minimum. Because the Earth’s orbit is embedded in the slow wind during solar minimum, changes in the geometry and/or content of the slow wind region (SWR) can have a direct influence on the solar wind parameters near the Earth. In this study, we analyze solar wind plasma and magnetic field data of hourly values acquired by Ulysses. It is found that the solar wind, when averaging over the first (1995.6?–?1995.8) and third (2006.9?–?2008.2) Ulysses’ perihelion (\({\sim}\,1.4~\mbox{AU}\)) crossings, was about the same speed, but significantly less dense (\({\sim}\,34~\%\)) and cooler (\({\sim}\,20~\%\)), and the total magnetic field was \({\sim}\,30~\%\) weaker during the third compared to the first crossing. It is also found that the SWR was \({\sim}\,50~\%\) wider in the third (\({\sim}\,68.5^{\circ}\) in heliographic latitude) than in the first (\({\sim}\,44.8^{\circ}\)) solar orbit. The observed latitudinal increase in the SWR is sufficient to explain the excessive decline in the near-Earth solar wind density during the recent solar minimum without speculating that the total solar output may have been decreasing. The observed SWR inflation is also consistent with a cooler solar wind in the SC-23/24 than in the SC-22/23 minimum. Furthermore, the ratio of the high-to-low latitude photospheric magnetic field (or equatorward magnetic pressure force), as observed by the Mountain Wilson Observatory, is smaller during the third than the first Ulysses’ perihelion orbit. These findings suggest that the smaller equatorward magnetic pressure at the Sun may have led to the latitudinally-wider SRW observed by Ulysses in SC-23/24 minimum.  相似文献   

11.
We aim to probe the dynamic structure of the extended Solar neighborhood by calculating the radial metallicity gradients from orbit properties, which are obtained for axisymmetric and non-axisymmetric potential models, of red clump (RC) stars selected from the RAdial Velocity Experiment’s Fourth Data Release. Distances are obtained by assuming a single absolute magnitude value in near-infrared, i.e. \(M_{Ks}=-1.54\pm0.04\) mag, for each RC star. Stellar orbit parameters are calculated by using the potential functions: (i) for the MWPotential2014 potential, (ii) for the same potential with perturbation functions of the Galactic bar and transient spiral arms. The stellar age is calculated with a method based on Bayesian statistics. The radial metallicity gradients are evaluated based on the maximum vertical distance (\(z_{max}\)) from the Galactic plane and the planar eccentricity (\(e_{p}\)) of RC stars for both of the potential models. The largest radial metallicity gradient in the \(0< z_{max} \leq0.5\) kpc distance interval is \(-0.065\pm0.005~\mbox{dex}\,\mbox{kpc}^{-1}\) for a subsample with \(e_{p}\leq0.1\), while the lowest value is \(-0.014\pm0.006~\mbox{dex}\,\mbox{kpc}^{-1}\) for the subsample with \(e_{p}\leq0.5\). We find that at \(z_{max}>1\) kpc, the radial metallicity gradients have zero or positive values and they do not depend on \(e_{p}\) subsamples. There is a large radial metallicity gradient for thin disc, but no radial gradient found for thick disc. Moreover, the largest radial metallicity gradients are obtained where the outer Lindblad resonance region is effective. We claim that this apparent change in radial metallicity gradients in the thin disc is a result of orbital perturbation originating from the existing resonance regions.  相似文献   

12.
We will show that the period T of a closed orbit of the planar circular restricted three body problem (viewed on rotating coordinates) depends on the region it encloses. Roughly speaking, we show that, \(2 T=k\pi +\int _\Omega g\) where k is an integer, \(\Omega \) is the region enclosed by the periodic orbit and \(g:{\mathbb {R}}^2\rightarrow {\mathbb {R}}\) is a function that only depends on the constant C known as the Jacobian constant; it does not depend on \(\Omega \). This theorem has a Keplerian flavor in the sense that it relates the period with the space “swept” by the orbit. As an application we prove that there is a neighborhood around \(L_4\) such that every periodic solution contained in this neighborhood must move clockwise. The same result holds true for \(L_5\).  相似文献   

13.
We consider a Yukawa-type gravitational potential combined with the Poynting-Robertson effect. Dust particles originating within the asteroid belt and moving on circular and elliptic trajectories are studied and expressions for the time rate of change of their orbital radii and semimajor axes, respectively, are obtained. These expressions are written in terms of basic particle parameters, namely their density and diameter. Then, they are applied to produce expressions for the time required by the dust particles to reach the orbit of Earth. For the Yukawa gravitational potential, dust particles of diameter \(10^{ - 3}\) m in circular orbits require times of the order of \(8.557 \times 10^{6}\) yr and for elliptic orbits of eccentricities \(e =0.1, 0.5\) require times of \(9.396 \times 10^{6}\) and \(2.129 \times 10^{6}\) yr respectively to reach Earth’s orbit. Finally, various cases of the Yukawa potential are studied and the corresponding particle times to reach Earth’s are derived per case along with numerical results for circular and various elliptical orbits.  相似文献   

14.
In this note, by using Smale’s \(\alpha \)-theorem on the convergence of Newton’s method, the \(\alpha \)-sets of convergence of some starters of solving the elliptic Kepler’s equation are derived. For each starter we compute the exact \(\alpha \)-set in the eccentricity-main anomaly \((e,M)\in [0,1)\times [0,\pi ]\), showing that these sets are larger than those derived by Avendaño et al. (Celest Mech Dyn Astron 119:27–44, 2014). Further, new convergence tests based on the Newton–Kantorowitch theorem are given comparing with the derived from Smale’s \(\alpha \)-test.  相似文献   

15.
In this work we consider the Kepler problem with linear drag, and prove the existence of a continuous vector-valued first integral, obtained taking the limit as \(t\rightarrow +\infty \) of the Runge–Lenz vector. The norm of this first integral can be interpreted as an asymptotic eccentricity \(e_{\infty }\) with \(0\le e_{\infty } \le 1\). The orbits satisfying \(e_{\infty } <1\) approach the singularity by an elliptic spiral and the corresponding solutions \(x(t)=r(t)e^{i\theta (t)}\) have a norm r(t) that goes to zero like a negative exponential and an argument \(\theta (t)\) that goes to infinity like a positive exponential. In particular, the difference between consecutive times of passage through the pericenter, say \(T_{n+1} -T_n\), goes to zero as \(\frac{1}{n}\).  相似文献   

16.
We consider an elliptic restricted four-body system including three primaries and a massless particle. The orbits of the primaries are elliptic, and the massless particle moves under the mutual gravitational attraction. From the dynamic equations, a quasi-integral is obtained, which is similar to the Jacobi integral in the circular restricted three-body problem (CRTBP). The energy constant \(C\) determines the topology of zero velocity surfaces, which bifurcate at the equilibrium point. We define the concept of Hill stability in this problem, and a criterion for stability is deduced. If the actual energy constant \(C_{\mathrm{ac}}\ ( {>} 0 ) \) is bigger than or equal to the critical energy constant \(C_{\mathrm{cr}}\), the particle will be Hill stable. The critical energy constant is determined by the mass and orbits of the primaries. The criterion provides a way to capture an asteroid into the Earth–Moon system.  相似文献   

17.
Taking into consideration a probe moving in an elliptical orbit around a celestial body, the possibility of determining conditions which lead to constant values on average of all the orbit elements has been investigated here, considering the influence of the planetary oblateness and the long-term effects deriving from the attraction of several perturbing bodies. To this end, three equations describing the variation of orbit eccentricity, apsidal line and angular momentum unit vector have been first retrieved, starting from a vectorial expression of the Lagrange planetary equations and considering for the third-body perturbation the gravity-gradient approximation, and then exploited to demonstrate the feasibility of achieving the above-mentioned goal. The study has led to the determination of two families of solutions at constant mean orbit elements, both characterised by a co-planarity condition between the eccentricity vector, the angular momentum and a vector resulting from the combination of the orbital poles of the perturbing bodies. As a practical case, the problem of a probe orbiting the Moon has been faced, taking into account the temporal evolution of the perturbing poles of the Sun and Earth, and frozen solutions at argument of pericentre 0\(^{\circ }\) or 180\(^{\circ }\) have been found.  相似文献   

18.
China is planning to land a spacecraft on the farside of the Moon, a premiere, by 2018. In essence, the traditional tracking modes, based on direct visibility, cannot operate for the lunar farside lander tracking, and therefore a relay satellite, visible at the same time by both the lander and the Earth, will be required, operating in the so-called four-way mode (Earth-relay satellite-lander-relay satellite-Earth). In this paper, we firstly give the mathematical formulation of the four-way relay tracking mode and of its partial derivatives with respect to the relevant parameters, implemented in our POD software WUDOGS (Wuhan University Deep-space Orbit determination and Gravity recovery System). In a second step, in simulation mode, we apply this relay mode to determining lander coordinates, which are absolutely needed for a sample return mission, or to add constraints on rotation models of the Moon. The results show that with Doppler measurements at a 0.1 mm/s error level, the positioning of the farside lander could be done at centimeters level (1-\(\delta\)) in the case of a circumlunar relay satellite; and at a 5 meters level (1-\(\delta\)) in the case of a Lagrange point (L2) Halo relay satellite.  相似文献   

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

20.
This paper considers disruption of triple close approaches with low initial velocities and equal masses in the framework of statistical escape theory in a three-dimensional space. The statistical escape theory is based on the assumption that the phase trajectory of a triple system is quasi-ergodic. This system is described by allowing for both energy and angular momentum conservation in the phase space. In this paper, “possibility of escape” is derived with the formation of a binary on the basis of relative distances of the participating bodies. The complete statistical solutions (i.e. the semi-major axis \(a\), the distributions of eccentricity \(e\) of the binary, binary energy \({E}_{{b}}\), escape energy \({E}_{{s}}\) of escaper, and its escape velocity \({v}_{{s}}\)) of the system are derived from the allowable phase space volumes and are in good agreement with the numerical results in the range of perturbing velocities \({v}_{{i}}\)(\(10^{ - 1} \le {v}_{{i}} \le 10^{ - 10}\)) and directions of \({v}_{{i}}(0 \le \alpha _{{i}},\beta _{{i}},\gamma _{{i}} \le \pi )\), \({i} = 1,2,3\). In this paper, the double limit process has been applied to approximate the escape probability. Through this process, it is observed that the perturbing velocity \({v}_{{i}} \to 0^{ +} \), as the product of the semi-major axis \(a\) of the final binary and the square of the escape velocity \({v}_{{s}}\) approach 2/3, i.e. \({a} {v}_{{s}}^{2} \to 2 / 3\), whatever direction of \(\mathbf{v}_{{i}}\) may be.  相似文献   

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

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