一类非规则网格上椭圆问题离散与求解研究

来源 :中国科学院软件研究所(北京) 中国科学院软件研究所 | 被引量 : 0次 | 上传用户:tjunu520
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文着眼于椭圆型偏微分方程的数值求解,重点研究了其在一类非规则结构化网格上的有限差分和有限元离散方法以及相应的快速解法。   经典有限元和有限差分方法在处理二维问题时多采用三角形和四边形网格。然而,平面正则剖分的方式除了三角形和四边形外还有六边形,并且与三角形和四边形相比六边形更接近于圆。六边形也广泛存在于自然界以及诸如材料科学、核工程等很多应用领域之中。这些因素启发我们研究六边形网格。本文研究的六边形网格主要有两类:一是对偶六边形网格,可以看作是平面三向三角形网格(由单一三角形剖分构成的网格)的对偶网格;另一类称作平行六边形网格,是由正六边形网格做仿射变换后得到的。   在本文第一部分,我们构造了对偶六边形网格上Laplace算子的成对四点差分格式,证明并且通过数值实验印证了该格式虽然仅有一阶截断误差,但能够达到二阶整体精度。我们证明的主要思想是利用矩阵变换,可以把对偶六边形网格解耦为两套粗三向三角形网格,而三角形网格上的七点差分格式则具有二阶精度。通过调用粗网格上的快速求解器,也可以把这种解耦思想用于求解四点格式的离散线性方程组。   本文第二部分的工作是研究平行六边形上的有限元。我们简要回顾了三种构建六边形元的思想,并重点研究了基于不完全多项式插值的方法。我们选择的不完全多项式空间是三向坐标意义下的三线性和旋转三线性空间,分别对应于基于边和基于顶点的六自由度六边形元。并且我们还提出对插值条件可以再额外增加单元上平均积分,从而获得两种新的七自由度六边形元。对这四种非协调元,我们都给出了先验误差估计,它们在L2模(或能量模)下都能达到二阶(或一阶)精度。这一误差估计也得到了数值实验的验证。   在第三部分,我们把目光放在典型的一阶椭圆方程组——Cauchy-Riemann方程上。受传统笛卡尔交错网格的启发,我们提出一类新的非规则交错网格,它由任意三向三角形网格及其对偶六边形网格组成。我们构造了这种网格上的一类三色差分格式,并且证明了这种格式与四点格式类似,都有二阶整体精度,虽然截断误差仅一阶。数值实验验证了差分格式的精度,并且检验了我们利用网格解耦思想提出的快速解法的效果。  
其他文献
学位
语义Web服务是利用语义知识对Web服务进行标注,从而实现机器和人都可理解;而语义Web服务组装则被视为领域应用快速有效开发的重要途径,也是国际上的研究热点。 首先,本文对当
本文首先介绍了Web Services的概念,对比分析和研究了WebServices技术在J2EE与.NET平台下的不同发展及应用,并详细介绍了.NET平台下Web Services的一些关键技术以及重要特点。深
从80年代开始税收电子化到1994年开始实施的金税工程,我国税收信息化逐步从如何管理数据发展到如何利用数据的层次。随着税收信息化的推进,各税务部门已建立了一系列税收征管系
数据仓库是一种有效的数据存储工具,OLAP是面向数据仓库的数据分析和知识提取的有效工具。本文以国家重大科学工程为背景,对LAMOST的天体光谱数据仓库系统进行了研究,对提高
目前上许多应用程序需要创建和管理网络会话。由于终端的移动性,多个命名,或者同时进行多个多媒体会话,所有这些使得应用程序变得很复杂。SIP(Session Initiation Protocol)协议
随着计算机技术与通信技术的迅速发展,计算机网络得到了广泛的应用,越来越多的政府机关、军事部门和商业机构通过计算机网络存储和处理有价值的、敏感的数据。以计算机网络为核
随着.NET平台的不断完善和多层体系结构的不断发展,基于.NET分布式多层体系结构在企业的开发中应用越来越广泛。三层C/S体系结构,克服传统二层C/S体系结构的缺陷,将表示层和业务逻
作为学术界和工业界关注的热点之一,数字水印被认为是一项在数字版权保护方面非常具有潜力的技术。人们对水印技术在标识所有权、跟踪使用、保证合法授权访问、阻止非法复制和
软件构件库是软件复用中管理可复用构件的基础设施,它为基于构件的软件开发中的构件生产和构件复用两方面工作之间建立了桥梁。构件库中的构件是来自软件生命周期各个阶段的可