四阶椭圆型方程若干有限元新方法和高效求解算法

来源 :上海交通大学 | 被引量 : 4次 | 上传用户:XIAOZHOU914
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
四阶椭圆型微分方程广泛应用于固体力学、材料科学和图像处理诸领域,因此对它的数值解研究不但具有重要的理论意义也具有直接应用价值.本文的主要工作是构造了重调和方程的一种基于Poisson求解子的快速求解器;对四阶椭圆偏微分方程的MWX元方法设计了两网格局部和并行求解算法并进行了误差分析;将C0间断有限元方法应用于四阶椭圆偏微分方程并给出求解由该方法离散得到的线性代数方程组的区域分解法;也提出了自适应有限元方法的一个抽象框架.首先,构造了重调和方程的一种基于Poisson求解子的快速求解器.通过位能极小原理,建立了重调和方程和Stokes方程之间的等价性,即重调和方程等价于一个Stokes方程和两个Possion方程.类似地,在有限元离散情形,建立了重调和方程Morley元方法和Stokes方程非协调P1-P0元方法之间的等价性,即重调和方程Morley元方法等价于一个Stokes方程非协调P1-P0元方法和两个Possion方程Morley元方法.然后基于这个等价性,结合代数多重网格法,构造了重调和方程的一种基于Poisson求解子的快速求解器.最后,用数值试验验证了这个求解器的有效性.其次,提出了三种用于求解任意空间维数的四阶椭圆偏微分方程的两重网格局部和并行的Morley-Wang-Xu (MWX)元方法.由于MWX元空间不是嵌套的,需要引进一些网格转移算子,以从全局的粗网格有限元解和局部的细网格修正来获得改进的全局有限元解.首先通过构造基于修正Argyris元的网格转移算子,提出了求解MWX元离散的四阶问题的三种局部和并行算法.然后,证明了相应的数值解的离散能量误差的阶为O(h + H2),其中H和h分别为有限元单形剖分的粗细网格的网格大小.同时,构造了基于算术平均网格转移算子的局部和并行算法,证明了相应的数值解的离散能量误差的阶为O(h + H2(H/h)(d?1)/2),其中d为求解区域的空间维数.进一步,我们提供了数值实验来说明这些算法的计算效果.再次,提出了求解薄板弯曲问题的一类新的C0间断有限元方法,并给出了数值解在离散能量误差和H1范数下的误差估计.首先将四阶偏微分方程写成二阶系统的形式,获得了构造求解原四阶问题的C0间断有限元方法的框架.然后,建立一个离散的稳定性等式,并由此选择可行的数值迹,得到求解薄板弯曲问题的一类稳定的C0间断有限元方法.根据二阶椭圆问题间断有限元方法误差分析中的思想和推导技巧,我们得到了薄板弯曲问题LCDG方法数值解在离散能量范数和H1范数下的最优阶误差估计以及CDG方法数值解在离散能量范数和H1范数下的误差估计.最后,用数值实验验证了理论的收敛阶.接着,构造了薄板弯曲问题的局部C0间断有限元方法的两水平加型Schwarz预条件子并估计了条件数.对于我们提出的求解薄板弯曲问题的LCDG方法,通过定义特殊的网格转移算子,给出了该方法的区域分解算法.进一步估计了经过预处理的线性代数方程组的系数矩阵的条件数,该上界估计在δ≈H时是最优的.在小部分重叠情形,即δ<< H时,利用更多的子区域的信息,改进了条件数的上界估计.然后,提出了基于标准循环的自适应有限元方法的抽象框架并得到了拟最优收敛性和最优复杂度.首先提出一些假设,基于这些假设得到自适应有限元方法的拟最优收敛性.在作了进一步的假设后,通过引进总误差,证明了自适应有限元方法的最优复杂度.最后通过前面提出的假设,将这个抽象框架应用于各种问题,包括一般二阶椭圆偏微分方程的高阶有限元方法, 2m阶椭圆偏微分方程的Morley型有限元方法,不定时谐Maxwell方程和H(div)方程.最后,对于薄板弯曲问题k = 0,1时的混合有限元方法,验证了自适应有限元方法抽象框架中提出的假设.利用Helmholtz分解,对所有的k≥0,构造了薄板弯曲问题的混合有限元方法(H-H-J方法)基于残差的后验误差估计子并证明了估计子的可靠性.利用泡函数技巧分析了误差估计子的有效性.而且根据力在内边上的跳跃可以被体积部分所控制的性质,改进了误差估计子.进一步,获得了离散的Helmoltz分解和离散的inf-sup条件,利用这两个工具和后验误差分析中得到的性质,建立了拟正交性和离散可靠性.对于其它的假设,在k = 0,1时容易得到.
其他文献
有机多孔聚合物制备方法的发展,使得科研工作者们能够在纳米尺度上构建孔结构。通过单体或者聚合物作为起始原料的模板法或利用嵌段共聚物自组装法通常可以制备具有特殊的规
语义分割是对图像中的每个像素点进行分类,将图像分割成若干个具有特定语义类别的区域,是计算机视觉领域的核心技术。近年来,卷积神经网络(CNN)不断取得突破性进展,利用深度
当前天文数据处理技术已经进入了数据密集型的天文信息学时代,大数据是比较典型的特征。在太阳观测中,具体表现为数据量庞大、数据采集速率高和数据连续性增长。传统的本地存储
脉宽调制(Pulse Width Modulation,PWM)器具有自然采样特性、边带效应和混叠效应等非线性特性,对DC-DC变换器的稳定性具有重要影响。而由PWM DC-DC变换器构成的级联系统,不仅
对于道安法师提出的"五失本三不易"之说中"失本"的内涵,学界一直存在分歧。代表性观点有二:一种认为"失本"即"改变原文形态"(形变),一种认为是"违背原文内容/精神"(义变)。本
光疗主要包括光动力疗法和光热疗法,其作为一类非侵入性、可实时操作的肿瘤治疗手段引起了科研工作者们广泛的研究兴趣。高效光疗试剂的构建是实现肿瘤有效治疗的关键环节。
目的探究糖尿病专科护士护理服务模式和应用效果。方法选择2014年12月—2015年12月来该院就诊的200例门诊住院2型糖尿病者为研究对象,确定糖尿病病专科护士,全面开展患者会员
本文文章分析了现代艺术和美学的命名,考察了艺术自主性的理论谱系,尤其是古典的自主论与现代自主论的差异,最后具体探讨了导致艺术自主性出现的主要社会文化动因。
随着数码类电子产品的广泛普及,清晰的高质量的图像在众多应用领域起着举足轻重的作用。然而由于成像系统自身的限制及成像过程中众多因素的影响,得到的图像普遍存在模糊和噪
大气重力波是大气中的基本波动形式之一,在中高层大气动力学过程中起着十分重要的作用。低层大气中激发的重力波携带着动量和能量向上传播。重力波在向上传播的过程中,随着高