TY - JOUR
T1 - Collisional disruption of planetesimals in the gravity regime with iSALE code
T2 - Comparison with SPH code for purely hydrodynamic bodies
AU - Suetsugu, Ryo
AU - Tanaka, Hidekazu
AU - Kobayashi, Hiroshi
AU - Genda, Hidenori
N1 - Funding Information:
We thank Thomas Davison and an anonymous reviewer for their useful comments on our manuscript. We gratefully acknowledge the developers of iSALE, including Gareth Collins, Kai Wünnemann, Boris Ivanov, Jay Melosh, Dirk Elbeshausen, and Thomas Davison. This work was supported by JSPS KAKENHI Grant Numbers JP26287101 and JP15H03716 . H.K. was supported by Grant-in-Aid for Scientific Research ( JP17K05632 , JP17H01105 , JP17H01103 ) and JSPS Core-to-Core Program “International Network of Planetary Sciences”. H.G. was also supported by JSPS KAKENHI Grant Number JP17H02990 .
PY - 2018/11/1
Y1 - 2018/11/1
N2 - In most of the previous studies related to collisional disruption of planetesimals in the gravity regime, Smoothed Particle Hydrodynamics (SPH) simulations have been used. On the other hand, impact simulations using grid-based hydrodynamic code have not been sufficiently performed. In the present study, we execute impact simulations in the gravity regime using the shock-physics code iSALE, which is a grid-based Eulerian hydrocode. We examine the dependence of the critical specific impact energy QRD * on impact conditions for a wide range of specific impact energy (QR) from disruptive collisions to erosive collisions, and compare our results with previous studies. We find that collision outcomes of the iSALE simulation agree well with those of the SPH simulation. Detailed analysis mainly gives three results. (1) The value of QRD * depends on numerical resolution, and is close to convergence with increasing numerical resolution. The difference in converged value of QRD * between the iSALE code and the SPH code is within 30%. (2) Ejected mass normalized by total mass (Mej/Mtot) generally depends on various impact conditions. However, when QR is normalized by QRD * that is calculated for each impact simulation, Mej/Mtot can be scaled by QR/QRD *, and is independent of numerical resolution, impact velocity and target size. (3) This similarity law for QR/QRD * is confirmed for a wide range of specific impact energy. We also derive a semi-analytic formula for QRD * based on the similarity law and the crater scaling law. We find that the semi-analytic formula for the case with a non-porous object is consistent with numerical results.
AB - In most of the previous studies related to collisional disruption of planetesimals in the gravity regime, Smoothed Particle Hydrodynamics (SPH) simulations have been used. On the other hand, impact simulations using grid-based hydrodynamic code have not been sufficiently performed. In the present study, we execute impact simulations in the gravity regime using the shock-physics code iSALE, which is a grid-based Eulerian hydrocode. We examine the dependence of the critical specific impact energy QRD * on impact conditions for a wide range of specific impact energy (QR) from disruptive collisions to erosive collisions, and compare our results with previous studies. We find that collision outcomes of the iSALE simulation agree well with those of the SPH simulation. Detailed analysis mainly gives three results. (1) The value of QRD * depends on numerical resolution, and is close to convergence with increasing numerical resolution. The difference in converged value of QRD * between the iSALE code and the SPH code is within 30%. (2) Ejected mass normalized by total mass (Mej/Mtot) generally depends on various impact conditions. However, when QR is normalized by QRD * that is calculated for each impact simulation, Mej/Mtot can be scaled by QR/QRD *, and is independent of numerical resolution, impact velocity and target size. (3) This similarity law for QR/QRD * is confirmed for a wide range of specific impact energy. We also derive a semi-analytic formula for QRD * based on the similarity law and the crater scaling law. We find that the semi-analytic formula for the case with a non-porous object is consistent with numerical results.
KW - Cratering
KW - Impact processes
KW - Planetary formation
UR - http://www.scopus.com/inward/record.url?scp=85048289177&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85048289177&partnerID=8YFLogxK
U2 - 10.1016/j.icarus.2018.05.027
DO - 10.1016/j.icarus.2018.05.027
M3 - Article
AN - SCOPUS:85048289177
VL - 314
SP - 121
EP - 132
JO - Icarus
JF - Icarus
SN - 0019-1035
ER -