# High-accuracy statistical simulation of planetary accretion: I. Test of the accuracy by comparison with the solution to the stochastic coagulation equation

Satoshi Inaba, Hidekazu Tanaka, Keiji Ohtsuki, Kiyoshi Nakazawa

Research output: Contribution to journalArticle

16 Citations (Scopus)

### Abstract

The object of this series of studies is to develop a highly accurate statistical code for describing the planetary accumulation process. In the present paper, as a first step, we check the validity of the method proposed by Wetherill and Stewart (1989) by comparing the results obtained by their method with the analytical solution to the stochastic coagulation equation (or to a well-evaluated numerical solution). As the collisional probability A(ij) between bodies with masses of im1 and jm1 (m1 being the unit mass), we consider the two cases: one is A(ij) Proportional to i x j and another is A(ij) Proportional to min(i, j)(i(1/3) + j(1/3))(i + j). In both cases, it is known that runaway growth occurs. The latter case corresponds to a simplified model of the planetesimal accumulation. We assumed that a collision of two bodies leads to their coalescence. Wetherill and Stewart's method contains some parameters controlling the practical numerical computation. Among these, two parameters are important: the mass division parameter δ, which determines the mass ratio of the adjacent mass batches, and the time division parameter ε, which controls the size of a time step in numerical integration. Through a number of numerical simulations for the case of A(ij) = i x j, we find that when δ ≤ 1.6 and ε ≤ 0.03 the numerical simulation can reproduce the analytical solution within a certain level of accuracy independently of the size of the body system. For the case of the planetesimal accumulation, it is shown that the simulation with δ ≤ 1.3 and ε ≤ 0.04 can describe precisely runaway growth. Because the accumulation process is stochastic, in order to obtain reliable mean values it is necessary to take the ensemble mean of the numerical results obtained with different random number generators. It is also found that the number of simulations, N(c), demanded to obtain the reliable mean value is about 500 and does not strongly depend on the functional form of A(ij). From the viewpoint of the numerical handling, the above value of δ(≤ 1.3) and N(c)(~ 500) are reasonable and, hence, we conclude that the numerical method proposed by Wetherill and Stewart is a valid and useful method for describing the planetary accumulation process. The real planetary accumulation process is more complex since it is coupled with the velocity evolution of the planetesimals. In the subsequent paper, we will complete the high-accuracy statistical code which simulate the accumulation process coupled with the velocity evolution and test the accuracy of the code by comparing with the results of N-body simulation.

Original language English 205-217 13 Earth, Planets and Space 51 3 Published - 1999 Yes

### Fingerprint

coagulation
accretion
protoplanets
planetesimal
simulation
division
random numbers
stochastic processes
stochasticity
coalescence
comparison
test
numerical integration
mass ratios
coalescing
numerical method
generators
collision
collisions
method

### ASJC Scopus subject areas

• Earth and Planetary Sciences (miscellaneous)

### Cite this

High-accuracy statistical simulation of planetary accretion : I. Test of the accuracy by comparison with the solution to the stochastic coagulation equation. / Inaba, Satoshi; Tanaka, Hidekazu; Ohtsuki, Keiji; Nakazawa, Kiyoshi.

In: Earth, Planets and Space, Vol. 51, No. 3, 1999, p. 205-217.

Research output: Contribution to journalArticle

@article{2d32a87ac5754bb58f5c1bc481900623,
title = "High-accuracy statistical simulation of planetary accretion: I. Test of the accuracy by comparison with the solution to the stochastic coagulation equation",
abstract = "The object of this series of studies is to develop a highly accurate statistical code for describing the planetary accumulation process. In the present paper, as a first step, we check the validity of the method proposed by Wetherill and Stewart (1989) by comparing the results obtained by their method with the analytical solution to the stochastic coagulation equation (or to a well-evaluated numerical solution). As the collisional probability A(ij) between bodies with masses of im1 and jm1 (m1 being the unit mass), we consider the two cases: one is A(ij) Proportional to i x j and another is A(ij) Proportional to min(i, j)(i(1/3) + j(1/3))(i + j). In both cases, it is known that runaway growth occurs. The latter case corresponds to a simplified model of the planetesimal accumulation. We assumed that a collision of two bodies leads to their coalescence. Wetherill and Stewart's method contains some parameters controlling the practical numerical computation. Among these, two parameters are important: the mass division parameter δ, which determines the mass ratio of the adjacent mass batches, and the time division parameter ε, which controls the size of a time step in numerical integration. Through a number of numerical simulations for the case of A(ij) = i x j, we find that when δ ≤ 1.6 and ε ≤ 0.03 the numerical simulation can reproduce the analytical solution within a certain level of accuracy independently of the size of the body system. For the case of the planetesimal accumulation, it is shown that the simulation with δ ≤ 1.3 and ε ≤ 0.04 can describe precisely runaway growth. Because the accumulation process is stochastic, in order to obtain reliable mean values it is necessary to take the ensemble mean of the numerical results obtained with different random number generators. It is also found that the number of simulations, N(c), demanded to obtain the reliable mean value is about 500 and does not strongly depend on the functional form of A(ij). From the viewpoint of the numerical handling, the above value of δ(≤ 1.3) and N(c)(~ 500) are reasonable and, hence, we conclude that the numerical method proposed by Wetherill and Stewart is a valid and useful method for describing the planetary accumulation process. The real planetary accumulation process is more complex since it is coupled with the velocity evolution of the planetesimals. In the subsequent paper, we will complete the high-accuracy statistical code which simulate the accumulation process coupled with the velocity evolution and test the accuracy of the code by comparing with the results of N-body simulation.",
author = "Satoshi Inaba and Hidekazu Tanaka and Keiji Ohtsuki and Kiyoshi Nakazawa",
year = "1999",
language = "English",
volume = "51",
pages = "205--217",
journal = "Earth, Planets and Space",
issn = "1343-8832",
publisher = "Terra Scientific Publishing Company",
number = "3",

}

