专利名称：二阶波动方程的非分裂完全匹配层的构造方法
技术领域：
本发明属于数值模拟仿真技术领域，具体为一种适用于二阶波动方程的非分裂完全吸收匹配层的构造方法。
背景技术：
基于波动方程的数值仿真中，需采用数值方法离散化波动方程。由于电脑内存和计算时间的限制，需截断计算区域。由此将带来边界反射的问题。通常，在截断边界处可设置吸收边界条件，以消除或减弱由于截断计算区域所引起的边界反射波的影响。完全匹配层(PML) [1]是目前仿真中运用最广的吸收边界之一。PML是一种特殊的介质层，它具有这样的性质，当入射波从仿真区域进入其内部区域时，不会发生反射，而在其内部的传播过程中迅速衰减，进而达到消除反射波的目的。由于PML最先是为一阶波动方程而设计，以往对PML的应用多局限于一阶形式的波动方程。对二阶方程的PML研究则较少。一般情况下，一 个二阶方程可分解为一阶方程组，然后便可使用PML。然而某些情况下必须使用二阶波动方程，如超声非线性成像模拟中运用广泛的全波模型和Westervelt模型[2]’[3]，因此找到一种有效的适用于二阶波动方程的PML对此类仿真具有重要意义。另外，二阶形式的方程由于只涉及一个主要变量，在迭代和计算上也具有优势。目前很多商业有限元仿真软件是以二阶方程为核心，二阶形式的方程更便于采用这些软件直接进行仿真。若将二阶方程化为一阶方程组，则会破坏原来方程的结构，增加计算的复杂度。现有的几种适用于二阶波动方程的PML主要有两种类型。第一种是基于场量分裂_[6]的方法。该方法将声压分为若干部分，然后分别处理各个部分。其主要缺点为所需的额外存储量较大，计算量大。第二种是一种混合形式(mixed form)的方法，它形式上是二阶的，但执行以一阶的方式进行。该方法的主要缺点是破坏了原二阶方程的执行结构，处理不便，计算量较大。本发明克服了已有方法的缺点，提供了一种执行效率高，存储量低的PML，以解决二阶波动方程仿真中的边界问题。
发明内容
本发明的目的是提出一种执行效率高，存储量低的构造二阶波动方程的完全匹配层(PML)的方法，以解决基于二阶波动方程的仿真中的吸收边界问题。本发明提供的适用于二阶波动方程完全匹配层的构造方法，其步骤如下
1、基于伸缩坐标变换，采用直接求导的方法获得二阶波动方程的PML的频域表达式；
2、以角频率为中心变量，进行因式分解，逐步引入辅助变量，构造辅助微分方程，得到时域表达式。下面对各步骤进行具体描述。步骤一、基于伸缩坐标变换，采用直接求导的方法获得二阶波动方程的PML的频域表达式；二阶声场波动方程的频域表达式为
权利要求
1.一种二阶波动方程的非分裂完全匹配层的构造方法，其特征在于具体步骤为(O基于伸缩坐标变换，采用直接求导的方法，构造二阶波动方程的PML的频域表达式； (2)以角频率为中心变量，进行因式分解，逐步引入辅助变量，构造辅助微分方程，得到时域表达式。
2.根据权利要求I所述的构造方法，其特征在于构造二阶波动方程的PML频域表达式的步骤为 对于空间坐标，通过直接微分求导的方法得到PML的频率表达式
3.根据权利要求2所述的构造方法，其特征在于得到时域表达式的具体步骤为 对方程(2)右边各项以角频率为中心进行因式分解
4.根据权利要求2或3所述的构造方法，其特征在于对坐标变换的其它算子
全文摘要
本发明属于数值仿真技术领域，具体为一种构造二阶波动方程的非分裂完全吸收匹配层的方法。该方法直接采用坐标变换得到二阶方程的完全匹配层的频域表达式；对得到的频域表达式以角频率为中心进行分式分解，调整其结构；通过引入辅助变量，构造辅助微分方程的方法，得到简洁的完全匹配层的时域表达式。在本方法中，引入的辅助微分方程都为形式相同的一阶微分方程，大幅降低仿真中的离散化难度，提高了执行效率，降低了匹配层执行所需的存储量。本发明能作为一种有效边界条件，应用于基于二阶方程的数值仿真中。
文档编号G06F17/14GK102880590SQ20121036035
公开日2013年1月16日 申请日期2012年9月25日 优先权日2012年9月25日
发明者余锦华, 马友能, 汪源源 申请人:复旦大学