Numerical Procedures for Obtaining P(e, i)

このQ&Aのポイント
  • This article presents numerical procedures for efficiently obtaining P(e, i) based on the framework of Hill's equations.
  • The procedures involve reducing the degrees of freedom of particle motion and simplifying the computation process to avoid difficulties.
  • Specific simplifications include neglecting certain orbits, using a two-body formula for collision determination, and omitting computations in certain regions.
回答を見る
  • ベストアンサー

この文章の和訳を教えてください。

4. Numerical procedures for obtaining <P(e, i)> Based on the results in Sect. 3, we develop numerical procedures for efficiently obtaining <P(e, i)>. By using the framework of Hill’s equations, one can reduce the degrees of freedom of particle motion, as described in Paper I. Nevertheless, they are still too numerous to compute orbits densely and uniformly over the five-dimensional phase space (e, i, b, τ_s, ω_s) of initial conditions. Furthermore, we frequently find orbits which comes very close to the protoplanet or which revolve complicatedly many times around the protoplanet. In these cases, an enormously long computation time is needed to calculate them. To avoid these difficulties, we introduce the following three simplifications: (i) We neglect n-recurrent orbits in evaluating <P(e, i)> when n≧3; that is, we stop an orbital calculation if the particle does not collide with the protoplanet even after two close encounters. (ii) When the planetesimal crosses the sphere surface of the two-body approximation, the two-body formula determines whether or not a collision occurs. (iii)We omit orbital computations in some regions of τ_s, which depend on a given set of e, i, and b; it is found empirically that there are no collision orbits in these regions. よろしくお願いします。

  • 英語
  • 回答数2
  • ありがとう数2

質問者が選んだベストアンサー

  • ベストアンサー
  • ddeana
  • ベストアンサー率74% (2976/4019)
回答No.1

4.<P(e, i)>を得る為の数値計算手法 3章における結果に基づき、効率的に<P(e, i)>得るための数値計算手法を開発する。論文Iで述べているように、ヒル方程式の枠組みを使うことにより粒子運動の自由度を縮小することが可能である。にもかかわらず、それらは初期条件の5次元位相空間(e, i, b, τ_s, ω_s)上で軌道を緻密で一様に計算するには、いまだ数が多すぎるのである。加えて、原始惑星の非常に近くまで接近するか、原始惑星の周りを何度も複雑に周回する軌道が頻繁にみつかる。こうしたケースではそれらを計算するのにコンピューターを使った尋常ではない長さの算出時間が必要となるのである。 これらの問題を避ける為に、次なる3つの単純化を紹介する。 (1)nが3と同等かそれよりも大きい場合、<P(e, i)>の評価において、n回帰性軌道を無視する。つまり、2つの接近遭遇後であっても粒子が原始惑星と衝突しなければ軌道計算をやめるということである。 (2)微惑星が、二体近似の球体表面を横切る時、衝突が起こるか否かにかかわらず2体公式を決定する。 (3)特定のe, i, そしてbに依存するτ_sのいくつかの領域での軌道計算を省略する。なぜならこうした領域に衝突軌道がないことが実験的にわかっているからだ。

mamomo3
質問者

お礼

どうもありがとうございました。 感謝いたします。

その他の回答 (1)

  • ddeana
  • ベストアンサー率74% (2976/4019)
回答No.2

ddeanaです。前回答の一部、simplificationの訳を「簡略化」に変更いたします。 理由:計算式における時間の短縮化などを念頭にしたものなので、「単純」よりも複雑な行程を減らす「簡略」の方が適していると考えた為。

mamomo3
質問者

お礼

了解いたしました。 補足いただきまして、ありがとうございます。

