[论文解读] Numerical Simulation of Microflows using Hermite Spectral Methods

Zhicheng Hu, Zhenning Cai+1 more|arXiv (Cornell University)|Jul 17, 2018

Gas Dynamics and Kinetic Theory参考文献 28被引用 4

一句话总结

本论文提出了一种用于求解空间非均匀玻尔兹曼方程的埃尔米特谱方法，采用广义近似二次碰撞算子来处理任意分布函数。通过引入一种快速变换以高效离散化碰撞算子，该方法实现了高精度和高收敛性，计算复杂度为 O(M₀⁹ + M⁴)，在1D和2D微流基准测试中与DSMC结果高度一致。

ABSTRACT

We propose a Hermite spectral method for the spatially inhomogeneous Boltzmann equation. For the inverse-power-law model, we generalize an approximate quadratic collision operator defined in the normalized and dimensionless setting to an operator for arbitrary distribution functions. An efficient algorithm with a fast transform is introduced to discretize this new collision operator. The method is tested for one-dimensional benchmark microflow problems.

研究动机与目标

• 开发一种用于稀薄气体动力学中空间非均匀玻尔兹曼方程的确定性、高阶数值求解器。
• 将埃尔米特谱方法从归一化、空间均匀情况扩展至一般分布函数。
• 设计一种高效的碰撞算子算法，保持低阶矩的二次型形式，同时实现快速计算。
• 在基准微流问题上验证该方法，包括不同克努森数下的1D和2D驱动腔流。
• 在水动力和动态区域中，实现比线性模型或DSMC更高的精度和更平滑的结果。

提出的方法

• 提出一种基于正交多项式展开的埃尔米特谱方法，采用麦克斯韦加权埃尔米特多项式来表示分布函数。
• 将文献[29]中的近似二次碰撞算子推广至任意分布函数，仅对前几阶矩保持二次型。
• 引入一种快速变换算法，以高效计算碰撞算子，将计算复杂度降低至 O(M₀⁹ + M⁴)。
• 在结构化网格（100×100）上采用显式时间推进并自适应时间步长进行2D模拟，自由度参数为 M₀ = 10，M = 25 或 35。
• 边界条件处理存在吉布斯现象问题，建议未来采用滤波器或改进的边界处理方法。
• 采用基于CPU的实现方式，在Intel Xeon E5-2680 v4处理器（28个线程）上运行，测量总CPU时间和每步CPU时间以评估性能。

实验结果

研究问题

• RQ1 基于广义近似碰撞算子的埃尔米特谱方法是否能在求解微流中的空间非均匀玻尔兹曼方程时实现高精度和高效率？
• RQ2 所提出的快速变换算法如何在保持物理解一致性的同时降低碰撞算子的计算成本？
• RQ3 该方法在1D和2D微流基准测试（如驱动腔流）中，对DSMC结果的再现程度如何？
• RQ4 当改变埃尔米特模态数（M）和碰撞近似阶数（M₀）时，精度与计算成本之间的权衡如何？
• RQ5 如吉布斯现象等数值伪影如何影响边界附近的解，可采取何种改进措施？

主要发现

• 在2D驱动腔流中，对于 Kn = 0.1 和 Kn = 1.0 的情况，该埃尔米特谱方法与DSMC结果均表现出良好一致性。
• 对于 Kn = 0.1，采用 M = 25 和 M₀ = 10，对应 3,276 个自由度，总CPU时间为 6.93 × 10⁶ 秒。
• 对于 Kn = 1.0，采用 M = 35 和 M₀ = 10，对应 8,436 个自由度，总CPU时间为 1.07 × 10⁷ 秒。
• 每自由度的CPU时间为 Kn = 0.1 时为 1.07 × 10⁻⁵ 秒，Kn = 1.0 时为 5.30 × 10⁻⁶ 秒，表明具有良好的可扩展性。
• 在边界附近观察到差异，可能源于分布函数不连续引起的吉布斯现象，提示需引入滤波或改进边界条件。
• 该方法表现出 O(M₀⁹ + M⁴) 的时间复杂度，使其在复杂多维及多原子气体问题中具有良好的计算前景。