时间分数阶热方程透明边界条件是否保持扩散和规则？

和规则决定透明性而非通量守恒

图结点上的透明边界条件，意在让入射脉冲像在无分支直线上那样进入出射键，不产生虚假反射。在三键星图上，此前针对 β = 1 的经典结果给出了实现这一点的条件：扩散系数须满足 κ1 = κ2 + κ3，其中键 1 是入射键，键 2 和键 3 是出射键。

有两件事容易混淆，值得分开。结点处的基尔霍夫通量平衡与和规则无关；论文通过数值验证表明，无论和规则是否满足，q1 = q2 + q3 都精确到机器精度（最大偏差约 10⁻¹⁷–10⁻¹⁸）。和规则管的是另一件事：入射通量如何在两条出射键之间分配，以及键 1 自身是否看到反射分量。破坏和规则不会改变通量恒等式，却会在键 1 上产生可见反射。

结点核的阶数改变但条件不变

作者用拉普拉斯变换求解外部键问题，证明入射键所看到的非局部结点算子是阶为 β/2 的 Abel 型卷积核。当 β = 1 时，它退化为此前经典构造中出现的经典半导数核 (t − τ)^(−1/2)。

该核仍是弱奇异的，但可积：由于 0 < β/2 ≤ 1/2 < 1，(t − τ)^(−β/2) 在 [0, t] 上的积分对每个有限 t 都有限。因此该算子是作用在边界迹上的真正卷积，而非形式对象。对读者而言，关键在于分数阶同时通过同一个核进入分支侧和自由空间侧的算子，这正是它从透明性条件中消去的原因。

和规则为何在分数阶推广下依然成立

论文的主要解析结果是：κ1 = κ2 + κ3 对 (0,1] 内的每一个 β 都消除回流，与 β = 1 时完全一样。作者严格证明了这一点，并证明了耦合结点问题的唯一性。

其机制是消去而非巧合。由于分数阶在分支侧和无分支侧以同一个 Abel 型核出现，透明性条件最终只取决于扩散系数。论文直接指出：结点的透明性只取决于网络的扩散系数，而不取决于输运律的记忆效应。这是读者应当记住的结果，它比仅仅观察到该规则在若干测试过的 β 值下碰巧成立要更强。

数值验证是如何设置的

解析结果是在耦合偏微分方程组层面验证的，而不是通过离散化非局部结点算子本身。每条键用 Caputo 导数的标准隐式 L1 格式离散，并配合全隐式空间拉普拉斯算子。L1 权重 bj = (j+1)^(1−β) − j^(1−β) 在每个时间步都乘以整个节点历史，因此计算代价和内存随时间步数 Nt 按 O(Nt²) 增长，而非经典隐式格式的 O(Nt)。

结点条件只涉及同一时间层上的值，按经典 Crank-Nicolson 格式离散：基尔霍夫通量用单侧有限差分，并用一个共享的结点未知量强制连续性。分数阶只影响内部节点方程，不影响结点耦合代数。作者指出，若要模拟显著更长的时间，需要快速卷积技术。

验证将键 1 上的星图解 u1(x, t) 与参考解比较，后者是在一条无分支直线上用相同扩散系数 κ1 求解同一时间分数阶热方程得到的，直线被对称截断且截断足够远，在所考虑的时间范围内不可见。两者使用相同格式、相同 Δx 和 Δt，以及键 1 上相同的高斯初始数据。若结点透明，两个解应在 x ≤ 0 上重合；任何偏差都量化了回流。

透明性误差在和规则处始终最小

论文在固定末时刻 T 测量相对 L2 透明性误差 E(κ1)，在 x ≤ 0 上比较键 1 的星图解与无分支参考解。在 β = 0.99、0.7 和 0.5，且 κ2 = 0.089、κ3 = 0.08 固定的情况下，E 在 κ1 = κ2 + κ3 = 0.169 处或紧邻处取得最小值。

表 2 给出了 κ1 = 0.169 与 κ1 = 0.35 在五个 β 值下的数值。满足和规则时，E 在 β = 0.99 为 3.74%，β = 0.90 为 3.03%，β = 0.70 为 1.60%，β = 0.50 为 0.68%，β = 0.30 为 0.23%。破坏和规则时，对应值为 20.07%、17.65%、12.34%、7.55% 和 3.87%。和规则成立时，相对误差始终小三到五倍。

