一维线性Helmholtz方程的紧致差分数值模拟

来源 :山东师范大学 | 被引量 : 0次 | 上传用户:qyyqyy202
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文考虑如下的一维Helmholtz方程(?),其中,k为波数,f为源函数,u为待求函数表示波场的压力。Helmholtz方程主要描述的是一类波传播现象,包括电磁波、声波、光辐射等,在工程实际和科学技术中有很重要应用.其数值求解是微分方程数值解法研究领域的热点问题之一.对于有界区域上的Helmholtz方程,高阶紧致方法大都采用的是紧致有限差分法.此方法缺点是在处理非一致网格时会降低格式的精度,并且不具有守恒性.因此,在本文第一章中我们采用紧致有限体积方法研究有界区域上Helmholtz方程相关的边值问题.第一章第二节中介绍了Helmholtz方程的由来及物理背景;第三、四节中分别提出了一维Helmholtz方程基于Dirichlet和周期边值问题的四阶紧致有限体积方法;第五节的两个数值实验表明给出的两种格式均是四阶格式,并且周期边值问题的格式对于大波数问题同样适用.对于无界区域上的Helmholtz方程,区域的无界性给此类问题的数值求解带来很大的困难.目前,完美匹配层(PML)是解决此困难的有效方法之一.引入完美匹配层使得波穿过内边界处不产生反射,在层中有阻尼,使得波进入吸收层后迅速衰减,并在达到吸收层外边界时衰减到零,以便在计算中可以使用简单的齐次边条件.从而,在有限计算区域上带人工边界条件的边值问题是原问题的很好地近似.基于这些考虑,本文第二章在无界区域上构造完美匹配层(PML),使大部分外行的波都能被吸收.第二章第二节中,我们采用复坐标延伸法,得到了Helmholtz方程的PML-Helmholtz方程.由此,无界区域上的Helmholtz方程就简化为有界计算区域上的边值问题.有限体积法和有限差分法是数值求解PML-Helmholtz方程边值问题的有效方法之一.本文在第二章第三节中对此问题提出了有限体积格式,得到了较好的数值结果,但精度不高.而有限差分方法往往产生严重的数值频散,高阶显式有限差分格式需要用较多的网格点,同时也不利于边界的处理.为此,我们分别采用紧致有限体积方法和紧致有限差分方法来处理PML-Helmholtz方程边值问题,使得两种格式在内部区域和人工边界上均能达到四阶精度.本文第二章第四节提出了紧致有限体积格式,第三章第二节给出了PML-Helmholtz方程边值问题的紧致有限差分格式,通过差分解的先验误差估计,证明了格式的稳定性.在这两章的最后,我们均进行了数值算例,证明了两种格式的有效性。
其他文献
目的探讨分析小儿急性中毒的常见病因,总结疾病的临床特点以及有效的治疗效果。方法方便选取该院2014年1月—2017年1月收治的52例急性中毒患儿作为研究对象,回顾性分析所有患
国家质检总局、国家标准委日前发布了《房间空气调节器能效限定值及能效等级》强制性国家标准。该标准将从6月1日起在全国范围内实施。这一标准规定了房间空调器产品新的能效
采用统计量的标准化方法,研究了均值的单侧限计量抽样检验的理论与方法,在方差σ^2已知的条件下,前人已给出严格理论,并依据该理论制定了GB8054表1;在方差未知的情况下,制定GB8054表2的依据是不严格
10月13日下午,江苏省第七届绿色建筑国际论坛暨江苏国际绿色建筑展览会组委会在南京召开新闻发布会,宣布该论坛与展览会将于11月5日至7日在南京国际展览中心举办。相关政府部
2009年2月3目,记者从有关部门了解到,去年,沈阳市结合“优化结构年”活动,把节能工作作为调整经济结构、转变发展方式的重要抓手和突破口,取得了积极进展。
旨在关注太阳能热水器“普通百姓式”的发展道路,着眼于广为普通居民所接受的直插式太阳能热水器,分析了其特性、优点与一般规格下的安装要求;并针对其在推广过程中遇到的问
近日,全国建设领域节能减排检查组对福建省建设领域节能减排工作进行检查。被查21个项目均按节能标准规范设计、施工和验收,节能措施合理,很多项目达到较高水平。检查组对福建的