一个时间分数阶扩散方程的源强识别反问题.docVIP

一个时间分数阶扩散方程的源强识别反问题 　　摘 要： 对于一类带有多个点源的分数阶维扩散方程问题，应用有限差分法给出了一个数值求解格式，并在已知点源个数及其位置的前提下，应用最佳摄动量正则化算法对源强进行了数值反演，并讨论了扰动数据和微分阶数的不同选取对反演算法的影响. 　　关键词： 分数阶扩散方程 多点源 最佳摄动量正则化算法 　　整数阶扩散方程是根据质量守恒定律和费克梯度扩散定律得出的，但是实际情况下许多扩散现象并不满足经典的费克梯度扩散定律，通常称之为“反常”扩散.分数阶导数与反常扩散过程本质上有一定的相似性，所以分数阶扩散方程应运而生，其研究已受到人们越来越多的关注.对于分数阶微分方程的求解，由于解析解一般都含有特殊函数或者复杂级数的形式，计算起来比较困难，因此越来越多的研究者讨论分数阶微分方程的数值解，而对于分数阶扩散方程的源项反演问题却少有研究.本文研究多点源分数阶对流-扩散方程的差分求解方法，进而应用最佳摄动量正则化算法对分数阶扩散方程中的多点源源强进行了数值反演. 　　考虑一个带多个点源的时间分数阶对流扩散方程的源强反演问题： 　　■-D■+v■=f（x，t） 0xl，0tT（1） 　　其中u=u（x，t）为t时刻x处污染物的浓度分布，D为扩散系数，D≥0，v为对流系数，v≥0，f（x，t）=■Q■δ（x-x■），x■为污染源的坐标位置，Q■为污染物排放强度（单位时间排放量），M为点污染源的个数，δ为狄拉克函数.■是Caputo意义下的分数阶导数，定义为 　　？鄣■■u（x，t）=■？蘩■■■■，0α1■，α=1 　　式中Γ表示Γ函数，Γ（x）=e■t■，且0α≤1. 　　初边值条件给定为： 　　u（x，0）=g（x），u（0，t）=u（l，t）=0（2） 　　假设知道污染源的个数，以及每个源的位置坐标，而需要确定各个污染源的强度值.这时，需要额外的关于污染物浓度分布的附加数据，联合模型（1）—（2）形成一个源强度识别的反问题.设在出流端t=T处观测到多处的浓度值为： 　　u（x■，T）=■■，k=1，2，…，K（3） 　　这样，由附加数据（3）式联合模型（1）—（2）式即构成一维多点源时间分数阶扩散模型的源强识别反问题. 　　1.正问题的求解 　　下面应用隐式差分方法给出问题的数值解.设x■=ih，h0；i=0，1，2，…，M；t■=nτ，τ0，n=0，1，2，…，N，其中h和τ分别是空间和时间步长，且Mh=l，Nτ=T.在方程（1）中采用如下有限差分近似： 　　？鄣■■u（x■，t■）=■■■？蘩■■■+O（τ）（4） 　　=■{u（x■，t■）-u（x■，t■）+■[u（x■，t■）-u（x■，t■）]}[（k+1）■-k■]+O（τ） 　　■|■=■+O（h■）（5） 　　■|■=■+O（h）（6） 　　将式（4）—（6）代入方程（1），可得 　　■{u（x■，t■）-u（x■，t■）+■[u（x■，t■）-u（x■，t■）]}[（k+1）■-k■]-D■+v■=f■（7） 　　其中i=1，2，…，M-1，f■=■Q■δ（ih-x■）.设p=■，q=■，u■■是u（x■，t■）的数值近似，即得 　　-pu■■+（1+2p-q）u■■-（q-p）u■■=u■■-■（u■■-u■■）[（k+1）■-k■]+τ■Γ（2-α）f■（8） 　　u■■=0，u■■=0，u■■=g（x■） 　　因此，可得隐式差分格式： 　　当n=0时， 　　-（q+p）u■■+（1+2p+q）u■■-pu■■=u■■+τ■Γ（2-α）f■（9） 　　当n0时， 　　-（q+p）u■■+（1+2p+q）u■■-pu■■=（2-2■）u■■-■u■■[2（k+1）■-（k+2）■-k■]+u■■[（n+1）■-n■]+τ■Γ（2-α）f■（10） 　　2.最佳摄动量正则化算法 　　最佳摄动量正则化算法是一种基于Tikhonov正则化、求解参数反演问题的迭代方法.对于源强反问题（1）—（2），设污染源的位置已知，即x■，x■，…，x■为，需要反演相应的源强度Q■，Q■，…，Q■.记R=[Q■，Q■，…，Q■.]，在R已知的前提下，借助上述求解正问题的方法，求出正问题（1）—（2）的解u（x，t；R），进而得到其在时刻t=T和x■处的值，记为u（x■，t；R）.这样求解反演问题（1）—（2）的一种优化方法就是使得计算值与观测值在某种误差意义下最小，通常化为下述极小问题的求解 　　min||u（x■，T；R）-■■||■+μ||R||■■，（11） 　　其中μ为正则参数.根据摄动算法，极小问题的求解又转化为对于给定的R■，求解最佳摄动量δR■，进而由下式确定出R■的一种迭代算法