κ1 = 0.169 处的残余误差并非恰好为零，论文将其归因于 L1/隐式格式的离散化伪影，而非和规则失效。作者还指出，这一比较是在耦合偏微分方程组及其标准基尔霍夫条件层面验证和规则的效果；它并非在解析论证已确立的内容之外，对特定算子恒等式提供独立的数值确认。

反射体现在温度剖面上

对于 β = 0.7 且满足和规则的情形，入射脉冲在键 2 和键 3 之间平滑分裂，键 1 上除普通亚扩散展宽外没有可见畸变，星图解与无分支参考解在整个演化过程中在视觉上无法区分。

当和规则被破坏（κ1 = 0.35）时，键 1 在 t = 15、30 和 60 时于结点附近出现一个可见的次级隆起，这在无反射情形中不存在。这是脉冲的反射分量，直接体现在空间剖面上，而不仅仅体现在聚合诊断量 E(κ1) 中。论文指出，结点通量图的两个面板在形状上定性相似，因此破坏和规则引起的反射并不直接体现在结点通量平衡本身。

收敛率揭示结点在何处损失精度

在无耦合参考直线上，β = 0.7 且 Δx = 0.025 固定时细化 Δt，观测到的时间收敛率约为 1.03–1.05，低于渐近率 2 − β = 1.3。论文将此差距归因于 t = 0 附近有限的时间正则性：即使初始数据光滑，时间分数阶方程的解在该处通常也只保留有限的时间正则性，而这一初始层奇异性主导了均匀网格上的观测率。Δt = 0.02 固定时细化 Δx，观测率约为 2.0–2.2，与预期的二阶空间精度相符。

在满足和规则的耦合星图系统上重复该研究，远离结点处（x = −2）的收敛率与无耦合参考直线接近：时间约 1.03–1.05，空间约 2.04–2.22。在结点本身（x = 0），空间收敛率降至约 1.21–1.35，证实用于基尔霍夫通量的单侧有限差分是结点附近空间误差的主要局部来源，而每条键的主体仍保持完整的二阶精度。

为确认均匀网格时间收敛率降低确实源于 t = 0 初始层，作者在分级网格 tk = (k/N)^r T 上重复时间细化，对 β = 0.7 取 r = (2 − β)/β ≈ 1.86。分级网格给出的观测率约为 1.28–1.29，与理论率 2 − β = 1.3 接近，而在相同时间步数下均匀网格为 1.03–1.06。作者还报告，早期实现中存在一个真实的 L1 历史配对索引错误，仅靠和规则比较无法发现，因为它会一致地降低参考解和星图解的精度，并在两者之差中大部分抵消。

不变性结果的适用范围止于何处

和规则不变性是在三键度量星图上、所有键共享单一分数阶 β 的条件下建立的。论文指出该结果可推广到 n 键星图，但树、环和多结点图未予处理。键依赖的阶数 βj 会破坏不变性结果背后的共同核论证；作者将其与已有的边依赖阶数模型和时空分数阶度量图框架联系起来，但并未用它们解决该问题。

论文证明了耦合结点问题的唯一性，但未给出与之匹配的存在性理论，完整的适定性与唯一性处理留待未来工作。真正离散的分数阶结点条件也仍未解决。这些是解析结果的边界，而非数值演示的边界；数值演示被明确界定为在耦合偏微分方程组层面验证和规则的效果，而非对算子恒等式的独立确认。

本文确立的实际优先级是一个比较关系，而非新阈值：对于共享分数阶的三键星图，κ1 的透明性最优选择在每个测试过的 β 下都是 κ1 = κ2 + κ3，偏离它的代价是相对透明性误差增大三到五倍，并在入射键上出现可见的反射隆起。读者若要决定如何在分支亚扩散网络中分配扩散系数，应把和规则作为目标，并视分数阶与该目标无关，同时记住证据覆盖的是共享单一 β 的星图。

本文引用的文献