时间分布阶扩散方程的快速算法研究

来源 :山东大学 | 被引量 : 0次 | 上传用户:hulin99
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
近年来,关于时间分布阶扩散方程的研究得到了广泛的关注,许多扩散指数随时间变化的复杂扩散过程,如减速的慢扩散和加速的超扩散、减速的超扩散和加速的慢扩散等,都可以用这类方程来描述。目前时间分布阶扩散方程已经被广泛应用于复合材料的流变特性、信号控制和处理、介电感应及扩散以及粘弹性材料的应力-应变行为等各项工程技术领域。  本文研究的主要内容是:时间分布阶扩散方程的两种有限差分格式的快速计算方法以及拟线性时间分布阶扩散方程的一种有限差分格式,模型如下:{Dwtu(x,t)=(e)2u/(e)x2+f(x,t),0<x<L,0<t≤T,u(x,0)=0,0<x<L,u(0,t)=ψ1(t),u(L,t)=ψ2(t),0≤t≤T,其中,ψ1(0)=0,ψ2(0)=0,且Dwtu(x,t)=∫10w(α)C0Dαtu(x,t)dα,C0Dαtu(x,t)={1/Γ(1-α)∫t0(t-ξ)-α(e)u/(e)ξ(x,ξ)dξ,0≤α<1,ut(x,t),α=1.w(α)>0,f10w(α)dα=c0>0,f,ψ1和ψ2均为已知函数。  因为时间分布阶导数的非局部性质,利用有限差分法求解该方程在当前时间层的数值解时需要用到之前所有时间层的函数值,这导致了右端项的计算量非常巨大,当求解结束时,计算量能够达到O(MN2J),这里M为空间剖分个数,N为时间剖分个数,J是对分布阶导数剖分的个数。当网格剖分加密时,计算量和存储量会越来越大,计算时间会越来越长。通过研究离散得到的系数矩阵的结构发现,将系数矩阵通过一种合理的方式重新转化后,可以得到由块状Toeplitz矩阵构成的矩阵,结合快速Fourier变换及逆变换,我们可以得到该模型的快速求解方法,将计算量减少到O(MNlogN)+O(NJ),从而大大降低了计算时间。在后面,我们将考虑右端项f为非线性项的情况,并设f局部有界且满足局部Lipschitz条件,提出了一种有限差分格式,证明了该格式的收敛阶为O(τ+h2+△α2)。  本文主要针对上述模型讨论了快速算法,并分析了右端项为非线性的情况。全文共分为四章:  第一章:介绍了时间分步阶扩散方程的背景、国内外研究现状以及模型知识。  第二章:结合时间分布阶扩散方程的两种有限差分格式,推导了这两种格式的快速算法,并通过数值算例验证了快速算法的有效性。  第三章:给出了右端项为非线性项的有限差分格式,并证明了该格式的收敛性,并用数值算例验证了理论结果。  第四章:对本文内容做了总结和展望。
其他文献
测验等值是教育学、心理学中的一项重要研究内容,它对于考试的公平性、可比性、题库建设、教学质量评价和计算机自适应测验都具有重要的意义。但是对于多维结构的测验,如果仍使
在各类工程,金融,数学等领域中,经常要遇到类似下面形式的方程dyt=f(yt)dxt,(0.0.1)其中x是一个多维的驱动信号,f是一列驱动的向量场。1如果x∈C1或者x∈C1-var那么这个方程可以理
网格曲面已经被广泛应用于计算机图形学和几何造型中。随着网格数量的快速增长和质量的不断提高,实际中产生了很多基于人类感知的网格应用,这预示着在处理网格的时候要将人类
近年来,复Banach空间几何理论的研究已经逐渐成为国内外数学工作者所关注的领域。复Banach空间几何性质的讨论起源于向量值解析函数相关性质方面的研究,在研究过程中学者们发现
本文考虑如下带五次项的更一般的非线性Schr(o)dinger方程的初边值问题{iut+uxx+q(|u|2)u-β|u|4u=f(x,t)u, xl≤x≤xr,0<t≤Tu(xl,t)=u(xr,t)=0u(x,0)=u0(x)其中q(s),f(x,t)为已知的
外汇期权作为一种出现最晚发展最快的金融衍生品,其良好的避险作用,受到各个国家和经济体的青睐,对于其定价模型的研究也在不断的发展与完善.传统的均值回归对数模型在波动率为
在现代物理学的研究中,出现了许多非线性发展方程.电报方程就是从研究电报线上电压和电流的变化规律推导出来的,它描述了均匀传输线上电压和电流的关系,所以又被称为传输线方程.
可靠性分析是统计学的一门重要分支,有着较强的应用性,因此受到了工程界等的高度重视。在实际的工作和生活中我们经常会遇到可靠性问题,例如,我们在购买电视机等家用电器时都想要
图染色理论和极值图论是图论中非常重要的两个研究课题,在计算机科学、信息安全以及模式匹配等领域有着广泛的应用。在本文中,我们主要讨论了以下几类图染色问题以及3-一致超图