一种改进的二阶弹性波动方程的最佳匹配层吸收边界条件

在分析前人有关二阶波动方程的最佳匹配层(PML)吸收边界条件的构建方式的基础上,讨论了这些PML吸收边界条件在计算效率和计算精度方面的不足,并给出了一种新的改进PML吸收边界条件.本文的核心思想是在频率域引入一个中间变量,避免了褶积运算,对该变量做傅里叶逆变换得到时间域的精确解.通过模型试算,把本文的改进算法与前人的算法的计算精度和计算时间进行比较.对比分析表明,本文提出的PML吸收边界条件计算量小且精度较高,并且是一种稳定的算法.

发展方程的完美匹配边界层模型

在有界区域上使用数值方法模拟无界区域上的偏微分方程时,要求在有界区域边界上添加人工边界条件。构造的人工边界条件必须能消除边界上不必要的反射,完美匹配边界层是非常有效的人工吸收边界条件之一。本研究采用复值坐标拉伸的方法,在不分离变量的前提下对发展方程推导出了完美匹配边界层模型。通过平面波分析可知,波在初始求解区域边界上无反射,且在有界的完美匹配边界层上振幅是指数衰减的。

完全匹配层吸收边界条件的研究

边界反射是数值模拟中的一个关键问题。构建完全匹配层(PML)吸收边界条件是消除边界反射的理想方法之一。此方法是在研究区域的边界上加入吸收层,使边界上传入吸收层的波随传播距离按指数规律衰减,不产生任何反射,达到消除边界反射的目的。首先给出了二维声波方程完全匹配层的公式,然后通过地面地震模型和井间地震模型数值算例与常规吸收边界条件进行了比较。结果表明,PML边界条件具有吸收更干净、相位对应关系好且能够吸收各种角度的边界反射等优点,从而验证了完全匹配层吸收效果的优越性。

完全匹配层吸收边界条件应用研究

在地震波正演和波动方程偏移的数值计算中,人工边界的引入限制了实际计算区域的大小,边界处反射波的能量严重影响数值计算结果.文中研究了完全匹配层(PML)吸收边界条件在地震波正演及逆时偏移中的有效应用,将PML边界条件的衰减因子引入到声波方程中,利用有限差分的数值方法求解声波方程,在保证差分近似稳定性的情况下,也使边界处的人工反射得到了最小化的处理.数值结果表明,文中所述的完全匹配层方法对地震波的吸收能力不受入射角度的限制,在地震波传播的数值模拟以及叠前逆时偏移算法测试中,均能有效地吸收地震波,是一种高效的边界吸收算法.

完全匹配层(PML)吸收边界条件的理论分析

本文以屏蔽微带线为例,应用谱域方法推导了时域差分(FDTD)法中所用的完全匹配层(PML)吸收边界条件的公式,并且严格证明了有关参数的选择原则.计算表明,在屏蔽微带线中,这种吸收边界条件的反射系数可以小于-80dB.本文所给出的推导方法对于平面分层介质结构是普遍适用的,因为谱域中的每一个波谱实际上都代表一个平面波.文中所给出的推导揭示了这种吸收边界条件的机理,它对于正确使用以及进一步改进这种吸收边界条件都是有意义的.

三维新型非分裂场完全匹配层吸收边界条件

研究了一种3D新型时域有限差分法(FDTD)非分裂场近似完全匹配层(NPML)吸收边界。在直角坐标系中,通过直接对麦克斯韦方程组中的空间偏导的变量进行空间坐标拉伸,得到三维非分裂场的NPML吸收边界,这种吸收边界相对于PML吸收边界编程实现简单,编程复杂度低。通过对三维算例的电磁散射分析,验证了NPML吸收边界的正确性。计算结果表明了该方法是一个完美的吸收边界,可以用于实际目标的电磁散射问题分析。

曲线坐标系下的完全匹配层吸收边界条件

在地震波数值模拟中,需要采用吸收边界条件以吸收人为边界反射。本文针对曲线坐标系下的二阶弹性波方程提出了一种完全匹配层(PML)吸收边界条件。与直角坐标系下的PML吸收边界条件类似,曲线坐标系下的PML吸收边界条件是一种在频率域中给出的人工边界条件,由相应的复坐标变换得到。在变换到时间域后,完全匹配层中将出现复杂的卷积运算。为了避免这些卷积运算,引入了4个中间变量。为了简化自由边界条件,采用正交贴体网格对起伏地表模型进行网格剖分。数值算例表明,该方法可以有效消除人为边界反射。

