微分求积法、Sinc方法中几个新方法的研究

来源 :同济大学 | 被引量 : 3次 | 上传用户:qingxu007
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
微分求积法(Differential Quadrature Method,英文简称DQM)是由Bellman及其同事(BELLMAN和CASTI 1971,BELLMAN,KASHEF and CASTI 1972)在70年代初期提出的求解非线性偏微分方程的一种新方法。自提出以来,DQM已经被成功地应用到许多工程物理问题中去。这个算法数学原理简单,计算精度高,计算量小,使用方便,不依赖泛函和变分原理,边界条件不必另外考虑。当问题具有全局性光滑解时,DQM因为其精度高的逼近及只需较少的网格点而成为优于有限元和有限差分的更佳选择。在最近的一份比较研究(Malik和Civan 1995)中,我们看到DQM在计算精度及计算量方面要明显优于有限差分及有限元方法。Bert和Malik(1996)给出了把微分求积法应用到几个机械工程问题中的示例。区域分裂法(Domain Decomposition Method,英文简称DDM)是上个世纪六十年代由德国数学家H.Schwarz为解复杂区域上的偏微分方程而提出的。DDM有三种类型:(1)带重叠的区域分裂法(DDM-1),(2)不带重叠的区域分裂法(DDM-2),(3)带虚拟成分的区域分裂法(DDM-3)。DDM有许多优点,如区域分裂方式的任意性、区域分裂后物理问题的数学描述的多样性以及问题求解方法的灵活性,因此它已成为大型科学计算问题的算法设计的重要方法。通过DDM,使乘法的复杂性大大的减少了,同时也使DQM的有效性大大的增加了。DDM发展到今天,它已不再是一个单纯的计算方法甚至一类计算方法,而是问题求解一整过程的思想与方法的集成。但是众所周知,在某些情况下这种方法并不适用。例如,考虑Dirichlet问题,传统的Lagrange插值并不能直接用于在无界区域上的问题,同时它也不适合逼近在空间方向振荡的函数。但是若使用Sinc函数作为基,可以很好地解决这些问题(Stenger 1983,Stenger 1993,Lund和Bowers 1992)。Sinc配点法是用Sinc函数(Stenger 1983,Stenger 1993,Lund和Bowers 1992)作为基函数来构造微分求积法并以此为基础来求解微分方程。多项式基是在Sobolev空间中讨论逼近性态的,而Sinc函数是在解析函数空间中研究逼近性态的。Sinc函数逼近方法在引进指数变换后其逼近真解的误差可达到指数阶收敛。Sinc函数对边界层问题与振荡问题具有很好的逼近效果。所以,Sinc配点法在求解微分方程时将达到真正的高效率与高精度,对微分方程数值解计算方法具有革命性的意义。Sinc数值方法要推广与应用遇到的一个主要瓶颈问题是如何处理有限区域的边界条件特别是边界点上的法向导数。Narasimham,Majdaln和Stenger(2002)以及Morlet,Lybeck和Bowers(1997)等在处理边界点法向导数时改用差分方法处理。这将大大降低精度。并且由于Sinc点在做了坐标变换后是向区间两端密集的,所以简单地用差分代替微商将使有效位数损失很多,误差增大。不是一个好办法。寻找处理计算区间的各种边界条件,特别是如何计算边界上的法向导数是推动Sinc数值方法发展的一个重要且极有意义的工作。Stokes方程和Navier-Stokes方程在力学与工程领域中有着广泛的应用。数值计算Stokes方程和Navier-Stokes方程在力学与工程中已占有越来越重要的地位。1996年美国能源部及其属下的三大国家实验室提出了”ASCI”计划(即加速战略计算创新计划),该计划的实质是将过去以实验为主的科学研究方法转化为以数值模拟为主的科学研究方法。自从这两个方程提出以来,人们一直努力地寻找一种计算简单,实施方便同时精度有比较高的算法。每年在各种计算数学期刊杂志上有成百上千篇的这方面的文献资料,它们分别从不同的角度分析研究这两个方程的数值算法。如何快速精确计算Stokes方程和Navier-Stokes方程,一直是微分方程数值计算领域中的热点问题。这一问题的主要难点在于:压力与速度耦合,没有明确地给出确定压力的方程,计算得到的速度如何确保其满足散度为零的条件等。环绕这些问题不同的解决方案已提出。特别是压力和速度分离的方案被认为是一种准确且有效的方案而备受关注。传统的分离方案要寻找满足divu=0或其弱形式的子空间的一组基底。但要做到这一点是极其困难与复杂的。寻找简单且有效的方法来使速度与压力分离从而高效准确地计算Stokes方程和Navier-Stokes方程将是一件极具意义的工作。本文的主要贡献和创新之处有如下几点。我们提出了一种增量计算的方法,并且将其与微分求积相结合,可以很有效的处理非线性方程;我们提出了一种微分求积区域分裂法,可以有效地处理一些带有奇异性的问题;提出了微分求积消元法,这个方法可以有效地求解Stokes方程以及Navier-Stokes方程;然后我们创新地提出了一种带有边界处理的Sinc collocation方法,可以方便地处理微分方程,而无论它带有什么样的边界条件;最后,我们提出了一种基于最高阶导数逼近的Sinc方法,它克服了数值微分所带来的固有缺点。本文的结构如下:第一章是引言;第二章是对于微分求积区域分裂法以及增量计算法研究;第三章是研究Stokes方程以及Navier-Stokes方程上的这些方法;第四章是对于Sinc方法的研究,并提出一种带有边界处理的新的Sinc配点法;第五章是在第四章的基础上,研究基于高阶导数逼近的Sinc方法;第六章是总结与展望。
其他文献
在科技迅速发展的今天,网络已经成为现代社会主要的信息交流方式,并且覆盖了社会生活的各个领域。而随着网络的渗透,大学生成为网民中的主要群体并且在飞速的增长,网络正以空
我们在实践中引入"生态教学"理念,贯彻"以生为本,以师为主"的教学思想,结合多维度课程资源,调动学生学习英语的主动性和积极性,从教师资源、生活资源、网络资源和研究性学习
本文介绍了日照市“北水南调”工程赔青工作的做法,赔青工作中存在的主要矛盾和问题,及解决方法。
本文结合王铺渡槽工程建设监理实践,介绍了监理工作的实施方法、过程与效果.
提问是课堂教学的重要手段之一,问得好,一方面能灵活地推动课堂教学进程,另一方面不仅能深化学生对知识的理解与掌握,还能培养学生的思维能力,同时也能起到活跃课堂气氛的效
为了寻求局部可分度量空间的映像,引入双点星sn覆盖这一概念,并借用近年来由传统的Ponomarev-system推广而得的ls-Ponomarev-system,给出了局部可分度量空间1-序列覆盖映像的
教育部制订的《普通高中英语课程标准(实验)》对高中英语必修(一)写作有明确的要求:(1)能用恰当的格式写便条和简单的信函;(2)能简单地描述人物或事件,并简单地表达自己的意
放着人人羡慕的百万富翁不当,却把自己价值300万的煤矿捐给了村子,并且在危难之时当起了这个“烂村”的书记。王振昌,一个土生土长的农民,凭着一腔热血和无私的奉献精神,带领
期刊
新课标下的课堂是鲜活的,教师要给学生营造一种民主、和谐、愉悦的学习氛围,让他们"有机可动".要以学生为中心,让学生在课堂上充分动耳、动口、动眼、动手,从而使学生成为主
“要不是参加新型农村合作医疗,我恐怕活不到现在,‘新农合’真是我的救命恩人啊!”近日,61岁的鹿邑县生铁冢乡丁怀村村民刘万善感慨地说。