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

来源 :山东大学 | 被引量 : 0次 | 上传用户:wzxisno2
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
近几十年来,分数阶微积分理论成功应用于工程科学的各个领域,比如电磁学、流体力学、粘弹性、反常扩散和信号处理,呈现出欣欣向荣之势。这表明分数阶微积分理论具有独特的优势,体现了其不可替代性,关于这方面的理论和应用研究已成为一个国际热点问题。相比较经典的整数阶微分算子,由于分数阶微分算子具有全局相关性或非局部的特性,故更适合于描述具有记忆性和遗传特性材料的复杂力学行为。本文主要介绍分数阶微积分理论在粘弹性材料、流体力学、生物组织传热和激光加热中的应用。为了更好的分析这些反常现象,通过合适的参数估计方法获得了模型中未知参数的最优估计结果。首先,利用由分数元组成的分数阶本构关系模型:分数元模型,分数阶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,0<y<L,(6)u(0,t)= u(L,t)= 0,t>0.(7)借助于有限差分的方法获得此问题在常压力梯度驱动下的速度分布,探讨分数阶导数和广义雷诺数对不定常粘性流动的影响,结果显示两个模型参数对两平行平板之间的流动均有较大影响。为了更好的分析两平行平板之间的压力驱动流动的特性,探讨两个模型参数的估计问题,并利用Levenberg-Marquardt(L-M)方法估计了分数阶导数和广义雷诺数,结果表明用Levenberg-Marquardt方法解决空间分数阶微分方程中的反问题是有效的。第四章,主要研究时间分数阶双相延迟模型及其在生物组织传热中的应用。首先基于经典傅里叶定律提出以下生物组织传热中的时间分数阶双相延迟模型(?)其中0<α,β<1,分数阶微分为Caputo型分数阶微分。生物组织热传导通常采用以下Pennes’方程(?)其中Q(r,t)=ωbCb(Tα-T(r,t))+qm+qr,ρ,c和T分别是生物组织的密度,比热和温度;cb是血液的比热,ωb是血液灌注率;Ta是动脉血的温度;qm是生物组织的新陈代谢产生的热量,qr是空间热源项。方程(8)和(9)联立可得带有两个分数阶参数α和β的生物传热方程,其中包括体现生物组织内在热性能的松弛时间τq和延迟时间τT。基于Mitra等的实验,提出以下初边值条件:(?)(?)(?)并借助于积分变换和H函数,得到了以上初边值条件下分数阶双相延迟热传导方程的精确解。进一步用非线性最小二乘的方法预测了 Mitra等给出的实验Ⅰ和Ⅲ的温度,结果显示分数阶双相延迟模型能较好地吻合测量数据且能很好的抓住温度快速上升的趋势。对比热波模型和双相延迟模型,分数阶双相延迟模型的L2误差最小。这说明分数阶双相延迟模型有助于提升对生物组织中热传导现象的认识。第五章,主要研究短脉冲激光加热一个半无穷介质。以时间分数阶双相延迟模型(?)为热传导模型,建立相应的包含体积热源项的分数阶热传导方程(?)这里g(x,t)=(1-rf)I0δf(t)e-δx.气考虑初边值条件(?)(?)(?)(?)T(x,t)= T0,x→,t>0,(16)下定解问题的解,通过拉普拉斯变换的方法得到了温度分布的半解析解。最后,利用数值逆拉普拉斯变换的方法数值剖析了分数阶阶数和延迟时间对温度分布产生的影响。第六章,给出本文的总结和未来可能的研究方向。
其他文献
以核心素养推进教育改革与发展是当今世界教育发展的热潮。中国学生发展核心素养以培养"全面发展的人"为核心,充分反映国家经济社会发展对人才培养所需要的关键能力和必备品
本文介绍了一起500kV智能变电站为了实现大二次遥控功能,进行保护程序升级时发现保护厂家提供的新ICD文件与保护装置的CID文件的GOOSE虚端子差异较大,防止了不合格下装保护装
股权结构是公司治理的产权基础,当企业的控制权发生转移,公司治理机制将会发生变动。本文以2004年标的为上市公司股权(标的属广东省)且控制权发生转移的股权交易案为研究样本
实施生态扶贫移民搬迁工作,是加快贵州经济社会发展的重要举措。到2020年,贵州省将完成对200万贫困群众的搬迁工作,规模为全国之最。但是,生态移民搬迁不仅仅是简单的生活地
传统KNN查询是一种稳定性和准确率性能均较好的算法,但是在样本规模过大时,算法的计算效率受到影响较大,对此提出一种基于聚类中心文本串联的并行(Mapreduce for KNN,MKNN)文
园林绿化是一座城市的重要组成部分,发挥着装点城市、净化城市空气和缓解噪音污染等作用。而园林中各类植物的种类与数量的配置则影响着园林绿化是否能较好的发挥作用,也决定
对220名学生第二年和第四年两年体质测试数据进行比较分析结果表明:速度和耐力素质在第四年明显下降,肺活量低于本省其他学校,力量性项目没有差异性变化。为此建议:加强速度和耐力素
农家书屋逐步成为农村社区文化娱乐中心,在农村留守儿童问题日益突出的今天,农家书屋应积极发挥自身优势,为农村留守儿童健康成长搭建平台。
<正>3月11日消息,来自ABI Research的最新研究报告显示,2013年全球LTE用户数已达到2.297亿,并且2013—2019年将以43.6%的年复合增长率继续增长,到2019年将超过20亿。截至2013
期刊
罐车是一种车体呈罐形,主要用来装载各种液体、液化气以及粉末状货物的车辆。油罐车罐体的截面形状一般有四种:圆形、椭圆形、圆矩形及异形。油罐车罐体选择哪种界面型式,是需要