[论文解读] A Numerical Scheme for the Quantum Boltzmann Equation Efficient in the Fluid Regime

Francis Filbet, Jingwei Hu+1 more|arXiv (Cornell University)|Sep 17, 2010

Gas Dynamics and Kinetic Theory参考文献 17被引用 10

一句话总结

该论文提出了一种新颖的渐近保持（asymptotic-preserving）数值格式，通过使用经典BGK算子对量子碰撞算子进行惩罚，而非计算量巨大的量子BGK算子。利用经典麦克斯韦分布（以内部能量作为温度）与量子麦克斯韦分布在前五阶矩上一致的特性，该方法避免了对逸度和温度的复杂非线性系统求解，从而在所有Knudsen数范围内实现了稳定、高效的模拟，且在动能与流体区域均保持一致的精度。

ABSTRACT

Numerically solving the Boltzmann kinetic equations with the small Knudsen number is challenging due to the stiff nonlinear collision term. A class of asymptotic preserving schemes was introduced in [6] to handle this kind of problems. The idea is to penalize the stiff collision term by a BGK type operator. This method, however, encounters its own difficulty when applied to the quantum Boltzmann equation. To define the quantum Maxwellian (Bose-Einstein or Fermi- Dirac distribution) at each time step and every mesh point, one has to invert a nonlinear equation that connects the macroscopic quantity fugacity with density and internal energy. Setting a good initial guess for the iterative method is troublesome in most cases because of the complexity of the quantum functions (Bose-Einstein or Fermi-Dirac function). In this paper, we propose to penalize the quantum collision term by a 'classical' BGK operator instead of the quantum one. This is based on the observation that the classical Maxwellian, with the temperature replaced by the internal energy, has the same first five moments as the quantum Maxwellian. The scheme so designed avoids the aforementioned difficulty, and one can show that the density distribution is still driven toward the quantum equilibrium. Numerical results are present to illustrate the efficiency of the new scheme in both the hydrodynamic and kinetic regimes. We also develop a spectral method for the quantum collision operator.

研究动机与目标

• 解决在渐近保持格式中求解量子麦克斯韦分布（玻色-爱因斯坦或费米-狄拉克分布）时非线性方程的计算困难问题。
• 开发一种在所有Knudsen数范围内均稳定且一致收敛的量子玻尔兹曼方程数值格式。
• 用经典BGK算子替代量子BGK算子，避免求解复杂矩系统的迭代过程。
• 确保即使在平均自由程未被解析时，该格式仍能捕捉流体动力学极限（Euler方程）。
• 将谱方法扩展至量子碰撞算子，以实现高精度模拟。

提出的方法

• 以经典麦克斯韦分布 $ M_c $（其中温度由内部能量替代）作为惩罚项，对量子碰撞算子 $ Q_q $ 进行惩罚，即 $ M_c - f $。
• 利用经典麦克斯韦分布的前五阶矩（密度、动量、能量）与量子麦克斯韦分布一致，确保流体极限下的相容性。
• 采用半拉格朗日或有限体积格式对惩罚后的方程进行半隐式离散化，对BGK项采用隐式处理以保证稳定性。
• 通过守恒性质显式求解隐式BGK项，避免在每个时间步求解非线性系统。
• 采用球谐函数与高斯-勒让德求积法实现量子碰撞算子的谱方法，其中 $ M=4 $, $ N=32 $, $ L=8 $（玻色气体取 $ L=6 $）。
• 通过数值积分计算费米-狄拉克函数，通过级数展开计算玻色-爱因斯坦函数，以获得平衡分布。

实验结果

研究问题

• RQ1 在流体区域中，经典BGK算子是否能有效近似量子碰撞算子，而无需求解非线性矩系统？
• RQ2 所提出的格式是否在所有Knudsen数范围内均保持一致的稳定性和精度？
• RQ3 该格式是否能避免求解量子麦克斯韦分布所需的逸度-温度系统的昂贵迭代过程？
• RQ4 当平均自由程未被解析时，该格式在捕捉流体动力学极限（Euler方程）方面表现如何？
• RQ5 谱方法在模拟量子气体时，对量子碰撞算子的性能如何？

主要发现

• 新格式成功避免了在每个时间步求解非线性系统 $ \rho = \rho(z,T) $, $ e = e(z,T) $（其中 $ z $ 为逸度，$ T $ 为温度），消除了主要的计算瓶颈。
• 数值结果与量子Euler方程的KFVS格式在流体区域（$ \varepsilon=1\text{e-}4 $）下高度一致，验证了渐近保持行为。
• 在动能区域（$ \varepsilon=1\text{e-}1 $ 或 $ 1\text{e-}2 $），该格式与显式前向欧拉格式结果一致，验证了其在不同尺度下的精度。
• 即使采用大时间步长，该格式仍保持稳定且准确，表现出在Knudsen数上的均匀稳定性。
• 对量子碰撞算子的谱方法在 $ M=4 $, $ N=32 $, $ L=8 $（玻色气体取 $ L=6 $）下实现了高精度，能够精确计算量子碰撞项。
• 对玻色气体与费米气体的模拟结果在不同逸度值（$ z=0.01 $ 至 $ z=9 $）下表现出一致行为，验证了方法的鲁棒性。