行星增生过程中金属-硅酸盐分离：中等尺度撞击诱发的铁乳化仍受限

《Icarus》：Metal-silicate segregation during planetary accretion: Limited iron emulsification through intermediate impacts

　　

夹带系数（entrainment coefficient）在高雷诺数（Reynolds number, Re）、韦伯数（Weber number, We）和邦德数（Bond number, Bo）条件下达到0.094的平台值。对于中等尺度撞击体，铁云（iron cloud）的几何分散仍然有限。研究人员的模型表明，初始撞击后碎片化（fragmentation）会持续更长时间。

研究背景与问题：行星核形成阶段，分异撞击体与原始地球碰撞形成岩浆洋（magma ocean, MO），其深度和几何受撞击体尺寸、速度和角度控制。铁核穿过MO时，MO湍流会使其变形、破碎成铁云并促进金属-硅酸盐混合与平衡；夹带系数决定铁云下陷过程中卷入硅酸盐的量，从而影响地幔-地核最终成分。已有理论、数值和实验研究对夹带系数的估计存在差异，且对中等尺度撞击体在撞击后长期碎片化、铁乳化（iron emulsification）和几何分散程度仍不清楚；已有研究多将弗劳德数（Froude number, Fr）限定于撞击瞬间，未约束撞击后长期碎片化过程，因此需要开展碰撞后铁核与MO相互作用的显式模拟。研究人员采用不混溶多相流模型，考察撞击核直径和长宽比（aspect ratio）对夹带系数及铁-硅酸盐混合的影响，以约束行星分异过程中的金属-硅酸盐分离和化学平衡。研究意义在于提供夹带系数基准和初始乳化阶段约束，支持MO/核化学模型将夹带系数视为常量，并提示真实MO中碎片化可持续至更小尺度、延长平衡窗口。

方法概要：研究采用二维D₂Q₉ Rothman–Keller多相晶格玻尔兹曼方法（Rothman–Keller multiphase lattice Boltzmann method, RKLBM），以recoloring算法维持尖锐相界面并解析铁核向碎片铁体的拓扑转变；设置铁与硅酸盐两相密度、黏度及表面张力，采用侧向周期、上下自由滑移边界。模拟直径152、252、352、452、552 km的撞击核，长宽比0.75–5、步长0.25，并扩展至1100 km；模型Re约3.7×10³–1.6×10⁵，We 6.68–1420，Bo 393–20616。未纳入热对流，通过热参数标定组分重力；无样本队列来源。

研究结果：
3.1 夹带系数：通过跟踪552 km等撞击核在MO中的演化并利用形态学包络计算面积A和周长p，按夹带假设拟合dA/dz=pα，得到α。结果显示撞击核先形成薄片、韧带并破碎；长宽比大者横向扩展更宽、碎片更多；小撞击核受背景流影响较大、拟合不确定性更高。
3.2 夹带系数与无量纲数：通过固定其他参数只增大D，以及固定D而改变域尺寸、重力和黏度两类模拟，发现α随D和无量纲数增大；当D≥452 km（Re≥1.9×10⁴、We≥5.9×10²、Bo≥3.5×10³、D/L_x≥7.5×10^?2、D/L_z≥1.51×10^?1）进入平台，18次运行平均α=0.094±0.011。研究人员指出，Re已足够高，但Bo和We仍低于真实MO多个量级，α是否最终对Bo、We平台仍需更高Bo、We研究。
3.3 碎片化：通过碎片半径时间序列、碎片数目和gamma分布拟合，发现撞击核在到达MO底部前逐步破碎至网格可解析最小尺度（半径可低于10 km），碎片尺寸呈分布而非单一特征半径；随D增大gamma形状参数k增大、尺度参数θ减小。模型仅捕捉片状/韧带状初级失稳，未解析毛细破碎；碎片下界由网格而非表面张力控制。不同D的碎片起始时间相近，且碎片化在初始撞击后持续，非单次事件。
3.4 碎片尺寸与碎片化时间的分辨率依赖性：通过2000²至8000²网格分辨率重复552 km、长宽比5的模拟，发现碎片数和平均半径依赖分辨率，平均半径约24个网格单元；但碎片化起始在最高两档分辨率收敛，后期破碎阶段无系统分辨率依赖。结论是破碎尺度为数值尺度，但破碎时序在测试范围内较稳健。
3.5 铁-硅酸盐界面分布分析：通过原始密度场和混合熵（mixing entropy, ME）指标分析界面像素和铁云空间分散，发现所有尺寸和长宽比下ME均较低（<0.1），铁云几何分散有限；较大撞击体峰值ME略高且混合状态分布更宽，长宽比影响随D增大而减弱。最大混合集中在MO中下部；界面像素比例随D增大而增加。该结果与大型撞击体不完全乳化的已有认识定性一致，但ME是网格尺度几何代理，不等同原子尺度热力学混合。

讨论部分总结：模型在较高Re下再现铁核破碎为铁云并持续碎片化，但未解析高效金属-硅酸盐平衡的小尺度子区间；其大尺度拉伸-减薄过程在高Re湍流中具有普适性。固定直径时，较长宽比撞击体横向分散更强、碎片更多。与实验相比，模型在更低Fr（0.085–0.67）下仍显示广泛碎片化，说明初始撞击动量并非唯一驱动，浮力驱动不稳定性可在低能条件下触发并维持碎片化；实验报道的Fr阈值和iSALE模拟与模型差异也反映了撞击后长期过程未被前者约束。模型α（0.00302–0.12426）低于实验值0.2–0.3，作者归因于实验管状几何与长宽比、MO域长宽比、未模拟超声速可压缩效应及二维几何抑制三维涡拉伸和湍流级联。D≥452 km后α平台支持以往将夹带系数视为常数的地幔-地核化学计算；真实MO的Re、We、Bo更高，α=0.094±0.011可适用于全尺寸范围。更细碎片沉降更慢、停留时间更长，有利于平衡；但研究未模拟撞击动力学、坑/射流、化学扩散和冲击加热。局限包括未解析毛细破碎、二维和未纳入热效应。

研究结论（翻译）：研究人员考查了分异撞击体铁核与约3000 km厚完全熔融MO相互作用的碰撞后时间演化，先模拟直径152–552 km（间隔100 km）、长宽比0.75–5（间隔0.25）且无热效应的铁核。模型显示，MO湍流强烈改造撞击核形态，使其破碎成铁云并在MO中分散；碎片化逐步发生而非单次事件，铁被拉伸减薄为细韧带，并在厚度接近网格尺度时分离。因此，表观碎片化级联由分辨率而非毛细破碎设定，模型不能将其与剥离物质的渐进侵蚀区分开。将模拟扩展至最大1100 km后，研究人员在充分发展的湍流区确定α=0.094±0.011，并将其作为MO条件的代表值。结果表明，细长撞击体显著增加铁-硅酸盐界面面积，为后续碎片再破碎创造条件；但显式追踪这些亚尺度再破碎超出当前分辨率。利用ME指标，研究人员发现所有考虑尺寸（包括中等范围和50 km、1100 km外围案例）的铁云分散均有限，与既有理论、数值和实验工作定性一致。结果还表明MO中下部铁云分散更宽，有利于金属-硅酸盐平衡和进一步碎片化，但捕捉这些子网格尺度破碎仍超出当前分辨率极限。