[论文解读] On deterministic approximation of the Boltzmann equation in a bounded domain

Francis Filbet|arXiv (Cornell University)|Jun 6, 2011

Gas Dynamics and Kinetic Theory参考文献 25被引用 27

一句话总结

本文提出了一种完全确定性的高阶数值方法，用于在具有周期性、镜面反射和扩散性（Maxwell型）边界条件的有界区域内求解玻尔兹曼方程。该方法在空间上采用二阶有限体积格式，在碰撞项上采用傅里叶谱方法，实现了速度空间中的谱精度，从而能够在广泛的克努森数范围内实现高效模拟，包括流体动力学极限和稀薄流区域。

ABSTRACT

In this paper we present a fully deterministic method for the numerical solution to the Boltzmann equation of rarefied gas dynamics in a bounded domain for multi-scale problems. Periodic, specular reflection and diffusive boundary conditions are discussed and investigated numerically. The collision operator is treated by a Fourier approximation of the collision integral, which guarantees spectral accuracy in velocity with a computational cost of $N\\,\\log(N)$, where $N$ is the number of degree of freedom in velocity space. This algorithm is coupled with a second order finite volume scheme in space and a time discretization allowing to deal for rarefied regimes as well as their hydrodynamic limit. Finally, several numerical tests illustrate the efficiency and accuracy of the method for unsteady flows (Poiseuille flows, ghost effects, trend to equilibrium).

研究动机与目标

• 开发一种完全确定性的数值方法，用于求解有界区域内具有多尺度行为的时空依赖玻尔兹曼方程。
• 准确处理各种边界条件，包括周期性、镜面反射和扩散性（Maxwell型）条件。
• 通过结合二阶有限体积格式与谱方法，在空间和速度方向上均实现高阶精度。
• 通过解耦时间离散化，实现从稀薄流到流体动力学区域的全克努森数范围模拟。
• 将该方法与已知的渐近极限进行验证，特别是当连续介质模型失效时的鬼效应和趋向平衡的趋势。

提出的方法

• 对玻尔兹曼方程的输运项采用二阶有限体积格式进行空间离散化。
• 对碰撞积分应用傅里叶谱方法，确保在速度空间中具有谱精度，计算成本为 $N\log N$。
• 通过算子分裂方法将输运与碰撞步骤解耦，实现高效的时序积分与并行化。
• 采用时间离散化方法，确保在包括流体动力学极限在内的多尺度行为下保持稳定性和精度。
• 通过混合Maxwell型模型处理边界条件：比例 $\alpha$ 的粒子以麦克斯韦分布重新发射，其余 $1-\alpha$ 的粒子发生镜面反射。
• 通过边界处的质量守恒求解壁面弛豫系数 $\mu(t,x)$，确保物理一致性。

实验结果

研究问题

• RQ1 能否通过完全确定性的方法在有界区域内实现玻尔兹曼方程在空间和速度方向上的高阶精度？
• RQ2 该方法在小克努森数条件下的周期性温度驱动流中，能否准确捕捉鬼效应？
• RQ3 在流体动力学极限下，数值解是否收敛于渐近理论，而非经典热传导方程的解？
• RQ4 该方法能否使用单一算法和最小网格分辨率，高效模拟稀薄流和流体动力学区域？
• RQ5 与基于标准求积的碰撞求解器相比，速度方向的谱方法在精度和计算成本方面表现如何？

主要发现

• 该方法在空间上达到二阶精度，在速度方向上达到谱精度，仅需较少的速度网格点即可获得高分辨率结果。
• 泊肃叶流的数值模拟结果与理论预测和基准解高度一致。
• 在鬼效应问题中，随着克努森数 $\varepsilon$ 减小，温度场升高，收敛于 Sone 等人的渐近理论，而非热传导方程的解。
• 在 $\varepsilon \to 0$ 的极限下，平均速度场趋于零，与渐近理论一致，证实了无虚假流动的存在。
• 该方法成功捕捉了周期性壁面温度问题中的非平衡温度分布，等温线偏离了热传导方程的预测结果。
• 碰撞步骤的计算成本为 $N\log N$，在高分辨率速度网格下具有高效性，且解耦结构支持并行实现。