带阻尼的多项Caputo-Fabrizio柯西问题的求解算子族

《Symmetry》：Resolvent Families for Multi-Term Caputo–Fabrizio Cauchy Problems with Damping

　　

研究人员构建了一个用于Banach空间中带阻尼的多项分数阶柯西问题的求解算子框架，其中分数阶导数为Caputo-Fabrizio（CF）型。利用拉普拉斯变换技术，研究人员获得了相关求解算子族生成的充要条件，并建立了广义解的存在唯一性。由于CF导数在原点为零，当????≠0时经典初始条件 u(0⁺)=x 无法达到；此处构造的求解算子产生了一个广义解，满足 u(0⁺)=T(0)x≠x。研究人员精确刻画了这一初始跳跃，证明其以 O(1?α) 的速率随 α→1^? 而消失，并指出其在物理上有意义的情形。该抽象理论被应用于具有结构阻尼的双调和板模型，产生了显式解析解和跳跃的定量界。

论文解读

分数阶微积分将整数阶微积分推广到任意阶，能够刻画具有记忆和遗传行为的复杂过程，已在反常扩散、粘弹性力学、控制理论和数学物理等领域得到广泛应用。抽象分数阶柯西问题将标量分数阶微分方程提升到Banach空间中的算子理论框架，Da Prato和Iannelli开创了求解算子族（resolvent family）方法，后经Lizama、Bazhlekova、Peng-Li等人发展完善。然而，基于Riemann-Liouville（RL）导数的框架虽然数学丰富，但其奇异核导致非经典的初始条件，并且物理意义不清晰。为克服此缺陷，Caputo和Fabrizio引入了基于非奇异指数核的CF导数，该导数不含原点奇异性，已在热传导和材料建模等领域得到应用。不过，Douaifia等人近期证明，采用CF导数等非奇异核的自治分数阶系统与经典初始条件在经典解意义下本质上不相容。具体而言，由于CF导数在原点为零，当方程右端A??≠0时，经典初始条件u(0⁺)=x无法满足，初始时刻必然发生跳跃。

针对上述问题，研究人员在Banach空间X中研究了带阻尼的多项CF抽象分数阶柯西问题（CFAFCP）：{多项CF导数组合} u(t)=Au(t)，t>0，其中α∈(0,1)，β≥0为阻尼参数，c为阻尼系数，A为闭稠定线性算子，x为初始状态。该问题中，CF核的非奇异结构使得经典解不存在，但可通过分布意义下的拉普拉斯变换将x理解为左极限，从而构造广义解。本文的核心贡献是构建了CF型α-求解算子族框架，证明其生成定理，并建立广义解的存在唯一性；该框架揭示出初始跳跃是CF内核的固有量化结果，并给出其随α→1^?消失的速率。

该研究的意义在于首次为含阻尼的CF型抽象柯西问题提供了统一的算子理论处理。与RL框架中求解算子在原点具有t^?α奇异性不同，CF求解算子在t>0上有界且收敛到非平凡初始算子T(0)，从而精确刻画了初始跳跃。该框架还澄清了先前文献中关于CF导数与经典初值不相容性的争议，表明该不相容性可通过广义解和初始跳跃得到定量解释。论文发表在《Symmetry》。

在技术方法上，研究人员主要运用拉普拉斯变换技术处理CF核的测度卷积结构。首先证明了与CF导数相关的卷积核k不是局部可积函数，而是局部有限测度，其中包含位于原点的Dirac原子，其拉普拉斯变换为有理函数，具有两个负实轴极点。基于此，定义了CF型α-求解算子族，并利用拉普拉斯变换刻画生成元的谱条件，得到充要的生成定理。解的性质方面，采用测度卷积意义下的Volterra积分方程和拉普拉斯变换方法，证明广义解的存在唯一性，并通过解析表示分离初始跳跃项。最后，以正方形域上的双调和板模型作为具体应用，通过特征函数展开和谐振条件验证求解算子族条件，获得显式解和跳跃的定量界。

研究结果分为三个部分：

（一）CF-求解算子族的基本性质（第3节）。通过定义CF型α-求解算子族{T(t)}_t≥0，研究人员证明了其生成元A的谱条件：σ(A)包含于截断区域之外，这源于核函数K的极点位置。生成定理（定理1）给出了A生成指数有界CF-求解算子族的充要条件：存在强连续且有界的一致有界算子族满足拉普拉斯变换恒等式。还证明了该求解算子族是唯一的，并推导出初始算子满足代数关系T(0)=I?(γ/α)A，其中γ为系统参数。特别地，当A??≠0时，初始跳跃T(0)x?x=?(γ/α)A??非零，体现了CF核的非平凡贡献。

（二）广义解的存在唯一性与初始跳跃（第4节）。通过定义满足分布意义下CF导数的广义解，研究人员建立了表示公式：u(t)=T(t)x=T(0)x+∫₀^tk₀(t?s)T(s)xds，其中k₀为核的光滑部分。该式明确显示解从T(0)x出发而非x，并通过光滑Volterra扰动演化。进一步证明了，对于任意初始数据x，u(t)=T(t)x是广义解且唯一，唯一性建立在拉普拉斯变换的唯一性和解的连续性上。定理7进一步证明，当α→1^?时，初始跳跃以O(1?α)速率消失，说明CF框架与经典一阶初值问题渐近相容。

（三）双调和板模型的应用（第5节）。在正方形域上考虑带结构阻尼的分数阶双调和板方程，其中A为双调和算子减去基础系数λ与恒等算子之和，施加Navier边界条件。特征函数为双正弦函数，特征值严格递减且为负。通过验证谱条件、算子交换性、拉普拉斯变换及初始值定理，证明了A生成指数有界CF-求解算子族。利用Routh-Hurwitz判据和紧致性论证，获得了谱参数θ上统一的指数衰减界，进而得到整体一致有界性。命题4给出了初始跳跃的定量界限，表明跳跃随α接近1、初始剖面光滑或其双调和能量减小而变小。数值实验以第一模态为例展示了跳跃随α、阻尼系数c和基础系数λ的变化，验证了跳跃以O(1?α)消失的速率。

讨论与结论：研究人员从算子理论视角总结了CF框架与RL框架的根本差异：RL求解算子在原点具有奇异性，初始条件以卷积意义施加；CF求解算子在原点有界且趋于非平凡算子T(0)，该有界性直接源于非奇异指数核。初始跳跃T(0)x?x是CF模型的结构特征，而非缺陷，它量化了非奇异核导致的初始数据非经典效应。在物理上，CF模型适用于具有瞬时响应的系统，如标准线性固体模型中的瞬时弹性应变；对于不存在瞬时响应的系统，非奇异核可能是不合适的建模选择。未来方向包括高阶正则性分析、非齐次项与时间相关算子的推广、与其他非奇异核（如Atangana-Baleanu）的系统比较，以及捕捉初始跳跃的数值方法开发。