[论文解读] Coarse-grained analysis of a lattice Boltzmann model for planar streamer fronts

Wim Vanroose, Giovanni Samaey+1 more|ArXiv.org|Jan 2, 2007

Lattice Boltzmann Simulation Studies参考文献 3被引用 4

一句话总结

本文提出一种粗粒度时间推进方法，用于在包含速度依赖电离反应的格子Boltzmann模型中计算平面流光前驱的行进波解。通过结合短时模拟与提升和限制步骤，该方法计算固定点解，从而获得最小物理波速，为具有陡峭梯度的模型提供了一种稳健的替代传统Chapman-Enskog分析的方法。

ABSTRACT

We study the traveling wave solutions of a lattice Boltzmann model for the planar streamer fronts that appear in the transport of electrons through a gas in a strong electrical field. To mimic the physical properties of the impact ionization reaction, we introduce a reaction matrix containing reaction rates that depend on the electron velocities. Via a Chapman--Enskog expansion, one is able to find only a rough approximation for a macroscopic evolution law that describes the traveling wave solution. We propose to compute these solutions with the help of a coarse-grained time-stepper, which is an effective evolution law for the macroscopic fields that only uses appropriately initialized simulations of the lattice Boltzmann model over short time intervals. The traveling wave solution is found as a fixed point of the sequential application of the coarse-grained time-stepper and a shift-back operator. The fixed point is then computed with a Newton-Krylov Solver. We compare the resulting solutions with those of the approximate PDE model, and propose a method to find the minimal physical wave speed.

研究动机与目标

• 通过基于速度分辨的格子Boltzmann模型引入微观电子撞击电离物理机制，扩展最小流光模型。
• 解决传统Chapman-Enskog分析在捕捉电子密度和电场中陡峭梯度方面的局限性。
• 开发并验证一种计算方法，用于在难以推导精确宏观PDE的非理想动力学模型中求解行进波解。
• 利用粗粒度时间推进方法与Newton-Krylov求解器，数值确定流光前缘的最小物理波速。
• 将粗粒度方法与通过Chapman-Enskog展开导出的近似PDE模型进行比较。

提出的方法

• 构建一个具有多速度的格子Boltzmann模型以模拟电子动力学，基于微观电离截面引入速度依赖的反应速率。
• 粗粒度时间推进方法（CGTS）定义为一系列步骤：提升（从宏观场初始化分布函数）、模拟（在粗时间步长ΔT内推进LBM）和限制（提取宏观场）。
• 通过将CGTS与回移算子结合，将行进波解表示为固定点，以考虑波的传播。
• 采用Newton-Krylov求解器计算固定点，实现对行进波形的收敛。
• 该方法避免显式推导宏观PDE，即使在存在陡峭梯度时也能实现精确模拟。
• 通过将计算波速与通过Chapman-Enskog展开导出的近似PDE模型结果对比，验证了该方法的有效性。

实验结果

研究问题

• RQ1 粗粒度时间推进方法能否准确计算具有速度依赖电离反应的格子Boltzmann模型中的行进波解？
• RQ2 最小物理波速如何依赖于微观模型中电离截面的强度？
• RQ3 与传统Chapman-Enskog展开相比，粗粒度时间推进方法在预测临界波速方面表现如何？
• RQ4 该方法能否定位共向传播和反向传播的流光前缘，包括正电荷前缘？
• RQ5 高阶提升过程对粗粒度时间推进方法的精度有何影响？

主要发现

• 粗粒度时间推进方法成功计算了具有速度依赖电离反应的格子Boltzmann模型中平面流光前缘的行进波解。
• 该方法通过结合CGTS与回移算子的固定点计算，实现了对最小物理波速的数值计算。
• 计算得到的波速依赖于电离截面的强度，且该依赖关系可通过粗粒度方法定量获取。
• 粗粒度方法的结果与通过Chapman-Enskog展开导出的近似PDE模型结果一致，验证了该方法的可靠性。
• 该方法足够稳健，能够定位共向传播和反向传播的流光前缘，包括正电荷前缘。
• 高阶提升过程有望通过更精确地近似初始分布函数，减少模拟步数并提高收敛性，从而提升精度。