[论文解读] The Consistent Multi-fluid Advection method

T. Plewa, Ewald Müller|ArXiv.org|Jul 23, 1998

Ionosphere and magnetosphere dynamics被引用 13

一句话总结

本文提出了一种一致多流体平流（CMA）方法，这是一种保守的高阶Godunov型格式，通过强制质量分数归一化并利用改进的插值剖面减少数值扩散，从而精确平流多流体流动。该方法保持了守恒性和单调性，在流体动力学模拟中显著提升了组分分辨率，特别是在II型超新星模型中，由于扩散性降低，重元素丰度可改变数倍。

ABSTRACT

Simple modifications for higher-order Godunov-type difference schemes are presented which allow for accurate advection of multi-fluid flows in hydrodynamic simulations. The constraint that the sum of all mass fractions has to be equal to one in every computational zone throughout the simulation is fulfilled by renormalizing the mass fractions during the advection step. The proposed modification is appropriate for any difference scheme written in conservation form. Unlike other commonly used methods it does not violate the conservative character of the advection method. A new steepening mechanism, which is based on modification of interpolation profiles, is used to reduce numerical diffusion across composition discontinuities. Additional procedures are described, which are necessary to enforce monotonicity. Several numerical experiments are presented which demonstrate the capability of our Consistent Multi-fluid Advection (CMA) method in case of smooth and discontinuous distributions of fluid phases and under different hydrodynamic conditions. It is shown that due to the reduced diffusivity of the proposed scheme the abundance of some heavy elements obtained from hydrodynamic simulations of type II supernova explosions can change by a factor of a few in the most extreme cases.

研究动机与目标

• 解决多流体流体动力学模拟中的数值扩散和质量分数不一致问题。
• 在平流多个流体组分时，保持高阶Godunov格式的守恒特性。
• 在不违反守恒性或单调性的情况下，提高组分间断面的分辨率。
• 实现对天体物理环境（如II型超新星）中流体混合的精确建模。
• 开发一种可推广至任意守恒形式差分格式的一般性方法。

提出的方法

• 通过在平流步骤中进行重归一化，强制在每个计算单元中所有质量分数之和等于1。
• 通过修改插值剖面，创建一种新的陡化机制，以减少组分界面处的数值扩散。
• 该格式与任何以守恒形式书写的高阶Godunov型有限体积方法兼容。
• 实施了额外的算法程序以确保单调性，并防止解中出现非物理振荡。
• 该方法保持了原始格式的守恒特性，避免引入人工源项。
• 在一系列数值实验中应用该方法，以测试其在光滑和不连续流条件下的性能。

实验结果

研究问题

• RQ1 如何在保持组分间断面处精度的同时，以守恒方式执行多流体平流？
• RQ2 对高阶Godunov格式进行何种修改，可减少数值扩散而不破坏守恒性？
• RQ3 质量分数的重归一化对多流体模拟的稳定性和精度有何影响？
• RQ4 基于插值的新陡化机制在多大程度上提升了流体界面的分辨率？
• RQ5 减少的数值扩散对超新星模拟中重元素丰度预测有何影响？

主要发现

• CMA方法在整个模拟域内成功将质量分数归一化保持在机器精度水平。
• 改进的插值剖面显著减少了在尖锐组分间断面处的数值扩散。
• 该方法保持了底层Godunov格式的守恒特性，避免了人工源项的引入。
• 通过额外的算法程序强制实现单调性，防止出现非物理振荡。
• 在II型超新星模拟中，由于扩散性降低，极端情况下重元素丰度的变化可达数倍。
• 数值实验证实，该方法在包括光滑流和不连续流在内的多种流体动力学条件下均表现出稳健性能。