几种有限体元格式及其在辐射热传导问题中的应用

来源 :湘潭大学 | 被引量 : 6次 | 上传用户:matlab_walker
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
有限体积法由于能保持某些物理量(如质量、能量等)的局部守恒性,它已成为一种广泛应用于科学与工程计算领域的重要数值方法.有限体元(FVE)法是一种重要的有限体积法,虽然对它已有大量的研究工作,但仍存在许多需要研究的问题.ICF问题描述的是一种聚变等离子体流体力学问题,其中三温辐射热传导方程的数值求解在该问题的数值模拟中占有主要工作量,因此,对该类问题的有限体元格式及其快速算法进行研究,是一项具有重要理论和实际应用价值的工作.本文针对定常扩散和二维辐射热传导等问题构造和分析了几种有限体元格式,获得了一批算法和理论结果.针对定常扩散问题,当扩散系数连续时,在一致四边形网格上,首次从理论上严格证明了逐点意义下等参双线性有限体元解的渐近展式和超收敛结果;当扩散系数间断时,在含混合单元的四边形网格上,构造了两种等参双线性有限体元格式;数值实验验证了理论结果的正确性,表明了混合网格上的第二种格式具有饱和收敛阶.针对线性抛物问题,设计了相应的等参双线性有限体元格式,数值结果表明有限体元解在逐点意义下存在超收敛性.针对连续系数和间断系数两种情形的定常扩散问题,在非结构四边形网格上,分别建立了两种保对称有限体元格式,对前一情形,当网格拟一致时,从理论上证明了有限体元解在L2和H1范数下均具有饱和收敛阶;针对一类具有非局部边界的二维椭圆问题,设计了一种保对称有限体元格式,并得到其最优L2模估计.数值实验验证了理论结果的正确性,表明了新格式具有一些特点:对扭曲网格、间断系数问题的较强适应性,对流函数在正交网格单元中心点的超逼近性,混合网格上的第二种格式具有饱和收敛阶.针对连续系数和间断系数两种情形的定常扩散问题,分别在四边形网格上建立了相应的二阶混合有限体元格式,通过建立并从理论上证明了分层基下该有限体元刚度矩阵和二次有限元刚度矩阵的谱等价关系,设计了一种高效GAMG预条件子;数值实验验证了理论结果的正确性和新预条件子的高效性与稳定性,表明了混合网格上的第二种格式具有饱和收敛阶.针对柱坐标系下多介质输运管单温辐射热传导模型问题,首先设计并分析了相应的分片线性保对称有限体元格式,接着将上述的二阶混合有限体元格式应用于该问题的数值求解,数值实验表明:经新格式计算得到的热传导过程与实际物理现象相符,且具有较小的能量守恒误差.针对含和不含混合单元两种情形对应的轴对称辐射热传导问题以及含混合单元情形的球截面辐射热传导问题,在Euler四边形网格上建立了相应的保对称有限体元格式和等参双线性有限体元格式,数值结果表明:经这些格式计算得到的温度分布和辐射热传导过程比较合理,且能量守恒误差较小,从而验证了这些格式在求解该类问题时的有效性.
其他文献
软件定义网络(Software-Defined Networking,SDN)是一种新型网络体系结构,它实现了网络设备中控制逻辑与数据转发功能的分离,采用集中式的控制方式控制整个网络,并向上层应用提供开放的可编程接口。这种设计模式极大地简化了网络策略的部署,缩短了网络应用的开发周期,在数据中心和云计算中得到了广泛的应用。但是,这种网络体系结构也引入了新的安全威胁。本文将对SDN体系结构中的若干安全
热障涂层(Thermal Barrier Coatings,TBCs)因为耐高温、高隔热、抗腐蚀等优异性能,已成为航空发动机涡轮叶片等热端部件不可缺少的热防护材料。然而,在含有杂质颗粒的高温、高速燃气环境下服役时,会因为界面氧化、颗粒冲蚀、CMAS腐蚀(钙镁铝硅等金属氧化物的混合物,简称CMAS)等多种因素导致涂层剥落。其中,CMAS腐蚀指高温下钙镁铝硅混合氧化物熔融、渗透到涂层中,引起涂层结构、
本文利用谱配置方法求解一类带弱奇异核的和非线性的Volterra型积分微分方程,并且构造高精度算法,着重分析该方法的误差估计和收敛性,并进行数值实验,验证所给方法的有效性。第二章,用Legendre谱配置方法对非线性Volterra-Fredholm-Hammerstein积分方程进行求解。并且对于所给的数值格式和误差进行了分析,即当核函数充分光滑时,计算所得数值格式的L2范数和L∞范数误差呈指数
半拉格朗日(半拉氏)方法没有CFL稳定性条件限制且比欧拉方法更加稳定。在文献当中,水平集重新初始化方程的半拉氏方法是复杂的,这是没有必要的。因为重新初始化过程是辅助的,我们提出了耦合有改善紧邻交界面网格点处精度的投影技术的一阶半拉氏格式。标准的二阶半拉氏格式被用于推进水平集对流方程。实现是简单的,包括在块结构自适应网格上。半拉氏方法的有效性通过大量的数值算例被证实,其中包括在给定速度场中做被动对流
复杂网络已经吸引了来自科学技术不同领域研究者的大量关注。继复杂网络中的小世界特性与无标度性质之后,自相似性已经成为复杂网络的第三大基本特征并且在近年来得到了广泛的研究。分形分析可以有效地揭示一些分形对象的自相似性。然而,重分形分析是一个用来系统性地刻画理论和实验复杂分形对象空间异质性的更加强大而有效的工具。尽管已有一些研究复杂网络自相似性的分形和重分形分析算法,但是这些算法的效率并不高而且已有的算
本文主要研究了阶数为α∈ (0,1)的具有Caputo导数的时间分数阶Navier-Stokes方程,这类方程可以用来模拟分形介质中的反常扩散现象.我们主要讨论它的适定性.本文的工作主要分成五个部分:在第二章中我们首先通过把Helmholtz投影子作用到目标方程上消去压强项,将其转化为抽象形式的发展方程;然后根据分数阶发展方程理论给出其相应的积分方程并定义适度解;利用半群理论和不动理论得到方程的全
本论文主要研究了M-模糊化凸结构理论,(L,M)-模糊凸结构理论,及(L,M)-模糊凸结构(子)范畴理论.在M-模糊化凸结构理论中,首先,我们引入并刻画了M-模糊化(弱)JHC凸空间、M-模糊化几何区间空间、M-模糊化(弱)Peano、Pasch、沙漏区间空间,并讨论了它们之间的关系。其次,我们引入了M-模糊化基点序的概念,并用其刻画了 M-模糊化几何性质、M-模糊化Peano性质、M-模糊化Pa
对于非线性扩散占优Volterra型偏泛函微分方程的时间离散化,以往最常用的方法是Volterra泛函微分方程(VFDE)的隐式Euler方法及梯形方法,有时候也使用隐式中点法及Bellen和Zennaro所建议的二级Lobatto IIIC Runge-Kutta方法。为了避免出现Bellen和Zennaro所指出的“order failure”和“stability failure”现象,本文
分数微分方程(Fractional Differential Equations)在许多学科领域有广泛的应用,这些学科领域的许多数学模型都是用分数微分方程来描述的.近十年来,分数微分方程边值问题得到了迅速的发展.随着这一方向研究的进一步深入,研究内容和研究方法不断得到丰富,非线性分数微分方程边值问题获得了许多研究成果.本文主要研究非线性分数微分方程边值问题解的存在性.本文由四章组成,主要内容如下:
本文研究二维各向异性Delaunay三角化网格的生成与优化算法及其应用.当物体在不同方向上的变化速率不同时,例如树木向上生长与横向生长的速率是不同的,我们称其具有各向异性性质.各向异性问题在各向异性晶体材料、多孔介质、石油工程、粘性流模拟、图像处理等科学工程计算中有着广泛的应用.给出一个各向异性问题,进行有限元数值求解时,合适的网格剖分至关重要,选择符合各向异性特征的三角化网格是自然而然的想法.各