共查询到20条相似文献，搜索用时 0 毫秒
1.
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. 相似文献
2.
Coronal mass ejections (CMEs) are large-scale eruptions of plasma from the Sun, which play an important role in space weather. Faraday rotation is the rotation of the plane of polarization that results when a linearly polarized signal passes through a magnetized plasma such as a CME. Faraday rotation is proportional to the path integral through the plasma of the electron density and the line-of-sight component of the magnetic field. Faraday-rotation observations of a source near the Sun can provide information on the plasma structure of a CME shortly after launch. We report on simultaneous white-light and radio observations made of three CMEs in August 2012. We made sensitive Very Large Array (VLA) full-polarization observations using 1?–?2 GHz frequencies of a constellation of radio sources through the solar corona at heliocentric distances that ranged from 6?–?\(15~\mathrm{R}_{\odot}\). Two sources (0842+1835 and 0900+1832) were occulted by a single CME, and one source (0843+1547) was occulted by two CMEs. In addition to our radioastronomical observations, which represent one of the first active hunts for CME Faraday rotation since Bird et al. (Solar Phys., 98, 341, 1985) and the first active hunt using the VLA, we obtained white-light coronagraph images from the Large Angle and Spectrometric Coronagraph (LASCO) C3 instrument to determine the Thomson-scattering brightness [\(\mathrm{B}_{\mathrm{T}}\)], providing a means to independently estimate the plasma density and determine its contribution to the observed Faraday rotation. A constant-density force-free flux rope embedded in the background corona was used to model the effects of the CMEs on \(\mathrm{B}_{\mathrm{T}}\) and Faraday rotation. The plasma densities (\(6\,\mbox{--}\,22\times10^{3}~\mbox{cm}^{-3}\)) and axial magnetic-field strengths (2?–?12 mG) inferred from our models are consistent with the modeling work of Liu et al. (Astrophys. J., 665, 1439, 2007) and Jensen and Russell (Geophys. Res. Lett., 35, L02103, 2008), as well as previous CME Faraday-rotation observations by Bird et al. (1985). 相似文献
3.
Fast Radio Bursts (FRBs) are short duration highly energetic dispersed radio pulses. We developed a generic formalism (Bera et al. 2016, MNRAS, 457, 2530) to estimate the FRB detection rate for any radio telescope with given parameters. By using this model, we estimated the FRB detection rate for two Indian radio telescope; the Ooty Wide Field Array (OWFA) (Bhattacharyya et al. 2017, J. Astrophys. Astr., 38, 17) and the upgraded Giant Metrewave Radio Telescope (uGMRT) (Bhattacharyya et al. 2018, J. Astrophys. Astr.) with three beam-forming modes. Here, we summarize these two works. We considered the energy spectrum of FRBs as a power law and the energy distribution of FRBs as a Dirac delta function and a Schechter luminosity function. We also considered two scattering models proposed by Bhat et al. (2004, Astrophys. J. Suppl. Series, 206, 1) and Macquart & Koay (2013, ApJ, 776, 125) for these works and we consider FRB pulse without scattering as a special case. We found that the future prospects of detecting FRBs by using these two Indian radio telescopes is good. They are capable to detect a significant number of FRBs per day. According to our prediction, we can detect \(\sim 10^5{-}10^8\), \(\sim 10^3{-}10^6\) and \(\sim 10^5{-}10^7\) FRBs per day by using OWFA, commensal systems of GMRT and uGMRT respectively. Even a non detection of the predicted events will be very useful in constraining FRB properties. 相似文献
4.
In this paper, we have studied the magnetized quark matter (QM) and strange quark matter (SQM) distributions in the presence of \(f(R,T)\) gravity in the background of Friedmann-Lemaître-Robertson-Walker (FLRW) metric. To get exact solutions of modified field equations we have used \(f(R,T ) = R + 2 f(T)\) model given by Harko et al. with two different parametrization of geometrical parameters i.e. the parametrization of the deceleration parameter \(q \), and the scale factor \(a \) in hybrid expansion form. Also, we have obtained Einstein Static Universe (ESU) solutions for QM and SQM distributions in \(f(R,T)\) gravity and General Relativity (GR). All models in \(f(R,T)\) gravity and GR for FRW and ESU Universes with QM also SQM distributions, we get zero magnetic field. These results agree with the solutions of Akta? and Aygün in \(f(R,T)\) gravity. However, we have also discussed the physical consequences of our obtained models. 相似文献
5.
IUE has made very successful long term and intense short time-scale monitoring spectroscopic study of NGC 4151, a Seyfert 1 galaxy for over nearly 18 years from its launch in 1978 to 1996. The long-term observations have been useful in understanding the complex relation between UV continuum and emission line variability Seyfert galaxies. In this paper, we present the results of our studies on the short-timescale intense monitoring campaign of NGC 4151 undertaken during December 1–15, 1993. A most intense monitoring observation of NGC 4151 was carried out by IUE in 1993, when the source was at its historical high flux state with a shortest interval of 70 min between two successive observations. We present our results on emission line and continuum variability amplitudes characterized by \(F_{\mathrm{var}}\) method. We found highest variability of nearly 8.3% at 1325 \(\AA \) continuum with a smallest amplitude of 4% at 2725 \(\AA \). The relative variability amplitudes (\(R_\mathrm{max}\)) have been found to be 1.372, 1.319, 1.302 and 1.182 at 1325, 1475, 1655 and 2725 \(\AA \) continuum respectively. The continuum and emission line variability characteristics obtained in the present analysis are in very good agreement with the results obtained by Edelson et al. (1996) and Crenshaw et al. (1996) from the analysis of the same observational spectral data. The large amplitude rapid variability characteristics obtained in our study have been attributed to the continuum reprocessing of X-rays absorbed by the material in the accretion disk as proposed by Shakura and Sunyaev (1973). The continuum and emission light curves have shown four distinct high amplitude events of flux maxima during the intense monitoring campaign of 15 days, providing a good limit on the amplitude of UV variability and the BLR size in low luminosity Seyfert galaxies and are useful for constraining the continuum emission models. The decreasing \(F_{\mathrm{var}}\) amplitude of UV continuum with respect to increasing wavelength obtained in the present study and consistent with similar observations by Edelson et al. (1996) and Crenshaw et al. (1996) is a significant result of the intense monitoring observations. 相似文献
6.
Here we report an in-depth reanalysis of an article by Vats et al. (Astrophys. J. 548, L87, 2001) that was based on measurements of differential rotation with altitude as a function of observing frequencies (as lower and higher frequencies indicate higher and lower heights, respectively) in the solar corona. The radial differential rotation of the solar corona is estimated from daily measurements of the disc-integrated solar radio flux at 11 frequencies: 275, 405, 670, 810, 925, 1080, 1215, 1350, 1620, 1755, and 2800 MHz. We use the same data as were used in Vats et al. (2001), but instead of the twelfth maxima of autocorrelograms used there, we use the first secondary maximum to derive the synodic rotation period. We estimate synodic rotation by Gaussian fit of the first secondary maximum. Vats et al. (2001) reported that the sidereal rotation period increases with increasing frequency. The variation found by them was from 23.6 to 24.15 days in this frequency range, with a difference of only 0.55 days. The present study finds that the sidereal rotation period increases with decreasing frequency. The variation range is from 24.4 to 22.5 days, and the difference is about three times larger (1.9 days). However, both studies give a similar rotation period at 925 MHz. In Vats et al. (2001) the Pearson’s factor with trend line was 0.86, whereas present analysis obtained a \({\sim}\,0.97\) Pearson’s factor with the trend line. Our study shows that the solar corona rotates more slowly at higher altitudes, which contradicts the findings reported in Vats et al. (2001). 相似文献
7.
This work is a continuation of our previous articles (Yermolaev et al. in J. Geophys. Res.120, 7094, 2015 and Yermolaev et al. in Solar Phys.292, 193, 2017), which describe the average temporal profiles of interplanetary plasma and field parameters in large-scale solar-wind (SW) streams: corotating interaction regions (CIRs), interplanetary coronal mass ejections (ICMEs, including both magnetic clouds (MCs) and ejecta), and sheaths as well as interplanetary shocks (ISs). Changes in the longitude angle, \(\varphi\), in CIRs from ?2 to \(2^{\circ}\) agree with earlier results (e.g. Gosling and Pizzo, 1999). We have also analyzed the average temporal profiles of the bulk velocity angles in sheaths and ICMEs. We have found that the angle \(\varphi\) in ICMEs changes from 2 to \(-2^{\circ}\), while in sheaths it changes from ?2 to \(2^{\circ}\) (similar to the change in CIRs), i.e. the angle in CIRs and sheaths deflects in the opposite sense to ICMEs. When averaging the latitude angle \(\vartheta\) on all the intervals of the chosen SW types, the angle \(\vartheta\) is almost constant at \({\sim}\,1^{\circ}\). We made for the first time a selection of SW events with increasing and decreasing \(\vartheta\) and found that the average \(\vartheta\) temporal profiles in the selected events have the same “integral-like” shape as for \(\varphi\). The difference in \(\varphi\) and \(\vartheta\) average profiles is explained by the fact that most events have increasing profiles for the angle in the ecliptic plane as a result of solar rotation, while for the angle in the meridional plane, the numbers of events with increasing and decreasing profiles are equal. 相似文献
8.
We investigate a method to test whether a numerically computed model coronal magnetic field \({\boldsymbol {B}}\) departs from the divergence-free condition (also known as the solenoidality condition). The test requires a potential field \({\boldsymbol {B}}_{0}\) to be calculated, subject to Neumann boundary conditions, given by the normal components of the model field \({\boldsymbol {B}}\) at the boundaries. The free energy of the model field may be calculated using \(\frac{1}{2\mu _{0}}\int ({\boldsymbol {B}}-{\boldsymbol {B}}_{0})^{2}\mathrm{d}V\), where the integral is over the computational volume of the model field. A second estimate of the free energy is provided by calculating \(\frac{1}{2\mu _{0}}\int {\boldsymbol {B}}^{2}\,\mathrm{d}V-\frac{1}{2\mu _{0}}\int {\boldsymbol {B}}_{0}^{2}\,\mathrm{d}V\). If \({\boldsymbol {B}}\) is divergence free, the two estimates of the free energy should be the same. A difference between the two estimates indicates a departure from \(\nabla \cdot {\boldsymbol {B}}=0\) in the volume. The test is an implementation of a procedure proposed by Moraitis et al. (Solar Phys.289, 4453, 2014) and is a simpler version of the Helmholtz decomposition procedure presented by Valori et al. (Astron. Astrophys.553, A38, 2013). We demonstrate the test in application to previously published nonlinear force-free model fields, and also investigate the influence on the results of the test of a departure from flux balance over the boundaries of the model field. Our results underline the fact that, to make meaningful statements about magnetic free energy in the corona, it is necessary to have model magnetic fields that satisfy the divergence-free condition to a good approximation. 相似文献
9.
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. 相似文献
10.
Abhas Mitra 《Astrophysics and Space Science》2011,332(1):43-48
By applying Birkhoff’s theorem to the problem of the general relativistic collapse of a uniform density dust, we directly show that the density of the dust ρ=0 even when its proper number density n would be assumed to be finite! The physical reason behind this exact result can be traced back to the observation of Arnowitt et al. (Phys. Rev. Lett. 4: 375, 1960) that the gravitational mass of a neutral point particle is zero: m=0. And since, a dust is a mere collection of neutral point particles, unlike a continuous hydrodynamic fluid, its density ρ=mn=0. It is nonetheless found that for k=?1, a homogeneous dust can collapse and expand special relativistically in the fashion of a Milne universe. Thus, in reality, general relativistic homogeneous dust collapse does not lead to the formation of any black hole in conformity of many previous studies (Logunov et al., Phys. Part. Nucl. 37: 317, 2006; Kiselev et al., Theor. Math. Phys. 164: 972, 2010; Mitra, J. Math. Phys. 50: 042502, 2009a; Suggett, J. Phys. A 12: 375 1979b). Interestingly, this result is in agreement with the intuition of Oppenheimer and Snyder (Phys. Rev. 56: 456, 1939) too:“Physically such a singularity would mean that the expressions used for the energy-momentum tensor does not take into account some essential physical fact which would really smooth the singularity out. Further, a star in its early stages of development would not possess a singular density or pressure, it is impossible for a singularity to develop in a finite time.” 相似文献
11.
We studied plane symmetric cosmological model in the presence of quark and strange quark matter with the help of f(R, T) theory. To decipher solutions of plane symmetric space-time, we used power law relation between scale factor and deceleration parameter. We considered the special law of variation of Hubble’s parameter proposed by Berman (Nuovo Cimento B74, 182, 1983) which yields constant deceleration parameter. We also discussed the physical behavior of the solutions by using some physical parameters. 相似文献
12.
D. O. Gough 《Solar physics》2017,292(5):70
It has been argued (Gough and McIntyre in Nature 394, 755, 1998) that the only way for the radiative interior of the Sun to be rotating uniformly in the face of the differentially rotating convection zone is for it to be pervaded by a large-scale magnetic field, a field which is responsible also for the thinness of the tachocline. It is most likely that this field is the predominantly dipolar residual component of a tangled primordial field that was present in the interstellar medium from which the Sun condensed (Braithwaite and Spruit in Nature 431, 819, 2004), and that advection by the meridional flow in the tachocline has caused the dipole axis to be inclined from the axis of rotation by about \(60^{\circ}\) (Gough in Geophys. Astrophys. Fluid Dyn., 106, 429, 2012). It is suggested here that, notwithstanding its turbulent passage through the convection zone, a vestige of that field is transmitted by the solar wind to Earth, where it modulates the geomagnetic field in a periodic way. The field variation reflects the inner rotation of the Sun, and, unlike turbulent-dynamo-generated fields, must maintain phase. I report here a new look at an earlier analysis of the geomagnetic field by Svalgaard and Wilcox (Solar Phys. 41, 461, 1975), which reveals evidence for appropriate phase coherence, thereby adding support to the tachocline theory. 相似文献
13.
The term “jumping” Trojan was introduced by Tsiganis et al. (Astron Astrophys 354:1091–1100, 2000) in their studies of long-term dynamics exhibited by the asteroid (1868) Thersites, which had been observed to jump from librations around \(L_4\) to librations around \(L_5\). Another example of a “jumping” Trojan was found by Connors et al. (Nature 475:481–483, 2011): librations of the asteroid 2010 TK7 around the Earth’s libration point \(L_4\) preceded by its librations around \(L_5\). We explore the dynamics of “jumping” Trojans under the scope of the restricted planar elliptical three-body problem. Via double numerical averaging we construct evolutionary equations, which allow analyzing transitions between different regimes of orbital motion. 相似文献
14.
We examine the average magnetic field magnitude (\(| \boldsymbol{B} | \equiv B\)) within magnetic clouds (MCs) observed by the Wind spacecraft from 1995 to July 2015 to understand the difference between this \(B\) and the ideal \(B\)-profiles expected from using the static, constant-\(\alpha\), force-free, cylindrically symmetric model for MCs of Lepping, Jones, and Burlaga (J. Geophys. Res. 95, 11957, 1990, denoted here as the LJB model). We classify all MCs according to an assigned quality, \(Q_{0}\) (\(= 1, 2, 3\), for excellent, good, and poor). There are a total of 209 MCs and 124 when only \(Q_{0} = 1\), 2 cases are considered. The average normalized field with respect to the closest approach (\(\mathit{CA}\)) is stressed, where we separate cases into four \(\mathit{CA}\) sets centered at 12.5 %, 37.5 %, 62.5 %, and 87.5 % of the average radius; the averaging is done on a percentage-duration basis to treat all cases the same. Normalized \(B\) means that before averaging, the \(B\) for each MC at each point is divided by the LJB model-estimated \(B\) for the MC axis, \(B_{0}\). The actual averages for the 209 and 124 MC sets are compared to the LJB model, after an adjustment for MC expansion (e.g. Lepping et al. in Ann. Geophys. 26, 1919, 2008). This provides four separate difference-relationships, each fitted with a quadratic (Quad) curve of very small \(\sigma\). Interpreting these Quad formulae should provide a comprehensive view of the variation in normalized \(B\) throughout the average MC, where we expect external front and rear compression to be part of its explanation. These formulae are also being considered for modifying the LJB model. This modification will be used in a scheme for forecasting the timing and magnitude of magnetic storms caused by MCs. Extensive testing of the Quad formulae shows that the formulae are quite useful in correcting individual MC \(B\)-profiles, especially for the first \({\approx\,}1/3\) of these MCs. However, the use of this type of \(B\) correction constitutes a (slight) violation of the force-free assumption used in the original LJB MC model. 相似文献
15.
The gravitational interaction between two objects on similar orbits can effect noticeable changes in the orbital evolution even if the ratio of their masses to that of the central body is vanishingly small. Christou (Icarus 174:215–229, 2005) observed an occasional resonant lock in the differential node \(\varDelta \varOmega \) between two members in the Himalia irregular satellite group of Jupiter in the N-body simulations (corresponding mass ratio \(\sim 10^{-9}\)). Using a semianalytical approach, we have reproduced this phenomenon. We also demonstrate the existence of two additional types of resonance, involving angle differences \(\varDelta \omega \) and \(\varDelta (\varOmega +\varpi )\) between two group members. These resonances cause secular oscillations in eccentricity and/or inclination on timescales \(\sim \)1 Myr. We locate these resonances in (a, e, i) space and analyse their topological structure. In subsequent N-body simulations, we confirm these three resonances and find a fourth one involving \(\varDelta \varpi \). In addition, we study the occurrence rates and the stability of the four resonances from a statistical perspective by integrating 1000 test particles for 100 Myr. We find \(\sim \)10 to 30 librators for each of the resonances. Particularly, the nodal resonance found by Christou is the most stable: 2 particles are observed to stay in libration for the entire integration. 相似文献
16.
Amit Mittal Rajiv Agarwal Md Sanam Suraj Monika Arora 《Astrophysics and Space Science》2018,363(5):109
This paper deals with the photo-gravitational restricted four-body problem (PR4BP) with variable mass. Following the procedure given by Gascheau (C. R. 16:393–394, 1843) and Routh (Proc. Lond. Math. Soc. 6:86–97, 1875), the conditions of linear stability of Lagrange triangle solution in the PR4BP are determined. The three radiating primaries having masses \(m_{1}\), \(m_{2}\) and \(m_{3}\) in an equilateral triangle with \(m_{2}=m_{3}\) will be stable as long as they satisfy the linear stability condition of the Lagrangian triangle solution. We have derived the equations of motion of the mentioned problem and observed that there exist eight libration points for a fixed value of parameters \(\gamma (\frac{m \ \text{at time} \ t}{m \ \text{at initial time}}, 0<\gamma\leq1 )\), \(\alpha\) (the proportionality constant in Jeans’ law (Astronomy and Cosmogony, Cambridge University Press, Cambridge, 1928), \(0\leq\alpha\leq2.2\)), the mass parameter \(\mu=0.005\) and radiation parameters \(q_{i}, (0< q_{i}\leq1, i=1, 2, 3)\). All the libration points are non-collinear if \(q_{2}\neq q_{3}\). It has been observed that the collinear and out-of-plane libration points also exist for \(q_{2}=q_{3}\). In all the cases, each libration point is found to be unstable. Further, zero velocity curves (ZVCs) and Newton–Raphson basins of attraction are also discussed. 相似文献
17.
We perform a statistical analysis on 157 M-class soft X-ray flares observed during 1997?–?2014 with and without deca-hectometric (DH) type II radio bursts aiming at the reasons for the non-occurrence of DH type II bursts in certain events. All the selected events are associated with halo Coronal Mass Ejections (CMEs) detected by the Solar and Heliospheric Observatory (SOHO) / Large Angle Spectrometric and COronograph (LASCO). Out of 157 events, 96 (61%; “Group I”) events are associated with a DH type II burst observed by the Radio and Plasma Wave (WAVES) experiment onboard the Wind spacecraft and 61 (39%; “Group II”) events occur without a DH type II burst. The mean CME speed of Group I is \(1022~\mbox{km}/\mbox{s}\) and that of Group II is \(647~\mbox{km}/\mbox{s}\). It is also found that the properties of the selected M-class flares such as flare intensity, rise time, duration and decay time are greater for the DH associated flares than the non-DH flares. Group I has a slightly larger number (56%) of western events than eastern events (44%), whereas Group II has a larger number of eastern events (62%) than western events (38%). We also compare this analysis with the previous study by Lawrance, Shanmugaraju, and Vr?nak (Solar Phys. 290, 3365L, 2015) concerning X-class flares and confirm that high-intensity flares (X-class and M-class) have the same trend in the CME and flare properties. Additionally we consider aspects like acceleration and the possibility of CME-streamer interaction. The average deceleration of CMEs with DH type II bursts is weaker (\(a = - 4.39\mbox{ m}/\mbox{s}^{2}\)) than that of CMEs without a type II burst (\(a = -12.21\mbox{ m}/\mbox{s}^{2}\)). We analyze the CME-streamer interactions for Group I events using the model proposed by Mancuso and Raymond (Astron. Astrophys. 413, 363, 2004) and find that the interaction regions are the most probable source regions for DH type II radio bursts. 相似文献
18.
H. A. Folonier S. Ferraz-Mello E. Andrade-Ines 《Celestial Mechanics and Dynamical Astronomy》2018,130(12):78
This paper deals with a new formulation of the creep tide theory (Ferraz-Mello in Celest Mech Dyn Astron 116:109, 2013—Paper I) and with the tidal dissipation predicted by the theory in the case of stiff bodies whose rotation is not synchronous but is oscillating around the synchronous state with a period equal to the orbital period. We show that the tidally forced libration influences the amount of energy dissipated in the body and the average perturbation of the orbital elements. This influence depends on the libration amplitude and is generally neglected in the study of planetary satellites. However, they may be responsible for a 27% increase in the dissipation of Enceladus. The relaxation factor necessary to explain the observed dissipation of Enceladus (\(\gamma =1.2{-}3.8\times 10^{-7}\ \mathrm{s}^{-1}\)) has the expected order of magnitude for planetary satellites and corresponds to the viscosity \(0.6{-}1.9 \times 10^{14}\) Pa s, which is in reasonable agreement with the value recently estimated by Efroimsky (Icarus 300:223, 2018) (\(0.24 \times 10^{14}\) Pa s) and with the value adopted by Roberts and Nimmo (Icarus 194:675, 2008) for the viscosity of the ice shell (\(10^{13}{-}10^{14}\) Pa s). For comparison purposes, the results are extended also to the case of Mimas and are consistent with the negligible dissipation and the absence of observed tectonic activity. The corrections of some mistakes and typos of paper II (Ferraz-Mello in Celest Mech Dyn Astron 122:359, 2015) are included at the end of the paper. 相似文献
19.
Kiyoshi Kawabata 《Astrophysics and Space Science》2018,363(1):1
We have established an iterative scheme to calculate with 15-digit accuracy the numerical values of Ambartsumian-Chandrasekhar’s \(H\)-functions for anisotropic scattering characterized by the four-term phase function: the method incorporates some advantageous features of the iterative procedure of Kawabata (Astrophys. Space Sci. 358:32, 2015) and the double-exponential integration formula (DE-formula) of Takahashi and Mori (Publ. Res. Inst. Math. Sci. Kyoto Univ. 9:721, 1974), which proved highly effective in Kawabata (Astrophys. Space Sci. 361:373, 2016). Actual calculations of the \(H\)-functions have been carried out employing 27 selected cases of the phase function, 56 values of the single scattering albedo \(\varpi_{0}\), and 36 values of an angular variable \(\mu(= \cos\theta)\), with \(\theta\) being the zenith angle specifying the direction of incidence and/or emergence of radiation. Partial results obtained for conservative isotropic scattering, Rayleigh scattering, and anisotropic scattering due to a full four-term phase function are presented. They indicate that it is important to simultaneously verify accuracy of the numerical values of the \(H\)-functions for \(\mu<0.05\), the domain often neglected in tabulation. As a sample application of the isotropic scattering \(H\)-function, an attempt is made in Appendix to simulate by iteratively solving the Ambartsumian equation the values of the plane and spherical albedos of a semi-infinite, homogeneous atmosphere calculated by Rogovtsov and Borovik (J. Quant. Spectrosc. Radiat. Transf. 183:128, 2016), who employed their analytical representations for these quantities and the single-term and two-term Henyey-Greenstein phase functions of appreciably high degrees of anisotropy. While our results are in satisfactory agreement with theirs, our procedure is in need of a faster algorithm to routinely deal with problems involving highly anisotropic phase functions giving rise to near-conservative scattering. 相似文献
20.
Rogerio Deienno Diogo Merguizo Sanchez Antonio Fernando Bertachini de Almeida Prado Georgi Smirnov 《Celestial Mechanics and Dynamical Astronomy》2016,124(4):433-456
We analyze the families of central configurations of the spatial 5-body problem with four masses equal to 1 when the fifth mass m varies from 0 to \(+\infty \). In particular we continue numerically, taking m as a parameter, the central configurations (which all are symmetric) of the restricted spatial (\(4+1\))-body problem with four equal masses and \(m=0\) to the spatial 5-body problem with equal masses (i.e. \(m=1\)), and viceversa we continue the symmetric central configurations of the spatial 5-body problem with five equal masses to the restricted (\(4+1\))-body problem with four equal masses. Additionally we continue numerically the symmetric central configurations of the spatial 5-body problem with four equal masses starting with \(m=1\) and ending in \(m=+\infty \), improving the results of Alvarez-Ramírez et al. (Discrete Contin Dyn Syst Ser S 1: 505–518, 2008). We find four bifurcation values of m where the number of central configuration changes. We note that the central configurations of all continued families varying m from 0 to \(+\infty \) are symmetric. 相似文献