2025-11-25T08:43:16.765271

Remarks on propagation of discontinuities in stationary radiative transfer

Kawagoe

We consider the stationary transport equation with the incoming boundary condition. We are interested in discontinuities of the solution. Under the generalized convexity condition, it is known that it has only boundary-induced discontinuities, which are discontinuities arising from discontinuous boundary data, they propagate along positive characteristic lines, and we can reconstruct the attenuation coefficient from boundary measurements by the inverse X-ray transform. In this article, we observe that coefficient-induced discontinuities, discontinuities of the solution arising from discontinuous coefficients, would also appear without the generalized convexity condition. If the set of discontinuous points of the coefficients contains at most finite number of flat parts, coefficient-induced discontinuities do not affect the inverse X-ray transform. We also remark that, under the generalized convexity condition, a three dimensional inverse problem can be reduced to the two dimensional one. A numerical experiment is exhibited.

academic

Remarks on propagation of discontinuities in stationary radiative transfer

基本信息

论文ID: 2510.12130
标题: Remarks on propagation of discontinuities in stationary radiative transfer
作者: Daisuke Kawagoe (京都大学信息学研究科)
分类: math.AP (数学分析)
发表时间: 2025年10月15日
论文链接: https://arxiv.org/abs/2510.12130

摘要

本文研究带有入射边界条件的稳态输运方程中解的不连续性传播问题。在广义凸性条件下，已知解仅存在边界诱导的不连续性，这些不连续性源于不连续的边界数据，沿正特征线传播，并可通过逆X射线变换从边界测量中重构衰减系数。本文观察到，在不满足广义凸性条件时，还会出现系数诱导的不连续性，即由不连续系数引起的解的不连续性。如果系数不连续点集合最多包含有限个平坦部分，系数诱导的不连续性不会影响逆X射线变换。文章还指出，在广义凸性条件下，三维逆问题可以化简为二维问题，并展示了数值实验。

研究背景与动机

问题背景

本文研究的核心问题是稳态辐射传输方程： $\xi \cdot \nabla_x I(x, \xi) + (\mu_a(x) + \mu_s(x))I(x, \xi) = \mu_s(x) \int_{S^{d-1}} p(x, \xi, \xi')I(x, \xi') d\sigma_{\xi'}$

这个方程描述了光子在浑浊介质（如生物组织）中的传播，其中：

$I(x, \xi)$ ：位置 $x$ 处沿方向 $\xi$ 的光子密度
$\mu_a$ ：吸收系数
$\mu_s$ ：散射系数
$p$ ：散射相函数
$\mu_t = \mu_a + \mu_s$ ：总衰减系数

研究动机

理论完善需求：现有理论主要在广义凸性条件下研究边界诱导的不连续性，但实际应用中这一条件可能不满足
实际应用驱动：光学层析成像等医学成像技术需要从边界测量重构内部参数
数值实验启发：文献4中的数值实验显示，即使违反广义凸性条件，重构仍可能成功
逆问题求解：理解不连续性传播对于逆X射线变换的稳定性至关重要

现有方法局限性

严格依赖广义凸性条件
仅考虑边界诱导的不连续性
对系数不连续性的影响缺乏系统分析
三维问题的计算复杂性高

核心贡献

扩展了不连续性理论：在不满足广义凸性条件时，系统分析了系数诱导不连续性的出现和传播
建立了新的分类框架：将解的不连续性分为边界诱导（boundary-induced）和系数诱导（coefficient-induced）两类
证明了逆变换的稳定性：当系数不连续点集合最多包含有限个平坦部分时，系数诱导不连续性不影响逆X射线变换
提出了维数约简方法：在广义凸性条件下，将三维逆问题化简为二维问题
提供了数值验证：通过具体实验验证了理论结果的有效性

方法详解

任务定义

研究稳态输运方程边值问题(1.1)-(1.2)解的不连续性：

输入：边界条件 $I_0(x,\xi)$ ，系数 $\mu_a, \mu_s, p$
输出：解 $I(x,\xi)$ 的不连续性集合
目标：从边界测量重构衰减系数 $\mu_t$

理论框架

1. 广义凸性条件回顾

对于区域划分 $\{Ω_j\}_{j=1}^N$ ，广义凸性条件要求：对所有 $(x,ξ) \in Ω \times S^{d-1}$ ，半直线 $\{x-tξ|t≥0\}$ 与 $∂Ω_0$ 的交点数有限。

2. 解的积分表示

通过沿特征线积分，得到解的积分方程： $I(x,ξ) = \exp(-M_t(x,ξ;τ_-(x,ξ)))I_0(P(x,ξ),ξ) + \int_0^{τ_-(x,ξ)} \mu_s(x-sξ)\exp(-M_t(x,ξ;s)) \times \int_{S^{d-1}} p(x-sξ,ξ,ξ')I(x-sξ,ξ')dσ_{ξ'}ds$

其中 $M_t(x,ξ;s) = \int_0^s \mu_t(x-rξ)dr$ 。

3. 不连续性分类

二维情况（定理3.1）： $\text{disc}(I) \subset \text{disc}_B \cup \text{disc}_C$

其中：

$\text{disc}_B$ ：边界诱导不连续性
$\text{disc}_C = \cup_{j=1}^K \{(x+tξ_±^{(j)}, ξ_±^{(j)}) | x \in L_j, 0 ≤ t < τ_+(x,ξ_±^{(j)})\}$

三维情况： $\text{disc}_{C,3d} = \cup_{j=1}^K \{(x+tξ,ξ) | x \in F_j, ξ \in S_j^2, 0 ≤ t < τ_+(x,ξ)\}$

技术创新点

迭代分解方法：将解分解为 $I = \sum_{n=0}^∞ I^{(n)}$ ，其中 $I^{(0)}$ 包含边界数据， $I^{(n≥1)}$ 表示散射贡献
不连续性传播分析：
- $I^{(0)}$ 的不连续性来源于边界数据和 $M_t$ 函数
- $I^{(n≥1)}$ 的不连续性仅来源于系数
跳跃量计算：对于边界诱导不连续性，跳跃量为： $[I](x^*,ξ^*) = C\exp\left(-\int_0^{τ_-(x^*,ξ^*)} \mu_t(x^*-rξ^*)dr\right)$

实验设置

三维数值实验配置

区域：单位球 $Ω = \{x \in \mathbb{R}^3 | |x| < 1\}$
区域划分：
- $Ω_1$ ：环形区域， $\mu_a = 0.3$
- $Ω_2$ ：球形区域， $\mu_a = 0.2$
- $Ω_3$ ：其余区域， $\mu_a = 0.1$
散射系数： $\mu_s = 0.3$ （常数）
相函数：Henyey-Greenberg核， $g = 0.9$