关于完全匹配层吸收边界条件中差分格式的讨论

对完全匹配层 (PML)中Maxwell方程差分格式的推导 ,本文认为文献 [1]求解微分方程时有不妥之处 ,给出了改正。根据PML介质中电磁波传播特性 ,我们还用直接方法导出差分格式的不同形式 ,在证明两者一致的同时指出后者更具精确性 ,并易于编程实现。

声波完全匹配层吸收边界条件的改进算法

传统位移形式的声波完全匹配层(PML)吸收边界条件以三分裂位移参量来构建,需要求解时间3阶导数,计算量大,计算时间长。为此,提出了一种改进的声波分裂PML吸收边界条件。阐述了传统和改进的PML吸收边界条件的构建原理,通过数值模拟讨论了改进和传统的PML吸收边界以及低阶Higdon吸收边界对边界反射的吸收效果。结果表明,低阶Higdon吸收边界条件吸收效果较差,在波场快照和模拟记录中存在较强的边界反射;改进的PML吸收边界条件与传统的PML吸收边界条件效果相当,均能有效地吸收衰减任意角度的边界反射波。改进算法用褶积近似运算代替时间3阶导数的求解,因此简化了计算过程,减少了计算量,是一种高效稳健的算法。

完全匹配层吸收边界条件研究

为了提高人为截断边界的吸收效率,提出了一种完全匹配层边界吸收衰减因子——正、余弦型吸收衰减因子。采用声波波动方程,并应用统一格式的高阶交错网格有限差分法来构建正演递推算子实现声波全波场正演计算(边界处差分阶数依次递减)。以单层均匀介质模型为例,研究不同吸收层厚度与介质速度和声波模拟主频的关系。数值试算结果表明,不同吸收层厚度所消耗的时间按一定的函数关系变化,边界吸收层厚度与速度和最大频率之间无明确的对应关系,边界吸收层厚度所占节点数不小于6个时就可以得到较好的边界吸收效果。以Marmousi模型为例研究了完全匹配层边界吸收衰减因子在复杂地质模型中的边界吸收效果,与解析解的对比可知,该边界吸收衰减因子能在复杂地质模型中取得满意的边界吸收效果。

有耗介质空间完全匹配层吸收边界条件及其应用

针对Gedney提出的完全匹配层(PML)无法用于有耗各向同性计算域的截断及其对倏逝波的衰减不理想等问题,提出了一种扩展方法。扩展的PML的主要思想是在各向异性的PML中引入与有耗介质空间相一致的复介电常数和复磁导率,使之可以与有耗介质计算域相匹配。通过给PML的张量介电常数、张量磁导率增加衰减因子以加速倏逝波的衰减。构造了PML吸收效果验证模型,数值结果证明了扩展的PML在处理有耗介质计算域截断问题中的有效性。利用该吸收边界条件,采用时域有限差分法计算了电磁脉冲作用下地面铺设电缆的电磁脉冲响应,计算结果和试验结果取得了较好的一致。

广义正交坐标系FDTD算法的准完全匹配层吸收边界条件

本文推导出了广义正交坐标系下的各向异性准完全匹配吸收层 (QPML)边界条件 ,给出了直角坐标、圆柱坐标、保角变换柱形坐标的各向异性准完全匹配吸收媒质 .

完全匹配层吸收边界条件在弹性波波场分离数值模拟中的应用

吸收边界条件是用有限区域问题代替求解无限区域问题来研究地震波传播规律的一种有效手段。为此,本文将完全匹配层吸收边界条件应用于求解等效的弹性波动方程以精确地实现波场分离数值模拟,推导了相应的统一格式的高阶交错网格有限差分计算格式,并对均匀各向同性介质模型和复杂模型进行数值试验,比较本文方法与非分离弹性波数值模拟方法在相同吸收厚度情况下的边界吸收能力,数值计算结果表明,本文方法和分离方法均取得了较理想的边界吸收效果,同时也成功合成了混合多分量波场和完全分离的纯纵横波波场,从而可以用来研究弹性波场的传播规律。

吸收边界条件和完美匹配层

证明了一个Maxwell方程的完美匹配层的的充要条件.条件的核心是吸收边界条件.这个判据可以用来检验文献中的各个模型,也可用来设计一些新的模型.最后给出了一些例子.

指数差分离散的完全匹配层吸收边界条件在ADI-FDTD中的应用

