[论文解读] On linear hypersingular Boltzmann transport equation and its variational formulation

Jouko Tervo|arXiv (Cornell University)|Aug 29, 2018

Numerical methods in engineering参考文献 18被引用 3

一句话总结

本文为带电粒子输运中出现的线性超奇异玻尔兹曼输运方程（BTE）提出了一种严格的变分公式，特别针对Møller型散射，其微分截面表现出超奇异特性。结果表明，精确的BTE包含与Hadamard有限部分积分算子耦合的一阶能量导数以及二阶角向导数，通过变分方法推导出一种CSDA-Fokker-Planck型近似算子，使得利用Lions-Lax-Milgram理论进行解的存在性证明成为可能。

ABSTRACT

For charged particle transport the linear Boltzmann transport equation (BTE) turns out to be a partial hyper-singular integro-differential operator. This is due to the fact that the related differential cross-sections $σ(x,ω',ω,E',E)$ may have hyper-singularities. In these cases the energy integral appearing in the collision terms must be interpreted as the Hadamard finite part integral leading to hyper-singular integral operators. The article considers a refined expression for the exact transport operator and related variational formulations of the inflow initial boundary value problem for one particle equation containing hyper-singularities. We find that the exact BTE contains the first-order partial derivatives with respect to energy combined by partial Hadamard (first-order) singular integral operators. In addition, it contains the second-order partial derivatives with respect to angle and some mixed terms. The analysis will be carried out only for the so called Møller-type interaction (scattering) which is a kind of prototype of hyper-singular interactions. The generalizations to other type of collisions, such as to Bremsstrahlung, go analogously. We also expose a weak form (the variational formulation) of the hyper-singular transport problem. Another variant variational formulation decreases the level of singularities in the integration (appearing in the due bilinear form) containing only singularities of order one that is, singularities like ${1\over{E'-E}}dE' dE$. The variational formulation is an essential step in order to show the existence of generalized solutions e.g. by Lions-Lax-Milgram Theorem based methods (proceedings for solution spaces and existence theory are omitted here). The corresponding approximative transport operator is deduced. It turns out to be a CSDA-Fokker-Planck type operator.

研究动机与目标

• 分析当微分截面包含超奇异特性时，线性玻尔兹曼输运算子的精确结构，特别是Møller散射情形。
• 为超奇异BTE的流入初始-边界值问题构建一个无近似假设的变分公式。
• 降低双线性形式中的奇异性阶数，以实现通过Galerkin有限元方法进行数值处理。
• 利用变分方法和Lions-Lax-Milgram定理，为广义解的存在性建立理论基础。
• 推导一种与CSDA-Fokker-Planck模型结构相符的近似输运算子，适用于数值计算。

提出的方法

• 本文指出，微分截面中的超奇异特性导致与能量相关的偏导数（一阶）与Hadamard有限部分积分算子耦合。
• 将精确输运算子表述为包含二阶角向导数和混合项的偏积分微分方程，其中碰撞积分被解释为Hadamard有限部分积分。
• 通过在适当的Sobolev-Slobodevskij空间中使用测试函数，推导出弱形式（变分公式），将双线性形式中的奇异性阶数降低至一阶。
• 证明该变分公式可通过Lions-Lax-Milgram定理应用于存在性理论，从而实现广义解的存在性证明。
• 通过主粒子的泰勒展开和角向近似，推导出近似输运算子，其形式为CSDA-Fokker-Planck型。
• 分析局限于Møller型散射，但可通过类似方法推广至其他相互作用（如轫致辐射）。

实验结果

研究问题

• RQ1 微分截面中存在超奇异特性时，如何改变线性玻尔兹曼输运方程的结构？
• RQ2 当碰撞积分被解释为Hadamard有限部分积分时，精确输运算子的精确形式是什么？
• RQ3 能否为超奇异BTE构建一个变分公式，以降低奇异性阶数并支持存在性证明？
• RQ4 从精确超奇异算子推导出的近似输运算子的结构是什么？
• RQ5 该变分公式如何支持Galerkin有限元方法等数值求解策略？

主要发现

• Møller散射的精确线性BTE包含与Hadamard有限部分积分算子耦合的一阶能量偏导数，形成类似伪微分算子的项。
• 输运算子还包含二阶角向偏导数和混合导数项，其角向部分类似于Fokker-Planck算子。
• 变分公式将双线性形式中的奇异性阶数降低至一阶，具体形式为 $ \frac{1}{E' - E} dE' dE $，从而实现数值稳定性与收敛性分析。
• 推导出的近似输运算子为CSDA-Fokker-Planck型，与放疗剂量计算中常用的常规近似一致。
• 变分公式为通过Lions-Lax-Milgram定理严格证明广义解的存在性提供了数学路径，即使对于超奇异算子亦成立。
• 分析表明，流入边界值问题具有可变边界重数，增加了适定性分析的复杂性，但变分框架为此类情况提供了可行的处理路径。