[论文解读] Simulating dust grain-radiation coupling on a moving mesh

R. McKinnon, Rahul Kannan+4 more|arXiv (Cornell University)|Dec 5, 2019

Astrophysics and Star Formation Studies参考文献 2被引用 6

一句话总结

本文提出了一种新型混合方法，应用于 arepo-rt 辐射流体动力学代码中，将实时尘埃模拟粒子与非结构化移动网格相结合，用于辐射传输。通过自洽地计算尘埃消光截面、辐射压力、光子吸收和热辐射，该模型实现了真实尘埃-气体-辐射之间的能量与动量交换，无需假设局部热动平衡，并通过解析解验证了速度与温度演化的准确性。

ABSTRACT

We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code AREPO-RT, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically-symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in AREPO-RT, which treats UV and optical photons as single-scattered and IR photons as multi-scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

研究动机与目标

• 开发一种在辐射流体动力学模拟中实现尘埃-辐射耦合的物理准确、自洽的模型。
• 克服子网格模型的局限性，后者假设固定的尘埃-气体比或局部热动平衡。
• 通过将尘埃视为响应辐射压力和阻力的实时模拟粒子，实现动态尘埃行为。
• 改进多频辐射场中尘埃加热、热辐射和光子吸收的建模。
• 通过涉及辐射驱动尘埃和薄层中光子吸收的解析测试问题，对框架进行验证。

提出的方法

• 尘埃被建模为代表颗粒群体的实时模拟粒子，每个粒子具有质量、速度以及与颗粒尺寸相关的截面。
• 通过将粒子的尘埃截面插值到气体单元，计算辐射压力力，从而实现通过阻力将动量传递给气体。
• 通过核光滑插值将尘埃消光截面映射到气体单元，计算光子吸收，并自洽地更新辐射场。
• arepo-rt 中的 M1 闭合方案处理多频辐射，可区分单次散射的紫外/光学光子与多次散射的红外光子。
• 尘埃和气体温度独立演化，通过碰撞过程实现能量交换，尘埃通过吸收辐射和热辐射实现加热与冷却。
• 该方法采用混合粒子-单元方法：尘埃粒子随流体动力学运动，而辐射在非结构化移动网格上求解。

实验结果

研究问题

• RQ1 混合粒子-单元方法在球对称辐射场中模拟尘埃颗粒所受辐射压力的精度如何？
• RQ2 在移动网格辐射流体动力学代码中，能否在不假设固定消光截面的前提下，自洽地模拟尘埃对光子的吸收？
• RQ3 非平衡尘埃-气体能量交换的引入如何影响尘埃和气体的动力学与热演化？
• RQ4 在辐射主导环境中，颗粒尺寸分布的变化对辐射压力和动量传递有何影响？
• RQ5 该模型在辐射压力和吸收作用下，对尘埃加速及与气体速度同步的解析解重现效果如何？

主要发现

• 该模型成功重现了恒定光度源驱动下尘埃加速的解析解，初始阶段由于辐射压力导致尘埃与气体速度分离，随后通过阻力实现速度同步。
• 在无光子吸收的情况下，尘埃与气体速度持续增加，速度差约为 108 km s⁻¹，与解析预测一致。
• 当包含光子吸收时，尘埃与气体速度在约 44 km s⁻¹ 处达到平衡，与动量传递和能量吸收的解析解一致。
• 该模型准确捕捉了通过薄尘埃层的光子吸收，光学深度依赖关系与解析预期相符。
• 未假设尘埃与气体处于局部热动平衡；分别演化其温度，通过碰撞过程实现能量交换。
• 该框架基于局部颗粒丰度自洽计算尘埃消光截面，避免了对尘埃-气体比或分布的假设。