0 引言
目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法。Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射。之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果。有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果。但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] 。为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用。如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演。同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收。然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟。
本文基于Maxwell方程推导了一种在时域有限差分算法中运用递推卷积求解电磁场的方法,该方法在离散网格条件下通过递推形式计算卷积,不分裂变量,直接计算卷积,避免了对卷积求解的复杂计算。并将该方法运用到多源并发下对电磁波的吸收,为阵列探地雷达数值仿真做了铺垫。
1 基本原理
1.1 坐标伸缩变化下Maxwell卷积方程
非递归卷积完全匹配层主要求解是在麦克斯韦方程的频率域进行,对连续和时谐场,麦克斯韦旋度支配方程可写为:
(1) ∇ × H = j ω ε E + σ E = j ω ε - j σ ω E = j ω ε - E ,
(2) $\begin{aligned}\nabla \times \boldsymbol{E} &=-j \omega \mu \boldsymbol{H}+\sigma_{m} \boldsymbol{H} \\&=-j \omega\left(\mu-j \frac{\sigma_{m}}{\omega}\right) \boldsymbol{H}=-j \omega \bar\mu{-} \boldsymbol{H}\end{aligned}$
在PML中电导率不是由物理空间模型的基本参数决定的,而是设置电导率参数使网格截断的反射最小。在此意义上,电导率σ 具有任意性,将一个从属于特定位置的相对介电常数来规范电导率,在分裂场PML条件下定义:
(3) S w = 1 + σ j ω ε 0 ,
当频率趋于0时,Sw 没有意义。这将导致低频引发数值异常,为修正这个问题,增加额外的参数来保证频率趋于0时Sw 值的有限性。引入更一般化的坐标伸缩因子Si :
(4) S i = K i + σ i α i + j ω ε 0 ,
式中:Ki 是改善PML对表面波吸收的参数;σi 为PML层内k 方向电导率参数;αi 为改善PML对低频分量的吸收参数;σi 和αi 都为大于零的正实数;Ki ≥1。在三维空间下,i ∈{x ,y ,z },每一个Si 项总是与i ∈{x ,y ,z }方向的道数是成对的。有耗介质中,在TM极化模式下(Hx ,Hy ,Ez ),三维空间下电场频域方程表达式为:
(5) j ω ε E z + σ E z = 1 S x ∂ H y ∂ x - 1 S y ∂ H x ∂ y 。
将方程(5)通过傅里叶逆变换转化到时域,由于坐标伸缩因子与频率无关,所以在时域内方程右边存在卷积,即:
(6) ε ∂ ∂ t E z + σ E z = s - x ( t ) * ∂ ∂ x H y - s - y ( t ) * ∂ ∂ y H x ,
式中:* 表示卷积运算;S - x 和S - y 分别为1 S x 和1 S y 的逆傅里叶变换,即S - i =F -1 1 S i 。坐标伸缩因子Si 的倒数为:
(7) 1 S i = 1 K i + σ i α i + j ω ε 0 = α i + j ω ε 0 α i K i + σ i + j ω ε 0 ,
(8) F - 1 b d + a c - b d 1 + j ω d c = F - 1 ε 0 K i ε 0 + α i α i K i + σ i - ε 0 K i ε 0 1 + j ω K i ε 0 α i K i + σ = 1 K i δ ( t ) - σ i K i 2 ε 0 e x p - t α i ε 0 + σ i K i ε 0 u ( t ) ,
其中δ (t )是Dirac冲击函数,u (t )是单位阶跃函数。定义:
(9) ζ i ( t ) = - σ i K i 2 ε 0 e x p - t α i ε 0 + σ i K i ε 0 u ( t ) , i = ( x , y , z ) 。
(10) S - i = δ ( t ) K i - σ i ε 0 K i 2 e - t α i ε 0 + σ i K i ε 0 u ( t ) = δ ( t ) K i + ζ i ( t ) 。
因为Dirac冲击函数和其他函数的卷积仍为原函数,即δ (t )*f (t )=f (t ),将式(10)代入式(6)有:
(11) ε r ε 0 ∂ E z ∂ t + σ E z = 1 K x ∂ ∂ x H y - 1 K y ∂ ∂ y H x + ζ x ( t ) * ∂ ∂ x H y - ζ y ( t ) * ∂ ∂ y H x 。
式(11)为电场z 分量的时域递推卷积方程,直接在递推方程(14)中计算时域卷积效率是很低的,为有效计算卷积,下面给出递归卷积求解原理和在时域有限差分方法中的实现步骤。
1.2 基于递归卷积的非分裂场PML的FDTD实现
时域有限差分方法是对时间和空间进行离散化,在时间相差半个步长上迭代求解电场和磁场值。对式(11)的卷积求解,定义函数ψ E u i q 为:
(12) $\begin{aligned} \psi_{E_{u} i}^{q} &=\left.\zeta_{i}(t) * \frac{\partial H_{v}}{\partial i}\right|_{t=q \Delta t} \\ &=\int_{0}^{q \Delta t} \zeta_{i}(\tau) \frac{\partial H_{v}(q \Delta t-\tau)}{\partial i} \mathrm{~d} \tau, \end{aligned}$
式中:Eu i 表明这个函数将会在电场Eu 分量上更新,并与i (i =x ,y ,z )方向的空间导数有关。假设积分变量τ 在式(12)中是连续变化的,由于时域有限差分方法的离散性,所以Hv 场是离散变化的。用连续变量t 表示Hv 时,只取离散值,即
(13) ∂ H v ( t ) ∂ i = ∑ i = 0 I m a x ∂ H v ( k Δ t ) ∂ i p k ( t ) 。
令f (t )=∂ H v ∂ i ,此函数的阶梯表示如图1 a所示。假设函数在第一个时间步时为0,在t =Δt 时刻函数值为f 1 并保持不变,直到t =2Δt 时函数值变为f 2 ,当t =3Δt 时函数值为f 3 ,以此类推得到在任意时刻的磁场随空间i 分量上的变化值。
图1
图1
函数f (t )阶梯表示法
Fig.1
Stepped representation of function f (t )
卷积包含f (q Δt -τ )函数。在时间步为q =0时,函数值为f (-τ ),如图1 b所示。在所有采样点fn 都关于原点反转对称,在-Δt ≤τ <0上函数值为f 1 ,在-2Δt ≤τ <-Δt 上函数值为f 2 ,在-3Δt ≤τ <-2Δt 上函数值为f 3 ,其余各时刻的函数值以此类推。
图1 c展示了在特定时刻q =4时f (q Δt -τ )的函数值。如图所示,扩展到τ =0右侧右边的第一个脉冲具有fq 的值,扩展到τ =Δt 右边的值为fq -1 ,扩展到τ =2Δt 右边的值为fq -2 ,以此类推。此时移动函数可表示为:
(14) f ( q Δ t - τ ) = ∑ k = 0 q - 1 f q - k p k ( τ ) ,
(15) ∂ H v ( q Δ t - τ ) ∂ i = ∑ k = 0 q - 1 ∂ H v q - k ∂ i p k ( τ ) ,
其中:H v q - k 表示磁场在时间步q -k 时刻的值,空间导数用空间有限差分来实现。在时间步q 时刻,ψ E u i q 变换积分和求和顺序可得到:
(16) $\begin{aligned}\psi_{E_{u} i}^{q} &=\sum_{k=0}^{q-1} \frac{\partial H_{v}^{q-k}}{\partial i} \int_{0}^{q \Delta t} \zeta_{i}(\tau) p_{k}(\tau) \mathrm{d} \tau \\&=\sum_{k=0}^{q-1} \frac{\partial H_{v}^{q-k}}{\partial i} \int_{k \Delta t}^{(k+1) \Delta t} \zeta_{i}(\tau) \mathrm{d} \tau\end{aligned}$
式(16)中,脉冲函数pk (t )用来确定积分上下限。考虑如下特殊积分:
(17) - ∫ k Δ t ( k + 1 ) Δ t e - a t d t = 1 a e - a t | k Δ t ( k + 1 ) Δ t = 1 a [ e - a ( k + 1 ) Δ t - e - a k Δ t ] = 1 a ( e - a Δ t - 1 ) ( e - a Δ t ) k ,
(18) $\begin{aligned}Z_{0 i}(k) &=\int_{k \Delta t}^{(k+1) \Delta t} \zeta_{i}(\tau) \mathrm{d} \tau \\&=-\frac{\sigma_{i}}{\varepsilon_{0} K_{i}^{2}} \int_{k \Delta t}^{(k+1) \Delta t} \mathrm{e}^{-\left(\frac{\sigma_{i}}{\varepsilon_{0} K_{i}}+\frac{\alpha}{\varepsilon_{0}}\right) \tau} \mathrm{d} \tau \\&=a_{i} \mathrm{e}^{-\left(\frac{\sigma_{i}}{K_{i}}+\alpha_{i}\right) \frac{k \Delta t}{\varepsilon_{0}}}=a_{i}\left(b_{i}\right)^{k},\end{aligned}$
(19) a i = σ i ( σ i K i + K i 2 α i ) ( e - σ i K i + α i Δ t ε 0 - 1.0 ) i = ( x , y , z ) ,
(20) b i = e x p - σ i ε 0 K i + α ε 0 Δ t i = ( x , y , z ) ,
将式(17)、(18)、(19)、(20)代入式(16)中有:
(21) ψ E u i q = ∑ k = 0 q - 1 ∂ H v q - k ∂ i a i ( b i ) k ,
(22) ψ E u i q = a i ∂ H v q ∂ i + ∑ k = 1 q - 1 ∂ H v q - k ∂ i a i ( b i ) k ,
(23) ψ E u i q = a i ∂ H v q ∂ i + ∑ n = 0 q - 2 ∂ H v q - n - 1 ∂ i a i ( b i ) n + 1 ,
(24) ψ E u i q = a i ∂ H v q ∂ i + ∑ k = 0 [ q - 1 ] - 1 ∂ H v q - 1 - k ∂ i a i ( b i ) k ,
(25) ψ E u i q = a i ∂ H v q ∂ i + b i ψ E u i q - 1 ,
根据式(25),ψ E u i q 在时间步q 是q -1时刻的函数,因此ψ E u i q 求解可以通过递推求解。
空间离散下的电场递推式,按照Yee氏网格将式(25)进行时间和空间上的离散,可得:
(26) ε r ε 0 E z n + 1 ( i , j , k ) - E z n ( i , j , k ) Δ t + σ E z n + 1 ( i , j , k ) - E z n ( i , j , k ) 2 = H y n + 1 / 2 ( i , j , k ) - H y n + 1 / 2 ( i - 1 , j , k ) K x Δ x - H x n + 1 / 2 ( i , j , k ) - H x n + 1 / 2 ( i - 1 , j , k ) K y Δ y + ∑ m = 0 N - 1 Z 0 x ( m ) H y n - m + 1 / 2 ( i , j , k ) - H y n - m + 1 / 2 ( i - 1 , j , k ) Δ x + ∑ m = 0 N - 1 Z 0 y ( m ) H x n - m + 1 / 2 ( i , j , k ) - H x n - m + 1 / 2 ( i , j - 1 , k ) Δ y
上式中Z 0 i (k )项含有卷积的计算,离散卷积计算十分复杂,但是同时Z 0 i (k )项是简单的指数形式,从而它们的和可以通过递归卷积来得到,引入一组新的辅助表达式φi :
(27) φ e z i n + 1 / 2 ( i , j , k ) = b i φ e z i n - 1 / 2 ( i , j , k ) + a i [ H i n + 1 / 2 ( i , j , k ) - H i n + 1 / 2 ( i - 1 , j , k ) ] / Δ i ( i = x , y )
将φi 代入式(26)经整理后得到电场Ez 分量的递推表达式:
(28) E z n + 1 ( i , j , k ) = 1 - σ ( m ) Δ t 2 ε ( m ) 1 + σ ( m ) Δ t 2 ε ( m ) E z n ( i , j , k ) + Δ t ε ( m ) 1 + σ ( m ) Δ t 2 ε ( m ) [ φ e z x n + 1 / 2 ( i , j , k ) - φ e z y n + 1 / 2 ( i , j , k ) ] + Δ t ε ( m ) 1 + σ ( m ) Δ t 2 ε ( m ) H y n + 1 / 2 ( i , j , k ) - H y n + 1 / 2 ( i - 1 , j , k ) K x Δ x - H x n + 1 / 2 ( i , j , k ) - H x n + 1 / 2 ( i , j - 1 , k ) K y Δ y ,
式(28)中:σ 为电导率,ε 为介质的介电常数,m =(i ,j ,k ),在仿真时令模拟区域的K 值为1,在PML层内使K 值渐进变化,PML层为有限的厚度单元,σ 、K 和α 在PML层中单调变化:
(29) σ z ( z ) = σ m a x | z - z 0 | m d m ,
(30) K z ( z ) = 1 + ( K m a x - 1 ) | z - z 0 | m d m ,
(31) α z ( z ) = α m a x | z - z 0 | d ,
式中:d 为PML层的厚度,z 0 为PML层靠近FDTD区域的界面位置,在仿真时取m =4。σ max 的最佳值可取为σ max =m + 1 ε r 150 π δ ,其中δ 为元胞尺寸。K max =5~ 11,α max =0~ 0. 05。根据探地雷达的频率段,K max 取5左右,α max 取0. 008为适当值。
1.3 Mur二阶吸收边界条件的FDTD实现
(32) f ( x , y , t ) = A e x p [ j ( ω t - α x x - α y y ) ] ,
假设在x =0的位置设置截断边界,则在x ≥0的区域会同时存在入射波和反射波。因此式(32)可写为:
(33) f ( x , y , t ) = A - e x p [ j ( ω t + α 2 - α y 2 x + α y y ) ] + A + e x p [ j ( ω t + α 2 - α y 2 x + α y y ) ] ,
可将上式分解为左行波f - (入射波)和右行波f + (反射波)之和:
(34) f - = A - e x p [ j ( ω t + α 2 - α y 2 x + α y y ) ] ,
(35) f + = A + e x p [ j ( ω t + α 2 - α y 2 x + α y y ) ] 。
(36) L ± = ∂ ∂ x ± 1 v 2 ∂ 2 ∂ t 2 - ∂ 2 ∂ y 2 ,
对左行波使得L - f - |x =0 中并令f =Ez ,化简后得到:
(37) 1 v ∂ 2 ∂ x ∂ t - 1 v 2 ∂ 2 ∂ t 2 + 1 2 ∂ 2 ∂ y 2 E z = 0
根据∂ E z ∂ y =-μ ∂ H x ∂ t ,对二维电磁场问题,Mur指出二阶吸收边界条件可降低为只含电场(E )和磁场(H )分量的一阶导数:
(38) ∂ E z ∂ x - 1 v ∂ E z ∂ t - v μ 2 ∂ H x ∂ y = 0 ,
在(i +1/ 2,j )点和t =(n +1/ 2)Δt 时刻对式(62)进行差分离散并进行线性插值可得到左截断边界处Mur的二阶吸收边界条件递推式:
(39) E z n + 1 ( i , j ) = E z n ( i + 1 , j ) + v Δ t - Δ x v Δ t + Δ x [ E z n + 1 ( i + 1 , j ) - E z n ( i , j ) ] - v 2 μ Δ t 2 ( v Δ t + Δ x ) · Δ x Δ y × [ H x n + 1 / 2 i , j + 1 2 - H x n + 1 / 2 i , j - 1 2 ] + H x n + 1 / 2 i + 1 , j + 1 2 - H x n + 1 / 2 i + 1 , j - 1 2 ,
式(39)中,v 为电磁波传播速度,μ 为磁导率,Δx 、Δy 、Δt 分别为x 、y 方向的空间步长和时间步长。式(39)给出的Mur二阶吸收边界条件递推公式没有在边角点处对电场值进行修正,因此,在4个边角位置处发生明显的波反射,如图2 所示,给出了激励源位于模型区域中心位置处,Ez 波场在600个时间步的变化,在第400个时间步时,由于没有对边角进行修正,此时4个角点位置处都发生了明显的反射,在500和600时间步时刻的研究区域内存在4个角点反射的波场。
图2
图2
未修正角点不同时间步的Ez 波场
Fig.2
Ez wave field at different time steps of uncorrected corners
1.4 Mur二阶吸收边界条件的角点修正
二维仿真时,研究区域是矩形,此时在4个边角点位置处出现场重合的情况,因此需要对角点进行处理,以TM波为例,假设角点处存在沿元胞对角线传播的行波。如图3 所示,在距离坐标(i 0 ,j 0 )为2v Δt 处的对角线上取一点P ,根据FDTD的离散稳定条件,有2v Δt ≤2 δ ,其中δ =Δx =Δy ,沿对角线传播的行波:
(40) E z n ( i 0 , j 0 ) = E z n - 2 ( P ) ,
创建新坐标系ηoξ ,新坐标系于原始坐标系逆时针旋转45° ,在新坐标系下Mur一阶近似公式可表示为:
(41) ∂ E z ∂ ξ - 1 v ∂ E z ∂ t = 0 。
图3
图3
TM波左下角点修正示意
Fig.3
Schematic diagram of TM wave lower left corner correction
由于P 点与原点在同一离散网格内,忽略波振幅衰减,利用线性插值可得到:
(42) E z n - 2 ( i 0 , j 0 ) - E z n - 2 ( P ) 2 v Δ t = E z n - 2 ( i 0 , j 0 ) - E z n - 2 ( i 0 + 1 , j 0 + 1 ) 2 δ 。
根据式(40)结合式(42),当Δx ≠Δy 时,有:
(43) E z n + 1 ( i 0 , j 0 ) = E z n ( i 0 + 1 , j 0 + 1 ) + v Δ t - ( Δ x ) 2 + ( Δ y ) 2 v Δ t + ( Δ x ) 2 + ( Δ y ) 2 · [ E z n + 1 ( i 0 + 1 , j 0 + 1 ) - E z n ( i 0 , j 0 ) ] ,
根据式(43),修正角点位置处的电场值取决于该角点相邻一个时间步位置处前一时刻的电场值和该点的速度值。
2 模型验证
2.1 均匀介质
为验证在均匀介质下Mur二阶近似吸收边界条件和基于递归卷积完全匹配层作为吸收边界条件在多源并发下对电磁波的吸收效果,设置模拟区域网格大小为200×200,四周PML层厚度为10个网格。模拟区域的物性参数为:εr =3. 5,ur =1,σ =2. 5 mS/m。两个激励源采用主频为900 MHz的布莱克曼—哈里斯脉冲,两个激励源的位置沿x 轴相距80个空间步长,空间步长为0.006 m,时间步长为 0.015 ns。不同时间步的Ez 波场如图4 所示。图4 为应用递归卷积完全匹配层作为吸收边界条件的Ez 波场快照,仿真时间持续400个时间步长。在第100时间步时(图4 a),两个并发激励源电磁波相遇,经相互干涉作用后继续扩散;在时间步200时(图4 b),部分电磁波进入PML层被吸收掉; 时间步300时(图4 c)两个激励源产生的扩散电磁波到达彼此对面PML层,此时部分波到达边界进入PML层无任何反射发生;时间步400时(图4 d),波扩散离开研究区域,在整个仿真持续时间段内,无任何明显反射发生,电磁场在网格截断位置处被PML层有效吸收掉了。
图4
图4
均匀介质递归卷积完全匹配层不同时间步Ez 波场快照
Fig.4
Snapshots of Ez wave field at different time steps of the homogeneous medium recursive convolution perfectly matched layer
图5 给出了Mur二阶吸收边界条件不同时间步Ez 波场快照,二阶Mur吸收边界条件是运用微分二阶近似来求解单向行波方程,它是将微分方程进行泰勒级数展开后去掉高阶项的近似算法,此时反射系数不为零,因此不能实现对波的良好吸收。在两个激励源并发的条件下,电磁波相互干涉,重构后的波场形成小方形包络,如图5 c箭头所示位置,随着时间推移,小方形包络在水平方向上扩展开来,以平面波束形式达到网格边界位置。这时,Mur二阶吸收边界对电磁波吸收效果较差,会在边界处发生反射。在前300个时间步内(图5 a,5b,5c),由单个激励源产生的电磁波在均匀介质中扩散,以球面的形式到达网格边界处, 在第400个时间步时(图5 d),经过干涉后的波到达边界,此时部分波返回研究区域;在时间步500时,波已经扩散离开研究区,然而因吸收不完全导致部分反射波返回研究区域,如图5 e、5f所示,此时波在边界位置发生多次反射,形成谐振,对下一时间段内的有效波场造成二次干扰。
图5
图5
均匀介质Mur二阶吸收边界不同时间步Ez 波场快照
Fig.5
Snapshots of Ez wave field at different time steps of Mur second-order absorbing boundary in homogeneous medium
2.2 低洼模型
在复杂地质构造条件下,在地下介质空间中扩散的电磁波会产生反射、绕射现象。绕射波会以不同的角度入射到匹配层中,为体现两种不同的吸收边界对波的吸收效果,构造了复杂低洼地质模型,模拟现实地质环境中的正、逆断层。设置模型区域网格大小为300×400,模拟区域的物性参数为:上层介质相对介电常数ε 1 =8,电导率σ 1 =0.003 S/m,下层介质相对介电常数ε 2 =5,电导率σ 2 =0.001 S/m。所有激励源采用主频为900 MHz的布莱克曼—哈里斯脉冲,空间步长为0.006 m,时间步长为 0.015 ns,26个激励源的位置沿x 轴分布,每个激励源相距5个空间步长,其模型的空间分布如图6 所示。对于递归卷积完全匹配层在四周设置10个网格大小的PML层厚度。为方便更加清晰显示波场信息,下面将使用灰度图显示不同时间步Ez 波场的信息。
图6
图6
低洼模型结构示意
Fig.6
Schematic diagram of low-lying model structure
在多个激励源并发的情况下,采用相同主频的脉冲电磁波,此时在近地表区域,由单个激励源产生的球面波在空间上进行波场重构,子波波前形成平面电磁波束,如图7 a,7b所示。平面电磁波束在地下空间扩散,在第400个时间步时遇到第二层介质产生绕射波,如图7 c箭头指示位置;第500个时间步时,低洼模型底部产生的反射波如图7 d标记所示,此时反射波的能量明显高于绕射波能量;两种吸收边界条件下,电磁波在介质区域内传播规律相同,有相同的波场快照;在第700个时间步时,平面电磁波束到达底部边界,基于递归卷积完全匹配层作为吸收边界条件的波场快照如图8 a所示,此时平面电磁波束进入完全匹配层被有效地吸收掉,而以二阶Mur作为吸收边界条件在底部网格截断处发生了反射,部分波返回介质区域,如图8 b圆圈位置,这将对能量相对较弱的绕射波造成了二次干扰,影响有效波场。
图7
图7
低洼模型不同时间步Ez 波场快照
Fig.7
Snapshots of the Ez wave field of the low-lying model at different time steps
图8
图8
低洼模型700时间步Ez 波场快照
a—递归卷积完全匹配层波场快照;b—Mur二阶吸收边界波场快照
Fig.8
Snapshot of the low-lying model at 700 time step Ez wave field
a—recursive convolution complete matching layer wave field snapshot;b—snapshot of Mur second-order absorbing boundary wave field
为了对比在不同偏移距下,两种吸收边界条件对多源并发电磁波的吸收效果,设置了宽角法观测低洼模型,将26个发射天线(激励源)置于低洼位置正上方,每个发射天线间隔两个空间步长,接收天线沿x 轴移动,采集200道雷达数据。天线参数依然采用900 MHz主频的布莱克曼—哈里斯脉冲,多个发射天线无延时同时发射相同信号脉冲电磁波,得到在多源并发下,以递归卷积完全匹配层作为吸收边界条件的数据记录如图9 a所示。此时,由于采用的是相同极化方向的激励信号,直达波、反射波与单个激励源雷达数据记录(如图10 a所示)相比加大了子波时宽,增强了来自地下界面的反射信号,此时在直达波左、右两翼无明显畸变,无虚假反射。而Mur二阶吸收边界条件在大偏移距处有强烈的虚假反射记录;这是由于偏移距越大,波入射到匹配层的入射角度越大,掠射情况越严重。而在多个激励源的情况下,在大偏移距下,虚假反射会更加明显;这是由于加大了子波的时宽和介质的不均一性导致的色散造成的。如图9 b箭头所指示的位置,在直达波的左右翼都发生了明显的波形畸变和虚假反射记录。虽然在单个天线时的雷达数据记录上,直达波左右两翼畸变不明显,但是也产生了明显的虚假反射信号,如图10 b箭头所指示的位置。
图9
图9
多源并发雷达数据记录
a—递归卷积完全匹配层数据记录;b—Mur二阶吸收边界数据记录
Fig.9
Multi-source concurrent radar data recording
a—recursive convolution complete matching layer data record;b—Mur second-order absorbing boundary data record
图10
图10
单个激励源雷达数据记录
a—递归卷积完全匹配层数据记录;b—Mur二阶吸收边界数据记录
Fig.10
Single excitation source radar data record
a—recursive convolution complete matching layer data record;b—Mur second-order absorbing boundary data record
3 结论与展望
1)在均匀介质下,Mur二阶吸收边界条件对多源并发下电磁波的吸收不完全,会在网格截断处产生反射回波返回波场模拟区域。而基于非分裂的递推卷积完全匹配层能很好地吸收多源并发下的电磁波,有效改善边界位置处对干涉后的电磁波吸收效果。
2)在复杂地质构造下,电磁波在地下介质中发生反射、绕射后以不同的入射角度到达网格截断处,传统的Mur二阶吸收边界条件在网格截断位置直接对微分方程进行近似求解,不能很好地消除强电磁反射,在边界网格截断位置处会有部分反射波返回研究区域对有效波造成二次干扰,特别是在大偏移距下,会造成直达波波形畸变,形成虚假反射记录;而基于非分裂递归卷积的完全匹配层作为吸收边界条件能很好地抑制直达波波形畸变现象,消除虚假反射记录。
3)基于非分裂递归卷积的完全匹配层作为吸收边界条件在不分裂波场的情况下,又避免了直接对卷积的求解的复杂运算,提升了计算的效率,改善了吸收的效果。
本文详细介绍了非分裂递归卷积完全匹配层和Mur二阶吸收边界结合时域有限差分的实现方法,并将它们运用到多个激励源并发下探地雷达正演模拟,基于C语言和matlab平台开发了相应的数值计算程序,为阵列探地雷达的数值模拟研究提供了帮助。
致谢
感谢煤炭资源与安全开采国家重点实验室提供的计算设备,感谢评审的专家提出的宝贵意见。
参考文献
View Option
[1]
Bayliss A , Turkel E . Radiation boundary conditions for wave-like equations
[J]. Communications on Pure and Applied Mathematics , 1980 , 33 (6 ):707 -725 .
[本文引用: 1]
[2]
Engquist B , Majda A . Absorbing boundary conditions for numerical simulation of waves
[J]. Proceedings of the National Academy of Sciences of the United States of America , 1977 , 74 (5 ):1765 -1766 .
[本文引用: 1]
[3]
Mur G . Absorbing boundary conditions for the finite-difference approximation of the time-domain electromagnetic-field equations
[J]. IEEE Transactions on Electromagnetic Compatibility , 1981 , 23 (4 ):377 -382 .
[本文引用: 1]
[4]
Berenger J P . Three-dimensional perfectly matched layer for the absorption of electromagnetic waves
[J]. Journal of Computational Physics , 1996 , 127 (2 ):363 -379 .
[本文引用: 1]
[5]
Chen Y H , Weng C C , Oristagilo M L . Application of perfectly matched layers to the transient modeling of subsurface EM problems
[J]. Geophysics , 1997 , 97 (6 ):1730 -1736 .
[本文引用: 1]
[6]
Festa G . Interaction between surface waves and absorbing boundaries for wave propagation in geological basins:2D numerical simulations
[J]. Geophysical Research Letters , 2005 , 32 :L20306 .
[本文引用: 1]
[7]
Qin Z , Lu M , Zheng X , et al . The implementation of an improved NPML absorbing boundary condition in elastic wave modeling
[J]. Applied Geophysics , 2009 , 6 (2 ):113 -121 .
[本文引用: 1]
[8]
Higdon R L . Absorbing boundary conditions for elastic waves
[J]. Siam Journal on Numerical Analysis , 2012 , 31 (2 ):64 -100 .
[本文引用: 1]
[9]
Liao Z P , Wong H L , Yang B P , et al . A transmitting boundary for transient wave analysis
[J]. Science in China Series A-Mathematics Physics Astronomy & Technological Science , 1984 , 27 (10 ):1063 -1076 .
[本文引用: 1]
[10]
Marfurt K J . Accuracy of finite difference and finite element modeling of the scalar and elastic wave equation
[J]. Geophysics , 1984 , 49 (4 ):533 -549 .
[本文引用: 1]
[11]
丁科 . PML吸收边界条件影响因素分析
[J]. 物探与化探 , 2012 , 36 (4 ):623 -627 .
[本文引用: 1]
Ding K . Analysis of influencing factors of PML absorption boundary conditions
[J]. Geophysical and Geochemical Exploration , 2012 , 36 (4 ):623 -627 .
[本文引用: 1]
[12]
Sacks S Z , Kingsland M D , Lee R , et al . A perfectly matched anisotropic absorber for use as an absorbing boundary condition
[J]. IEEE Transactions on Antennas & Propagation , 1995 , 43 (12 ):1460 -1463 .
[本文引用: 1]
[13]
Gedney D S . An anisotropic perfectly matched layer-absorbing medium for the truncation of FDTD lattices
[J]. IEEE Transactions on Antennas & Propagation , 2002 , 44 (12 ):1630 -1639 .
[本文引用: 1]
[14]
肖明顺 , 昌彦君 , 曹中林 , 等 . 探地雷达数值模拟的吸收边界条件研究
[J]. 工程地球物理学报 , 2008 , 5 (3 ):315 -320 .
[本文引用: 1]
Xiao M S , Chang Y J , Cao Z L , et al . Study on absorption boundary conditions for numerical simulation of ground penetrating radar
[J]. Journal of Engineering Geophysics , 2008 , 5 (3 ):315 -320 .
[本文引用: 1]
[15]
詹应林 , 昌彦君 , 曹中林 . 基于UPML吸收边界的探地雷达数值模拟研究
[J]. 资源环境与工程 , 2008 , 22 (2 ):235 -238 .
[本文引用: 1]
Zhan Y L , Chang Y J , Cao Z L . Numerical simulation of ground penetrating radar based on UPML absorption boundary
[J]. Resources Environment and Engineering , 2008 , 22 (2 ):235 -238 .
[本文引用: 1]
[16]
冯德山 , 谢源 . 基于单轴各向异性完全匹配层边界条件的ADI-FDTD三维GPR全波场正演
[J]. 中南大学学报:自然科学版 , 2011 , 42 (8 ):202 -210 .
[本文引用: 1]
Feng D S , Xie Y . Full Wave field forward modeling of 3D GPR based on boundary conditions of uniaxial anisotropic completely matched layer
[J]. Journal of Central South University:Science and Technology , 2011 , 42 (8 ):202 -210 .
[本文引用: 1]
[17]
冯德山 , 陈承申 , 戴前伟 . 基于UPML边界条件的交替方向隐式有限差分法GPR全波场数值模拟
[J]. 地球物理学报 , 2010 , 53 (10 ):2484 -2496 .
[本文引用: 1]
Feng D S , Chen C S , Dai Q W . Numerical simulation of GPR full wave field with alternate direction implicit finite difference method based on UPML boundary condition
[J]. Chinese Journal of Geophysics , 2010 , 53 (10 ):2484 -2496 .
[本文引用: 1]
[18]
冯德山 , 王向宇 . 基于卷积完全匹配层的旋转交错网格高阶差分法模拟弹性波传播
[J]. 物探与化探 , 2018 , 42 (4 ):766 -776 .
[本文引用: 1]
Feng D S , Wang X Y . Elastic wave propagation simulation in anisotropic media and random media using heig-order difference method of rotation staggered grids based on convolutional perfectly matched layer
[J]. Geophysical and Geochemical Exploration , 2018 , 42 (4 ):766 -776 .
[本文引用: 1]
[19]
李静 , 曾昭发 , 吴丰收 , 等 . 探地雷达三维高阶时域有限差分法模拟研究
[J]. 地球物理学报 , 2010 , 53 (4 ):974 -981 .
[本文引用: 1]
LI J , Zeng Z F , Wu F S , et al . Study of three dimension high-order FDTD simulation for GPR
[J]. Chinese Journal of Geophysics , 2010 , 53 (4 ):974 -981 .
[本文引用: 1]
[20]
Chew W C , Weedon W H . A 3D perfectly matched medium from modified maxwell’s equations with stretched coordinates
[J]. Microwave & Optical Technology Letters , 2010 , 7 (13 ):599 -604 .
[本文引用: 1]
[21]
Kuzuoglu M , Mittra R . Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers
[J]. IEEE Microwave & Guided Wave Lett , 1996 , 6 (12 ):447 -449 .
[本文引用: 1]
Radiation boundary conditions for wave-like equations
1
1980
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
Absorbing boundary conditions for numerical simulation of waves
1
1977
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
Absorbing boundary conditions for the finite-difference approximation of the time-domain electromagnetic-field equations
1
1981
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
Three-dimensional perfectly matched layer for the absorption of electromagnetic waves
1
1996
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
Application of perfectly matched layers to the transient modeling of subsurface EM problems
1
1997
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
Interaction between surface waves and absorbing boundaries for wave propagation in geological basins:2D numerical simulations
1
2005
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
The implementation of an improved NPML absorbing boundary condition in elastic wave modeling
1
2009
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
Absorbing boundary conditions for elastic waves
1
2012
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
A transmitting boundary for transient wave analysis
1
1984
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
Accuracy of finite difference and finite element modeling of the scalar and elastic wave equation
1
1984
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
PML吸收边界条件影响因素分析
1
2012
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
PML吸收边界条件影响因素分析
1
2012
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
A perfectly matched anisotropic absorber for use as an absorbing boundary condition
1
1995
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
An anisotropic perfectly matched layer-absorbing medium for the truncation of FDTD lattices
1
2002
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
探地雷达数值模拟的吸收边界条件研究
1
2008
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
探地雷达数值模拟的吸收边界条件研究
1
2008
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
基于UPML吸收边界的探地雷达数值模拟研究
1
2008
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
基于UPML吸收边界的探地雷达数值模拟研究
1
2008
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
基于单轴各向异性完全匹配层边界条件的ADI-FDTD三维GPR全波场正演
1
2011
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
基于单轴各向异性完全匹配层边界条件的ADI-FDTD三维GPR全波场正演
1
2011
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
基于UPML边界条件的交替方向隐式有限差分法GPR全波场数值模拟
1
2010
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
基于UPML边界条件的交替方向隐式有限差分法GPR全波场数值模拟
1
2010
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
基于卷积完全匹配层的旋转交错网格高阶差分法模拟弹性波传播
1
2018
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
基于卷积完全匹配层的旋转交错网格高阶差分法模拟弹性波传播
1
2018
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
探地雷达三维高阶时域有限差分法模拟研究
1
2010
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
探地雷达三维高阶时域有限差分法模拟研究
1
2010
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
A 3D perfectly matched medium from modified maxwell’s equations with stretched coordinates
1
2010
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...
Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers
1
1996
... 目前,在时域有限差分方法中运用边界条件来截断网格存在很多方法.Bayliss和Turkel[1 ] 采用外行波的模拟法作为吸收边界条件,实现了对电磁波的简单吸收;Engquist和Majda [2 ] 提出了基于单向波动方程的Engquist-Majda吸收边界条件;Mur[3 ] 给出了波动方程的各阶近似及差分形式,表明了在高阶近似时,对电磁波的吸收效果较好,然而对于各向异性介质和色散介质,以二阶近似Mur作为吸收边界条件会在网格截断处发生强反射.之后,Berenger[4 ] 提出了完全匹配层吸收边界条件,他将电磁波在吸收边界区域进行波场分裂,并对各个分裂场赋予不同的损耗值,这相当于在FDTD网格外加入一层吸收媒介,它具有不依赖于外向传播的电磁波入射角及频率的波阻抗,因此,进入完全匹配层内的电磁波快速地衰减,实现了对不同频率、不同角度的入射波较好的吸收效果.有很多研究发现[5 ⇓ -7 ] ,Berenger完全匹配层作为吸收边界条件比旁轴近似吸收边界条件、Higdon吸收边界条件[8 ] 、廖氏吸收边界条件[9 ] 和指数衰减吸收边界条件[10 ] 具有更好地吸收效果.但是,Berenger的PML理论违背了波的折射原理,在Maxwell方程中是不能实现的,同时,电磁场的分裂增加了数值计算的难度,计算效率低,而且只对行波具有吸收效果,对空域中衰减的隐失波、低频波以及入射角度较小的掠射波吸收效果差[11 ] .为了改善边界条件的吸收效果,Sacks等[12 ] ,Gedney等[13 ] 提出了单轴各向异性完全匹配层作为吸收边界条件,该方法在数学上与Berenger提出的经典完全匹配层等价,它是基于Maxwell方程的,不分裂波场,而是在吸收系数中引入线性吸收因子,实现了对隐失波、低频波以及在PML层内凋落波等干扰波的吸收,该方法作为良好的吸收边界条件在差分算法中被广泛使用.如肖明顺等[14 ] 、詹应林等[15 ] 将各向异性完全匹配层应用于二维探地雷达的正演中;中南大学冯德山等[16 ⇓ -18 ] 实现了在二维和三维空间中将各向异性完全匹配层引入交替方向隐式时域有限差分算法中(ADI-FDTD);吉林大学的李静等[19 ] 人开展了三维各向异性完全匹配层作为吸收边界条件的高阶时域有限差分正演.同一时间段内,Chew等[20 ] 提出了拉伸坐标完全匹配层理论;Kuzuoglu等[21 ] 在复频域内的完全匹配层变量Sk 中引进低频分量吸收参数,将复平面的极点从实轴移动到虚轴,加强了对低频波和隐失波的吸收.然而到现研究阶段,所有的吸收边界条件都是在单个激励源的条件下对电磁波进行数值模拟. ...