若干非局部和非线性问题的数值方法

发布时间：2020-11-07 06:53

　　 随着科学和工程技术的发展,越来越多的领域提出了对复杂问题高性能数值计算的要求.比如对于机械工程、材料科学、航空航天和卫星遥感等在各种实际条件下的形态以及行为做出更精确的模拟及进一步的预测.对这些复杂系统做更精确的模拟则需要更复杂的如非局部和非线性的抽象模型,比如不满足Fick定律的反常扩散现象在自然科学和社会科学中是大量存在的,需要采用具有历史依赖和非局部性质的分数阶微分算子来进行模拟.这些非局部和非线性复杂模型在数值求解时又往往会形成一个大规模且稠密的离散系统,如何快速求解这样的大规模问题或者降低问题的规模成为当前日益关注的问题.下面介绍本文的工作.第一部分、分数阶微分方程的相关模型及其数值方法与应用在自然现象中,不混溶的液相、液相与气相、气相与固相、以及液相与固相是普遍存在的.不混溶相的界面动力学在多相系统的形成与相变机制以及演化中起着非常重要的作用,但是在建模和数值求解计算上都是一项非常困难的任务.有效描述多相系统的可计算数学模型的发展和相应的有效且高效的数值方法已经成为多相系统研究中的主要挑战.另一方面,由于分数阶微积分具有的历史依赖与非局部的特性比较适合诸如反常扩散中的记忆和非局部等性质,因此分数阶方程比整数阶方程更能有效的描述这些复杂系统.随着涉及的应用领域越来越多,分数阶微积分方程的研究逐渐成为一个新的活跃领域.在理论分析和数值模拟等多方面,分数阶方程为描述复杂现象提供新的视角和工具,但同时带来了很多新的挑战.我们考虑了如下的几个问题:1.相场模型的分数阶建模与快速算法·时空分数阶Allen-Cahn模型及其数值算法.Allen-Cahn模型可以用于描述如凝固及结晶等物理问题中多相变化过程,通过将经典模型中的拉普拉斯算子替换为分数阶拉普拉斯算子可以在界面厚度参数较大的情况下能够灵活地控制调节界面锐度来更精确地跟踪界面演化,同时还引入时间分数阶导数来描述模型的长时间记忆或延迟效应.我们提出了一个预条件的快速及节约存储的数值格式,计算量由直接解法的O(NV3)降至O(NVlogN),存储量由O(NV2)降至O(N),可以快速地求解大规模系统.并且初步研究了时空分数阶参数的变化对该模型的锐度和延迟效应的影响,数值实验验证了该模型具有良好的建模能力可以更好地跟踪界面演化.·空间分数阶Allen-Cahn模型及其二阶无条件能量稳定的数值算法.由于模型中包含的非线性势函数项和分数阶微分算子导致的复杂性,尚未有文献提出具有严格证明的无条件能量稳定的高效且精确的数值计算格式.而我们采用不变能量二次化方法,引入辅助变量来保证自由能密度是一个不变的二次泛函,得到等价模型,通过半显式处理非线性项,构造了一阶和二阶半离散格式,其方程组由每个时刻上的线性分数阶拉普拉斯方程构成,由于系数矩阵算子是对称正定的,从而可以有效求解半离散系统.我们进一步证明采用这种方法所提出的数值格式是无条件能量稳定的,在采用节点配置法进行空间离散得到全离散格式后,通过数值实验验证了该数值算法的稳定性和准确性.2.形状记忆聚合物分数阶建模与算法·形状记忆聚合物的变阶分数阶微分方程模型及其数值算法.形状记忆聚合物能够记忆其原始形状,可以在形变时获得暂时形状并且响应外部刺激而恢复到其永久形状.描述其形状记忆效应的早期模型是结合描述理想固体的胡克定律和描述理想流体的牛顿粘性定律得到的,而考虑到形状记忆聚合物的形状可以响应外部刺激温度的波动而产生巨大的变化,并反应在其微观网络结构上,因此变阶分数阶微分方程模型更适合描述形状记忆效应.对于该模型我们构造了数值格式求解给定可变阶数的正问题,而由于实际中可变阶数的未知性,我们还研究了其反问题,即通过物理数据来确定可变阶数,并通过自适应方法保证其精度.数值实验验证了模型的模拟效果及数值格式的有效性.第二部分、最优控制问题及其数值方法最优控制问题可以概括为:对一个受控的系统,从已有的控制方案中找出一个最优的控制方案,使得目标性能指标值达到最优.在一定的条件下,投入最少的成本,获得最大的收入和利益,这普遍存在于工程、金融、医药等众多实际应用领域的模型中.本文中我们研究了带有点态控制受限的受椭圆偏微分方程约束的最优控制问题的自适应有限元方法.·自适应有限元方法具有诸多优点,能够根据局部误差指示子自动判断并在误差最大处进行加密,提高了计算精度和效率,以及通过误差估计给出计算结果的误差范围等.最优控制问题自适应有限元方法的收敛率理论分析较少,而且已有算法需要外层自适应加密迭代和内层求非线性问题迭代两层迭代,对于内迭代的误差及大量的计算工作在于反复求解非线性问题等方面没有考虑.我们采用变分离散方式离散控制变量以及分片线性连续函数逼近状态变量,提出了一个收敛并节省计算量的自适应有限元算法,该算法仅含一层加密迭代,即内层非线性迭代仅需迭代一次,有效地节省了计算量.基于对该问题的后验误差估计以及数据振荡的缩减,得到了算法的整体收敛性并证明其具有线性收敛速率,并通过一些数值实验验证了理论结果.
【学位单位】：华东师范大学
【学位级别】：博士
【学位年份】：2018
【中图分类】：O241.8
【部分图文】：

