[论文解读] Analysis of an asymptotic preserving scheme for stochastic linear kinetic equations in the diffusion limit

Nathalie Ayi, Erwan Faou|arXiv (Cornell University)|Mar 16, 2018

Gas Dynamics and Kinetic Theory参考文献 18被引用 9

一句话总结

本文提出了一种用于带乘性噪声的随机线性动能方程的新型渐近保持（AP）格式，通过微观-宏观分解将分布函数分解为宏观与微观分量。该格式在扩散模型和一般情况下均被证明在显式CFL条件下的均匀稳定性，确保在动能与扩散区之间的精确高效模拟，即使当平均自由程 ε → 0 时亦然。

ABSTRACT

We present an asymptotic preserving scheme based on a micro-macro decomposition for stochastic linear transport equations in kinetic and diffusive regimes. We perfom a mathematical analysis and prove that the scheme is uniformly stable with respect to the mean free path of the particles in the simple telegraph model and in the general case. We present several numerical tests which validate our scheme.

研究动机与目标

• 开发一种数值格式，使其在带乘性噪声的随机线性输运方程中，于动能与扩散区均保持稳定与准确。
• 将渐近保持（AP）框架扩展至带乘性维纳噪声的随机偏微分方程，这是AP文献中的新贡献。
• 证明关于尺度参数 ε（平均自由程）的均匀稳定性，确保在 ε → 0 时的鲁棒性。
• 通过与显式格式和Crank-Nicolson格式在不同 ε 值下的比较，对格式进行数值验证。

提出的方法

• 基于碰撞算子的守恒性与平衡性质，将分布函数 f 分解为宏观分量 ρ 和微观分量 g，实现微观-宏观分解。
• 对耦合的微观-宏观系统采用半隐式时间离散化，确保在刚性区域的稳定性。
• 通过柱状维纳过程 dWt 引入乘性噪声，通过带有自伴算子 Q 的随机积分建模，其中 Q 满足可和性条件。
• 利用伊藤微积分与随机格拉沃尔德型方法推导能量估计，建立在CFL条件下的稳定性。
• 实现一种随机AP格式（SMM），在极限 ε → 0 时退化为扩散求解器，保持渐近行为。
• 采用具有周期性边界条件的有限差分空间离散化，并通过100次随机实现计算系综平均。

实验结果

研究问题

• RQ1 能否为带乘性噪声的随机线性动能方程构造一种渐近保持格式？
• RQ2 所提出的格式在扩散模型和一般情况下是否对平均自由程 ε 具有均匀稳定性？
• RQ3 该格式在动能区（ε = 1）与扩散区（ε = 10⁻²）是否均保持准确与稳定？
• RQ4 在一般情况下，确保稳定性的CFL条件是什么？其如何依赖于 ε 与空间分辨率？
• RQ5 在平均密度演化方面，SMM格式与显式格式和Crank-Nicolson格式相比表现如何？

主要发现

• 在扩散模型中，该格式对 ε 具有均匀稳定性，稳定性通过能量估计与伊藤微积分证明。
• 在一般情况下，稳定性在显式CFL条件下得以建立：Δt ≤ (σ̃Δx²)/(2(2+ε)) + (εΔx)/(2+ε)，确保能量范数的有界增长。
• 能量估计表明，宏观密度 ρ 与微观涨落 g 的期望平方 L² 范数随时间有界，其增长受依赖于噪声系数的常数控制。
• 数值测试证实，SMM 格式在动能区（ε=1）与扩散区（ε=10⁻²）均能准确捕捉平均密度 ρ，与参考解一致。
• 该格式在极限 ε → 0 时成功退化为扩散求解器，确认其在随机设置下的渐近保持性质。
• 该方法在有限时间区间内保持稳定，但该证明不适用于无限时间，原因在于随机能量估计的性质。