[论文解读] 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 more|arXiv (Cornell University)|Dec 30, 2013

Astro and Planetary Science参考文献 1被引用 3

一句话总结

本研究对比了光滑粒子动力学（SPH）模拟中固体体（弹塑性）模型与流体动力学模型在Ceres质量星子与冰壳碰撞中的表现，表明材料强度显著增加了碎片化和水损失，尤其是在低至中等能量碰撞中，同时实现了可检测的水转移，挑战了早期行星系统模型中完美合并的假设。

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体模拟早期行星系统中使用计算成本较高的固体模型的可行性与必要性。

提出的方法

• 采用自定义SPH代码，实现无强度流体动力学模型与完整弹塑性材料模型（含脆性破坏），基于Libersky & Petschek（1991）与Benz & Asphaug（1994, 1995）的理论。
• 模拟Ceres质量目标体（表面含30 wt%水冰）与固体岩石撞击体在1至6倍逃逸速度范围内的碰撞。
• 改变撞击参数，从正碰至近距离飞越，以探索广泛的碰撞几何与能量范围。
• 追踪物质演化过程，包括质量损失、碎片形成以及幸存体与碎片云之间水含量的再分布。
• 模拟持续2000分钟以上，以评估最终系统状态及碎片的长期稳定性。
• 采用GPU加速SPH计算，实现相较于CPU实现约50倍的性能提升，支持大规模参数空间探索。

实验结果

研究问题

• RQ1 与无强度流体动力学模型相比，包含材料强度如何影响星子碰撞中的碎片化程度与碎片云扩散？
• RQ2 在低能碰撞中，固体模型相较于流体模型是否显著预测更高的水损失，尤其是对冰质星子而言？
• RQ3 材料强度是否能实现可检测的水（冰）在碰撞的Ceres大小天体之间的转移，特别是在“击中逃逸”或部分吸积情形下？
• RQ4 是否存在特定碰撞模式（如低能或近中心撞击），使得固体模型与流体模型产生显著不同的结果？
• RQ5 在大规模n体模拟早期行星系统中，使用固体体模型的计算成本与收益之间存在何种权衡？

主要发现

• 固体模型产生的碎片数量约为流体动力学模型的10倍，且碎片云更广泛、更接近粉尘状。
• 固体模型中幸存体的平均质量更低，源于增强的碎片化与材料损失。
• 水损失在固体模型中始终更高，且在低至中等能量碰撞中差异最为显著。
• 模拟结果明确显示了天体之间水冰的转移，尤其在“击中逃逸”与部分吸积情景中，这是流体动力学模型无法捕捉的特征。
• 低能与近中心碰撞表现出固体模型与流体模型之间最显著的差异，表明这些情形需依赖固体体物理才能实现准确模拟。
• 结果表明，为真实模拟早期行星系统中水输送与小行星家族形成，尤其当物质转移与碎片化为关键因素时，固体体模型是必要的。