基于EEP法的三维有限元超收敛计算与自适应分析

来源 :清华大学 | 被引量 : 2次 | 上传用户:huangguohao123
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
有限元法(FEM)超收敛计算的重要性体现在两个层面。其一,超收敛计算可以在相对稀疏的有限元网格上面获得较高精度的解答;其二,超收敛解可以在有限元自适应分析当中用于构造后验误差的估计量,此即本文主要的研究目标。单元能量投影(EEP)法是有限元超收敛计算的有效方法,已经在许多一维和二维问题中取得成功,但在尝试处理三维问题的时候遇到严重的阻碍。本文重新研究EEP法处理多维问题的理论和算法,实现了三维问题的超收敛计算,并取得了三维有限元自适应分析的初步成功。本文主要工作包括:(1)提出了广义一维问题理论框架。将多维问题当做广义一维问题,参照已经成熟的一维有限元EEP法建立各维度统一的公式,把求解n维问题等效地看成依次求解n步广义一维问题。(2)基于对二维问题各类超收敛计算方案的再研究,提出了适用于广义一维问题的EEP超收敛算法。仔细审查现有方案中输入的一维问题解答对输出的二维问题的影响,总结出保证其超收敛性和连续性的基本要求。基于以上要求,系统性地提出二维有限元的超收敛算法。同时,还计算了二维问题的二阶导数,为进一步的三维问题研究铺平了道路。(3)实现了三维有限元在非规则网格上的EEP超收敛计算。参考二维问题,总结出基于EEP法的三维有限元超收敛算法,可以处理Poisson方程和弹性力学问题等,适用于非规则网格和各类边界条件,突破了本文之前只适用于规则网格的局限。首先在三维六面体单元构成的非规则网格下,求出Poisson方程的超收敛位移和各向导数。然后应用于三维弹性力学问题,推导与三维有限元等效的逐维离散公式、逐维恢复EEP公式以及程序实施方案,求出超收敛位移和各向导数。(4)实现了本文典型问题的三维有限元自适应分析。采用超收敛位移代替精确解进行误差估计,利用六面体单元的棱边均差法进行单元划分和网格加密,直至求解域内有限元解的误差最大模满足用户给定的误差限,最终输出自适应的有限元网格和解答。经有精确解的算例检验,本文算法均能逐点满足误差限,成功实现自适应目标。本文算法适用于三维Poisson方程和弹性力学问题,在工程计算中具备较广泛的应用价值。
其他文献
我们引入(a,η)型区域的概念,对区域边界的凸性作了更细致的分类。通过选取适当的辅助函数,得到了解的先验估计,从而证明了有界区域上一类退化或奇异Monge-Ampere型方程Dirichlet问题解的存在唯一性,并得到解相应的边界H¨older模估计。从而,我们发现了解的边界正则性和区域的凸性之间的对应关系。我们所研究方程的一种特殊情形,例如双曲仿射球方程,和一些重要的几何问题相关。在这些情况下,
这篇论文是基于几篇我和导师丘成栋教授的合作文章。主要目的在于研究孤立奇点的导子李代数.孤立奇点的导子李代数的定义如下:令V是为原点附近的一个孤立奇点,它由解析函数f:(Cn,0)→(C,0)定义.L(V)定义为模代数A(V)的导子李代数.它是一个有限维可解代数,在研究奇点中起到重要作用.L(V)被称作丘代数,而λ(V)表示L(V)的维数,被称为丘数.有一类新的k-阶丘代数Lk(V),定义为模代数A
黎曼流形的分类问题一直是微分几何中一类重要的问题。本文给出了一些特殊的黎曼流形上的刚性定理。主要内容如下:1.令(Mn,g)(n≥ 4)为n维紧致局部共形平坦黎曼流形,有常数量曲率和常Ricci曲率张量的平方和。运用活动标架法,我们证明了 Ricci曲率张量有三个不同特征值的黎曼流形是不存在的。2.我们证明了一个n维(n≥ 4)紧致Bach平坦流形若它的数量曲率为正且σ2是正常数,且Weyl张量满
本文研究了蒙日-安培型方程的一些性质,包含一类蒙日-安培型方程解的径向对称性,以及一类以蒙日-安培型方程为特例的非线性奇异椭圆方程解的边界H(?)lder估计。我们先对一类来自于一些几何问题的蒙日-安培型方程解的对称性进行了讨论,在适当的结构性假设条件下,使用一种新的变换分析了解在无穷远处的渐近行为。进而结合移动平面法,证明了方程凸解的径向对称性。其次,我们研究了一类包含蒙日-安培方程、K-海森方
外尔半金属是一种新奇的拓扑物态,其低能激发与高能物理中外尔费米子遵循相同的规律。由于凝聚态系统更为多样的结构对称性,和丰富的相互作用,在一类正交相的过渡族金属二硫化物系统中还存在违背Lorentz不变量的第二类外尔费米子,并且这种新奇粒子没有标准模型粒子与之对应。尽管过去几十年凝聚态物理学家对过渡族金属二硫化物体系中的谷电子学、能隙可调半导体、电荷密度波以及超导的研究取得了巨大的进展,然而实验上对
实数或复数的超越性是数论的基本问题之一。虽然我们知道几乎所有的实数或复数都是超越数,但要判断一个给定的实数或复数是否为超越数则通常极为困难。现代数论给我们的启示是:同样的问题放在有理数域或正特征函数域上时,在处理技巧上会呈现出许多共性和差异。本文从函数域的角度出发来研究形式幂级数的线性相关性、超越性以及代数独立性,主要包括以下四个方面的内容:一.线性无关性判别准则:我们在正特征函数域上给出了判断形
强相互作用下的量子多体系统可以在低能下演生出丰富的强关联物相与相变现象。在传统理论中,物相与相变由对称性来刻画,然而有一大类新发现的物相,其根本描述在于拓扑而非对称。在本文中,我将围绕高温超导与量子磁性聚焦到几个典型的零温量子强关联多体系统来探讨一些新奇的拓扑物相与相变。受到铜基高温超导实验的启发,我们在描述铜氧面低能物理的t-J模型中引入反铁磁外场耦合,研究超导配对对称性与能隙所受到的影响。在电
Landau-Ginzburg模型一直以来同时受到数学家们与物理学家们的双重关注。围绕Landau-Ginzburg模型的数学研究,将奇点理论与非交换几何、Hodge理论、形变理论和量子上同调等多个不同的数学理论紧密关联起来,并提供了诸多重要的研究课题。其中,Landau-Ginzburg模型之间的镜像对称问题是相关的林林总总的研究方向中最为重要也最有丰富的课题之一。但围绕这个课题的相关研究远未充
颗粒撞击液面是自然界和工业过程的常见现象,也是流体力学和颗粒动力学等学科的基础问题。本文对微米级颗粒撞击液面过程的颗粒和流体运动行为开展研究,揭示颗粒撞击液面的动力学和能量转化机制,为相关自然现象的理解和工业技术的开发提供理论支撑。首先,数值模拟研究了球形颗粒零速接触液面后的运动行为和漂浮条件。颗粒零速接触液面后的运动由邦德数、接触角和密度比控制。给出了颗粒撞击液面后能够漂浮的极限密度比,与实验结
近年来,由于在数据存储系统、通信系统和消费电子产品等方面的应用,具有很少重量的线性码,被专家学者们广泛地研究。文献[1]提出了线性码的一般构造,即由定义集来构造线性码。通过选择合适的定义集,可以生成许多已知的线性码。基于这种构造,目前已经构造出了许多线性码。在本文中,对奇素数p和正整数m≥2,我们通过选取定义集D={x∈F*pm:Tr(x2+x)=0}、D0={x∈Fpm:Tr(x2+x)∈C0(