基于反射系数和离散误差的PML参数优化

来源 :浙江大学 | 被引量 : 0次 | 上传用户:socks2010
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
完全匹配层(PML)是由J.P Berenger最先提出的一种有限差分时间域(FDTD)方法的吸收边界条件。理论上,在无穷计算区域加入PML后,任意角度和频率的外行波都能够完全被吸收。该方法首先基于Maxwell模型提出,但是在近年被大量的研究应用于各种波动传播控制方程中,包括Helmholtz方程,并对其进行误差分析和优化等。一般来说,PML方法的误差主要来源于两个方面:截断处产生的反射误差,离散产生的误差。为了减小反射误差,可以选取合适的材料使得吸收函数很大,阻尼很大,但是这种情况下会使得离散误差变大;如果使得离散步长尽量小以减小离散误差,此时不仅计算量是非常难以解决的问题,并且反射误差也会受到影响。由此看出两者呈现此消彼长的趋势,所以必须要寻求两者的平衡。目前由于误差的复杂性以及计算量等问题,大多数的工作都假设离散误差很小,只关注反射产生的误差而很少聚焦离散误差。本文首先推导了PML在数学上的表现是一种复坐标的变换,然后基于Helmholtz方程模型,在有限差分格式下得到反射系数和离散误差的表达式,以赋予相同权重的反射系数和离散误差之和为目标函数,选取合适的吸收函数形式使得目标函数最小,更接近外行波在无穷测度传播区域的行为。  对于PML最重要的指标反射系数,它作为入射角θ的函数,是不确定的值。下文中将采用反射系数模的平均值,即平均反射系数来表现,整体的展现了反射系数的大小,用|R|来表示。对于离散误差,主要考虑FDTD方法产生的离散误差损失,用error来表示,error是每一个离散区间上产生的误差的和。  然后选取目标函数M=|R|+|error|,当M最小时,我们得到的参数作为PML最优设计的参数,最终在数值试验的验证下,得到一些较优的参数选取。
其他文献
十六大以来,在党中央、国务院的坚强领导下,我国反腐倡廉法规制度建设取得新的明显成效。——以党章为核心、以监督条例为主干,以配套规定和其他监督规范为重要补充的党内监
学位
当前,在论述科学发展观时,必然会联系干部政绩观以及干部考核体系。从某种程度上讲,科学发展观能否真正落实,取决于政绩观能否真正确立,以及干部考核体系能否真正得到完善。
加法组合是组合数论的主要内容,也是近年来国际数学界研究的热点之一,主要包括整数Ramsey型问题,和集问题,零和问题等方面.在本篇论文中,我们主要应用代数方法与调和分析方法研究
学位
学位
就连多年从事新闻史研究的许多学者可能也未必清楚,《人民日报》在初创时期曾经出版过两份子报,一份是《人民晚报》,另一份是《人民週报》。可能因为出版时间较短、发行数量
本文主要研究了几类非线性方程的动力学行为,主要讨论了其解的全局稳定性,收敛性,及Lotka-Volterra方程的定性永久生存性.  首先考虑了二次Logic差分方程xn+1=xner(1+xn-ax2n)
可压缩多介质辐射(磁)流体力学问题的数值方法研究是惯性约束聚变、Z箍缩(Z-pinch)、武器物理等领域的重要研究课题.目前,基于网格的拉氏方法和ALE方法是数值求解此类问题的主要
论述了如何有效地回收细粒尾矿,使排往尾矿库的尾矿减少,进而减小或消除尾矿库,减少占地,利于环保。这对我国的矿山企业具有很重要的参考价值。 It discusses how to effectivel