関連するQ&A

  • この文章の和訳をよろしくお願いします。

    3. Examples of orbital calculations To find efficient numerical procedures for obtaining <P(e ,i)>, we have made detailed orbital calculations for two typical cases, (e, i)=(0, 0) and (1, 0.5). The former is the simplest case with a single parameter b, and corresponds to that studied by Giuli (1968), Nishida (1983) and Petit and Hénon (1986). In the latter, the orbit changes with three parameters b, τ, and ω in a complicated manner. From this example, we can see the characteristic features of orbits in the three-dimensional case. In the present examples, it is supposed that a protoplanet with a mean mass density 3gcm^-3, orbits in the Earth’s region. Furthermore, we adopt 1% as the limiting accuracy criterion for use of the two-body approximation. These give 0.005 and 0.03 for the radius of the protoplanet and that of the two-body sphere, respectively (see Eqs. (12) and (13)). よろしくお願いします。

  • この文章の和訳をお願いします。

                    3. Results of the orbital calculations   In this paper, we describe our numerical results only for the special cases e_i~=0 and e_i~=4 in order to find the characteristic features of the two-body encounters of Keplerian particles. 3.1. For the case e_i~=0   In this case the orbital element, δ_i, loses its meaning because the particle orbits have no periheria and , hence, we can actually assign the initial condition by one parameter, b_i~. The orbital calculations are performed for about 3800 cases with various values of b_i~ from -10 to 10. よろしくお願いします。

  • この文章の和訳を教えてください。

    3.1. Case of e=0 and i=0 We have first calculated 6000 orbits in the parameter range of b from b_min=1.9 to b_max=2.5 at intervals of 0.0001. It is already known from previous studies (Nishida, and Petit and Hénon) that no collision orbits exist outside this region.    The orbits vary in a complicated way with the value of parameter b (see Petit and Hénon, 1986). In spite of the complex behavior of the orbits, we can classify them in terms of the number of encounters with the two-body sphere, from the standpoint of finding collision orbits. The classes are: (a) non-encounter orbit, (b) n-recurrent non-collision orbit, and (c) n-recurrent collision orbit, where the term “n-recurrent orbit” means the particle encounters n-times with the two-body sphere. That is, n-recurrent non-collision (or collision) orbits are those which fly off to infinity (or collide with the protoplanet ) after n-times encounters with the two-body sphere, while non-recurrent orbits are those which fly off without penetrating the two-body sphere. Examples of orbits in the classes (a), (b), and (c) are illustrated in Figs. 2,3, and 4, respectively. The above classification of orbits will be utilized for developing numerical procedures for obtaining <P(e, i)>, as described in the next section. よろしくお願いします。

  • この文章の和訳を教えてください。

    Collisional probability of planetesimals revolving in the solar gravitational field.III Summary. We have calculated the collisional rate of planetesimals upon the protoplanet, taking fully into account the effect of solar gravity. Our numerical scheme is based on Hill’s equations describing approximately the three –body problem. By the adoption of Hill’s equations , we can reduce the degrees of freedom of orbital motion. Furthermore, we made some simplifications: First, an orbital motion is determined by the formula of the two-body approximation when the distance between the protoplanet and a planetesimal is smaller than a certain critical length. Second, collision orbits in the chaotic zones are neglected in evaluating the collisional rate because of their very small measure. These simplifications enable us to save considerable computation time of orbital integration and, hence, to find numerically the phase volume occupied by collision orbits over wide ranges of orbital initial conditions. よろしくお願いします。

  • この文章の和訳を教えてください。

    Collisional probability of planetesimals revolving in the solar gravitational field.III Summary. We have calculated the collisional rate of planetesimals upon the protoplanet, taking fully into account the effect of solar gravity. Our numerical scheme is based on Hill’s equations describing approximately the three –body problem. By the adoption of Hill’s equations , we can reduce the degrees of freedom of orbital motion. Furthermore, we made some simplifications: First, an orbital motion is determined by the formula of the two-body approximation when the distance between the protoplanet and a planetesimal is smaller than a certain critical length. Second, collision orbits in the chaotic zones are neglected in evaluating the collisional rate because of their very small measure. These simplifications enable us to save considerable computation time of orbital integration and, hence, to find numerically the phase volume occupied by collision orbits over wide ranges of orbital initial conditions. よろしくお願いします。

  • 和訳してください。よろしくお願いします。

    In the present paper, we first describe initial conditions for orbital calculations in the framework of Hill’s equations (Sect.2). In Sect.3, orbital calculations are made for particular sets of e and i, i.e., (e,i)=(0,0) and (1,0.5). Based on these calculations and previous works: Papers I and II, Nishida (1983), Hénon and Petit (1986), and Petit and Hénon (1986), we develop an efficient numerical procedure obtaining <P(e, i)> (Sect. 4).According to this procedure, we systemically integrate the orbits of relative motion, to find collision orbits. Furthermore, compiling results of these orbital calculations, we find <P(e, i)> for various sets of e and i and compare them with those in the two-body approximation <P(e, i)>_2B (Sects. 5 to 7). The results will be compared with those of Nishida (1983) and Wetherill and Cox (1985) who also studied the collisional rate taking account of the effect of solar gravity (but their studies were restricted, as mentioned later).

  • この英文の和訳をお願いします。

      Before a detailed description of our numerical procedures, we comment on above simplifications. As seen in an example of e=1.0 and i=0.5, n-recurrent (n≧3) collision orbits contribute by only 1% or less to the collision rate. On the one hand, from calculations with other e and i, the degree of contribution by n-recurrent (n≧3) collision orbits is found to be largest in the case e≒1. Hence, the error in <P(e, i)> introduced by simplification (i) is of the order of 1% or less. As for the applicability of the two-body approximation, we have confirmed in PaperII that the orbit are well described by the two-body formula inside the two-body sphere, whose radius is given by Eq. (13). No appreciable error in <P(e, i)> comes from simplification (ii). Simplification (iii) follows the discussion in the last section. Using above simplifications, we have developed numerical procedures for obtaining <P(e, i)> efficiently; their flow chart is illustrated in Fig. 10. Choosing initial values of orbital elements (e, i, b, τ_s, ω_s), we start to compute numerically Hill’s equations (6) by an ordinary fourth-order Runge-Kutta method from a starting point given by Eqs. (21) and (22). The distance r is checked at every time step of the numerical integration. If the particle flies off to a sufficient distance from the protoplanet after approaching it, i.e., if |y|>y_0+2e,                     (26) then, the orbital computation is stopped. If a particle approaches the protoplanet and crosses the two-body sphere surface, i.e., if r≦r_cr, the two-body formula is employed to predict whether or not a collision occurs. When no collision occurs at the first encounter, the numerical integration of Hill’s equations is continued. Since a particle which enters the two-body sphere inevitably escapes from the sphere (see Paper II), the particle follows alternatives: one is that it departs to such a distance that Eq. (26) is satisfied, and the other is that it crosses the two-body sphere surface again. In the former case, we stop the computation, considering that the orbit is non-collisional. In the latter case, the occurrence of collision is checked in the same way as earlier by means of the two-body formula, and the orbital calculation is terminated. Using the numerical procedures developed in this way, we obtain <P(e, i)> in many sets of (e, i); the results are presented in the preceeding sections. Fig. 10. Flow chart of orbital calculation for finding collision orbits. Fig. 10. 拡大画像↓ http://www.fastpic.jp/images.php?file=2113240192.jpg 長文になりますが、よろしくお願いします。

  • この文章の和訳をお願いします。

                     2. Adopted assumptions and basic equations We consider two planetesimals revolving around the proto-Sun (being called the Sun). Here we assume that the mass of the one of these, which hereafter is called a protoplanet or simply a planet, is much larger than that of the other (being called a particle). We also assume that the planet moves circularly around the Sun without the influence of gravity of the particle. Furthermore, each orbit of the particle is limited in the ecliptic plane of the planetary orbit. Under these assumptions, the particle motion is simply given by a solution to the plane circular RTB problem. よろしくお願いします。

  • 和訳をお願いします。

    In order to obtain the total collisional rate <Γ(e_1, i_1)>, we have to find <P(e,i)> numerically by computing orbits of relative motion between the protoplanet and the planetesimal. Since <P(e,i)> should be provided for wide ranges of e and i with a sufficient accuracy, we are obliged to compute a very large number of orbits. In practice, it is an important problem to find an efficient method for numerical computation. In the second paper (Nakazawa et al., 1989b, referred to as Paper II), we have studied the validity of the two-body approximation and found that within the sphere of the two-body approximation (hereafter referred to as the two-body sphere), the relative motion can be well described by a solution to the two-body problem: the sphere radius has been found to be r_cr=0.03(a_0*/1AU)^(-1/4)(ε/0.01)^(1/2). ・・・・・(13) Within the sphere, the nearest distance can be predicted with an accuracy εby the well-known formula of the two-body encounter. We can expect the above result to be useful to reduce computation time for obtaining <P(e,i)> numerically. よろしくお願いします。

  • この文章の和訳をよろしくお願いします。

    According to Paper I, the total collisional rate <Γ(e_1,i_1)> of planetesimals upon the protoplanet with the heliocentric orbital elements e_1 and i_1 is given by <Γ(e_1,i_1)>=2π^2∫<n_2>e_i<P(e,i)>dedi, ・・・・・・・(9) where <n_2> is the distribution function of planetesimals averaged by the phase angles τ_1 and ω_1 of the protoplanet, and e and i are the orbital elements of relative motion between the protoplanet and the planetesimal at infinity. よろしくお願いします。