移动网格方法及其应用

来源 :北京大学 | 被引量 : 12次 | 上传用户:woai6672690
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
在广泛的实际应用问题中,往往出现解的性质相对恶劣,方程在求解区域的局部变化非常剧烈,或者是求解区域整体相对较大,却又要对其中小部分上解的细节信息要求很高的情况。对于这样的问题,在均匀的网格上求解是不现实的,尤其是高维的问题,计算量远远超出硬件的能力。自适应方法是解决这种问题的一个途径。移动网格方法,作为自适应方法的一种,主要是为了解决发展方程的计算问题而设计的方法,已经有二十多年的研究历史。本文主要是阐述了一种基于调和映射的移动网格方法。在移动网格方法中,需要引进一个逻辑区域作为参考,网格的移动往往通过一个区域变换来实现,网格可能发生缠绕的问题是移动网格方法中一个一直没有解决的问题,我们利用调和映射来构造区域间的变换,使得变换的存在唯一性有了理论的保证,这就为避免网格缠绕打下了基础。我们进而引进了一个迭代的过程来实现网格的移动,避免了数值原因导致的网格缠绕,彻底地解决了网格缠绕的问题。我们将网格移动和方程求解完全分开,从而使得移动网格方法在各类不同的问题中的应用,被完全归结为构造控制函数的问题,并且有利于程序开发。我们还设计了一个在不同的网格之间进行插值的格式,用来将解函数在不同的网格间过渡,数值结果表明这个格式和多项式插值方法相比,具有比较好的性质。我们用一系列比较难于处理的典型发展问题进行了计算,都取得了比较好的效果。计算结果说明,这样的方法产生的网格,变化的幅度可以很大,灵活性高,能够处理一些很复杂的情况,鲁棒性强。在网格移动以后,解的误差获得了有效的减少。我们发展的方法可以自然的在三维情况下实现,并且给出了数值算例。三维问题的移动网格方法,在以前是很少看见结果的。边界网格和内部网格偶合移动的问题,是一个实践中要求非常迫切,却又一直没有结果的问题。长期以来,只能对边界当成一维问题进行特殊处理。在二维的情况下,为了能够比较系统的处理边界节点的移动问题,我们扩大区域变换的求解空间,给出了一个将内部节点移动和边界节点移动偶合的格式,在一个扩大的空间里构造一个调和映射,将问题转换为一个边界控制问题,并且有存在唯一性的保证。在一定的情况下,这样的方法还能够在三维实现,并能推广到更一般的情况,将边界特殊处理的方法是不能做这样的推广的。在边界固定的情况下,构造网格移动的方向是比较容易实现的,但是在将边界和内部偶合的情况下,构造出网格移动的方向,我们的这个工作是有创新意义的。和边界固定的方法相比,使用这种偶合模式产生的网格,网格的质量更高,进一步的减小了解的误差,得到的网格的性态,显示出一种内蕴的合理性。我们还将移动网格方法应用于静态问题,并在其中研究使用后验误差估计来构造控制函数的方法。尽管加密的方法几乎总是根据后验误差估计来进行自适应,但对于移动网格方法来说,我们还没有看到使用后验误差估计构造控制函数的文献。我们得到的静态问题的数值结果,说明了用于构造控制函数的后验误差估计必须具有足够的精度才能够构造出能够获得满意的网格的控制函数,这些问题包括变分不等式问题和椭圆的最优控制问题,我们对这两个问题给出了能够用于构造控制函数的后验误差估计。本文后面的内容组织如下:第零章介绍需要用到的数值计算方法的一些基本知识,第一章介绍几种自适应方法的关系和历史,第二章介绍基于调和映射的移动网格方法的基本格式,第三章介绍二维带边界的移动网格方法,第四章介绍变分不等式和最优控制问题的后验误差估计,并使用后验误差估计构造控制函数的技术,第五章是二维发展问题的算例,第六章是变分不等式和最优控制问题的算例,第七章是带边界二维问题的算例,第八章是三维问题的算例,第九章是总结和下一步工作的计划,附录A中是一个关于移动网格方法的误差估计的结果。
其他文献
该文主要从两方面讨论连续Domain理论中的几个问题.一方面在经典Domain理论中,对连续Domain本身的性质进行了深入研究,得到连续Domain的许多类似连续格的性质,给出了连续Doma
在评估软件可靠性之前,必须要进行软件测试.该文采用遗传算法来进行测试算例的自动生成,重点对于两类典型的遗传算法--杂乱遗传算法和遗传K-均值算法进行较为详细的分析,从而
该文首先回顾了基于内容图像检索的理论及其应用研究的整个发展过程,全面综述了在图像检索,尤其是基于内容的图像检索领域的技术和现状,探讨了图像检索系统中的关键技术,并介
Dempster,Laird和Rubin于1977年提出的EM算法是一种用于不完全数据情形下进行最大似然估计或最大后验估计的方法。由于具有很好的收敛性质和广泛的应用背景,EM算法一直受到统
随着信息技术及互联网的发展,为图书情报档案管理工作带去挑战的同时,也带去了新的发展途径,图书情报档案管理工作一体化的进程己经势在必行.本文对图书情报档案信息一体化管
全文共分八章,第一章综述了金融工程的内涵及其研究对象.第二章通过分析股指期货的产生背景,探讨了股指期货产生的原因、难产的原因和各国发展的情况;并分析了股指期货所经受
高层建筑的防火问题已越来越引起政府和社会的关注.目前,建筑防火安全评估通常是根据建筑防火安全规范对构成系统的各因子进行校验和评价.通过对影响建筑防火安全因子的控制,
该文首先运用分形理论的R/S方法来判断个股是否存在明显的主力投资行为.对于存在主力投资行为的个股,对相关时间序列进行对数线性去趋势化处理,以大盘指数为特征,消去宏观面
该文对用于模拟紧离散动力系统的形似KdV方程的Rosenau方程提出了混合有限元格式,在空间上使用时间连续离散,在时间方向上我们考虑使用Euler向后差分离散.并且给出了半离散和
该文主要研究延迟微分方程的并行迭代的Runge-Kutta方法.利用隐式Runge-Kutta方法较好地稳定性,采用"分而治之"的思想,构造具有良好并行性的迭代格式,其目的是提高计算效率.