广义时间分数阶扩散方程解的估计.pdfVIP

第 50卷 第 6期 复旦 学 报 (自然科学版) Vo1．50No．6 2011年 l2月 JournalofFudanUniversity(NaturalScience) Dec．2011 文章编号：0427—7104(2011)06—0705-09 广义时间分数阶扩散方程解 的估计 孙立博，徐 翔，张 瀛 (复旦大学 数学科学学院，上海 200433) 摘 要：采用特征函数展开的方法对一类广义时间分数阶扩散方程进行研究，导出解的L。。范数的估计，也可 以看作解的“强”极值原理．Luchko指出在适当的条件下，方程的解的极值可以在区域的底部和侧面取到．作为 与Luchko已得的极值原理对比，新获得的结果排除了区域内部极值点的存在，因此更具代表性和精确性．最后 得到了在连续和离散情况下的理论结果，并通过数值实验进行了验证． 关键词：分数阶扩散方程；离散极值原理；极值原理 中图分类号：0175．2 文献标志码：A 众所周知，用扩散方程来描述长时间和非局部依赖的不规则运动的规律是不恰当的．例如 Fokker (1914)、Planck(1917)和Kolmogorov(1938)得到的著名的FPK方程，它只能刻画局部的马尔可夫运动， 或者确切地说布朗运动．因为在布朗运动中，粒子下一时刻的状态只依赖于当前时刻状态，而与从初始时 刻起的历史状态无关．然而事实上，在 自然界中经常产生长记忆和非局部依赖的反常扩散，并且和某些随 机过程，例如分数布朗运动、L6vyflight、连续时间随机游走等关系紧密．半全局依赖的分数阶扩散方程 已经被证明能够很好的描述上述复杂的过程． 经过将近二十年的发展，带有分数阶的微分和积分理论现在在扩散方程领域扮演着非常重要的角 色．目前已经有很多强调分数阶扩散方程应用的理论结果，例如文献[1—9]．最近文献El0-1证明了广义时间 分数阶扩散方程的极值原理，由此很容易推出古典解的惟一性和对初值条件的连续依赖性．然而这个结 果没有排除极值可以同时在区域内部和边界上取到的情况，因此具有其 自身的局限性，例如在证明某些 反问题解的惟一性方面．另一方面有很多作者在研究分数阶扩散方程的数值实现，例如文献[11-13]．因此 离散的分数阶扩散方程的性质引起了很多数学家的兴趣． 本文我们主要推导一个 “强”极值条件，这个条件排除了区域内部的点存在极值的可能性．首先我们在理 论上获得了一个离散的极值条件，然后进行了数值实验．我们的工作是基于研究如下的广义时间分数阶扩 散方程 ： D ( ，￡)=== (z，) ( ，)∈ (0，L)×(0，T)， (1) Bu(x，)：==0 Xon8I， (2) u(x，O)一 ()． (‘3) 这里，a∈(0，1)是分数阶，B是边界算子， ·(n(z) )--k(x)是一个广义的二阶椭圆算子，口()0， k()≥O．Caputo意义下的分数阶的微分算子定义如下： ’ rt r，， 、 厂一而1 J。 d· (4) 这里 I、(·)为伽马函数． 1 连续情形下解的估计 引理 1 1)如果厂∈C1[O，T]在to∈(0， 达到极大值(或者极小值)，则D~f(to)≥O(或者≤0)．2)等号 收稿 日期：2009—11—11 作者简介：孙立博 (1983一)，男 ，硕士研究生，E-mail：sunlibo_comp_fdu@hotmail．CO1TI． 706 复 旦 学 报 (自然科学版) 第 5O卷 成立当且仅当f(t