将指数差分离散的完全匹配层(PML)吸收边界条件应用于交替方向隐式的时域有限差分法(ADI-FDTD)中,可以更好地反映PML中电磁波衰减快的特点。数值实验表明,采用这种计算方法在ADI-FDTD中有很好的稳定性和吸收效果。

非线性欧拉方程组的完美匹配层吸收边界条件(英文)

对于非线性Euler方程,提出一类基于完美匹配层(PML)技术的吸收边界条件。首先对线性化的Euler方程设计出PML公式,然后将线性化Euler方程中的通量函数替换成相对应的非线性通量函数,得到非线性的PML方程。考虑到PML方程中包含有一个刚性的源项,文中采用一种隐显Runge-Kutta方法来求解空间半离散后得到的ODE系统。数值实验表明设计的非线性PML吸收边界条件优于传统的特征边界条件。

三维间断伽辽金玻尔兹曼方法的完全匹配层吸收边界条件研究

为解决计算气动声学数值模拟中远场边界反射问题,在三维间断伽辽金有限元玻尔兹曼方法框架下应用完全匹配层算法构造了吸收边界条件,并采用了两种不同的离散求解公式。利用三维高斯脉动源数值案例测试了两种公式的有效性和稳定性,发现仅有一种公式既能有效衰减反射波又具有良好的稳定性;接着研究了影响完全匹配层算法无反射性能的若干因素,其中无反射性能采用归一化误差范数来衡量。研究表明:衰减因子总是存在最优值,且最优值仅与吸收层的参数有关,而与高斯脉动源自身参数无关,故所定义的归一化误差范数具有一定的普适性;当衰减因子采用幂函数分布律时,幂指数取2时无反射性能更好,取4时无反射性能随衰减因子变化更平缓;吸收层厚度和区域半径越大,无反射性能越好,但相应的计算量也显著提高。研究工作可为构造更加实用有效的无反射边界条件提供参考。

完全匹配层吸收边界在孔隙介质弹性波模拟中的应用

模拟弹性波在孔隙介质中传播,需要稳定有效的吸收边界来消除或尽可能的减小由人工边界引起的虚假反射.本文在前人工作基础上,首次建立了弹性孔隙介质情况下完全匹配层吸收边界的高阶速度-应力交错网格有限差分算法,并详细讨论了完全匹配层的构建及其有限差分算法实现.首先,本文通过均匀孔隙模型的数值解与解析解的对比,验证所提出的数值方法的正确性;然后,本文考察了完全匹配层对不同入射角度入射波和自由表面上的瑞利波的吸收性能,将完全匹配层与廖氏和阻尼吸收边界进行了对比,研究了这三种吸收边界在不同吸收厚度情况下对弹性波吸收能力.数值结果表明,在孔隙介质中,完全匹配层作为吸收边界能十分有效地吸收衰减外行波,无论对体波还是面波,是一种高效边界吸收算法.

声波方程完全匹配层吸收边界

本文针对声波方程 ,给出了完全匹配层吸收边界方法。其基本思想是在所研究区域的边界上引入吸收层 ,波由研究区域边界传到吸收层时不产生任何反射 ,在吸收层内按传播距离的指数规律衰减 ,不产生反射 ,从而达到吸收边界的效果。文中给出了完全匹配层方法的基本原理 ,导出了声波方程完全匹配层的控制方程 ,并给出了令人满意的数值模拟结果。

黏弹性多层介质中SH波动的一种吸收边界条件

提出一种高精度时域吸收边界条件,与有限元法结合用于模拟瞬态标量SH波在达朗贝尔黏弹性多层介质中传播问题.建立时域吸收边界条件的过程是:首先将半无限域沿着竖向半离散得到半离散的位移方程以及人工边界处的力-位移关系,再通过引入模态分解,将物理空间下的量转化到模态空间,从而获得半无限域模态空间下的频域动力刚度;其次采用一种在全频范围内收敛的连分式精确逼近单层介质模态空间下标量形式的频域动力刚度,将标量连分式扩展为矩阵形式用来表示多层介质的频域动力刚度;最后通过引入辅助变量技术,将模态空间下基于连分式的频域动力刚度关系转化为时域吸收边界条件,进一步转换到物理空间后得到物理空间下的时域吸收边界条件.单层介质和五层介质的数值算例表明,建立的高精度时域吸收边界条件对于达朗贝尔黏弹性单层介质是精确且稳定的;对于达朗贝尔黏弹性多层介质,为了保证其高精度特性,需要将人工边界放置在距离感兴趣区域约为0.5倍无限域高度的位置处.