[论文解读] Well-balanced and asymptotic preserving schemes for kinetic models

Casimir Emako, Min Tang|arXiv (Cornell University)|Mar 10, 2016

Mathematical Biology Tumor Growth参考文献 17被引用 4

一句话总结

本文提出了一种通用框架，用于设计在两个参数区间下均具备良好平衡（WB）与渐近保持（AP）特性的数值格式：当 ε → 0 时的扩散极限，以及 ε = O(1) 时的稳态行为。该方法结合了 AP 时间演化步骤与局部稳态求解，以实现 WB 性质，成功在 ε 的全范围内以统一精度捕捉渐近极限与稳态解。

ABSTRACT

In this paper, we propose a general framework for designing numerical schemes that have both well-balanced (WB) and asymptotic preserving (AP) properties, for various kinds of kinetic models. We are interested in two different parameter regimes, 1) When the ratio between the mean free path and the characteristic macroscopic length $ε$ tends to zero, the density can be described by (advection) diffusion type (linear or nonlinear) macroscopic models; 2) When $ε$ = O(1), the models behave like hyperbolic equations with source terms and we are interested in their steady states. We apply the framework to three different kinetic models: neutron transport equation and its diffusion limit, the transport equation for chemotaxis and its Keller-Segel limit, and grey radiative transfer equation and its nonlinear diffusion limit. Numerical examples are given to demonstrate the properties of the schemes.

研究动机与目标

• 开发一种统一的数值框架，使动能方程同时满足良好平衡（WB）与渐近保持（AP）特性。
• 解决在两个不同参数区间中准确模拟动能模型的挑战：扩散极限（ε ≪ 1）与刚性源项区间（ε = O(1)）。
• 确保数值格式在 ε 的所有取值下均保持稳定与准确，尤其在标准格式失效的小 ε 区域。
• 在 ε → 0 时，于离散层面正确保持宏观极限（如 Keller-Segel 方程或非线性扩散方程）。
• 通过准确捕捉具有刚性源项的双曲系统的稳态解，维持长时间模拟中的高精度。

提出的方法

• 该框架采用两步法：首先，利用 AP 格式求解时间依赖的动能方程，生成预测解。
• 其次，利用预测解在每个计算单元内设定局部稳态问题的流入边界条件。
• 通过专为稳态方程设计的 AP 格式求解该稳态问题，确保与宏观极限的一致性。
• 将所得单元内流出量用于校正数值通量，从而在离散层面强制实现良好平衡性质。
• 该方法应用于三个模型：中子输运、趋化（Othmer-Dunbar-Alt）与灰辐射输运，提供了详细的离散化方案。
• 该方法具有灵活性，可兼容任何现有的 AP 格式用于时间演化与稳态问题，支持高效实现。

实验结果

研究问题

• RQ1 是否能通过单一数值框架，在不同参数区间下同时实现动能模型的优良平衡与渐近保持特性？
• RQ2 如何使离散格式在 ε → 0 时保持正确的宏观极限（如扩散或非线性扩散），而与 ε 的大小无关？
• RQ3 当 ε = O(1) 时，该格式能否准确捕捉具有刚性源项的双曲系统之稳态解？
• RQ4 该格式在网格尺寸 Δx 与参数 ε 上的收敛行为如何？其收敛性是否在 ε 的全范围内保持一致？
• RQ5 如何在不依赖非线性系统迭代求解器的前提下，在离散层面强制实现良好平衡性质？

主要发现

• 所提格式在所有 ε 取值下均表现出与网格尺寸 Δx 的一致收敛性，证实其渐近保持（AP）特性，如图 7 所示，收敛阶保持稳定。
• 当 ε → 0 时，格式正确捕捉了极限行为：辐射强度 Iε 趋近于 ψ⁰/2，其中 ψ⁰ = ac(T⁰)⁴，ε 的收敛阶接近 0.6，如图 8 所示。
• 当 ε = 1 时，格式成功捕捉了精确稳态解 I(x,v) = |v|x 与 ρ = x，残差在 T = 1 时减小并趋于稳定，证实了良好平衡（WB）性质（图 9 与图 10）。
• 数值结果表明，该格式在扩散极限（ε ≪ 1）与刚性源项区间（ε = O(1)）下均保持精度与稳定性，且随着 ε 减小，性能无退化。
• 该框架计算高效，因其允许在预测步骤中使用显式或简单隐式 AP 格式，并避免对非线性系统使用昂贵的牛顿型求解器。
• 该方法具有通用性与可扩展性，可推广至其他动能模型，包括具有松弛项的双曲系统，在多尺度模拟中具有广阔应用前景。