TY - JOUR

T1 - High-accuracy statistical simulation of planetary accretion

T2 - I. Test of the accuracy by comparison with the solution to the stochastic coagulation equation

AU - Inaba, Satoshi

AU - Tanaka, Hidekazu

AU - Ohtsuki, Keiji

AU - Nakazawa, Kiyoshi

PY - 1999

Y1 - 1999

N2 - The object of this series of studies is to develop a highly accurate statistical code for describing the planetary accumulation process. In the present paper, as a first step, we check the validity of the method proposed by Wetherill and Stewart (1989) by comparing the results obtained by their method with the analytical solution to the stochastic coagulation equation (or to a well-evaluated numerical solution). As the collisional probability A(ij) between bodies with masses of im1 and jm1 (m1 being the unit mass), we consider the two cases: one is A(ij) Proportional to i x j and another is A(ij) Proportional to min(i, j)(i(1/3) + j(1/3))(i + j). In both cases, it is known that runaway growth occurs. The latter case corresponds to a simplified model of the planetesimal accumulation. We assumed that a collision of two bodies leads to their coalescence. Wetherill and Stewart's method contains some parameters controlling the practical numerical computation. Among these, two parameters are important: the mass division parameter δ, which determines the mass ratio of the adjacent mass batches, and the time division parameter ε, which controls the size of a time step in numerical integration. Through a number of numerical simulations for the case of A(ij) = i x j, we find that when δ ≤ 1.6 and ε ≤ 0.03 the numerical simulation can reproduce the analytical solution within a certain level of accuracy independently of the size of the body system. For the case of the planetesimal accumulation, it is shown that the simulation with δ ≤ 1.3 and ε ≤ 0.04 can describe precisely runaway growth. Because the accumulation process is stochastic, in order to obtain reliable mean values it is necessary to take the ensemble mean of the numerical results obtained with different random number generators. It is also found that the number of simulations, N(c), demanded to obtain the reliable mean value is about 500 and does not strongly depend on the functional form of A(ij). From the viewpoint of the numerical handling, the above value of δ(≤ 1.3) and N(c)(~ 500) are reasonable and, hence, we conclude that the numerical method proposed by Wetherill and Stewart is a valid and useful method for describing the planetary accumulation process. The real planetary accumulation process is more complex since it is coupled with the velocity evolution of the planetesimals. In the subsequent paper, we will complete the high-accuracy statistical code which simulate the accumulation process coupled with the velocity evolution and test the accuracy of the code by comparing with the results of N-body simulation.

AB - The object of this series of studies is to develop a highly accurate statistical code for describing the planetary accumulation process. In the present paper, as a first step, we check the validity of the method proposed by Wetherill and Stewart (1989) by comparing the results obtained by their method with the analytical solution to the stochastic coagulation equation (or to a well-evaluated numerical solution). As the collisional probability A(ij) between bodies with masses of im1 and jm1 (m1 being the unit mass), we consider the two cases: one is A(ij) Proportional to i x j and another is A(ij) Proportional to min(i, j)(i(1/3) + j(1/3))(i + j). In both cases, it is known that runaway growth occurs. The latter case corresponds to a simplified model of the planetesimal accumulation. We assumed that a collision of two bodies leads to their coalescence. Wetherill and Stewart's method contains some parameters controlling the practical numerical computation. Among these, two parameters are important: the mass division parameter δ, which determines the mass ratio of the adjacent mass batches, and the time division parameter ε, which controls the size of a time step in numerical integration. Through a number of numerical simulations for the case of A(ij) = i x j, we find that when δ ≤ 1.6 and ε ≤ 0.03 the numerical simulation can reproduce the analytical solution within a certain level of accuracy independently of the size of the body system. For the case of the planetesimal accumulation, it is shown that the simulation with δ ≤ 1.3 and ε ≤ 0.04 can describe precisely runaway growth. Because the accumulation process is stochastic, in order to obtain reliable mean values it is necessary to take the ensemble mean of the numerical results obtained with different random number generators. It is also found that the number of simulations, N(c), demanded to obtain the reliable mean value is about 500 and does not strongly depend on the functional form of A(ij). From the viewpoint of the numerical handling, the above value of δ(≤ 1.3) and N(c)(~ 500) are reasonable and, hence, we conclude that the numerical method proposed by Wetherill and Stewart is a valid and useful method for describing the planetary accumulation process. The real planetary accumulation process is more complex since it is coupled with the velocity evolution of the planetesimals. In the subsequent paper, we will complete the high-accuracy statistical code which simulate the accumulation process coupled with the velocity evolution and test the accuracy of the code by comparing with the results of N-body simulation.

UR - http://www.scopus.com/inward/record.url?scp=0032773141&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0032773141&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:0032773141

VL - 51

SP - 205

EP - 217

JO - Earth, Planets and Space

JF - Earth, Planets and Space

SN - 1343-8832

IS - 3

ER -