[论文解读] Implicit-explicit, realizability-preserving first-order scheme for moment models with Lipschitz-continuous source terms

Florian Schneider|arXiv (Cornell University)|Nov 4, 2016

Gas Dynamics and Kinetic Theory参考文献 30被引用 3

一句话总结

本文提出了一种隐式-显式（IMEX）的、保持可实现性的的一阶数值格式，用于具有利普希茨连续源项的矩模型，特别应用于最小熵（M_N）模型。通过隐式处理刚性源项并显式处理非刚性通量，该格式消除了对物理参数的时间步长限制，从而实现稳定且高效的模拟，同时确保矩向量保持可实现性——即与非负的动能分布一致。

ABSTRACT

We derive an implicit-explicit (IMEX), realizability-preserving first-order scheme for moment models with Lipschitz-continuous source terms. In contrast to fully-explicit schemes the time step does not depend on the physical parameters, removing the stiffness from the system. Furthermore, a wider class of collision operators (e.g. the Laplace-Beltrami operator) can be used. The derived scheme is applied to minimum-entropy models.

研究动机与目标

• 开发一种数值格式，以在刚性源项下保持最小熵（M_N）模型中矩向量的可实现性。
• 克服完全显式格式的时间步长限制，这些限制依赖于吸收系数和散射系数等物理参数。
• 使更广泛的碰撞算子类（包括拉普拉斯-贝尔特拉米算子）能够应用于矩模型。
• 确保格式在通过隐式处理刚性源项的同时保持一阶精度，且维持可实现性。
• 为未来工作提供基础，将高阶 IMEX 格式扩展至完全显式、高阶的动能和间断伽辽金方法。

提出的方法

• 采用隐式-显式（IMEX）时间离散化，将刚性源项隐式处理，非刚性双曲通量显式处理。
• 对简化的、空间均匀的矩系统应用保持可实现性的隐式离散化，确保重构的动能密度非负。
• 在每个时空单元内局部求解隐式系统，实现并行化和计算效率。
• 依赖于碰撞算子满足质量守恒、局部熵耗散以及时间推进下非负性保持的假设。
• 采用拉普拉斯-贝尔特拉米算子作为满足所需假设并实现更广泛应用的代表性碰撞算子。
• 通过理论分析和数值测试（包括构造解和基准平面源问题）验证了可实现性的保持。

实验结果

研究问题

• RQ1 能否构建一种 IMEX 格式，使其在具有刚性、利普希茨连续源项的矩模型中保持可实现性？
• RQ2 对源项进行隐式处理是否能消除完全显式格式中固有的时间步长限制，同时不损失稳定性？
• RQ3 该格式能否应用于更广泛的碰撞算子类（如拉普拉斯-贝尔特拉米算子），同时保持可实现性？
• RQ4 在实际应用中，该格式能多接近可实现性边界，特别是对于高阶 M_N 模型？
• RQ5 当理论框架的假设在实际中仅被近似满足时，可实现性保持特性是否依然稳健？

主要发现

• 该 IMEX 格式在标准 CFL 条件下保持可实现性，即使在强吸收或强散射等刚性源项下亦然。
• 时间步长不再受物理参数限制，相比完全显式格式可采用更大、更高效的步长。
• 通过与构造解的收敛性测试验证了格式的一阶精度。
• 在平面源基准测试中，解始终接近可实现性边界，相对距离最小值分别为 0.0039（M₁）、2.2147×10⁻⁵（M₂）和 9.0981×10⁻⁷（M₃），表明随着阶数升高，矩问题求解难度增加。
• 相空间中 M₂ 和 M₃ 解的直方图显示，粒子集中在可实现集的下部，避免远离边界的区域，证实了该格式的鲁棒性。
• 即使在理论假设可能未完全成立的情况下（如使用拉普拉斯-贝尔特拉米算子的混合矩模型），该方法仍表现良好，表明其有潜力扩展至更弱假设或更广泛类别的模型。