分数阶相场问题的有限元数值模拟

来源 :山东师范大学 | 被引量 : 0次 | 上传用户:semitic
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
相场数学模型利用引入的纯量相函数,刻画了不同相物质组分的扩散运移与界面的形成与发展过程,在多相流体力学[17,85,86]、晶体生长[70.102.104]、固体力学[24,20,121]等重要自然过程中具有广泛的应用.经典的相场模型由Cahn-Hilliard对具双阱位势的Ginzburg-Landau自由能二次泛函通过变分形式获得,主要表达形式为四阶Cahn-Hilliard方程与二阶 Allen-Cahn方程.在经典相场模型的数学推导过程中,粒子间的相互作用关系由空间卷积的形式表达.这是一个非局部的积分形式,主要用来描述粒子间长程力的作用结果,但是为了数学推导的方便,被简化成了整数阶微分算子进行处理.因此,借助于分数阶算子刻画粒子间的长程力作用,从而建立分数阶相场模型并对其进行有效的数值模拟是十分必要的.本学位论文旨在对分数阶相场模型构造高效Galerkin有限元数值模拟方法并对其建立严谨的数值分析理论,为实际的相场模拟提供理论与决策依据.主要内容分为两个部分:一、分数阶扩散方程的有限元方法与分数阶导数指数的辨识基于对相场模型中分数阶Laplace算子可视为一类特殊的分数阶扩散算子的认知,我们首先考虑了变系数分数阶扩散问题的有限元方法.同时,为适应现场试验中须明确分数阶指数的需要,我们对非稳态分数阶扩散方程构造了关于分数阶指数的有限元辨识算法.通过这两个内容的研究,深入了解对分数阶问题乃至非局部问题数值模拟的难点与关键技术,积累经验,以期更好的实现对分数阶相场问题的高效有限元数值模拟.(1)具变系数分数阶扩散方程的最小二乘混合元数值模拟.由于变系数的影响,直接对该问题建立的Galerkin变分形式不能保证其强制性,从而不能保证解的存在唯一性.对此,我们通过引入扩散通量作为中间变量,将原方程转化为由变系数一阶方程与常系数分数阶方程构成的方程组,并借助最小二乘的思想,构造了具唯一解的最小二乘混合变分形式,证明了一个新的关于解的正则性结论.进一步,构造了可分别求解未知函数、伴随扩散通量的最小二乘混合有限元格式,有效降低了计算量,并证明了格式解对真解以及伴随扩散通量具有最优逼近精度.数值实验验证了理论分析结论.(2)分数阶扩散方程指数辨识的有限元-Armijo算法.与二阶扩散方程相比较,我们仅仅知道分数阶指数存在于某一个区间,如区间(1,2)中,但在实际的数值模拟应用中,需要清晰地给出确定这个分数阶指数的有限元算法与严谨的数学理论支持.对此,我们基于傅里叶变换理论,证明分数阶Riemann-Liouville导数算子满足某种意义下的半群性质,且关于指数在L2范数意义下满足强左连续性和弱右连续性;进一步,将指数辨识问题归结为某一个凸泛函在L2范数意义的最优解问题,并利用已获得的半群性质以及分数阶扩散方程解关于指数一致有界等结论,证明了分数阶扩散方程的解集构成了某一 Sobolev空间的弱闭子集,并证明了极小泛函实际上是一个定义在这个弱闭子集上的弱下半连续泛函,此极小泛函在该弱闭子集上必获极小值,从而证明最优解的存在唯一性.基于上述分析,我们采用普通的有限元空间设计了分数阶指数辨识的有限元-Armijo算法,并提供相应的数值实验佐证方法的可靠性.二、分数阶相场问题的有限元数值模拟及其理论分析本部分重点讨论一维分数阶Cahn-Hilliard相场模型的有限元数值模拟,并将该方法推广到二维分数阶Allen-Cahn相场模型与二维分数阶Cahn-Hilliard相场模型.(1)一维分数阶Cahn-Hilliard方程的混合有限元数值模拟与理论分析.我们通过引入中间变量,将原问题转化为由两个分数阶Laplace算子表示的非线性系统,构造了与之等价的混合变分形式;对时间导数采用向后Euler格式离散,而对空间采用混合有限元进行离散,构造了相应的凸分裂-混合有限元格式;利用Brouwer’s不动点原理,证明了混合元格式解的存在唯一性,进一步证明了格式解保持原始问题的能量衰减律以及在新能量定义下的能量守恒性质;利用离散Ehrling不等式,证明了格式在能量模意义下的收敛性质.为降低由分数阶Laplace算子非局部性导致的有限元刚度矩阵非稀疏以及带来的计算困难,我们将刚度矩阵的Toeplitz块结构、快速傅里叶变换(FFT)与传统的双共轭梯度(BiCG)巧妙结合,构造了高效求解该类问题的快速双共轭梯度算法(FBiCG),计算量由BiCG的O(M2)降低到O(M log M)数值实验也表明,本文中提出的凸分裂-混合元格式不仅具有良好的收敛精度与计算效率,而且保持了原始问题的能量衰减律、新能量下的守恒律,可清晰的刻画出物质各组分的粗化过程以及分数阶指数、扩散系数等参数对各组分界面发展的影响.(2)二维分数阶Allen-Cahn模型有限元方法与理论分析.我们在Galerkin框架下构造了相应的变分形式与有限元格式,证明格式的可解性、稳定性、能量定律与收敛性结论.数值算例验证了格式的收敛性、能量定律和各组分的粗化过程.(3)二维分数阶Cahn-Hilliard模型的有限元算法与实现.本章中,我们将第三章中关于一维问题的结论推广到二维问题,构造了凸分裂-混合元离散格式,证明格式解的存在唯一性、能量定律与收敛性.数值实验表明文中提出算法的有效性.有必要指出的是,与一维问题比较,对高维分数阶相场问题的计算难度显著增加,其最大困难在于刚度矩阵形成过程中的关于多重奇异积分的数值计算.
其他文献
学位
近年来,细菌耐药性的出现和传播给临床抗感染治疗带来了巨大挑战,超级耐药细菌对于人类健康的威胁日益得到重视,能治疗耐药细菌感染的噬菌体疗法重新引起了人们的广泛关注。然而,在噬菌体治疗中,细菌很容易对噬菌体产生抗性,极大地阻碍了噬菌体治疗耐药菌的发展。因此,探究细菌快速产生烈性噬菌体抗性的机制并克服细菌对噬菌体的耐受成为了噬菌体治疗细菌感染领域的首要问题。理论上,与抗生素不同,噬菌体是一种病毒,它可以
学位
发酵床网上养鸭作为一种新型的养殖模式越来越受养殖场户欢迎,但新模式在技术应用方面还有许多细节需要完善熟化,针对模式应用过程中发酵床承载力过度、设施设备不匹配的相关问题,本文开展了发酵床网上养鸭棚舍及设施现状调查并提出相关改进优化建议,以供参考。
生成对抗网络(generative adversarial network, GAN)快速发展,并在图像生成和图像编辑技术等多个方面取得成功应用。然而,若将上述技术用于伪造身份或制作虚假新闻,则会造成严重的安全隐患。多媒体取证领域的研究者面向GAN生成图像已提出了多种被动取证与反取证方法,但现阶段缺乏相关系统性综述。针对上述问题,本文首先阐述本领域的研究背景和研究意义,然后分析自然图像采集与GAN
学位
个人信息最小化原则为世界各地广泛采用。但其实践成效,日益受到质疑。对此,须对原则作更精细、多元分析。"精细"层面,将原则分解为"处理目的""信息类型""类型与目的间关联""处理对个人权益的影响"四部分,阐明各部分对应的基础性困难;"多元"层面,考虑数字经济中"信息作为价格"的现实,阐明最小化原则与备受争议的限价规制的相似性。在此基础上,就分析信息类型与目的间是否"直接相关",提出以下细化原则理解和
自1977年恢复高考以来,特别是经过40多年的改革开放,我国高等教育事业蓬勃发展,取得了令人瞩目的成绩:人才培养规模不断扩大,办学层次更加多元,科研创新水平不断提高,支撑服务经济社会发展的能力明显加强。特别是1999年开始,我国高等教育特别是地方本科高校通过扩大招生人数、扩大办学规模的外延式发展模式,实现我国高等教育的迅速扩张。这种外延式的发展模式在当时的历史阶段起到了一定的积极作用。目前我国已经
目的:观察分析面神经炎患者糖皮质激素联合神经节苷酯治疗的有效性。方法:选择2018年11月~2020年11月某院收治的面神经炎患者62例,随机分为对照组和研究组各31例,对照组应用糖皮质激素治疗,研究组在此基础上应用神经节苷酯治疗。观察两组患者的面神经功能恢复时间和面神经运动恢复时间以及不良反应发生率,观察两组患者的治疗效果和治疗前后的生活质量评分情况和面神经分级情况。结果:研究组面神经功能恢复时
2021年9月29日,2020/2021中国城市规划年会暨2021中国城市规划学术季专题会议二十"活力交通与空间营造——线上交通专场(青年托举)"顺利召开。本次会议由中国城市规划学会主办,中国城市规划学会城市交通规划学术委员会承办,同济大学交通运输工程学院副教授李健主持。会议以作者宣讲论文、行业专家点评讨论的形式共议城市活力交通发展,就大数据环境下活力交通跨学科研究、活力交通的空间保障与品质提升、