１．１．２相场模型??在自然界中，不混溶的液相、液相与气相、气相与固相、以及液相与固相是??普遍存在的（图１．１）．不混溶相的界面动力学在多相系统的形成与相变机制以及演??化中起着非常重要的作用丨３６，?６１，?９６Ｕ旦是在求解计算上是一项非常困难的任务．??有效描述多相系统的可计算数学模型的发展和相应的有效且高效的数值方法已经??成为多相系统研究中的主要挑战．??比如在描述凝固及结晶等过程中常用的物理界面模型可以分为两类：自由边??界模型和相场模型．自由边界模型采用尖锐界面近似，即假定不同相的交界面厚度??为零，在整个区域内非此即彼，相变在瞬间完成，每一个相都遵守自己的守恒定律??和本构关系，不同相的相互作用由界面上施加的边界条件决定，不同相的物理量在??界面上可能是不连续的丨１７，?７３，?１６５］．关于自由边界模型已经有多种不同的数值方??法来数值模拟，包括前沿跟踪方法丨７３，?１７０］和移动网格法丨５１，?８６］等．在数值求解??时，自由边界模型需要将不同的相分开处理，并且还要跟踪界面的演化，这在数值??上是很困难的

图１．４：尖锐界而与扩散界而对比．??

当，一o.i,T二
【相似文献】

相关期刊论文 前10条

1 卢曰海;丘东元;张波;陈艳峰;;大功率分数阶电感的电路实现[J];电源学报;年期

2 王娴;李东;;分数阶混沌系统的间歇控制同步[J];重庆工商大学学报(自然科学版);2018年04期

3 李明;陈旭;郑永爱;;一类分数阶混沌系统的自适应滑模同步[J];扬州大学学报(自然科学版);2016年03期

4 程春蕊;朱军辉;毛北行;;一类分数阶混沌系统的投影同步[J];河南科学;2016年11期

5 郝孟丽;任勤;;标量控制下的分数阶Lü系统的参数辨识和自适应同步[J];河南理工大学学报(自然科学版);2017年01期

6 王志兰;;分数阶电路阶跃响应特性研究[J];电子测试;2016年24期

7 李亮;李庆宾;毛北行;;分数阶同步发电机系统的混沌同步[J];河南科学;2017年03期

8 李特;袁建宝;吴莹;;一类不确定分数阶混沌系统同步的自适应滑模控制方法[J];动力学与控制学报;2017年02期

9 赖满丰;金玲玉;;分数阶Klein-Gordon-Schr?dinger方程弱解的存在性[J];佛山科学技术学院学报(自然科学版);2017年03期

10 田魏巍;;非线性分数阶动力系统的控制研究[J];教育现代化;2017年22期

相关博士学位论文 前10条

1 Azmat Ullah Khan Niazi;分数阶中立型微分方程的稳定性和可控性[D];安徽大学;2018年

2 赵猛;时空分数阶偏微分方程的快速算法及其应用[D];山东大学;2018年

3 李政;若干非局部和非线性问题的数值方法[D];华东师范大学;2018年

4 练婷婷;Banach空间中分数阶发展系统的能控性与优化控制问题[D];扬州大学;2018年

5 姜翠美;分数阶复动力系统的定性分析与同步研究[D];山东大学;2017年

6 孙红;基于多项式插值逼近的分数阶偏微分方程高精度差分方法[D];东南大学;2016年

7 吴艳萍;分数阶非线性系统同步与非线性电路理论若干问题研究[D];西北农林科技大学;2015年

8 王乔;分数阶混沌系统控制与同步理论研究[D];浙江大学;2015年

9 丁策;机载光电稳定平台的分数阶控制研究[D];中国科学院研究生院(长春光学精密机械与物理研究所);2013年

10 张金营;分数阶混沌系统的鲁棒同步研究[D];华北电力大学;2014年

相关硕士学位论文 前10条

1 张薇;分数阶积分微分方程数值算法研究[D];电子科技大学;2018年

2 蒋伟;几类带边值条件的分数阶微分方程解的存在性研究[D];安徽大学;2018年

3 杜瑞连;时间分数阶扩散波方程的高精度数值解法[D];集美大学;2018年

4 刘娜;分数阶长短波方程组的吸引子和近似惯性流形[D];鲁东大学;2018年

5 李丹妮;分数阶热弹理论下材料特性温度相关一维问题的非局部动态响应[D];兰州理工大学;2018年

6 王西丽;具有积分及无穷多点的分数阶微分方程解的存在性[D];安徽大学;2018年

7 马璐娉;一类次线性分数阶Schr?dinger-Poisson系统解的存在性[D];兰州理工大学;2018年

8 李琪;材料特性参数变化的分数阶广义热弹问题研究[D];兰州理工大学;2018年

9 陈瑶;回火分数阶布朗-朗之万运动的扩散行为：局部化以及弹道扩散[D];兰州大学;2018年

10 孙晶;二维回火分数阶拉普拉斯的算法实现及数值分析[D];兰州大学;2018年

本文编号：2873596