[论文解读] Existence and Regularity of solutions to 1-D Fractional Order Diffusion Equations

Lueling Jia, Huanzhen Chen+1 more|arXiv (Cornell University)|Aug 30, 2018

Fractional Differential Equations Solutions被引用 11

一句话总结

本文利用黎曼-刘维尔算子与黎曼-刘维尔-卡波托算子，建立了1D稳态分数阶扩散方程解的存在性与正则性。研究结果表明，解的正则性关键取决于右端函数与边界条件；黎曼-刘维尔算子的核具有代数奇异性，而黎曼-刘维尔-卡波托算子的核则为代数正则，从而明确了在狄利克雷与诺伊曼边界条件下模型的适定性。

ABSTRACT

In this article we investigate the existence and regularity of 1-d steady state fractional order diffusion equations. Two models are investigated:the Riemann-Liouville fractional diffusion equation, and the Riemann-Liouville-Caputo fractional diffusion equation. For these models we explicitly show how the regularity of the solution depends upon the right hand side (rhs) function. We also establish for which Dirichlet and Neumann boundary conditions the models are well posed.

研究动机与目标

• 研究由黎曼-刘维尔与黎曼-刘维尔-卡波托算子控制的1D稳态分数阶扩散方程解的存在性与正则性。
• 确定解的正则性如何依赖于右端函数与边界条件的选择。
• 建立模型适定性的条件，特别是针对狄利克雷与诺伊曼边界条件。
• 为真实世界应用中选择合适的分数阶扩散模型提供物理与数学上的洞见。

提出的方法

• 分析采用基于勒让德多项式展开的谱方法，研究加权L2空间中解的结构。
• 黎曼-刘维尔-卡波托算子定义为 $\,_{RLC}\mathcal{D}_{r}^{\alpha}u(x) := -D\left(rD^{-(2-\alpha)} + (1-r)D^{-(2-\alpha)*}\right)Du(x) = f(x)$，其中 $D$ 为微分算子，$D^{-\beta}, D^{-\beta*}$ 分别为左、右分数阶积分。
• 黎曼-刘维尔算子定义为 $\,_{RL}\mathcal{D}_{r}^{\alpha}u(x) := -D^{2}\left(rD^{-(2-\alpha)} + (1-r)D^{-(2-\alpha}*\right)u(x) = f(x)$，其与前者的主要区别在于微分算子的位置。
• 论文使用加权Sobolev空间 $H^{\alpha/2}_{0}(0,1)$ 与 $H^{-\alpha/2}(0,1)$ 来定义弱解并分析正则性。
• 关键估计通过雅可比多项式及其范数的性质推导，包括利用斯特林公式获得的渐近行为。
• 通过关于系数衰减与级数收敛性的引理，确保了解在适当函数空间下各项的有界性。

实验结果

研究问题

• RQ1 对于黎曼-刘维尔-卡波托模型，1D分数阶扩散方程解的正则性如何依赖于右端函数？
• RQ2 黎曼-刘维尔-卡波托算子的核结构是怎样的？其如何影响解的正则性？
• RQ3 在何种狄利克雷与诺伊曼边界条件下，黎曼-刘维尔与黎曼-刘维尔-卡波托模型是适定的？
• RQ4 两种分数阶扩散模型在能量通量的物理解释上存在哪些差异？
• RQ5 在何种源项条件下，可保证解在加权Sobolev空间中存在且具有正则性？

主要发现

• 黎曼-刘维尔-卡波托算子的核是代数正则的，意味着解的正则性可由右端函数继承。
• 黎曼-刘维尔算子的核具有代数奇异性，即使右端函数光滑，解在边界附近仍表现出奇异行为。
• 当 $f \in H^{-\alpha/2}(0,1)$ 时，在狄利克雷条件下，两种模型在 $H^{\alpha/2}_{0}(0,1)$ 中均存在弱解。
• 黎曼-刘维尔模型在诺伊曼条件下可能不适定，具体取决于边界条件类型，此结论已有前期研究支持。
• 谱方法确保了解系数的收敛性与有界性，衰减速率通过斯特林近似推导得出。
• 在 $f \in L^2_{\rho^{(\beta,\alpha-\beta)}}(I)$ 条件下，级数 $\sum_{i=0}^{\infty} \mu_i c_i G_{i+1}^{(\beta-1,\alpha-\beta-1)}(0)$ 收敛，从而保证边界条件的一致性。