[论文解读] Half-space Kinetic Equations with General Boundary Conditions

Qin Li, Jianfeng Lu+1 more|arXiv (Cornell University)|Sep 10, 2015

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

一句话总结

本文提出了一种针对具有通用边界条件（包括入射数据和多种反射类型，如漫反射、镜面反射、麦克斯韦反射）的线性半空间动力方程的数值方法。通过引入一种加阻尼-去阻尼的处理程序，建立了该方法的适定性，并证明了其逼近的次最优性。该方法在多组分、多频率及线性化BGK模型上得到了验证，精度极高。

ABSTRACT

We study half-space linear kinetic equations with general boundary conditions that consist of both given incoming data and various type of reflections, extending our previous work [LLS14] on half-space equations with incoming boundary conditions. As in [LLS14], the main technique is a damping adding-removing procedure. We establish the well-posedness of linear (or linearized) half-space equations with general boundary conditions and quasi-optimality of the numerical scheme. The numerical method is validated by examples including a two-species transport equation, a multi-frequency transport equation, and the linearized BGK equation in 2D velocity space.

研究动机与目标

• 将半空间动力方程的分析与数值解法扩展至包含入射数据与多种反射类型的通用边界条件。
• 在包含麦克斯韦反射与Cercignani-Lampis型反射的通用边界条件下，建立线性化半空间方程的适定性。
• 基于加阻尼-去阻尼程序，开发一种稳健且次最优的数值格式，用于求解此类方程。
• 在具有物理相关性的模型上验证该方法，包括双组分输运、多频率输运以及二维线性化BGK方程。

提出的方法

• 提出一种加阻尼-去阻尼程序，以处理具有复杂边界条件的半空间动力方程的适定性与数值求解问题。
• 应用一种加入耗散项的阻尼型动力方程，以稳定求解过程并支持数值计算。
• 构造一个特殊解 $ g_0 $，使其满足齐次方程与边界条件，从而使得精确解可表示为 $ f_h = f - c_h g_0 + c_h X_0 $，其中 $ f $ 为阻尼系统之解。
• 使用 $ L^2 $-投影算子 $ ilde{ ho} $ 与 $ ilde{ ho}^ot $ 投影到 $ ilde{ ho} $ 的零空间上，并定义 $ ilde{ ho}_1 $，以将零空间分解为正、负与零特征值子空间。
• 在频率变量 $ ilde{ ho} $ 上应用梯形积分规则，将问题降维为适合基函数构造的多组分系统。
• 处理边界算子 $ ilde{ ho} $ 与 $ ilde{ ho}_d $ 以模拟漫反射与镜面反射，并推导出 $ (I + ar{ ilde{ ho}})^{-1} $ 与 $ (I + ar{ ilde{ ho}})^{-1}(I - ar{ ilde{ ho}}) $ 的可逆形式。

实验结果

研究问题

• RQ1 如何为包含入射数据与多种反射机制的线性半空间动力方程建立适定性？
• RQ2 零空间分解 $ H^+ \bigoplus H^- \bigoplus H^0 $ 在确定无穷远处正确边界条件中的作用是什么？
• RQ3 加阻尼-去阻尼程序能否系统性地应用于构造此类方程的次最优数值格式？
• RQ4 该方法在具有多组分、多频率或多维速度结构的模型上的表现如何？
• RQ5 当精确解已知时（如 $ f_h = X_0 $ 的情形），该方法在多大程度上能恢复精确解？

主要发现

• 该方法在数值格式中实现了次最优性，确保了收敛速率在对数因子范围内最优。
• 在双组分输运方程中，当 $ h = X_0 $ 时，数值解与解析解完全一致，$ \tilde{\rho} $ 方向的误差低于 $ 10^{-10} $。
• 在麦克斯韦边界条件下，当 $ (\tilde{\rho}_d, \tilde{\rho}_s) = (0.3, 0.4) $ 时，方法成功恢复了精确解 $ f_h = X_0 $，且差值 $ f_h - X_0 $ 在机器精度范围内数值为零。
• 该方法成功处理了二维速度空间中的线性化BGK方程，在不同物理参数下均表现出强鲁棒性。
• 通过构造特殊解 $ g_0 $ 及使用系数 $ c_h = \frac{\braket{\tilde{\rho} X_0, f}}{\braket{\tilde{\rho} X_0, g_0}} $，实现了从阻尼系统中精确恢复解。
• 边界算子 $ \tilde{\rho} $ 被证明满足所需的强制性条件 $ (\tilde{\rho} f, f) \neq 0 $，从而保证了 $ (I + \bar{\tilde{\rho}}) $ 的可逆性。