分数阶偏微分方程高效数值算法及其参数估计

来源 :山东大学 | 被引量 : 1次 | 上传用户:crazyliang
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
分数阶微积分是一个古老而新鲜的话题,在早期,由于缺乏物理机理解释、应用背景研究等原因,分数阶微积分发展缓慢.后来随着科学技术的发展,关于分数阶微积分的研究内容由理论研究逐步转向应用研究,比如其在高能物理、反常扩散、生物医学工程、系统控制及复杂粘弹性材料力学本构关系中的应用.近年来,许多学者对分数阶微分方程的求解进行研究,但由于分数阶导数定义中的微分-积分卷积算子的存在,使得其具有历史依赖性和全局相关性的性质,因此直接求解分数阶偏微分方程的解析解比较困难.同时,根据实际问题所构建的分数阶微分方程模型,关于方程中各个参数的变化对模型的影响的研究,对于分析问题内在物理机理及力学机制等具有重要的研究意义.但是,目前对于分数阶模型中的参数估计问题,主要采用的方法是“试差法”,这种方法效率低下同时相对不稳定.因此,对分数阶微分方程提出高效数值求解方法,并基于正问题的数值方法,提出合适的模型参数估计方法是一项有意义的研究.本文主要研究分数阶微分方程的高效数值求解方法及其参数估计问题.首先,对于具有周期性边界条件的时空分数阶电报方程提出时间方向Crank-Nicolson有限差分、空间方向傅里叶谱方法的数值格式,并对所提格式给出严格的收敛性与稳定性理论分析,数值实现过程中,借助快速傅里叶变换的技巧,极大地提高了运算效率.其次,我们对于具有齐次边界条件的二维分数阶加热广义二阶流体Stokes’第一问题提出时间方向L1有限差分,空间方向勒让德谱方法的数值格式,并考虑到方程具有非光滑解的情形,通过提出添加校正项的方法得出合适的数值解,同时给出所提方法严格的稳定性与收敛性理论证明.第三,借助于分数阶Maxwell本构方程,基于人体血液实验数据,采用非线性共轭梯度的方法给出本构方程中的参数估计值,极好的拟合实验数据.同时,结合连续性方程与柯西动量方程,构建了分数阶Maxwell流体电渗流流速的控制方程,并对控制方程提出一种高效的数值解法.第四,基于煤基质中甲烷气体吸附与解吸过程实验数据,借助分数阶菲克定律构建了多孔介质中的气体反常扩散模型,并提出模型的高效数值解法.将BFGS方法运用到分数阶模型的参数估计中,同时对比分析了BFGS方法、非线性共轭梯度方法及贝叶斯方法在参数估计中的优劣,研究结果表明,分数阶分形扩散模型可以很好的表示煤炭颗粒中甲烷气体的解吸过程.第五,发展了关于时间分数阶Cable方程的新的数值模拟技术,提出了基于勒让德多项式及其运算矩阵的高效时空谱tau算法,并给出了数值格式的收敛性证明.同时,利用非线性共轭梯度法建立方程的分数阶参数估计模型,并针对分数阶阶数的估计问题进行稳定性与收敛性分析,完善了分数阶模型参数估计反问题的理论框架.具体地:第一章,我们简要介绍分数阶微积分,并给出本文中用到的分数阶微积分的几种定义.然后,简要概述本文的主要研究内容.第二章,我们研究了具有周期性边界条件的时空分数阶电报方程的数值求解方法及其参数估计问题.对于这一问题,我们首先构造了时间方向Crank-Nicolson(CN)有限差分,空间方向傅里叶谱方法的数值离散格式,并在数值实现过程中,采用了快速傅里叶变换(FFT)的技巧来提高计算效率,同时给出离散格式的稳定性和收敛性的严格证明.对于参数估计问题,我们借助于正问题所提的CN有限差分傅里叶谱方法,提出Levenberg-Marquardt(LM)方法进行参数估计,通过添加随机扰动项与否来数值模拟实验数据的方式,对模型中分数阶阶数进行估计.最后,我们给出数值算例,数值结果表明,我们所提数值方法是无条件稳定的,该格式在时间方向可达到2-α阶精度,空间方向可以达到谱精度.同时,对于不同的初始猜测值和随机扰动误差,采用LM方法均可得到合理的参数估计值.这验证了所提方法的有效性和理论分析的合理性.第三章,我们对具有齐次边界条件的二维分数阶加热广义二阶流体Stokes’第一问题进行研究.对于这一问题,我们构造了时间方向采用L1有限差分,空间方向采用勒让德谱Galerkin方法的数值离散格式.借助所采用的L1有限差分格式对时间分数阶导数进行逼近时系数的性质关系,给出离散格式的稳定性与收敛性分析.同时,考虑了方程具有非光滑解的情形,并提出添加校正项的方法,将低正则项进行校正处理,给出合理的数值解.最后,对于不同的分数阶阶数γ,通过数值算例,我们验证了数值解与精确解之间分别关于时间方向和空间方向的误差、收敛阶和计算所需的CPU时间,结果表明我们所提格式在时间方向收敛阶为1+γ,在空间方向可达到谱精度,这也验证了理论分析的合理性.第四章,我们研究了矩形微通道中分数阶Maxwell流体的电渗流动.非牛顿流体在微通道中的应用主要与生物流体的研究有关,其中,最常用的应用是对血液的研究.在这一章节,我们基于人体血液耗散与存储模量实验数据,借助分数阶Maxwell本构方程,提出非线性共轭梯度的方法对方程中分数阶参数α,β,松弛时间τ和剪切模量E进行参数估计,为了得到精确的参数估计值,我们提出一种新型对数误差函数表示形式.所得结果与实验数据高度吻合,表明分数阶Maxwell模型可以极好的用来模拟人体血液的特性.其次,我们以上述模型为基础,建立了电渗驱动下的粘弹性流体动力学模型,构建了关于分数阶Maxwell流体电渗流流速的控制方程,并采用时间方向有限差分,空间方向勒让德谱配点的数值求解方法求出流体的速度场.研究结果对于微流控装置的设计及微通道中流体粘弹性行为的预测具有重要的指导意义.第五章,我们基于实际实验,研究了煤基质中甲烷气体吸附与解吸过程.由于煤炭颗粒对甲烷气体具有吸附性,同时煤基质孔隙具有非均匀性,其孔隙大的有几百微米,小的只有几十纳米,甚至几埃米.在这样的多孔介质中,气体分子的扩散路径受到扭曲,分子与壁面也会发生碰撞,而且在吸附或者解吸过程中,气体的压力随时间发生改变,压力变化引起平均自由程的改变,自由程与扩散路径有关,所以扩散过程通常呈现反常扩散.在这里,我们基于煤炭颗粒的这些特性,构建了煤基质中甲烷气体解吸过程的分数阶分形扩散模型,并提出时间方向有限差分、空间方向勒让德谱配点法的方法给出模型的数值解.然后,我们基于正问题的数值方法,对模型中分数阶阶数α、分形维数df、结构参数θ进行估计,并将BFGS方法运用到分数阶模型的参数估计中,同时对比分析了 BFGS方法、非线性共轭梯度方法及贝叶斯方法在参数估计中的优劣,所研究结果极好的拟合了煤炭颗粒中甲烷气体的解吸过程,为地下煤矿中瓦斯气体的开采及煤矿安全生产的研究提供理论支撑.第六章,我们研究了时间分数阶Cable方程正反问题的新的数值模拟技术.首先,对于时间分数阶Cable方程的正问题,我们提出基于转换的勒让德多项式及其运算矩阵的时空勒让德谱tau方法,给出了正问题的数值解,并证明了数值算法的收敛性.其次,以分数阶导数阶数为待估计参数,利用非线性共轭梯度法建立时间分数阶Cable方程的分数阶参数估计模型,并基于正问题的数值方法,通过仿真试验对分数阶参数进行了估计,讨论了待估计参数的初始猜测值、测量误差等对参数估计结果的影响.同时,对参数估计方法的稳定性和收敛性给出理论证明,完善了分数阶模型参数估计的理论框架,并为前面几个章节中参数估计在模型中的应用提供理论支撑.最后,给出数值算例验证所提方法的有效性及理论分析的合理性,并验证了时空勒让德谱tau方法对于模型中边界条件处理的灵活性.第七章,给出本文的总结和未来的研究工作展望.
其他文献
1.研究背景动脉粥样斑块钙化(atherosclerotic calcification,AC)是在动脉粥样硬化负荷斑块内发生的钙化,是一种最常见的钙化性血管病变,在炎症因子和氧化脂质等病理性刺激作用下,斑块内血管细胞亚群发生成骨样分化而形成的。动脉粥样斑块钙化与斑块的不稳定性密切相关,是心血管疾病的独立危险因素。糖尿病微环境中存在着普遍的促钙化刺激,如炎症、氧化应激损伤、高糖和无机磷酸盐等,可促
角膜是一层透明无血管的纤维膜,在组织学结构上角膜共分为5层,从外向内依次为:上皮细胞层,前弹力层,基质层,后弹力层和内皮细胞层。角膜内皮细胞为多边形单细胞层,是维持角膜相对脱水状态,保持角膜透明性和维持角膜厚度的重要因素。在人体内角膜内皮细胞难以再生,若被损伤只能依赖残存细胞形态的改变和移行来进行修复。当损伤超过一定限度造成内皮细胞的密度低于约500-800/mm2时,残存的内皮细胞便不能完全代偿
纤维素是自然界中丰富的可再生资源,利用纤维素酶降解纤维素,转化获得生物燃料和化工产品,能够保护环境、促进经济可持续发展。在纤维素酶生产菌株中,丝状真菌能够分泌完整的纤维素降解酶组分。探索其纤维素酶合成调控的关键因子,阐明纤维素酶合成的调控网络,将为菌株的遗传改造提供合理的靶点,对于提高酶的产量和降低成本具有重要的现实意义。据报道,假定的甲基转移酶LaeA/LAE1调控丝状真菌发育、纤维素酶基因及次
丙二酰辅酶A(malonyl-coenzyme A)是脂肪酸类衍产物,黄酮类化合物,聚酮类化合物合成的关键前体化合物,这些化合物在医药类,绿色能源及化工领域有重要的开发价值,因此构建丙二酰辅酶A平台菌株在工业生产中具有巨大潜力。酿酒酵母的遗传背景研究透彻,安全性高,基因操作手段成熟,是微生物工业合成的常用菌株,以酿酒酵母为细胞工厂合成丙二酰辅酶A衍生产品具有很大发展潜力。在酿酒酵母中,中间代谢物丙
研究背景动脉粥样硬化(Atherosclerosis,AS)是一类重要的动脉硬化血管病变,其病变从动脉内膜开始,先后出现脂质积聚、出血和血栓形成、纤维组织增生和钙质沉着、动脉中层退变和钙化。免疫炎症反应、氧化应激反应、血流动力改变等多个方面影响AS发生发展,机制尚未完全明确。探讨影响AS发生发展的分子机制,明确潜在治疗靶点,对于早期干预斑块形成及稳定易损斑块、降低不良心血管事件的发生具有重大意义。
硫是组成生命有机体的重要元素之一,在自然界中主要以硫化氢和硫酸盐的形式存在。硫化氢(H2S)作为一种气体信号分子可以参与生物体中的很多生理活动,先后有研究指出硫化氢在血管舒张、心脏保护、抵抗氧化应激、神经传递和抗炎症作用等过程中发挥重要作用。然而,过量的硫化物对细胞是有毒的,可以通过毒害细胞色素c氧化酶来抑制细胞呼吸或通过产生不同的活性硫物质来提高氧化应激。生物细胞可以自身调节胞内硫化物含量,把硫
酯类水解酶能够催化酯键的水解和合成,其广泛存在于植物、动物和微生物中。由于具有广泛的底物特异性、精确的化学选择性和区域选择性、反应不需要辅因子以及在有机溶剂中高度稳定等特性,酯类水解酶在工业上具有十分广泛的应用。糖酯酶是一大类糖类活性酶,能够催化酯类取代基从多糖的糖链上脱离,这些多糖来自于不同生物体,包括微生物多糖肽聚糖、动物多糖几丁质以及植物多糖木聚糖等。一些酯类水解酶还具有糖酯酶的活性。海洋中
肽聚糖是细菌细胞壁特有的基本组成成分,在维持细胞特定形状、保护细胞免受内部膨胀压破坏和参与细菌的生长和分裂过程中起重要作用。阐明细菌细胞壁肽聚糖的组织结构,是微生物学研究中的基本问题。解释清楚其结构对于研究细菌的生长和分裂机制、胞外聚合物与细胞壁的锚接形式、抗生素靶点的选择以及仿生学等都有重要意义。肽聚糖是由聚糖链和肽链组成的种类繁多、变化复杂的超分子复合物。到目前为止,细菌细胞壁肽聚糖主要的拓扑
自夸克模型和量子色动力学(QCD)理论的相继建立,直至标准模型的完善,人们在认识微观粒子的深层结构和性质上已经取得了极大成功。由于QCD理论在低能区域的非微扰性质,至今仍不能从第一性原理出发直接研究强子结构和强子谱等相关问题,尤其对于奇特强子而言,非微扰问题变得更加突出,因此很多研究中会使用基于QCD精神的有效理论以及相关唯象模型以分析这些问题。近些年来,实验上陆续发现了一系列性质奇特的介子和重子
第一部分人副流感病毒3型F蛋白HRB连接区结构与功能研究研究背景:人副流感病毒3型(human parainfluenza virus type3,HPIV3)为副粘病毒科(Paramyxoviridae)副粘病毒属(Paramyxovirus)病毒家族中的一种重要病毒,可引起5岁以下婴幼儿下呼吸道感染,导致支气管炎和肺炎。在病毒性呼吸道感染患者中,HPIV3的临床阳性检出率仅次于呼吸道合胞病毒(