时间分数阶Allen-Cahn方程紧格式及其快速算法

来源 :山东大学 | 被引量 : 0次 | 上传用户:q363342684
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
近年来,关于分数阶相场模型的研究得到了广泛的关注,由于分数阶算子的非局部性,它可以更准确地描述相变机理,例如连续时间随机游走的粒子运动、粒子之间的长距离相互作用和结构损伤、疲劳的演化过程。本文研究的主要内容是:时间分数阶Allen-Cahn方程在非均匀时间步长下的紧致有限差分格式,分析它的收敛性和稳定性,并介绍了三类加速方法来提高计算性能。模型如下:(?)其中Ω=(0,L)2的闭包为(?),非线性体积力f(u)=u3-u,常数0<ε(?)1是描述材料之间过渡边界厚度的相互作用长度,?tα:=0CDtα表示α阶的Caputo导数:(?)时间分数阶Allen-Cahn方程的解具有初始奇点,而且它的数值解允许采用多个时间尺度,即初始动力学演变以快速的时间尺度发展,而随后的粗化阶段以非常慢的时间尺度发展。因此,采用非均匀时间步长(例如渐变网格)是合理且必要的。由于时间分数阶导数的历史依赖性,当我们考虑使用有限差分方法求解该方程在某一时间层的数值解时,需要用到之前所有时间层的函数值,致使计算工作量变得非常巨大。因为分数阶微分方程可以通过在时间或空间上用非局部关系的幂律记忆核得到,我们考虑在某一时间段上逼近核函数,即采用指数求和(SOE)算法求解时间分数阶Allen-Cahn方程。当我们求解某一时间层的数值解时,只需使用前面部分时间层的函数值,进而减少计算量。此外,还可以考虑采用自适应时间步长策略和减少自由度的方法来对方程的计算过程进行加速。本文主要研究了时间分数阶Allen-Cahn方程的紧致有限差分格式的收敛性和稳定性,讨论了关于该模型的加速方法。全文共分为五章:第一章:介绍了时间分数阶Allen-Cahn方程的背景、国内外研究现状以及问题模型。第二章:提出了非均匀时间网格下的Alikhanov型紧致有限差分格式,并证明了解的存在唯一性。第三章:利用卷积结构分析了时间分数阶Allen-Cahn方程紧致有限差分格式的收敛性和稳定性。第四章:结合差分格式,推导出三类加速方法,并通过数值实验验证了理论结果的合理性。第五章:对本文的工作做了总结和展望。
其他文献
研究背景:随着我国居民生活方式的改变以及生活水平不断提高,消化系统疾病发病率居高不下。目前在我国胃食管肿瘤、消化道溃疡、结直肠肿瘤等各类消化系统疾病发病率越来越高。电子胃镜及内镜下活组织检查是诊断各类消化道疾病的重要手段,但是因普通电子胃镜检查痛苦且技术水平要求高和部分人群接受度较差等局限性,一定程度上限制了其在人群中大规模应用,同时电子胃肠镜难以对小肠行进一步检查。胶襄内镜在2000年被研究发明
学位
20世纪80年代白羽肉鸡引入中国,40多年的行业发展发生了巨大变化,白羽肉禽从最初的满足出口到今天逐渐适应满足国人消费,从最初的手工养殖到现在全自动化操作,从最初的小规模养殖、屠宰到今天的大型化、全产业链模式。2014年中国鸡肉产量已经跃居全球第二。受利益追逐本能驱使,且行业本身进入门槛相对较低,2018年起,大量资金、技术等涌入到白羽禽肉行业,行业进入高周期运转,肉禽企业间的竞争也持续加剧。LH
学位
正Markov跳变线性系统是一类特殊的随机切换系统,可用来描述由于突然的环境干扰,不受控制的配置更改,子系统连接方式改变等引起的正线性子系统之间的随机跳变.本文分别讨论了连续时间和离散时间时滞正Markov跳变线性系统的性能与动态输出反馈控制器设计问题,采用线性Lyapunov泛函的方法给出了时滞正Markov跳变线性系统γ-指数均值内稳定意义下的L1(l1)增益性能的充分条件,及其正动态输出反馈
学位
2014年10月份,国家颁布新型显示产业创新发展行动计划,彰显了新型显示在信息产业中的地位。受到国家相关政策的影响,以京东方、深圳天马和华星光电为代表的大陆液晶显示屏生产企业在AMOLED显示技术领域的加快研发速度,2017年逐步开始进入AMOLED显示屏市场,实现批量生产。LY公司作为该市场的开拓者,正在面临市场新加入者的竞争,随着越来越多新竞争者的加入,AMOLED显示屏正在由蓝海市场转向红海
学位
本文所探讨的是(1+n)形中心构型共形解的稳定性问题.中心构型是天体力学中的重要概念,利用中心构型构造平面N体问题的共形解是已知的构造周期解的唯一方法.在这类周期解中,所有质点沿离心率相同的Kepler轨道运动并时刻保持相同的构型.如果Kepler轨道是椭圆的,这种解被称为椭圆相对平衡解,离心率e=0时,则称为相对平衡解.一直以来,相对平衡解的稳定性问题受到了广泛的关注.其中,(1+n)形中心构型
学位
华林-哥德巴赫问题一直是堆垒素数论的经典问题之一,其研究具有重大的理论意义.华林-哥德巴赫问题主要研究对于满足必要同余条件的正整数n,方程n=p1k+p2k+…+Psk是否有解,其中k为正整数,p1,p2,…,ps均为素数.记H(k)表示使得上述方程可解的最小的s.著名的哥德巴赫猜想可表述为H(1)=2,这个猜想至今仍未被证明.1937年,Vinogradov[2]利用圆法及素变量三角和估计证明了
学位
金融危机的过程中往往伴随着市场流动性的匮乏,金融危机在市场上一个最直接的体现就是资产价格的急剧波动,而市场流动性的匮乏会导致交易使得价格出现明显变化。期权作为金融市场的重要组成部分,常被投资者用于对冲风险提高收益,同时也稳定了市场经济。然而当期权的标的资产市场流动性短缺的时候,对应期权的价格也会受到影响,因此,将流动性风险因素加入期权定价的模型中,对于让期权定价模型更贴近真实的金融市场具有重要的参
学位
本文考虑半导体问题的Drift-Diffusion(DD)模型,它是由浓度方程和电子位势方程构成的一个方程组,浓度方程中不仅包含一阶导数对流项,还包含二阶导数扩散项。本文给出了 DD模型一维和二维问题的局部间断Galerkin(LDG)方法,并进行了数值模拟。主要技术困难包括由LDG方法的不连续性引起的单元边界处理,以及浓度方程中浓度与电场强度的非线性耦合。模拟一维问题时,在浓度变化剧烈的区域采用
学位
在本文中,我们研究的是应用二重网格有限元方法求解非线性多孔弹性体模型,问题模型如下:(?)对于非线性多孔弹性模型,我们一般采用Picard迭代法或者Newton迭代法进行求解,但是在剖分比较密的时候,计算效率非常低,为提高计算效率,本文提出应用二重网格有限元方法求解非线性多孔弹性模型。二重网格算法的主要思想是将一个在细网格上难于求解的非线性问题转化为一个在粗网格上易于求解的非线性问题和一个在细网格
学位
目的 观察结直肠侧向发育型肿瘤(LST)组织中RNF6蛋白的表达变化,并探讨其意义。方法 LST患者33例、结直肠隆起型腺瘤(PA)患者37例、结直肠癌患者35例均接受手术治疗,手术过程中留取切除的肿瘤组织分别记为LST组、PA组、CRC组,留取结直肠癌患者的癌旁正常组织(距癌组织>5 cm)记为Control组。采用免疫组织化学染色法检测各组织中RNF6蛋白,采用半定量评分法评估蛋白表达情况,包
期刊