各向异性三角化网格生成与优化研究及应用

来源 :湘潭大学 | 被引量 : 0次 | 上传用户:ln86119
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文研究二维各向异性Delaunay三角化网格的生成与优化算法及其应用.当物体在不同方向上的变化速率不同时,例如树木向上生长与横向生长的速率是不同的,我们称其具有各向异性性质.各向异性问题在各向异性晶体材料、多孔介质、石油工程、粘性流模拟、图像处理等科学工程计算中有着广泛的应用.给出一个各向异性问题,进行有限元数值求解时,合适的网格剖分至关重要,选择符合各向异性特征的三角化网格是自然而然的想法.各向异性三角化网格由对称正定度量来反映各向异性特征.已知各向异性问题,如何定义合适的度量并自动生成高质量的各向异性三角化网格是本文的主要研究内容.本文首先提出了常度量下各向异性三角化网格生成与优化算法,该算法主要通过波前法由区域边界向内部逐次插入网格结点,同时基于Anisotropic Centroidal Voronoi Tessellation(ACVT)进行波前局部优化和网格整体优化,由此生成的常度量下各向异性三角化网格不仅可以很好地符合区域边界,而且其三角形单元在该度量意义下几乎为正三角形.我们将此网格应用于求解常系数各向异性椭圆偏微分方程(PDEs).通过大量的数值实验,找到了该类方程的匹配度量为系数矩阵的逆,基于匹配度量的各向异性网格上求解,我们发现方程离散系统的刚度矩阵条件数远小于其它非匹配网格,所得的数值解具有高精度,并且在网格结点上有超收敛现象,收敛阶为O(h2+α),α≈0.5,而在其它非匹配网格上不具有超收敛特性.因此,常系数各向异性椭圆偏微分方程的合适度量恰为系数矩阵的逆.对于度量随点变化而变化的情况,本文提出了新的各向异性变度量网格生成与优化算法,该算法通过在具有各向异性特征的背景网格上,进行结点的插入与删除,同时基于力平衡原理移动网格结点,生成了符合目标尺寸的高质量各向异性Delaunay网格.然后我们将变度量各向异性网格应用于变系数各向异性椭圆偏微分方程有限元求解中,将常系数时所得到的结论,即系数矩阵的逆为匹配度量,推广到变系数情况.通过数值实验发现,在变系数情况下,匹配的变度量各向异性网格上,同样具有好的离散系统、高精度的数值解以及结点上l2误差的超收敛现象.因此我们的结论适用于一般的各向异性椭圆偏微分方程:当系数矩阵为各向异性,方程右端项为各向同性时,选择系数矩阵的逆作为度量所生成的各向异性网格是方程相匹配的网格.另外我们根据显式多项式恢复(EPR)技术,对有限元解进行恢复,通过数值实验发现,EPR对匹配网格的作用明显优于其它网格.并基于EPR后验误差指示子进行各向异性网格自适应,数值结果证实了该指示子的有效性.最后我们探索了一个实际应用问题――具有强各向异性材料的电池热传导问题,通过在各向异性匹配网格上有限元求解,并与一致网格上的数值解进行分析比较,发现各向异性网格优于一致网格.
其他文献
铁电材料是一种功能材料,其具有的优异电学和光学性能孕育出了它广阔的应用前景。电畴翻转是铁电材料显示宏观非线性本构行为的微观物理机制,铁电畴的取向直接决定着铁电材料的物理性质和应用方向,而电畴翻转后不同取向铁电畴体积分数的变化则直接与铁电材料及其器件的效率和稳定性挂钩。因此,确定铁电材料的畴取向及其体积分数(分布特性)对铁电器件工程至关重要。多铁性材料作为一种特殊的铁电材料,除具有铁电性外,还兼具有
软件定义网络(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)在许多学科领域有广泛的应用,这些学科领域的许多数学模型都是用分数微分方程来描述的.近十年来,分数微分方程边值问题得到了迅速的发展.随着这一方向研究的进一步深入,研究内容和研究方法不断得到丰富,非线性分数微分方程边值问题获得了许多研究成果.本文主要研究非线性分数微分方程边值问题解的存在性.本文由四章组成,主要内容如下: