[논문 리뷰] The role of material strength in collisions -- Comparing solid body and hydrodynamic physics for simulating collisions of planetesimals with icy shells

Thomas I. Maindl, R. Dvořák외 2명|arXiv (Cornell University)|2013. 12. 30.

Astro and Planetary Science참고 문헌 1인용 수 3

한 줄 요약

이 연구는 Ceres 질량의 원형행성형체가 얼음 껍질과 충돌할 때, 탄성-플라스틱 고체 모델과 유체역학 모델을 스플래터드 파티클 유체역학(SPH) 시뮬레이션에서 비교하여, 물질 강도가 분열과 물질 손실을 크게 증가시킨다는 것을 보여주며, 특히 저에너지에서 중간 에너지 충돌에서 그러한 영향이 두드러지며, 물체 간에 검출 가능한 물질 이동을 가능하게 하여 초기 행성계 모델에서의 완전한 융합 가정을 도전한다.

ABSTRACT

Context. We investigate the effects of including material strength in multi-material planetesimal collisions. Aims. The differences between strengthless material models and including the full elasto-plastic model for solid bodies with brittle failure and fragmentation when treating collisions of asteroid-sized bodies as they occur frequently in early planetary systems are demonstrated. Methods. We study impacts of bodies of Ceres-mass with a solid rock impactor and a target with 30 weight-% water content as surface ice.} The initial impact velocities and impact parameters are varied between the escape velocity $v_\mathrm{esc}$ to about 6 $v_\mathrm{esc}$ and from head-on collisions to close fly-bys, respectively. We simulate the collisions using our own SPH code utilizing both strengthless material and the full elasto-plastic material model including brittle failure. Results. One of the most prominent differences is the higher degree of fragmentation and shattered debris clouds in the solid model. In most collision scenarios however, the final outcomes are very similar and differ primarily by the about one order of magnitude higher degree of fragmentation in the solid case. Also, the survivors tend to be of less mass in the solid case which also predicts a higher water loss than the strengthless hydro model. This may be an effect of the relatively low-energy impacts that cannot destroy the solid material instantly. As opposed to giant impacts we also observe an indication that some water ice gets transferred between the bodies.

연구 동기 및 목표

• 탄성-플라스틱 고체 모델을 사용하여 Ceres 질량의 천체에서 물질 강도가 충돌 결과에 미치는 영향을 평가함. 특히 얼음 표면을 가진 천체에 초점을 맞춤.
• 초기 행성계 시뮬레이션에서 흔히 사용되는 물질 강도가 없는 유체역학 모델의 가정을 도전함. 이는 물질 이동과 분열을 단순화할 수 있음.
• 탄성-플라스틱 거동와 파손 파손 메커니즘을 포함한 고체 모델이, 유체역학 모델과 비교하여 잔해 분포, 질량 손실, 물질 이동 등의 결과에 어떤 영향을 미치는지 조사함.
• 특히 저에너지에서 중간 에너지 충돌과 중심 충돌에 가까운 경우와 같이 고체 모델이 정확한 충돌 결과 예측에 필수적인 조건을 규명함.
• 초기 행성계의 대규모 n체 시뮬레이션에서 계산 비용이 높은 고체 모델을 사용하는 것이 실현 가능하고 필수적인지 평가함.

제안 방법

• Libersky & Petschek(1991), Benz & Asphaug(1994, 1995) 기반으로, 물질 강도가 없는 유체역학 모델과 완전한 탄성-플라스틱 물질 모델을 구현한 고유의 SPH 코드를 사용함.
• Ceres 질량의 목표체(표면에 30 wt%의 물의 얼음 포함)와 고체 암석 충돌체 간의 충돌을, 탈출 속도의 1배에서 6배의 속도로 시뮬레이션함.
• 정면 충돌에서 가까운 프라이 베이를 포함한 다양한 충돌 기하학과 에너지를 탐색하기 위해 충돌 매개변수를 다양화함.
• 생존체와 잔해 구름 내에서 질량 손실, 조각 형성, 물 함량 재분포 등의 물질 변화를 추적함.
• 충돌 후 약 2000분 동안 시뮬레이션을 수행하여 최종 시스템 상태와 잔해의 장기적 안정성을 평가함.
• 성능 향상을 위해 GPU 가속을 사용하여 CPU 기반 구현 대비 약 50배의 속도 향상을 달성함. 이는 대규모 파arameter 공간 탐색을 가능하게 함.

실험 결과

연구 질문

• RQ1 물질 강도를 고려한 고체 모델이 물질 강도가 없는 유체역학 모델과 비교해 원형행성형체 충돌에서 분열 정도와 잔해 구름 확산 정도에 어떤 영향을 미치는가?
• RQ2 저에너지 충돌에서 특히 고체 모델이 유체역학 모델보다 얼음으로 이루어진 원형행성형체에서 더 높은 물질 손실을 예측하는 정도는 어느 정도인가?
• RQ3 물질 강도가 Ceres 크기의 천체 간에 검출 가능한 물(ice) 이동을 가능하게 하는가? 특히 히트-앤드-런 또는 부분적 응집 영역에서 그러한 현상이 나타나는가?
• RQ4 저에너지 충돌이나 중심에 가까운 충돌 등 특정 충돌 영역에서 고체 모델이 유체 모델과 유의미하게 다른 결과를 도출하는가?
• RQ5 초기 행성계의 대규모 n체 시뮬레이션에서 고체 모델을 사용할 경우의 계산 비용 대비 성능 향상의 균형은 어떠한가?

주요 결과

• 고체 모델은 유체역학 모델 대비 약 10배 더 많은 조각과 더 광범위하고 먼지 같은 잔해 구름을 생성함.
• 고체 모델에서 생존체의 평균 질량은 더 낮은 편이 되며, 이는 강화된 분열과 물질 손실 때문임.
• 물질 손실은 고체 모델에서 일관되게 더 높으며, 특히 저에너지에서 중간 에너지 충돌에서 그 차이가 가장 두드러짐.
• 시뮬레이션 결과에서 히트-앤드-런 및 부분적 응집 상황에서 천체 간에 물의 얼음 이동이 명확하게 관측됨. 이는 유체역학 모델에서는 포착되지 않은 특성임.
• 저에너지 및 중심에 가까운 충돌에서 고체 모델과 유체 모델 간의 결과 차이가 가장 두드러지며, 이러한 영역에서는 고체 모델 물리학이 정확한 시뮬레이션을 위해 필수적임을 시사함.
• 결과적으로 고체 모델은 초기 행성계에서 물질 공급과 소행성 군 형성의 현실적인 모델링을 위해 필수적이며, 특히 물질 이동과 분열이 핵심 고려 사항일 경우 그러한 모델이 필요함을 시사함.