分数阶微积分在反常输运过程中的应用研究

作     者：续焕英 

作者单位：山东大学 

学位级别：博士

导师姓名：蒋晓芸

授予年度：2017年

学科分类：07[理学] 070102[理学-计算数学] 0701[理学-数学] 

主      题：解析解 有限差分算法 分数阶微积分 分数阶双相延迟模型 分数阶热传导模型 分数阶Poynting-Thomson模型 分数阶粘弹性模型 内点算法 激光加热 积分变换 Levenberg-Marquardt算法 非傅里叶热传导 参数估计 Riesz分数阶导数 粘弹性材料 

摘      要：近几十年来,分数阶微积分理论成功应用于工程科学的各个领域,比如电磁学、流体力学、粘弹性、反常扩散和信号处理,呈现出欣欣向荣之势。这表明分数阶微积分理论具有独特的优势,体现了其不可替代性,关于这方面的理论和应用研究已成为一个国际热点问题。相比较经典的整数阶微分算子,由于分数阶微分算子具有全局相关性或非局部的特性,故更适合于描述具有记忆性和遗传特性材料的复杂力学行为。本文主要介绍分数阶微积分理论在粘弹性材料、流体力学、生物组织传热和激光加热中的应用。为了更好的分析这些反常现象,通过合适的参数估计方法获得了模型中未知参数的最优估计结果。首先,利用由分数元组成的分数阶本构关系模型:分数元模型,分数阶Maxwell模型,分数阶Kelvin-Voigt模型和分数阶Poynting-Thomson模型,描述粘弹性材料的时间依赖蠕变行为。借助于聚合物和岩石的三组蠕变实验数据对比这些分数阶本构关系模型的有效性,并通过内点算法求得这些模型参数的最优估计结果,借助图形将三组材料的蠕变数据与计算所得结果进行对比,分析表明分数阶Poynting-Thomson模型在描述材料的蠕变行为时是最优的。其次,考虑到空间分数阶微分算子可以准确的描述反常力学行为的路径依赖、长程相关等特性,故将流体力学中的经典Navier-Stokes方程中拉普拉斯算子替换为Riesz分数阶微分算子,得到空间分数阶Navier-Stokes方程。借助分数阶微分方程的有限差分算法研究两平行平板之间的压力驱动流动,分析模型参数对流体流动的影响;并运用Levenberg-M arquardt算法获得模型中未知参数的最优估计值。结果显示,两个模型参数对速度场均有较强的影响,且Levenberg-Marquardt方法在空间分数阶微分方程反问题的研究中是有效的。再次,基于生物组织建立分数阶双相延迟模型和相应的生物传热方程,以此解释预处理肉中的传热现象。借助于傅里叶变换和拉普拉斯变换获得由H函数表示的解析解。模型中的两个松弛时间和分数阶导数的阶数均由非线性最小二乘法获得最优估计值。通过图形对比测量温度和计算所得温度,得到较好的拟合效果,并详细分析了模型参数的影响。最后,基于分数阶Taylor级数和Tzon提出的双相延迟模型,建立分数阶双相延迟热传导模型;针对短脉冲激光加热半无穷介质问题建立分数阶双相延迟热传导方程,利用拉普拉斯变换给出问题的半解析解,并利用数值求逆拉普拉斯变换的方法深入剖析了短脉冲激光加热的热传输过程。具体来讲:第一章,主要介绍分数阶微积分的发展历史,本文要研究的主要问题以及可能用到的一些预备知识。第二章,基于分数阶本构关系模型研究粘弹性材料的时间依赖蠕变行为。基于分数阶微积分的粘弹性材料本构关系模型不仅具有拟合实验数据好、使用参数少的优点,而且可以合理地描述具有记忆和时间依赖的物理现象,故研究人员在粘弹性材料本构模型的研究中开始采用分数阶微积分理论。本章借助于分数阶本构关系模型:分数元模型,分数阶Maxwell模型,分数阶Kelvin-Voigt模型和分数阶Poynting-Thomson模型研究粘弹性材料(聚合物和岩石)的蠕变行为,这些模型的蠕变函数分别为(?)(?)(?)(?)这里Ep,q(z)是Mittag-Leffler函数。通过包括HDPE,PEEK和岩石在内的三组蠕变实验数据说明分数阶粘弹性模型的有效性,并分析了所得拟合结果。结果显示,分数阶Poynting-Thomsom模型的拟合效果是本章涉及的分数阶微分模型中最好的,其可以很好的抓住粘弹性固体短时和长时的蠕变行为。模型参数的识别是分数阶微分模型研究中的一个重要问题,本章将采用内点算法估计了未知的模型参数,得到未知参数的最优估计值。内点算法对于解决分数阶微分模型参数估计的反问题是可行的。第三章,主要研究空间分数阶Navier-Stokes方程,并用分数阶微分方程有限差分的方法获得两平行平板之间压力驱动流动的速度分布。将Navier-Stokes方程中的拉普拉斯算子替换为Riesz分数阶微分算子,可得以下形式的空间分数阶Navier-Stokes方程(?)(?)以上述方程为运动方程,考虑垂直距离是L的两平行平板之间流体的不定常压力驱动流动,其初边值条件分别为u(y,0)= 0,00.(7)借助于有限差分的方法获得此问题在常压力梯度驱动下的速度分布,探讨分数阶导数和广义雷诺数对不定常粘性流动的影响,结果显示两个模型参数对两平行平板之间的流动均有较大影响。为了更好的分析两平行平板之间的压力驱动流动的特性,探讨两个模型参数的估计问题,并利用Levenberg-Marquardt(L-M)方法估计了分数阶导数和广义雷诺数,结果表明用Levenberg-Marquardt方法解决空间分数阶微分方程中的反问题是有效的。第四章